# Radial Vibration Characteristics of Piezoelectric Ceramic Composite Ultrasonic Transducer

Authors: Chundong Yao; Chao Zhang; Jiani Zhang; Yang Zhang
Journal: Shock and Vibration
 (2020)

---

## Abstract

Different from the existing equivalent circuit analysis method of the transducer, based on the vibration theory of the mechanical system and combined with the constitutive equation, this paper analyzes the radial vibration characteristics of the transducer. The piezoelectric ceramic composite ultrasonic transducer is simplified as a mechanical model of a composite thick wall tube composed of a piezoelectric ceramic tube and a metal prestressed tube. The mathematical model of radial vibration of the transducer is established, which consists of the wave equation of radial coupling vibration of the piezoelectric ceramic tube and the metal prestressed tube, the continuity conditions, and the boundary conditions of radial vibration of composite thick wall tube. The characteristic equation and the mode function of radial vibration are derived. The calculated results of natural frequency are in good agreement with the existing experimental results. Based on the analytical method and the difference method, the numerical simulation models of radial vibration are established, and the amplitude-frequency characteristic curves and the displacement responses are given. The simulation results show that the amplitude-frequency characteristic curves and the displacement responses of the two methods are the same, which verifies the correctness of simulation results. Through the simulation analysis, the influence rule of the transducer’s structure sizes on its radial vibration natural frequency is given: when the thickness of the metal prestressed tube and the piezoelectric ceramic tube are constant, the natural frequency decreases with the increase of the inner diameter of the piezoelectric ceramic tube; when the outer diameter of the metal prestressed tube and the inner diameter of the piezoelectric ceramic tube are constant, the natural frequency decreases with the increase of the thickness-to-wall ratio. The calculation method of natural frequency based on elastic vibration theory is clear in concept and simple in calculation, and the simulation models can analyze the mechanical vibration of the transducer.


---

## Body

## 1. Introduction

Ultrasonic transducer is one of the most important parts of an ultrasonic vibration system. Since the Curie brothers discovered the piezoelectric effect in 1880 and the French physicist Langevin invented piezoelectric transducers in 1916 [1, 2], piezoelectric ceramic transducer has gradually become the most widely used type of transducer. Piezoelectric ceramic composite ultrasonic transducer has been the research focus of experts and scholars in recent years. Because of its outstanding advantages of stable performance, large acoustic radiation area, high radiation efficiency, uniform directivity, high sensitivity, and so on [3], it is widely used in underwater acoustic receiving or transmitting [4–6].Piezoelectric ceramic composite ultrasonic transducer is composed of a piezoelectric ceramic tube and a metal prestressed tube along radial direction, and the assembling method is the temperature difference method [7, 8]. In order to obtain the optimal working state of the transducer, the operating frequency of the transducer needs to be its natural frequency [9, 10]. Therefore, it is extremely important to study the radial vibration characteristics of the transducer.In recent years, many scholars have done a lot of researches on the radial vibration characteristics of this kind of transducers. Based on the theory of elasticity, the vibration and stress in the length and radius direction of the thin walled piezoelectric ceramic short tube are ignored; references [11, 12] discussed the natural frequency of transducer through its equivalent circuit. For a transducer whose length is much larger than its radial dimension, reference [13] regarded it as the plane strain problem in mechanics approximately, its pure radial vibration is analyzed, and its equivalent circuit and frequency equation are obtained. According to the principle of electromechanical analogy, reference [14] studied a kind of composite piezoelectricity ultrasonic transducer and obtained the equivalent circuit and frequency equation of the system. Reference [15] proposed a tangential polarized composite cylindrical transducer, which is connected by a tangential polarized piezoelectric tube and an outer metal circular tube, and its natural frequency is analyzed. Reference [16] developed a transducer composed of two radially reinforced composite vibrators, derived the electromechanical equivalent circuit, and obtained the input impedance of the transducer. Although the structures of those transducers studied in the above literatures are different, the essence is the analysis of the coupling vibration between piezoelectric ceramic ring (tube) and metal circular ring (tube). The main research method is to derive the equivalent circuit of transducer combined in the constitutive equation, so as to analyze the natural frequency and performance parameters of the transducer in the resonance state. The calculation process of the above method is complicated, which can only describe the performance parameters of the transducer in the resonance state, not the mechanical vibration. However, dynamic performance is closely related to its sensitivity analysis and structural optimization, which causes certain limitations to the equivalent circuit method.The natural frequency is determined by the structure, the boundary conditions, and the continuous conditions of the vibration system [17]. Based on the vibration theory of the mechanical system and combined with the constitutive equation, the natural frequency and displacement response are obtained. So far, there is no relevant report. Compared with the equivalent circuit method, the calculation process of this method is simpler, and it can describe the mechanical vibration of transducers, which also provides theoretical support for the sensitivity analysis and optimal design of transducers. This paper studies this method. The piezoelectric ceramic composite ultrasonic transducer is simplified as a mechanical model of a composite thick wall tube composed of a piezoelectric ceramic tube and a metal prestressed tube. The mathematical model of radial vibration of the transducer is established, which consists of the wave equation of radial coupling vibration of the piezoelectric ceramic tube and the metal prestressed tube, the continuity conditions, and the boundary conditions of radial vibration of composite thick wall tube. The characteristic equation and mode function of the radial vibration of the transducer are established. Based on the analytical method and the difference method, the numerical simulation models of the radial vibration of the transducer are established, and the amplitude-frequency characteristic curves and the displacement responses of the radial vibration are given. The relationship between the natural frequency of radial vibration and structure sizes was discussed. This paper will provide some guidance for the design and application of the transducer.


## 2. Mechanical and Mathematical Model of Radial Vibration of the Transducer

### 2.1. Mechanical Model

To facilitate research, the following assumptions are made: (1) the axial dimension of the transducer largely outweigh its radial dimension. According to the problem of plane strain in elastic mechanics [18], its axial vibration is ignored, and it is analyzed as a pure radial vibration; (2) the influence of tangential stress is ignored; (3) the influence of the shrink range between the piezoelectric ceramic tube and the metal prestressed tube is ignored.Figure1 shows a schematic diagram of a piezoelectric ceramic composite ultrasonic transducer. It is composed of a piezoelectric ceramic tube and a metal prestressed tube with a negative tolerance by the temperature difference method. A radial excitation voltage is applied to the piezoelectric ceramic tube, which causes the metal prestressed tube to produce radial mechanical vibration by the inverse piezoelectric effect, so that the transducer can produce radial radiation of sound.Figure 1
Schematic diagram of the transducer.Figure2 shows the mechanical model of the transducer. ra, rb are the inside and the outside diameters of the piezoelectric ceramic tube. rb, rc are the inside and the outside diameters of the metal prestressed tube. Pe is the environmental load.Figure 2
Mechanical model of the transducer.


