# Output-Only Modal Analysis Based on Improved Empirical Mode Decomposition Method

Authors: Shiqiang Qin; Qiuping Wang; Juntao Kang
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

The output-only modal analysis for bridge structures based on improved empirical mode decomposition (EMD) is investigated in this study. First, a bandwidth restricted EMD is proposed for decomposing nonstationary output measurements with close frequency components. The advantage of bandwidth restricted EMD to standard EMD is illustrated by a numerical simulation. Next, the modal parameters are extracted from intrinsic mode function obtained from the improved EMD by both random decrement technique and stochastic subspace identification. Finally, output-only modal analysis of a railway bridge is presented. The study demonstrates the mode mixing issues of standard EMD can be restrained by introducing bandwidth restricted signal. Further, with the improved EMD method, band-pass filter is no longer needed for separating the closely spaced frequency components. The modal parameters extracted based on the improved EMD method show good agreement with those extracted by conventional modal identification algorithms.


---

## Body

## 1. Introduction

Monitoring, maintenance, and management at operation stage play an important role in the whole life span of large-scale bridge structures. Damage detection in time and proper repairing measures prolong the service life of bridge. Thus, the structural health monitoring (SHM) of bridge structures has gained more and more concern by both researchers and practicing engineers. By installing various sensors on bridge, mechanical behaviors and dynamic responses of bridge are recorded under normal load conditions (ambient, vehicles, etc.) and extreme load conditions (hurricane, earthquake, etc.). In this way, structural health condition can be assessed resulting in good advices for maintenance and management. Besides, the recorded data can be further used to validate and improve the design theory of bridge with nonconventional structure system. Successful applications of SHM for bridge structures can be found in [1–6] and the references therein.A relevant research topic in SHM is to experimentally determine the dynamic characteristics with high accuracy. For bridge structures, the basic dynamic properties are the modal parameters, that is, natural frequencies, damping ratios, and mode shapes. These modal properties are the inherent properties of structure, only related with mass distribution, stiffness, and boundary condition. Variation of modal parameter implies change in the inherent properties of structure, by which structural damage can be detected. Therefore, modal parameters are often served as indicators in modal-based damage detection techniques [7, 8]. In addition, experimentally determined modal parameters are regarded as key parameters for updating and calibrating the preliminary finite element (FE) model developed from design documents [9–12]. Served as baseline for SHM, the updated FE model can be used to predict the dynamic responses and damage detection of bridge. Furthermore, the modal parameters are important for condition assessment; see [13], for example. Traditionally, modal analysis is often performed based on both input and output measurements of structures, meaning special excitation devices, such as drop weights, shakers, are needed. However, these requirements are difficult to be fulfilled for large-scale bridge structures. This is because for one thing, bridge structure cannot be fully excited by only one or several excitation devices; moreover, these devices are very expensive for large-scale application. For the other, the operation service needs to be interrupted. Therefore, it is more meaningful to extract modal parameters from operational vibration measurements, namely, output-only modal analysis or operational modal analysis (OMA). At present, there are several sophisticated approaches for output-only modal analysis of bridge structures. These OMA approaches can be grouped into two categories, namely, frequency domain methods and time domain methods. The typical frequency domain methods include Pick-Peaking (PP) from the power spectra density [14], frequency domain decomposition (FFD) [15], and the operational PolyMax [16]. The common used time domain methods contain auto regressive moving averaging (ARMA) method [17], the natural excitation method (NExT) [18], the least square complex exponential (LSCE) method [19], the Ibrahim time domain (ITD) method [20], and the stochastic subspace identification (SSI) [21].Two commonly used assumptions of output-only modal analysis algorithms are the white noise assumption of unmeasured input excitations as well as the stationary process assumption of the outputs. For example, the ambient excitations are assumed to be white noise in SSI algorithm so that they can be combined with test inaccuracies and modelling errors. Thus, the stochastic state space model is obtained from the discrete state space model. However, operational vibration measurements are not always stationary signals. Although the existing output-only modal analysis algorithms yield satisfactory results, it is still worthy investigating alternative approaches for modal identification with nonstationary output-only responses, especially from theoretical perspective.Several examples of output-only modal analysis performed in large civil structures have been reported in technical literature, of which the most relevant cases regarding nonstationary operational vibration measurements are presented herein. Yang et al. [22, 23] proposed a modal identification method based on Hilbert-Huang transform (HHT) in conjunction with band-pass filter using measured free vibration time histories. Yang’s study also pointed out that limited by poor frequency resolution of wavelet transform (WT), the identified results from HHT seem to be more accurate than those from WT. He et al. [24] extracted the modal parameters of Nanjing Railway Bridge based on empirical mode decomposition (EMD) and random decrement technique (RDT). Yu and Ren [25] proposed an EMD based SSI modal identification method, which reduced the postprocessing work in SSI algorithm (picking the stable poles). Band-pass filter was also employed to treat the closely spaced modes. Yan and Miyamoto [26] made a comparative study between HHT and WT for modal analysis. The results found that for well-separated modes both methods are applicable, whereas for close spaced modes, the WT method is more effective. Abovementioned literatures proved that it is possible to extract modal parameters from nonstationary operational vibration measurements, but more attention should be paid to closely spaced modes. The bridge structures, often large in scale, fall into the category has low frequencies and close modes. Although some useful algorithms such as ensemble EMD (EEMD) [27] and the masking signal method [28] have been proposed to deal with close spaced modes, it is still worthy to investigate this problem, especially that EMD is often combined with time domain algorithms to identify modal parameters.This paper focuses on output-only modal analysis of bridge structure based on an improved EMD method. The bandwidth restricted EMD is proposed to deal with closely spaced modes. Without band-pass filter, the closely spaced modes can be well separated with bandwidth restricted EMD. Both RDT and SSI algorithm are employed to identify modal parameters when the single frequency component of original signal is obtained. The text of this paper is organized as follows: in Section2, the improved EMD method is introduced. A numerical simulation is used to illustrate the effectiveness of bandwidth restricted EMD in decomposing signals with close frequency components. Section 3 explains the basic principles of extracting modal parameters from intrinsic mode functions by RDT and SSI. The output-only modal analysis of a railway bridge based on improved EMD is accomplished in Section 4. The PP and SSI algorithms are used for comparison. Conclusions are made at last.


