Dynamics Analysis on Piezoelectric Laminated Vibrator and Optimization of PZT Position

Piezoelectric laminated structure is widely used as actuator’s drive part. The different position of PZT on a piezoelectric vibrator causes different incentive effects. Therefore, seeking an optimal PZT position is of great significance to improve actuator’s drive forces and electromechanical conversion efficiency. In this research, the optimization of PZT position was studied using the approximate solution of piezoelectric vibrator mode shape with mutation sections. The vibration mode function was expressed as a linear superposition of the admissible function according to Rayleigh Ritz method. Then solving of functional variation was converted into the solving of the coefficient matrix of the admissible function by Hamilton’s principle.Through analyzing the forms of admissible functions, the admissible functions that satisfied the boundary conditions of displacement were chosen. For a given vibrator, approximate function for natural frequency and vibration mode was calculated in given admissible functions. Calculated values and experimental results were compared. Results showed that themore items an admissible function sequence had the closer the calculated results were to the experimental values.The errors of calculations were analyzed based on the selection of admissible functions and electromechanical coupling effect. Optimization of PZT position was achieved by analyzing the mode forces of the piezoelectric laminated vibrator.


Introduction
Piezoelectric laminated structure is widely used as actuator's drive part as it is well characterized with simple structure, convenient control, quick response, and no electromagnetic pollution [1].Researchers from all over the world have been involved in the study on theoretical models of such structures extensively [2][3][4][5].Research on how to paste PZT on traveling wave ultrasonic motor indicated that the electromechanical conversion efficiency of piezoelectric ceramics was closely related to the force coefficient which depended on the size and position of the ceramic transducer as well as the mode shape of the stator [6].Thus, to improve the electromechanical conversion efficiency of the vibrator theoretically, the first thing that must be done is to obtain the vibrator's mode shape and then to determine optimal size and location of a piezoelectric ceramics transducer accordingly.Koplow et al. conducted dynamic analysis on the free vibration of discontinuous beam with one-step change in cross section and predicted the dynamic response of the beam [7].Sadri et al. established a theoretical model for a thin plate covered with piezoelectric sheets and derived control equations for the plate based on Rayleigh-Ritz method; based on genetic algorithms the optimal placement of piezoelectric actuators was determined, but follow-up experiment verification has not been reported [8].Bashash et al. divided a variable cross section stepped homogeneous material beam with arbitrary boundary into several segments and built models for them.As a result, characteristic matrix was formed by taking boundary conditions and continuity conditions into consideration [9].For nonhomogeneous piezoelectric beam, however, this approach created a high-order algebraic equation, which is difficult to solve.Crawley and de Luis studied the strain transmitting mechanism from piezoelectric actuator to substructure and conducted researches on PZT optimization position of constant cross section piezoelectric beam and ignored the influence of PZT on its vibration mode [10].From the perspective of vibration suppression, Ren and Jiang presented a systematic approach for the free vibration analysis and forced response of the beam bonded with PZT employing the travelling wave method [11].They found that the PZT bonded position near the fixed end in beam had the powerful actuated capability.However, experiments demonstrated that this conclusion was only applicable to vibration suppression.Cho et al. established an equivalent electric circuit for piezoelectric bimorph beam by using piezoelectric PZT as equivalent to an impedance element [12], and the researches provided the new idea for mathematical model of piezoelectric laminated structure.However, the PZT position optimization was not mentioned.Xu and Zhou presented a two-dimensional analysis for piezoelectric beam with variable thickness under the boundary condition of the simple support [13].Mao proposed a solution for the problem of finding the shape of piezoelectric mode sensors for nonuniform Euler-Bernoulli beams with rectangular cross sections [14].Cupiał discussed a perturbation solution of the natural frequencies and mode shapes of a piezoelectric rectangular plate [15].
Daraji and Hale applied piezoelectric elements in active vibration reduction and proposed fitness and objective functions to seek optimal distribution of segmented sensors on a flexible plate [16,17].These researchers also developed a model for complex isotropic plate using FEM and Hamilton's principle [18].
Cazzulani et al. used the independent mode space control logics to achieve the optimal actuator and sensor placement [19].This most popular model control logic used independently controlled modes' damping and stiffness.In this method, Ambrosio et al. proposed a negative derivative feedback strategy, which feeds velocity back to the mode and also generates a control input proportional to the compensator velocity.The feedback was able to filter out both the higher and lower uncontrolled components.In active vibration control of smart structures, the actuator and sensor placement was a key point of the control system design.Ambrosio et al. [20] also proposed an H2 norm approach to solve actuator and sensor placement.This approach was examined using a finite element model of a square plate.Piezoelectric patch actuators and acceleration sensors were mounted on the three sides of the square plate.To seek for optimal actuator and sensor placement, researchers have realized that more research effort should be made to driving forces and deformation of a piezoelectric vibrator in piezoelectric motor, micro pump, and scanner.
PZT's excitation to the vibrator is achieved through the substrate deformation caused by the piezoelectric effect.Thus, the matching relations between PZT's position and mode shape of the vibrator have great influences upon driving forces and electromechanical conversion efficiency.From the point of mode shape, this paper describes the position optimization of a PZT on a piezoelectric laminate vibrator.First, the vibrator is divided into several segments according to the cross section mutation, and then the dynamic analysis for the piezoelectric laminate structure was carried out.Then, approximate solutions of vibration modes were obtained using Rayleigh-Ritz method.These approximate solutions were compared to experimental results at their resonant frequencies.Furthermore, force coefficients at different PZT positions were discussed.These results will provide a theoretical foundation for optimizing the design of the piezoelectric laminated beam.In this paper the research objectives were to (1) conduct dynamic analysis of the piezoelectric laminated structure, (2) derive approximate solutions of vibration modes using Rayleigh-Ritz method, and (3) propose an optimization method for the design of a piezoelectric laminated beam.

