You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
Department of Engineering Mechanics, Northwestern Polytechnical University, China
2
3
4
==Abstract==
5
6
A new damage feature index is presented for the structural health monitoring based on Hilbert-Huang transform (HHT). The energy marginal spectrum of the dynamic signal is used to construct damage characteristic parameter, which can reflect the signal energy variation and benefit the structural damage detection. A sinusoidal wave with frequency change and a composite plate vibration experiment with pre-defined damage are designed to verify the effectiveness of characteristic parameter in damage detection. Results obtained from simulation and test show that the extracted non-model-based damage feature index is available and sensitive in damage detection of time-varying system.
7
8
==1.  Introduction==
9
10
Structural damage induces the change of structural mechanical properties which can affect the stability and lifetime of structures in service. This is why the structural health monitoring (SHM) technology, especially for the critical unit in structure with large scale, has attracted an increasing attention in recent years. In addition, the composites need to pay more attention since the structures are extensively used in engineering field while prone to be damaged [1].
11
12
Damage detection technology is the major part in SHM which can early warn in the initial damage condition to avoid failure of structure, in which the signal processing method is playing an increasingly important role. Among the signal processing methods, the wavelet transform [2-4] and the Hilbert-Huang transform (HHT)[5, 6] are the most typical representatives. W.L. Bayissa [7] developed a structural damage method based on energy density function using wavelet transform. However, the wavelet transform relies on the choice of wavelet transform which limits the application in engineering.
13
14
HHT is composed of empirical mode decomposition (EMD) and Hilbert transform. Compared with the Fourier transform and the wavelet transform, EMD in HHT does not need prior functions [8]. Owing to the adaptive nature, HHT can be used in analyzing nonlinear and non-stationary signals [9]. R. Jha [10] processed a multi-level structural health monitoring based on HHT and succeeded in detecting and locating the damage. C.X. Bao [11] improved the HHT algorithm to identify the time-varying systems. The results of simulation and composite beam test showed that the improved algorithm was effective for the analysis of nonlinear signals. M. Bransch [12] calculated elastic-plastic building ground deformations and failure based on FEM, nonlinear HHT-alpha method and scaled boundary FEM. Pines [13] used the empirical mode decomposition and the Hilbert phase information to detect structural damage which is proved to be feasible. Through the scholars’ achievements, signal processing methods based on HHT and improved HHT show that HHT has broad prospects whether in system detection [14] or in structural damage detection [15-18]. However, there still remains great potential of HHT in SHM.
15
16
In parallel with the damage detection in SHM, the data acquisition system and experiment system [19] are also important. N. Roveri [20] performed the structural health monitoring of railway track and train wheels with the help of a FBG(Fiber Bragg grating) sensor array system. Z. Xia [21] investigated a multi-agent system based SHM to achieve the compatibility of different kinds of sensors in monitoring of large scale engineering structures. K. Masayuki [22] evaluated the feasibility of outside strain monitoring method by finite element analyses and experiments. With the advantage of high sensitivity and accuracy, FBG sensor has more and more extensive application in experiment systems and SHM of large-span structures.
17
18
HHT is employed here to get the energy marginal spectrum of signal which can reflect the change of structural mechanical properties. A new damage characteristic parameter is constructed based on the energy marginal spectrum. To verify the proposed damage feature index, it is applied in processing a sinusoidal wave signal with frequency change and identifying pre-defined damage in composite plate vibration experiment. Meanwhile, FBG sensors are used in composites experiment to collect the acceleration data. It is proved that the constructed damage feature index is effective and reliable for detecting structural damage of composites.
19
20
==2.  The Hilbert-Huang Transform==
21
22
The HHT method is introduced here in which the concept of instantaneous frequency is applied. According to this theory, any arbitrary time series are composed of intrinsic mode functions (IMFs), which are different with each other, simple and non-sine. Through the empirical mode decomposition (EMD), an intricate time series can be decomposed into several basic time series from high frequency to low frequency, termed IMFs. After the IMFs of Hilbert transforming, the time-frequency spectrum of original signal is got and the analysis in time domain can proceed. To make sure the above process progress smoothly, each IMF have to satisfy the following conditions: first, the number of extremum and zero-crossings of signal can differ by no more than one; second, the mean value of the upper and lower envelopes defined by the local maxima and the local minima should be zero at any point.
23
24
The EMD method is the first step of HHT to decompose non-stationary signal. For a given time series <math display="inline">x(t)</math>, the process of EMD is introduced as follows. (a)All local extremums are detected firstly, and produce the upper envelope <math display="inline">{x}_{max}</math> and the lower envelope <math display="inline">{x}_{min}</math> by spline approximations using the local maxima and the local minima; (b)obtain <math display="inline">{h}_{1}=</math><math>x\left( t\right) -({x}_{max}+{x}_{min})/2</math>, and restart the process by replacing <math display="inline">{h}_{1}</math> with the original signal if <math display="inline">{h}_{1}</math> does not meet the above IMF criteria; (c) if <math display="inline">{h}_{1}</math> is an IMF, regard <math display="inline">{r=x\left( t\right) -h}_{1}</math> as a new <math display="inline">x(t)</math> and repeat the previous steps till <math display="inline">r</math> becomes a monotonic function or cannot decompose into another IMF any more. Based on the decomposition of the original signal so far, <math display="inline">x(t)</math> can be expressed as Eq.(1).
25
26
{| class="formulaSCP" style="width: 100%; text-align: center;" 
27
|-
28
| 
29
{| style="text-align: center; margin:auto;" 
30
|-
31
| <math display="inline">x\left( t\right) =\sum _{i=1}^{n}{h}_{i}+r</math>
32
|}
33
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
34
|}
35
36
where <math display="inline">\sum _{i=1}^{n}{h}_{i}</math> is the sum of all IMFs, and <math display="inline">r</math> is the residual component. The IMF expresses a simple oscillation style relative to the harmonic function, which is more general.
37
38
After the EMD of <math display="inline">x(t)</math>, Hilbert Transform(HT) is the next step. The HT of <math display="inline">x(t)</math> is defined as the convolution of signal and <math display="inline">\frac{1}{\pi }</math>, which is
39
40
{| class="formulaSCP" style="width: 100%; text-align: center;" 
41
|-
42
| 
43
{| style="text-align: center; margin:auto;" 
44
|-
45
| <math display="inline">y\left( t\right) =\frac{1}{\pi }P\int_{-\infty }^{+\infty }\frac{x(\tau )}{t-\tau }{d}_{\tau }</math>
46
|}
47
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
48
|}
49
50
Then the analytical signal is obtained
51
52
{| class="formulaSCP" style="width: 100%; text-align: center;" 
53
|-
54
| 
55
{| style="text-align: center; margin:auto;" 
56
|-
57
| <math display="inline">z\left( t\right) =x(t)+iy(t)=A(t){e}^{i\theta (t)}</math>
58
|}
59
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
60
|}
61
62
in which <math display="inline">A(t)</math> is the instantaneous amplitude and <math display="inline">\theta (t)</math> is the phase angle as shown in Eq.(4-5).
63
64
{| class="formulaSCP" style="width: 100%; text-align: center;" 
65
|-
66
| 
67
{| style="text-align: center; margin:auto;" 
68
|-
69
| <math display="inline">A\left( t\right) =\sqrt{{x}^{2}\left( t\right) +{y}^{2}(t)}</math>
70
|}
71
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
72
|}
73
74
{| class="formulaSCP" style="width: 100%; text-align: center;" 
75
|-
76
| 
77
{| style="text-align: center; margin:auto;" 
78
|-
79
| <math display="inline">\theta \left( t\right) =\mathrm{arctan}\,\left( \frac{y(t)}{x(t)}\right)</math>
80
|}
81
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
82
|}
83
84
The instantaneous frequency is defined as
85
86
{| class="formulaSCP" style="width: 100%; text-align: center;" 
87
|-
88
| 
89
{| style="text-align: center; margin:auto;" 
90
|-
91
| <math display="inline">\omega \left( t\right) =\frac{d\theta (t)}{dt}</math>
92
|}
93
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
94
|}
95
96
Intuitively, the instantaneous frequency is the differential coefficient of the phase of an analytic signal. The signal is decomposed to the components with local symmetric firstly to get the right instantaneous frequency in HHT.
97
98
Through applying the Hilbert transform to each IMF component, the Hilbert amplitude spectrum is designated as
99
100
{| class="formulaSCP" style="width: 100%; text-align: center;" 
101
|-
102
| 
103
{| style="text-align: center; margin:auto;" 
104
|-
105
| <math display="inline">H\left( \omega ,t\right) =Re\sum _{j=1}^{n}{A}_{j}(t)\mathrm{exp}\,(i\int_{}^{}{\omega }_{j}(t)dt)</math>
106
|}
107
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
108
|}
109
110
which represents the time-frequency distribution of the amplitude of signal. The instantaneous frequency can also be computed according to Eq. (6). Also, the Hilbert marginal spectrum is defined as
111
112
{| class="formulaSCP" style="width: 100%; text-align: center;" 
113
|-
114
| 
115
{| style="text-align: center; margin:auto;" 
116
|-
117
| <math display="inline">h\left( \omega \right) =\int_{0}^{T}H(\omega ,t)dt</math>
118
|}
119
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
120
|}
121
122
and the energy marginal spectrum is
123
124
{| class="formulaSCP" style="width: 100%; text-align: center;" 
125
|-
126
| 
127
{| style="text-align: center; margin:auto;" 
128
|-
129
| <math display="inline">E\left( \omega \right) =\int_{0}^{T}{H}^{2}(\omega ,t)dt</math>
130
|}
131
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
132
|}
133
134
Compared with the Fourier expansion, the HHT marginal spectrum represents the distribution of signal components in different frequencies. The spectrum value of amplitude and the energy of signal in each frequency band are <math display="inline">h\left( {\omega }_{i}\right)</math>  and <math display="inline">E\left( {\omega }_{i}\right)</math> , <math display="inline">i=</math><math>1,n</math>, respectively.
135
136
<span id='OLE_LINK1'></span><span id='OLE_LINK2'></span>
137
138
==3.  Damage detection method==
139
140
Whether the structure is damaged or not can be reflected by the change of response energy in different states. Compared with the amplitude of signal in marginal spectrum, the energy marginal spectrum magnifies the change of response energy generated by structural damage, which is more favorable for the small damage detection.
141
142
{| class="formulaSCP" style="width: 100%; text-align: center;" 
143
|-
144
| 
145
{| style="text-align: center; margin:auto;" 
146
|-
147
| <math display="inline">V=\sum _{i=1}^{n}\left| \frac{{E}^{u}\left( {\omega }_{i}\right) -{E}^{d}\left( {\omega }_{i}\right) }{\mathrm{max}\,[{E}^{u}\left( {\omega }_{i}\right) ]}\right| \times 100\%</math>
148
|}
149
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
150
|}
151
152
In which <math display="inline">{E}^{u}\left( {\omega }_{i}\right)</math>  and <math display="inline">{E}^{d}\left( {\omega }_{i}\right)</math>  are the Hilbert energy of dynamic signal before and after the structure are damaged, respectively. Through Eq. (10), we can see that there is no change before and after the signal energy when the value of <math display="inline">V</math> is 0, which means that structural damage is nonexistent. However, the structural damage occurred when the value does not equal to 0. The bigger value means the greater difference between the signal energy of intact and testing structure, and the damage is more likely to happen. So the structural damage can be indicated by the value of signal energy.
153
154
==4. Numerical simulation==
155
156
In this simulation, a 10s sinusoidal wave with a sudden frequency change at 5s from 1Hz to 1.1Hz is analyzed. Through EMD, the original signal is resolved into 4 IMFs and residue. Original signal and IMFs are shown in Fig. 1. The instantaneous frequency of IMFs and Hilbert spectrum are obtained and given in Figs. 2-3 through HT. As shown in Fig. 2, the frequency shift at 5s can be directly captured from instantaneous frequency of the first IMF while the shifts of other IMFs are not obvious. The change of amplitude in Hilbert spectrum also reflects the frequency shift of signal at 5s.
157
158
After the Hilbert energy is obtained, the value of signal energy can be calculated through Eq.10 by comparing two signals in different conditions. In the first case, the original signal is divided into two parts, which is the signal of front 5s and last 5s, respectively. The value of <math display="inline">V</math> is 9.23%, which means that the signal parameter of last 5s has changed compared with the signal of front 5s.In the second case, the 10s signal of 1Hz is compared with the original signal. In this condition the value of <math display="inline">V</math> is 3.7%, which also means that the character of two signals is different and structural damage occurs.
159
160
161
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
162
|-
163
| [[Image:draft_Hua_323164655-image2-c.jpeg|510px]]
164
|- style="text-align: center; font-size: 75%;"
165
| colspan="1" | '''Figure 1.''' Original sine wave and IMFs.
166
|}
167
168
169
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
170
|-
171
| [[Image:draft_Hua_323164655-image3-c.png|510px]]
172
|- style="text-align: center; font-size: 75%;"
173
| colspan="1" | '''Figure 2.''' Instantaneous frequency of the IMFs.
174
|}
175
176
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
177
|-
178
| [[Image:draft_Hua_323164655-image4-c.jpeg|498px]]
179
|- style="text-align: center; font-size: 75%;"
180
| colspan="1" | '''Figure 3.''' Hilbert spectrum.
181
|}
182
183
==5.  Experimental verification==
184
185
In Section 4, the parameter based on HHT was used to identify the change of a sinusoidal wave which was proved to be effective. To test the structural damage detection method, a vibration experiment of composite plate form wing skin was designed and performed.
186
187
<span id='_GoBack'></span>Composite material used in test was T300/QY8911 and the material property is shown in Table 1. The thickness of single layer was 0.18mm and stacking sequence of composite laminates was <math display="inline">{[45/-45/0/-45/45/0]}_{S}</math>. One side of the <math display="inline">200\times 2000 \times  2.16</math> mm<sup>3</sup> composite plate was clamped to vibrating table. There were two symmetrical isosceles triangular breaches near the clamped end. The acceleration of vibrating table was 4g and the frequency of vibration was 74Hz.
188
189
190
<div class="center" style="font-size: 75%;">'''Table 1'''. Material properties of T300/QY8911 (MPa)
191
</div>
192
193
{| style="width: 100%;border-collapse: collapse;font-size:85%;text-align: center;" 
194
|-
195
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{E}_{XX}</math>
196
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{E}_{YY}</math>
197
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{E}_{ZZ}</math>
198
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{G}_{XY}</math>
199
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{G}_{XZ}</math>
200
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{G}_{YZ}</math>
201
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{v}_{XY}</math>
202
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{v}_{XZ}</math>
203
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"| <math>{v}_{YZ}</math>
204
|-
205
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|135000
206
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|8800
207
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|8800
208
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|4470
209
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|4470
210
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|3200
211
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|0.33
212
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;border-right: 1pt solid black;"|0.33
213
|  style="border-left: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;text-align: center;border-right: 1pt solid black;"|0.48
214
|}
215
216
217
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
218
|-
219
| style="padding-top:10px;" |[[Image:draft_Hua_323164655-image5-c.jpeg|center|600px]]
220
|- style="text-align: center; font-size: 75%;"
221
| colspan="1" style="padding-bottom:10px;"| '''Figure 4.''' Composite plate with pre-defined damage.
222
|}
223
224
225
As shown in Fig. 4, a crack in one breach was made to be pre-defined damage. Structural damage was detected and located by comparing the signal of damaged structure with the intact one. The sample frequency was 1000Hz. Acceleration sensors were distributed as shown in Fig. 5. The wavelengths of FBG sensors are 1541nm, 1555nm, 1538nm and 1558nm, respectively.
226
227
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width:auto;max-width: auto;"
228
|-
229
| style="padding-top:10px;" |[[Image:draft_Hua_323164655-image6.png|500px]]
230
|- style="text-align: center; font-size: 75%;"
231
| colspan="1" style="padding-bottom:10px;"| '''Figure 5.''' Sensor set and damage location.
232
|}
233
234
235
Taking the acceleration data of sensor 4 for example, HHT is employed in signal processing of intact structure and damaged one. The results of 0.5s signal in two cases are plotted in Figs. 6-9. As shown in Fig. 6, two original signals were decomposed into 4 IMFs respectively. From the instantaneous frequency given in Fig. 7, it can be found that there is a little end-effect around boundary area in IMF of low frequency which is ignorable since the energy of these frequencies is quite low. It is hard to find the difference between the two HHT spectrums from Fig. 8. The marginal energy of intact laminate ranges from 40Hz to 100Hz, while the damaged one ranges from 50Hz to 100Hz are shown in Fig. 9.
236
237
238
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
239
|-
240
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image7-c.jpeg|404px]]
241
|-
242
| style="text-align: center; font-size:75%;| (a) IMFs of intact laminate
243
|-
244
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image8-c.jpeg|400px]]
245
|-
246
| style="text-align: center; font-size:75%;| (b) IMFs of damaged laminate
247
|- style="text-align: center; font-size: 75%;"
248
| colspan="1" style="padding-bottom:10px;padding-top:10px;"| '''Figure 6.''' EMD results of structural signals.
249
|}
250
251
252
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
253
|-
254
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image9-c.gif|498px]]
255
|-
256
| style="text-align: center; font-size:75%;| (a) Instantaneous frequency of intact laminate
257
|-
258
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image10-c.gif|498px]] 
259
|-
260
| style="text-align: center; font-size:75%;| (b) Instantaneous frequency of damaged laminate
261
|- style="text-align: center; font-size: 75%;"
262
| colspan="1" style="padding-bottom:10px;padding-top:10px;"| '''Figure 7.''' Instantaneous frequency of structural signals.
263
|}
264
265
266
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
267
|-
268
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image11-c.jpeg|498px]]
269
|-
270
| style="text-align: center; font-size:75%;| (a) HHT spectrum of intact laminate
271
|-
272
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image12-c.jpeg|498px]]
273
|-
274
| style="text-align: center; font-size:75%;| (b) HHT spectrum of damaged laminate
275
|- style="text-align: center; font-size: 75%;"
276
| colspan="1" style="padding-bottom:10px;padding-top:10px;"| '''Figure 8.''' HHT spectrum of structural signals.
277
|}
278
279
280
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;max-width: 100%;"
281
|-
282
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image13-c.gif|498px]]
283
|-
284
| style="text-align: center; font-size:75%;| (a) Energy marginal spectrum of intact laminate
285
|-
286
| style="padding-top:10px;padding-bottom:10px;" |[[Image:draft_Hua_323164655-image14-c.gif|498px]]
287
|-
288
| style="text-align: center; font-size:75%;| (b) Energy marginal spectrum of damaged laminate
289
|- style="text-align: center; font-size: 75%;"
290
| colspan="1" style="padding-bottom:10px;padding-top:10px;"| '''Figure 9.''' Energy marginal spectrum of structural signals.
291
|}
292
293
294
Based on Eq. (10), the change of signal energy can be obtained as 6.63%. Since <math display="inline">V</math> of sensor 4 signal is greater than 0, structural damage has occurred. However, the location of damage has not been captured yet. <math display="inline">V</math> of structural signals are shown in Table 2. From comparing the change of signal energy, <math display="inline">V</math> of sensor 3 signal is bigger which means that the structural damage is more likely to happen around this sensor. The damage detection results are in consistent with the pre-defined damage condition.
295
296
<div class="center" style="font-size: 75%;">'''Table 2'''. V of structural response signals
297
</div>
298
299
{| style="width: auto;border-collapse: collapse;font-size:85%;" 
300
|-
301
|  style="border-left: 1pt solid black;border-top: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
302
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Sensor 1
303
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Sensor 2
304
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Sensor 3
305
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Sensor 4
306
|-
307
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|V
308
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|3.02%
309
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|1.47%
310
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|13.7%
311
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|6.63%
312
|}
313
314
==6.  Conclusions==
315
316
The Hilbert-Huang transformation is an effective method for analyzing nonlinear and non-stationary signal. A new structural damage detection parameter based on the Hilbert-Huang transformation is developed here. The damage characteristic value obtained from Hilbert energy spectrum is used to identify structural damage in a sinusoidal wave with a sudden frequency change and a composite plate with pre-defined damage. The change of signal energy is got through comparing the signal energy of intact and damaged structural which can implement the detection and localization of structural damage. Through the results of simulation and experiment, the following conclusions can be drawn: 1) the Hilbert spectrum energy can reflect the distribution of signal energy and the constructed structural damage feature index can be used in identifying structural damage and localizing the damage; 2) signal energy variation is sensitive to identify small damage without specific structural parameter which satisfies the requirement of structural health monitoring of dynamic and real-time system; 3) the size of FBG acceleration sensor used in composites vibration test is quite relative to test sample, which can be substituted by sensor with smaller volume in future research such as strain sensor. After the detection and localization of structural damage, the assessment of damage level is the future goal. Meanwhile, this parameter needs to be further improved which includes eliminating the end-effect in EMD.
317
318
==Acknowledgements==
319
320
This study is supported by the National Natural Science Foundation (No: 51475369).
321
322
==References==
323
324
[1] X.J. Dang, Statistic strategy of damage detection for composite structure using the correlation function amplitude vector. Procedia Engineering, 99(2015)1395–1406.
325
326
[2] M.M. Alamdari, J. Li, B. Samali, Damage detection using 2-D discrete wavelet transform on extended operational mode shapes. Archives of Civil and Mechanical Engineering, 15(2015)698-710.
327
328
[3] L. Montanari, B. Basu, et al., A padding method to reduce edge effects for enhanced damage detection using wavelet analysis. Mechanical Systems and Signal Processing, 52(2015)264-277.
329
330
[4] Y. Song, C.R. Bowen, H. Alicia Kim, et al., Non-invasive damage detection in beams using marker and wavelets. Mechanical Systems and Signal Processing, 49(2014)13-23.
331
332
[5] N.E. Huang, Shen, et al., The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings of Royal Society London, 454(1998)903-995.
333
334
[6] L. Law, J. Kim, et al., An approach based on wavelet packet decomposition and Hilbert-Huang transform (WPD-HHT) for spindle bearing condition monitoring. Mechanical Systems and Signal Processing, 33(2012)197-211.
335
336
[7] W. L. Bayissa, N. Haritos, S. Thelandersson, Vibration-based structural damage detection using wavelet transform. Mechanical Systems and Signal Processing, 22(2008)1194–1215.
337
338
[8] D.J. Huang, J.P. Zhao, et al., Practical implementation of Hilbert-Huang Transform algorithm. Acta Oceanologiea Sinica, 22(2003)1-14.
339
340
[9] Z. Wang, G. Chen, Analytical mode decomposition with Hilbert transform for modal parameter detection of building under ambient vibration. Engineering Structures, 59(2014)173-184.
341
342
[10] R. Jha, S. Xu, G. Ahmadi. Health monitoring of a multi-level structure based on empirical mode decomposition and hilbert spectral analysis. In: Chang FK, editor. Fifth international workshop on structural health monitoring, CA, USA: Stanford University; 2005.
343
344
[11] C.X. Bao, H. Hong, et al., Time-varying system detection using a newly improved HHT algorithm. Computers and Structures, 87(2009)1611–1623.
345
346
[12] M. Bransch, L. Lehmann, A nonlinear HHT-a method with elastic–plastic soil–structure interaction in a coupled SBFEM/FEM approach. Computers and Geotechnics, 38(2011)80–87.
347
348
[13] D. Pines, L. Salvino. Structural health monitoring using empirical mode decomposition and the Hilbert phase. Journal of Sound and Vibration, 294(2006)97–124.
349
350
[14] T.L. Huang, M.L. Lou. System detection of nonlinear structures based on HHT. Earthq Eng Eng Vib, 26(2006)80–3.
351
352
[15] Z.K. Peng, W.T. Peter, F.L. Chu. An improved Hilbert–Huang transform and its application in vibration signal analysis. Journal of Sound and Vibration, 286(2004)187–205.
353
354
[16] P. F. Pai, A.N. Palazotto. HHT-based nonlinear signal processing method for parametric and non-parametric detection of dynamical systems. International Journal of Mechanical Sciences, 50(2008)1619-1635.
355
356
[17] S.T. Quek, P.S. Tua, Q. Wang. Detecting anomalies in beams and plate based on the Hilbert-Huang transform of real signals. Smart Mater. Struct., 12 (2003)447–460.
357
358
[18] Q.Zhang, Y.Wang, et al., Hilbert-Huang Transform based method for monitoring the crack of concerete arch by using FBG sensors. Optik, 127(2016)3417-3422.
359
360
[19] C. Heng, M. Kurt, S.L. Young, et al., Experimental system detection of the dynamics of a vibro-impact beam with a view towards structural health monitoring and damage detection. Mechanical Systems and Signal Processing, 46(2014)91–113.
361
362
[20] N. Roveri, A. Carcaterra, A. Sestieri, Real-time monitoring of railway infrastructures using fibre Bragg grating sensors. Mechanical Systems and Signal Processing, 60-61(2015)14–28.
363
364
[21] Z. Xia, Shenfang Yuan, Zhenhua Yu, et al., Designing strategy for multi-agent system based large structural health monitoring. Expert Systems with Applications, 34(2008)1154–1168.
365
366
[22] K. Masayuki, M. Koji, Monitoring of inside surface crack growth by strain measurement of the outside surface: A feasibility study. Nuclear Engineering and Design, 241(2011)1–11.
367

Return to Hua et al 2017a.

Back to Top

Document information

Published on 03/01/18
Accepted on 26/10/17
Submitted on 12/07/17

Volume 34, Issue 1, 2018
DOI: 10.23967/j.rimni.2017.9.002
Licence: CC BY-NC-SA license

Document Score

0

Views 289
Recommendations 0

Share this document