## 2. An Improved EMD Method

### 2.1. Standard EMD

The EMD, first proposed by Huang et al. [29], is an effective approach to decompose nonstationary, multicomponent signal. It can decompose multicomponent signal into a series of intrinsic mode functions (IMFs) and a residual. By definition an IMF has to satisfy the following two conditions:(
1
) In the whole data set, the number of extremes and the number of zero crossings may differ but no more than one; (
2
) At any point, the local average is zero.For an arbitrary time signal, the basic steps of EMD are as follows.(a)
Sifting process: Find all the extremes ofx
t. Define the upper and lower envelope of x
t, which can be obtained through connecting all the local maxima and local minima with cubic spline curve, respectively. Calculate the mean value of the two envelopes. The difference between the x
t and the mean value is designated as h
1
t.
(b)
Checking process: Check whetherh
1
t satisfies the two conditions of IMF. If h
1
t is an IMF, then take c
1
t
=
h
1
t as the first IMF of x
t. Otherwise, let x
t
=
h
1
t and repeat the sifting process until the first IMF is found out.
(c)
Looping process: Definer
1
t
=
x
t
-
c
1
t as the new signal and repeat (a), (b) until r
n
t is smaller than a predetermined threshold. Then, the original signal x
t can be rewritten as(1)
x
t
=
∑
i
=
1
n
c
i
t
+
r
n
t
,where c
i
(
t
)
(
i
=
1,2
,
3
,
…
,
n
) represents the IMFs of the x
(
t
) from high frequency to low frequency component. r
n
(
t
) is a single value function and represents the residue of x
(
t
). The flow chart of standard EMD is shown in Figure 1. Ideally, each IMF should only contain one frequency component, but in practice, it is often the case that one IMF contains several frequency components. This phenomenon is called mode mixing. It will be further discussed in the next section.Figure 1
The flow chart of standard EMD.


### 2.2. Bandwidth Restricted EMD

As discussed in Section2.1, EMD is an effective signal processing method especially when dealing with nonlinear and nonstationary signals. But it still faces some big challenges. Mode mixing is one of these challenges. It means that one IMF contains multifrequency components or one frequency component is decomposed into several IMFs. Under the following two conditions, it may lead to mode mixing during EMD:(1)
The intermittency of signal: This means that the signal has intermittency frequency components. In this case, several different frequency components will be decomposed into one IMF. Figure2 gives an example of mode mixing because of intermittency of signal. x
1
(
t
) and x
2
(
t
) are two frequency components of signal x
(
t
), and imf1 and imf2 are the IMFs obtained by applying standard EMD to x
(
t
). It is obvious that because of intermittency of x
1
(
t
), both imf1 and imf2 contain multifrequency components.
(2)
The closely spaced frequency components of signal: In this case, one frequency component will be decomposed into several IMFs. Figure3 shows an example of this kind of mode mixing. x
1
(
t
) and x
2
(
t
) are two sinusoidal waves with the frequency of 1.8 Hz and 2.5 Hz, respectively. imf1 and imf2 are two IMFs obtained from standard EMD. It is obvious that the two frequency components are decomposed into both IMFs.Figure 2
Mode mixing due to intermittency of signal.x
1
(
t
) and x
2
(
t
) are the two frequency components of x
(
t
). imf1 and imf2 are the two IMFs obtained from standard EMD of x
(
t
).Figure 3
Mode mixing due to closely spaced frequency components of signal.x
1
(
t
) and x
2
(
t
) are the two frequency components of x
(
t
). imf1 and imf2 are the two IMFs obtained from standard EMD of x
(
t
).In this paper, a bandwidth restricted EMD is proposed to deal with the problem of mode mixing. It has the advantage of easy implementation and does not need the a priori knowledge of the signal to be analyzed. From the theory background of EMD, it is known that each IMF is obtained by subtracting the mean value of upper bound and lower bound of original signal, meaning that EMD has a property similar to binary filter. The proposed method is based on this property and can be realized as the following steps:(1)
Supposex
(
t
) is a time history signal. Apply standard EMD to x
(
t
), set stop criterion options to make it decompose x
(
t
) into onlyone IMF and a residual; that is, (2)
x
t
=
I
M
F
1
t
+
r
1
t
.
(2)
Calculate the instant amplitudea
1
(
i
) and instant frequency f
1
(
i
) of I
M
F
1
(
t
) based on Hilbert transform. The frequency content of bandwidth restricted signal is set as(3)
f
=
A
0
×
∑
i
=
1
k
a
1
i
f
1
2
i
∑
i
=
1
k
a
1
i
f
1
i
,where A
0 is bandwidth restricted signal coefficient.
(3)
Construct the bandwidth restricted signals
(
t
) by(4)
s
t
=
B
sin
⁡
2
π
f
t
,where B is a constant; it can be selected according to the amplitude of x
(
t
).
(4)
Apply standard EMD toy
(
t
)
=
x
(
t
)
+
s
(
t
), making that it only decomposes one IMF and a residual(5)
y
t
=
I
M
F
1,1
t
+
r
1,1
t
.
(5)
Check whetherh
1
(
t
)
=
I
M
F
1,1
(
t
)
-
s
(
t
) leads to mode mixing. If not, take h
1
(
t
) as the first IMF of x
(
t
). Otherwise, take h
1
(
t
) as x
(
t
) and repeat (1)~(4) for n times until h
1
,
n
(
t
) do not lead to mode mixing any more. In order to check whether an IMF leads to mode mixing, the frequency distribution of the IMF obtained from Hilbert transform can be employed.
(6)
Repeat (1)~(5) until all IMFs of x
(
t
) are extracted.Figure4 shows the flow chart of bandwidth restricted EMD.Figure 4
The flow chart of bandwidth restricted EMD.Numerical Simulation. To evaluate the effect of bandwidth restricted EMD, a simulation example is presented here. The simulation signal y
(
t
) is composed of three sinusoidal waves with frequency of 10 Hz, 15 Hz, and 18 Hz. It can be expressed as y
t
=
sin
⁡
(
20
π
t
)
+
sin
⁡
(
30
π
t
)
+
sin
⁡
(
36
π
t
). The sampling frequency is 200 Hz and the sampling time is 2.5 s. Figure 5 shows the three IMFs obtained from EMD (C
1~C
3) and from bandwidth restricted EMD (BR-C
1~BR-C
3). It is clearly shown in the figure that C
1~C
3  lead to mode mixing while BR-C
1~BR-C
3 contain only one frequency component. The simulation example shows that the bandwidth restricted EMD effectively solved the problem of mode mixing.Figure 5
The comparison of IMFs obtained from EMD (C
1~C
3) and bandwidth restricted EMD (BR-C
1~BR-C
3). y
(
t
) is the simulated original signal.


