A Frequency-Dependent Dynamic Electric–Mechanical Network for Thin-Wafer Piezoelectric Transducers Polarized in the Thickness Direction: Physical Model and Experimental Confirmation

This paper is concerned with electric–acoustic/acoustic–electric conversions of thin-wafer piezoelectric transducers polarized in the thickness direction. By introducing two mechanical components with frequency-dependent values, i.e., radiation resistance and radiation mass, into the equivalent circuit of the thin-wafer piezoelectric transducer, we established a frequency-dependent dynamic mechanic-electric equivalent network with four terminals for an arbitrary given frequency, an enhancement from the conventional circuit networks. We derived the analytic expressions of its electric–acoustic and acoustic–electric conversion impulse responses using the four-terminal equivalent circuit to replace the traditional six-terminal equivalent circuit for a thin-wafer transducer with harmonic vibrational motion. For multifrequency electrical/acoustic signals acting on the transducer, we established parallel electric–acoustic/acoustic–electric conversion transmission networks. These two transmission network models have simple structures and clear physical and mathematical descriptions of thin-wafer transducers for electric–acoustic/acoustic–electric conversion when excited by a multifrequency electric/acoustic signal wavelet. The calculated results showed that the transducer’s center frequency shift relates to its mechanical load and vibration state. The method reported in this paper can be applied to conventional-sized and small-sized piezoelectric transducers with universal applicability.


Introduction
Vibration is the source of waves; a wave is a vibration state propagating through a medium.The two are closely connected and are the basis of acoustics.Thus, acoustic measurements are essential for measuring particle vibration states inside a medium.
Acoustic measurements are ubiquitous in the industry and in daily life.As an essential component in acoustic measurements, a transducer is fundamental for achieving the mutual conversion of acoustic and electric signals.A wide variety of piezoelectric transducers that are conventionally sized and small-sized are applied in many different fields, and some accompanying acoustic phenomena are investigated.In medical and biomedical fields, scientists distinguish normal and pathological issues to prevent, discover, and treat patients carrying some diseases by analyzing the interaction between human tissues and acoustic signals measured by small-sized transducers [1,2].Liu et al. used an arrayed flexible piezoelectric micromachined ultrasonic transducer with a sandwich structure for adjuvant treatment of bone injury [3].Rong et al. applied a piezoelectric micromachined ultrasonic transducer to wireless power transfer for implantable biomedical devices [4].By observing and analyzing acoustic signals measured by the small-sized transducer, Hong et al. found that acoustic stress and wave resonance play a crucial role in plasma bubbles, and relevant acoustic studies offer new ways to achieve sustainable chemistry [5].Ogawa pointed out that in measuring biological properties, acoustic effects play an essential role in designing miniaturized electronic instruments [6].Pop demonstrated a high-data-rate intrabody communication link based on Lithium Niobate piezoelectric micromachined ultrasonic transducers [7].For fluid acoustic sensors, some conventional-sized and smallsized piezoelectric transducers are also widely investigated for efficiently and precisely detecting the fluid's physical properties, e.g., impurities in oil, toxic substances in sewage, viscidity variation of blood, and so on [8][9][10][11].Zhao also researched the frequency shift induced by fluid viscosity for the frequency property of the measured acoustic signal [12].Wang et al. performed a study on the high-frequency vibration of piezoelectric plates excited by lateral electric fields produced by surface electrodes under viscous liquid loadings for sensing [13].For the industrial application of nondestructive detection, the acoustic signals measured by conventional-sized and small-sized transducers can provide valuable information for detecting the level of liquid carbon dioxide in containers or the oil level in a pipeline [14][15][16][17] and the quality of mechanical components [18][19][20][21][22].The vibration responses of the objects (such as bridges, aero engines, and so on), which are measured using conventional-sized and small-sized transducers, may provide data for safety inspections, early warning against structural damage, and any types of anomalies [23,24].For the fields of environmental protection and geological disaster prediction, using a conventional-sized transducer to perform acoustic measurements of buildings and surrounding environments can help improve hearing conditions [25][26][27].Ambient noise measured by the conventionalsized transducer can be used to invert some valuable information, such as the inversion of underground rock dynamic velocity structure from imaging background noise for earthquake prediction [28][29][30]; moreover, because ocean noise carries some valuable marine environmental information, its acoustic inversion is essential for ocean research [31][32][33].Conventional-sized piezoelectric transducers have also been widely used in petroleum logging and seismic exploration instruments [34][35][36].In underwater acoustics applications, larger-sized transducers are used to detect acoustic signals propagating over a long distance for real-time tracking and the high-precision positioning of underwater targets [37][38][39].
For different purposes, either application or scientific exploration, scientists have built physical models of many different transducers with various structures, shapes, polarization modes, etc., to study their physical properties using distinct methods.Wang et al. studied interdigital transducers and provided accurate modeling of the piezoelectric effect of coupled structures [40][41][42].Zhao et al. performed analyses on the lateral edge effect and segmented boundary conditions at a solid-liquid interface of the transducer [43].Hara et al. proposed a structural system identification method that uses semiactive inputs generated from piezoelectric transducers [44].Hur studied planar ultrasonic transducers based on a metasurface piezoelectric ring array for subwavelength acoustic focusing [45].Antonino et al. developed and investigated a particular geometry of transducers, emulating the shape of bats' cochlea, to transmit and receive ultrasounds in the air [46].Meanwhile, Yan et al. also studied piezoelectric material's electric-mechanical conversion coefficient to improve the piezoelectric transducer's electric-acoustic conversion efficiency [47].To generate precise and guided ultrasonic wave pulses for nondestructive testing purposes, Upadhyay and Schaal conducted research on an experimental model for commonly used guided ultrasound transducers [48].Sherrit et al. compared Mason and KLM's equivalent circuits with six terminals for piezoelectric resonators in the thickness mode [49].Mel'kanovich derived electric-acoustic six-terminal network equations for the transducer [50], and Sherrit and Konovalov provided a method for all transfer functions of a piezoelectric transducer with allowance for its electrical load [51].
The transient response of an object, e.g., a transducer, provides information on its physical property and measurement quality.Theoretical and experimental studies reported by Piquette showed that, under the excitation of a sinusoidal voltage signal, the dynamical state on the surface of a piezoelectric transducer went through a transient transition process from a static state to a steady state sinusoidal vibration [52,53].Moon et al. found that the resonant frequency of the piezoelectric transducer changed with the load's fluctuations and environmental conditions and proposed a method to track the resonant frequency using the transducer's transient characteristics [54].Wang et al. investigated the relationship between eigenfrequency and the early strength development of cement mortar by using a piezoelectric transducer [55].Their study indicated that the eigenfrequency is related to the mechanical loading (i.e., the property of the cement mortar around the transducer) and is an integrated response from the piezoelectric transducer's electric-acoustic (or acoustic-electric) conversions.Redwood qualitatively described the transient properties of piezoelectric transducers: "In the generation of ultrasonic waves a transient electrical signal is applied to the transducer and this produces a transient mechanical vibration, while in detection, the application of a transient mechanical signal produces an electrical vibration" [56].On the basis of Fa and Zhao's works [57][58][59], this paper performs a study on the physical characteristics of a thin-wafer piezoelectric transducer.The studies reported in this paper continue the earlier work with new developments and provide a mathematical description of the physical mechanism of Redwood's research [56].
In studies on isotropic solids, Fa and Zhao et al. predicted that the frequency spectrum corresponding to the temporary progression of dynamical oscillations is closely related to the intrinsic noise of the transmission medium, which is related to the internal structure and physical properties of the medium, and the maximum of this frequency spectrum is determined by the solid medium's resonant frequency and the property of external acoustic disturbance [60].
With our understanding of the physical mechanism of electric-acoustic (acousticelectric) conversion, this paper reports a frequency-dependent dynamic electric-mechanical network for a single thin-wafer piezoelectric transducer polarized in the thickness direction by introducing radiation resistance and radiation mass [61].The acoustic measurement goes through three processes: (i) The source transducer (referred to as the source) converts the electric-driving signal into an outward-radiated acoustic signal, and this electric-acoustic conversion process causes the vibration of particles inside the source to generate intrinsic noise; (ii) the radiated acoustic wave causes the vibration of particles inside the medium, and the interaction between the particles inside the transmission medium leads to the propagation of the acoustic wave in the medium [62], so it would also generate the intrinsic noise inside the medium [60]; (iii) the acoustic wave propagating in the medium is converted into an electric signal by the receiving transducer (referred to as the receiver), and the intrinsic noise generated by the acoustic-electric conversion will also appear.Therefore, the measured acoustic signal, i.e., the electric signal at the electric terminals of the receiver, contains intrinsic noise from the transient transition of the electric-acoustic conversion of the source: that of the vibration of particle inside the medium and that of the acousticelectric conversion of the receiver.Specifically, to describe the physical process of acoustic measurements clearly, we established a parallel transmission network consisting of many mechanical-electric equivalent circuits for the electric-acoustic/acoustic-electric conversions of the thin-wafer transducer under the excitation of multifrequency electric/acoustic signal wavelets.The impulse responses and system functions of electric-acoustic and acoustic-electric exchanges were derived from the residue theorem for a given sinusoidal frequency.Furthermore, a new method was proposed to extract the intrinsic noise generated by the propagation medium from the integrated frequency response for transitioning from the source to the receiver.
This study provides a theoretical basis for applying intrinsic noise to acoustic measurements.Based on the analysis of the frequency response characteristics under the excitation of multifrequency-driving signal wavelet, we discuss the physical mechanism of intrinsic noise generated by the electric-acoustic and acoustic-electric conversions of the transducer.Because any object with mass has inertia, there are transient transition processes for electric-acoustic and acoustic-electric conversions as well as the vibrations of particles inside the media.The signals obtained from acoustic measurements contain three kinds of intrinsic noises.The electric-acoustic (and acoustic-electric) conversion properties are vital to extracting noise generated inside the medium for the accurate inversion of its physical properties and internal structure.
Studying the transient characteristics of piezoelectric transducers as a general transducer or as a fluid acoustic sensor is essential.One potential application of the research reported in this paper is that by comparing the difference between the integrated frequency responses of source-medium-receiver and source-receiver processes, the physical properties and internal structure of the medium should be inverted accurately from the measured acoustic signals.