#### 2.1.1. Mechanical Model of the Piezoelectric Ceramic Tube

Figure3 shows the mechanical model of the piezoelectric ceramic tube. Here, l1 is the height of the piezoelectric ceramic tube, Frb1 is the radial force exerted on the piezoelectric ceramic tube by the metal prestressed tube.Figure 3
Mechanical model of the piezoelectric ceramic tube: (a) overall stress of piezoelectric ceramic tube, (b) infinitesimal element, (c)x−y plane force, and (d) r−z plane force.
(a)(b)(c)(d)


#### 2.1.2. Mechanical Model of the Metal Prestressed Tube

Figure4 shows the mechanical model of the metal prestressed tube. Here, l2 is the height of the metal prestressed tube, Frb2 is the radial force exerted on the metal prestressed tube by the piezoelectric ceramic tube, Pe is the environmental load.Figure 4
Mechanical model of the metal prestressed tube: (a) overall stress of metal prestressed tube, (b) infinitesimal element, (c)x−y plane force, and (d) r−z plane force.
(a)(b)(c)(d)


### 2.2. Mathematical Model

The mathematical model of radial vibration of the transducer is composed of the wave equations of radial vibration of the piezoelectric ceramic tube and the metal prestressed tube, the continuity conditions, and the boundary conditions [19].In polar coordinates, the wave equation of radial vibration of the piezoelectric ceramic tube is as follows [20]:(1)ρ1∂2u1∂t2=∂σr∂r+σr−σθr,where u1 is the radial displacement of the piezoelectric ceramic tube, and ρ1 is the density of the piezoelectric ceramic.According to the three-dimensional theory of piezoelectric elasticity, when the axial normal stress and the axial shear stress are not considered, the piezoelectric equation is as follows [21]:(2)εr=S33Eσr+S13Eσθ+d33Er,εθ=S13Eσr+S11Eσθ+d31Er,where S11E, S12ES13E and S33E are elastic compliance coefficients, d31 and d33 are permittivities, and Er is the external excitation voltage.The relationship between strains and displacement is as follows [22]:(3)εr=∂u1∂r,εθ=u1r.The relationship between stresses and displacement can be obtained by combining (2) and (3) are as follows:(4)σr=e11∂u1∂r+e12u1r+e13Er,σθ=e12∂u1∂r+e22u1r+e23Er.Here,(5)e11=S11ES11ES33E−S13ES13E,e12=−S13ES11ES33E−S13ES13E,e13=S13Ed31−S11Ed33S11ES33E−S13ES13E,e22=S33ES11ES33E−S13ES13E,e23=S13Ed33−S33Ed31S11ES33E−S13ES13E.By substituting (4) into (1), the wave equation of radial vibration of the piezoelectric ceramic tube can be expressed as follows:(6)1V12∂2u1∂t2=∂2u1∂r2+1r∂u1∂r−e22e111r2u1+e13−e23e111rEr,where V1=e11/ρ1 is the radial propagation velocity of the sound wave in piezoelectric ceramics.The analysis process of the wave equation of the radial vibration of the metal prestressed tube is the same as that of the piezoelectric ceramic tube. In polar coordinates, the wave equation of radial vibration of the metal prestressed tube is as follows [23]:(7)ρ2∂2u2∂t2=∂σr∂r+σr−σθr,where u2 is the radial displacement of the piezoelectric ceramic tube, and ρ2 is the density of metallic material.According to the thick tube theory [24], the relationship between strains and stresses is as follows:(8)εr=1E2σr−γ2σθ,εθ=1E2σθ−γ2σr,where E2 and γ2 are the elastic modulus and the Poisson ratio of metallic material.The relationship between strains and displacement is as follows:(9)εr=∂u2∂r,εθ=u2r.The relationship between stresses and displacement can be obtained by combining (8) and (9) as follows:(10)σr=c11∂u2∂r+c12u2r,σθ=c12∂u2∂r+c11u2r.Here,(11)c11=E21−γ22,c12=E2γ21−γ22.By substituting (10) into (7), the wave equation of radial vibration of the metal prestressed tube can be expressed as follows:(12)1V22∂2u2∂t2=∂2u2∂r2+1r∂u2∂r−1r2u2,where V2=c11/ρ2 is the radial propagation velocity of the sound wave in metallic materials.The external boundary condition of the metal prestressed tube is a stress boundary condition. Its external load is environmental pressure:(13)σrr=rc=−Pe.The inner boundary condition of the piezoelectric ceramic tube is as follows:(14)σrr=ra=0.The displacement continuity condition and the force continuity condition are as follows:(15)vrb,t=urb,t,S1σrb1=S2σrb2,where S1=2πrbl1 and S2=2πrbl2 are the area of piezoelectric ceramic tube outer wall and the metal prestressed tube inner wall, σrb1 and σrb2 are the radial stress on the outer surface of the piezoelectric ceramic tube and the inner surface of the metal prestressed tube.In conclusion, the mathematical model of radial vibration of the transducer can be obtained as follows:(16)1V12∂2u1∂t2=∂2u1∂r2+1r∂u1∂r−e22e111r2u1+e13−e23e111rEr,1V22∂2u2∂t2=∂2u2∂r2+1r∂u2∂r−1r2u2,σrr=ra=0,σrr=rc=−Pe,u1rb,t=u2rb,t,S1σrb1=S2σrb2.


## 2.1. Mechanical Model

To facilitate research, the following assumptions are made: (1) the axial dimension of the transducer largely outweigh its radial dimension. According to the problem of plane strain in elastic mechanics [18], its axial vibration is ignored, and it is analyzed as a pure radial vibration; (2) the influence of tangential stress is ignored; (3) the influence of the shrink range between the piezoelectric ceramic tube and the metal prestressed tube is ignored.Figure1 shows a schematic diagram of a piezoelectric ceramic composite ultrasonic transducer. It is composed of a piezoelectric ceramic tube and a metal prestressed tube with a negative tolerance by the temperature difference method. A radial excitation voltage is applied to the piezoelectric ceramic tube, which causes the metal prestressed tube to produce radial mechanical vibration by the inverse piezoelectric effect, so that the transducer can produce radial radiation of sound.Figure 1
Schematic diagram of the transducer.Figure2 shows the mechanical model of the transducer. ra, rb are the inside and the outside diameters of the piezoelectric ceramic tube. rb, rc are the inside and the outside diameters of the metal prestressed tube. Pe is the environmental load.Figure 2
Mechanical model of the transducer.