## 2.1. Standard EMD

The EMD, first proposed by Huang et al. [29], is an effective approach to decompose nonstationary, multicomponent signal. It can decompose multicomponent signal into a series of intrinsic mode functions (IMFs) and a residual. By definition an IMF has to satisfy the following two conditions:(
1
) In the whole data set, the number of extremes and the number of zero crossings may differ but no more than one; (
2
) At any point, the local average is zero.For an arbitrary time signal, the basic steps of EMD are as follows.(a)
Sifting process: Find all the extremes ofx
t. Define the upper and lower envelope of x
t, which can be obtained through connecting all the local maxima and local minima with cubic spline curve, respectively. Calculate the mean value of the two envelopes. The difference between the x
t and the mean value is designated as h
1
t.
(b)
Checking process: Check whetherh
1
t satisfies the two conditions of IMF. If h
1
t is an IMF, then take c
1
t
=
h
1
t as the first IMF of x
t. Otherwise, let x
t
=
h
1
t and repeat the sifting process until the first IMF is found out.
(c)
Looping process: Definer
1
t
=
x
t
-
c
1
t as the new signal and repeat (a), (b) until r
n
t is smaller than a predetermined threshold. Then, the original signal x
t can be rewritten as(1)
x
t
=
∑
i
=
1
n
c
i
t
+
r
n
t
,where c
i
(
t
)
(
i
=
1,2
,
3
,
…
,
n
) represents the IMFs of the x
(
t
) from high frequency to low frequency component. r
n
(
t
) is a single value function and represents the residue of x
(
t
). The flow chart of standard EMD is shown in Figure 1. Ideally, each IMF should only contain one frequency component, but in practice, it is often the case that one IMF contains several frequency components. This phenomenon is called mode mixing. It will be further discussed in the next section.Figure 1
The flow chart of standard EMD.


## 2.2. Bandwidth Restricted EMD

As discussed in Section2.1, EMD is an effective signal processing method especially when dealing with nonlinear and nonstationary signals. But it still faces some big challenges. Mode mixing is one of these challenges. It means that one IMF contains multifrequency components or one frequency component is decomposed into several IMFs. Under the following two conditions, it may lead to mode mixing during EMD:(1)
The intermittency of signal: This means that the signal has intermittency frequency components. In this case, several different frequency components will be decomposed into one IMF. Figure2 gives an example of mode mixing because of intermittency of signal. x
1
(
t
) and x
2
(
t
) are two frequency components of signal x
(
t
), and imf1 and imf2 are the IMFs obtained by applying standard EMD to x
(
t
). It is obvious that because of intermittency of x
1
(
t
), both imf1 and imf2 contain multifrequency components.
(2)
The closely spaced frequency components of signal: In this case, one frequency component will be decomposed into several IMFs. Figure3 shows an example of this kind of mode mixing. x
1
(
t
) and x
2
(
t
) are two sinusoidal waves with the frequency of 1.8 Hz and 2.5 Hz, respectively. imf1 and imf2 are two IMFs obtained from standard EMD. It is obvious that the two frequency components are decomposed into both IMFs.Figure 2
Mode mixing due to intermittency of signal.x
1
(
t
) and x
2
(
t
) are the two frequency components of x
(
t
). imf1 and imf2 are the two IMFs obtained from standard EMD of x
(
t
).Figure 3
Mode mixing due to closely spaced frequency components of signal.x
1
(
t
) and x
2
(
t
) are the two frequency components of x
(
t
). imf1 and imf2 are the two IMFs obtained from standard EMD of x
(
t
).In this paper, a bandwidth restricted EMD is proposed to deal with the problem of mode mixing. It has the advantage of easy implementation and does not need the a priori knowledge of the signal to be analyzed. From the theory background of EMD, it is known that each IMF is obtained by subtracting the mean value of upper bound and lower bound of original signal, meaning that EMD has a property similar to binary filter. The proposed method is based on this property and can be realized as the following steps:(1)
Supposex
(
t
) is a time history signal. Apply standard EMD to x
(
t
), set stop criterion options to make it decompose x
(
t
) into onlyone IMF and a residual; that is, (2)
x
t
=
I
M
F
1
t
+
r
1
t
.
(2)
Calculate the instant amplitudea
1
(
i
) and instant frequency f
1
(
i
) of I
M
F
1
(
t
) based on Hilbert transform. The frequency content of bandwidth restricted signal is set as(3)
f
=
A
0
×
∑
i
=
1
k
a
1
i
f
1
2
i
∑
i
=
1
k
a
1
i
f
1
i
,where A
0 is bandwidth restricted signal coefficient.
(3)
Construct the bandwidth restricted signals
(
t
) by(4)
s
t
=
B
sin
⁡
2
π
f
t
,where B is a constant; it can be selected according to the amplitude of x
(
t
).
(4)
Apply standard EMD toy
(
t
)
=
x
(
t
)
+
s
(
t
), making that it only decomposes one IMF and a residual(5)
y
t
=
I
M
F
1,1
t
+
r
1,1
t
.
(5)
Check whetherh
1
(
t
)
=
I
M
F
1,1
(
t
)
-
s
(
t
) leads to mode mixing. If not, take h
1
(
t
) as the first IMF of x
(
t
). Otherwise, take h
1
(
t
) as x
(
t
) and repeat (1)~(4) for n times until h
1
,
n
(
t
) do not lead to mode mixing any more. In order to check whether an IMF leads to mode mixing, the frequency distribution of the IMF obtained from Hilbert transform can be employed.
(6)
Repeat (1)~(5) until all IMFs of x
(
t
) are extracted.Figure4 shows the flow chart of bandwidth restricted EMD.Figure 4
The flow chart of bandwidth restricted EMD.Numerical Simulation. To evaluate the effect of bandwidth restricted EMD, a simulation example is presented here. The simulation signal y
(
t
) is composed of three sinusoidal waves with frequency of 10 Hz, 15 Hz, and 18 Hz. It can be expressed as y
t
=
sin
⁡
(
20
π
t
)
+
sin
⁡
(
30
π
t
)
+
sin
⁡
(
36
π
t
). The sampling frequency is 200 Hz and the sampling time is 2.5 s. Figure 5 shows the three IMFs obtained from EMD (C
1~C
3) and from bandwidth restricted EMD (BR-C
1~BR-C
3). It is clearly shown in the figure that C
1~C
3  lead to mode mixing while BR-C
1~BR-C
3 contain only one frequency component. The simulation example shows that the bandwidth restricted EMD effectively solved the problem of mode mixing.Figure 5
The comparison of IMFs obtained from EMD (C
1~C
3) and bandwidth restricted EMD (BR-C
1~BR-C
3). y
(
t
) is the simulated original signal.


