Structural nonlinear damage detection using improved DempsterShafer theory and time domain model
Huiyong Guo^{1} , Rong Zhou^{2} , Feng Zhang^{3}
^{1, 2, 3}School of Civil Engineering, Chongqing University, Chongqing, 400045, P. R. China
^{1, 2, 3}Laboratory of New Technology for Construction of Cities in Mountain Area, Ministry of Education, Chongqing, 400045, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 21, Issue 6, 2019, p. 16791693.
https://doi.org/10.21595/jve.2019.20858
Received 14 June 2019; received in revised form 6 September 2019; accepted 16 September 2019; published 30 September 2019
JVE Conferences
In the service period, a crack may appear in some engineering structures. The development of accurate and effective methods for crack damage detection has become a topic of great importance. In this paper, a nonlinear damage detection method based on the improved DempsterShafer (DS) theory and time domain model is presented. First, acceleration responses in the undamaged and damaged states are measured by using accelerometers. Then, acceleration responses are utilized to establish an autoregressive (AR) model, and residual time series of acceleration responses are used to establish an autoregressive conditional heteroskedasticity (ARCH) model. A cepstral metric conversion (CMC) method based on the AR model is employed to obtain local damage solution and an autoregressive conditional heteroskedasticity conversion (ARCHC) method based on ARCH model is presented to acquire another local damage solution. Finally, the DS theory is applied to detect damages by integrating these local damage solutions, and an improved DS theory is further presented to enhance the detection accuracy. The numerical and experimental examples show that the improved DS theory has high detection accuracy and good performance.
 An efficient method based on DempsterShafer theory and time domain model is proposed to detect crack damage.
 A cepstral metric conversion method based on the AR model is employed to obtain local damage solution.
 An autoregressive conditional heteroskedasticity conversion method based on ARCH model is presented to acquire another local damage solution.
 An improved probability quantity method based on different importance is presented to enhance the detection accuracy of the DS theory.