Dynamic Analysis on the Shape of Vibration Modes
Assume that coordinates of cross section change positions of the stepped beam of piezoelectric laminate structure in  direction are  1 ,  2 , and  3 , as shown in Figure 1.The piezoelectric laminate structure was divided into several segments with equal cross section and differential equations were given for them, respectively.A characteristic matrix for the whole structure was assembled based on the continuity relationship between these segments [9].Differential equation of the segment  is where   (), ()  , and   are the mode shape, bending rigidity, and mass of a unit long beam, respectively.Define Then (1) can be expressed as The general solution for (3) is where   ,   ,   , and   are undetermined coefficients.For a cantilever beam, deflection and rotation angle on the left end are zero.That is to say, In this paper, this variable cross section beam is divided into three segments.At the right end of the beam, the boundary condition is free, and both bending moment and shear force are zero.Consider At the mutation sections of the beam, the deflection, rotation angle, bending moment, and shear force are continuous.That is, The characteristic equation of the system is obtained by substituting (5)∼( 10) into (4).Obviously,   ( = 1, 2, 3) is a function of the whole beam's natural frequency .  can be expressed as where The characteristic equation can also be expressed in the form of a matrix: where J is a characteristic matrix and P is a feature vector which means a coefficient vector of mode shape.Consider ) J consists of three parts, as shown in (14).The first part J1 and the third part J3 represent the boundary conditions at the positions  = 0 and  =  3 , respectively.Part J2 represents the continuous relationship between these segments at the mutation sections.Consider where J1 and J3 are 2 × 12 order matrices and J2 is an 8 × 12 order matrix.The determinant of J is zero since the feature vector P has a nonzero solution.Consider     J 12×12     = 0.
Analytical expression for  is obtained by solving (15), and then the coefficients of mode function in (4) can be solved.It is difficult to get a nonzero solution since the equation is a higher-order transcendental equation after being unfolded.In practice, the equation is often approximately solved by making some reasonable simplifications.