## 3. Modal Parameter Identification Using IMFs

By introducing bandwidth restricted EMD, the mode mixing of standard EMD is effectively restrained. The obtained IMF is single frequency component, which can be seen as the dynamic response of a single degree of freedom (SDOF) system. Thus, the remaining problem is to identify modal parameters from SDOF system. Both RDT and SSI can be used to extract modal parameters from IMFs.


### 3.1. Random Decrement Technique

RDT is a time domain signal processing approach that is developed for extracting free decay response from stationary linear structure vibration response. Like NExT and LSCE method mentioned in Section1, RDT is also a time domain technique. These modal identification techniques, including RDT, NExT, and LSCE, are also capable of dealing with closely spaced modes [18, 19]. However, EMD based RDT modal identification technique is more suitable for large-scale structures since both excitations and responses are nonstationary processes.The basic principles of RDT are described as follows. For a SDOF linear vibration system, the displacement responsex
(
t
) is composed of three parts. The first part results from initial displacement, the second is from initial velocity, and the third is from random excitations. Thus, x
(
t
) can be expressed as(6)
x
t
=
x
0
D
t
+
x
˙
0
V
t
+
∫
0
t
h
t
-
τ
f
τ
d
τ
,where x
(
0
) and x
˙
(
0
) are the initial displacement and initial velocity, respectively. D
(
t
) is the free vibration response of the system with initial displacement being 1 and initial velocity being 0. V
(
t
) is the free vibration response with initial displacement being 0 and initial velocity being 1. f
(
τ
) is the random excitations at time segment τ. h
(
t
) is unit impulse response function of SDOF system. If the stochastic part can be separated and eliminated from x
(
t
), one can expect the free vibration part of x
(
t
). Select a reasonable value A to intercept the displacement response x
(
t
), the crossing points are denoted as (
t
i
,
x
(
t
i
)
)
(
i
=
1,2
,
3
,
…
). A series of subsample functions x
(
t
-
t
i
) can be obtained through moving the start time to the crossing points,(7)
x
t
-
t
i
=
x
t
i
D
t
-
t
i
+
x
˙
t
i
V
t
-
t
i
+
∫
t
i
t
h
t
-
τ
f
τ
d
τ
.In output-only modal analysis, the random excitations are assumed to be stationary Gaussian white noise. Therefore, the change of start time will not change the stochastic characteristic of random excitation functions. Thus, a series of random process can be obtained from(8)
x
i
t
=
A
D
t
+
x
˙
t
i
V
t
+
∫
0
t
h
t
-
τ
f
τ
d
τ
.Sincex
˙
t
i is a zero-mean stationary Gaussian process, hence, E
x
˙
t
i
=
0 and E
f
t
=
0. The mathematical expectation E
x
i
t is obtained by(9)
E
x
i
t
=
A
D
t
+
E
x
˙
t
i
V
t
+
∫
0
t
h
t
-
τ
E
f
τ
=
A
D
t
,which represents the free vibration response with initial displacement being A and initial velocity being 0.After obtaining the free decays of IMF, the modal parameter can be estimated by Hilbert transform. Without loss of generality, the free vibration response in (9) can be rewritten as(10)
x
t
=
A
0
e
-
ξ
ω
0
t
cos
⁡
ω
d
t
+
ϕ
0
,where ω
0 is the circular frequency and ω
d is the damping frequency. ξ is the damping ratio and A
0 is a constant. The corresponding analytical expression of x
(
t
) is obtained by(11)
z
t
=
x
t
+
i
x
~
t
=
A
t
e
-
i
θ
t
,where x
~
(
t
) is the Hilbert transform of x
(
t
). The amplitude function A
(
t
) and phase angle function θ
t can be expressed as(12)
A
t
=
A
0
e
-
ξ
ω
0
t
θ
t
=
ω
d
t
+
ϕ
0
.By introducing logarithm and differential operator, (12) yields(13)
ln
⁡
A
t
=
ln
⁡
A
0
-
ξ
ω
0
t
ω
t
=
d
θ
t
d
t
=
ω
d
.Equation (13) shows that the damping frequency is calculated from the slope of phase angle function. The slope of amplitude function ξ
ω
0 can also be obtained. Thus, the frequency ω
0 and damping ratio ξ are calculated from ω
d
=
ω
0
1
-
ξ
2.


### 3.2. Stochastic Subspace Identification