Background Information and New Physical Model
Let us consider a thin-wafer piezoelectric transducer in a stationary state with a given radius (r 0 ) and thickness (l t ), and it is polarized along the thickness direction.The two surfaces of the wafer disc of the transducer comprise two electrodes, as shown in Figure 1.Sur0, Sur1, and Sur2 refer to the side, rear, and front surfaces, respectively.
sponse for transitioning from the source to the receiver.
This study provides a theoretical basis for applying intrinsic noise to acoustic measurements.Based on the analysis of the frequency response characteristics under the excitation of multifrequency-driving signal wavelet, we discuss the physical mechanism of intrinsic noise generated by the electric-acoustic and acoustic-electric conversions of the transducer.Because any object with mass has inertia, there are transient transition processes for electric-acoustic and acoustic-electric conversions as well as the vibrations of particles inside the media.The signals obtained from acoustic measurements contain three kinds of intrinsic noises.The electric-acoustic (and acoustic-electric) conversion properties are vital to extracting noise generated inside the medium for the accurate inversion of its physical properties and internal structure.
Studying the transient characteristics of piezoelectric transducers as a general transducer or as a fluid acoustic sensor is essential.One potential application of the research reported in this paper is that by comparing the difference between the integrated frequency responses of source-medium-receiver and source-receiver processes, the physical properties and internal structure of the medium should be inverted accurately from the measured acoustic signals.

Background Information and New Physical Model
Let us consider a thin-wafer piezoelectric transducer in a stationary state with a given radius ( 0 r ) and thickness ( t l ), and it is polarized along the thickness direction.The two surfaces of the wafer disc of the transducer comprise two electrodes, as shown in Figure 1.Sur0, Sur1, and Sur2 refer to the side, rear, and front surfaces, respectively.Because the thickness of the transducer is much thinner than its radius and the intrinsic vibration frequency in the thickness direction is much higher than that in the radial direction, a simplified electric field applied to the transducer in the thickness direction may be helpful.Such a transducer is in a clipped state, with no time to deform in the radial direction, and only the strain component (S3) in the Z direction on the Z surface is nonzero.Let us label the Cartesian coordinates (X, Y, Z) as (1, 2, 3).In that case, we may label the component of electric field strength in the Z-axis (thickness) direction ( 3 E ) and the electric displacement vector in the Z-axis direction ( 3 D ).Now, the electric field applied to the transducer is in the same direction as the polarization, and its longitudinal piezoelectric equations are Because the thickness of the transducer is much thinner than its radius and the intrinsic vibration frequency in the thickness direction is much higher than that in the radial direction, a simplified electric field applied to the transducer in the thickness direction may be helpful.Such a transducer is in a clipped state, with no time to deform in the radial direction, and only the strain component (S 3 ) in the Z direction on the Z surface is nonzero.
Let us label the Cartesian coordinates (X, Y, Z) as (1, 2, 3).In that case, we may label the component of electric field strength in the Z-axis (thickness) direction (E 3 ) and the electric displacement vector in the Z-axis direction (D 3 ).Now, the electric field applied to the transducer is in the same direction as the polarization, and its longitudinal piezoelectric equations are Subscript "3" indicates the component of a physical quantity in the Z direction, such as the stress (T 3 = T ZZ ) and strain (S 3 = S ZZ ) and the elastic matrix element under constant electric displacement (C D 33 = C D ZZZZ ).Subscript "33" means "ZZ" for the dielectric isolation ratio matrix element under constant strain (β S 33 = β S ZZ ) and "ZZZ" for the piezoelectric stiffness matrix element (h 33 = h ZZZ ), where the first "3" (in h 33 ) stands for "Z" and the second "3" stands for "ZZ".Suppose the insulation performance of a transducer is ideal (no defects nor malfunctions).Then, only the Z-component of the electric displacement exists (D 3 = 0 and ∂D 3 /∂Z = 0), other components will be zeroes (D 1 = D 2 = 0, i.e., D X = D Y = 0), and there is no free charge inside the transducer.For a given density of piezoelectric material (ρ) and the particle displacement component (u Z ), we obtain the vibrational equation of motion along the thickness (Z-axis) direction as Multiplying both sides of Equation ( 3) by the volume of the transducer (V = πr 0 2 l t ) and noting the mass of the transducer (m), we obtain By placing the transducer in coupling liquid and vibrating in the Z-direction, the surrounding medium expands and contracts alternately, and the transducer radiates acoustic waves outward, which will generate a reaction force that is exerted on the transducer.
In Equation ( 5), we have several related parameters, namely, radiation resistance , the radiation reactance of the transducer (X r = ρ 0 c 0 [kr 0 /(1 + k 2 r 0 2 )]S 0 ), the density of the coupling medium (ρ 0 ), the acoustic velocity (c 0 ), and wave number (k = ω/c 0 ) in the coupling medium.Because the radius of the wafer is much larger than its thickness (r 0 >> l t ) and its side and rear surfaces, i.e., Sur0 and Sur1 are in the clamped state, we only need to consider the contribution of the front surface area to radiation resistance (R r ) and radiation reactance (X r ), which means that we choose to use an effective surface (S 0 ≈ πr 0 2 ).Due to the viscosity of the piezoelectric material, there is frictional force resistance (R m ) for the particle inside the transducer during its vibration, leading to the corresponding frictional force.
The magnitude of the frictional force is related to the viscosity of the piezoelectric material, the contact area of the vibrating surface of the transducer with the coupling fluid, and the physical property of the coupling fluid.
Overall, the resultant external (electric or mechanical) disturbance acting on the transducer when vibrating in coupling liquid is During vibrations, the transducer is affected by both external forces and self-generated internal stress.Then, we have a complete equation of motion instead of Equation (4).
Meanwhile, from Equations ( 1) and ( 2), we have Micromachines 2023, 14, 1641 6 of 27 By the substitution of Equation ( 9) into Equation ( 8), we obtain Combined with actual acoustic sensors, we may set the rear surface (Sur1) of the wafer transducer so that it is tightly cemented to the backing layer (i.e., this surface is considered stationary), and the front surface side (Sur2) is coated with a skinny layer of insulating material; i.e., the front surface is considered approximately free and radiates acoustic waves outward).Because the thickness of the thin-wafer transducer is much smaller than its radius, the edge effect on the establishment of the machine-electric transmission network model is negligible, and we can apply the compression/tension deformation of a uniform rod to the dynamic state of the transducer, as shown in Figure 2.
  Meanwhile, from Equations ( 1) and (2), we have By the substitution of Equation ( 9) into Equation ( 8), we obtain Combined with actual acoustic sensors, we may set the rear surface (Sur1) of the wafer transducer so that it is tightly cemented to the backing layer (i.e., this surface is considered stationary), and the front surface side (Sur2) is coated with a skinny layer of insulating material; i.e., the front surface is considered approximately free and radiates acoustic waves outward).Because the thickness of the thin-wafer transducer is much smaller than its radius, the edge effect on the establishment of the machine-electric transmission network model is negligible, and we can apply the compression/tension deformation of a uni- form rod to the dynamic state of the transducer, as shown in Figure 2. The micro-volume element (or particle) displacement in the thickness direction (i.e., Z-axis direction in Figure 2) inside the thin-wafer piezoelectric transducer can be written as follows: where l is the thickness of the transducer during its vibration process, and Z is the position of the micro-volume element inside the transducer.The relation of the strain versus particle displacement The micro-volume element (or particle) displacement in the thickness direction (i.e., Z-axis direction in Figure 2) inside the thin-wafer piezoelectric transducer can be written as follows: where l is the thickness of the transducer during its vibration process, and Z is the position of the micro-volume element inside the transducer.The relation of the strain versus particle displacement leads to We only need to consider the vibration state of the radiation surface, Sur2, and the displacement (l − l t ) created by this surface; thus, in the range inside the transducer near Sur2, Formula (13) can be rewritten as By substituting Equation ( 14) into Equation (10), the state equation of the transducer vibrating in the coupling liquid can be obtained.