### 2.1.1. Mechanical Model of the Piezoelectric Ceramic Tube

Figure3 shows the mechanical model of the piezoelectric ceramic tube. Here, l1 is the height of the piezoelectric ceramic tube, Frb1 is the radial force exerted on the piezoelectric ceramic tube by the metal prestressed tube.Figure 3
Mechanical model of the piezoelectric ceramic tube: (a) overall stress of piezoelectric ceramic tube, (b) infinitesimal element, (c)x−y plane force, and (d) r−z plane force.
(a)(b)(c)(d)


### 2.1.2. Mechanical Model of the Metal Prestressed Tube

Figure4 shows the mechanical model of the metal prestressed tube. Here, l2 is the height of the metal prestressed tube, Frb2 is the radial force exerted on the metal prestressed tube by the piezoelectric ceramic tube, Pe is the environmental load.Figure 4
Mechanical model of the metal prestressed tube: (a) overall stress of metal prestressed tube, (b) infinitesimal element, (c)x−y plane force, and (d) r−z plane force.
(a)(b)(c)(d)


## 2.1.1. Mechanical Model of the Piezoelectric Ceramic Tube

Figure3 shows the mechanical model of the piezoelectric ceramic tube. Here, l1 is the height of the piezoelectric ceramic tube, Frb1 is the radial force exerted on the piezoelectric ceramic tube by the metal prestressed tube.Figure 3
Mechanical model of the piezoelectric ceramic tube: (a) overall stress of piezoelectric ceramic tube, (b) infinitesimal element, (c)x−y plane force, and (d) r−z plane force.
(a)(b)(c)(d)


## 2.1.2. Mechanical Model of the Metal Prestressed Tube

Figure4 shows the mechanical model of the metal prestressed tube. Here, l2 is the height of the metal prestressed tube, Frb2 is the radial force exerted on the metal prestressed tube by the piezoelectric ceramic tube, Pe is the environmental load.Figure 4
Mechanical model of the metal prestressed tube: (a) overall stress of metal prestressed tube, (b) infinitesimal element, (c)x−y plane force, and (d) r−z plane force.
(a)(b)(c)(d)


## 2.2. Mathematical Model

The mathematical model of radial vibration of the transducer is composed of the wave equations of radial vibration of the piezoelectric ceramic tube and the metal prestressed tube, the continuity conditions, and the boundary conditions [19].In polar coordinates, the wave equation of radial vibration of the piezoelectric ceramic tube is as follows [20]:(1)ρ1∂2u1∂t2=∂σr∂r+σr−σθr,where u1 is the radial displacement of the piezoelectric ceramic tube, and ρ1 is the density of the piezoelectric ceramic.According to the three-dimensional theory of piezoelectric elasticity, when the axial normal stress and the axial shear stress are not considered, the piezoelectric equation is as follows [21]:(2)εr=S33Eσr+S13Eσθ+d33Er,εθ=S13Eσr+S11Eσθ+d31Er,where S11E, S12ES13E and S33E are elastic compliance coefficients, d31 and d33 are permittivities, and Er is the external excitation voltage.The relationship between strains and displacement is as follows [22]:(3)εr=∂u1∂r,εθ=u1r.The relationship between stresses and displacement can be obtained by combining (2) and (3) are as follows:(4)σr=e11∂u1∂r+e12u1r+e13Er,σθ=e12∂u1∂r+e22u1r+e23Er.Here,(5)e11=S11ES11ES33E−S13ES13E,e12=−S13ES11ES33E−S13ES13E,e13=S13Ed31−S11Ed33S11ES33E−S13ES13E,e22=S33ES11ES33E−S13ES13E,e23=S13Ed33−S33Ed31S11ES33E−S13ES13E.By substituting (4) into (1), the wave equation of radial vibration of the piezoelectric ceramic tube can be expressed as follows:(6)1V12∂2u1∂t2=∂2u1∂r2+1r∂u1∂r−e22e111r2u1+e13−e23e111rEr,where V1=e11/ρ1 is the radial propagation velocity of the sound wave in piezoelectric ceramics.The analysis process of the wave equation of the radial vibration of the metal prestressed tube is the same as that of the piezoelectric ceramic tube. In polar coordinates, the wave equation of radial vibration of the metal prestressed tube is as follows [23]:(7)ρ2∂2u2∂t2=∂σr∂r+σr−σθr,where u2 is the radial displacement of the piezoelectric ceramic tube, and ρ2 is the density of metallic material.According to the thick tube theory [24], the relationship between strains and stresses is as follows:(8)εr=1E2σr−γ2σθ,εθ=1E2σθ−γ2σr,where E2 and γ2 are the elastic modulus and the Poisson ratio of metallic material.The relationship between strains and displacement is as follows:(9)εr=∂u2∂r,εθ=u2r.The relationship between stresses and displacement can be obtained by combining (8) and (9) as follows:(10)σr=c11∂u2∂r+c12u2r,σθ=c12∂u2∂r+c11u2r.Here,(11)c11=E21−γ22,c12=E2γ21−γ22.By substituting (10) into (7), the wave equation of radial vibration of the metal prestressed tube can be expressed as follows:(12)1V22∂2u2∂t2=∂2u2∂r2+1r∂u2∂r−1r2u2,where V2=c11/ρ2 is the radial propagation velocity of the sound wave in metallic materials.The external boundary condition of the metal prestressed tube is a stress boundary condition. Its external load is environmental pressure:(13)σrr=rc=−Pe.The inner boundary condition of the piezoelectric ceramic tube is as follows:(14)σrr=ra=0.The displacement continuity condition and the force continuity condition are as follows:(15)vrb,t=urb,t,S1σrb1=S2σrb2,where S1=2πrbl1 and S2=2πrbl2 are the area of piezoelectric ceramic tube outer wall and the metal prestressed tube inner wall, σrb1 and σrb2 are the radial stress on the outer surface of the piezoelectric ceramic tube and the inner surface of the metal prestressed tube.In conclusion, the mathematical model of radial vibration of the transducer can be obtained as follows:(16)1V12∂2u1∂t2=∂2u1∂r2+1r∂u1∂r−e22e111r2u1+e13−e23e111rEr,1V22∂2u2∂t2=∂2u2∂r2+1r∂u2∂r−1r2u2,σrr=ra=0,σrr=rc=−Pe,u1rb,t=u2rb,t,S1σrb1=S2σrb2.


