Ultrasonic piezoelectric transducers used in classical nondestructive testing are producing in general longitudinal vibrations in the MHz range. A simple mechanical model of these transducers would be very useful for wave propagation numerical simulations, avoiding the existing complicated models in which the real components of the transducer are modeled by finite elements. The classical model for longitudinal vibrations is not adequate because the generated longitudinal wave is not dispersive, the velocity being the same at any frequency. We have adopted the Rayleigh-Bishop model, which avoids these limitations, even if it is not converging to the first but to the second exact longitudinal mode in an elastic rod, as obtained from the complicated Pochhammer-Chree equations. Since real transducers have significant vibrations damping, we have introduced a damping term in the Rayleigh-Bishop model, increasing the imaginary part and keeping almost identical real part of the wavenumber. Common transducers produce amplitude modulated signals, completely attenuated after several periods. This can be modeled by two close frequencies, producing a “beat” phenomenon, superposed on the high damping. For this reason, we introduce a two-rod Rayleigh-Bishop model with damping. Agreement with measured normal velocity on the transducer free surface is encouraging for continuation of the research.
1. Introduction
Longitudinal vibrations of short rods are the basic phenomenon involved in most ultrasonic transducers, since shear wave transducers are seldom used. From the classical theory of longitudinal vibrations of rods to the exact dynamic elasticity equation of motion, all deal with undamped motions. However, the conflicting requirements of high amplitude resonance displacements and short signal duration used for feature detection are implying more elaborate mechanical models which include damping of the elastic motion.
A review of fundamental models for longitudinal vibrations of rods is given by Shatalov et al. [1] and Fedotov et al. [2]. The classical approximate theory of longitudinal vibration of rods was developed during the 18th century by J. D’Alembert, D. Bernoulli, L. Euler, and J. Lagrange. The classical theory is generally accepted due to its simplicity and is called a unimodal theory since only the fundamental symmetric guided mode is approximated. It is broadly used for the design of low frequency waveguides such as ultrasonic transducers and mechanical filters. Other unimodal theories are presented, such as the Rayleigh-Bishop theories, with two, three, or four terms in radial and longitudinal displacements fields or the Rayleigh-Love theory neglecting the radial stress. All these unimodal theories are compared against the exact Pochhammer-Chree frequency equation in the mentioned reference.
Many researchers have elaborated transducers models including the piezoelectric component, the backing layer, and casing [3]. The transient response of 2 MHz, 12.7 mm diameter transducer is analyzed using an excitation signal of 2 MHz sine-wave half cycle. The model was implemented using commercial software using the finite element method (FEM) by Medina et al.
Coupled oscillator models of resonant transducers, finite element modeling issues in transducers and arrays were developed by Abboud et al. [4]. Accuracy of FEM analysis is dependent upon the accuracy of the dielectric, piezoelectric, and elastic properties used to represent the constituent materials of the model. The work of Ali [5] on a discrete time model of the impulse response of matched and backed piezoelectric transducers has led to the design of sensitive, short-response transducers. The model is obtained from the analysis of the electrical equivalent circuit of a piezoelectric device.
Guelaz et al. [6] present an ultrasonic piezoceramic transducer model integrated in a global measurement cell, a model dedicated to biological ultrasound characterization. Characterization of different mediums, like water and ethanol, by a comparison method is based on the ultrasonic parameter estimation. The good results obtained for temporal simulation and nonlinear parameter estimation compared to measurements and can be extended to multilayer characterization with an ultrasonic nonlinear imaging probe model.
The model used in [7] by Dragan et al. for the ultrasonic sandwich transducer is obtained by appropriate connecting of corresponding mechanical input ports of the transducer component parts, taking into account the actual existing contacts of the surfaces represented by these inputs. Under the condition that the outer contour surfaces of transducer component parts are loaded by surrounding environment acoustic impedances, which are small for the unloaded surface, and with alternating voltage applied to the electric inputs, it is possible to determine any transfer function of interest for the given system.
The FEM was proven by Nygren [8] to be a useful way of simulating piezoelectric ultrasonic transducers. The method is versatile as it can be used to emulate both one-dimensional models and more realistic three-dimensional transducer models. However, the suggested minimal requirement of 6 elements per wavelength would cause the solution of a complete three-dimensional model to become practically unattainable. Thankfully, many models can be significantly reduced in size, when mirror and axial symmetry is taken into consideration.
The work of Kocbach [9] presents an important basis for the first step of the design process of piezoelectric transducers. This has been done through an analysis of the radiated sound field, response functions, and vibration of the basic parts of piezoelectric disk transducers (the piezoelectric disk and the front layer) and an analysis of the influence of varying geometry and material parameters of the transducer on these quantities. The analysis has also been targeted at understanding the influence of the different design parameters of the simple piezoelectric transducer on the response functions of the transducer, in order to provide guidelines on how to vary the design parameters in the second step of the design process. The FEM has been used to simulate the 3D numerical models developed by Porto [10]. With the numerical optimization of new ultrasonic transducer designs in prospect, the computation time is a key problem. The different existing or developed optimization methods, absolutely valid with the correspondent analytical models, have shown their limitations. An optimization method based on genetic algorithms (GA) has been envisioned, thanks to its ease of programming and its intrinsic modularity. The developed and implemented new algorithms have been tested using many test functions and have shown satisfactory results. Consequently, new 2D models of ultrasonic transducers have been imagined and optimized, with accurate results and fast enough computations.
Schubert and Lamek [11] present the elastodynamic finite integration technique (EFIT) combined with point-source-synthesis (PSS) in order to efficiently calculate 3D ultrasonic wave fields in the time domain. The technique can be applied to arbitrarily layer locally isotropic media and includes the complete spectrum of linear elastic wave phenomena like surface and interface waves, head waves, mode conversion, and multiple scattering. The proposed technique is not restricted to axisymmetric isotropic problems since the basic point source data can also be calculated by using a single 3D calculation of an arbitrarily layered anisotropic medium.
An accurate simulation of the dynamic response of a longitudinal wave transducer requires parametric sweep of the damping characteristics. Damped longitudinal vibrations can be modeled in several ways. Pintelon et al. [12] used frequency dependent complex Young modulus, which can be can be experimentally determined and is capable of measuring low level of damping in the frequency domain. Wave propagation in viscoelastic structures, including rods, has also been investigated by Castaings et al. [13], using finite elements and complex moduli for the materials.
In the present work, a unimodal model for the longitudinal vibrations of a rod is investigated in the second section. The classical model for longitudinal vibrations is proved to be inadequate because the generated longitudinal wave is not dispersive, the velocity being the same at any frequency. The slightly more complicated Rayleigh-Love model is dispersive, but, beyond a critical frequency, the wave will not propagate. We have adopted the Rayleigh-Bishop model, which avoids these limitations, even if it is not converging to the first but to the second exact longitudinal mode in an elastic rod, which can be obtained from the complicated Pochhammer-Chree equations. Since real transducer has significant vibrations damping, we have introduced, in the third section, a damping term in the Rayleigh-Bishop model, proving negligible changes in the real part of the wavenumbers but significant frequency increasing imaginary parts absolute values. The real transducers produce an amplitude modulated signal, completely attenuated after several periods. Such behavior can be modeled by two close valued frequencies, producing a “beat” phenomenon, superposed on high damping. For this reason, we introduce a two-rod Rayleigh-Bishop model with damping in the fourth section. Agreement with measured normal velocity, shown in the fifth section, is encouraging.
2. Comparison of Theoretical Models in the High Frequency Range
The mechanical model considers a circular rod of length L, cross-section area S, mass density ρ, Young modulus E, and Poisson coefficient ν. According to the Rayleigh-Love theory, the longitudinal vibrations of the rod defined in Figure 1 are described by the following displacements field:
(1)u(x,r,t)=u(x,t),w(x,r,t)=-νru,x.
The radial displacement w is proportional to the local radius r(0≤r≤R) in which R is the outer radius of the rod and u,x=∂u(x,t)/∂x. The strains in cylindrical frame, using (1), are as follows:
(2)εxx=u,x;εrr=w,r=-νu,x;εxr=u,r+w,x=-νu,xx;εθθ=wr=-νu,x.
The relevant stress components, using the Lamé constants λ=Eν/(1-2ν)(1+ν),μ=E/2(1+ν), are as follows:
(3)σxx=(λ+2μ)εxx+λ(εrr+εθθ)=Eu,x,σxr=μεxr=-νμru,xx.
It can be proven that the radial displacements defined cancel the radial stress:
(4)σrr=(λ+2μ)εrr+λ(εxx+eθθ)=-ν(λ+2μ)u,x+λ(u,x+wr)=-ν(λ+2μ)u,x+λ(1-ν)u,x=0.
Geometry of the rod with local displacements in cylindrical frame (a). Two rods (b).
2.1. The Rayleigh-Love Model
The kinetic energy (T) and the strain potential energy (V) for this model are as follows:
(5)T=∫0Lρ2∫S(u˙2+ν2r2u˙,x2)dx=∫0Lρ2(Su˙2+ν2Ip2u˙,x2)dxV=12∫0L[∫S(σxxεxx+σrrεrr+σxrεxr)]dx≅12∫0LSEu,x2dx.
The contribution of the shear energy was neglected in the last expression, according to the Rayleigh-Love theory. The polar moment of inertia for a circular cross-section is Ip=πR4/2. The Lagrangian is L=T-V=∫0LΛdx in which the Lagrangian density Λ is
(6)Λ=S2(ρu˙2+ρν2rg2u˙,x2-Eu,x2),
with rg=R/2, which is the radius of gyration for a circular cross-section or radius R. The Lagrangian density must satisfy the Euler-Lagrange differential equation [14]:
(7)∂∂t(∂Λ∂u˙)+∂∂x(∂Λ∂u,x)-∂2∂x∂t(∂Λ∂u˙,x)=0,
which leads to the differential equation of motion:
(8)ρu¨-Eu,xx-ρν2rg2u¨,xx=0.
This theory presented in [1] proved to be more accurate in predicting the variation of the longitudinal wave velocity with the frequency than the classical theory, especially for the high frequency range, in which we are interested. The classical formulation used in any textbook of longitudinal vibrations of cylindrical rods (e.g., [15, 16], etc.), can be recovered by neglecting the last term.
A harmonic plane wave propagating with velocity c, wavenumber k=ω/c, and angular frequency ω=2πf, in which f is the frequency, as
(9)u(x,t)=Aexp[i(kx-ωt)]
can be injected in (8) to get
(10)[-ρω2+k2E-ρν2rg2ω2k2]u(x,t)=0.
The dispersion equation for an elastic rod according to the Rayleigh-Love theory is then
(11)k2(E-ρν2rg2ω2)=ρω2.
The dispersion curve can be plotted in several ways. Considering a rod of steel (E=210 GPa, ν=0.31, and ρ=7.8 g/cm3) of radius R=12.5 mm, the dispersion curves in physical quantities are shown on Figure 2. For comparison, the classical model solution is also plotted, which is a straight line, as follows:
(12)k=ωc0;c0=Eρ,
and the first two exact solutions of the Pochhammer-Chree dispersion equation for an elastic cylinder [17]. Solving the exact solution is however not a very simple task, involving Bessel functions of complex arguments.
Dispersion curve for a rod using the Rayleigh-Love theory, compared to the classical and exact theories.
The Pochhammer-Chree dispersion equation is considerably more complicated than (8) of the Rayleigh-Love model, which is a reason for which simpler models to be pursued by many authors. We have developed a numerical root-finding algorithm for the Pochhammer-Chree dispersion equation, and the obtained wavenumbers of the first two modes are used for comparison against the approximate theories. These are called “exact,” in the following, since they correspond to the “exact” elasticity theory but are obtained by numerical procedures with finite accuracy.
At relatively low frequencies (below 0.1 MHz in this case) the Rayleigh-Love, classical, and first exact modes are coinciding. Between 0.1 and 0.25 MHz, the Rayleigh-Love solution is better than the classical solution, being closer to the exact longitudinal mode, called L(0,1) of the exact solution. However, the Rayleigh-Love model has a major shortcoming for high frequencies: the wave number becomes purely imaginary, if
(13)f>12πEρν2rg2,
which means that f<300 kHz in this case. The Rayleigh-Love model is thus inadequate for longitudinal vibrations in the MHz range.
2.2. The Rayleigh-Bishop Model
The Rayleigh-Bishop model is also a unimodal model, meaning that only one longitudinal wave mode is approximated. The kinetic energy T and the strain potential energy V are as follows:
(14)T=∫0Lρ2∫S(u˙2+ν2r2u˙,x2)dx=∫0Lρ2(Su˙2+ν2Ip2u˙,x2)dx,V=12∫0L[∫S(σxxεxx+σrrεrr+σxrεxr)]dx≅12∫0L[SEu,x2+ν2μIpu,xx2]dx.
The Lagrangian density Λ is
(15)Λ=S2(ρu˙2+ρν2rg2u˙,x2-Eu,x2-ν2μrg2u,xx2),
in which rg=R/2 is the radius of gyration for a circular cross section or radius R.
The Lagrangian density must satisfy [14] the Euler-Lagrange differential equation:
(16)∂∂t(∂Λ∂u˙)+∂∂x(∂Λ∂u,x)-∂2∂x∂t(∂Λ∂u˙,x)-∂2∂x2(∂Λ∂u,xx)=0,
which leads to the differential equation of motion:
(17)ρu¨-Eu,xx-ρν2rg2u¨,xx+ν2μrg2u,xxxx=0.
A harmonic plane wave solution can be also used in this case:
(18)u(x,t)=Aexp[i(kx-ωt)].
Equation (17) becomes
(19)[-ρω2+k2E-ρν2rg2ω2k2+ν2μrg2k4]u(x,t)=0.
The dispersion equation for an elastic rod according to the Rayleigh-Bishop theory is then
(20)k2(E-ρν2rg2ω2+ν2μrg2k2)=ρω2.
The same rod is used in the numerical simulation and the results from our numerical analysis are shown on Figure 3. The model proves to be considerably more accurate than the classical model (dotted line). The dispersion curve of the Rayleigh-Bishop model tends towards the first exact mode of the Pochhammer-Chree theory in case of low frequencies and towards the second mode as the frequencies increase. There is no critical frequency, beyond which there are no solutions, as in the Rayleigh-Love model, so that a wavenumber k can be obtained for any frequency. This model is better adapted to model high frequency longitudinal vibrations, as required for ultrasonic transducers.
Dispersion curve for a rod using the Rayleigh-Bishop theory, compared to the classical and exact theories.
3. Damping of Longitudinal Vibrations
As shown above, the Rayleigh-Bishop unimodal model is the best in approximating the dispersion curves. The vibrations are always subjected to damping. In this paper, we are including damping, which is assumed to be proportional to the strain velocity variation of the selected model. Under this assumption, using a proportionality damping constant β, the differential equation of motion (17) becomes
(21)ρu¨-Eu,xx-ρν2rg2u¨,xx+ν2μrg2u,xxxx-βEu˙,xx=0.
Looking for damped vibrations of rods, one can select the displacement function as follows:
(22)u(x,t)=Aexp(ikx)exp(-iωt),
in which the angular frequency is expected to be complex valued. The dispersion equation becomes
(23)k2(E-ρν2rg2ω2+ν2μrg2k2-iβω)=ρω2.
For a given real wavenumber k, the complex angular frequency can be obtained as ω=ωr+iωi, by solving a simple polynomial equation of third degree, for which numerical and even closed form solutions exist.
The real and imaginary parts of the angular frequency as functions of the real wavenumber, for the same cylindrical rod are shown on Figures 4 and 5, respectively. The influence on the real part of the angular frequency is negligible for the selected frequency range, as can be seen in Figure 4. For values of β up to 104 kg·s-1·m-1, the influence on damping, represented by the imaginary part, is less than 3% of the real part, but, for β=105 kg·s-1·m-1, the imaginary part is increasing up to between 15 and 30% of the real part. Damping is increasing in absolute value with frequency, but, for increasing frequency, the imaginary part tends to a eliminated value depending on β.
Angular frequency (real part) versus frequency for three damping coefficients β.
Angular frequency (imaginary part) versus frequency for three damping coefficients β.
4. Longitudinal Damped Vibrations of a Short Rod
The foreseen application of the longitudinal vibrations is in ultrasonic transducers. An analytic model based on the Rayleigh-Bishop model without damping was presented in [18]. The damped vibrations mechanical model of such a transducer, based on the formulation of the previous section, can be represented as two coaxial rods of identical diameter connected by the piezoelectric crystal, as in [15].
We consider now two connected rods, made of different materials E1, E2, cross sections S1 and S2, and lengths L1 and L2, in order to follow the derivation of equivalent equations in Section 2.3 of [15]. A practical particular case will be deduced for E1=E2 and S1=S2, joined at x=0 (Figure 1(b)).
The damped harmonic motion in the two rods can be expressed as
(24)um(x,t)=[Amcos(kmx)+Bmsin(kmx)]exp(-iωmt);iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiim=1,2,
in which the complex angular frequencies ωm are solutions of (23) for given wavenumbers km. The boundary conditions for the first rod are represented by the longitudinal displacement and stress continuity at x=0 and longitudinal stress cancelling at x=L1:
(25)u1(0,t)=u2(0,t);E1S1u1,x(0,t)=E2S2u2,x(0,t);u1,x(L1,t)=0.
For the second rod, the longitudinal stress cancels at the free end
(26)u2,x(-L2,t)=0.
The exact boundary conditions for the free end in the Rayleigh-Bishop case can be found in [1]. Here, the shear stress is not taken into account. These four equations can be written as
(27)A1exp(-iω1t)=A2exp(-iω2t)E1S1k1B1exp(-iω1t)=E2S2k2B2exp(-iω2t)k1[-A1sin(k1L1)+B1cos(k1L1)]=0k2[A2sin(k1L1)+B2cos(-k2L2)]=0.
It follows that
(28)A1A2=B2B1k2E2S2k1E1S1;B1A1=tan(k1L1);B2A2=-tan(k2L2)
or
(29)k2E2S2k1E1S1sin(k2L2)cos(k1L1)+sin(k1L1)cos(k2L2)=0.
This equation resembles (2.3.16) from [15], except for the ratio of wavenumbers, which is missing by mistake in the cited reference. In the particular case of identical materials and cross sections, a relation between the two wavenumbers can be thus obtained. In practical cases of ultrasonic transducers, the first antisymmetric mode of vibration, with respect to x=0, is used. This means that the condition cos(k1L1)=0 will provide the wavenumber for the first rod and, from (29), it follows that cos(k2L2)=0 provides the smallest wavenumber for the second rod. Consequently,
(30)k1=π2L1;k2=π2L2.
If the two rods have different lengths, according to (23), the assembly will have two resonance frequencies and two damping ratios. For MHz range frequencies, the cylinders lengths are less than 1 mm for the considered steel.
5. Experimental Validation
A 1 MHz central frequency longitudinal wave ultrasonic transducer has received a short pulse of 200 V from the pulser-receiver. The free response of the transducer was measured by a laser velocimeter and recorded on a digital oscilloscope (Figure 6). The simplest model of the transducer can be imagined as two connected short steel rods, according to the previous section. The material choice is made in order to minimize acoustic impedance mismatch in case of nondestructive testing of steel materials.
Experimental setup used to measure free response of the ultrasonic transducer.
By iterative testing of rod lengths and attenuations, the following parameters have been obtained for a best fit of the recorded normal velocity of the transducer: L1=0.84 mm and L2=0.64 mm and attenuation: Imag (ω)=-6e5 rad/s.
The agreement (Figure 7) is satisfactory for most applications, proving the adequacy of the developed simple and efficient mechanical model.
Computed data superposed on the measured normal velocity on transducer surface.
6. Conclusions
Longitudinal vibrations of rods at high frequencies, in the MHz range, have been investigated using the classical and the Rayleigh-Love models, which both proved inadequate since the dispersion curves are too far from the exact values or simply constants in this frequency range. The Rayleigh-Bishop model was adopted and we have introduced a damping term in this model. As a consequence of this improved model, wavenumbers are almost unaffected by damping. On the other hand, the wave damping is increasing with frequency in an asymptotic manner, which can be accepted in a simple model.
A free vibrating ultrasonic transducer was modeled as two unequal rods acted by negligible thin piezoelectric element between them. The purpose is to mimic the behavior of the real transducer by the simplest mechanical model. By an iterative procedure, the optimal choice of parameters has led to an accurate representation of measured normal velocity of a 1 MHz longitudinal ultrasonic transducer.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The present work was supported by the Romanian UEFISCDI Exploratory Research Project PN-II-ID-PCE-2011-3-0512.
ShatalovM.MaraisJ.FedotovI.TenkamM. D.SchmidtM.M. SchmidtLongitudinal vibration of isotropic solid rods: from classical to modern theories2011InTech187214FedotovI.MaraisJ.ShatalovM.TenkamH. M.Hyperbolic models arising in the theory of longitudinal vibration of elastic bars201172118MR2794509MedinaJ. E. M.BuiochiF.AdamowskiJ.Numerical modeling of a circular piezoelectric ultrasonic transducer radiating in water2Proceedings of the ABCM Symposium Series in Mechatronics2006458464AbboudN. N.WojcikG. L.VaughanD. K.MouldJ.PowellD. J.NikodymL.Finite element modeling for ultrasonic transducers3341International Symposium on Medical Imaging, Ultrasonic Transducer Engineering ConferenceFebruary 1998San Diego, Calif, USA1942Proceedings of SPIE10.1117/12.3080152-s2.0-0032401685AliM.Analysis of broadband piezoelectric transducers by discrete time model2000232287295GuelazR.KourticheD.NadiM.HerveY.Ultrasonic piezoceramic transducer modeling with VHDL-AMS IEEE 1076.1Proceedings of the 3rd IEEE Conference on SensorsOctober 2004Vienna, Austria87902-s2.0-2794448090810.1109/ICSENS.2004.1426106DraganM.MilanR.ZoranP.GoranS.Influence of ultrasonic transducer acoustic impedances and dimensions on its input electrical impedance2008515972NygrenM. W.2011Trondheim, NorwayNorwegian University of Science and TechnologyKocbachJ.2000Bergen, NorwayDepartment of Physics, University of BergenPortoD.2009Lausanne, SwitzerlandÉcole Polytechnique Fédérale de LausanneSchubertF.LamekB.3-D ultrasonic transducer modeling using the elastodynamic finite integration technique in combination with point-source-synthesisProceedings of 4th International Workshop NDT in ProgressNovember 2007Prague, Czech RepublicPintelonR.GuillaumeP.VanlanduitS.de BelderK.RolainY.Identification of Young's modulus from broadband modal analysis experiments200418469972610.1016/S0888-3270(03)00045-12-s2.0-2442459945CastaingsM.BaconC.HostenB.PredoiM. V.Finite element predictions for the dynamic response of thermo-viscoelasicity material structures200411531125113310.1121/1.16393322-s2.0-1542376817CourantR.HilbertD.1953New York, NY, USAJohn Wiley & SonsGraffK. F.1975New York, NY, USADoverde SilvaC. W.2000Boca Raton, Fla, USACRC Press LLCAchenbachJ. D.1973Amsterdam, The NetherlandsNorth HollandTenkamH.A theory of piezoelectric thickness vibrator transducer based on the Rayleigh and Bishop model1Proceedings of the World Congress on Engineering (WCE '12)July 2012London, UK48