Under harmonic vibrations along the thickness direction with the amplitude, u 0 , of the response, the temporal particle displacement of the transducer's radiation surface Sur2 would be The substitution of Equation ( 16) into Equation ( 15) yields the following expression.
) ) depends only on the intrinsic properties of the transducer.Nevertheless, radiation resistance (R r ) and radiation quality (m r = X r /ω) are related to the geometric size of the transducer, the physical property of the coupling medium, and the frequency of the sinusoidal signal (electric or acoustic) that excites the transducer.
Based on Gauss's theorem, the total charge on each electrode is as follows [63].
A combination of Equations ( 2) and ( 14) yields The instantaneous current comes from the first-order derivation of the accumulated charge on each electrode surface concerning time.The substitution of Equations ( 17) and (19) into Equation (18) yields the instantaneous current.
Equation ( 20) is also known as the mechanical-electric state equation of the transducer, where several parameters are the static capacitance of the transducer (C 0 = πr 0 2 /β S 33 l t ), the voltage applied between the two electrodes (V = E 3 l t ), and the mechanical-electric conversion coefficient (φ = πr 0 2 h 33 /β S 33 l t ).We can establish the mechanical-electric equivalent circuits of the electric-acoustic and acoustic-electric conversions of the transducer based on Equation (20), as shown in Figure 3.The defined physical quantities in Figure 3 are the sinusoidal driving-voltage signal (U 1 (t)), the output resistance of the driving-voltage source (R o ), the voltage signal at the electric terminal of the receiver (U 3 (t)), the input resistance of the measurement circuit (R i ), and the particle displacement velocity on the surface of the transducer (v(t)).
We must pay attention to substantial differences in the characteristics and the analysis methods between conventional circuits and the mechanical-electric equivalent circuit of the transducer.The values of all components, e.g., resistance and inductance, etc., are independent of frequency in conventional circuits.In contrast, some mechanical elements in the mechanical-electric equivalent electric circuits of the transducer are functions of frequency, e.g., radiation resistance and radiation quality, as shown in Figure 3.
The mechanical-electric equivalent circuits in the time domain, as shown in Figure 3, were transformed into the s domain, as shown in Figure 4. From there, we analyze the inherent characteristics of the transducer and solve its electric-acoustic (and acoustic-electric) impulse responses and the corresponding electric-acoustic (and acoustic-electric) system functions.
of the transducer.The values of all components, e.g., resistance and inductance, etc., are independent of frequency in conventional circuits.In contrast, some mechanical elements in the mechanical-electric equivalent electric circuits of the transducer are functions of frequency, e.g., radiation resistance and radiation quality, as shown in Figure 3.
The mechanical-electric equivalent circuits in the time domain, as shown in Figure 3, were transformed into the s domain, as shown in Figure 4. From there, we analyze the inherent characteristics of the transducer and solve its electric-acoustic (and acoustic-electric) impulse responses and the corresponding electric-acoustic (and acoustic-electric) system functions.From the equivalent mechanical-electric circuit of the source, as shown in Figure 4a, and based on Kirchhoff's law, we have sis methods between conventional circuits and the mechanical-electric equivalent circuit of the transducer.The values of all components, e.g., resistance and inductance, etc., are independent of frequency in conventional circuits.In contrast, some mechanical elements in the mechanical-electric equivalent electric circuits of the transducer are functions of frequency, e.g., radiation resistance and radiation quality, as shown in Figure 3.
The mechanical-electric equivalent circuits in the time domain, as shown in Figure 3, were transformed into the s domain, as shown in Figure 4. From there, we analyze the inherent characteristics of the transducer and solve its electric-acoustic (and acoustic-electric) impulse responses and the corresponding electric-acoustic (and acoustic-electric) system functions.From the equivalent mechanical-electric circuit of the source, as shown in Figure 4a, and based on Kirchhoff's law, we have From the equivalent mechanical-electric circuit of the source, as shown in Figure 4a, and based on Kirchhoff's law, we have Since the initial state of the transducer is at equilibrium, the initial energy storage of the relevant electric and mechanical (or dynamical) components (such as capacitance, mass, radiation mass, compliance, etc.) in the mechanical-electric circuit is zero.Any response would have to be from external excitation, and the processes of electric-acoustic (and acoustic-electric) conversions would define the zero-state response.Therefore, in the s domain, the system function of electric-acoustic conversions is a simple ratio of the surface vibration velocity (v(s)) of the transducer to the driving voltage (U 1 (s)).And from Equations ( 21) and ( 22), we have where Let us define a few parameters to simplify the discussion and mathematical derivation: Then, we can write Equation ( 23) as Its discriminant (D = p 3 + q 2 = (b/3 + a 2 /9) 3 + (a 3 /27 − ab/6 + c/2) 2 ) controls the mode of dynamical motion.
Based on the residue theorem, we would have the following cases of the impulse response for the electric-acoustic conversion: (i) The discriminant has a negative value (D < 0), and system function (H 1 (s)) has three non-equal real poles.
The relevant parameters are α , and y = (−q − √ D) 1/3 .In this case, employing a unit step function (ε(t)), the impulse response from the electric-acoustic conversion of the source is This presents a state of over damping, where The discriminant is zero (D = 0), and the system function (H 1 (s)) has three real poles, two of which are identical.
The electric-acoustic conversion impulse response is where 2 .Equation (32) shows that the transducer is in a critical damping state.
(iii) The discriminant has a positive value (D > 0), and the system function (H 1 (s)) has a real-valued pole and a pair of complex conjugate poles.
The impulse response of the electric-acoustic conversion is in an oscillatory mode with an oscillation frequency, f 1s = ω 1s /2π, which can be written as ) where ).In this case, the source is in an under-damped state, and the radiated acoustic waves proceed outwards under the excitation of the driving-voltage signal, in line with the actual physical condition.The source's oscillation frequency is related to the radiation resistance and mass, so it is an important physical quantity for describing the characteristics of the transient process generated by the transducer under the excitation of any sinusoidal voltage signal with a given frequency.
Since Equation ( 36) is integrable with an existing relation (s = iω), the electric-acoustic conversion system function of a source, shown in Equation ( 24), yields Similarly, acoustic-electric conversion is the reverse of electric-acoustic transformation for a transducer.The system function of acoustic-electric conversion for the receiver is the ratio of voltages between its two electrodes to the particle displacement velocity on its surface in the frequency domain.
Since the transducer works in the oscillation mode, we only need to discuss the underdamped state of the receiver.Similarly to the source's electric-acoustic conversion, we can obtain the acoustic-electric impulse response, the system function of the acoustic-electric conversion, and the oscillation angular frequency of the receiver. ) Here, we have ), and θ = arctan(N /M ).

Calculation and Analysis
In the following discussion and calculations, we selected water as the coupling liquid around the transducer with the known acoustic velocity (c 0 = 1500 m/s) and density (ρ 0 = 1000 kg/m 3 ).Table 1 displays the physical parameters of the piezoelectric material PZT-5H used for the transducer [64].Symbol (ε 0 ) is the dielectric constant in a vacuum.

Structure Relationships of Oscillation and Center Frequencies versus Geometric Size
We selected the output resistance (R o ) of the driving-voltage signal source as 50 Ω, the input resistance (R i ) of the measuring circuit as 200 Ω, the free-load transducer is in a vacuum (R r = m r = R m = 0), and the loaded transducer is in water (R m = 0.2 × πr 2 b ρ m ).We also defined the transducer's center frequency as the frequency corresponding to its system function's maximum modulus.For the case of free-load, the center frequency (a source or a receiver) is irrelevant to the sinusoidal frequency of the excitation (electric or acoustic) exerted on the transducer, and it is only determined by the transducer's intrinsic physical and geometrical properties due to R r = 0 and m r = 0.The loading center frequency relies on the transducer's physical and geometrical parameters, the coupling medium's physical nature, and the sinusoidal frequency of the excitation signal (either electric or acoustic).
Taking the radius (r 0 ) of the transducer to be 8 mm as an example, Figure 5a displays the free-load center, load center, and free-load oscillation frequencies of the transducer versus its thickness.Figure 5a shows that regardless of whether the transducer is a source or a receiver, (i) the above three kinds of frequencies decrease significantly with increasing thickness.(ii) The free-load center, load center, and oscillation frequencies of the receiver are slightly higher than those of the source.(iii) The free-load center frequency is slightly higher than the load center frequency.(iv) The oscillation frequency can be greater than the free-load center frequency within a radius range, and it can also be lower than the latter in another radius range; this is because the effect of the radius on the direct current of the impulse response (electric-acoustic or acoustic-electric) is smaller than that of the oscillation frequency for the first case, and the effect of the radius on the direct current of the impulse response is more significant than the oscillation frequency for the second case.At l t = 2 mm and r 0 = 8 mm, f 10 = 1.001409MHz, f 1d = 1.000575MHz, f 1s0 = 1.002307MHz, f 30 = 1.006017MHz, f 3d = 1.005189MHz, and Similarly, in selecting the thickness of the transducer to be 2 mm, Figure 5b displays the free-load center, load center, and free-load oscillation frequencies of the transducer versus its radius.Figure 5b shows: (i) as a source, the free-load center, load center, and free-load oscillation frequencies of the transducer increase slightly with the radius; (ii) as a receiver, the above three kinds of frequencies of the transducer increase very slowly with increasing its radius.
Figure 6 is the calculated relationships of load oscillation frequencies ( f 1sd and f 3sd ) of both source and receiver versus the sinusoidal frequency (f ) of the excitation signal exerted on the transducer (electric or acoustic).It displays f 1sd and f 3sd has a monotonic increasing tendency with both the frequency of the sinusoidal driving-voltage signal and that of the acoustic signal arriving at the receiver.It also shows that the sinusoidal signal frequency has a more significant effect on the load oscillation frequency in the lower frequency range and has a minimal impact on the load oscillation frequency in the higher frequency range.The loading oscillation frequency of the transducer (as a source or as a receiver) approaches its free-load center frequency as the frequency of the applied sinusoidal signal increases.The load oscillation frequency of the source signal is slightly smaller than that of the receiver.

Electric-Acoustic and Acoustic-Electric Impulse Responses and Corresponding System Functions
The piezoelectric physical nature and geometric parameters of the transducers guarantee the system is always in an oscillatory mode in practical applications.Figures 7 and 8 show the calculated electric-acoustic (and acoustic-electric) impulse responses and the corresponding amplitude spectra for the piezoelectric transducers (both free-load and load), which are the same as discussed in the Section 5.  f and 3sd f has a monotonic increasing tendency with both the frequency of the sinusoidal driving-voltage signal and that of the acoustic signal arriving at the receiver.It also shows that the sinusoidal signal frequency has a more significant effect on the load oscillation frequency in the lower frequency range and has a minimal impact on the load oscillation frequency in the higher frequency range.The loading oscillation frequency of the transducer (as a source or as a receiver) approaches its free-load center frequency as the frequency of the applied sinusoidal signal increases.The load oscillation frequency of the source signal is slightly smaller than that of the receiver.f and 3sd f has a monotonic increasing tendency with both the frequency of the sinusoidal driving-voltage signal and that of the acoustic signal arriving at the receiver.It also shows that the sinusoidal signal frequency has a more significant effect on the load oscillation frequency in the lower frequency range and has a minimal impact on the load oscillation frequency in the higher frequency range.The loading oscillation frequency of the transducer (as a source or as a receiver) approaches its free-load center frequency as the frequency of the applied sinusoidal signal increases.The load oscillation frequency of the source signal is slightly smaller than that of the receiver.The radiation resistance and mass in the equivalent circuits of Figure 3 are frequencydependent.Figures 7 and 8 show the calculated electric-acoustic (and acoustic-electric) impulse responses and corresponding amplitude spectra at four selected frequencies with radiation resistance and radiation mass.Interestingly, the electric-acoustic (and acousticelectric) conversion impulse responses, the corresponding amplitude spectra, and the center frequencies are unvaried with respect to frequency for a free-load transducer but vary with the frequency for loaded transducers.It is plausible that the radiation resistance and mass of the transducer are indeed functions of frequency.The frequency spectrum corresponding to the transient process mainly reflects the inherent physical characteristics of the transducer, even though the sinusoidal signal (electric or acoustic) exerted on the transducer has some influence on it.The frequency corresponding to the amplitude spectrum's maximum stands for the central frequency of the transducer.