## 3. Analytical Solution of Radial Vibration of Transducer

### 3.1. Vibration Mode Function and Natural Frequency

The natural frequency is the inherent property of the object, which is independent of the electric field excitation and the external environment pressure [25]. When the frequency equation needs to be derived, the transducer should be in a free vibration state, and there is no voltage excitation and external environmental pressure.The separation of displacement functions is done by splitting space and temporal variables [26]. Let u1=ur1ejωtu2=ur2ejωt, equation (16) can be expressed as follows:(17)r2d2ur1dr2+rdur1dr+k12r2−m2ur1=0,r2d2ur2dr2+rdur2dr+k22r2−1ur2=0,where ω is the angular frequency, k1=ω/V1 and k2=ω/V2 are the wave numbers of radial vibration of the piezoelectric ceramic tube vibrator and the metal prestressed tube vibrator, and m=e22/e11.Equation (18) is a Bessel equation of order 1 and order m [27]. Its general solution is as follows:(18)ur1=A1Jmk1r+B1Ymk1r,ur2=A2J1k2r+B2Y1k2r.By substituting equation (19) into the boundary conditions and the continuity conditions, the following conclusions can be drawn.At the inner diameter of the piezoelectric ceramic tube, the stress boundary condition is as follows:(19)e11A1Jm′k1ra+B1Ym′k1ra+e12raA1Jmk1ra+B1Ymk1ra=0.At the outer diameter of the metal prestressed tube, the stress boundary condition is as follows:(20)c11A2J1′k2rc+B2Y1′k2rc+c12rcA2J1k2rc+B2Y1k2rc=0.The displacement continuity condition and the force continuity condition is as follows:(21)A1Jmk1rb+B1Ymk1rb=A2J1k2rb+B2Y1k2rb,S1e11A1Jm′k1rb+B1Ym′k1rb+e12rbA1Jmk1rb+B1Ymk1rb,=S2c11A2J1′k2rb+B2Y1′k2rb+S2c12rbA2J1k2rb+B2Y1k2rb.The characteristic equation of radial vibration of the transducer can be obtained by combining (20), (21), and (22):(22)S1e11Ym′k1rb+e12/rbYmk1rb−e11Jm′k1rb+e12/rbJmk1rbe11Ym′k1ra+e12/rbYmk1ra/e11Jm′k1ra+e12/rbJmk1raYmk1rb−Jmk1rbe11Ym′k1rb+e12/rbYmk1ra/e11Jm′k1ra+e12/rbJmk1ra=S2c11Y1′k2rb+c12/rbY1k2rb−c11J1′k2rb+c12/rbJ1k2rbc11Y1′k2rc+c12/rbY1k2rc/c11J1′k2rc+c12/rbJ1k2rcY1k2rb−J1k2rbc11Y1′k2rc+c12/rbY1k2rc/c11J1′k2rc+c12/rbJvk2rc.The natural frequency can be calculated by the characteristic equation [28]. The computing method is as follows:(23)fn=ωn/2π,n=1,2,3,….It can be obtained by combining (20)–(22):(24)A1=−e11Ym′k1ra+e12/raYmk1rae11Jmk1ra+e12/raJmk1raB1,A2=−c11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rcB2,B1B2=Y1k2rb−J1k2rbc11Y1′k2rc+c12/rcY1k2rc/c11J1′k2rc+k1raJvk2rcYmk1rb−Jmk1rbe11Ym′k1ra+e12/raYmk1ra/e11Jm′k1ra+e12/raJmk1ra=S2c11Y1′k2rb+c12/rbY1k2rb−c11J1′k2rb+c12/rbJ1k2rbc11Y1′k2rc+c12/rcY1k2rc/c11J1′k2rc+c12/rcJ1k2rcS1e11Ym′k1rb+e12/rbYmk1rb−e11Jm′k1rb+e12/rbJmk1rbe11Ym′k1ra+e12/raYmk1ra/e11Jm′k1ra+e12/raJmk1ra.Therefore, the mode function is as follows:(25)urn1=Y1k2nrb−J1k2nrbc11Y1′k2nrc+c12/rcY1k2nrc/c11J1′k2nrc+c12/rcJ1k2nrcYmk1nrb−Jmk1nrbe11Ym′k1nra+e12/raYmk1nra/e11Jm′k1nra+e12/raJmk1nrae11Ym′k1nra+e12/rbYmk1nrae11Jm′k1nra+e12/rbJmk1nraJmk1nr+Ymk1nr,urn2=−c11Y1′k2nrc+c12/rcY1k2nrcc11J1′k2nrc+c12/rcJ1k2nrcJ1k2nr+Y1k2nr.


## 3.1. Vibration Mode Function and Natural Frequency

