# Damage Identification in Structures Based on Energy Curvature Difference of Wavelet Packet Transform

Authors: Pengbo Wang; Qinghe Shi
Journal: Shock and Vibration
 (2018)

---

## Abstract

Damage identification is of tremendous significance in engineering structures. One key issue in damage identification is to determine an index that is sensitive to the structural damage. Current damage identification indices are generally focused on dynamic characteristics such as the natural frequencies, modal shapes, frequency responses, or their mathematical combinations. In this study, based on the wavelet packet transform, we propose a novel index, the energy curvature difference (ECD) index, to identify the damage in structures. The ECD index is the summation of component energy curvature differences after a signal is decomposed using WPT. Moreover, two numerical examples are used to demonstrate the feasibility and validity of the proposed ECD index for damage identification. Stiffness reduction is employed to simulate the structural damage. The damage can be identified by the ECD index curve plot. The results of the examples indicate that the proposed ECD index is sensitive to low damage levels because even 5% stiffness reduction can be apparently identified. The proposed ECD index can be employed to effectively identify structural damage.


---

## Body

## 1. Introduction

Damage is a serious threat to structures during their service life, and a damage identification technology is of tremendous significance in engineering applications. Once a structure is damaged, the dynamic characteristics of the structure, such as the modal shapes, natural frequencies, and frequency responses, will change [1]. Various methods have been proposed for the use of dynamic responses for damage identification [2–6]. These methods are regarded as vibration-based damage identification methods. Among these methods, the approach based on wavelet transform is playing an increasingly important role because it can focus on any detail of a signal in the time or frequency domain [7, 8].The wavelet transform is an extension of the traditional Fourier transform and is capable of performing local and self-adaptive time-frequency analyses. Therefore, the wavelet transform can reveal some hidden phenomena of a signal that other signal processing techniques fail to observe [9, 10]. This property has been introduced in the damage identification application. Numerous studies have been published regarding the use of wavelet transform for damage identification [11–15]. However, an obvious disadvantage of the wavelet transform is that the frequency resolution in the higher frequency domain is not high. It may be difficult to differentiate a signal containing high-frequency components that are very close to each other. Thus, scholars developed the wavelet packet transform (WPT) technique.The WPT can be viewed as an extension of the traditional wavelet transform. It is capable of executing a complete level-by-level decomposition on a signal [16]. By performing the WPT on a signal, the signal can be decomposed into a series of wavelet packet components with a certain decomposition level, and the component energies can be obtained. The energy of the original signal is the summation of the component energies corresponding to different frequency bands. The WPT has been used for the damage identification of structures [17–20]. Wavelet packet energy can be obtained in different frequency bands after a signal is decomposed using WPT. The wavelet packet energy has been employed for the damage identification in ancient wood structures, beams, cable-stayed bridges, and so on [21–25]. Curvature difference can reflect the small changes in a function or a signal and has been introduced in damage identification together with wavelet packet energy. Some researchers have used curvature difference based on wavelet packet energy to locate the damage in ancient wood structures, concrete frames, and bridges [26–29].One key issue in damage identification is to determine an index that is sensitive to the structural damage. Current damage identification indices are generally focused on dynamic characteristics such as the natural frequencies, modal shapes, frequency responses, or their mathematical combinations. In this study, based on the WPT, we propose a novel index, the energy curvature difference (ECD) index, to identify the damage in structures. The ECD index is the summation of component energy curvature differences after a signal is decomposed using WPT. It takes into account the spatial distribution of the collected signals and is sensitive to low damage levels because even a 5% stiffness reduction can be apparently identified. Moreover, two numerical examples are used to demonstrate the feasibility and validity of the proposed ECD index for damage identification. Mutation on the ECD index curve can identify the damage accurately.The remainder of this paper is organised as follows. In Section2, some basics of the WPT are briefly introduced. In Section 3, an ECD index based on the WPT is proposed for damage identification. In Section 4, two numerical examples are provided to illustrate the applicability of the proposed ECD index. Section 5 presents the conclusions of this study.


## 2. Basics of Wavelet Packet Transform