Electric-Acoustic and Acoustic-Electric Impulse Responses and Corresponding System Functions
and the center frequencies are unvaried with respect to frequency for a free-load transducer but vary with the frequency for loaded transducers.It is plausible that the radiation resistance and mass of the transducer are indeed functions of frequency.The frequency spectrum corresponding to the transient process mainly reflects the inherent physical characteristics of the transducer, even though the sinusoidal signal (electric or acoustic) exerted on the transducer has some influence on it.The frequency corresponding to the amplitude spectrum's maximum stands for the central frequency of the transducer.From Figures 7 and 8, we observed the following: (i) The center frequencies corresponding to the electric-acoustic and acoustic-electric conversions of the transducer increase with the sinusoidal frequency and tend toward the corresponding free-load centerfrequency, but the magnitudes of the corresponding amplitude spectra decrease with increasing frequency, namely the transducer's electric-acoustic and acoustic-electric conversion efficiencies decrease with increasing frequency.(ii) For the same transducer, although its acoustic-electric conversion is the inverse process of the electric-acoustic conversion and the difference between the center frequency corresponding to electric-acoustic conversion and that corresponding to acoustic-electric conversion is also insignificant, the characteristics of its two kinds of conversions are different.For example, the bandwidth of the amplitude spectrum of the electric-acoustic conversion is greater than that of its acoustic-electric conversion; that is to say that compared with the filtering effect of the transducer on the frequency component of the driving-voltage signal far from its center frequency, the acoustic-electric conversion has a more substantial filtering effect on the frequency component of the received acoustic signal far from its center frequency.We have normalized acoustic-electric impulse responses and amplitude spectra using the maximum value of the impulse response and amplitude spectrum for the case of free load, respectively.The cyan curve is the acoustic-electric impulse response and the corresponding amplitude spectrum for the free load.The other colors are those with mechanical load: the lines in magenta, blue, red, and black are the transducer acoustic-electric impulse responses, and the amplitude spectra of the corresponding system functions at selected frequencies of 0.1 3d f , 0.5 3d f , 1.0 3d f , and 10 3d f .Here, 3d f = 1.005189MHz.

Acoustic Signals Radiated by Transducers under the Excitation of Sinusoidal Driving-Voltage Signals with Different Frequencies
Four sinusoidal driving-voltage signals with an amplitude of 1 V are selected, and their four frequencies are selected as 0.1 1d f , 0. 5 1d f , 1 1d f , and 1.3 1d f to excite the source ( 1d f = 1.000575MHz). Figure 9 shows the calculated vibration velocity of the source's sur- face.We have normalized acousticelectric impulse responses and amplitude spectra using the maximum value of the impulse response and amplitude spectrum for the case of free load, respectively.The cyan curve is the acoustic-electric impulse response and the corresponding amplitude spectrum for the free load.The other colors are those with mechanical load: the lines in magenta, blue, red, and black are the transducer acousticelectric impulse responses, and the amplitude spectra of the corresponding system functions at selected frequencies of 0.1 f 3d , 0.5 f 3d , 1.0 f 3d , and 10 f 3d .Here, f 3d = 1.005189MHz.
From Figures 7 and 8, we observed the following: (i) The center frequencies corresponding to the electric-acoustic and acoustic-electric conversions of the transducer increase with the sinusoidal frequency and tend toward the corresponding free-load center-frequency, but the magnitudes of the corresponding amplitude spectra decrease with increasing frequency, namely the transducer's electric-acoustic and acoustic-electric conversion efficiencies decrease with increasing frequency.(ii) For the same transducer, although its acoustic-electric conversion is the inverse process of the electric-acoustic conversion and the difference between the center frequency corresponding to electric-acoustic conversion and that corresponding to acoustic-electric conversion is also insignificant, the characteristics of its two kinds of conversions are different.For example, the bandwidth of the amplitude spectrum of the electric-acoustic conversion is greater than that of its acousticelectric conversion; that is to say that compared with the filtering effect of the transducer on the frequency component of the driving-voltage signal far from its center frequency, the acoustic-electric conversion has a more substantial filtering effect on the frequency component of the received acoustic signal far from its center frequency.

Acoustic Signals Radiated by Transducers under the Excitation of Sinusoidal Driving-Voltage Signals with Different Frequencies
Four sinusoidal driving-voltage signals with an amplitude of 1 V are selected, and their four frequencies are selected as 0.1 f 1d , 0.5 f 1d , 1 f 1d , and 1.3 f 1d to excite the source ( f 1d = 1.000575MHz). Figure 9 shows the calculated vibration velocity of the source's surface.
The calculated results show that when excited by any sinusoidal driving-voltage signal with a given frequency, the transducer will go through a transient process from a static state to a steady sinusoidal vibration state due to the inertia of the particles with mass inside the transducer.When the frequency of the sinusoidal driving-voltage signal is equal to the load center frequency of the source, the source resonates with the sinusoidal driving-voltage signal, and the amplitude of its surface vibration is the largest.This also means that the amplitude of the acoustic wave signal radiated outward by the source is the largest (see Figure 9e).In this case, the frequency spectrum of the sinusoidal driving-voltage signal is stacked with the frequency spectrum of the electric-acoustic conversion precisely at the load center frequency of the source (see Figure 9f).When the frequency of the sinusoidal driving-voltage signal is not equal to (greater than or less than) the load center frequency of the source, the amplitude of the radiated acoustic signal decreases.We have normalized acoustic-electric impulse responses and amplitude spectra using the maximum value of the impulse response and amplitude spectrum for the case of free load, respectively.The cyan curve is the acoustic-electric impulse response and the corresponding amplitude spectrum for the free load.The other colors are those with mechanical load: the lines in magenta, blue, red, and black are the transducer acoustic-electric impulse responses, and the amplitude spectra of the corresponding system functions at selected frequencies of 0.1 3d f , 0.5 3d f , 1.0 3d f , and 10 3d f .Here, 3d f = 1.005189MHz.

Acoustic Signals Radiated by Transducers under the Excitation of Sinusoidal Driving-Voltage Signals with Different Frequencies
Four sinusoidal driving-voltage signals with an amplitude of 1 V are selected, and their four frequencies are selected as 0.1 1d f , 0. 5 1d f , 1 1d f , and 1.3 1d f to excite the source ( 1d f = 1.000575MHz). Figure 9 shows the calculated vibration velocity of the source's sur- face.The calculated results show that when excited by any sinusoidal driving-voltage signal with a given frequency, the transducer will go through a transient process from a static state to a steady sinusoidal vibration state due to the inertia of the particles with mass inside the transducer.When the frequency of the sinusoidal driving-voltage signal is equal to the load center frequency of the source, the source resonates with the sinusoidal driving-voltage signal, and the amplitude of its surface vibration is the largest.This also means that the amplitude of the acoustic wave signal radiated outward by the source is the largest (see Figure 9e).In this case, the frequency spectrum of the sinusoidal drivingvoltage signal is stacked with the frequency spectrum of the electric-acoustic conversion precisely at the load center frequency of the source (see Figure 9f).When the frequency of

The Electric Signal Is Converted from an Acoustic Signal by the Receiver Transducer
As mentioned in the earlier section, the acoustic-electric conversion of the transducer is an inverse process of the electric-acoustic transformation but differs in some characteristics.Below, we will analyze and discuss two cases.

Case 1
To analyze only the acoustic-electric conversion characteristics of the transducer, sinusoidal acoustic signals with an amplitude of 1 V and frequencies of 0.1 f 3d , 0.5 f 3d , f 3d , and 1.3 f 3d are directly input into the mechanical terminals of the receiver.Here, the central frequency is f 3d = 1.005189MHz. Figure 10 shows the voltage signal outputted by the electric terminals of the receiver after acoustic-electric conversion.
central frequency is 3d f = 1.005189MHz. Figure 10 shows the voltage signal outputted by the electric terminals of the receiver after acoustic-electric conversion.
The calculated results showed a transient process from a zero state to a steady sinusoidal voltage signal at the receiver's electric terminals.The temporary duration of the acoustic-electric conversion is longer than the electric-acoustic conversion because the amplitude spectrum's bandwidth corresponding to the transducer's acoustic-electric conversion is narrower than the corresponding electric-acoustic conversion.

Case 2
For the integrated response characteristics of the source-receiver, we set both the source and the receiver in an ideal liquid (viscosity and acoustic attenuation coefficients are zeroes) with an acoustic velocity of 1500 m per second, and the distance from the source transducer to the receiver is 0.15 m.
When excited by sinusoidal driving-voltage signals with an amplitude of 1 V and various frequencies (0.1 1d f , 0.5 1d f , 1.0 1d f , and 1.3 1d f ), the source radiates acoustic waves that propagate to the receiver without distortion and attenuation through the ideal liquid.
Figure 11 shows the voltage signals at the receiver's electric terminals.
Figure 11 shows the following: (i) Despite the fact that the bandwidth of the receiver's frequency response (amplitude spectrum) is narrower than the source's frequency response, the center frequency of the source is slightly lower than that of the receiver.(ii) When the frequency of the sinusoidal driving-voltage signal is close to the center frequency of the integrated frequency response of the source-receiver, the transient process of the electric signal at the receiver's electrical terminals experiences the lowest number of oscillation periods, and the amplitude spectrum of the sinusoidal driving-voltage signal is superimposed on the amplitude spectrum corresponding to the integrated frequency The calculated results showed a transient process from a zero state to a steady sinusoidal voltage signal at the receiver's electric terminals.The temporary duration of the acoustic-electric conversion is longer than the electric-acoustic conversion because the amplitude spectrum's bandwidth corresponding to the transducer's acoustic-electric conversion is narrower than the corresponding electric-acoustic conversion.