Approximate Solution of the Mode Shape
The essence for solving mode shape is to obtain the real vibrational track of the beam under a certain order mode shape.According to Hamilton's principle, the real motion's variation for Hamilton's processing object equals zero: where  is a variational symbol,  is Hamilton's magnitude, and   is a Lagrange function defined as   =   − , where   is system kinetic energy,  is potential energy, and  is generalized coordinates.Hamilton's principle gives a way to find a true movement from all possible movements and it is generally applicable in mechanics.It has become a basic principle in mechanics field.Rayleigh-Ritz method is a direct solution for this kind of variation problems, which extends Hamilton's principle into elastic dynamics.Its basic concept is to assume the approximate solution of a vibration mode function as follows: where   () is defined as admissible function,   is undermined coefficient, and  is the item number of the admissible function sequence.Equation (17) shows that approximate solution of a vibration mode function can be expressed by a linear combination of function sequences  1 , Supposing that the beam's deflection is (, ), then the kinetic energy of the system is Potential energy of the system can be written as The fonctionelle of Hamilton's processing object for the beam's bending vibration can be expressed as where ()  is the mass of a beam in unit length and ()  is flexural rigidity of the beam.
When the piezoelectric beam is doing a certain-order mode vibration,  (, ) =  () sin ( + ) . ( Substitute formula (22) into formula (21): Equation ( 23) is the fonctionelle of the function depending on the independent variable  and .In one cycle, which is in a range of [0, 2/], (23) can be simplified as a fonctionelle of () if, integrating on , Substitute ( 17) into (24); then the fonctionelle () can be changed into a function of ( 1 ,  2 , . . .,   ).Consider Equation (25) can also be written in a quadratic format about  1 ,  2 , . . .,   .Consider ], it is called generalized stiffness.Equation (26) can be written in the form of a quadratic matrix: where Calculate the derivatives on  1 ,  2 , . . .,   in formula (27), respectively, and let them be equal to zero.It can be expressed by a matrix as follows: In order to get the nonzero solutions of , its determinant of feature matrix should be equal to zero.Consider Unfold the elements of matrix M and K into a segmented form: where ()  and ()  are mass of the unit length beam and flexural rigidity of the segment  of the uniform beam ( = 1, 2, 3).And here Then generalized mass matrix and generalized stiffness matrix can be solved using (31).
The choice for an admissible function   (), which is shown in (17), has significant influence on calculating accuracy and convergence rate.Ideal admissible function should not only meet the displacement boundary conditions of structural vibration but also meet the force boundary conditions.However, it is very difficult to obtain an admissible function which meets both boundary conditions.In this research, an admissible function can be acceptable as long as it meets the displacement boundary conditions.The larger the number of items  is the more accurate the approximate solution becomes.When  → ∞, the approximate solution is its exact solution.
An admissible function sequence that meets the displacement boundary conditions can be assumed as When  = 1,  = 2, and  = 3, approximate vibration mode functions are as follows:

