 Open Access
 Total Downloads : 87
 Authors : Nguyen Van Thinh , Hoang Thi My Le
 Paper ID : IJERTV8IS060150
 Volume & Issue : Volume 08, Issue 06 (June 2019)
 Published (First Online): 10062019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Researching to Determine the Characteristic Parameters of the Power Ultrasonic Transducer by Finite Element Method and COMSOLMULTIPHYSICS Program
Nguyen Van Thinh
The University of Technology and Education The University of Danang
Danang, Vietnam
Hoang Thi My Le
The University of Technology and Education The University of Danang
Danang, Vietnam
Abstract The power ultrasonic transdurce operates in range of frequencies from 18 KHz to 45 kHz [9]. Depending on the material, the tranducer structure has the characteristic parameters that suitables for each application. The research results in this paper: building the power ultrasound transducer model from the PZT hard piezoelectric materials and using finite element method (FEM – Finite Element Method) and COMSOL Multiphysics (CM) prgram determine characteristic parameters. The research result has identified resonant oscillation frequency, the transducer displacement and the characteristics of the sound pressure emanating from the ultrasound transducer. These results are fundamental for the designing and applying in the high power ultrasound techniques..

Continuous equations in piezoelectric environment
{T} = [cE]{S} – [e]{E} (1)
{D} = [e]{S} + [S]{E} (2)
In which, {T}, {D}, {E} are the stress vector, the electric displacement vector and the electric field intensity vector respectively. [e], [S], [cE] are the piezoelectric coefficient matrix, the dielectric constant matrix and the stiffness matrix respectively. Index symbolizes for matrix or shift vector. Index "" symbolizes for matrix or shift vector. The mechanical terms in equation (1) follow 2 Newton's law, according to (3).
2{ } = {}
KeywordsLangevin trÃ¡nducer; PZT hard piezoelectric
(3)
ceramic ; sound pressure field; FEM; COMSOLMultiphysics
With = { u v w } is mechanical displacement vector
along axes x , y and z. The {S} strain vector relates the

INTRODUCTION
The Langevin ultrasound transducer is a key element using to make high power ultrasonic devices. Depending on the material and the design model, the ultrasonic transducer has different characteristic parameters. In order to achieve the
{R} the displacement vector according to the expression according to the expression:
{S} = [B ]{ } (4)
In which:
specific application purpose, it is necessary to determine the ultrasonic transducer characteristic parameters. By
0 0 0
0
0
constructing PZT doped hard piezoelectric ceramic models and Langevin transducer, using finite element method in the
[B ] =0
0
(5)
0
0
PZT hard piezoelectric material environment, COMSOL
Multiphysics program survey the dependence of characteristic
[ 0 0]
parameters according to frequency. In this papper, we simulated to determinate the oscillation characteristic, the resonance frequency of the PZT hard piezoelectric ceramic, the Langevin transducer displacement and the sound pressure field to emanate rom the front of the Langevin transducer. Research results are the basis for selecting materials, designing plans and applications.

FINITE ELEMENT METHOD IN PIEZOELECTRIC ENVIRONMENT [4], [6], [9]
The electric terms in (1) follow Gauss's law, with the assumption that piezoelectric materials are insulating materials, no charge current running in the transducer.
.{D} = 0 (6)
The E electric field is related to the voltage by
{E} = {B} (7)
Finite element method is applied to analyse the characteristic piezoelectric transducer parameters with
with, {B} = {
} (8)
different design models. In order to achieve accurate analytical results, it is necessary to set up mathematical and physical equations and impose boundary condition on each case.
In which, {G} = {T D}, {G} = {SE} are the stress vector and the general deformation vector respectively. The [c]
general elastic matrix form:
[c] = [ ](9)

SURVEYING THE LANGEVIN TRANSDUCER CHARACTERISTIC PARAMETERS BY COMSOL –
MUTIPHYSICS PROGRAM
The stress and strain relationships are represented by the equation:
{G} = [BG]{ G} (10)
In which, { G}= { } = { u v w } is the general displacement vector and the general matrix with form:
[BG] = [ 0 ] (11)0


Dynamic equations in piezoelectric environment
With piezoelectric materials that have satisfied the continuous equations, the G general displacement of any approximated point by [ ] displacement function and the {} vector containing the finite displacement nodes.
{ G } = [ ]{} (12)
i vi bÃ i toÃ¡n nÃ y, ti nÃºt th h cÃ³ 4 bc t do. Trong h ta Descartes 3 chiu vector ca nÃºt th h cÃ³ dng.
For this problem, at the h node has the 4 freedom degrees. In the 3dimensional Descartes coordinate system, the
{h}vector of the hth node is the form:
{h} = { h h} = { uh vh wh } (13)
Mathematical, physical basis for determining piezoelectric material parameters are described by general equations (1)Ã·(13).
An arbitrary problem is solved by imposing boundary conditions, algorithm of finite element analysis based on input data as matrix of material coefficients: elasticity coefficient, piezoelectric coefficient and dielectric constant. The resonant oscillation frequency and characteristics are obtained by combining individual solutions.