Case 2
For the integrated response characteristics of the source-receiver, we set both the source and the receiver in an ideal liquid (viscosity and acoustic attenuation coefficients are zeroes) with an acoustic velocity of 1500 m per second, and the distance from the source transducer to the receiver is 0.15 m.
When excited by sinusoidal driving-voltage signals with an amplitude of 1 V and various frequencies (0.1 f 1d , 0.5 f 1d , 1.0 f 1d , and 1.3 f 1d ), the source radiates acoustic waves that propagate to the receiver without distortion and attenuation through the ideal liquid.

Multifrequency Driving-Voltage Signal Exciting Transducer
To understand the comprehensive response characteristics of the source-receiver under the excitation of the wavelet of the driving-voltage signal with many different frequency components, we ignore the influence of the medium on the radiated sound wave.Assuming that the medium is ideal and elastic, the outer surfaces of the source and receiver transducer are tightly bonded to the medium, and the propagation process of the acoustic wave does not exhibit attenuation and dispersion.
In most acoustic measurements, both the driving-voltage signal exciting the source and the acoustic signal arriving at the receiver are wavelets containing many frequency components with different amplitudes and initial phases.Based on the Fourier transform, Figure 11 shows the following: (i) Despite the fact that the bandwidth of the receiver's frequency response (amplitude spectrum) is narrower than the source's frequency response, the center frequency of the source is slightly lower than that of the receiver.(ii) When the frequency of the sinusoidal driving-voltage signal is close to the center frequency of the integrated frequency response of the source-receiver, the transient process of the electric signal at the receiver's electrical terminals experiences the lowest number of oscillation periods, and the amplitude spectrum of the sinusoidal driving-voltage signal is superimposed on the amplitude spectrum corresponding to the integrated frequency response of the source-receiver.(iii) If the frequency of the sinusoidal driving-voltage signal is far from the center frequency of the integrated frequency response of the source-receiver, the number of periods of the electric signal at the receiver's electric terminals from the zero state to the stable sinusoidal vibration state increases, and the amplitude of the time domain waveform and the magnitude of amplitude spectrum decreases.

Multifrequency Driving-Voltage Signal Exciting Transducer
To understand the comprehensive response characteristics of the source-receiver under the excitation of the wavelet of the driving-voltage signal with many different frequency components, we ignore the influence of the medium on the radiated sound wave.Assuming that the medium is ideal and elastic, the outer surfaces of the source and receiver transducer are tightly bonded to the medium, and the propagation process of the acoustic wave does not exhibit attenuation and dispersion.
In most acoustic measurements, both the driving-voltage signal exciting the source and the acoustic signal arriving at the receiver are wavelets containing many frequency components with different amplitudes and initial phases.Based on the Fourier transform, an electric (or acoustic) signal used to excite a transducer can be a linear composition of sinusoidal components with different amplitudes, frequencies, and initial phases.Because radiation impedance and mass are frequency functions, each frequency component corresponds to a different electric-acoustic (or acoustic-electric) equivalent circuit within a circuit network.The electric-acoustic (and acoustic-electric) conversion can be achieved through a parallel transmission network consisting of many electric-acoustic (and acousticelectric) conversion circuits, as shown in Figure 12 (Parts I and III).Each electric circuit within the network would contribute its unique electric-acoustic (or acoustic-electric) impulse response due to its self-radiation resistance and mass.In this network, the related physical parameters include the j th frequency component (U 1j ) in the driving-voltage signal (U 1 ) (j = 1, 2, . ... . . N); the ee electric-acoustic impulse response (h 1j (t)) of the source corresponding to the j th sinusoidal frequency component in the driving-voltage signal (U 1 ); the particle displacement velocity (v 1j ) when excited by the sinusoidal component (U 1j ) in the driving-voltage signal; the particle displacement velocity (v 2 ) when excited by the driving-voltage signal (U 1 ); the impulse response (h 2 (t − ∆t)) of the medium, the propagation time (∆t) of the acoustic wave from the source to the receiver; the acoustic signal (v 3 ) reaching the receiving transducer; the k th frequency component (v 3k ) (k = 1, 2, . ... . . M) of a sound wave signal; the acoustic-electricity impulse response (h 3k (t − ∆t)) of the receiving transducer corresponding to the k th sinusoidal frequency (ω g ), and time window (t 0 ).We define the Heaviside function as (H(•)), and we use the signal to perform some calculations.
The driving-voltage signal in the frequency domain is For the driving-voltage signal at a set of given parameters (ω g = ω 1d , U 0 = 1 V, and t 0 = 6π/ω g = 6π/ω 1d ), Figure 13 shows the time-domain waveform, amplitude spectrum, and phase spectrum.The frequency corresponding to the maximum amplitude spectrum (0.9832 MHz) is slightly lower than that of the gate-selection sinusoidal signal (1.000575 MHz).
When writing the gate-selection sinusoidal signal into a series of sinusoidal components, each sinusoidal component is taken as an independent excitation source corresponding to the parallel equivalent circuit in Part I of Figure 12.
All frequency components that excite the source transducer would make it radiate acoustic signals outward, for which their accumulation forms the acoustic wavelet emitted from the source transducer, as shown in Figure 14e.The calculated results show that when each frequency component in the driving-voltage signal excites the source, the radiated acoustic signal quickly transitions from the zero states to the steady sinusoidal vibration.The corresponding amplitude spectrum shown in Figure 14f shows that the center frequency of the radiated acoustic signal wavelet is 0.997 MHz, which is slightly lower than the center frequency of the gate sinusoidal driving-voltage signal.When writing the gate-selection sinusoidal signal into a series of sinusoidal components, each sinusoidal component is taken as an independent excitation source corresponding to the parallel equivalent circuit in Part I of Figure 12.
All frequency components that excite the source transducer would make it radiate acoustic signals outward, for which their accumulation forms the acoustic wavelet emitted from the source transducer, as shown in Figure 14e.The calculated results show that when each frequency component in the driving-voltage signal excites the source, the radiated acoustic signal quickly transitions from the zero states to the steady sinusoidal vibration.The corresponding amplitude spectrum shown in Figure 14f shows that the center frequency of the radiated acoustic signal wavelet is 0.997 MHz, which is slightly lower than the center frequency of the gate sinusoidal driving-voltage signal.With respect to the integrated response characteristics of the source-receiver, Figure 15 shows the electric signal output at the receiver's electric terminals when the source is excited by a gate-selection sinusoidal driving-voltage signal.
The calculated results show that compared with the acoustic signals of the radiated acoustic signal when the source is excited by each frequency component in the gate sinusoidal drive-voltage signal, the duration of the transient transition process of the corresponding output electric signal at the electric terminals of the receiver (see each circuit of Part III of the parallel network shown in Figure 12) becomes more prolonged, and the time domain waveform becomes more disorderly, as shown in Figure 14a-d.acoustic signals outward, for which their accumulation forms the acoustic wavelet emitted from the source transducer, as shown in Figure 14e.The calculated results show that when each frequency component in the driving-voltage signal excites the source, the radiated acoustic signal quickly transitions from the zero states to the steady sinusoidal vibration.The corresponding amplitude spectrum shown in Figure 14f shows that the center frequency of the radiated acoustic signal wavelet is 0.997 MHz, which is slightly lower than the center frequency of the gate sinusoidal driving-voltage signal.With respect to the integrated response characteristics of the source-receiver, Figure 15 shows the electric signal output at the receiver's electric terminals when the source is excited by a gate-selection sinusoidal driving-voltage signal.a-d) display the convolutions for several frequency components (0.4 f 1d , 0.8 f 1d , 1.2 f 1d , and 1.6 f 1d ) and the corresponding electric-acoustic impulse response.Plot (e) shows the normalized accumulative output from the circuit network.And (f) is the normalized amplitude spectrum of the accumulative signal output.
The cumulative output of all parallel circuits in Part III in Figure 12, namely the electric signal outputted by the electric terminals at the receiver, is shown in Figure 15e.Compared with the acoustic signal wavelet radiated by the source, the oscillation duration of the electric signal at the receiver's electric terminals becomes more extended, and the bandwidth of the corresponding amplitude spectrum becomes narrower.In addition, the amplitude of the frequency component far away from the center frequency of the receiver decreases significantly or disappears.The center frequency of the electric signal at the receiver's electric terminals is 1.003 MHz, slightly higher than the acoustic signal radiated by the source due to the receiver's acoustic-electric filtering effect on the received acoustic signal.The calculated results show that compared with the acoustic signals of the radiated acoustic signal when the source is excited by each frequency component in the gate sinusoidal drive-voltage signal, the duration of the transient transition process of the corresponding output electric signal at the electric terminals of the receiver (see each circuit of Part III of the parallel network shown in Figure 12) becomes more prolonged, and the time domain waveform becomes more disorderly, as shown in Figure 14a-d.
The cumulative output of all parallel circuits in Part III in Figure 12, namely the electric signal outputted by the electric terminals at the receiver, is shown in Figure 15e.Compared with the acoustic signal wavelet radiated by the source, the oscillation duration of the electric signal at the receiver's electric terminals becomes more extended, and the bandwidth of the corresponding amplitude spectrum becomes narrower.In addition, the