The natural frequency is the inherent property of the object, which is independent of the electric field excitation and the external environment pressure [25]. When the frequency equation needs to be derived, the transducer should be in a free vibration state, and there is no voltage excitation and external environmental pressure.The separation of displacement functions is done by splitting space and temporal variables [26]. Let u1=ur1ejωtu2=ur2ejωt, equation (16) can be expressed as follows:(17)r2d2ur1dr2+rdur1dr+k12r2−m2ur1=0,r2d2ur2dr2+rdur2dr+k22r2−1ur2=0,where ω is the angular frequency, k1=ω/V1 and k2=ω/V2 are the wave numbers of radial vibration of the piezoelectric ceramic tube vibrator and the metal prestressed tube vibrator, and m=e22/e11.Equation (18) is a Bessel equation of order 1 and order m [27]. Its general solution is as follows:(18)ur1=A1Jmk1r+B1Ymk1r,ur2=A2J1k2r+B2Y1k2r.By substituting equation (19) into the boundary conditions and the continuity conditions, the following conclusions can be drawn.At the inner diameter of the piezoelectric ceramic tube, the stress boundary condition is as follows:(19)e11A1Jm′k1ra+B1Ym′k1ra+e12raA1Jmk1ra+B1Ymk1ra=0.At the outer diameter of the metal prestressed tube, the stress boundary condition is as follows:(20)c11A2J1′k2rc+B2Y1′k2rc+c12rcA2J1k2rc+B2Y1k2rc=0.The displacement continuity condition and the force continuity condition is as follows:(21)A1Jmk1rb+B1Ymk1rb=A2J1k2rb+B2Y1k2rb,S1e11A1Jm′k1rb+B1Ym′k1rb+e12rbA1Jmk1rb+B1Ymk1rb,=S2c11A2J1′k2rb+B2Y1′k2rb+S2c12rbA2J1k2rb+B2Y1k2rb.The characteristic equation of radial vibration of the transducer can be obtained by combining (20), (21), and (22):(22)S1e11Ym′k1rb+e12/rbYmk1rb−e11Jm′k1rb+e12/rbJmk1rbe11Ym′k1ra+e12/rbYmk1ra/e11Jm′k1ra+e12/rbJmk1raYmk1rb−Jmk1rbe11Ym′k1rb+e12/rbYmk1ra/e11Jm′k1ra+e12/rbJmk1ra=S2c11Y1′k2rb+c12/rbY1k2rb−c11J1′k2rb+c12/rbJ1k2rbc11Y1′k2rc+c12/rbY1k2rc/c11J1′k2rc+c12/rbJ1k2rcY1k2rb−J1k2rbc11Y1′k2rc+c12/rbY1k2rc/c11J1′k2rc+c12/rbJvk2rc.The natural frequency can be calculated by the characteristic equation [28]. The computing method is as follows:(23)fn=ωn/2π,n=1,2,3,….It can be obtained by combining (20)–(22):(24)A1=−e11Ym′k1ra+e12/raYmk1rae11Jmk1ra+e12/raJmk1raB1,A2=−c11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rcB2,B1B2=Y1k2rb−J1k2rbc11Y1′k2rc+c12/rcY1k2rc/c11J1′k2rc+k1raJvk2rcYmk1rb−Jmk1rbe11Ym′k1ra+e12/raYmk1ra/e11Jm′k1ra+e12/raJmk1ra=S2c11Y1′k2rb+c12/rbY1k2rb−c11J1′k2rb+c12/rbJ1k2rbc11Y1′k2rc+c12/rcY1k2rc/c11J1′k2rc+c12/rcJ1k2rcS1e11Ym′k1rb+e12/rbYmk1rb−e11Jm′k1rb+e12/rbJmk1rbe11Ym′k1ra+e12/raYmk1ra/e11Jm′k1ra+e12/raJmk1ra.Therefore, the mode function is as follows:(25)urn1=Y1k2nrb−J1k2nrbc11Y1′k2nrc+c12/rcY1k2nrc/c11J1′k2nrc+c12/rcJ1k2nrcYmk1nrb−Jmk1nrbe11Ym′k1nra+e12/raYmk1nra/e11Jm′k1nra+e12/raJmk1nrae11Ym′k1nra+e12/rbYmk1nrae11Jm′k1nra+e12/rbJmk1nraJmk1nr+Ymk1nr,urn2=−c11Y1′k2nrc+c12/rcY1k2nrcc11J1′k2nrc+c12/rcJ1k2nrcJ1k2nr+Y1k2nr.


## 4. Displacement

### 4.1. Analytical Solution

Letu1=ur1ejωt,u2=ur2ejωt, and Er=E0ejωt; equation (17) can be expressed as follows:(26)r2d2ur1dr2+rdur1dr+k12r2−m2ur1+e13−e23e11E0=0,r2d2ur2dr2+rdur2dr+k22r2−1ur2=0.Equation (26) is the Bessel equation and Lommel equation. Its general solution is as follows:(27)ur1=H1Jmk1r+G1Ymk1r+e23−e13e11E0S0,mk1r,ur2=H2J1k2r+G2Y1k2r.The boundary conditions and continuity conditions can be expressed as follows:(28)e11H1Jm′k1ra+G1Ym′k1ra+e12raH1Jmk1ra+G1Ymk1ra+e23−e13S0,v′k1ra+e12rae23−e13e11S0,mk1ra+e13E0=0,c11H2J1′k2rc+G2Y1′k2rc+c12rcH2J1k2rc+G2Y1k2rc=−Pe,H1Jmk1rb+G1Ymk1rb+e23−e13e11E0S0,mk1rb=H2J1k2rb+G2Y1k2rb,S1e11H2Jm′k1rb+G1Ym′k1rb+e12rbH1Jmk1rb+G1Ymk1rb+S1e23−e13S0,m′k1ra+e12rae23−e13e11S0,mk1ra+e13E0=S2c11H2J1′k2rb+G2Y1′k2rb+S2c12rbH2J1k2rb+G2Y1k2rb.According to the boundary conditions and the continuity conditions, the response of the system is as follows:(29)u1r,t=H1Jmk1r+G1Ymk1r+e23−e13e11E0S0,mk1rejwt,u2r,t=H2J1k2r+G2Y1k2rejwt.Here,(30)G1=F22F21G2+F23F21,G2=F23/F21−F13/F11F12/F11−F22/F21,H1=−e11Ym′k1ra+e12/rbYmk1rae11Jm′k1ra+e12/rbJmk1raG1−e23−e13E0S0,m′k1ra+e12/e11e23−e13/raE0S0,mk1ra+e13E0e11Jm′k1ra+e12/rbJmk1ra,H2=−c11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rcG2−Pec11J1′k2rc+c12/rbJ1k2rc,F11=Ymk1rb−Jmk1rbe11Ym′k1ra+e12/rbYmk1rae11Jm′k1ra+e12/rbJmk1ra,F21=S1e11Ym′k1rb+e12/rbYmk1rb−S1e11Jm′k1rb+e12/rbJmk1rbe11Ym′k1ra+e12/rbYmk1rae11Jm′k1ra+e12/rbJmk1ra,F12=Y1k2rb−J1k2rbc11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rc,F22=S2c11Y1′k2rb+c12/rbY1k2rb−S2c11J1′k2rb+c12/rbJ1k2rbc11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rc,F13=e23−e13e11E0S0,mk1ra−e23−e130S0,m′k1ra+e12/e11e23−e13/raE0S0,mk1ra+e13E0e11Jm′k1ra+e12/rbJmk1ra+J1k2rbPec11J1′k2rc+c12/rbJ1k2rc,F23=S1e23−e13E0S0,m′k1rb+e12e11e23−e13rbE0S0,mk1rb+e13E0−S1e11Jm′k1rb+e12/rbJmk1rbe11Jm′k1rb+e12/rbJmk1rbS2e23−e13E0S0,m′k1ra+e12e11e23−e13rbE0S0,mk1ra+e13E0+S2c11J1′k1rb+c12rbJ1k1rbPec11J1′k2rc+c12/rbJ1k2rc.


### 4.2. Difference Solution

