Electromechanical Impedance Analysis on Piezoelectric Smart Beam with a Crack Based on Spectral Element Method

An electromechanical impedance (EMI) analysis of a piezoelectric smart beam with a crack is implemented in this paper. Spectral element method (SEM) is used to analyze the EMI response of the piezoelectric smart beam. In this analysis, the spectral element stiffness matrices of different beam segments are derived in this paper. The crack is simulated using spring models, and the EMI signatures of piezoelectric smart beamwith and without crack are calculated using SEM, respectively. From the analysis results, it is found that the peak position and amplitude of the EMI signatures have significant changes with the change in crack depth, especially in higher frequency ranges. Different vibration modes of the piezoelectric smart beam are analyzed, and the effect of thickness of the adhesive layer on the admittance is also researched. An experimental study is also implemented to verify the validity of the analysis results using SEM.


Introduction
The presence of structural damage will result in the changes in the properties of the structure and will also induce the change in the mechanical impedance of the structure.Damage identification can be achieved by comparing the structural mechanical impedances before and after damage.However, structural mechanical impedance is difficult to measure in practice.The coupled electromechanical impedance (EMI) of structure with PZT patches bonded to it is convenient to obtain by using impedance analyzer, even in the high frequency ranges, which is quite sensitive to structural incipient damage.In recent years, the EMI-based technique has been applied to structural health monitoring.It has been successfully used to analyze various engineering structures, including aerospace structures, mechanical structures, and civil structures [1][2][3][4][5][6][7][8][9][10][11].The EMI technique is based on the direct and converse piezoelectric effects of the PZT patch.In order to analyze the EMI of coupled structures, the model characterizing electromechanical interaction between the host structure and PZT patch should be firstly established.Liang et al. [12] proposed the first EMI model for a PZT patch bonded onto an intact one-dimensional structure.Zhou et al. [13] extended the impedance method to model a twodimensional EMI structure.Bhalla and Soh also presented an improved 2D impedance model to characterize the PZTstructure electroelastic interactions based on the concept of "effective impedance" [14].Annamdas and Soh have examined the three-dimensional interaction of a PZT transducer or multiple PZT transducers with the host structure based on directional sum impedance formulation [15,16].Wang et al. also proposed an embedded 3D electromechanical impedance model for an embedded PZT transducer by considering the interaction between a square PZT patch and a host structure [17].
Spectral element method (SEM) has been proposed to the structural dynamics community as a numerical tool to analyze the dynamic responses of rods, beam, and plates [18].Spectral element is represented by an exponential interpolation function, which is an exact solution of the wave equation for representative structural element.Therefore, other than classical finite element method, SEM provides dynamic stiffness matrix of a structure, which is more accurate and efficient than the traditional static stiffness matrix, and allows an increased accuracy in modeling the structural dynamic behavior.Subsequently, Palacz and Krawczuk [19] analyzed the dynamic responses of a cracked rod using SEM.The crack is simulated by a rotational spring model, and the node displacement and shape function of the cracked rod are derived.Ostachowicz et al. further summarized the researches on spectral element method and derived the spectral element stiffness matrices of Euler beam, Timoshenko beam, damaged plate, and thin-walled shell structures [20][21][22][23].Samaratunga et al. presented a new 2D wavelet spectral finite element model for studying wave propagation in thin to moderately thick anisotropic composite laminates [24].Choi and Inman recently presented modeling of a cableharnessed structure by means of SEM, and a double beam model was formulated.The presented modeling was applied and compared with the conventional FEM to emulate a cableharnessed structure [25].
Some researchers also focused on identifying the damage from the EMI signatures of the cracked beam with surfacebonded PZT patches using SEM. Park et al. implemented damage identification study on a simple one-dimensional structure by combining SEM and EMI techniques [26,27].Ritdumrongkul et al. succeeded in quantitative identification of the structure damage based on the SEM and PZT active sensors.In the experiment, a bolted aluminum beam is studied, and the damage is simulated by looseness of the bolt [28].Wang and Tang further deduced the stiffness matrix of smart Timoshenko beam and successfully analyzed electric admittances of piezoelectric smart beam using the SEM [29].Combined with EMI method and nonlinear optimization technique, the spectral element model of a simply supported rod was simulated and numerical study on localization and quantitative identification of damage was made by Guo and Sun [30].Ostachowicz et al. developed their own numerical procedures using SEM in order to calculate damage indexes, which were used both for damage detection and for localization.The proposed methods were applied for structural health monitoring of metallic and fibre reinforced structural elements [31].Kim and Wang proposed an improved impedance-based damage identification method by incorporating a tunable piezoelectric circuitry with the structure to enrich the impedance measurements.Numerical case study on localizing damage in a fixed-fixed beam using SEM was performed to demonstrate the effectiveness of the new method for structural damage identification [32].
In this paper, a new spectral element model of a cracked Timoshenko beam is proposed.In the model, the crack is modeled accurately by three massless springs, that is, a shear spring, a translational spring, and a rotational spring.According to the relationship between dynamic load and displacement, the spectral element stiffness matrices of different beam segments are obtained, and then the EMI of the piezoelectric coupled beam with a crack is formulized.Subsequently, the proposed spectral element model is validated by experimental results.Finally, numerical investigations are addressed to analyze the EMI responses of the piezoelectric Timoshenko beam with a crack based on the spectral element model.
The coefficients   ( = 1, . . ., 6) can be calculated as a function of the nodal spectral displacements using the boundary conditions.The nodal spectral forces can be determined by differentiating the spectral displacements with respect to .Finally, using relation between the nodal forces and nodal displacements the spectral element stiffness matrix of an intact Timoshenko beam can be obtained as shown in [18].

Translational
spring Shear spring

Spectral Element Stiffness Matrix of a Cracked Timoshenko
Beam.The spectral element stiffness matrices of beam segments with and without a crack are similar in form.A spectral beam element model with a transverse open and nonpropagating crack is presented as shown in Figure 1.The crack is simulated by three massless springs, namely, one rotational spring, one translational spring, and one shear spring in this paper, whose flexibilities are calculated using Castigliano's theorem and the laws of the fracture mechanics.The length of the beam element is  and the crack is located at  1 .
The spectral solutions for the left and right parts of the beam divided by the crack are similar to the form of ( 5), which can be expressed using a local coordinate system as Node displacements at the two ends of the spectral cracked beam element are The continuity conditions at the crack location are where  11 ,  13 ,  22 , and  33 are stiffness coefficientsof the three springs, respectively.Their values can be obtained by Castigliano's theorem [33]: where  11 ,  13 ,  31 , and  33 are the flexibility coefficients of the three springs, respectively, which can be expressed as where  is the width of the beam and  = ℎ  /ℎ is the ratio between the crack depth ℎ  and the height of the beam ℎ.
In addition, the boundary conditions of the free-free beam are expressed as  (1)  1 =  (1) After the above analysis, the stiffness matrix of the crack can be obtained.By substituting (7) into ( 8)-( 9), the spectral nodal displacements can be written in a matrix form as follows: where By substituting ( 7) into (4), the spectral nodal forces can be expressed in the matrix form as follows: where The spectral force-displacement relation can be derived by eliminating the coefficient matrix A in ( 13) and ( 14) as follows: K is the frequency-dependent spectral element stiffness matrix of the Timoshenko beam with a crack.Crossing out 4 to 9 columns, the spectral stiffness matrixcan be written as a square matrix (6 × 6).

Spectral Element Stiffness Matrix of Beam Segment
with Surface-Bonded PZT.The vibration properties of beam segment with surface-bonded PZT, which is a PZT-adhesive layer-beam coupled structural element, as shown in Figure 2, must be considered synthetically.Below are the equivalent parameters of the coupled structural element: where   ,   ,   , and   denote the equivalent cross-sectional area, density, elasticity modulus, and inertia moment of the coupled structure, respectively; , ℎ, , , and  are the width, height, density, elasticity modulus, and inertia moment of the base beam, respectively;   , ℎ  ,   ,   , and   are the width, height, density, elasticity modulus, and inertia moment of the adhesive layer, respectively;   , ℎ  ,   ,  11 , and   are the width, height, density, elasticity modulus, and inertia moment of the PZT patch, respectively.Substituting these equivalent parameters into (7), the displacement equations of the beam segment with surfacebonded PZT can be obtained.According to the analysis method of cracked beam segment, the spectral element stiffness matrix of the beam segment with surface-bonded PZT can be also derived, whose form is similar to the cracked beam.
As shown in Figure 2, a single PZT is bonded onto the surface of the base beam.The piezoelectric smart beam consists of three layers: piezoelectric layer, adhesive layer, and base beam.The PZT will produce longitudinal expansion and shrink under voltage excitation.There exist interaction forces between the three layers.The base beam will produce both axial and bending vibration when there is only a single PZT.  ,   , and   denote the forces and moments applied to the base beam, which are generated by the driving voltage., , and  are the corresponding displacements.These forces and moments can be expressed in the spectral form as follows: As shown in Figures 3 and 4, two PZT patches are attached to the top and bottom surfaces of the beam  symmetrically.The vibration in this case is different to the case that only a single PZT is boned onto the top surface of the base beam.
Under the same size and same phase of the electrical excitation, bending vibration caused by the two PZT patches will be offset; the base beam only produces axial vibration.In this case, forces and moments applied to the base beam are only 2  .Details of the forces can be seen in Figure 3.
On the contrary, when the PZT patches are under the same size but opposite phase of the electrical excitation, axial vibration caused by the two PZT patches will be offset; the base beam only produces bending vibration.In this case, forces and moments applied to the base beam are 2  and 2  , which can be concluded from Figure 4 clearly.

Formation of Global Stiffness Matrix.
The PZT-adhesive layer-cracked beam coupled structure can be divided into different segments due to the location of PZT patch.The spectral element stiffness matrices of beam segments with and without a crack are similar in form.The spectral element stiffness matrix of beam segments with PZT is also derived.So the stiffness matrix of the beam can be determined.
Continuity conditions between different beam segments can be expressed as where superscript "" denotes the number of the beam segments, subscript "1" denotes the left end of the beam segment, and subscript "2" denotes the right end of the beam segment.By assembling stiffness matrices of different beam segments through the above continuity conditions, the global spectral element stiffness matrix of the piezoelectric smart beam can be obtained.
In order to analyze the displacement of the piezoelectric smart beam under load, restraint conditions at the two ends need to be considered.When one end of the beam is free, to the left end as an example, there exists  1 =  1 =  1 = 0; when it is simply supported, then  1 =  1 =  1 = 0; when it is fixed, then  1 =  1 =  1 = 0.By taking a different combination of the above restraints, different boundary conditions of the piezoelectric smart beam can be simulated.

EMI Analysis on Piezoelectric Smart Beam
A single PZT patch is bonded onto the top surface of the piezoelectric beam.Relative to the base beam, the size of piezoelectric patch is very small.The PZT is considered as a one-dimensional model and the force generated by piezoelectric patches can be simplified as a pair of axial force.The constitutive equations of the PZT patch are where  1 ,  1 are the strain and stress, respectively;  11 =  11 (1 + ) is complex Young's modulus at zero electric field with  denoting the mechanical loss factor;  3 = ()/ℎ  is the electric field intensity;  3 is the electric displacement along  direction, respectively;  33 =  33 (1−) is the complex dielectric constant at zero stress with  denoting the dielectric loss factor;  31 is the piezoelectric strain constant.From (19), the electric displacement and the electric displacement can be defined as The electric current passing through the PZT patch can be obtained from the electric displacement as follows: where   denotes the length of the PZT patch; , , and  are the change of axial displacement, rotational angle, and shear displacement of two ends of the bonded-PZT beam, which can be expressed as The coupled EMI of the piezoelectric smart beam can be described as For the piezoelectric smart beam, the PZT measured impedance curve is the function of the dynamic characteristic of PZT patch and impedance value of the base beam.It can be found that the impedance of the piezoelectric smart beam has relation with parameters of the beam.When the beam is damaged, the impedance curve of the piezoelectric smart beam will change.So change of the electromechanical impedance of the beam can reflect damage of the beam.

Experimental Verification
In this study, a free-free Timoshenko beam with a single PZT patch attached to its top surface is investigated.The type of PZT patch is PZT5A, and the PZT is located at  1 = 136.5 mm,  2 = 163 mm from the left end of the beam, and the crack is located at  3 = 182 mm.Here we consider intact beam and three damage cases, which are corresponding to the relative crack depths of 10%, 20%, and 30%, respectively.The detailed information is as shown in Figure 5; a single PZT patch is bonded onto the top surfaces of a Timoshenko beam with a crack.The PZT-adhesive layer-cracked beam coupled structure is divided into three segments due to the location of PZT patch.The material constants of the beam, adhesive layer, and PZT patch are listed in Table 1.
The PZT will produce longitudinal expansion and shrink under 1 voltage excitation, and then the base beam will also produce vibration.The realimpedances of intact smart beam are calculated using SEM in two different frequency ranges of 10-30 kHz and 80-100 kHz.In order to verify the validity of the analysis results using SEM, some experimental results for intact smart beam are also used for comparison.Numerical and experimental real impedance curves of smart intact beam and cracked beam are shown in Figures 6 and 7, respectively.
As shown in Figures 6 and 7, numerical and experimental results are in agreement by and large, and the impedance magnitudes are within an order of magnitude.Parts of the wave locations do not match.There are many reasons for this case.There exists error usingfoam to simulate free boundary conditions.Numerical models of piezoelectric patch and adhesive layer are considered too simple, which cannot fully reflect the actual vibration conditions.From the analysis results, it is found that the peak position and amplitude of the electromechanical impedance curve have significant changes with the crack, especially in higher frequency ranges.The experimental results are coincident with those using SEM.Therefore, it can be concluded that the SEM can be used to analyze the electromechanical impedance responses of a piezoelectric smart beam effectively.

Different Damage Cases.
The crack is simulated by three massless springs, namely, one rotational spring, one translational spring, and one shear spring in this paper.In order to certify the effectiveness of the spring model, piezoelectric beam under different damage conditions is analyzed.Because the process of vibration with a single PZT is too complicated, and the number of axial vibration modes is relatively small, therefore, the bending mode is taken to analyze the effect of crack depth on the admittance of the beam in this paper.On the other hand, the peak of bending vibration mode is more significant than the other two modes, so it is more suitable to be applied in structural health monitoring.
The real impedances for piezoelectric smart beam with a crack are calculated using SEM in two different frequency ranges of 10-30 kHz and 80-100 kHz.Numerical EMI signatures of smart beam with different crack sizes are shown in Figure 8. From Figure 8 it is found that the peak positions and amplitudes of the real impedance curves have significant left shifting trend with the variations of crack size, especially in higher frequency ranges.This is because the extension of cracks reduced the overall stiffness of the beam.That is, the greater the damage is, the more the curve shifts to the left.Therefore, better results can be achieved using high frequency for damage identification in the health monitoring.

Influence of Adhesive
Layer.The impact of the adhesive layer has been considered when analyzing the impedance of the piezoelectric smart beam.When the thickness of the adhesive layer changes, the impedance of the piezoelectric beam will also change.In order to analyze the effect of thickness of adhesive layer on the impedance, free-free beam without crack is taken as the object.The cumulative admittance shifts (CAS) and the cross correlation coefficients  (CC) between two real admittance data are defined as where Re(  ) and Re( 0 ) denote a real part of electric admittance with and without an adhesive layer, respectively;  denotes the number of measurements; Re() represents the real impedance; subscript "0" represents initial structural status without adhesive layer; superscript " " represents the mean of a Re();  Re(  ) and  Re( 0 ) are the standard deviations of the real parts of   and  0 .As shown in Figure 9, with increasing thickness of the adhesive, the index CC is becoming smaller, and the change of correlation coefficients in the high frequency is more obvious than in the low frequency.When the adhesive layer thickness increases, force and moment applied to the beam become larger; the magnitude of the admittance will also increase.The cumulative admittance shifts (CAS) are used to analyze the change of magnitude of the admittance.The analysis results in two different frequency ranges of 10-30 kHz and 80-100 kHz can be seen in Figure 9(b).

Different Vibration Modes.
When the PZT patches are pasted with different forms, the beam will produce different vibration modes.When pasted with a single PZT patch, the beam will produce both bending vibration and axial vibration.When two PZT patches attach to the top and bottom surfaces of the beam symmetrically, the vibration is different to the vibration of beam with a single PZT bonded onto it.Detail forces of the three vibration modes of the beam have been analyzed in the previous text.Under the same size and same phase of the electrical excitation, the base beam only produces axial vibration.Under the same size but opposite phase of the electrical excitation, the base beam only produces bending vibration.The three graphs in Figure 10 describe the vibration of the base beam with a single PZT bonded onto it, the vibration of the beam with collocated PZT patches under the same phase of electrical excitation, and that under the opposite phase of electrical excitation, respectively.The parameters of the beam are the same as those shown in Table 1.
As shown in Figure 11, the number of modes causing vibration of the base beam with a single PZT is equal to the number of bending vibration plus the number of axial vibration with collocated PZT patches in frequency ranges of 10-30 kHz and 80-100 kHz.This is because the beam will produce both bending vibration and axial vibration when only a single PZT patch vibrates under the voltage   excitation; the modes reflect information of bending vibration and axial vibration.On the other hand, the number of modes caused by the bending vibration is more than that caused by the axial vibration.This is because the value of the area moment of inertia is much smaller than the sectional area, so in a certain frequency band the interval between bending modes is smaller than the axial vibration modes.

Smart Beam with Dual PZT Patches.
In addition to the above analysis of a free-free beam, another beam with different damage conditions and patch pasted forms can be also analyzed.As shown in Figure 12, a cantilever Timoshenko beam with two PZT patches attached to its top surface is investigated.The PZT-adhesive layer-cracked beam coupled structure is divided into five segments due to the location of PZT patch.The length, width, and height of the piezoelectric smart beam are 400 mm, 20 mm, and 15 mm, respectively.PZT1 is located at  1 = 73.5 mm and  2 = 100 mm from the left end of the beam, PZT2 is located at  1 = 300 mm and  2 = 326.5 mm from the left end of the beam, and the crack is located at  3 = 150 mm.Here we consider intact beam and three damage cases, which are corresponding to the relative crack depths of 10%, 20%, and 30%, respectively.
As shown in Figures 13 and 14, the real impedance curves of PZT1 and PZT2 have significant left shifting trend with the variations of crack size, and the left shifting trend of impedance curves of PZT1 is more significant than that of PZT2.That is because the position of PZT1 is closer to crack than PZT2, so the results are more sensitive to damage.

Conclusions
An impedance analysis of a piezoelectric smart beam with a crack and a single PZT patch attached to it is implemented in this paper.The crack is simulated by spring models, and the influence of adhesive layer is also considered.
Different vibration modes of the piezoelectric beam are analyzed in this paper, and the effect of thickness of the adhesive layer on the admittance is also researched.In order to depict the effects of the depth of crack and thickness of adhesive layer on impedance curves of PZT patches, the damage index, cumulative admittance shifts (CAS), and cross correlation coefficient (CC) are defined.An experiment of a cracked beam with free boundary conditions at both ends is presented to verify the validity of the analysis results using SEM.The experimental and numerical results are in agreement by and large.It illustrates that using SEM to analyze the admittance of piezoelectric smart beam is feasible.Structural health monitoring technology based on electromechanical impedance using SEM has a certain application value.

Figure 1 :
Figure 1: Spectral beam element model of cracked beam segment.

Figure 3 :
Figure 3: Smart beam element with collocated PZT patches under the same phase of electrical excitation.

Figure 4 :
Figure 4: Smart beam element with collocated PZT patches under the opposite phase of electrical excitation.

Figure 5 :
Figure 5: Piezoelectric smart beam with a crack and a single PZT patch.

Figure 9 :Figure 10 :Figure 11 :
Figure 9: CC and CAS values of the piezoelectric smart beam in different frequency ranges: (a) CC index and (b) CAS index.

2 Figure 12 :Figure 13 :
Figure 12: Piezoelectric smart beam with a crack and two PZT patches.

Figure 14 :
Figure 14: CC of two PZT patches with different crack sizes in the frequency range of 10-30 kHz.