SSI is a robust time domain modal identification approach. Many successful applications based on SSI have been reported. The stochastic state space model is described as(14)
x
k
+
1
=
A
x
k
+
w
k
y
k
=
C
x
k
+
v
k
,where x
k is the state vector; y
k is the output vector; w
k and v
k denote the modelling inaccuracies and measuring errors, respectively; A is the state matrix and C output matrix. The main objective of SSI is to identify A and C from output-only response. The modal parameters can be obtained by applying eigenvalue decomposition to A. Based on the state space model, there are mainly two implementations of SSI algorithm: the first one directly starts from the sampled vibration data, namely, data driven SSI (SSI-data) and the other one from the covariance of the vibration data which is called covariance driven SSI (SSI-cov). Both SSI-data and SSI-cov reveal similar results with satisfying accuracy, while SSI-cov provides a better way to estimate the variance of the identified modal parameters.The idea of combining EMD and SSI for modal identification was first proposed by Peeters and De Roeck [21]. Since the IMFs are single component signals, the stabilization diagram obtained from EMD-SSI method only contains one stable pole, which will greatly reduce the human interaction work during postprocessing. However, the EMD-SSI method depends largely on band-pass filter, which is introduced to deal with the mode mixing problem of EMD. This means one needs a priori knowledge (frequency components) of the vibration response before modal analysis. With the improved EMD, the problem of mode mixing can be solved without the so-called a priori knowledge of the signals. This paper will use SSI-cov algorithm together with IMFs extracted by bandwidth restricted EMD to identify modal parameters.


## 3.1. Random Decrement Technique

RDT is a time domain signal processing approach that is developed for extracting free decay response from stationary linear structure vibration response. Like NExT and LSCE method mentioned in Section1, RDT is also a time domain technique. These modal identification techniques, including RDT, NExT, and LSCE, are also capable of dealing with closely spaced modes [18, 19]. However, EMD based RDT modal identification technique is more suitable for large-scale structures since both excitations and responses are nonstationary processes.The basic principles of RDT are described as follows. For a SDOF linear vibration system, the displacement responsex
(
t
) is composed of three parts. The first part results from initial displacement, the second is from initial velocity, and the third is from random excitations. Thus, x
(
t
) can be expressed as(6)
x
t
=
x
0
D
t
+
x
˙
0
V
t
+
∫
0
t
h
t
-
τ
f
τ
d
τ
,where x
(
0
) and x
˙
(
0
) are the initial displacement and initial velocity, respectively. D
(
t
) is the free vibration response of the system with initial displacement being 1 and initial velocity being 0. V
(
t
) is the free vibration response with initial displacement being 0 and initial velocity being 1. f
(
τ
) is the random excitations at time segment τ. h
(
t
) is unit impulse response function of SDOF system. If the stochastic part can be separated and eliminated from x
(
t
), one can expect the free vibration part of x
(
t
). Select a reasonable value A to intercept the displacement response x
(
t
), the crossing points are denoted as (
t
i
,
x
(
t
i
)
)
(
i
=
1,2
,
3
,
…
). A series of subsample functions x
(
t
-
t
i
) can be obtained through moving the start time to the crossing points,(7)
x
t
-
t
i
=
x
t
i
D
t
-
t
i
+
x
˙
t
i
V
t
-
t
i
+
∫
t
i
t
h
t
-
τ
f
τ
d
τ
.In output-only modal analysis, the random excitations are assumed to be stationary Gaussian white noise. Therefore, the change of start time will not change the stochastic characteristic of random excitation functions. Thus, a series of random process can be obtained from(8)
x
i
t
=
A
D
t
+
x
˙
t
i
V
t
+
∫
0
t
h
t
-
τ
f
τ
d
τ
.Sincex
˙
t
i is a zero-mean stationary Gaussian process, hence, E
x
˙
t
i
=
0 and E
f
t
=
0. The mathematical expectation E
x
i
t is obtained by(9)
E
x
i
t
=
A
D
t
+
E
x
˙
t
i
V
t
+
∫
0
t
h
t
-
τ
E
f
τ
=
A
D
t
,which represents the free vibration response with initial displacement being A and initial velocity being 0.After obtaining the free decays of IMF, the modal parameter can be estimated by Hilbert transform. Without loss of generality, the free vibration response in (9) can be rewritten as(10)
x
t
=
A
0
e
-
ξ
ω
0
t
cos
⁡
ω
d
t
+
ϕ
0
,where ω
0 is the circular frequency and ω
d is the damping frequency. ξ is the damping ratio and A
0 is a constant. The corresponding analytical expression of x
(
t
) is obtained by(11)
z
t
=
x
t
+
i
x
~
t
=
A
t
e
-
i
θ
t
,where x
~
(
t
) is the Hilbert transform of x
(
t
). The amplitude function A
(
t
) and phase angle function θ
t can be expressed as(12)
A
t
=
A
0
e
-
ξ
ω
0
t
θ
t
=
ω
d
t
+
ϕ
0
.By introducing logarithm and differential operator, (12) yields(13)
ln
⁡
A
t
=
ln
⁡
A
0
-
ξ
ω
0
t
ω
t
=
d
θ
t
d
t
=
ω
d
.Equation (13) shows that the damping frequency is calculated from the slope of phase angle function. The slope of amplitude function ξ
ω
0 can also be obtained. Thus, the frequency ω
0 and damping ratio ξ are calculated from ω
d
=
ω
0
1
-
ξ
2.


## 3.2. Stochastic Subspace Identification

SSI is a robust time domain modal identification approach. Many successful applications based on SSI have been reported. The stochastic state space model is described as(14)
x
k
+
1
=
A
x
k
+
w
k
y
k
=
C
x
k
+
v
k
,where x
k is the state vector; y
k is the output vector; w
k and v
k denote the modelling inaccuracies and measuring errors, respectively; A is the state matrix and C output matrix. The main objective of SSI is to identify A and C from output-only response. The modal parameters can be obtained by applying eigenvalue decomposition to A. Based on the state space model, there are mainly two implementations of SSI algorithm: the first one directly starts from the sampled vibration data, namely, data driven SSI (SSI-data) and the other one from the covariance of the vibration data which is called covariance driven SSI (SSI-cov). Both SSI-data and SSI-cov reveal similar results with satisfying accuracy, while SSI-cov provides a better way to estimate the variance of the identified modal parameters.The idea of combining EMD and SSI for modal identification was first proposed by Peeters and De Roeck [21]. Since the IMFs are single component signals, the stabilization diagram obtained from EMD-SSI method only contains one stable pole, which will greatly reduce the human interaction work during postprocessing. However, the EMD-SSI method depends largely on band-pass filter, which is introduced to deal with the mode mixing problem of EMD. This means one needs a priori knowledge (frequency components) of the vibration response before modal analysis. With the improved EMD, the problem of mode mixing can be solved without the so-called a priori knowledge of the signals. This paper will use SSI-cov algorithm together with IMFs extracted by bandwidth restricted EMD to identify modal parameters.