Keywords: damage detection, DempsterShafer theory, time domain model, nonlinearity, acceleration response.
1. Introduction
Vibrationbased damage diagnosis and detection have received considerable attention in recent years [1]. The feasibility of vibrationbased damage detection is that the damage usually changes the structural characteristics, such as stiffness, which are reflected in the measured vibration characteristics of the structure [2]. The vibration characteristics include measured frequency, mode shape and strain energy, which are often used in damage detection [35]. One category of vibrationbased damage detection methods is the application of time domain analysis by constructing a time series model, which has been widely used in feature exaction to realize structural damage diagnosis and detection. Muślewski et al. [6] analyzed the vibration characteristics using vibration time histories in the time domain for propulsion systems of ships. Grządziela et al. [7] proposed a noncoaxiality identification method for shaft line systems of engines propelling the ships. Das et al. [8] made a comparative study among the different vibrationbased damage detection methods: modal examination, local diagnostic approach, nonprobabilistic methodology and time series method, and the results showed that the time series method results were more credible in damage detection than the other methods.
Generally, the time domain damage can be divided into two categories: linear damage and nonlinear damage. The linear damage got more attention and AR, ARX and ARMA models were widely applied in damage detection by fitting the response data of linear systems [912]. Zheng et al. [13] proposed a damage indicator based on the distance between ARMA models to detect structural damage including its location and severity. Xing et al. [14] utilized loworder AR models and cepstral metric to find damage locations. The principle of the method is that the distance between the models of damaged and undamaged states is correlated with the damage location information. Liu et al [15] proposed a nonprobabilistic damage detection method, which employed the coefficient matrix of vector autoregressive model to extract diagnosis elements and used Mahalanobis distance of these diagnosis elements to detect damages. Mei et al. [12] presented a substructurebased damage identification method to find damage locations and severity, and autoregressive moving average with exogenous variables model residual was utilized to correct the damage indicator.
Most of the damages in engineering structures actually possesses nonlinear characteristic, such as fatigue cracks that open and close due to dynamic loading. Although some methods and theories have been utilized for structural nonlinear damage detection [1618], there are no damage detection methods with high reliability and sufficient research depth. The ARCH model was proposed by Engle [19], and subsequently was studied and applied by many researchers in the field of economics and finance. Pham et al [20] used a hybrid ARMA/GARCH model to diagnose the machine health state based on vibration signals. Chen et al. [21] also utilized the hybrid ARMA/GARCH model to diagnose the nonlinear damage of a threestorey structure and compared this model with the linear ARMA model. Yu et al. [22] presented a combination method based on the time domain model and statistical moments of structural responses to extract nonlinear damage information. Cheng et al [23] proposed the second order variance indicator based on the AR/ARCH model to identify nonlinear damage. These methods are mainly based on the degree of freedom (DOF) to identify structural damage locations, and it is still necessary to improve the identification reliability of nonlinear damage. Therefore, in this paper, a conversion technique is proposed to detect interstorey stiffness damage and an improved DS theory based on time domain model is presented to raise the identification reliability.
In this study, a nonlinear damage detection method based on the DS theory and time domain model is proposed. First, the ARCH model is analyzed, and an autoregressive conditional heteroskedasticity conversion method is proposed. Then, the cepstral metric theory based on the autoregressive model is given, and a cepstral metric conversion method is described. Finally, the DS theory is used to realize damage detection, and an improved strategy is presented. This paper consists of seven sections. In Section 2, the ARCH model is described, and the ARCHC method is proposed to localize nonlinear damage. In Section 3, the cepstral metric is analyzed, and the CMC method is given. The improved DS theory and its application to damage detection are proposed in Section 4. Numerical and experimental examples are provided in Sections 5 and 6, respectively. Finally, Section 7 summarizes the conclusions.
2. ARCH model and ARCHC method
2.1. Nonlinear ARCH model
ARCH model is a volatility model that uses observed values of returns or residuals as volatility shocks. The basic theory of ARCH model is: Let $\left\{{y}_{t}\right\}$ denote a discrete time series set and let ${\psi}_{t1}$ denote the information set available at time $t1$. The prediction value at time $t$ can be acquired by using the conditional expectation $E\left\{{y}_{t}{\psi}_{t1}\right\}$. Thus, the prediction residual or error ${\epsilon}_{t}$ can be expressed as:
Then, the conditional variance ${v}_{t}$ at time $t$ can be written as:
The values of $r$ past squared errors can be used to form the conditional variance function as follows:
In general, the error ${\epsilon}_{t}$ can be acquired from measured nonlinear time series $\left\{{y}_{t}\right\}$, and the conditional expectation can be obtained by an autoregressive analysis. Thus, the Eq. (1) can be written as:
where $\delta $ is a constant, ${\phi}_{i}$ is the autoregressive coefficient, and $p$ is the number of autoregressive order. Similarly, the conditional variance ${v}_{t}$ of $r$ past squared errors can be specified as a linear function as follows:
where ${\varsigma}_{0}$ and ${\gamma}_{j}$ ($j=$1,…, $r$) are the parameters of ARCH model, ${z}_{t}$ is a zeromean unitvariance white noise. The values of the parameters are constrained by ${\varsigma}_{0}>0$, ${\gamma}_{j}\ge 0$, $j=1,\cdots ,r$, ${\sum}_{j=1}^{r}{\gamma}_{j}<1$.
The first two constraints ensure that conditional variances are positive, and the third constraint ensures covariance stationarity.
It is necessary to determine the order of ARCH model and estimate model parameters. The maximum likelihood estimation method can be used for estimating the parameters of an ARCH model. Akaike information criterion (AIC) can be used to determine the order of ARCH model. It utilizes a function based on the likelihood to choose the best order. The mathematical form of AIC can be expressed as:
where $L$ is the value of likelihood function, and $m$ is the number of estimated parameters. AIC represents a tradeoff between the loglikelihood value and the number of parameters. In practice, we prefer to select the models with smaller AIC values.
2.2. ARCHC method
Generally, the structural damage identification can be implemented by comparing some features of a structure in different states, (i.e., undamaged state and damaged state). The structural acceleration time series in the undamaged and damaged states must be acquired from simulated or experimental vibration structure. The main steps of ARCH model application to establish the damage index are summarized as follows:
(1) Acquire the acceleration time series data in undamaged state and damaged state, so that the acquired time series meet stationarity.
(2) Extract residual sequences from the above acceleration data, and check if the residual sequences have ARCH effects.
(3) Estimate model parameters and determine the order of ARCH models. And then construct the ARCH models using the residual sequences $\left\{{\epsilon}_{t}^{B}\right\}$ and $\left\{{\epsilon}_{t}^{T}\right\}$ in the undamaged and damaged states.
(4) Acquire the conditional variance sequences $\left\{{{v}_{t}}^{B}\right\}$ and $\left\{{{v}_{t}}^{T}\right\}$ of the residual sequences $\left\{{\epsilon}_{t}^{B}\right\}$ and $\left\{{\epsilon}_{t}^{T}\right\}$ in the undamaged and damaged states, and calculate the corresponding variance $s\left({{v}_{t}}^{B}\right)$ and $s\left({{v}_{t}}^{T}\right)$.
(5) Establish autoregressive conditional heteroskedasticity conversion (ARCHC) method based on calculated values in the undamaged and damaged states to identify the nonlinear damage. The ARCHC method is described as follows:
For an $n$storey structure with nonlinear damage, the residual sequences $\left\{{\epsilon}_{t}\right\}$ of acceleration responses are no longer white noise processes, i.e., conditional variances sequences $\left\{{v}_{t}\right\}$ vary with time. The variance of conditional variances sequences $\left\{{v}_{t}\right\}$ is much greater than zero. On the other hand, for an undamaged structure, the conditional variance sequences are stable. Therefore, nonlinear damage can be diagnosed by using some methods based on the ARCH model. Here, the ARCHC method is proposed to directly identify damage locations. The damage probability value of the $j$th stiffness storey is given by:
where:
When $i=$1, ${R}_{i1}={R}_{0}=0.5{R}_{1}$. Thus, the bottom effect is considered. $n$ is the stiffness storey number of shear structure. The ARCHC method can directly identify structural nonlinear damage location or provide local decision information for the following fusion.
3. Cepstral metric and CMC method
3.1. Cepstral metric of AR model
Martin [24] proposed cepstral metric of ARMA models. Xing et al [14] further developed the Martin’s method and proposed cepstral metric of AR model. An autoregressive model can be written as:
where ${\phi}_{i}$ is the coefficient of AR model, $p$ is the order of AR model, and ${\epsilon}_{t}$ is a white noise time sequence. The cepstral metric of AR model is described as follows [14]:
where $N$ and $N\mathrm{\text{'}}$ are AR models before and after damage, respectively. $p$ and $p\mathrm{\text{'}}$ are the order of models $N$ and $N\mathrm{\text{'}}$. ${\alpha}_{i}$ and $\alpha \mathrm{\text{'}}$ are the pole parameters of models $N$ and $N\mathrm{\text{'}}$.
3.2. CMC method
In general, it is difficult to identify stiffness damage through the cepstral metric based on DOFs. Therefore, it is necessary to establish the damage detection method based on the stiffness storey. Here, the CMC method is proposed to directly identify damage locations. The damage probability value of the $j$th stiffness storey can be written as:
where:
When $i=$1, ${Q}_{i1}={Q}_{0}=0.5{Q}_{1}$. Thus, the bottom effect is considered. The CMC method can directly identify structural damage location.
4. Improved DS theory and its application to damage detection
4.1. DS theory and its application to damage detection
DempsterShafer theory is an effective information fusion approach for combining multiple data sources [25]. It has been widely utilized in the fields of pattern recognition.
Let $\mathrm{\Theta}$ be a finite set including $n$ hypotheses. The set $\mathrm{\Theta}$ is the frame of discernment. Basic probability quantity $m$ of focal element should meet the following constraint:
where $A$ expresses the focal element or evidence. The sum of basic probability quantity is 1. Dempster’s combination rule is the key part of DS theory, which can combine multisource data. Here, two sets of data sources ${A}_{1}$ and ${A}_{2}$ are considered. Let ${m}_{1}\left({A}_{1}\right)$ and ${m}_{2}\left({A}_{2}\right)$ be basic probability quantity given by sources ${A}_{1}$ and ${A}_{2}$, respectively. The combination rule can be written as:
where:
where ${A}_{1}$ and ${A}_{2}$ are focal elements. $r$ is conflict factor between the data sources ${A}_{1}$ and ${A}_{2}$. The item ($1r$) is utilized to compensate for probability loss caused by nonintersecting subsets. Thus, the generated probability quantity $m\left(B\right)$ is normalized to one.
In this study, two different data sources, the residual time series data and autoregressive time series data of acceleration responses are applied to damage detection. The local decisions can be acquired by using ARCHC method and CMC method, respectively. Then, the basic DS theory can be used to obtain the global solution and localize damages.
4.2. Improved DS theory and its application to damage detection
Generally, different information sources shall have different importance or reliability, but the basic DS theory uses the equal importance strategy to deal with the fusion problem. Therefore, an improved strategy, which considers the effect of different importance, is proposed.
Consider two information sources ${A}_{h}$ and ${A}_{l}$ with different importance. Let ${m}_{h}$ and ${m}_{l}$ be the basic probability quantity given by sources ${A}_{h}$ and ${A}_{l}$ and let ${\alpha}_{h}$ and ${\alpha}_{l}$ be the probability importance coefficients which are the probability importance degree of sources ${A}_{h}$ and ${A}_{l}$, respectively. Here, ${A}_{h}$ has a higher importance degree than ${A}_{l}$, i.e., ${\alpha}_{h}\ge {\alpha}_{l}$. In addition, the probability importance coefficients should meet the following constraint:
The ${\alpha}_{h}{m}_{h}$ and ${\alpha}_{l}{m}_{l}$ could not be directly utilized as modified probability quantity of sources ${A}_{h}$ and ${A}_{l}$. The reason is the ${\alpha}_{h}{m}_{h}$ or ${\alpha}_{l}{m}_{l}$ could not meet the probability condition that the sum is 1. Here, an improved probability quantity method based on different importance is given as follows:
The generated probability quantity ${m}_{h}^{\text{'}}$ and ${m}_{l}^{\text{'}}$ all meet the normalization condition. Then, Dempster’s combination rule can be used by:
where $r=\sum _{{A}_{h}\cap {A}_{l}=\phi}{m}_{h}^{\text{'}}\left({A}_{h}\right)\times {m}_{l}^{\text{'}}\left({A}_{l}\right).$
If the importance degrees of sources ${A}_{h}$ and ${A}_{l}$ are equal (i.e., ${\alpha}_{h}={\alpha}_{l}=\text{0.5}$), the improved DS theory can automatically become basic DS theory.
In this study, residual time series data and autoregressive time series data of acceleration response are considered as two sources. In general, the residual time series data contain more nonlinear information than autoregressive time series data. Therefore, residual time series data is the source ${A}_{h}$ and autoregressive time series data is the source ${A}_{l}$. ${m}_{h}$ can be acquired by using the ARCHC method and ${m}_{l}$ can be obtained by using the CMC method. Nonlinear ARCH model of ARCHC method can perfectly deal with nonlinear signals. Therefore, for nonlinear damage problem, the probability importance coefficient ${\alpha}_{h}$ of ARCHC method is greater than ${\alpha}_{l}$ of CMC method. Probability importance coefficients can be acquired by the relative recognition rate of the damage detection method. Of course, the coefficients can also be obtained in practice.
5. Numerical examples
An eightstorey shear structure, which is shown in Fig. 1, is utilized to verify the effectiveness of the improved DS theory. This structure can be simplified as an eight degreeoffreedom system. The parameters of the shear structure are: mass of every storey is 100 kg, lateral stiffness is 1 MN/m, and Rayleigh damping parameters $\alpha $ and $\beta $ are 0.9842 and 0.00036. Nonlinear damage is simulated by means of the bilinear stiffness approach [26], which is given as follows:
where ${x}_{i}\left(t\right)$ and ${x}_{i1}\left(t\right)$ are the displacement of the $i$th and $\left(i1\right)$th storey, and ${k}_{i}\left[{x}_{i}\left(t\right)\right]$ denotes the $i$th interstorey stiffness. If $i=$1, ${x}_{i1}\left(t\right)={x}_{0}\left(t\right)=0$, expressing the displacement of structure basement. $c$ is damage coefficient. The bilinear stiffness approach is used to realize the simulation of breathing crack.
In this study, the damage coefficient has been assigned with two degrees of damage, i.e., 20 % and 30 %. Single damage occurs on the column from the 1st storey to 8th storey in turns. Thus, in total there are 16 damage scenarios. We can use the ‘stiffness storey’ to denote the interstorey stiffness damage location. White noise which is filtered by Butterworth lowpass filter has the similar features as wind load or small earthquake. Therefore, the filtered white noise data is utilized as excitation loads, and the Wilson$\theta $ method is adopted for calculating structural response. In addition, the disturbance of measurement noise cannot be ignored. Random noise, whose amplitude is adjusted to the same level as structural response, is added into the acceleration response time series as measurement noise. Thus, the acceleration time series with measurement noise can be written as:
where $\left\{{y}_{t}\right\}$ is acceleration time series with measurement noise, $\left\{{y}_{t}\mathrm{\text{'}}\right\}$ is acceleration time series calculated by the Wilson$\theta $ method, $S$ is noise level, and $\left\{{w}_{t}\right\}$ is a random noise which amplitude has been adjusted to the same level as structural response. The noise level $S$ is 5 %.
Fig. 1. Eightstorey shear structure
The ARCHC and CMC methods are utilized for damage detection. For the ARCHC method, the order of ARCH model should be determined. AIC is used for order determination of ARCH model, and the maximum likelihood estimation method is adopted for parameter estimation. Here, the acquired order $r$ of ARCH model is 5. For the CMC method, the Xing’s method [14] is used to obtain the order of AR model. The literature [14] already described the order determination method of AR model in detail. For the DS theory, residual time series data is considered as source ${A}_{h}$ and autoregressive time series data is considered as source ${A}_{l}$. So ${m}_{h}$ can be acquired by using the ARCHC method and ${m}_{l}$ can be obtained by using the CMC method. Then, the basic DS theory and improved DS theory can be utilized to realize damage detection. For the improved DS theory, the probability importance coefficients ${\alpha}_{h}$ and ${\alpha}_{l}$ are 0.7 and 0.3, respectively.
First, 20 % damage is considered. Single damage with 20 % stiffness reduction occurs from the 1st to 8th storey in turns. The CMC method, ARCHC method, DS theory (or basic DS theory) and improved DS theory are applied for damage detection. The greater the calculation value of stiffness storey is, the higher the probability of damage is. The detection results from the 1st storey to the 8th storey with 20 % damage are shown in Figs. 2 and 3. From Figs. 2 and 3, it can be observed that the two methods based on the DS theory exhibit better results than the CMC and ARCHC methods, except for the 1st storey. For the 1ststorey damage, the value of basic DS theory is lower than that of ARCHC method, but the result of improved DS theory is still higher than that of ARCHC method. Of course, for nonlinear damage problem, the ARCHC method is better than the CMC method and the calculated results of the improved DS theory are superior to those of the basic DS theory. Therefore, the improved DS theory is better than the other three methods.
Fig. 2. Detection results from the 1ststorey to 4thstorey damage with 20 % stiffness reduction
Fig. 3. Detection results from the 5thstorey to 8thstorey damage with 20 % stiffness reduction
Then, 30 % damage is also considered. Single damage with 30 % stiffness reduction occurs from the 1st to 8th storey in turns. The CMC method, ARCHC method, DS theory and improved DS theory are used to realize nonlinear damage detection. The detection results from the 1st storey to the 8th storey with 30 % damage are depicted in Fig. 4 and 5. From Figs. 4 and 5, it can be seen that the two methods based on the DS theory exhibit better results than both the CMC method and ARCHC method, except for the 1st storey. For the 1ststorey damage, the result of the basic DS theory is inferior to that of the ARCHC method, but the result of improved DS theory is obviously better than that of the ARCHC method. Therefore, the improved DS theory is the best one among the four damage detection methods.
Fig. 4. Detection results from the 1ststorey to 4thstorey damage with 30 % stiffness reduction
Fig. 5. Detection results from the 5thstorey to 8thstorey damage with 30 % stiffness reduction
For nonlinear damage problem, the CMC method utilizes autoregressive time series data to identify damage, and the ARCHC method mainly utilizes residual time series data to identify damage. In general, the residual time series data contains more nonlinear information than autoregressive time series data. Therefore, the identification results of the ARCHC method are better than those of CMC method. Then, DS theory is applied to further improve the identification results.
6. Experimental examples
A threestorey frame structure is utilized to verify the proposed method. The experimental frame structure with sensors is shown in Fig. 6 [27]. Nonlinear damage is modeled by a mechanism constituted of a bumper and a center column, which are used to generate nonlinear characteristic such as breathing crack prediction. The mechanism simulates the opening or closing of crack due to dynamic loads. The distance between the tip of the center column and the bumper is adjustable. In general, the greater the gap distance is, the weaker the nonlinear damage degree becomes. Therefore, in the experiment, nonlinear damage degree can be changed through adjusting the gap distance.
Fig. 6. Threestorey frame
a) Threestorey frame
b) Adjustable bumper and suspended column
Table 1. Sensor channels and parameters
Channel Number

Transducer

Reference

Nominal Sensitivity

Channel 1

Load cell

PCB 208 C03 SN 22569

2.2 mV/N

Channel 2

Accelerometer

PCB 336C SN 10099

1000 mV/g

Channel 3

Accelerometer

PCB 336C SN 10120

1000 mV/g

Channel 4

Accelerometer

PCB 336C SN 9916

1000 mV/g

Channel 5

Accelerometer

PCB 336C SN 10100

1000 mV/g

Table 2. Damaged states of threestorey frame structure
Damaged states

Description

State 1

Gap (0.05 mm)

State 2

Gap (0.10 mm)

State 3

Gap (0.13 mm)

State 4

Gap (0.15 mm)

State 5

Gap (0.20 mm) (weak nonlinearity)

Accelerometers are installed on each floor including the bottom floor to record the time domain acceleration responses. The sensor channels and parameters are depicted in Table 1. The acceleration responses of the 1st storey, 2nd storey and 3rd storey of the structure are the measured data of Channels 3, 4 and 5. The damaged states are listed in Table 2. The greater the gap distance is, the weaker the nonlinear damage degree becomes. If the gap distance is great enough, the bumper will not contact the center column and nonlinear characteristic caused by stiffness change will disappear. The acceleration responses in undamaged state are depicted in Fig. 7, and the acceleration responses in damaged state 1 are depicted in Fig. 8. It can be seen that their response curves are apparently different.
Fig. 7. Acceleration responses in undamaged State
Fig. 8. Acceleration responses in damaged State 1
Then, the order of time series model should be determined. Akaike information criterion is utilized to determine the order of ARCH model, and the maximum likelihood estimation method is adopted for parameter estimation. Here, the acquired order $r$ of ARCH model is 5. For the CMC method, the Xing’s method [14] is utilized to obtain the order of AR model. The order of AR model is 30. For the DS theory, residual time series data is considered as source ${A}_{h}$, and autoregressive time series data is considered as source ${A}_{l}$. And ${m}_{h}$ can be acquired by using the ARCHC method, and ${m}_{l}$ can be obtained by using the CMC method. Then, the basic DS theory and improved DS theory can be used to detect structural damage. For the improved DS theory, the probability importance coefficients ${\alpha}_{h}$ and ${\alpha}_{l}$ are 0.7 and 0.3, respectively.
For damaged state 1, measured acceleration response signals are utilized to realize nonlinear damage detection. The third stiffness storey is the damaged location. The CMC method, ARCHC method, DS theory (or basic DS theory) and improved DS theory are used to identify the structural damage location. The detection results for damaged state 1 are depicted in Fig. 9. From Fig. 9, it can be seen that all four methods possess the capacity of nonlinear damage detection. The two methods based on the DS theory have relatively higher values than the CMC and ARCHC methods do for the damaged third storey. The reason is that the methods based on the DS theory combine multiple source information. In addition, the calculated value of the improved DS theory is higher than that of the basic DS theory for the damaged third storey. Therefore, the improved DS theory is the best one among all four damage detection methods.
Fig. 9. Nonlinear damage detection results in damaged State 1
For damaged states 24, the damage is located in the third stiffness storey. The CMC method, ARCHC method, basic DS theory and improved DS theory are also utilized to detect structural damage location. The damage detection results with stiffness storey are depicted in Figs. 1012. From Figs. 1012, it can be seen that the two methods based on the DS theory exhibit better results than both CMC method and ARCHC method do. In addition, the calculated value of the improved DS theory is obviously higher than those of the other three methods for the damaged third storey. Although the other three methods can also approximately find the damage location, the improved DS theory is obviously superior to the other three methods.
Fig. 10. Nonlinear damage detection results in damaged State 2
Fig. 11. Nonlinear damage detection results in damaged State 3
Generally, the greater the gap distance is, the weaker the nonlinear damage degree becomes. Therefore, the damaged state 5 can be considered as a weak nonlinearity problem. For damaged state 5, the damage detection results are shown in Fig. 13. For the weak nonlinearity problem, the results of CMC and ARCHC methods are not good enough, and the nonlinear identification results of DS theory are still superior to those of both CMC method and ARCHC method. The value of improved DS theory is slightly higher than that of the basic DS theory for the damaged third storey. Therefore, the improved DS theory is still the best one among the four methods for the weak nonlinearity damage problem.
From numerical and experimental examples, for the nonlinear damage problem with measured error, the CMC method, ARCHC method, basic DS theory and improved DS theory have the capacity of detecting the damage location. The ARCHC method based on the ARCH model is better than the CMC method based on the cepstral metric and AR model. The detection results of two DS theories are superior to those of both CMC method and ARCHC method. The reason is that the DS theory integrates the damage information of two sources. In addition, the results of improved DS theory are better than those of the basic DS theory. So, it is obvious that the improved DS theory is the best one among the four methods.
Fig. 12. Nonlinear damage detection results in damaged State 4
Fig. 13. Nonlinear damage detection results in damaged State 5 (weak nonlinearity)
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
7. Conclusions
In this study, an improved DS theory based on the time series model has been proposed to raise the detection accuracy of nonlinear damage. First, acceleration responses in the undamaged and damaged states are measured by using accelerometers. Then, acceleration responses are utilized to establish an AR model and residual time series of acceleration responses are used to establish an ARCH model. The CMC method based on the AR model is proposed to obtain local damage solution, and the ARCHC method based on the ARCH model is also proposed to acquire another local damage solution. Finally, the DS theory is applied to damage detection by integrating these local damage solutions from the ARCHC and CMC methods. An improved DS theory is further presented to enhance the detection accuracy. The main advantage of the improved DS theory is that the proposed theory can effectively solve the problem that different information sources have different importance. By taking into account the results of simulation and experiment, the following conclusions can be drawn: (1) For the nonlinear damage problem, the ARCHC method based on the ARCH model is better than the CMC method based on the AR model. (2) The detection results of the improved DS theory are not only superior to those of the basic DS theory, but also better than those of the ARCHC and CMC methods. In a word, the improved DS theory exhibits better performance and higher accuracy than other methods.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant No. 51578094).
References
 Sarrafi A., Mao Z., Niezrecki C., et al. Vibrationbased damage detection in wind turbine blades using phasebased motion estimation and motion magnification. Journal of Sound and Vibration, Vol. 421, Issue 12, 2018, p. 300318. [Publisher]
 Li P. J., Xu D. W., Zhang J. Probabilitybased structural health monitoring through Markov chain Monte Carlo sampling. International Journal of Structural Stability and Dynamics, Vol. 16, Issue 7, 2016, p. 1550039. [Publisher]
 Xu Y., Qian Y., Chen J., et al. Probabilitybased damage detection using model updating with efficient uncertainty propagation. Mechanical Systems and Signal Processing, Vol. 60, Issue 8, 2015, p. 958970. [Publisher]
 Tatar A., Niousha A., Rofooei F. R. Damage detection in existing reinforced concrete building using forced vibration test based on mode shape data. Journal of Civil Structural Health Monitoring, Vol. 7, Issue 1, 2017, p. 123135. [Publisher]
 Guo H. Y., Yuan H. F., Huang Q. Structural damage identification based on gray cloud rule generator algorithm. Advances in Mechanical Engineering, Vol. 11, Issue 1, 2019, https://doi.org/10.1177/1687814018819904. [Publisher]
 Muślewski L., Pająk M., Grządziela A., et al. Analysis of vibration time histories in the time domain for propulsion systems of minesweepers. Journal of Vibroengineering, Vol. 17, Issue 3, 2015, p. 13091316. [CrossRef]
 Grządziela A., Musiał J., Muślewski L., et al. A method for identification of noncoaxiality in engine shaft lines of a selected type of naval ships. Polish Maritime Research, Vol. 22, Issue 1, 2015, p. 6571. [Publisher]
 Das S., Saha P., Patro S. K. Vibrationbased damage detection techniques used for health monitoring of structures: a review. Journal of Civil Structural Health Monitoring, Vol. 6, Issue 3, 2016, p. 477507. [Publisher]
 Eloi F., Joaquim F., Gyuhae P., et al. Influence of the autoregressive model order on damage detection. Computer – Aided Civil and Infrastructure Engineering, Vol. 26, Issue 3, 2011, p. 225238. [Publisher]
 Bernal D., Zonta D., Pozzi M. ARX residuals in damage detection. Structural Control and Health Monitoring, Vol. 19, Issue 4, 2012, p. 535547. [Publisher]
 Bao C., Hao H., Li Z. X. Integrated ARMA model method for damage detection of subsea pipeline system. Engineering Structures, Vol. 48, Issue 3, 2013, p. 176192. [Publisher]
 Mei L., Mita A., Zhou J. An improved substructural damage detection approach of shear structure based on ARMAX model residual. Structural Control and Health Monitoring, Vol. 23, Issue 2, 2016, p. 218236. [Publisher]
 Zheng H., Mita A. Damage indicator defined as the distance between ARMA models for structural health monitoring. Structural Control and Health Monitoring, Vol. 15, Issue 7, 2008, p. 9921005. [Publisher]
 Xing Z., Mita A. Locating the damaged storey of a building using distance measures of loworder AR models. Smart Structures and Systems, Vol. 6, Issue 9, 2010, p. 9911005. [Publisher]
 Liu G., Mao Z., Luo J. Damage detection with interval analysis for uncertainties quantification. International Conference on Performancebased and Lifecycle Structural Engineering, Queensland, Australia, 2015, p. 272281. [Publisher]
 Prawin J., Rao A. R. M. Nonlinear structural damage detection based on adaptive Volterra filter model. International Journal of Structural Stability and Dynamics, Vol. 18, Issue 2, 2018, p. 1871003. [Publisher]
 Ebrahimian H., Astroza R., Conte J. P. Nonlinear finite element model updating for damage identification of civil structures using batch Bayesian estimation. Mechanical Systems and Signal Processing, Vol. 84, Issue 1, 2017, p. 194222. [Publisher]
 Farrar C. R., Worden K., Todd M. D., et al. Nonlinear System Identification for Damage Detection. Technical Report, Los Alamos National Laboratory (LANL), Los Alamos, 2007. [CrossRef]
 Engle R. F. Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica, Vol. 50, Issue 4, 1982, p. 9871007. [Publisher]
 Pham H. T., Yang B. S. Estimation and forecasting of machine health condition using ARMA/GARCH model. Mechanical Systems and Signal Processing, Vol. 24, Issue 2, 2010, p. 546558. [Publisher]
 Chen L. J., Yu L. Structural nonlinear damage identification algorithm based on time series ARMA/GARCH model. Advances in Structural Engineering, Vol. 16, Issue 9, 2013, p. 15971609. [Publisher]
 Yu L., Zhu J. H. Nonlinear damage detection using higher statistical moments of structural responses. Structural Engineering and Mechanics, Vol. 54, Issue 2, 2015, p. 221237. [Publisher]
 Cheng J. J., Guo H. Y., Wang Y. S. Structural nonlinear damage detection method using AR/ARCH model. International Journal of Structural Stability and Dynamics, Vol. 17, Issue 8, 2017, p. 1750083. [Publisher]
 Martin R. J. A metric for ARMA processes. IEEE Transactions on Signal Processing, Vol. 48, Issue 4, 2000, p. 11641170. [Publisher]
 Sun L., Srivastava R. P., Mock T. J. An information systems security risk assessment model under the DempsterShafer theory of belief functions. Journal of Management Information Systems, Vol. 22, Issue 4, 2006, p. 109142. [Publisher]
 Rébillat M., Hajrya R., Mechbal N. Nonlinear structural damage detection based on cascade of Hammerstein models. Mechanical Systems and Signal Processing, Vol. 48, Issue 1, 2014, p. 247259. [Publisher]
 Figueiredo E., Park G., Figueiras J., et al. Structural Health Monitoring Algorithm Comparisons Using Standard Data Sets. Report No.: LA14393, Los Alamos National Laboratory (LANL), Los Alamos, United States, 2009. [Publisher]