COMSOL Forums: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/
Most recent forum messagesThu, 03 Jan 2013 21:12:36 +0000COMSOL Forums: Problems generating lamb waves in an Aluminum platehttp://www.comsol.com/shared/images/logos/comsol_logo.gif
http://www.comsol.com/community/forums/general/thread/24175/
Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p93103
Rohan Repale,I just find I have the same problem as you ,could you share your model with me about this issue ,I am exceedly need your help.because I want to do PZT transducer to excite lamb in 3D model.what's more,I need to control the mode of lamb wave.could you help me ?you can send me related materials as more as possible.I am really appreciate your help.you can contact me via maochengtao@shu.edu.cn<br />
<br />
God bless you! Thank you!Thu, 03 Jan 2013 21:12:36 +000024175.1357247556.93103Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p72842
Hi Glenston<br />
<br />
Hope everything went fine with you.<br />
<br />
I am working on the same problem, did you find something helpful ?!<br />
<br />
WaelFri, 09 Mar 2012 13:14:33 +000024175.1331298873.72842Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69257
Hi<br />
<br />
that KB article is related first of all to RF where the meshing is probably the most delicate, because of the short wavelength of the EM waves (particularly in dense medium. <br />
<br />
But the principle remains valid for all FEM (just as for digital signal sampling where most engineers accepted this now, without even thinking), <br />
i,.e. in HT its the heat diffusivity alpha=k/rho/Cp that relates to the CFL number (mesh size time step physical property link), in CFD among others the visosity and Reynolds number for turbulence aspects. In structural p waves its the sqt(E/rho)=c speed of sound that links in material propewrties to the wavelength, hence the mesh density to resolve these ...<br />
<br />
--<br />
Good luck<br />
IvarFri, 27 Jan 2012 12:27:45 +000024175.1327667265.69257Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69255
Hi<br />
<br />
Thanks a lot for your reply and effort.<br />
This will definitely help me a lot.<br />
Will refer to the user guide for more information.7<br />
<br />
Thanks a lot again.<br />
<br />
Regards<br />
GlenstonFri, 27 Jan 2012 12:10:24 +000024175.1327666224.69255Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69253
Hi<br />
<br />
the first is related to the variable name, you define a parameter "w" but COMSOL uses "w" for the thrd space displacement (dependent variable). You may argue that you are in 2D, but COMSOL works always in 3D so it's dragging internally all variables, also for 3D, even if they are mostly set to "0" (but I do not believe its solving or uneccesary variables)<br />
<br />
Try to avoid using COMSOL internal names for your own use, add a letter or a "_" or whatever<br />
<br />
Read carefully the doc on the solvers, it's scattered but you have a lot in the user guide and more in the reference guide.<br />
<br />
When you define a time step range, these are the values COMSOL will extract from the time stepping it decides (in free mode), but COMSOL might decide by its own to do two steps and extract your values by interpolation in between the two steps. This is mostly OK for diffusion models or critical damped models, but for any periodical signal it's a catastrophy, COMSOL might fully overlook the oscillatory behaviour of a BC condition. therefore you should use "intermediate" or "strict", there Comsol ensure it's taking at least one step inbetween your define time steps, or strict when it takes at least a step AT your defined value (you see that by the - out" log messages).<br />
<br />
Depending on the solver BFD, gen-alpha ... COMSOL has different ways to estimate what is reasonably the "next step" via the Jacobian or other resolution parameter (i.e. Newton method ...) This chapter is rather lengthy to explain, and I'm not that comfortable with it, as there are many buttons to tweak and I have not had time to go systematically through them all ;) <br />
<br />
COMSOL has a rather good course on solvers, but even a full day is by far not enough to go through all the buttons<br />
<br />
--<br />
Good luck<br />
IvarFri, 27 Jan 2012 11:59:18 +000024175.1327665558.69253Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69249
Hi Ivar,<br />
<br />
When i wanted to mention the time steps.<br />
I used it referring to the article in the knowledge base<br />
<br />
<a href="http://www.comsol.com/support/knowledgebase/1118/" title="www.comsol.com/support/knowledgebase/1118/">www.comsol.com/support/knowledgebase/1118/</a><br />
This uses a manual time stepping and is based on the CFL criteria.<br />
Which is right?<br />
<br />
Should i be getting a -out for each and every step?<br />
Does that mean my output is right?<br />
<br />
Regard<br />
Glenston Fri, 27 Jan 2012 11:45:33 +000024175.1327664733.69249Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69236
<div class="quote"><br />
Hi<br />
<br />
Two important things to make your model work:<br />
<br />
as I hadn't open your latest models (I'm not always by my WS, often on travel or in meetings) I hadn't noticed <br />
<br />
1) your w for the 2*pi*f is already defined by COMSOL, for the dependent variable 3rd dimension "w"<br />
2) more important, you left the time stepping in "free" mode, you should use "strict" or "intermediate" then you will get far more interesting results ;)<br />
<br />
Note the difference in time stepping, you see it in your log file: there are several " - out" without any steps in between, this is not compatible with a sinusoidal BC entry, te default free is for time decaying processes such as diffusion<br />
<br />
Rule: if any BC is peridic & you are in time stepping => always use strict or intermediate (never "free")<br />
<br />
--<br />
Good luck<br />
Ivar<br />
</div><br />
<br />
Hi Ivar,<br />
<br />
Thanks a lot for your reply.<br />
I din not quite understand your first suggestion. Could you please elaborate<br />
<br />
I have a simple question whats the significance of time stepping and time steps set in the <br />
time edit field.<br />
<br />
How does the internal time step have a role to play in the formation of a solution.<br />
Wont manual time stepping with a very small value suffice.?<br />
<br />
Thanks a lot for your time and effort.<br />
<br />
Regards<br />
GlenstonFri, 27 Jan 2012 02:03:44 +000024175.1327629824.69236Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69201
Hi<br />
<br />
Two important things to make your model work:<br />
<br />
as I hadn't open your latest models (I'm not always by my WS, often on travel or in meetings) I hadn't noticed <br />
<br />
1) your w for the 2*pi*f is already defined by COMSOL, for the dependent variable 3rd dimension "w"<br />
2) more important, you left the time stepping in "free" mode, you should use "strict" or "intermediate" then you will get far more interesting results ;)<br />
<br />
Note the difference in time stepping, you see it in your log file: there are several " - out" without any steps in between, this is not compatible with a sinusoidal BC entry, te default free is for time decaying processes such as diffusion<br />
<br />
Rule: if any BC is peridic & you are in time stepping => always use strict or intermediate (never "free")<br />
<br />
--<br />
Good luck<br />
IvarThu, 26 Jan 2012 18:59:33 +000024175.1327604373.69201Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69159
Hi<br />
<br />
you are right that first of all you cannot give a time pulse in harmonic sweep solver, it's a constant driving sinus excitation of a given frequency and of a to be defined amplitude and phase (the latter in option if required)<br />
<br />
some physics do not have the time derived equation written up, these cannot be solved in "time". One way to better understand what is happening, is to look at the equations COMSOL uses when you switch the different solver, set up a simple physics model and add a stationary, time and frequency harmonic solver, then look at the physics equation while you switch the solver tab, in stationary there is no d/dt^2 present (we are in a quasi static equation mode, with time dependent they appear, and in harmonic they are replaced by the "omega^2" (it's similar to the Fourier like "s" or "z" functions in data analysis)<br />
<br />
Now if you mix some ACDC and mechanical physics you might have a time dependent for one and none for the other, the last version of COMSOL have some new frequency - transient and frequency stationary solver set up, these evaluate a period (no pulse variation) in temporal mode and then apply the rms values to the time decay or transient analysis. Mostly ACDC EM waves and physics are "instantaneous" w.r.t. mechanical vibrations frequencies, or thermal diffusion, (not always fully true in the MEMS dimensions !) so as they have completely different time scales, the ACDC effect is a mean RMS response applied to the slower mechanics or HT.<br />
<br />
But one can always imagine a case where you do not find (yet) the full combination desired in the Comsol solvers, in these cases one must decide for the closes, and apply some "good engineering practice" to your results<br />
<br />
--<br />
Good luck<br />
IvarThu, 26 Jan 2012 01:02:45 +000024175.1327539765.69159Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p69155
<div class="quote"><br />
Hi<br />
<br />
take a look at the main "equations" of your physics when you choose a time solver and when you have a frequency solver (and compare to a steady state solver).<br />
<br />
Time solver adds the time derivatives to the equations, (compared to a stationary case)<br />
The frequency sweep solvers replaces the time derivatives by a "omega" (and adds a phase term on the sources)<br />
This represents a harmonic development of the solution.<br />
<br />
Consider your FEM with the parallels you probably know better: as when you need to record an analogue signal with a digital DAC: the Nyquist criteria tells you that you need a minimum number of samples (>2) to get a given frequency through your analysis, this corresponds to the FEM meshing and meshing density w.r.t. the dependent variables, and often the gradients of the dependent variables that need to be solved<br />
<br />
In the same way when you do digital data analysis like Fourier transforms you work in the time domain or in the frequency domain and you use the "s" parameters to solve your filtering equations or feedback loop calculation, with harmonic developments. This is a closes parallel I can give you for time solving, and frequency sweep harmonic solving (take car it's "analogies" one cannot apply everything 1:1 ;)<br />
<br />
So for me mixing both solvers for a simple physics like that seems strange and not really physical<br />
<br />
--<br />
Good luck<br />
Ivar<br />
</div><br />
<br />
Hi Ivar,<br />
<br />
This is very elementary but it still is bothering me.<br />
I used a transient physics and time dependent study. <br />
I gave a sine pulse as the input of the form sin(2*pi*1000000[Hz]*t)*(t<=0.000001).<br />
This works fine for me.<br />
<br />
Now i want to use another physics which is not transient so therefore im in the frequency domain right?<br />
My first question is how can i give a pulse as an input ?<br />
In the transient case t was my variable ie time, now i dont know how i can make it work?<br />
By giving the same input as above i dont get the same result as in the transient case.<br />
<br />
My second question is a more general question could seem very elementary but have to ask.<br />
If i use frequency domain as a study with giving different frequencies in the frequency sweep.<br />
Are these frequency values going to be used as the frequency of the sine excitation at the application of the excitation with user defined amplitude? Or are the results checked for the steady state response to the frequency values .<br />
<br />
I don know if ive made it completely clear.<br />
I'm totally confused.<br />
Any kind of reply will help me a lot !<br />
<br />
Thanks a lot for your time and effort!!<br />
<br />
Regards<br />
Glenston<br />
Thu, 26 Jan 2012 00:37:18 +000024175.1327538238.69155Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p68711
Hi<br />
<br />
take a look at the main "equations" of your physics when you choose a time solver and when you have a frequency solver (and compare to a steady state solver).<br />
<br />
Time solver adds the time derivatives to the equations, (compared to a stationary case)<br />
The frequency sweep solvers replaces the time derivatives by a "omega" (and adds a phase term on the sources)<br />
This represents a harmonic development of the solution.<br />
<br />
Consider your FEM with the parallels you probably know better: as when you need to record an analogue signal with a digital DAC: the Nyquist criteria tells you that you need a minimum number of samples (>2) to get a given frequency through your analysis, this corresponds to the FEM meshing and meshing density w.r.t. the dependent variables, and often the gradients of the dependent variables that need to be solved<br />
<br />
In the same way when you do digital data analysis like Fourier transforms you work in the time domain or in the frequency domain and you use the "s" parameters to solve your filtering equations or feedback loop calculation, with harmonic developments. This is a closes parallel I can give you for time solving, and frequency sweep harmonic solving (take car it's "analogies" one cannot apply everything 1:1 ;)<br />
<br />
So for me mixing both solvers for a simple physics like that seems strange and not really physical<br />
<br />
--<br />
Good luck<br />
IvarThu, 19 Jan 2012 11:20:24 +000024175.1326972024.68711Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p68693
Hi Ivar,<br />
<br />
Thanks a lot for your reply.<br />
Your comments are always valuable. <br />
Is it possible to run a frequency sweep using structural mechanics physics and time dependent study?<br />
<br />
Regards<br />
GlenstonThu, 19 Jan 2012 02:31:06 +000024175.1326940266.68693Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p68691
<div class="quote"><br />
Hi Glenston,<br />
I tried opening your model, but I couldn't, So I will simply tell you, how I was able to generate lamb waves using point-force model.<br />
<br />
I used Solid Mechanics module on a 2D cross-sectional aluminum plate, with small thickness, and considerably large width.<br />
<br />
The only boundary condition, I applied is the edge fixed constraint BC, on the edge opposite to the application of the force (that is opposite to the direction of wave propagation)<br />
<br />
The mesh size is considerably small (smaller the better) with minimum element size = 1/10th of lambda.<br />
<br />
The study used is time-dependent, with period = 1/freq and time step = 1/10th of period<br />
<br />
Results must be checked for solid.mises (Von-Mises Stress)<br />
Simulated result should give two modes (S0 and A0) where S0 precedes A0 in space.<br />
Note: S0 will not be as visible as A0, you will have to play with color range<br />
<br />
For further information, please refer to the following paper:<br />
<a href="http://www.ece.cmu.edu/%7Edwg/research/Waves25rev.pdf" title="www.ece.cmu.edu/%7Edwg/research/Waves25rev.pdf">www.ece.cmu.edu/%7Edwg/research/Waves25rev.pdf</a><br />
<br />
Best,<br />
Rohan<br />
</div><br />
<br />
Hi Rohan,<br />
<br />
I went through the paper that you listed in detail and then tried modeling it exactly according to the <br />
specifications in the paper. The metal i have taken is Aluminium with specs that of the paper. <br />
E=70GPa, den=2.7gm/cm^3, v=0.33, thickness=1.59mm and applied force 3.2mm from the y axis.<br />
<br />
I've used Structural Mechanics as the Physics and time dependent study.<br />
A point load is used to give the point force as F*sin(wt)*(sin(wt/10))^2 *(t<10*pi/w).<br />
A fixed constraint point as a BC to the origin. The remaining geometry i have used free boundary condition.<br />
<br />
I ran this for 100kHz frequency.<br />
<br />
My question is it possible to give a frequency sweep in a time dependent study?<br />
If no, is how can u simulate the graphs for the x component of velocity and y component of velocity as a function of frequency?<br />
<br />
Ive not quite got the simulations results that 100kHz should be showing compared to the what the paper shows.<br />
Ive attached the model that ive created in 4.2.<br />
Also attached is the pic of the geometry that i have created.<br />
I would be very helpful if you could upload the model that u created.<br />
<br />
Thanks a lot for your time and effort.<br />
<br />
Regards<br />
Glenston<br />
<br />
Thu, 19 Jan 2012 02:27:14 +000024175.1326940034.68691Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p68044
Hi Glenston,<br />
I tried opening your model, but I couldn't, So I will simply tell you, how I was able to generate lamb waves using point-force model.<br />
<br />
I used Solid Mechanics module on a 2D cross-sectional aluminum plate, with small thickness, and considerably large width.<br />
<br />
The only boundary condition, I applied is the edge fixed constraint BC, on the edge opposite to the application of the force (that is opposite to the direction of wave propagation)<br />
<br />
The mesh size is considerably small (smaller the better) with minimum element size = 1/10th of lambda.<br />
<br />
The study used is time-dependent, with period = 1/freq and time step = 1/10th of period<br />
<br />
Results must be checked for solid.mises (Von-Mises Stress)<br />
Simulated result should give two modes (S0 and A0) where S0 precedes A0 in space.<br />
Note: S0 will not be as visible as A0, you will have to play with color range<br />
<br />
For further information, please refer to the following paper:<br />
<a href="http://www.ece.cmu.edu/%7Edwg/research/Waves25rev.pdf" title="www.ece.cmu.edu/%7Edwg/research/Waves25rev.pdf">www.ece.cmu.edu/%7Edwg/research/Waves25rev.pdf</a><br />
<br />
Best,<br />
RohanTue, 10 Jan 2012 17:35:39 +000024175.1326216939.68044Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p68039
Hi<br />
<br />
I'm not equipped with an acoustic module so I cannot open your file. <br />
In general: using point constraints is rather nasty (stress concentration effects) but if you are not interested in the stress, and you can ignore what is happening in the vicinity of your point load, then the rest should be rather correct. <br />
<br />
You have the choice:i.e. use a harmonic solver (harmonic development of the time based equations) that is much simpler (computational wise) than a sinus time excitation.<br />
<br />
The best is to check the physics equations of COMSOL and that of the Lamb wave theory and see how they match, this to help you select the correct physics<br />
<br />
--<br />
Good luck<br />
IvarTue, 10 Jan 2012 17:17:27 +000024175.1326215847.68039Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p68033
<div class="quote"><br />
Hi<br />
<br />
I do not have 4.0 and I'm not sure my 4.2a is reading correctly the original model, but I have a few comments:<br />
<br />
Solid: <br />
your PZT material isn't really linear isotropic, you must add a 2nd node of linear Elastic material and anisotropic for the PZT material (but this is not enough, see lower)<br />
your device is floating in space, no fixed nodes, you could fix at least one point ans apply a symmetry to the left boundary or "fixed" (careful this would sort your modes to symmetric ones, one need to rerun with antisymmetry BC's, a 2nd time to catch all)<br />
<br />
ES:<br />
there is no need to include the full model in ES the aluminium can be dropped, no ?<br />
<br />
You are in 2D so you are modelling "per meter depth", be aware (nothing wrong with that)<br />
<br />
The main point is: there is no PZT link effects (ES V => PZT forces => solid strain => ES V changes) in solid+ ES you have to write them out yourself<br />
You should rather use PZD physics to have this made precooked by COMSOL (not sure which modules adds PZD though, MEMS certainly, check the COMSOL matrix)<br />
<br />
--<br />
Good luck<br />
Ivar<br />
</div><br />
<br />
Hi Ivar,<br />
<br />
Im trying to do that same thing a little different.<br />
I'm exciting a point source using a single sine pulse and thereafter checking for <br />
the lamb waves in the steel surface.<br />
<br />
Am i applying the right boundary conditions.<br />
Is my physics right?<br />
Acoustic Structure should be fine i'm guessing.<br />
<br />
I've attached a mph file for your further review.<br />
Could you please have a look at it as i have been struggling with it for weeks now!!<br />
<br />
Thanks a lot for your time and effort!<br />
<br />
Regards<br />
Glenston<br />
<br />
Tue, 10 Jan 2012 16:35:26 +000024175.1326213326.68033Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p65001
HI Ivar,<br />
According to your suggestion, I switched to Electromechanics physics module which is added new to the 4.2a version and redesigned the model considering all the BCs as suggested. Hopefully there were no errors, Can you take a look and see If I am going right?<br />
Also can you tell me how to plot 1D plot from the surface plot, thank you.<br />
Wed, 23 Nov 2011 17:21:09 +000024175.1322068869.65001Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p64999
HI Ivar,<br />
According to your suggestion, I switched to 'Electromechanics' physics module which is added new to the 4.2a version and redesigned the model considering all the BCs as suggested. Hopefully there were no errors, Can you take a look and see If I am going right?<br />
Thanks a lotWed, 23 Nov 2011 17:19:42 +000024175.1322068782.64999Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p64997
Yes It is possible, but I am not sure how, and which mode, I should be gettting by my approach.Wed, 23 Nov 2011 17:16:34 +000024175.1322068594.64997Re: Problems generating lamb waves in an Aluminum plate
http://www.comsol.com/community/forums/general/thread/24175/#p64965
<br />
<br />
Is this something like what you are trying to get?? Thoseare antisymmetrical lamb modes if im right.<br />
<br />
<a href="http://www.iaeng.org/publication/WCECS2011/WCECS2011_pp166-171.pdf" title="www.iaeng.org/publication/WCECS2011/WCECS2011_pp166-171.pdf">www.iaeng.org/publication/WCECS2011/WCECS2011_pp166-171.pdf</a><br />
<br />
Guessing from that, lamb modes are possible in comsol...Wed, 23 Nov 2011 12:08:35 +000024175.1322050115.64965