## 4. Application: Output-Only Modal Analysis of a Railway Bridge

To illustrate the effectiveness of the proposed method, output-only modal analysis of a railway bridge is performed in this section. The modal identification results from Peak-Picking and SSI-cov are also listed for comparison. Before modal analysis, some basic information about the railway bridge and the ambient vibration test are introduced.


### 4.1. Bridge Description

Songtoujiang bridge (Figure6) is a four span prestressed concrete bridge located at the mountain area in Jiangxi province, China. It is a continuous box girder bridge with a span arrangement of (
60
+
2
×
100
+
60
) m. One of its piers is as high as 99 m. The width of the bridge deck is 12.5 m, on which located two opposite parallel railway lines.Figure 6
Elevation, plan, and layout of measuring points in ambient vibration test. The circles indicate reference nodes (dimensions are in centimeter).


### 4.2. Ambient Vibration Test

In order to know its dynamic properties and also gather useful information for future establishment of long-term SHM system, an ambient vibration test campaign of Songtoujiang bridge was performed. The test campaign was based on a low-cost data acquisition system with 6 channels. In total, 3 accelerometers were used for a multi-setups test with 2 reference nodes and 47 roving measuring points. As shown in Figure6, all measuring points were placed at the center of the bridge deck, and the number in circles shows the location of reference sensors.In each setup, vertical and lateral accelerations of bridge deck at two reference sensors and one roving sensor were recorded. The sampling frequency is 20 Hz and the sampling time is around 5 minutes for every setup. Figure7 shows typical measured acceleration time histories at test point 5. Before output-only modal analysis was performed, the test signals were first preprocessed with (
1
) DC components of the test signal were removed. (
2
) Low-pass filter was applied to test signals in order to remove high frequency noises.Figure 7
Vertical and lateral acceleration responses recorded at measuring point 5.


### 4.3. Output-Only Modal Analysis

In this section, the output-only modal analysis of Songtoujiang bridge is performed. In total, three methods are used for comparison: Peak-Picking, SSI-cov, and the improved EMD based modal identification method.Peak-Picking is a frequency domain method aimed at determining the Frequency Response Function (FRF) of a structure. The peaks in the FRF figure represent the eigenfrequencies and the damping ratios are estimated by half-power spectrum bandwidth method. FRF is estimated by the input and output measurements of the structure. In the output-only case, the FRF cannot be determined but only instead by the Fourier spectrum of the output response. Figure8 shows the vertical and lateral spectrum of acceleration responses recorded at measuring point 5.Figure 8
The Fourier spectrum of the accelerations shown in Figure4. (a) Vertical spectrum; (b) lateral spectrum.
(a)
(b)Figure9 plots the stabilization diagram of setup 5 obtained from SSI-cov algorithm. There are many stable poles in the figure. One need to decide which stable poles represent system modes, which are spurious modes. This kind of human interaction could be very time-consuming; especially there are 47 stabilization diagrams in total. Besides some stable poles are quite close, it is hard to choose system modes. For example, there are two stable poles near 1.5 Hz; both of them could be system poles.Figure 9
The stabilization diagram obtained from conventional SSI-cov. The stars denote the stable points with following stable criterion:d
f
i
≤
1
%, d
ξ
i
≤
5
%, and d
M
A
C
i
≤
1
%. MAC is the Modal Assurance Criterion value of two adjacent mode shapes.For the improved EMD method, first the acceleration responses should be decomposed into IMFs. Take the lateral acceleration response at measuring point 5 as an example. Figure10 shows the first six IMFs of the lateral acceleration response. It is shown that the frequencies of the IMFs are decreasing. After the IMFs were extracted, the modal parameter can be identified by either RDT or SSI-cov.Figure 10
The first six IMFs obtained by applying bandwidth restricted EMD to the lateral acceleration at measuring point 5.Figure11 shows the free decay response extracted from IMF by RDT. The corresponding amplitude function and phase angle function are shown in Figure 12. Through least square fitting, the slopes of the amplitude function and phase angle function were calculated. The slope of amplitude function -
ξ
ω
0
=
-
9.473 rad/s. The slope of phase angle function ω
d
=
ω
0
1
-
ξ
2
=
0.102. Solving the two equations simultaneously yields the frequency ω
0
=
1.509Hz and the damping ratio ξ
=
0.068. This mode is the fourth lateral bending mode. Following the same rules, the frequencies and damping ratios of the left modes can be obtained.Figure 11
The free decay response extracted from IMF (BR-C
4 in Figure 10) by random decrement technique.Figure 12
Phase angle and amplitude function of the free decay response in Figure6. (a) Phase angle and its least square fitting. (b) The amplitude logarithm and its least square fitting.
(a)
(b)Assemble the IMFs contain the same frequency component from different setups and through SSI-cov algorithm, the stabilization diagrams in Figure13 were obtained. In these stabilization diagrams, there is only one stable pole. The spurious modes are separated into some discrete stable points. The system modes can be easily determined. The Fourier spectrum of each IMF is also shown in the stabilization diagram which further illustrated that IMF obtained from the improved EMD is single frequency component signal.Figure 13
The stabilization diagrams obtained from bandwidth restricted EMD based SSI-cov algorithm. (a) Mode 1:f
=
0
.
583 Hz; (b) mode 2: f
=
0
.
716 Hz; (c) mode 3: f
=
0
.
915 Hz; (d) mode 5: f
=
1
.
512 Hz. The stars denote the stable points with following stable criterion: d
f
i
≤
1
%, d
ξ
i
≤
5
%, and d
M
A
C
i
≤
1
%. The fourier spectrums are scaled to sever as comparison.
(a)
(b)
(c)
(d)