The transducer is discretized along the radial direction. The planer−z is the research plane, the transducer is divided into I nodes along the radial direction, and the unit length is dr. The inner wall of the piezoelectric ceramic tube is the first node, the outer wall of the metal prestressed tube is the node I, and the composite position is the node i0. Figure 5 shows the discrete model of the radial vibration of the transducer. Time is divided into J+1 time nodes. ui,j,i=0,1,2,3,…,I represents the radial vibration displacement of the i space node at the j time node.Figure 5
Discrete model of radial vibration of the transducer.The numerical simulation model can be obtained by using the forward difference formula and the central difference formula [29].The difference equations are as follows:(31)ui,j+11=v22Δt21/Δr2ui+1,j−2ui,j+ui−1,j+1/ra+iΔrΔrui+1,j−ui,j−e22/e111/ra+iΔr2ui,j+e13−e23/e111/ra+iΔrErj+1+2ui,j−ui,j−1,ui,j+1=v22Δt21/Δr2ui+1,j−2ui,j+ui−1,j+1/ra+iΔrΔrui+1,j−ui,j−e22/e111/ra+iΔr2ui,j+e13−e23/e111/ra+iΔrErj+1+2ui,j−ui,j−1.The difference form of the boundary conditions is as follows:(32)u0,j+1=1e11/Δr−e12/rae11Δru1,j+1+e13Erj+1,uI,j+1=1c11/Δr−c12/rcc11ΔruI−1,j+1.The difference equation of the composite position is as follows:(33)ui0,j+1=S1c111/Δr−1/2rbvi0+1,j+S1c111/2rb−1/Δr+Δr/2rb2−S2e111/Δr−1/2rb+Δre22/2e11rb2−S2e12/rb+S1c12/rbui0,j−S2e111/2rb−1/Δrui0−1,j−S2e13Erj+11/2Δr1/Δt2S2e111/v12−S1c111/v22.The initial displacement and velocity are 0. The difference form of the initial conditions is as follows:(34)ui,01=0,ui,02=0,ui,11=12ui,2,ui,12=12ui,22.Based on the numerical simulation model established above, the difference grid as shown in Figure6 can be established.Figure 6
Difference grid.As shown in Figure6, rectangle “a” represents the numerical calculation path of the inner boundary (the first node) of the difference method: the displacement of the first node is calculated from the second node by the difference form of the inner condition. Rectangle “b” represents the numerical calculation path of the outer boundary (the last node): the displacement of the last node is calculated from the previous node by the difference form of the outer condition. Rectangle “c” and “d” represent the numerical calculation path of internal nodes: the displacement of an internal node is calculated by the displacement of three nodes at the previous time node and one node at the previous two-time node.


## 4.1. Analytical Solution

Letu1=ur1ejωt,u2=ur2ejωt, and Er=E0ejωt; equation (17) can be expressed as follows:(26)r2d2ur1dr2+rdur1dr+k12r2−m2ur1+e13−e23e11E0=0,r2d2ur2dr2+rdur2dr+k22r2−1ur2=0.Equation (26) is the Bessel equation and Lommel equation. Its general solution is as follows:(27)ur1=H1Jmk1r+G1Ymk1r+e23−e13e11E0S0,mk1r,ur2=H2J1k2r+G2Y1k2r.The boundary conditions and continuity conditions can be expressed as follows:(28)e11H1Jm′k1ra+G1Ym′k1ra+e12raH1Jmk1ra+G1Ymk1ra+e23−e13S0,v′k1ra+e12rae23−e13e11S0,mk1ra+e13E0=0,c11H2J1′k2rc+G2Y1′k2rc+c12rcH2J1k2rc+G2Y1k2rc=−Pe,H1Jmk1rb+G1Ymk1rb+e23−e13e11E0S0,mk1rb=H2J1k2rb+G2Y1k2rb,S1e11H2Jm′k1rb+G1Ym′k1rb+e12rbH1Jmk1rb+G1Ymk1rb+S1e23−e13S0,m′k1ra+e12rae23−e13e11S0,mk1ra+e13E0=S2c11H2J1′k2rb+G2Y1′k2rb+S2c12rbH2J1k2rb+G2Y1k2rb.According to the boundary conditions and the continuity conditions, the response of the system is as follows:(29)u1r,t=H1Jmk1r+G1Ymk1r+e23−e13e11E0S0,mk1rejwt,u2r,t=H2J1k2r+G2Y1k2rejwt.Here,(30)G1=F22F21G2+F23F21,G2=F23/F21−F13/F11F12/F11−F22/F21,H1=−e11Ym′k1ra+e12/rbYmk1rae11Jm′k1ra+e12/rbJmk1raG1−e23−e13E0S0,m′k1ra+e12/e11e23−e13/raE0S0,mk1ra+e13E0e11Jm′k1ra+e12/rbJmk1ra,H2=−c11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rcG2−Pec11J1′k2rc+c12/rbJ1k2rc,F11=Ymk1rb−Jmk1rbe11Ym′k1ra+e12/rbYmk1rae11Jm′k1ra+e12/rbJmk1ra,F21=S1e11Ym′k1rb+e12/rbYmk1rb−S1e11Jm′k1rb+e12/rbJmk1rbe11Ym′k1ra+e12/rbYmk1rae11Jm′k1ra+e12/rbJmk1ra,F12=Y1k2rb−J1k2rbc11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rc,F22=S2c11Y1′k2rb+c12/rbY1k2rb−S2c11J1′k2rb+c12/rbJ1k2rbc11Y1′k2rc+c12/rbY1k2rcc11J1′k2rc+c12/rbJ1k2rc,F13=e23−e13e11E0S0,mk1ra−e23−e130S0,m′k1ra+e12/e11e23−e13/raE0S0,mk1ra+e13E0e11Jm′k1ra+e12/rbJmk1ra+J1k2rbPec11J1′k2rc+c12/rbJ1k2rc,F23=S1e23−e13E0S0,m′k1rb+e12e11e23−e13rbE0S0,mk1rb+e13E0−S1e11Jm′k1rb+e12/rbJmk1rbe11Jm′k1rb+e12/rbJmk1rbS2e23−e13E0S0,m′k1ra+e12e11e23−e13rbE0S0,mk1ra+e13E0+S2c11J1′k1rb+c12rbJ1k1rbPec11J1′k2rc+c12/rbJ1k2rc.


## 4.2. Difference Solution