Comparison of Theoretical and Experimental Results
In this section, the approximation method for the mode used in earlier sections was verified through vibration tests.As The structure of the piezoelectric laminated vibrator in the experiment is shown in Figure 1.The PZT material was P81, and substrate material was phosphor bronze.The properties of the materials and the structural parameters are shown in Tables 1 and 2. The width of the vibrator was 16 mm.
These dimensions indicate that the midpoint of the ceramic sheet was exactly coincident with the midpoint of the substrate.When the vibrator works in higher orders of vibration modes (over the third order), the amplitude was very small.Thus, only the first two modes were considered in the paper.By taking data of Tables 1 and 2 into formulas (31)-(32), elements of matrix M and K were obtained as shown in (30).Solving (30) value  was then obtained while  brings the information of resonant frequency ( =  2 ).Then, taking the value  into (29), coefficients of mode function () can be solved as  1 ,  2 , and  3 .The calculation results showed that only resonance frequency of the first-order bending vibration is obtained as  1 = 310 Hz when () took a combination of admissible function with the first term in (34).When () took a combination of admissible function with the first two terms in (35), resonance frequencies of the first-order and second-order bending vibrations were obtained as  1 = 230 Hz and  2 = 234 Hz.When () was a combination of admissible function with the first three items in (36), resonance frequencies of the first-order and second-order bending vibrations were  1 = 216 Hz and  2 = 1450 Hz as shown in Table 3.
To verify the results from the method described above, a laser based noncontact instrument, Doppler vibrometer, was used to perform the scan test on the vibrator.This laser Doppler vibrometer system's components and working principle are shown in Figure 2. Operating parameters (e.g., scan range and excitation voltage) are specified on the terminal device.In this experiment, scanning frequency ranges from 10 to 3000 Hz.Excitation voltage was 60 V, and the scanning signal was assigned by the master controller, which motivates the vibrator to work.Scanning laser head (OFV-56) sends out frequency-stabilized laser to scan the surface of the vibrator and collects the reflected light.Doppler's frequency shift signal, which is proportional to the target's velocity, is generated through interference phenomenon.After the signal is processed by a decoded controller, velocity and displacement analog values for the measured object are received.In addition, the resonant frequencies and vibration modes of the vibrator can also be obtained.All of these data will go through a high-speed A/D converter to be processed and displayed on a PC.The calculated values and experimental results are shown in Table 3.
It can be seen from Table 3 that with just taking the first term of the admissible function into consideration ( = 1), the error of the theoretical value relative to the experimental value is 46.7% and that with just taking the first two terms of the admissible function sequence to form a basic function ( = 2), the error of the theoretical value relative to the experimental value is 28.3% for the first-order bending vibration.Similarly, it is 39.5% for the second-order bending vibration.With taking the first three terms to form a function ( = 3), the relative error between them is 23.6% for the first-order bending vibration and 2.1% for the secondorder bending vibration.These results indicate that with more terms the approximate solution becomes more accurate.The calculated values are always greater than experimental values.These errors are mainly due to following factors.
(1) Approximation of the Admissible Function.The admissible functions only meet the displacement boundary conditions, not the force boundary conditions.The beam's ends are free in its physical model, but there are bending moment and shear force on the right end of the beam when deriving the admissible function, which is equivalent to the imposed additional constraints.
(2) Influence Caused by the External Electric Field.With the action of an electric field, the stress-strain relationship of the PZT beam can be approximately expressed using following equation: where    is the stress of PZT in  direction;    is the elastic modulus of PZT in  direction;    is strain of PZT in  direction;  31 is the piezoelectric constant of PZT;  is the voltage of external electric field; and ℎ  is the thickness of PZT.
Equation (37) can be simplified as follows: Equations ( 37) and (38) indicate that the stress of the PZT generated from a certain strain not only is associated with the strain itself but also is related to the piezoelectric effect induced from external electric field.Under the same strain, the stress caused by PZT with external electric field became smaller than that with no external electric field.That is to say that the effective elastic modulus  effective is smaller than that when the external electric field is open.The higher the applied voltage in external electric field is the faster the stress decreases.
The theoretical calculation in this paper was carried out without consideration of electromechanical coupling effect from the PZT.While measuring these values at the driving voltage of 60 V, the effective elastic modulus reduces due to the action of external electric field.Thus, measured mode frequencies were less than the theoretical values.
When  = 3, for example, the first-order and secondorder approximate vibration mode functions can be written as  (1)   (2) For the piezoelectric bimorph vibrator described in Figure 1, using the second-order mode can get higher electromechanical coupling effect [21].Thus, optimized analysis of later part is aimed at the second-order mode in this study.After the length  3 of piezoelectric beam is normalized, the curve of the second-order vibration mode shape drawn using (39) is shown in Figure 3.When comparing this curve with the vibration curve of a uniform beam, it indicates that (1) when the midpoint of the PZT beam is exactly coincident with that of the substrate the second-order vibration curve's extreme point of multiple-step beam is almost coincident with that of the uniform beam and that (2) after the vibration amplitudes on the two ends are normalized the deformation of multiplestep beam is smaller than that of uniform beam due to the influence caused by PZT.