### 4.4. Discussion

Table1 shows the modal parameters of Songtoujiang Bridge identified from different approaches. The experimental results from PP, SSI-cov, and standard EMD are listed for comparison. The analytical results from a preliminary Finite Element Model (FEM) are also listed for comparison. In total, four lateral bending modes and two vertical bending modes can be extracted from operational vibration measurements. Because of the relatively low sampling frequency and short sampling time of the low-cost data acquisition system, there might exist some inaccuracies in the higher modes. From Table 1, it is shown that the experimental results from different methods show good agreement, while the analytical results show some deviations especially in 2nd and 3rd lateral bending mode. This means that the preliminary FE model need to be updated. A higher accuracy of frequencies is found than the damping ratios. Due to mode mixing, the third lateral bending mode cannot be extracted from standard EMD. With the improved EMD method, the natural frequencies and damping ratios can be extracted with satisfactory accuracy.Table 1
The comparison of the modal parameters obtained from different algorithms.
fem
PP
SSI-cov
Std. EMD-RDT
Improved EMD
Mode
with RDT
with SSI-cov
f
f (Hz)
f
PP (Hz)
f
SSI (Hz)
ξ
SSI (%)
f
e (Hz)
ξ
e (%)
f
e
r (Hz)
ξ
e
r (%)
f
e
s (Hz)
ξ
e
s (%)
1st↔
0.662
0.58
0.583
2.74
0.530
7.41
0.530
3.12
0.583
2.28
2nd↔
0.789
0.71
0.721
3.63
0.714
7.23
0.714
2.85
0.716
3.52
3rd↔
1.137
0.89
0.914
3.91
n/a
n/a
0.916
1.96
0.915
4.33
1st ↕
1.333
1.32
1.302
1.54
1.319
1.28
1.331
1.68
1.325
2.48
4th↔
1.570
1.50
1.481/1.508
1.65/1.26
1.505
2.84
1.509
1.37
1.512
1.03
2nd ↕
1.975
2.09
1.996
4.82
2.097
0.83
1.99
4.23
1.983
3.68
Note:↔ denotes lateral bending mode; ↕ denotes vertical bending mode; n/a: not a number.The main contribution of the proposed improved EMD based modal identification technique is that it provide more reasonable theory explanations to deal with nonstationary excitations (ambient loads) and measurements (dynamic responses) without assuming them to be stationary processes. Although some subspace algorithms were reported to have constant estimates of eigenstructures even under nonstationary excitations [30, 31], the theory derivation of subspace algorithms is based on stationary inputs assumption. However, the EMD is proved to be very suitable for nonstationary signals, even nonlinear signals. Another contribution of the proposed method is that band-pass filter, which is often employed to deal with closely spaced modes, is no longer needed with bandwidth restricted EMD. The application example demonstrates that quite accurate modal parameters can be obtained by the proposed method, meaning that the improved EMD based modal identification technique can serve as an alternative to existing output-only modal analysis methods.


## 4.1. Bridge Description

Songtoujiang bridge (Figure6) is a four span prestressed concrete bridge located at the mountain area in Jiangxi province, China. It is a continuous box girder bridge with a span arrangement of (
60
+
2
×
100
+
60
) m. One of its piers is as high as 99 m. The width of the bridge deck is 12.5 m, on which located two opposite parallel railway lines.Figure 6
Elevation, plan, and layout of measuring points in ambient vibration test. The circles indicate reference nodes (dimensions are in centimeter).


## 4.2. Ambient Vibration Test

In order to know its dynamic properties and also gather useful information for future establishment of long-term SHM system, an ambient vibration test campaign of Songtoujiang bridge was performed. The test campaign was based on a low-cost data acquisition system with 6 channels. In total, 3 accelerometers were used for a multi-setups test with 2 reference nodes and 47 roving measuring points. As shown in Figure6, all measuring points were placed at the center of the bridge deck, and the number in circles shows the location of reference sensors.In each setup, vertical and lateral accelerations of bridge deck at two reference sensors and one roving sensor were recorded. The sampling frequency is 20 Hz and the sampling time is around 5 minutes for every setup. Figure7 shows typical measured acceleration time histories at test point 5. Before output-only modal analysis was performed, the test signals were first preprocessed with (
1
) DC components of the test signal were removed. (
2
) Low-pass filter was applied to test signals in order to remove high frequency noises.Figure 7
Vertical and lateral acceleration responses recorded at measuring point 5.


## 4.3. Output-Only Modal Analysis