The transducer is discretized along the radial direction. The planer−z is the research plane, the transducer is divided into I nodes along the radial direction, and the unit length is dr. The inner wall of the piezoelectric ceramic tube is the first node, the outer wall of the metal prestressed tube is the node I, and the composite position is the node i0. Figure 5 shows the discrete model of the radial vibration of the transducer. Time is divided into J+1 time nodes. ui,j,i=0,1,2,3,…,I represents the radial vibration displacement of the i space node at the j time node.Figure 5
Discrete model of radial vibration of the transducer.The numerical simulation model can be obtained by using the forward difference formula and the central difference formula [29].The difference equations are as follows:(31)ui,j+11=v22Δt21/Δr2ui+1,j−2ui,j+ui−1,j+1/ra+iΔrΔrui+1,j−ui,j−e22/e111/ra+iΔr2ui,j+e13−e23/e111/ra+iΔrErj+1+2ui,j−ui,j−1,ui,j+1=v22Δt21/Δr2ui+1,j−2ui,j+ui−1,j+1/ra+iΔrΔrui+1,j−ui,j−e22/e111/ra+iΔr2ui,j+e13−e23/e111/ra+iΔrErj+1+2ui,j−ui,j−1.The difference form of the boundary conditions is as follows:(32)u0,j+1=1e11/Δr−e12/rae11Δru1,j+1+e13Erj+1,uI,j+1=1c11/Δr−c12/rcc11ΔruI−1,j+1.The difference equation of the composite position is as follows:(33)ui0,j+1=S1c111/Δr−1/2rbvi0+1,j+S1c111/2rb−1/Δr+Δr/2rb2−S2e111/Δr−1/2rb+Δre22/2e11rb2−S2e12/rb+S1c12/rbui0,j−S2e111/2rb−1/Δrui0−1,j−S2e13Erj+11/2Δr1/Δt2S2e111/v12−S1c111/v22.The initial displacement and velocity are 0. The difference form of the initial conditions is as follows:(34)ui,01=0,ui,02=0,ui,11=12ui,2,ui,12=12ui,22.Based on the numerical simulation model established above, the difference grid as shown in Figure6 can be established.Figure 6
Difference grid.As shown in Figure6, rectangle “a” represents the numerical calculation path of the inner boundary (the first node) of the difference method: the displacement of the first node is calculated from the second node by the difference form of the inner condition. Rectangle “b” represents the numerical calculation path of the outer boundary (the last node): the displacement of the last node is calculated from the previous node by the difference form of the outer condition. Rectangle “c” and “d” represent the numerical calculation path of internal nodes: the displacement of an internal node is calculated by the displacement of three nodes at the previous time node and one node at the previous two-time node.


## 5. Example Analysis

The material of the piezoelectric ceramic tube is PZT-4. Its material parameters are as follows:ρ1 = 7500 kg/m3, S11E = 12.3 × 10−12 m2/N, S12E = −4.05 × 10−12 m2/N, S13E = −5.31 × 10−12 m2/N, and S33E = 1.55 × 10−12 m2/N. The material of the metal prestressed tube is aluminum alloy. Its material parameters are as follows: ρ2 = 2700 kg/m3, E2 = 7.15 × 1010 N/m3, and γ2=0.34.In order to verify the correctness of the simulation models, the above transducer is calculated as an example.


### 5.1. Natural Frequency Calculation

The natural frequencies of the three groups of transducers are solved by using the characteristic equation (23), which is also compared with the existing experimental results in reference [13]. As shown in Table 1, the theoretical calculating values ft are in good agreement with the existing experimental results fe in the reference, and the accuracy of the obtained characteristic equation meets the requirements of engineering design.Table 1
Theoretical and existing values in the reference.
No.ra (mm)rb (mm)rc (mm)l1, l2 (mm)ft (Hz)fe (Hz)Δ (%)11115213442833430420.4822126313424597256213.9935058638410704108511.88


### 5.2. Amplitude-Frequency Characteristics and Displacement Response

When the excitation frequency of the external voltage is consistent with the natural frequency of the transducer, the transducer will resonate, and its displacement response will suddenly increase [30]. Based on the resonance method, the amplitude-frequency characteristics of the analytical solution and the difference solution are given, respectively.The structure sizes of the transducer are as follows:ra = 21 mm,rb = 26 mm,rc = 31 mm,l1=l2 = 34 mm. A 50 V AC voltage is applied to this transducer. For the analytical solution and the difference solution, the amplitude-frequency characteristics of the transducer at 1.2 × 10−5s are obtained by iteration in 1 Hz step.Figure7 shows the frequency characteristic curve of the analytical solution. It can be observed from Figure 7 that the displacement suddenly increases when the frequency of the voltage of the transducer is 24597 Hz. It can be concluded that the natural frequency of the radial vibration of the transducer is 24597 Hz, which meets the value obtained by the characteristic equation.Figure 7
Amplitude-frequency characteristics of the analytical solution.Limited by the time length and the accuracy of computer calculation, letdt = 3 × 10−9 s, dr = 5 × 10−5 m. The amplitude-frequency characteristic of the difference solution at 1.2 × 10−5 s can be obtained by calculating the 4000th time node and the 101st space node. Figure 8 shows the frequency characteristic curve of the difference solution. It can be concluded that the natural frequency of the radial vibration of the transducer is 24597 Hz, which is consistent with the conclusion of the characteristic equation and the analytical solution.Figure 8
Amplitude-frequency characteristics of the difference solution: (a) there-dimensional diagram; (b) the amplitude-frequency characteristics.
(a)(b)For the transducer of the same structure sizes, when a radial excitation voltage of 20000 Hz and 50 V is applied to the transducer, the displacement responses of the analytical solution and the differential solution are as shown in Figure9.Figure 9
Displacement.It can be observed from Figure9 that the radial vibration displacement curves obtained by the analytical method, and the difference method are same when the transducer with same structure sizes and same external excitation voltage, which realizes the mutual verification between the analytical solution and the differential solution of radial vibration of the transducer and proves the correctness of the simulation result.


## 5.1. Natural Frequency Calculation

The natural frequencies of the three groups of transducers are solved by using the characteristic equation (23), which is also compared with the existing experimental results in reference [13]. As shown in Table 1, the theoretical calculating values ft are in good agreement with the existing experimental results fe in the reference, and the accuracy of the obtained characteristic equation meets the requirements of engineering design.Table 1
Theoretical and existing values in the reference.
No.ra (mm)rb (mm)rc (mm)l1, l2 (mm)ft (Hz)fe (Hz)Δ (%)11115213442833430420.4822126313424597256213.9935058638410704108511.88


## 5.2. Amplitude-Frequency Characteristics and Displacement Response

