Printer Friendly

High frequency longitudinal damped vibrations of a cylindrical ultrasonic transducer.

I. 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 [rho], Young modulus E, and Poisson coefficient v. According to the Rayleigh-Love theory, the longitudinal vibrations of the rod defined in Figure 1 are described by the following displacements field:

u (x, r,t) = u (x, t), w (x, r, t) = -[vru,.sub.x]. (1)

The radial displacement w is proportional to the local radius r (0 [less than or equal to] r [less than or equal to] R) in which R is the outer radius of the rod and [u,.sub.x] = [partial derivative]u(x,t)/[partial derivative]x. The strains in cylindrical frame, using (1), are as follows:

[[epsilon].sub.xx] = [u,.sub.x]; [[epsilon].sub.rr] = [w,.sub.r] = -[vu,.sub.x]; [[epsilon].sub.xr] = [u,.sub.r] + [w,.sub.x] = -[vu,.sub.xx]; [[epsilon].sub.[theta][theta]] = w/r =-[vu,.sub.x]. (2)

The relevant stress components, using the Lame constants [lambda] = Ev/(1 - 2v)(1 + v), [mu] = E/2(1 + v), are as follows:

[[sigma].sub.xx] = ([lambda] + 2[mu])[[epsilon].sub.xx] + [lambda]([[epsilon].sub.rr] + [[epsilon].sub.[theta][theta]) = [Eu,.sub.x], [[sigma].sub.xr] = [mu][[epsilon].sub.xr] = -[v[mu]ru,.sub.xx]. (3)

It can be proven that the radial displacements defined cancel the radial stress:


2.1. The Rayleigh-Love Model. The kinetic energy (T) and the strain potential energy (V) for this model are as follows:


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 [I.sub.p] = [pi][R.sup.4]/2.

The Lagrangian is L = T - V = [[integral].sup.L.sub.0] [LAMBDA] dx in which the Lagrangian density [LAMBDA] is


with [r.sub.g] = R/[square root of 2], which is the radius of gyration for a circular cross-section or radius R. The Lagrangian densitymust satisfy the Euler-Lagrange differential equation [14]:


which leads to the differential equation of motion:


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 = [omega]/c, and angular frequency [omega] = 2[pi]f, in which f is the frequency, as

u(x,t) = A exp [i(kx - [omega]t)] (9)

can be injected in (8) to get

[-[rho][[omega].sup.2] + [k.sup.2]E - [rho][v.sup.2][r.sup.2.sub.g][[omega].sup.2][k.sup.2]] u (x, t) = 0. (10)

The dispersion equation for an elastic rod according to the Rayleigh-Love theory is then

[k.sup.2] (E - [rho][v.sup.2][r.sup.2.sub.g][[omega].sup.2]) = [rho][[omega].sup.2]. (11)

The dispersion curve can be plotted in several ways. Considering a rod of steel (E = 210 GPa, v = 0.31, and [rho] = 7.8 g/[cm.sup.3]) 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:

k = [omega]/[c.sub.0]; [c.sub.0] = [square root of E/[rho], (12)

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.

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

f > 1/2[pi] [square root of (e/rho][v.sup.2][r.sup.2.sub.g]], (13)

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:


The Lagrangian density [LAMBDA] is


in which [r.sub.g] = R/[square root of 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:


which leads to the differential equation of motion:

[[rho][??] - [Eu,.sub.xx] - [rho][v.sup.2][r.sup.2.sub.g][[??].sub.xx] + [v.sup.2][mu][r.sup.2.sub.g][u,.sub.xxxx] = 0 (17)

A harmonic plane wave solution can be also used in this case:

u(x,t) = A exp [i(kx-[omega]t)]. (18)

Equation (17) becomes

[-[rho][[omega].sup.2] + [k.sup.2]E - [rho][v.sup.2][r.sup.2.sub.g][[omega].sup.2][k.sup.2] + [v.sup.2][mu][r.sup.2.sub.g][k.sup.4]] u (x, t) = 0. (19)

The dispersion equation for an elastic rod according to the Rayleigh-Bishop theory is then

[k.sup.2] (E - [rho][v.sup.2][r.sup.2.sub.g][[omega].sup.2] + [v.sup.2][mu][r.sup.2][k.sup.2]) = [rho][[omega].sup.2]. (20)

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 theoryin 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.

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 [beta], the differential equation of motion (17) becomes


Looking for damped vibrations of rods, one can select the displacement function as follows:

u (x, t) = A exp (ikx) exp (-i[omega]t), (22)

in which the angular frequency is expected to be complex valued. The dispersion equation becomes

[k.sup.2](E - [rho][v.sup.2][r.sup.2.sub.g][[omega].sup.2] + [v.sup.2][mu][r.sup.2.sub.g] [k.sup.2] - i[beta][omega]) = [rho][[omega].sup.2]. (23)

For a given real wavenumber k, the complex angular frequency can be obtained as [omega] = [[omega].sub.r] + i[[omega].sub.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 [beta] up to [10.sup.4] kg x [s.sup.-1] x [m.sup.-1], the influence on damping, represented by the imaginary part, is less than 3% of the real part, but, for [beta] = [10.sup.5] kg x [s.sup.-1] x [m.sup.-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 [beta].

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 [E.sub.1], [E.sub.2], cross sections [S.sub.1] and [S.sub.2], and lengths [L.sub.1] and [L.sub.2], in order to follow the derivation of equivalent equations in Section 2.3 of [15]. A practical particular case will be deduced for [E.sub.1] = [E.sub.2] and [S.sub.1] = [S.sub.2], joined at x = 0 (Figure 1(b)).

The damped harmonic motion in the two rods can be expressed as

[u.sub.m] (x, t) = [[A.sub.m] cos ([k.sub.m]x) + [B.sub.m] sin ([k.sub.m]x)] exp (-i[[omega].sub.m]t); m = 1,2, (24)

in which the complex angular frequencies [[omega].sub.m] are solutions of (23) for given wavenumbers [k.sub.m]. 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 = [L.sub.1]:

[u.sub.1] (0,t) = [u.sub.2] (0,t); [E.sub.1][S.sub.1] [u.sub.1,x] (0,t) = [E.sub.2][S.sub.2][u.sub.2,x] (0,t); [u.sub.1,x] ([L.sub.1], t) = 0. (25)

For the second rod, the longitudinal stress cancels at the free end

[u.sub.2,x] (-[L.sub.2],t)=0. (26)

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


It follows that



[k.sub.2][E.sub.2][S.sub.2]/[k.sub.1][E.sub.1][S.sub.1]sin([k.sub.2][L.sub.2])cos ([k.sub.1][L.sub.1])cos([k.sub.1][L.sub.1]) + sin ([k.sub.1][L.sub.1])cos ([k.sub.2][L.sub.2]) = 0. (29)

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([k.sub.1][L.sub.1]) = 0 will provide the wavenumber for the first rod and, from (29), it follows that cos([k.sub.2][L.sub.2]) = 0 provides the smallest wavenumber for the second rod. Consequently,

[k.sub.1] = [pi]/2[L.sub.1]; [k.sub.2] = [pi]/2[L.sub.2]. (30)

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 1MHz 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.

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: [L.sub.1] = 0.84 mm and [L.sub.2] = 0.64 mm and attenuation: Imag ([omega]) = -6e5 rad/s.

The agreement (Figure 7) is satisfactory for most applications, proving the adequacy of the developed simple and efficient mechanical model.

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.


The present work was supported by the Romanian UEFISCDI Exploratory Research Project PN-II-ID-PCE-2011-3-0512.


[1] M. Shatalov, J. Marais, I. Fedotov, M. D. Tenkam, and M. Schmidt, "Longitudinal vibration of isotropic solid rods: from classical to modern theories," in Advances in Computer Science and Engineering, M. Schmidt, Ed., pp. 187-214, InTech, 2011.

[2] I. Fedotov, J. Marais, M. Shatalov, and H. M. Tenkam, "Hyperbolic models arising in the theory of longitudinal vibration of elastic bars," The Australian Journal of Mathematical Analysis and Applications, vol. 7, no. 2, pp. 1-18, 2011.

[3] J. E. M. Medina, F. Buiochi, and J. Adamowski, "Numerical modeling of a circular piezoelectric ultrasonic transducer radiating in water," in Proceedings of the ABCM Symposium Series in Mechatronics, vol. 2, pp. 458-464, 2006.

[4] N. N. Abboud, G. L. Wojcik, D. K. Vaughan, J. Mould, D. J. Powell, and L. Nikodym, "Finite element modeling for ultrasonic transducers," in International Symposium on Medical Imaging, Ultrasonic Transducer Engineering Conference, vol. 3341 of Proceedings of SPIE, pp. 19-42, San Diego, Calif, USA, February 1998.

[5] M. Ali, "Analysis of broadband piezoelectric transducers by discrete time model," Egyptian Journal of Solids, vol. 23, no. 2, pp. 287-295, 2000.

[6] R. Guelaz, D. Kourtiche, M. Nadi, and Y. Herve, "Ultrasonic piezoceramic transducer modeling with VHDL-AMS IEEE 1076.1," in Proceedings of the 3rd IEEE Conference on Sensors, pp. 87-90, Vienna, Austria, October 2004.

[7] M. Dragan, R. Milan, P. Zoran, and S. Goran, "Influence of ultrasonic transducer acoustic impedances and dimensions on its input electrical impedance," Working and Living Environmental Protection, vol. 5, no. 1, pp. 59-72, 2008.

[8] M. W. Nygren, Finite Element Modeling of Piezoelectric Ultrasonic Transducers, Norwegian University of Science and Technology, Trondheim, Norway, 2011.

[9] J. Kocbach, Finite Element Modeling of Ultrasonic Piezoelectric Transducers, Department of Physics, University of Bergen, Bergen, Norway, 2000.

[10] D. Porto, Design Methodology and Numerical Optimization of Ultrasonic Transducers for Spinal Surgery, Ecole Polytechnique Federale de Lausanne, Lausanne, Switzerland, 2009.

[11] F. Schubert and B. Lamek, "3-D ultrasonic transducer modeling using the elastodynamic finite integration technique in combination with point-source-synthesis," in Proceedings of 4th International Workshop NDT in Progress, Prague, Czech Republic, November 2007.

[12] R. Pintelon, P. Guillaume, S. Vanlanduit, K. de Belder, and Y. Rolain, "Identification of Young's modulus from broadband modal analysis experiments," Mechanical Systems and Signal Processing, vol. 18, no. 4, pp. 699-726, 2004.

[13] M. Castaings, C. Bacon, B. Hosten, and M. V Predoi, "Finite element predictions for the dynamic response of thermoviscoelasicity material structures," Journal of the Acoustical Society of America, vol. 115, no. 3, pp. 1125-1133, 2004.

[14] R. Courant and D. Hilbert, Methods of Mathematical Physics, John Wiley & Sons, New York, NY, USA, 1953.

[15] K. F. Graff, Wave Motion in Elastic Solids, Dover, NewYork, NY, USA, 1975.

[16] C. W. de Silva, Vibration. Fundamentals and Practice, CRC Press LLC, Boca Raton, Fla, USA, 2000.

[17] J. D. Achenbach, Wave Propagation in Elastic Solids, North Holland, Amsterdam, The Netherlands, 1973.

[18] H. Tenkam, "A theory of piezoelectric thickness vibrator transducer based on the Rayleigh and Bishop model," in Proceedings of the World Congress on Engineering (WCE '12), vol. 1, p. 48, London, UK, July 2012.

Mihai Valentin Predoi, (1) Cristian Catalin Petre, (2) Ovidiu Vasile, (1) and Mihail Boiangiu (1)

(1) Department of Mechanics, Polytechnic University of Bucharest, Splaiul Independentei 313, 60048 Bucharest, Romania

(2) Department of Strength of Materials, University Politehnica of Bucharest, Splaiul Independent 313, 60048 Bucharest, Romania

Correspondence should be addressed to Mihai Valentin Predoi;

Received 27 February 2014; Accepted 16 June 2014; Published 24 June 2014

Academic Editor: Xinjie Zhang
COPYRIGHT 2014 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Predoi, Mihai Valentin; Petre, Cristian Catalin; Vasile, Ovidiu; Boiangiu, Mihail
Publication:Shock and Vibration
Article Type:Report
Date:Jan 1, 2014
Previous Article:A stochastic wavelet finite element method for 1D and 2D structures analysis.
Next Article:Feedback linearisation for nonlinear vibration problems.

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |