Piezoelectric Electro-Mechanical Impedance (EMI) Based Structural Crack Monitoring

: To detect small cracks in plate like structures, the high frequency characteristics of local dynamics were studied with the piezoelectric electro-mechanical impedance (EMI) method, and damages were monitored by the changes of the EMI. The ﬁnite element simulation model of EMI was established, and numerical analysis was conducted. The simulation results indicated that the peak frequency of the piezoelectric admittance signal is a certain order resonance frequency of the structure, and the piezoelectric impedance method could e ﬀ ectively detect the dynamic characteristics of the structure. The piezoelectric admittance simulation and experimental study of aluminum beams with di ﬀ erent crack sizes were performed. Simulation and experimental results revealed that the peak admittance frequency decreases with the increase of crack size, and the higher resonance frequency is more sensitive to the small-scale damage. The proposed method has good repeatability and strong signal-to-noise ratio to monitor the occurrence and development of small-scale crack damage, and it has an important application prospect.


Introduction
During the life cycle of some engineering structures, cracks will occur due to the effects of dynamic load [1,2], varying temperature, environmental corrosion [3][4][5] and other factors [6]. If the cracks are not detected and restrained after the initial occurrence, the crack growth will accelerate, leading to eventual structural failure with potentially significant losses. To avoid catastrophic consequences caused by cracks in structures, it is of particularly importance to detect initial cracks in these structures. At present, the techniques of structural crack detection include magnetic particle detection, ultrasonic detection and radiographic detection, among others [7,8]. Most of the above methods have their shortcomings, such as limited applicability, low detection efficiency and inability to detect crack growth in real time.
Compared with the traditional detection methods mentioned, structural health monitoring (SHM) [9][10][11] is advancing rapidly because of its real time capacity. In SHM, vibration-based structural damage detection is a widely used technology. Since for any structure, the structural dynamic characteristic is a function of its structural parameters such as stiffness, mass and damping, structural damage means the change of structural parameters, and the change of the structural parameters will inevitably lead to the change of structural dynamic characteristics, such as resonance frequency [12,13] technique combined with guided wave, Sohn and Park [48] utilized the piezoelectric arrays to conduct bridge damage detection, and proposed a symmetric index (SYMI) to evaluate the damage degree.
The piezoelectric impedance method is an active sensing method that uses only one transducer and is applied in the bolt looseness detection [49], stress monitoring [50] and structural damage diagnosis [51]. It is sensitive to structural small-scale damage due to its high frequency detection band, and can realize on-line monitoring under certain conditions. Zagrai et al. [52] verified that the electro-mechanical impedance method can identify the presence of damage in a thin plate by theoretical and experimental study. Xu and Giurgiutiu [53] presented a low-cost impedance analyzer for active structural health monitoring, and succeeded detecting a disbonding on a spacecraft panel of aluminum. Bhalla et al. [54] and Yang et al. [55] interrogated the damage condition of the structure by extracting the changes of equivalent mass, stiffness and damping of the structure through the piezoelectric impedance signature before and after damage, establishing a damage interrogation method. Ai et al. [56] applied the piezoelectric impedance technology to monitor the corrosion process of metals, and the mechanical impedance of the host structure was extracted from the electromechanical coupling impedance by the united mechanical impedances (UMI) method. It found that the mechanical impedance damage index extracted by the proposed UMI method was better than that by the traditional method.
In this paper, the piezoelectric electro-mechanical impedance (EMI) was adopted for the detection of the small crack in a beam. By analyzing the piezoelectric admittance formula of the electromechanical coupling system and combining with the harmonic response analysis of the finite element method (FEM), it found that the local dynamic characteristics of the structure would change when cracks occur and propagate in the beam. This dynamic variation is represented by the shift of peak frequency in the imaginary part of the piezoelectric admittance spectrum. The method of piezoelectric admittance peak frequency shift was proposed for the detection of the small-scale crack detection with a relatively higher monitoring frequency band. An experimental device was constructed to measure the admittance signal of a piezoelectric patch installed on the structure by a precision impedance meter. The relationship between the peak frequency in the piezoelectric admittance spectrum and the change of the crack size was analyzed. The finite element simulation and experimental results show that the small crack in the beam can be detected by the proposed method.

Principle of Piezoelectric Electro-Mechanical Impedance Method
Because of its strong piezoelectric effect, a PZT (lead-zirconate-titanate) type piezoelectric transducer has been successfully used in SHM and is adopted in this research. For a PZT patch of length l a , width b a , and thickness h a is bonded on a structure, based on the assumption of elastic connection (represented by a pair of stiffness springs 2K str ) between piezoelectric material and host structure, as shown in Figure 1, Giurgiutiu and Zagrai [57] proposed an EMI model and the piezoelectric admittance expression is given as: where ω is the angular frequency of excitation signal, ρ is the density of the PZT patch, κ 2 31 is the electromechanical coupling factor, ϕ = γ · l a /2, while γ = ω/c is the wave number, c = 1/ρs E 11 is the wave speed, s E 11 is the mechanical compliance at zero field, C = l a b a ε T 33 h −1 a is the conventional capacitance of the PZT patch, K str is the stiffness of the host structure, K PZT = A a / s E 11 l a is the stiffness of the PZT patch, A a = h a × b a is the cross-sectional area of PZT patch, and ε T 33 is dielectric constant at zero stress. The expressions in Equations (1) cover the complete frequency spectrum and encompass both structure and sensor dynamics, overcoming low-frequency structure-focused analysis of Liang's model [58]. connection (represented by a pair of stiffness springs 2Kstr) between piezoelectric material and host structure, as shown in Figure 1, Giurgiutiu and Zagrai [57] proposed an EMI model and the piezoelectric admittance expression is given as: Equation (1) reveals that the change of admittance Y of the coupled system is only related on the host structural stiffness K str while the parameters of piezoelectric material and host structural material are determined. In order to find out the effect of the stiffness of the host structure on the admittance signal, Equation (1) is derived with ω simultaneously on both sides. Here, we expand the cotangent function by a Taylor series and substitute the first two terms in the Taylor series into Equation (1): Some peak values in the admittance signature represent the resonant frequency of the host structure. In this case, the peak frequency is a certain order resonance frequency, ω n , of the coupled system of the host structure and the piezoelectric patch. Thus, in the position of the resonance frequency, ω n , the derivative of piezoelectric admittance signal equals zero, i.e., dY/dω n = 0. In Equation (2), both c and K PZT are complex.
Considering that the imaginary part of c and K PZT depend on the material damping and the effect of material damping is small, the imaginary part is neglected in order to simplify the analysis process. Meanwhile, let κ 2 31 = a + bi in Equation (2). By extracting the imaginary part of Equation (2) and making it equal to zero, the expressions about K str and ω n can be obtained as Equation (4) shows that the peak frequency ω n in the piezoelectric admittance signature depends on the stiffness of the host structure. When a crack occurs and propagates in the host structure, the local stiffness of the host structure will reduce. While the local stiffness of the host structure decreases, the peak frequency in the piezoelectric admittance signature also decreases [59]. The piezoelectric patch bonded near the crack can transform this stiffness variation into peak frequency changes in the piezoelectric admittance signature. By analyzing the peak frequency changes in the piezoelectric admittance signature, the structural crack can be detected and thus it provides a theoretical basis for monitoring the crack status of a host structure with the peak frequency changes of the piezoelectric impedance.

Finite Element Analysis of Piezoelectric Impedance
To verify the above proposed method, the finite element analysis (FEA) is adopted to simulate the piezoelectric impedance and to find out the influence of crack and crack size on the peak frequency in the piezoelectric admittance signature, i.e., the local resonance frequency of the coupled system. The host structure is an aluminum alloy beam with a dimension of 250 × 30 × 5 mm. A PZT patch, with a size of 8 × 7 × 1 mm, is bonded on the center of the beam surface with epoxy resin adhesive. The adhesive layer thickness is kept 0.125 mm. The properties of the piezoelectric patch are listed in Table 1 [60], and the material constants of the adhesive layer and the aluminum alloy beam are listed in Table 2. Cracks are simulated by setting rectangular grooves on a beam. The width of each rectangular groove is 0.5 and the depth is 1 mm. Only the length of groove is changed to simulate the different crack sizes. The length of the groove is set as 0, 6 and 12 mm, respectively. The location of the piezoelectric patch and the crack on the beam is shown in Figure 2. In the numerical analysis, free boundaries are set at both ends of the beam. The system damping ratio is set to 0.005.
Appl. Sci. 2020, 10, x 5 of 14 frequency in the piezoelectric admittance signature, i.e., the local resonance frequency of the coupled system. The host structure is an aluminum alloy beam with a dimension of 250 × 30 × 5 mm. A PZT patch, with a size of 8 × 7 × 1 mm, is bonded on the center of the beam surface with epoxy resin adhesive. The adhesive layer thickness is kept 0.125 mm. The properties of the piezoelectric patch are listed in Table 1 [60], and the material constants of the adhesive layer and the aluminum alloy beam are listed in Table 2. Cracks are simulated by setting rectangular grooves on a beam. The width of each rectangular groove is 0.5 and the depth is 1 mm. Only the length of groove is changed to simulate the different crack sizes. The length of the groove is set as 0, 6 and 12 mm, respectively. The location of the piezoelectric patch and the crack on the beam is shown in Figure 2. In the numerical analysis, free boundaries are set at both ends of the beam. The system damping ratio is set to 0.005.  The finite element model of the cracked beam with piezoelectric patch is built, as shown in Figure 3a. The model is meshed by a mapping grid with an element size of 2 mm, except where the crack generates and the piezoelectric patch is bonded. For the region where the crack is located, the mapping mesh with the element size of 0.25 mm is used. In the vicinity of the crack, the free mesh is used to gradually transfer the specialized element size of 0.25 mm to the overall model element size of 2 mm. The mesh at the crack is detailed in Figure 3b. The piezoelectric patch is discretized with an element size of 0.5 mm by mapping, and the adhesive layer under the piezoelectric patch is discretized by free mesh with a 1 mm element size gradually changing to a 2 mm element size.  The finite element model of the cracked beam with piezoelectric patch is built, as shown in Figure 3a. The model is meshed by a mapping grid with an element size of 2 mm, except where the crack generates and the piezoelectric patch is bonded. For the region where the crack is located, the mapping mesh with the element size of 0.25 mm is used. In the vicinity of the crack, the free mesh is used to gradually transfer the specialized element size of 0.25 mm to the overall model element size of 2 mm. The mesh at the crack is detailed in Figure 3b. The piezoelectric patch is discretized with an element size of 0.5 mm by mapping, and the adhesive layer under the piezoelectric patch is discretized by free mesh with a 1 mm element size gradually changing to a 2 mm element size.
Appl. Sci. 2020, 10, x 5 of 14 frequency in the piezoelectric admittance signature, i.e., the local resonance frequency of the coupled system. The host structure is an aluminum alloy beam with a dimension of 250 × 30 × 5 mm. A PZT patch, with a size of 8 × 7 × 1 mm, is bonded on the center of the beam surface with epoxy resin adhesive. The adhesive layer thickness is kept 0.125 mm. The properties of the piezoelectric patch are listed in Table 1 [60], and the material constants of the adhesive layer and the aluminum alloy beam are listed in Table 2. Cracks are simulated by setting rectangular grooves on a beam. The width of each rectangular groove is 0.5 and the depth is 1 mm. Only the length of groove is changed to simulate the different crack sizes. The length of the groove is set as 0, 6 and 12 mm, respectively. The location of the piezoelectric patch and the crack on the beam is shown in Figure 2. In the numerical analysis, free boundaries are set at both ends of the beam. The system damping ratio is set to 0.005.  The finite element model of the cracked beam with piezoelectric patch is built, as shown in Figure 3a. The model is meshed by a mapping grid with an element size of 2 mm, except where the crack generates and the piezoelectric patch is bonded. For the region where the crack is located, the mapping mesh with the element size of 0.25 mm is used. In the vicinity of the crack, the free mesh is used to gradually transfer the specialized element size of 0.25 mm to the overall model element size of 2 mm. The mesh at the crack is detailed in Figure 3b. The piezoelectric patch is discretized with an element size of 0.5 mm by mapping, and the adhesive layer under the piezoelectric patch is discretized by free mesh with a 1 mm element size gradually changing to a 2 mm element size.   In the simulation, the modal analysis is firstly carried out to obtain the resonance frequency, and then the piezoelectric harmonic analysis is performed. In the piezoelectric harmonic analysis, Appl. Sci. 2020, 10, 4648 6 of 14 the top and bottom surface piezoelectric degrees of the piezoelectric patch is coupled, and the voltage of 1 V and 0 V are applied on the top and bottom surfaces of the piezoelectric patch respectively to simulate 1 V AC voltage. The selected frequency range in the harmonic analysis is near the modal resonant frequencies obtained from the modal analysis. The time-varying charge on the piezoelectric surface is obtained, and then the piezoelectric impedance, piezoelectric admittance and other information are calculated. Comparing the results of the modal analysis with the piezoelectric admittance analysis, it is found that the modal resonant frequencies correspond approximately to the peak frequencies in the admittance signature, indicating that some of the peak admittance frequencies are the resonant frequencies of the coupled system. The peak frequencies in the admittance signature can reflect the local dynamic characteristics of the host structure with the piezoelectric patch.
After comparative analysis, the piezoelectric admittance around some resonant frequencies is simulated and analyzed with different crack sizes. The scanning frequency bands are selected around 19, 62 and 174 kHz, which is near the coupled structure resonance frequencies.
The detailed frequency bands of 19.2-19.3, 62.5-63 and 174.5-175 kHz, covering from relatively low to relatively high frequencies, are chosen as scanning frequency bands for the finite element harmonic response analysis. Then the piezoelectric impedance is computed. The piezoelectric admittance spectra and peak frequencies in the admittance spectrum with different cracks lengths (0 mm crack length means primitive state) on the beam are shown from Figures 4-6, the vertical coordinate is an imaginary part of the piezoelectric admittance spectrum, whose unit is siemens (S).
Appl. Sci. 2020, 10, x 6 of 14 In the simulation, the modal analysis is firstly carried out to obtain the resonance frequency, and then the piezoelectric harmonic analysis is performed. In the piezoelectric harmonic analysis, the top and bottom surface piezoelectric degrees of the piezoelectric patch is coupled, and the voltage of 1 V and 0 V are applied on the top and bottom surfaces of the piezoelectric patch respectively to simulate 1 V AC voltage. The selected frequency range in the harmonic analysis is near the modal resonant frequencies obtained from the modal analysis. The time-varying charge on the piezoelectric surface is obtained, and then the piezoelectric impedance, piezoelectric admittance and other information are calculated. Comparing the results of the modal analysis with the piezoelectric admittance analysis, it is found that the modal resonant frequencies correspond approximately to the peak frequencies in the admittance signature, indicating that some of the peak admittance frequencies are the resonant frequencies of the coupled system. The peak frequencies in the admittance signature can reflect the local dynamic characteristics of the host structure with the piezoelectric patch.
After comparative analysis, the piezoelectric admittance around some resonant frequencies is simulated and analyzed with different crack sizes. The scanning frequency bands are selected around 19, 62 and 174 kHz, which is near the coupled structure resonance frequencies.
The detailed frequency bands of 19.2-19.3, 62.5-63 and 174.5-175 kHz, covering from relatively low to relatively high frequencies, are chosen as scanning frequency bands for the finite element harmonic response analysis. Then the piezoelectric impedance is computed. The piezoelectric admittance spectra and peak frequencies in the admittance spectrum with different cracks lengths (0 mm crack length means primitive state) on the beam are shown from Figures 4-6, the vertical coordinate is an imaginary part of the piezoelectric admittance spectrum, whose unit is siemens (S).     From Figures 4a, 5a and 6a, it is clear that as the crack appears and grows, in both selected frequency bands, the peak frequencies in the piezoelectric admittance spectra shift to the left. It means that the peak frequency decreases when the crack grows, as shown in Figures 4b, 5b and 6b  From Figures 4a, 5a and 6a, it is clear that as the crack appears and grows, in both selected frequency bands, the peak frequencies in the piezoelectric admittance spectra shift to the left. It means that the peak frequency decreases when the crack grows, as shown in Figures 4b, 5b and 6b  From Figures 4a, 5a and 6a, it is clear that as the crack appears and grows, in both selected frequency bands, the peak frequencies in the piezoelectric admittance spectra shift to the left. It means that the peak frequency decreases when the crack grows, as shown in Figures 4b, 5b and 6b. Based on the above theoretical analysis, the local stiffness of the host structure around the crack decreases when the crack occurs and propagates. Once the local stiffness decreases, the resonant frequency of the host structure will decrease, shifting the peak frequency to left in the piezoelectric admittance spectrum. The numerical simulation results are consistent with the theoretical analysis.
Based on the above findings, the peak frequency in the piezoelectric admittance is chosen as the feature parameter to characterize the crack propagation. The peak frequency shift between crack-free and 6 mm crack is defined as ∆f 1 , while the peak frequency offset between crack-free and 12 mm crack is defined as ∆f 2 . The peak frequency shifts at different frequency bands with different crack sizes are extracted as shown in Table 3. It shows that the higher the scanning frequency band, the greater the peak frequency shift, thus the higher the frequency band, the higher the sensitivity of the proposed method.

Experiments and Discussion
To verify the above theoretical model and finite element analysis, an experimental device was set up, as shown in Figure 7, and experiments were conducted. The dimension of the test specimen was consistent with the simulation model and the piezoelectric patch was bonded on the center of the one surface of the beam, as shown in Figure 8. When the piezoelectric patch was bonded to the host structure with epoxy resin, two optical fibers of 0.125 mm diameter were placed parallel on the surface of the beam with the epoxy resin to ensure the consistent adhesive layer thickness for all test pieces. First, the PZT patch, bonded on the without crack beam, was performed the piezoelectric impedance test by the WK6530B impedance analyzer (Wayne Kerr, West Sussex, UK). Then, a crack with length of 6 mm was machined and the test was performed again. At last, the piezoelectric impedance test was carried out with crack length of 12 mm machined on the beam. The piezoelectric admittance spectra of the beam with different crack lengths (0, 6 and 12 mm) were measured and the peak frequencies in the admittance spectra are extracted, as shown in Figures 9-11. As the figures indicate, as the crack grows, the peak frequency, i.e., the resonant frequency of the coupled system shifts to the left. That means the frequency decreases as the crack grows. There may be a slight difference between the simulation and experimental results, which may be resulted from an error between the actual and the designed positions of the damage and PZT. The experimental results have the same trend with the theoretical analysis and finite element simulation results. The tested resonant frequencies are also consistent with the simulation results; therefore, the finite element simulation results can be used to guide the experiment to select an appropriate scanning frequency band of the piezoelectric impedance measurement.
Appl. Sci. 2020, 10, x 8 of 14 on the above theoretical analysis, the local stiffness of the host structure around the crack decreases when the crack occurs and propagates. Once the local stiffness decreases, the resonant frequency of the host structure will decrease, shifting the peak frequency to left in the piezoelectric admittance spectrum. The numerical simulation results are consistent with the theoretical analysis. Based on the above findings, the peak frequency in the piezoelectric admittance is chosen as the feature parameter to characterize the crack propagation. The peak frequency shift between crack-free and 6 mm crack is defined as Δf1, while the peak frequency offset between crack-free and 12 mm crack is defined as Δf2. The peak frequency shifts at different frequency bands with different crack sizes are extracted as shown in Table 3. It shows that the higher the scanning frequency band, the greater the peak frequency shift, thus the higher the frequency band, the higher the sensitivity of the proposed method.

Experiments and Discussion
To verify the above theoretical model and finite element analysis, an experimental device was set up, as shown in Figure 7, and experiments were conducted. The dimension of the test specimen was consistent with the simulation model and the piezoelectric patch was bonded on the center of the one surface of the beam, as shown in Figure 8. When the piezoelectric patch was bonded to the host structure with epoxy resin, two optical fibers of 0.125 mm diameter were placed parallel on the surface of the beam with the epoxy resin to ensure the consistent adhesive layer thickness for all test pieces. First, the PZT patch, bonded on the without crack beam, was performed the piezoelectric impedance test by the WK6530B impedance analyzer (Wayne Kerr, West Sussex, UK). Then, a crack with length of 6 mm was machined and the test was performed again. At last, the piezoelectric impedance test was carried out with crack length of 12 mm machined on the beam. The piezoelectric admittance spectra of the beam with different crack lengths (0, 6 and 12 mm) were measured and the peak frequencies in the admittance spectra are extracted, as shown in Figures 9-11. As the figures indicate, as the crack grows, the peak frequency, i.e., the resonant frequency of the coupled system shifts to the left. That means the frequency decreases as the crack grows. There may be a slight difference between the simulation and experimental results, which may be resulted from an error between the actual and the designed positions of the damage and PZT. The experimental results have the same trend with the theoretical analysis and finite element simulation results. The tested resonant frequencies are also consistent with the simulation results; therefore, the finite element simulation results can be used to guide the experiment to select an appropriate scanning frequency band of the piezoelectric impedance measurement.       Hz Figure 10. Experimental results (a) piezoelectric admittance spectrum and (b) peak frequency tendencies in 62.2 kHz-62.6 Hz with different crack sizes.
(b) Figure 10. Experimental results (a) piezoelectric admittance spectrum and (b) peak frequency tendencies in 62.2 kHz-62.6 Hz with different crack sizes.
(a) (b) Figure 11. Experimental results (a) piezoelectric admittance spectrum and (b) peak frequency tendencies in 174.1 kHz-174.6 kHz with different crack sizes.
The peak frequency shifts extracted from the piezoelectric admittance spectra at different frequency bands are shown in Table 4. The peak frequencies in the piezoelectric admittance spectra decrease with the increase of crack length in each frequency band. Moreover, at the same crack length, the peak frequency shifts increase with the increase of the scanning frequency band, which means increasing the scanning frequency band can increase the sensitivity of the proposed method. Therefore, the crack occurrence and propagation can be interrogated by the peak frequency shift in the imaginary part of the piezoelectric admittance spectrum. With the increase of the scanning frequency, the higher the frequency band is, and the larger the shift of peak frequency is, i.e., the higher the frequency band is, the more sensitive is the method to damage. Peak Frequency /kHz Crack size Figure 11. Experimental results (a) piezoelectric admittance spectrum and (b) peak frequency tendencies in 174.1 kHz-174.6 kHz with different crack sizes.
The peak frequency shifts extracted from the piezoelectric admittance spectra at different frequency bands are shown in Table 4. The peak frequencies in the piezoelectric admittance spectra decrease with the increase of crack length in each frequency band. Moreover, at the same crack length, the peak frequency shifts increase with the increase of the scanning frequency band, which means increasing the scanning frequency band can increase the sensitivity of the proposed method. Therefore, the crack occurrence and propagation can be interrogated by the peak frequency shift in the imaginary part of the piezoelectric admittance spectrum. With the increase of the scanning frequency, the higher the frequency band is, and the larger the shift of peak frequency is, i.e., the higher the frequency band is, the more sensitive is the method to damage. The experimental studies reveal that the peak frequency in the piezoelectric admittance spectrum, i.e., the resonant frequency of the coupled system of the beam and piezoelectric patch, decreases with the increase of crack length, which indicates that the peak frequency shift can well reflect the size of cracks. Therefore, the peak frequency can be used as a feature parameter to characterize the crack length. With the increase of the scanning frequency band, the shift of peak frequency tends to increase. In application, the detection sensitivity can be improved by properly increasing the monitoring frequency band, and small cracks can be detected and crack propagation can be monitored.
At the same time, the experimental results are basically consistent with the finite element numerical analysis results, which shows the accuracy and validity of the finite element model. That indicates the established finite element model can effectively help choosing the sensitive frequency band and the installation location of the piezoelectric patch in practical applications.