When the excitation frequency of the external voltage is consistent with the natural frequency of the transducer, the transducer will resonate, and its displacement response will suddenly increase [30]. Based on the resonance method, the amplitude-frequency characteristics of the analytical solution and the difference solution are given, respectively.The structure sizes of the transducer are as follows:ra = 21 mm,rb = 26 mm,rc = 31 mm,l1=l2 = 34 mm. A 50 V AC voltage is applied to this transducer. For the analytical solution and the difference solution, the amplitude-frequency characteristics of the transducer at 1.2 × 10−5s are obtained by iteration in 1 Hz step.Figure7 shows the frequency characteristic curve of the analytical solution. It can be observed from Figure 7 that the displacement suddenly increases when the frequency of the voltage of the transducer is 24597 Hz. It can be concluded that the natural frequency of the radial vibration of the transducer is 24597 Hz, which meets the value obtained by the characteristic equation.Figure 7
Amplitude-frequency characteristics of the analytical solution.Limited by the time length and the accuracy of computer calculation, letdt = 3 × 10−9 s, dr = 5 × 10−5 m. The amplitude-frequency characteristic of the difference solution at 1.2 × 10−5 s can be obtained by calculating the 4000th time node and the 101st space node. Figure 8 shows the frequency characteristic curve of the difference solution. It can be concluded that the natural frequency of the radial vibration of the transducer is 24597 Hz, which is consistent with the conclusion of the characteristic equation and the analytical solution.Figure 8
Amplitude-frequency characteristics of the difference solution: (a) there-dimensional diagram; (b) the amplitude-frequency characteristics.
(a)(b)For the transducer of the same structure sizes, when a radial excitation voltage of 20000 Hz and 50 V is applied to the transducer, the displacement responses of the analytical solution and the differential solution are as shown in Figure9.Figure 9
Displacement.It can be observed from Figure9 that the radial vibration displacement curves obtained by the analytical method, and the difference method are same when the transducer with same structure sizes and same external excitation voltage, which realizes the mutual verification between the analytical solution and the differential solution of radial vibration of the transducer and proves the correctness of the simulation result.


## 6. Relationship between Natural Frequencies and Structural Dimensions

Continuum mechanics models depicted that thickness plays a significant role for vibrational analysis of tubes [31, 32]. Taking the transducer of metal prestressed tube as steel and piezoelectric ceramic tube as PZT-4 as an example, the relationship between natural frequency and structure sizes is given by programming with MATLAB.Its material parameters are as follows:ρ1 = 7500 kg/m3, ρ2 = 7930 kg/m3, E2 = 209 × 1010 N/m3, γ2 = 0.28, S11E =12.3 × 10−11 m2/N, S12E = −4.05 × 10−11 m2/N, S13E = −5.31 ×10−12 m2/N, and S33E = 1.55 × 10−12 m2/N.Keep the wall thickness of the piezoelectric ceramic tube and the metal prestressed tube as 5 mm, and take the inner diameter of the piezoelectric ceramic tube as 18–30 mm. Other structure sizes are as follows:l1 = 40 mm, l2 = 30 mm. For the above group of transducers, the relationship curve between the natural frequency and the inner diameter obtained by using the derived characteristic equation is as shown in Figure 10. It can be observed from Figure 10, when the wall thickness of the metal prestressed tube and the piezoelectric ceramic tube is kept constant, the larger the inner diameter of the piezoelectric ceramic tube, the smaller the natural frequency of radial vibration of the transducer.Figure 10
Relationship between natural frequency and inner diameterKeep the wall thickness of the piezoelectric ceramic tube and the metal prestressed tube, the outer diameter of the metal prestressed tube, and the inner diameter of the piezoelectric ceramic tube unchanged. The only change is the inner diameter of the metal prestressed tube or the outer diameter of the piezoelectric ceramic tube. The thickness-to-wall ratioλ=rb−ra/rc−ra is introduced, which represents the proportion of the thickness of the piezoelectric ceramic tube to the total thickness of the transducer. When λ is 0, the transducer becomes a metal prestressed tube, and when λ is 1, the transducer becomes a piezoelectric ceramic tube. The relationship between the natural frequency and the thickness-to-wall ratio can be obtained by using the derived characteristic equation.Figure11 shows the relationship between the natural frequency and the wall thickness-to-wall ratio. It can be observed from Figure 11 that the natural frequency decreases with the increase of the thickness-to-wall ratio. The reason is that the elastic modulus of piezoelectric ceramic material is smaller than that of metal material. The larger thickness of the piezoelectric ceramic tube equates with the larger proportion of piezoelectric ceramic material in the transducer, so the natural frequency will decrease.Figure 11
Relationship between natural frequency and wall thickness ratio.


## 7. Conclusions

(1)
The focus in the paper is on the radial vibration characteristics of the piezoelectric ceramic composite ultrasonic transducer. The transducer is simplified as a mechanical model of a composite thick wall tube composed of a piezoelectric ceramic tube and a metal prestressed tube. The mathematical model of radial vibration of the transducer is established, which consists of the wave equation of radial coupling vibration of the piezoelectric ceramic tube and the metal prestressed tube, the continuity conditions, and the boundary conditions of radial vibration of composite thick wall tube. The characteristic equation and the mode function of the radial vibration of the transducer are obtained. The results of natural frequency calculations are consistent with the existing experimental results, which proves the correctness of the simulation models and the accuracy of the obtained characteristic equation meets the requirement of engineering design.(2)
Based on the analytical method and the difference method, the numerical simulation models of radial vibration of the transducer are established, and the amplitude-frequency characteristic curves and the displacement responses of radial vibration of the transducer are given. The simulation results show that the amplitude-frequency characteristic curves and the displacement responses of the two methods are the same, which verifies the correctness of the simulation results.(3)
The relationship between the natural frequency and the transducer’s structural sizes are obtained: when the thickness of the metal prestressed tube and the piezoelectric ceramic tube are constant, the natural frequency decreases with the increase of the inner diameter of the piezoelectric ceramic tube. And when the outer diameter of the metal prestressed tube and the inner diameter of the piezoelectric ceramic tube are constant, the natural frequency decreases with the increase of the thickness-to-wall ratio.(4)
The theory in this paper is based on the vibration theory of the mechanical system and combined with the constitutive equation. This method can be used to analyze the radial vibration characteristics of the transducer composed of a piezoelectric ceramic ring (tube) and an outer metal ring (tube). Compared with the current analysis method of radial vibration characteristics of the transducer based on the equivalent circuit, the calculation method of natural frequency based on elastic vibration theory is clear in concept and simple in calculation, and the simulation model can analyze the mechanical vibration of the transducer. This method has a guiding significance for the design of the piezoelectric ceramic composite ultrasonic transducer.



---

*Source: 9834324-2020-07-16.xml*