Experimental Measurement
To confirm the quality of our theoretical prediction, we selected two thin-wafer piezoelectric transducers with a thickness of 2 mm and a radius of 8 mm for experimental measurements in this paper.One acts as the source transducer, and another one works as the receiver.Table 1 illustrates that the material constituting the transducer is the manufactured PZT-5H and provides the corresponding physical parameters.Figure 16 shows the physical image of the transducer.The system software based on the LabVIEW platform (Version of 2011) is used for measurement control, which consists of an electric-signal module, a display/storage module, and a slippage/rotation module.We placed the source and receiver in a tank filled with water.The virtual panel on a computer screen gave instructions to place the two transducer electrodes opposite of one another at 15 cm apart.The control panel selects the gate-selection sinusoidal drivingvoltage signal with an amplitude of 1 V, a gate width of 3 cycles, and frequencies of 0.3, 1.0, and 1.7 1d f to excite the source transducer.The measurement system collects and    The system software based on the LabVIEW platform (Version of 2011) is used for measurement control, which consists of an electric-signal module, a display/storage module, and a slippage/rotation module.We placed the source and receiver in a tank filled with water.The virtual panel on a computer screen gave instructions to place the two transducer electrodes opposite of one another at 15 cm apart.The control panel selects the gate-selection sinusoidal drivingvoltage signal with an amplitude of 1 V, a gate width of 3 cycles, and frequencies of 0.3, 1.0, and 1.7 1d f to excite the source transducer.The measurement system collects and The system software based on the LabVIEW platform (Version of 2011) is used for measurement control, which consists of an electric-signal module, a display/storage module, and a slippage/rotation module.
We placed the source and receiver in a tank filled with water.The virtual panel on a computer screen gave instructions to place the two transducer electrodes opposite of one another at 15 cm apart.The control panel selects the gate-selection sinusoidal drivingvoltage signal with an amplitude of 1 V, a gate width of 3 cycles, and frequencies of 0.3, 1.0, and 1.7 f 1d to excite the source transducer.The measurement system collects and stores the signals from the receiver's electric terminals.Finally, the calculated and measured data are compared on the MATLAB platform to verify the correctness of the model established in this paper.
Figure 18 shows the time-domain waves and amplitude spectra for the three selected driving-voltage signals.Figure 19 shows the calculated and measured time-domain waveforms and the corresponding spectra at the receiver's electric terminals.Figure 19a,c Figure 18 shows the time-domain waves and amplitude spectra for the three selected driving-voltage signals.Figure 19 shows the calculated and measured time-domain waveforms and the corresponding spectra at the receiver's electric terminals.Figure 19a,c,e show the normalized time-domain curve, and Figure 19b,d,f show the normalized frequency domain curve.By comparing the red dotted line and solid blue line in Figure 18, the theoretically calculated results in this paper agree well with the experimental measurements.When the driving-voltage signal frequency's center frequency is equal to or close to the transducer's load center frequency, the signal amplitude at the receiver's electric terminals reaches the maximum.When the center frequency of the driving-voltage signal is substantially different from the load center frequency of the source, the signal amplitude at the receiver's electric terminals becomes smaller, and two peaks appear in the corresponding amplitude spectrum curves.One peak (maximum) is near the driving-voltage signal frequency's center frequency, and the other is near the transducer's load center frequency.
Figure 19 also shows some minor differences between the theoretical calculations and experimental measurements, especially for the difference between sub-main peaks in the amplitude spectrumʹs curves.Figure 18 shows the time-domain waves and amplitude spectra for the three selected driving-voltage signals.Figure 19 shows the calculated and measured time-domain waveforms and the corresponding spectra at the receiver's electric terminals.Figure 19a,c,e show the normalized time-domain curve, and Figure 19b,d,f show the normalized frequency domain curve.By comparing the red dotted line and solid blue line in Figure 18, the theoretically calculated results in this paper agree well with the experimental measurements.When the driving-voltage signal frequency's center frequency is equal to or close to the transducer's load center frequency, the signal amplitude at the receiver's electric terminals reaches the maximum.When the center frequency of the driving-voltage signal is substantially different from the load center frequency of the source, the signal amplitude at the receiver's electric terminals becomes smaller, and two peaks appear in the corresponding amplitude spectrum curves.One peak (maximum) is near the driving-voltage signal frequency's center frequency, and the other is near the transducer's load center frequency.
Figure 19 also shows some minor differences between the theoretical calculations and experimental measurements, especially for the difference between sub-main peaks in the amplitude spectrumʹs curves.By comparing the red dotted line and solid blue line in Figure 18, the theoretically calculated results in this paper agree well with the experimental measurements.When the driving-voltage signal frequency's center frequency is equal to or close to the transducer's load center frequency, the signal amplitude at the receiver's electric terminals reaches the maximum.When the center frequency of the driving-voltage signal is substantially different from the load center frequency of the source, the signal amplitude at the receiver's electric terminals becomes smaller, and two peaks appear in the corresponding amplitude spectrum curves.One peak (maximum) is near the driving-voltage signal frequency's center frequency, and the other is near the transducer's load center frequency.
Figure 19 also shows some minor differences between the theoretical calculations and experimental measurements, especially for the difference between sub-main peaks in the amplitude spectrum's curves.
The measured signal at the receiver's electric terminals is from the contributions of the driving-voltage signal, the electric-acoustic conversion of the source, the physical property of the medium (water), and the acoustic-electric conversion of the receiver.At the same time, the calculated signal at the receiver's electric terminals is from the contributions of the driving-voltage signal, the source's electric-acoustic conversion, and the receiver's acoustic-electric conversion.The most prominent peak reflects the source, receiver, or medium's comprehensive frequency response.Meanwhile, the sub-main peak reflects the frequency property of the driving voltage signal used to excite the source.