Optimization of PZT Position
The second-order vibration mode of the vibrator discussed above was the case when the length  piezo of PZT was nearly equal to the length of the substrate and the position of PZT was located at the center of the substrate.Next scenario to be discussed is when the length of  piezo is smaller than the half wavelength of the vibration mode curve.In addition, optimal pasting position for a PZT with ideal paste condition will be discussed.
Under the ideal paste condition of the PZT, the thickness of adhesive layer between the PZT and the substrate tends to zero, and then the force exerted on the substrate by the PZT focuses on two endpoints of the ceramic sheet [10]; the relationship between the force  and the parameters of PZT is expressed as At the left end of the ceramic sheet, the force  + is expressed as At the right end of the ceramic sheet, the force  − is expressed as In ( 40), (41), and (42),  + and  − are concentrated force exerted on the substrate by PZT;   is the elastic modulus of substrate;  is the width of PZT and ℎ  is the thickness of the substrate;  is the stiffness ratio between the substrate and PZT;    is the surface strain of the substrate;  is a constant that reflects the strain distribution of the beam, and, for a linear Euler-Bernoulli beam,  = 6; Λ =  31 /ℎ  , where  31 is the piezoelectric constant,  is the voltage of external electric field, and ℎ  is the thickness of PZT.
The action force and moment of PZT on the piezoelectric vibrator are shown in Figure 4, in which  is the distance from fixed end to the midpoint of PZT,  0 is the moment exerted from PZT to the vibrator.Generalized force  of the piezoelectric laminated vibrator is expressed as follows: where  is vibration mode of the vibrator and ℎ  is the thickness of the substrate.
When the system is vibrating, the strain on top and bottom surfaces of the substrate is Substituting (44) into (41) and (42), respectively, and then substituting (41) and ( 42) into (43), respectively, thus generalized vibrating force  can be achieved as follows: For the first half of (45), except for , it can be expressed as Then (45) can be simplified as where When the piezoelectric vibrator works near the resonant frequency, it can be treated as a single degree of freedom vibration, and the vibration equation is Substituting (45) into (49), then (49) became Equations ( 50) and (48) indicate that the mode force not only is associated with piezoelectric material properties  31 , voltage  of external electric field, and structural parameter ℎ  but also related to the position of PZT relative to the mode shapes of the vibrator.In this study, the part in brackets [] of (48) is defined as the positional force coefficient Π of PZT.Consider  According to Figures 1 and 4, Then (51) can be expressed as From ( 51) and (53) it can be seen that, with the material and geometry sizes fixed, the vibration shape curve of the piezoelectric vibrator and the positional force coefficient Π of PZT will change along with the change of PZT position.Therefore, the position of PZT on the piezoelectric vibrator will have optimal solution ever if the influence of PZT on the piezoelectric vibrator is considered.
According to the conclusion above, the first three admissible functions in sequence of the function shown in (33) were selected to study the variation of the positional force coefficient of PZT.
The structure of the vibrator is shown in Figure 1 and the parameters are given in Table 4.The width of the vibrator is 16 mm.
Through calculation, it indicates that the positional force coefficient Π of the second-order bending vibration changes along with the changing of PZT position on the vibrator, as shown in Figure 5. From the changing curve, it can be seen that when  1 = 4 mm the value of Π reaches its minimum and PZT's incentive effect is lowest and when  1 = 21 mm the value of Π reaches its maximum, and then the PZT's incentive effect is the best.

Conclusions
In this paper the researches on the influence laws of PZT position and mode shape on performance of piezoelectric vibrator are conducted.And a new optimization method for the design of a piezoelectric laminated beam is proposed.The outcomes of this study are summarized as follows.
(1) After dividing and assembling the piezoelectric laminated vibrator with mutation sections, the dynamics equation is deduced based on the continuity boundary conditions of the vibrator.
(2) The form of vibration mode function was assumed based on the Hamilton principle and Rayleigh-Ritz method, and a sequence of admissible functions was chosen which met the displacement boundary conditions.
(3) Approximate functions of frequency and mode shape were calculated by giving different values to the item number  of the admissible function sequence, and the calculated results were compared with experimental results.It was concluded that with more items the approximate solution became more accurate.The cause of these errors was analyzed based on the selection of admissible function and the external electric field.
(4) Mode force of the piezoelectric laminated vibrator was analyzed with the action of PZT's inverse piezoelectric effect.The result indicates that mode force not only was associated with piezoelectric material properties  31 , voltage  of external electric field, and structural parameter but also related to the position of PZT relative to the mode shape curve of the vibrator.

Figure 3 :
Figure 3: The second-order vibration curve of the vibrator with the PZT pasted in the center of the substrate.

Figure 5 :
Figure 5: The variation curve of positional force coefficient Π along with the changing of PZT position on the vibrator.

Table 1 :
Properties of the materials used.

Table 2 :
Dimension parameters of the vibrator (mm).

Table 3 :
Calculated and experimental values of resonant frequency under the driving voltage of 60 V.

Table 4 :
Dimension parameters of the vibrator (mm).