The WPT is capable of accomplishing a complete level-by-level decomposition of the signals. Wavelet packetψj,ki(t) can be expressed as follows: (1)ψj,kit=2-j/2ψi2-jt-k,i,j∈Z+,k∈Z,where i indicates the modulation parameter, j indicates the scale parameter, and k indicates the translation parameter. The wavelet packet functions possess the orthogonality property, i.e., (2)ψj,kmt·ψj,knt=0,m≠n.The first wavelet function is the mother wavelet function, i.e.,(3)ψ1t=ψt.Wavelet functionsψi(t) can be calculated using the recursive equation (4)ψ2it=2∑k=-∞+∞hkψi2t-kψ2i+1t=2∑k=-∞+∞gkψi2t-k,where h(k) and g(k) are the quadrature mirror filters. h(k) is related to the scaling function and g(k) is related to the mother wavelet function.Assume thatf(t) is a time signal. If f(t) is decomposed to the jth level, it can be expressed as (5)ft=∑i=02j-1fjit,where fji represents the ith order wavelet packet component signal at level j. fji can be expressed by a linear combination of wavelet packet functions ψj,ki(t) in the following manner: (6)fjit=∑k=-∞+∞cj,kit·ψj,kit,where cj,ki(t) represents the wavelet packet coefficient, expressed as (7)cj,kit=∫-∞+∞ft·ψj,kitdt.WPT is a generalization of wavelet decomposition that offers a richer range of possibilities for signal analysis. In wavelet transform, a signal is split into an approximation and a detail. The approximation is then itself split into a second-level approximation and detail, and the process is repeated. In WPT, the details as well as the approximations can be split. This yields more different ways to encode the signal. The decomposition tree of WPT is plotted in Figure1, where A represents the approximations and D represents the details, and j denotes the decomposition level.Figure 1
Decomposition tree of WPT (j is the decomposition level).The energy of signalf(t) at the jth decomposition level is defined as (8)Ej=∫-∞+∞f2tdt=∑m=02j-1∑n=02j-1∫-∞+∞fjmt·fjntdt.Substituting (2) and (6) into (8) yields (9)Ej=∑i=02j-1Eji=∑i=02j-1∫-∞+∞fjit2dt,where Eji represents the signal energy stored in the ith order component signal at the jth level. Eji is regarded as the component energy of the WPT, and it can be expressed as follows: (10)Eji=∫-∞+∞fjit2dt.The physical implication of (10) is that the total energy stored in the original signal can be decomposed into a summation of the component energies, where each component corresponds to an individual frequency band.


## 3. Damage Identification Based on Energy Curvature Difference