Discissions and Concluding Remarks
By introducing two components (namely, radiation resistance and radiation mass) into the equivalent electric-mechanical circuit developed by predecessors, this paper established a new frequency-dependent dynamic electric-mechanical network for thinwafer piezoelectric transducers.Because the newly introduced mechanical components are frequency functions, we extended the electric-mechanical circuit model to a parallel circuit network to treat multiple frequency signal transmissions.
For the first time, we derived analytical expressions of electric-acoustic/acousticelectric impulse responses and their transmission functions under the sinusoidal electric/acoustic signals excitation with any given frequency for thin-wafer piezoelectric transducers polarized in the thickness direction.
The driving-voltage signal used to excite the source (in electric-acoustic conversion) is a signal wavelet with many frequency components.The particles inside the source transducer have mass and inertia.The vibration of particles inside the receiver and the electric signal at its electric terminals would still go through a transient process, even for a sinusoidal acoustic signal arriving at the receiver.The measured acoustic signal, i.e., the electric signal at the receiver's electric terminals, will also not be a single-frequency sinusoidal signal.
By carrying out analytical derivations, calculations, analyses, and experiment measurements, we conclude with the following remarks: (a) The effects of the transducer's geometric parameters on various frequencies related to the transducer's properties presented in Figures 5-8 are significant.Load and freeload center frequencies and free-load oscillation frequencies decrease significantly with increased thicknesses for all given radius ranges and a minimum change is observed with respect to the radius after a range of initial increases (Figure 5).The free-load center frequency is slightly higher than the load center frequency.The free-load center and oscillation frequencies are constant for a thin-wafer transducer with fixed physical and geometrical parameters.
In contrast, the loading center frequency varies slightly with the transducer's steadystate vibration frequency (Figures 7 and 8).With respect to the variation tendency of the frequency, the loading oscillation frequency increases with the (electric or acoustic) sinusoidal signal frequency (f) exerted on the transducer.The sinusoidal signal frequency has a more significant effect on the load oscillation frequency in the low-frequency range.There is a minimal impact on the load oscillation frequency in the higher frequency range (Figure 6).(b) Concerning the frequency bandwidth, the source's frequency response bandwidth is greater than the receiver's bandwidth.Due to the electric-acoustic filtering of the source signal, the frequency bandwidth of the radiated acoustic signal is narrower than the driving-voltage signal.Meanwhile, due to the acoustic-electric filtering of the receiver, the frequency bandwidth of the electric signal at the receiver's electric terminals, i.e., the measured acoustic signal, is narrower than that of the acoustic signal arriving at the receiver.(c) Extending the single-frequency circuit network (Figure 3) to parallel multiple-circuit networks (Parts I and III of Figure 12) is necessary for the transitions of the multifrequency signal wavelet.Since thin-wafer transducers' radiation resistance and mass are frequency functions, we may express the electric-acoustic (acoustic-electric) conversion excited by a multifrequency signal wavelet using a parallel transmission network, i.e., an equivalent mechanical-electric circuit network.The cumulative output from the parallel-circuit network represents the acoustic signal wavelet that is radiated outward by the source (and the measured acoustic signal wavelet, i.e., the electric signal wavelet at the receiver's terminals).(d) Comparing the integrated frequency response of the source-receiver with that of the source-medium-receiver, we can use the measured acoustic signal wavelet to obtain the medium's physical properties and internal structure accurately.(e) The center frequency shift of the transducer varies with the property mechanical loading (the medium around the transducer) and that of the driving-voltage signal to excite the transducer.If using free loading as a reference, the transducer exhibits higher sensitivity corresponding to lower frequency radiation resistance and radiation mass.(f) Compared with Fa and Zhao's previous works on thin spherical shells and cylindrical piezoelectric transducers [7,8], the time domain waveform obtained by a thin-wafer transducer vibrates for a longer period, and the corresponding amplitude spectrum has a narrower bandwidth, which is likely from the transducer's different mechanical sizes, shapes, and polarization modes.
The new advancement discussed in this paper on the electric-acoustic conversion of the piezoelectric transducer circuit network should be worth some additional points considering the relevant published work in the literature.
Piquette's research has shown that there is a transient transition process in the electricacoustic conversion of the piezoelectric transducer under the excitation of a sinusoidal electrical signal.The equivalent circuit given by Piquette missed two mechanical components, i.e., radiation resistance and radiation mass, for which their values are frequency functions.Without radiation resistance and radiation mass, the capacitive and inductive reactance generated by the capacitive and inductive components in the traditional equivalent circuit are frequency functions.Still, the values of the capacitive and inductive elements are independent of frequency.Therefore, that work missed carrying out further analyses on the transient transition process of the electric-acoustic conversion of multifrequency electrical signal (wavelet) excitation transducers [52,53].Our paper introduced the radiation resistance and radiation mass into the equivalent circuit.These components are frequency functions, and the related and generated radiation force resistance and mechanical reactance are also frequency functions.
Although Fa et al. introduced radiation resistance and radiation mass into the equivalent circuit of radially polarized thin spherical shell piezoelectric transducers [61,65] when multifrequency electrical signals stimulated the transducers, however, that work adjusted the transducer's electric-acoustic conversion impulse response artificially and the values of radiation resistance and radiation mass constantly.Therefore, the obtained results from calculations and analyses have human subjective factors [65].In the current paper, based on linear superpositions and the Fourier transform, the amplitude and phase spectra of the multifrequency electric-driving signal are discretized into many frequency components with different amplitude and initial phases.The equivalent circuits of a thin-wafer transducer polarized in the thickness direction form a parallel electric-acoustic conversion network corresponding to each electric driving-signal frequency (with different radiation resistances and radiation masses).Yet, the transmission network model proposed in this paper is the continuation and enhancement of the model of acoustic source transducers excited by multifrequency electrical signals reported by Fa et al. [65].
The method described in this paper differs from those of Mason and KLM for thin-disc piezoelectric transducers with respect to thickness polarization [49], and the transitional transient process generated by the electric-acoustic/acoustic-electric conversion of the transducer proposed in this paper is conceptually different from the transient characteristics of the transducer reported by Redwood [56].In this paper, the spectrum corresponding to the abovementioned process is the intrinsic noise generated by the particles inside the transducer during the electric-acoustic/acoustic-electric conversion.We noted that the equivalent circuits of Mason and KLM used conceptual physical knowledge to describe the characteristics of thickness polarization piezoelectric transducers.In contrast, our new piezoelectric transducer model combines physics and signal transmission/processing knowledge to describe the characteristics of thin-wafer piezoelectric transducers.
Meanwhile, in the equivalent circuit models of Mason and KLM, some physical parameters of piezoelectric materials are assumed to be numerically complex, as shown in Tables 2 and 3 of reference [49], which is against physical nature.Even the propagation speed of the acoustic wave included an imaginary part, e.g., v = 2700(1 + 0.005i).In our new model, we applied all factual physical parameters (real numbers) of piezoelectric materials for calculations and analyses, which is physically more meaningful.
Considering the work of Sherrit et al. [49] on Mason and KLM Equivalent Circuits, the values of specific electrical components and mechanical components are independent of frequency, even though the impedance of these components can be frequency dependent.Similarly to ordinary equivalent circuits, the impedance values corresponding to the capacitive/inductive elements in the equivalent circuit are indeed frequency-dependent.In our new model, the two newly introduced mechanical components, i.e., radiation resistance and radiation mass, are numerically frequency functions [61,65]; the mechanical values of other electrical components are fixed and unchanged, which leads to obtaining a circuit of a different nature than the equivalent circuits reported in the previous literature; this is the basis for the establishment of a parallel electric-acoustic conversion transmission network, and it reflects the realistic situation of transducers better.
Additionally, the equivalent circuit of Mason and KLM is a six-terminal mechanicelectric network, i.e., four mechanical and two electrical terminals corresponding to two acoustic ports and one electrical port.Two acoustic ports correspond to the front and back surfaces of the thickness polarization transducer, and one electrical port corresponds to the two electrodes on the transducer's front and back disk surfaces [49].The disk's back surface was tightly and entirely cemented with a thick rigid solid medium (assumed approximated that the back surface does not vibrate).Its front surface comprised a radiating surface coated with a thin insulating layer, considering that the thickness of the insulating layer was thin relative to the wavelength of the radiated acoustic wave, and its influence on the acoustic wave emitted by the transducer is ignorable.In our new physical model, we use the compression/tension deformation of a uniform rod (sheet) to describe the rear surface fixation and the front surface radiation of the thickness polarization thin-wafer transducer.By carrying this out, the above six-terminal machine-electrical equivalent circuit of Mason and KLM was simplified into a four-terminal equivalent circuit similar to but different from Camp's thin ring piezoelectric transducer, noting that Camp's equivalent circuit does not contain two frequency-related mechanical components, namely radiation resistance and radiation mass [63].
The transient transition process of particle vibrations inside the transducer reported in this paper also conceptually differs from the transient characteristics reported by Redwood [56].Redwood's qualitative description of the transient properties of piezoelectric transducers is as follows: "In the generation of ultrasonic waves a transient electrical signal is applied to the transducer and this produces a transient mechanical vibration, while in detection, the application of a transient mechanical signal produces an electrical vibration" [56].If so, the transient electrical signal (applied to the transducer) and the transient mechanical vibration are essentially the electric driving signals of the acoustic source transducer and the acoustic signal radiated from the transducer.Then, both transient processes are signal wavelets within waveforms in the time domain with a specific time duration.In our model, when a transient mechanical signal is applied to the transducer, it generates a dynamic vibration.The concept of intrinsic noise generated from vibration particles inside viscous solids is introduced into the electric-acoustic and acoustic-electric conversion processes of thick-polarized thin-wafer piezoelectric transducers, which provides a theoretical basis for both extracting the intrinsic noise generated by the measured medium from the acoustic signal obtained in acoustic measurements and using the intrinsic noise to invert the physical characteristics and internal microstructure of the measured medium.It is essentially a vibrational wave propagating to the observation point to make its electrical terminals output a time-domain waveform with a specific duration with respect to the electric signal wavelet.
We conclude that the updated parallel transmission network reported in this paper with the two new physical factors is more realistic relative to the physical nature of a multifrequency signal wavelet excitation transducer.It is applicable not only to conventionalsized but also to small-sized piezoelectric transducers, and it is universal in the descriptions of transducers.

Figure 1 .
Figure 1.A schematic presentation of a thin-wafer transducer polarized in its thickness direction.

Figure 1 .
Figure 1.A schematic presentation of a thin-wafer transducer polarized in its thickness direction.

Figure 2 .
Figure 2. Compression/tension deformation of a uniform rod.

Figure 2 .
Figure 2. Compression/tension deformation of a uniform rod.

Figure 3 .
Figure 3.The time-domain mechanical-electric equivalent circuits of a transducer: (a) the transducer as a source; (b) the transducer as a receiver.

Figure 4 .
Figure 4.The s-domain mechanical-electric equivalent circuits of the transducer: (a) source; (b) receiver.

Figure 3 .
Figure 3.The time-domain mechanical-electric equivalent circuits of a transducer: (a) the transducer as a source; (b) the transducer as a receiver.

Figure 3 .
Figure 3.The time-domain mechanical-electric equivalent circuits of a transducer: (a) the transducer as a source; (b) the transducer as a receiver.

Figure 4 .
Figure 4.The s-domain mechanical-electric equivalent circuits of the transducer: (a) source; (b) receiver.

Figure 4 .
Figure 4.The s-domain mechanical-electric equivalent circuits of the transducer: (a) source; (b) receiver.

Figure 5 .
Figure 5. Relationships of the free-load center, load center, and free-load oscillation frequencies versus the geometric dimension of the transducer: (a) frequencies and thickness; (b) frequencies and radius.

Figure 6
Figure 6 is the calculated relationships of load oscillation frequencies ( 1sd f and 3sd f ) of both source and receiver versus the sinusoidal frequency (f) of the excitation signal exerted on the transducer (electric or acoustic).It displays 1sd

Figure 6 .
Figure 6.Relationships of the load oscillation frequency versus the frequency of sinusoidal (electric or acoustic) signal exerted on the transducer.

Figure 5 .
Figure 5. Relationships of the free-load center, load center, and free-load oscillation frequencies versus the geometric dimension of the transducer: (a) frequencies and thickness; (b) frequencies and radius.

Figure 5 .
Figure 5. Relationships of the free-load center, load center, and free-load oscillation frequencies versus the geometric dimension of the transducer: (a) frequencies and thickness; (b) frequencies and radius.

Figure 6
Figure 6 is the calculated relationships of load oscillation frequencies ( 1sd f and 3sd f ) of both source and receiver versus the sinusoidal frequency (f) of the excitation signal exerted on the transducer (electric or acoustic).It displays 1sd

Figure 6 .
Figure 6.Relationships of the load oscillation frequency versus the frequency of sinusoidal (electric or acoustic) signal exerted on the transducer.

Figure 6 .
Figure 6.Relationships of the load oscillation frequency versus the frequency of sinusoidal (electric or acoustic) signal exerted on the transducer.

Figure 7 .
Figure 7.The electric-acoustic impulse responses of the source and the corresponding amplitude spectra: (a) electric-acoustic impulse responses; (b) amplitude spectra.We normalized the electricacoustic impulse responses and amplitude spectra using the maximum values of the impulse response and amplitude spectrum for the case of free-loading, respectively.The cyan curve is the electric-acoustic impulse response and the corresponding amplitude spectrum for free-loading.The other colors are those with mechanical loads: the lines in magenta, blue, red, and black are the transducer's electric-acoustic impulse responses and the corresponding amplitude spectra at selected frequencies of 0.1 1d f , 0.5 1d f , 1.0 1d f , and 10 1d f , respectively.Here, 1d f =1.000575 MHz.

Figure 7 . 28 Figure 8 .
Figure 7.The electric-acoustic impulse responses of the source and the corresponding amplitude spectra: (a) electric-acoustic impulse responses; (b) amplitude spectra.We normalized the electricacoustic impulse responses and amplitude spectra using the maximum values of the impulse response and amplitude spectrum for the case of free-loading, respectively.The cyan curve is the electricacoustic impulse response and the corresponding amplitude spectrum for free-loading.The other colors are those with mechanical loads: the lines in magenta, blue, red, and black are the transducer's electric-acoustic impulse responses and the corresponding amplitude spectra at selected frequencies of 0.1 f 1d , 0.5 f 1d , 1.0 f 1d , and 10 f 1d , respectively.Here, f 1d = 1.000575MHz.Micromachines 2023, 14, 1641 14 of 28

Figure 8 .
Figure 8.The acoustic-electric impulse responses of the receiver and the corresponding amplitude spectra: (a) acoustic-electric impulse responses; (b) amplitude spectra.We have normalized acousticelectric impulse responses and amplitude spectra using the maximum value of the impulse response and amplitude spectrum for the case of free load, respectively.The cyan curve is the acoustic-electric impulse response and the corresponding amplitude spectrum for the free load.The other colors are those with mechanical load: the lines in magenta, blue, red, and black are the transducer acousticelectric impulse responses, and the amplitude spectra of the corresponding system functions at selected frequencies of 0.1 f 3d , 0.5 f 3d , 1.0 f 3d , and 10 f 3d .Here, f 3d = 1.005189MHz.