In this section, the output-only modal analysis of Songtoujiang bridge is performed. In total, three methods are used for comparison: Peak-Picking, SSI-cov, and the improved EMD based modal identification method.Peak-Picking is a frequency domain method aimed at determining the Frequency Response Function (FRF) of a structure. The peaks in the FRF figure represent the eigenfrequencies and the damping ratios are estimated by half-power spectrum bandwidth method. FRF is estimated by the input and output measurements of the structure. In the output-only case, the FRF cannot be determined but only instead by the Fourier spectrum of the output response. Figure8 shows the vertical and lateral spectrum of acceleration responses recorded at measuring point 5.Figure 8
The Fourier spectrum of the accelerations shown in Figure4. (a) Vertical spectrum; (b) lateral spectrum.
(a)
(b)Figure9 plots the stabilization diagram of setup 5 obtained from SSI-cov algorithm. There are many stable poles in the figure. One need to decide which stable poles represent system modes, which are spurious modes. This kind of human interaction could be very time-consuming; especially there are 47 stabilization diagrams in total. Besides some stable poles are quite close, it is hard to choose system modes. For example, there are two stable poles near 1.5 Hz; both of them could be system poles.Figure 9
The stabilization diagram obtained from conventional SSI-cov. The stars denote the stable points with following stable criterion:d
f
i
≤
1
%, d
ξ
i
≤
5
%, and d
M
A
C
i
≤
1
%. MAC is the Modal Assurance Criterion value of two adjacent mode shapes.For the improved EMD method, first the acceleration responses should be decomposed into IMFs. Take the lateral acceleration response at measuring point 5 as an example. Figure10 shows the first six IMFs of the lateral acceleration response. It is shown that the frequencies of the IMFs are decreasing. After the IMFs were extracted, the modal parameter can be identified by either RDT or SSI-cov.Figure 10
The first six IMFs obtained by applying bandwidth restricted EMD to the lateral acceleration at measuring point 5.Figure11 shows the free decay response extracted from IMF by RDT. The corresponding amplitude function and phase angle function are shown in Figure 12. Through least square fitting, the slopes of the amplitude function and phase angle function were calculated. The slope of amplitude function -
ξ
ω
0
=
-
9.473 rad/s. The slope of phase angle function ω
d
=
ω
0
1
-
ξ
2
=
0.102. Solving the two equations simultaneously yields the frequency ω
0
=
1.509Hz and the damping ratio ξ
=
0.068. This mode is the fourth lateral bending mode. Following the same rules, the frequencies and damping ratios of the left modes can be obtained.Figure 11
The free decay response extracted from IMF (BR-C
4 in Figure 10) by random decrement technique.Figure 12
Phase angle and amplitude function of the free decay response in Figure6. (a) Phase angle and its least square fitting. (b) The amplitude logarithm and its least square fitting.
(a)
(b)Assemble the IMFs contain the same frequency component from different setups and through SSI-cov algorithm, the stabilization diagrams in Figure13 were obtained. In these stabilization diagrams, there is only one stable pole. The spurious modes are separated into some discrete stable points. The system modes can be easily determined. The Fourier spectrum of each IMF is also shown in the stabilization diagram which further illustrated that IMF obtained from the improved EMD is single frequency component signal.Figure 13
The stabilization diagrams obtained from bandwidth restricted EMD based SSI-cov algorithm. (a) Mode 1:f
=
0
.
583 Hz; (b) mode 2: f
=
0
.
716 Hz; (c) mode 3: f
=
0
.
915 Hz; (d) mode 5: f
=
1
.
512 Hz. The stars denote the stable points with following stable criterion: d
f
i
≤
1
%, d
ξ
i
≤
5
%, and d
M
A
C
i
≤
1
%. The fourier spectrums are scaled to sever as comparison.
(a)
(b)
(c)
(d)


## 4.4. Discussion

Table1 shows the modal parameters of Songtoujiang Bridge identified from different approaches. The experimental results from PP, SSI-cov, and standard EMD are listed for comparison. The analytical results from a preliminary Finite Element Model (FEM) are also listed for comparison. In total, four lateral bending modes and two vertical bending modes can be extracted from operational vibration measurements. Because of the relatively low sampling frequency and short sampling time of the low-cost data acquisition system, there might exist some inaccuracies in the higher modes. From Table 1, it is shown that the experimental results from different methods show good agreement, while the analytical results show some deviations especially in 2nd and 3rd lateral bending mode. This means that the preliminary FE model need to be updated. A higher accuracy of frequencies is found than the damping ratios. Due to mode mixing, the third lateral bending mode cannot be extracted from standard EMD. With the improved EMD method, the natural frequencies and damping ratios can be extracted with satisfactory accuracy.Table 1
The comparison of the modal parameters obtained from different algorithms.
fem
PP
SSI-cov
Std. EMD-RDT
Improved EMD
Mode
with RDT
with SSI-cov
f
f (Hz)
f
PP (Hz)
f
SSI (Hz)
ξ
SSI (%)
f
e (Hz)
ξ
e (%)
f
e
r (Hz)
ξ
e
r (%)
f
e
s (Hz)
ξ
e
s (%)
1st↔
0.662
0.58
0.583
2.74
0.530
7.41
0.530
3.12
0.583
2.28
2nd↔
0.789
0.71
0.721
3.63
0.714
7.23
0.714
2.85
0.716
3.52
3rd↔
1.137
0.89
0.914
3.91
n/a
n/a
0.916
1.96
0.915
4.33
1st ↕
1.333
1.32
1.302
1.54
1.319
1.28
1.331
1.68
1.325
2.48
4th↔
1.570
1.50
1.481/1.508
1.65/1.26
1.505
2.84
1.509
1.37
1.512
1.03
2nd ↕
1.975
2.09
1.996
4.82
2.097
0.83
1.99
4.23
1.983
3.68
Note:↔ denotes lateral bending mode; ↕ denotes vertical bending mode; n/a: not a number.The main contribution of the proposed improved EMD based modal identification technique is that it provide more reasonable theory explanations to deal with nonstationary excitations (ambient loads) and measurements (dynamic responses) without assuming them to be stationary processes. Although some subspace algorithms were reported to have constant estimates of eigenstructures even under nonstationary excitations [30, 31], the theory derivation of subspace algorithms is based on stationary inputs assumption. However, the EMD is proved to be very suitable for nonstationary signals, even nonlinear signals. Another contribution of the proposed method is that band-pass filter, which is often employed to deal with closely spaced modes, is no longer needed with bandwidth restricted EMD. The application example demonstrates that quite accurate modal parameters can be obtained by the proposed method, meaning that the improved EMD based modal identification technique can serve as an alternative to existing output-only modal analysis methods.


## 5. Conclusions

Starting from the nonstationary vibration measurements and closely spaced mode of bridge structures, this paper has developed a bandwidth restricted EMD for output-only modal analysis. Specifically, the mode mixing problem of standard EMD when dealing with closely spaced modes has been effectively restrained by the improved EMD. Through numerical simulation and application in a railway bridge, the following conclusions have been made:(1)
The bandwidth restricted EMD have restrained the problem of mode mixing, making sure that each IMF obtained from EMD contains only one frequency component. It is based on the binary filter property of EMD and can be easily implemented, without introducing too much computational complexities. The numerical simulation have illustrated the effectiveness of the improved EMD.
(2)
The modal parameters obtained through improved EMD algorithm show good agreement with those extracted from PP and SSI-cov. Compared with conventional EMD based RDT or SSI-cov modal identification algorithms, the improved EMD based method does not rely on band-pass filter when dealing with closely spaced modes, meaning that the a priori knowledge (frequency components) of nonstationary output measurements is not needed.



---

*Source: 945862-2015-10-04.xml*