Building the Langevin transducer model [2], [3], [5], [8]
Figure 1, describing the Langevin transducer structure consists of two PZT hard piezoelectric ceramic plates. The front is an aluminum metal block. The back is a metal block of steel. The whole transducer is linked by steel bolt.
Fig. 1. The Langevin transducer structure [3]

Setting the piezoelectric ceramic plate model and sound pressure field
Figure 2 is the oscillation survey model according to the different frequencies of the PZT piezoelectric hard ceramic to be built from the COMSOL – Multiphysics program, the meshes and nodes distribution with FEM. Sizes: the 40mm outer diameter, the 16mm internal diameter and the 5mm thickness that are also the actual sizes to manufacture the Langevin transducer.
Fig. 2. The PZT piezoelectric ceramic plate model distributing the meshes
and buttons with FEM
The model surveys the displacement and the emitted sound pressure field from the front of the Langevin transducer according to different frequencies, in Figure 3.
Fig. 3. The surveying model of the displacement and the sound pressure
field
After the built model, we simulate the hard PZT piezoelectric ceramic operation and the Langevin transducer
according to frequency. The simulation process purpose is to determine the resonant frequency, the displacement dependence according to the resonant frequency, the emitted sound pressure field characteristic from the Langevin transducer front.

The simulation results to determine the PZT piezoelectric ceramic plate resonant frequency [7], [2]
Determining the resonance frequency, we surveyed the PZT hard piezoelectric ceramic plate oscillation state in the frequency range from 16.855 kHz to 44.334 kHz.
Figure 4 shows the resonant oscillation dependence of the PZT piezoelectric ceramic plae by frequency. The simulation result at the 40 kHz oscillation frequency of the ceramic plate according to the radius in the Figure 4.c is the best match. The ceramic plate oscillates according to the thickness from 16.855 kHz to 27.5 kHz frequencies (Figure 4.a, b) and to be twisted at the 44.334 kHz frequency.
Fig. 4. The PZT piezoelectric ceramic plate resonance oscillation dependence according to frequency
The oscillation spectrum represents the PZT hard piezoelectric impedancedependence by frequency in the Figure 5. The result shows that the impedance reaches the minimum value (Zmin) at the 40 kHz frequency and the resonance phenomenon starts happening. This result is completely suitable for the PZT piezoelectric ceramic plate resonance oscillation at 40 kHz frequency (Figure 4.c).
Fig. 5. The PZT piezoelectric ceramic plate impedance dependence by
frequency

The simulation result to determine the displacement of the Langevin transducer [2], [1]
The results in Figure 6.c and Figure 7 show that the transducer displacement is maximum at 40 kHz frequency.
Fig. 6. The transducer displacement dependence by frequency
Fig. 7. The oscillation spectrum represents the transducer displacement by
frequency

The survey result to determine the sound pressure field of the Langevin transducer front [1], [2]
Figure 8 is the simulation result of the emitted sound pressure field from the Langevin transducer front. The sound pressure field clearly shows the compression – expansion process, the spread of the ultrasonic wave to the frequency dependence environment. At the 40 kHz frequency, the sound pressure is maximun at the center and widened the boundary with the best orientation, Figure 8c.
Fig. 8. The sound pressure field of the Langevin transducer front
CONCLUSION
By finite element method, we have established the continuous equations, kinetic equations, the elastic coefficient matrix , piezoelectric coefficient, dielectric constants of piezoelectric environment.
Using the COMSOLMultiphysics program builds the model to survey the characteristic parameters of the PZT hard piezoelectric ceramic plate and the Langevin ultrasonic transducer.
The resonant frequency of the PZT hardened piezoelectric ceramic is determined at 40 kHz. At this frequency the displacement and the sound pressure field are maximum and the best orientability.
REFERENCES

C. H. Sherman, J. L. Butler, Transducers and Arrays for Underwater Sound.: Springer Science+Business Media, LLC, 2016.

COMSOL Inc. (2011), Introduction to Comsol Multiphysics.: 1998 2017 COMSOL

Dong Chen, Sanjay K. Sharma, Ackmez Mudhoo (2012), Handbook on Applications of Ultrasound, Taylor & Francis Group, LLC

H. Allik and T. J. R. Hughes (1970), "Finite element method for piezoelectric vibration," International Journal of Numerical Methods in Engineering, vol. 2, pp. 151l 57.

Juan A. GallegoJuarez and Karl F. Graff (2015), Power Ultrasonics Applications of Highintensity Ultrasound, Woodhead Publishing.

M. Naillon, R. H. Coursant and F. Besnier (1983), "Analysis of piezoelectric structures by a finite element method," Acta Electronica, vol. 25, pp. 341362.

N. Guo and P. Cawley (1992), "The Finite Element Analysis Of The Vibration Characteristics Of Piezoelectric Discs," Journal of Sound and Vibration, vol. 159, no. 1, pp. 115l 38.

R. S. Woollett, "Theory of the piezoelectric flexural disk transducer with applications to underwater sound," USL research report No.490 S F001 03041, U. S. Navy Underwater Sound Laboratory 1960.

Y. Kagawa and T. Yamabuchi (1976), "A finite element approach to electromechanical problems with an application to energytrapped and surfacewave devices," Transactions of Sonics and Ultrasonics SUU, vol. 3, pp. 263272.