Conclusions
For small crack detection and crack propagation monitoring in structures, the method based on the peak frequency shift in the piezoelectric admittance spectrum is proposed. The feasibility of the method is analyzed theoretically. The FEA and experimental results verified the effectiveness of the method. The following conclusions are drawn: (1) The peak frequency in the imaginary part of the piezoelectric admittance spectrum can sensitively detect a crack in a structure, and the peak frequency decreases with the increase of the crack length.
(2) Appropriately increasing the detection frequency band can increase the peak frequency shift at the same crack size and improve the detection sensitivity to the crack, therefore, making the detection of a small crack in a structure possible.
(3) The experimental results are consistent with the FEM simulation results, which show the accuracy and validity of the finite element model. In future research, the resonant frequencies of complex structures can be found by a finite element model firstly, then the scanning frequency bands of piezoelectric admittance can be set to include the resonant frequencies, and finally the damage degree of structures can be judged by the peak frequency shift in the piezoelectric admittance signature.
(4) Since the frequency is selected as a feature parameter for damage monitoring, the proposed method has a better repeatability and signal-to-noise ratio than the root mean square deviation based damage index.
Since the method proposed in this paper is based on the change of local dynamic characteristics to monitor the small crack in a structure, there still needs to be more in-depth study on the local dynamic modeling and analysis. Further research is still needed on how to extract the feature parameter to characterize the crack propagation quantitatively, determine the optimal detection frequency band, and select the sensitive mode frequency bands, etc. Moreover, the influence of the damage size including the width and depth, the damage position, different excitation positions, the structural load and the damage identification of composite materials will be investigated in-depth in our future work.