One key issue of damage identification is to determine an index that is sensitive to structural damage. In this study, we propose an energy curvature difference (ECD) index to identify the damage in structures. The wavelet packet energy of a signal can be expressed as(11)E=Ej0+Ej1+Ej2+⋯+Eji+⋯+Ej2j-1,where i represents the order of the component signal and j represents the decomposition level. We use the node acceleration responses of the structures as the signals to be analysed. These acceleration signals can be conveniently obtained by using acceleration sensors.In practice, the curvature of a functionf(x) is usually obtained by using the second-order difference equation: (12)f′′xi=fxi+1-fxi/xi+1-xi-fxi-fxi-1/xi-xi-1xi+1-xi+xi-xi-1/2,where xi-1, xi, and xi+1 represent three different points of variable x.If the distance betweenxi-1 and xiequals the distance between xi and xi+1, i.e.,(13)xi-xi-1=xi+1-xi=h,then in this case, (12) can be simplified as (14)f′′xi=fxi+1-2fxi+fxi-1h2.The curvature is essentially a second-order difference of a function. Compared with the original function, the curvature is more sensitive to minor changes in the variable. This is an advantage of the curvature concept. Now, we introduce the concept of curvature into the wavelet packet energy for achieving damage identification. The component energy curvature is defined as(15)Ej,ki′′=Ej,k+1i-Ej,ki/dk,k+1-Ej,ki-Ej,k-1i/dk-1,kdk,k+1+dk-1,k/2,where Ej,ki denotes the energy related to the collected signal component fji(t) of the kth node and Ej,ki′′ represents the curvature of Ej,ki. Symbol dk,k+1 represents the distance between the node k and node k+1, and dk-1,k represents the distance between the node k-1 and node k. In particular, if dk,k+1 is equal to dk-1,k, i.e.,(16)dk,k+1=dk-1,k=d,then in this situation, (15) can be simplified as (17)Ej,ki′′=Ej,k+1i-2Ej,ki+Ej,k-1id2.From (17), it can be observed that, by using the curvature technique, the component energy change of a single node is transformed to the distribution change of the curvature values of three adjacent nodes. The spatial distribution of the component energies has been taken into account in (17). This technique can amplify the changes in a signal.The component energy curvature difference is defined as(18)ΔEj,ki′′=Ej,ki′′d-Ej,ki′′u,where subscript d denotes a damaged structure and subscript u denotes an undamaged structure.Now we perform a summation of the component energy curvature differences, and the summation is defined as the ECD index, expressed as(19)ECD=∑i=02j-1ΔEj,ki′′,where i is the order of the wavelet packet frequency band, j is the wavelet packet decomposition level, and k is the signal collection node number.By substituting (18) into (19), the ECD index can also be expressed as follows: (20)ECD=∑i=02j-1Ej,ki′′d-Ej,ki′′u.The numerical value of the ECD index is regarded as an indicator of the damage identification in structures. In practice, we collect the node acceleration response signals of the undamaged structure and damaged structure. Then, we perform WPT on the collected signals. With the calculation of the wavelet packet component energies, we can obtain the ECD indices. Subsequently, we can plot the ECD indices of different nodes in one figure. A mutation on the plot is indicative of damage in a structure. A sudden change in the ECD plot is indicative of a damage occurrence. The location and severity of the damage can be reflected by the ECD index plot.The proposed method requires signals of both damaged and undamaged structure, which sometimes could be difficult to obtain. In practice, with acceleration sensors installed appropriately, the signals of both damaged and undamaged structure are accessible. Besides, finite elements method is commonly used in engineering and we can simulate the damage in a fine meshed finite element model. The proposed method can give accurate results for a finite element model with damage. This will be very useful in the initial design stage of a structure. Thus, the proposed method has extent of engineering values.


## 4. Numerical Examples