Figure 8 .
Figure 8.The acoustic-electric impulse responses of the receiver and the corresponding amplitude spectra: (a) acoustic-electric impulse responses; (b) amplitude spectra.We have normalized acoustic-electric impulse responses and amplitude spectra using the maximum value of the impulse response and amplitude spectrum for the case of free load, respectively.The cyan curve is the acoustic-electric impulse response and the corresponding amplitude spectrum for the free load.The other colors are those with mechanical load: the lines in magenta, blue, red, and black are the transducer acoustic-electric impulse responses, and the amplitude spectra of the corresponding system functions at selected frequencies of 0.1 3d f , 0.5 3d f , 1.0 3d f , and 10 3d f .Here, 3d f = 1.005189MHz.

Figure 9 .
Figure 9. Surface vibration velocity and corresponding amplitude spectrum of the source excited by sinusoidal driving-voltage signals with different frequencies: (a,c,e,g) are the waveforms in the time domain; (b,d,f,h) are the corresponding amplitude spectra.

Figure 9 .
Figure 9. Surface vibration velocity and corresponding amplitude spectrum of the source excited by sinusoidal driving-voltage signals with different frequencies: (a,c,e,g) are the waveforms in the time domain; (b,d,f,h) are the corresponding amplitude spectra.

Figure 10 .
Figure 10.Electric signals converted from sinusoidal acoustic signals with different frequencies: (a,c,e,g) are the waveforms in the time domain; (b,d,f,h) are the corresponding amplitude spectra.

Figure 10 .
Figure 10.Electric signals converted from sinusoidal acoustic signals with different frequencies: (a,c,e,g) are the waveforms in the time domain; (b,d,f,h) are the corresponding amplitude spectra.

Figure 11
shows the voltage signals at the receiver's electric terminals.Micromachines 2023, 14, 1641 16 of 28 response of the source-receiver.(iii) If the frequency of the sinusoidal driving-voltage signal is far from the center frequency of the integrated frequency response of the sourcereceiver, the number of periods of the electric signal at the receiver's electric terminals from the zero state to the stable sinusoidal vibration state increases, and the amplitude of the time domain waveform and the magnitude of amplitude spectrum decreases.

Figure 11 .
Figure 11.Waveform and amplitude spectrum of electric signals at the receiver's electric terminals for the integrated response of the source-receiver: (a,c,e,g) are the waveforms in the time domain; (b,d,f,h) are the corresponding amplitude spectrum in the frequency domain.

Figure 11 .
Figure 11.Waveform and amplitude spectrum of electric signals at the receiver's electric terminals for the integrated response of the source-receiver: (a,c,e,g) are the waveforms in the time domain; (b,d,f,h) are the corresponding amplitude spectrum in the frequency domain.

Figure 12 . 2 v
Figure 12.A schematic representation of a parallel transmission network of an acoustic measurement process.In this network, the related physical parameters include the j th frequency component ( 1j U ) in the driving-voltage signal ( 1 U ) (j = 1, 2, …… N); the electric-acoustic impulse  ) of the medium, the propagation time ( t  ) of the acoustic wave from the source to the receiver; the acoustic signal ( 3 v ) reaching the receiving transducer; the k th frequency component ( 3 v k ) (k = 1, 2, …… M) of a sound wave signal; the acoustic-electricity impulse response (  ) of the receiving transducer corresponding to the k th sinusoidal fre- quency component in the acoustic signal ( 3 v ); output electric signal ( 3 k U ) of the receiving

Figure 12 .
Figure 12.A schematic representation of a parallel transmission network of an acoustic measurement process.

Figure 13 .
Figure 13.The gate-selection sinusoidal driving-voltage signal with three periods: (a) the waveform; (b) the amplitude spectrum; (c) the phase spectrum.

Figure 13 .
Figure 13.The gate-selection sinusoidal driving-voltage signal with three periods: (a) the waveform; (b) the amplitude spectrum; (c) the phase spectrum.

Figure 14 .
Figure 14.Figures (a-d) display the convolutions for several frequency components (0.4 1d f , 0.8 1d f , 1.2 1d f , and 1.6 1d f ) and the corresponding electric-acoustic impulse response.Plot (e) shows the normalized accumulative output from the circuit network.And (f) is the normalized amplitude spectrum of the accumulative signal output.

Figure 14 .
Figure 14.Figures (a-d) display the convolutions for several frequency components (0.4 f 1d , 0.8 f 1d , 1.2 f 1d , and 1.6 f 1d ) and the corresponding electric-acoustic impulse response.Plot (e) shows the normalized accumulative output from the circuit network.And (f) is the normalized amplitude spectrum of the accumulative signal output.

Figure 15 .
Figure 15.Plots (a-d) are the electric signals when the source is excited by several frequency components (0.4 1d f , 0.8 1d f , 1.2 1d f , and 1.6 1d f ) of the gate sinusoidal driving-voltage signal.Plot (e) is the normalized time-domain waveform of the accumulative output.Plot (f) is the normalized amplitude spectrum of the accumulative signal output.

Figure 15 .
Figure 15.Plots (a-d) are the electric signals when the source is excited by several frequency components (0.4 f 1d , 0.8 f 1d , 1.2 f 1d , and 1.6 f 1d ) of the gate sinusoidal driving-voltage signal.Plot (e) is the normalized time-domain waveform of the accumulative output.Plot (f) is the normalized amplitude spectrum of the accumulative signal output.

Figure 16 .
Figure 16.The physical image of the transducer.

Figure 17
Figure 17 is the structural chart of hardware and software modules used in experimental testing.The hardware is composed of three modules: the mechanical module, the electronic hardware module, and the systemʹs software module.The mechanical module includes a stepper motor, a moving slide rail, a tank, and a steering gear.The electronic hardware module consists of a signal generator (AFG-125P, made by Gwinstek in China) with up to a bandwidth of 25 MHz and a resolution of 1 µHz, a microcontroller (STC12C5A60S2, made by STC MCU Limited in China) used to control the spatial position of the source and receiver, and a data collector (PXI-5922 made by NI in USA) with a sampling rate of 5-15 MHz and 16-24 bits.The system software based on the LabVIEW platform (Version of 2011) is used for measurement control, which consists of an electric-signal module, a display/storage module, and a slippage/rotation module.

Figure 17 .
Figure 17.The structural chart of the experimental testing apparatus.

Figure 16 .
Figure 16.The physical image of the transducer.

Figure 17
Figure 17 is the structural chart of hardware and software modules used in experimental testing.The hardware is composed of three modules: the mechanical module, the electronic hardware module, and the system's software module.The mechanical module includes a stepper motor, a moving slide rail, a tank, and a steering gear.The electronic hardware module consists of a signal generator (AFG-125P, made by Gwinstek in China) with up to a bandwidth of 25 MHz and a resolution of 1 µHz, a microcontroller (STC12C5A60S2, made by STC MCU Limited in China) used to control the spatial position of the source and receiver, and a data collector (PXI-5922 made by NI in USA) with a sampling rate of 5-15 MHz and 16-24 bits.

Figure 16 .
Figure 16.The physical image of the transducer.

Figure 17
Figure 17 is the structural chart of hardware and software modules used in experimental testing.The hardware is composed of three modules: the mechanical module, the electronic hardware module, and the systemʹs software module.The mechanical module includes a stepper motor, a moving slide rail, a tank, and a steering gear.The electronic hardware module consists of a signal generator (AFG-125P, made by Gwinstek in China) with up to a bandwidth of 25 MHz and a resolution of 1 µHz, a microcontroller (STC12C5A60S2, made by STC MCU Limited in China) used to control the spatial position of the source and receiver, and a data collector (PXI-5922 made by NI in USA) with a sampling rate of 5-15 MHz and 16-24 bits.The system software based on the LabVIEW platform (Version of 2011) is used for measurement control, which consists of an electric-signal module, a display/storage module, and a slippage/rotation module.

Figure 17 .
Figure 17.The structural chart of the experimental testing apparatus.

Figure 17 .
Figure 17.The structural chart of the experimental testing apparatus.
,e show the normalized time-domain curve, and Figure 19b,d,f show the normalized frequency domain curve.Micromachines 2023, 14, 1641 22 of 28 stores the signals from the receiver's electric terminals.Finally, the calculated and measured data are compared on the MATLAB platform to verify the correctness of the model established in this paper.

Figure 18 .
Figure 18.Time domain waveforms and amplitude spectra of three gate-selected sinusoidal voltage signals with different frequencies: (a,c,e) are the normalized time domain waveform; (b,d,f) are the corresponding amplitude spectra.

Figure 19 .
Figure 19.The electric signals at the receiver's electric terminals under the excitations of the above three different gated sinusoidal electric signals: (a,c,e) waveforms; (b,d,f) amplitude spectra.Blue solid lines come from the theoretical calculation, and red dotted lines come from the experimental measurement.

Figure 18 .
Figure 18.Time domain waveforms and amplitude spectra of three gate-selected sinusoidal voltage signals with different frequencies: (a,c,e) are the normalized time domain waveform; (b,d,f) are the corresponding amplitude spectra.

Figure 18 .
Figure 18.Time domain waveforms and amplitude spectra of three gate-selected sinusoidal voltage signals with different frequencies: (a,c,e) are the normalized time domain waveform; (b,d,f) are the corresponding amplitude spectra.

Figure 19 .
Figure 19.The electric signals at the receiver's electric terminals under the excitations of the above three different gated sinusoidal electric signals: (a,c,e) waveforms; (b,d,f) amplitude spectra.Blue solid lines come from the theoretical calculation, and red dotted lines come from the experimental measurement.

Figure 19 .
Figure 19.The electric signals at the receiver's electric terminals under the excitations of the above three different gated sinusoidal electric signals: (a,c,e) waveforms; (b,d,f) amplitude spectra.Blue solid lines come from the theoretical calculation, and red dotted lines come from the experimental measurement.