To demonstrate the feasibility of the proposed ECD index for damage identification of structures, the following two numerical examples are considered. Both the location and severity of the damage are investigated with the proposed damage identification index.Example 1.
Consider a cantilever beam under the action of an impact force, as shown in Figure2. The length of the beam is L=1 m. The section moment of inertia of the beam is I = 1.6×10−7 m4. Young’s modulus of the material is E=210GPa, the mass density isρ=7850kg/m3, and Poisson’s ratio is ν=0.3. It is assumed that there is an impact force F acting at the endpoint of the beam in the vertical direction. The force-time history is plotted in Figure 3. The beam is discretized by 21 nodes and meshed with 20 elements. The node acceleration responses are calculated using the software ANSYS 15.0 and are regarded as the collected signals. The sampling frequency employed in this example is 105 Hz. The signals are decomposed to level 4 with Db10 wavelet, and 16 component energies are generated in total.
The element stiffness reduction can be used to simulate the structural damage. In practice, it is convenient to use Young’s modulus reduction to characterize damage. In this example, we use two test cases to investigate the application of the proposed ECD index in damage location identification, as presented in Table1. The locations of the damage elements in the finite element model are illustrated in Figures 4 and 6, and the corresponding ECD index curves are plotted in Figures 5 and 7, respectively. The damaged elements are displayed using a different color in the finite element model that it is employed for the verification of the proposed index.
We use four test cases to investigate the application of the proposed ECD index in damage severity identification, as presented in Table2. The damage level of specific elements increases from 5% to 20%. The ECD curves corresponding to all the four cases are plotted in Figure 8.
The decomposition tree of WPT in Example1 is plotted in Figure 9. The tree consists of several knots. Knot (j,i) represents the ith order wavelet packet component signal at decomposition level j.
For graphical illustration, the wavelet packet coefficients corresponding to the 16th node in case S1 of knots (1,0), (1,1), (2,0), (2,1), (3,0), (3,1), (4,0), and (4,1) are presented in Figure10.
From the plots in Figures5 and 7, it can be observed that the proposed ECD index can identify the damage locations accurately. There is an obvious mutation on the ECD index curve at the location where the damage is present. Moreover, the ECD index is sensitive to low damage levels because even 5% stiffness reduction can be apparently reflected by the ECD index curve. From the plot in Figure 8, it can be observed that the damage severity can be reflected by the amplitude of the mutation point. For the same damage locations with different damage levels, the amplitudes of the mutation are higher for the cases that are more severely damaged. The results of this example indicate that the proposed ECD index can be employed to identify the damage location and damage severity effectively.Table 1
Test cases for damage location identification in Example1.
Case No. Damage element Damage location Damage level Young’s modulus L1 Element 9 Between node 9 and node 10 5% 210 GPa×95% L2 Element 9 Between node 9 and node 10 5% 210 GPa×95% Element 16 Between node 16 and node 17 5% 210 GPa×95%Table 2
Test cases for damage severity identification in Example1.
Case No. Damage element Damage location Damage level Young’s modulus S1 Element 9 Between node 9 and node 10 5% 210 GPa×95% Element 16 Between node 16 and node 17 5% 210 GPa×95% S2 Element 9 Between node 9 and node 10 10% 210 GPa×90% Element 16 Between node 16 and node 17 10% 210 GPa×90% S3 Element 9 Between node 9 and node 10 15% 210 GPa×85% Element 16 Between node 16 and node 17 15% 210 GPa×85% S4 Element 9 Between node 9 and node 10 20% 210 GPa×80% Element 16 Between node 16 and node 17 20% 210 GPa×80%Figure 2
Cantilever beam under the action of an impact force.Figure 3
Force-time history curve for Example1.Figure 4
Location of damage element in case L1 in Example1.Figure 5
ECD index curve for damage location identification in case L1 in Example1.Figure 6
Locations of damage elements in case L2 in Example1.Figure 7
ECD index curve for damage location identification in case L2 in Example1.Figure 8
ECD index curve for damage severity identification in cases S1 to S4 in Example1.Figure 9
Decomposition tree of WPT in Example1.Figure 10
Wavelet packet coefficients corresponding to the 16th node in case S1 in Example1.
(a)
Knot (1,0) (b)
Knot (1,1) (c)
Knot (2,0) (d)
Knot (2,1) (e)
Knot (3,0) (f)
Knot (3,1) (g)
Knot (4,0) (h)
Knot (4,1)Example 2.
Now, we consider a spatial frame structure under the action of an impact force, as shown in Figure11. The dimensions are as follows: h1 = 1.6 m, h2 = 1.5 m, l = 1.8 m, and w = 0.9 m. The frame section is plotted in Figure 12. The section dimensions are b = 0.1 m and t = 0.009 m. Young’s modulus of the material is E=210 GPa, the mass density is ρ=7850kg/m3, and Poisson’s ratio is ν=0.3. It is assumed that there is an impact force acting at point D in the z direction. The force-time history is plotted in Figure 13. The structure is meshed with 760 beam elements. The node acceleration responses are calculated using the software ANSYS 15.0 and are regarded as the collected signals. The sampling frequency employed in this example is 2000 Hz. The signals are decomposed to level 4 with Db10 wavelet, and 16 component energies are generated in total.
It is assumed that the damage occurs in the EF part. In the finite element model, the EF part is discretized by 36 elements. The element close to point E is numbered 1, and the element close to point F is numbered 36. The length of each element is 0.05 m. There are totally 37 nodes located in the EF part. The node at point E is numbered as the 1st node, and the node at point F is numbered as the 37th node. The element numbers and node numbers in the EF part are illustrated in Figure14.
The element stiffness reduction can be used to simulate the structural damage. In practice, it is convenient to use Young’s modulus reduction to characterize damage. In this example, we use two test cases to investigate the application of the proposed ECD index in damage location identification, as presented in Table3. The locations of the damage elements in the finite element model are illustrated in Figures 15 and 17, and the corresponding ECD index curves are plotted in Figures 16 and 18, respectively. The damaged elements are displayed using a different color in the finite element model that it is employed for the verification of the proposed index.
We use four test cases to investigate the application of the proposed ECD index in damage severity identification, as presented in Table4. The damage level of specific elements increases from 5% to 20%. The ECD curves corresponding to all the four cases are plotted in Figure 19.
The decomposition tree of WPT in Example1 is plotted in Figure 20. The tree consists of several knots. Knot (j,i) represents the ith order wavelet packet component signal at decomposition level j.
For graphical illustration, the wavelet packet coefficients corresponding to the 24th node in case S1 of knots (1,0), (1,1), (2,0), (2,1), (3,0), (3,1), (4,0), and (4,1) are presented in Figure21.
From the plots in Figures16 and 18, it can be observed that the proposed ECD index can identify the damage locations accurately. There is an obvious mutation on the ECD index curve at the location where the damage is present. Moreover, the ECD index is sensitive to low damage levels because even 5% stiffness reduction can be apparently reflected by the ECD index curve. From the plot in Figure 19, it can be observed that the damage severity can be reflected by the amplitude of the mutation point. For the same damage locations with different damage levels, the amplitudes of the mutation are higher for the cases that are more severely damaged. The results of this example indicate that the proposed ECD index can be employed to identify the damage location and damage severity effectively.Table 3
Test cases for damage location identification in Example2.
Case No. Damage element Damage location Damage level Young’s modulus L1 Element 18 Between node 18 and node 19 5% 210 GPa×95% L2 Element 12 Between node 12 and node 13 5% 210 GPa×95% Element 24 Between node 24 and node 25 5% 210 GPa×95%Table 4
Test cases for damage severity identification in Example2.
Case No. Damage element Damage location Damage level Young’s modulus S1 Element 12 Between node 12 and node 13 5% 210 GPa×95% Element 24 Between node 24 and node 25 5% 210 GPa×95% S2 Element 12 Between node 12 and node 13 10% 210 GPa×90% Element 24 Between node 24 and node 25 10% 210 GPa×90% S3 Element 12 Between node 12 and node 13 15% 210 GPa×85% Element 24 Between node 24 and node 25 15% 210 GPa×85% S4 Element 12 Between node 12 and node 13 20% 210 GPa×80% Element 24 Between node 24 and node 25 20% 210 GPa×80%Figure 11
Spatial frame structure under the action of an impact force.Figure 12
Illustration of the frame section.Figure 13
Force-time history curve for Example2.Figure 14
Element numbers and node numbers of the EF part.Figure 15
Location of damage element in case L1 in Example2.Figure 16
ECD index curve for damage location identification in case L1 in Example2.Figure 17
Locations of damage elements in case L2 in Example2.Figure 18
ECD index curve for damage location identification in case L2 in Example2.Figure 19
ECD index curve for damage severity identification in cases S1 to S4 in Example2.Figure 20
Decomposition tree of WPT in Example2.Figure 21
Wavelet packet coefficients corresponding to the 24th node in case S1 in Example2.
(a)
Knot (1,0) (b)
Knot (1,1) (c)
Knot (2,0) (d)
Knot (2,1) (e)
Knot (3,0) (f)
Knot (3,1) (g)
Knot (4,0) (h)
Knot (4,1)


## 5. Conclusions

In this study, we propose an ECD index based on the WPT to identify a structural damage. The node acceleration responses of the undamaged and damaged structures are collected. We perform WPT on the signals and compute the component energy of each decomposed component signal. Then the ECD index can be obtained. With the illustration of the ECD index curve, the damage location and damage severity can be identified effectively. Two numerical examples are used to demonstrate the applicability of the proposed ECD index for damage identification. The ECD index is sensitive to low damage levels because even 5% stiffness reduction can be apparently identified by it. The mutation on the ECD index curve can identify the damage accurately, which is demonstrated by the illustration of damage elements in the finite element model. The results of this study indicate that the proposed ECD index can be employed to effectively identify the damage to structures for engineering applications.



---

*Source: 4830391-2018-05-28.xml*