Technology Developments in Structural Health Monitoring and Integrity Maintenance
View this Special IssueResearch Article  Open Access
Damage Detection of Structures Identified with DeterministicStochastic Models Using Seismic Data
Abstract
A deterministicstochastic subspace identification method is adopted and experimentally verified in this study to identify the equivalent singleinputmultipleoutput system parameters of the discretetime state equation. The method of damage locating vector (DLV) is then considered for damage detection. A series of shaking table tests using a fivestorey steel frame has been conducted. Both single and multiple damage conditions at various locations have been considered. In the system identification analysis, either full or partial observation conditions have been taken into account. It has been shown that the damaged stories can be identified from global responses of the structure to earthquakes if sufficiently observed. In addition to detecting damage(s) with respect to the intact structure, identification of new or extended damages of the asdamaged counterpart has also been studied. This study gives further insights into the scheme in terms of effectiveness, robustness, and limitation for damage localization of frame systems.
1. Introduction
Structural health monitoring of civil engineering structures has received considerable attention over the last two decades with the progress in digital signal processing and system identification techniques [1–6]. Damage to a structure alters its dynamic characteristics in terms of the modal frequencies, damping ratios, and mode shapes, so as the physical parameters in terms of the stiffness or flexibility matrices. Despite that damage of a structure might be revealed from the changes of frequencies, it is difficult to locate the damages only with this information. The mode shapes perhaps provide a preferable basis for damage detection as they are spatially specific and reflective of local structural behavior. However, information of mode shapes alone is not sensitive enough for damage localization from the global dynamic vibration of the structure. The physical parameters indeed are more useful as far as damage localization is concerned. The study by Zhao and DeWolf [7] comparing various vibrational parameters for damage detection of springmass systems concluded that the model flexibilities were more sensitive to structural damages than either the natural frequencies or mode shapes. Stiffness is intuitively the most direct physical parameter related to damages in structures. However, the sensitivity analysis developed for damage detection based on stiffness matrix requires an accurate analytical model of the intact structure as a reference. Unfortunately, obtaining such an analytical model is in itself a difficult task. Moreover, synthesis of the stiffness matrix requires contributions of the higher modes that are practically difficult to identify. In contrast, the flexibility matrix can be sufficiently synthesized with a limited number of the lower modes as the modal contribution to the flexibility matrix decreases with the square of the corresponding modal frequency. Therefore, flexibilitybased techniques have been considered of great potential in damage localization of structures from the global dynamic responses.
Another branch on the development of SHM techniques based on wavelet transformation has also gained a great deal of progress recently. A review by Li et al. [8] gives an insight of the uptodate innovations on the SHM of infrastructures. In this paper, related theories and methodologies including sensing technology, sensor placement, signal processing, data fusion, system identification, and damage detection have been thoroughly discussed. Li et al. [9] propose a novel wavelet approach integrated with an autoregressive moving average (ARMA) model for the damage detection of structures with progressive damage in time. The methodology has been validated both numerically and experimentally through shaking table tests of a reinforced concrete frame. The concept of wavelet transformation has also been adopted by Li et al. [10] to develop a method for determining the critical incidence of the seismic wave interpreted via energy principles. Yi et al. [11] proposes an energybased multistage structural damage diagnosis method by adopting the schemes of wavelet packet transform (WPT) and artificial neural network (ANN). Using the benchmark structure of the American Society of Civil Engineers (ASCE) as the target, the authors demonstrate via numerical simulations that the proposed method is able to detect structural damage of various extents. Yet, localization of damages requires other measures.
Among the damage localization techniques based on variation of the physical parameters as the structure deteriorates, the flexibilitybased approaches have been shown to be very promising and computationally efficient. Pandey and Biswas [12] demonstrated that the damage locations of a wideflange steel beam could be successfully identified by interrogating the change in the flexibility. This technique has been further extended for damage detection of the plane trusses [13]. The flexibilitybased damage localization technique referred to as the method of damage locating vectors (DLV) has been advanced further by Bernal [14] and Bernal and Gunes [15]. This methodology has also been adopted for damage localization of space trusses or plates by Gao et al. [16, 17] and Huynh et al. [18]. The concept of the DLV method is to identify the members with zero stress under some specific loading patterns, namely, the DLVs that span the null space of the change in flexibility matrix of the structure before and after the damage state. The loading vectors are obtained by performing the singular value decomposition (SVD) on the flexibility differential matrix. Those structural elements resulting with zero stresses (or internal forces) under the static loads of the DLVs are considered potentially damaged. The DLV technique is capable of identifying scenarios of multiple damages in the structure via a truncated modal basis without a predetermined reference model. This null space approach has recently been further extended by Bernal [19, 20] as the dynamic damage locating vector (DDLV) theorem that connects changes in the transfer matrices to the spatial location of stiffness related damage. Bernal [21] revised the DLV method further in the context of linear statespace models to comply with the stochastic realization of the systems by acknowledging that the DLVs in the null space of the flexibility differential can be extracted without a priori knowledge of the flexibility matrix. This scheme is adopted in this study.
System identification problems are commonly concerned with determining, from the relation of input/output sequences, the parameters of either polynomial models (e.g., ARMA, autoregressive model with external input (ARX)) using the prediction error method [1] or statespace models by subspace system identification (SSI) [22], numerical algorithms for subspace state space system identification (N4SID) [23] for deterministic/stochastic systems, or SRIM for deterministic systems [24]. In contrast to the prediction error method that identifies the system parameters by minimizing the errors between the observations and those predicted by the model in a least squares sense, the stochastic realization theory initiated by Akaike [25] and Faure [26] is not based on the concept of optimization and leads to noniterative, convergenceguaranteed and efficient numerical algorithms. The covariance matrix is first constructed from block Hankle matrix [27] of shifted process sequences. The statespace model is in turn realized from the observability and/or controllability matrix via singular value decomposition of the covariance matrix. The theory of covariancedriven subspace method has been unified by van Overschee and de Moor [23] for deterministic, stochastic, and combined systems by defining the estimated state sequences as the projection of inputoutput data. The projected state sequences turn out to be the outputs of nonsteadystate Kalman filter banks. To facilitate implementation of the DLV method, system identification technique of subspace technique or SRIM that identifies the system matrix and output (observation) matrix of a statespace model is considered a preferable alternative. The method of system realization using information matrix (SRIM) proposed by Juang [24] is based on a deterministic statespace system. The equivalent system matrix is identified from the covariance matrix of the input and output signals. This simple and elegant approach works well for systems in response to transient excitation if the noise level is negligible. The performance degrades, however, as the noise becomes pronounced. On the contrary, the stochastic subspace identification (SSI) method based on a stochastic model [23] without knowing the input is less sensitive to noise and works well if the excitation is Gaussian white noise. This scheme is not sufficient, however, for systems under transient excitations such as earthquakes. Therefore, a mixed deterministicstochastic model is considered more robust to transient systems with nonnegligible noise contamination.
In this study, the feasibility of DLV method for damage detection of frame structures using seismic response data is explored in association with the subspace identification algorithm developed by van Overschee and de Moor [23] for mixed deterministicstochastic systems. As an effort to verify the robustness of the deterministicstochastic model against noise, system identification analysis by the N4SID algorithm is conducted with various noise levels and compared with those obtained by the SRIM. The effectiveness of the DLV method for damage detection using system parameters identified from contaminated signals is further investigated. Moreover, a series of shaking table tests using a fivestorey steel frame has been conducted in National Center for Research on Earthquake Engineering (NCREE), Taiwan. Damage condition is simulated by reducing the crosssectional area of some of the columns at the bottom. Both single and combinations of multiple damage conditions at various locations have been considered. In the system identification analysis, either full or partial observation conditions have been taken into account. It has been shown that the damaged stories can be identified from global responses of the structure to earthquakes if sufficiently observed. In addition to detecting damage(s) with respect to the intact structure, identification of new or extended damages of the asdamaged (illconditioned) counterpart has also been studied. The proposed scheme proves to be effective. This study gives further insights into the scheme in terms of effectiveness, robustness, and limitation for damage localization of frame systems.
2. DeterministicStochastic Subspace System Identification
A deterministicstochastic linear timeinvariant system is represented in a discretetime statespace model as with where and are, respectively, the state and output vectors and is the input vector at time instant with , , and being, respectively, the dimension of the dynamic system, observational outputs, and external inputs. is the system matrix, is the input influence matrix, is the observation matrix, and is the direct transmission matrix. and are the unmeasurable vector signals assumed to be zeromean, stationary white noise vector sequences. is the Kronecker .
The state and output vectors are split into deterministic and stochastic parts as and satisfying, respectively, the deterministic subsystem and the stochastic subsystem By recursive substitution into the statespace equations of consecutive shifted processes, it leads to [28] where is the output block Hankle matrix [27] of the past where is the number of steps to be considered in the Markov process and represents the length of the data to be included in the analysis; is the input block Hankle matrix of the past; is the stochastic output block Hankle matrix of the past; is the deterministic state matrix; is the observability matrix; is the controllability matrix; is the triangular Toeplits matrix [27], and
The stochastic covariance equations of the subspace can be defined as in which and represent, respectively, the autocovariance matrix of the state and output vectors and is the covariance matrix of the state and output vectors.
The block Toeplitz covariance matrix of the stochastic output is defined as where and the block Toeplitz cross covariance matrix is defined similarly as where and represent the data length. Equation (18) is the stochastic controllability matrix. Projecting the future output state matrix onto the input matrix and the past output state matrix , the outputinput block Hankel matrices are constructed as where denotes the pseudoinverse of the matrix . To further simplify (19) and (20), the deterministic subspace and stochastic subspace are utilized to define what follows: Assuming that the deterministic input and the deterministic state are independent of the stochastic output , then (19) and (20) can be simplified as where where in which It follows immediately from (22) that where or denotes the pseudoinverse of its counterpart. The sequences and proved to be the states of a bank of nonsteadystate Kalman filters [23]. The nonsteadystate Kalman filter state can be defined in a recursive form as [23] If we collect consecutive columns of the Kalman states in parallel, (27) can then be extended as where . It is trivial that where . Combining (30) and (31), one gets Substituting (26) for and into (32), it leads to in which By forcing the noise terms in (33) to be zero, the coefficient matrices may be resolved as where is the pseudoinverse of . A numerically stable and efficient algorithm referred to as the N4SID devised by van Overschee and de Moor [22] is adopted in this study to solve for the system matrix.
3. N4SID Algorithm
To facilitate implementation of the numerical algorithm, van Overschee and de Moor [22] proposed a numerically stable procedure utilizing the decomposition of the block Hankle matrix where and is a lower triangular matrix. It can be expressed in the following partitioned form as which can be written in a more condensed manner as or with Now if we express and , then (19) can be written as By the same token, (20) may also be written as It has been shown by van Overschee and de Moor [22] that the observability matrix can be obtained from SVD of The rank is determined from the dominant singular values of the decomposition, and can be determined as and can be extracted directly from without the last rows. With and determined, one can calculate and without difficulty and in turn solve (35).
4. Review of the DLV Method
Bernal and Gunes [15] proposed that the structure subjected to the damage locating vectors, , would undergo the same deformation before and after the damaged state. This statement immediately leads to where is the flexibility differential of the structure before and after being damaged. When rank ( is the degree of freedom of the structure), the basis corresponding to the null space of is the damage locating vectors, , which can be derived from singular value decomposition of the flexibility differential of the structure before and after the damage state. Members with nearly zero stress under the loadings of DLVs are considered potentially damaged.
The flexibility matrix of the structure can be expressed with the system matrices of the continuoustime statespace representation as where is the continuoustime system matrix; ; ; ( being the mass matrix of the system).
With (46), the flexibility differential can be expressed as where and since the mass matrix is unchanged. By taking the singular value decomposition of , the eigenvectors corresponding to the singular eigenvalues are the damage locating vector .
The weighted stress index is defined as where is the normalized stress index of the th member or d.o.f. subjected to the th DLV. Member (or storey ) is considered seriously damaged when the normalized weighted stress index in which , whereas it is considered moderately damaged as .
5. Numerical Illustration of the Proposed Schemes
As an effort to verify the robustness of the deterministicstochastic model against noise, system identification analysis by the N4SID algorithm is conducted with various noise levels and compared with those obtained by the SRIM. A fivestory diagonally braced shear frame (Figure 1) with parameters shown in Table 1 is considered in the numerical example. The 1940 El Centro earthquake with peak ground acceleration (PGA) scaled to 0.1 g is considered as the input. To simulate noisy conditions, white noises of various noisetosignal ratios (NSR) are added to the dynamic response signals. The noisetosignal ratio at the th storey, , is defined as where and represent, respectively, the root mean squares (RMS) of the added white noise and the original response signals at the th storey. In the simulation, the noise levels at all floors are assumed to be of the same amount in each case with %, 5% and 10%, respectively.

To assess the accuracy of the identification results, an error index of the th mode is defined as where is the total number of modes, in which and are, respectively, the th element of the th mode shape derived from the eigen analysis and system identification.
Figure 2(a) shows the comparison of mode shapes identified by the N4SID algorithm at various NSR levels with the analytical solutions. When the signals are noisefree (NSR = 0%), all the identified mode shapes coincide with their corresponding analytical counterparts. Deviation of the mode shapes increases with the noise level and becomes more pronounced for higher modes as observed from Figure 2(a). This is also revealed quantitatively from the error index bar charts as illustrated in Figure 3(a). When the noise level achieves , only the first three modes can be identified with acceptable accuracy. The results by the SRIM exhibit similar trends, as shown in Figures 2(b) and 3(b), but deviate in a larger extent from the analytical solutions. When the noise level achieves NSR = 10%, moreover, only the first mode is identified with fidelity. The N4SID algorithm proves more robust to noises than the SRIM, as expected.
(a)
(b)
(a)
(b)
In order to examine if the damage detection scheme works sufficiently when the dynamic response signals are contaminated with noises, numerical simulations are next conducted using the same analytical model. Damage conditions of the structure are simulated by removing some of the diagonal bracings. Both single and multiple damages have been considered, which includes Cases Dl, D2, D3, D4, D5, D13, D15, and D135 where the number(s) denote the damaged stories. Only the N4SID algorithm will be adopted and a noise level of NSR = 10% is considered. Analysis for the noisefree condition has also been carried out as a reference. The results are summarized and compared in Table 2. The shaded area in the table with the normalized weighted stress index indicates those identified as potentially damaged. It is observed that in both cases with NSR = 0% and 10%, the damaged stories can be successfully identified without exception, despite that results of the noisefree condition are shown to be more sensitive in the sense of getting smaller at the damaged locations. The proposed scheme proves effective even if there exists a certain level of noise and only the first three modes are accurately identified.

6. Experimental Verifications
As a further step in verifying the feasibility of the proposed scheme experimentally, a series of shaking table tests has been carried out in the National Center for Research on Earthquake Engineering (NCREE), Taiwan, using a benchmark model (Figure 4). Weighing 22.5 tons, the benchmark model consisting of steel Ibeams is a fivestorey frame of 7.5 m in height and in plane, as illustrated in Figure 5. Accelerometers have been implemented at each floor and the base to monitor the dynamic responses of the model to ground motion which together serve as the basis for system identification. In order to get more insight of the structural behavior, additional velocity meters, LVDT, and load cells have also been implemented in the tests, as illustrated in Figure 6. Only the acceleration information is utilized, however, in this study.
Damage of the structure is simulated by cutting out a small portion of the flange near the bottom of the column(s), as shown in Figure 7. In order to sufficiently examine the damage at various extents, the damages were progressively enforced on one side of the frame from the first storey up to the third storey from one test to another. It is meant to represent a moderate damage condition as a single column is damaged in the same storey while representing a serious damage condition as two columns are damaged. In total, six damage conditions have been considered in the test program. The 1940 El Centro earthquake has been adopted as the input with the PGA scaled to 0.1 g.
The test results are analyzed under considerations of full observation (utilizing acceleration responses of all stories), partial observation (ignoring accelerations at some stories), and the illconditioned condition where the reference structure has been wounded in previous tests.
6.1. Full Observation
The six damaged conditions simulated in the tests are designated as:M1:single column damaged at the first storey, representing a moderate damage condition of the first storey;S1:two columns damaged at the first storey, representing a serious damage condition of the first storey;S1M2:two columns damaged at the first storey and single column damaged at the second storey;S12:two columns damaged at both the first and second stories;S12M3:two columns damaged at both the first and second stories and single column damaged at the third storey;S123:two columns damaged from storey 1 to 3.
The assessment results of various damage conditions based on full observation data of the structure are summarized in Table 3 where the shaded area corresponds to those being screened out as potentially damaged stories. Take case M1 for example, only the normalized weighted stress index of 1F is below 0.1, indicating damage of the first storey. This agrees with the actual damage condition. And in case S1M2, the index is 0.04 for 1F and 0.12 for 2F, agreeing with the actual condition that the first storey is seriously damaged and the second storey is moderately damaged. It is evident that, under a full observation condition, the damaged location(s) are successfully identified, regardless of single or multiple damage conditions.
 
indicates serious damage; indicates moderate damage. *Good indicates the damaged location(s) being identified without missjudgment. Fair indicates the damaged storey being identified but the extent might be underestimated. Poor indicates failing to identify one of the damaged stories. Fail indicates failing to identify more than one of the damaged stories. 
6.2. Partial Observation
The 12 cases analyzed in a partial observation condition are designated below as follows: M1/135: single column damaged at the first storey, representing a moderate damage condition of the first storey with only the first, third, and fifth floors observed; S1/135: two columns damaged at the first storey, representing a serious damage condition of the first storey with only the first, third, and fifth floors observed; M1/124: single column damaged at the first storey with only the first, second, and fourth floors observed; S1/124: two columns damaged at the first storey with only the first, second, and fourth floors observed; S1M2/124: two columns damaged at the first storey and single column damaged at the second storey with only the first, second, and fourth floors observed; S1M2/125: two columns damaged at the first storey and single column damaged at the second storey with only the first, second, and fifth floors observed; S1M2/1235: two columns damaged at the first storey and single column damaged at the second storey without observing the fourth floor; S12/124: two columns damaged at both the first and second stories with only the first, second, and fourth floors observed; S12/125: two columns damaged at both the first and second stories with only the first, second, and fifth floors observed; S12/1235: two columns damaged at both the first and second stories without observing the fourth floor; S12M3/1235: two columns damaged at both the first and second stories and single column damaged at the third storey without observing the fourth floor; S123/1235: two columns damaged at the first, second and third stories without observing the fourth floor.
The assessment results of various damage conditions based on partial observation data of the structure are summarized in Table 4 where the shaded area corresponds to those being screened out as potentially damaged stories. The performance index (PMC) is similarly defined as in Table 3. The results show that the scheme with partial observation remains effective for single damage conditions, as in the cases of M1/135, S1/135, M1/124 and S1/124 where the index corresponding to the first storey is below 0.2. The scheme fails, however, to locate the damaged stories in multiple damage conditions except for case S12/1235 where the first two stories are seriously damaged and 4 out of 5 stories are observed.

6.3. IllConditioned Structures
At times the structural health monitoring system might be introduced after the target building has been previously damaged. It is of interest to verify if the scheme is able to identify new or extended damage(s) in an earthquake event of an illconditioned structure that has been previously damaged. The system identification analysis will be based on full observation data as it provides more reliable structural information for damage assessment. The 10 cases considered for damage detection analysis of illconditioned structures are designated below as follows: S1/M1: structure seriously damaged at the first storey of the current state versus the one with its first storey moderately damaged earlier; S1M2/S1: structure damaged seriously at the first storey and moderately at the second storey of the current state versus the one with its first storey seriously damaged earlier; S12/S1M2: structure damaged seriously at both the first and second stories of the current state versus the one seriously damaged at its first storey and moderately at the second in advance; S12M3/S12: structure damaged seriously at both the first and second stories and moderately at the third storey of the current state versus the one with its both first and second stories seriously damaged earlier; S123/S12M3: structure seriously damaged at the first, second, and third stories of the current state versus the one seriously damaged at its both first and second stories and moderately at the third in advance; S12/M1: structure seriously damaged at both the first and second stories of the current state versus the one with its first storey moderately damaged earlier; S12M3/M1: structure seriously damaged at both the first and second stories and moderately at the third of the current state versus the one with its first storey moderately damaged earlier; S123/M1: structure seriously damaged at the first, second, and third stories of the current state versus the one with its first storey moderately damaged earlier; S12/S1: structure seriously damaged at both the first and second stories of the current state versus the one with its first storey seriously damaged earlier; S123/S1: structure seriously damaged at the first, second, and third stories of the current state versus the one with its first storey seriously damaged earlier.
The assessment results of various damage conditions in respect to an illconditioned structure damaged earlier are summarized in Table 5. The performance index (PMC) is similarly defined as in Table 3. In all the cases considered, the scheme proves sufficient in identifying new or extended damage(s) without exception under a full observation condition in the system identification process, as indicated from the indices.

7. Summary and Conclusion
A scheme integrated with deterministicstochastic subspace system identification and the method of damage localization vector for damage detection of structures is considered and verified based on seismic response data in this study. Both numerical simulation and experimental verification have been carried out to verify the feasibility of the proposed scheme. A series of shaking table tests using a fivestorey steel frame has been conducted in National Center for Research on Earthquake Engineering (NCREE), Taiwan. Both single and combinations of multiple damage conditions at various locations have been considered. Either a full or partial observation conditions have been taken into account in the system identification process. This study gives further insights into the scheme in terms of effectiveness, robustness, and limitation for damage localization of frame systems. Based on the numerical simulation and test results, the conclusions are drawn as the following.(a)The N4SID algorithm for deterministicstochastic models is more robust to noise contamination than the SRIM algorithm for purely deterministic models. With the first three modes of the sheartype frame system accurately identified from noisecontaminated signals with 10% NSR by using the N4SID and the flexibilitybased DLV method, all damaged conditions considered can be successfully predicted.(b)Damage localization utilizing seismic response data, in particular the floor accelerations, proves feasible not only analytically but also experimentally.(c)Under a full observation condition where all floor responses are observed, the damaged location(s) can be successfully identified, regardless of single or multiple damage conditions.(d)Under a partial observation condition where 3 out of 5 floor responses are observed, only those with single damage can be identified if the damaged storey is colocated with one of the observed floors. The scheme fails, however, to fully locate the damages in multiple damage conditions in general.(e)The scheme proves to be sufficient in identifying new or extended damage(s) without exception under a full observation condition.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work is supported by the National Science Council of Republic of China under contracts NSC 1022221E009086 and NSC1022625M344001.
References
 L. Ljung, System Identification: Theory for the Users, Prentice Hall PTR, Upper Saddle River, NJ, USA, 1987.
 E. Safak, “Adaptive modeling, identification, and control of dynamic structural systems. I: Theory,” Journal of Engineering Mechanics, ASCE, vol. 115, no. 11, pp. 2386–2405, 1989. View at: Publisher Site  Google Scholar
 E. Safak, “Adaptive modeling, identification, and control of dynamic structural systems. II: applications,” Journal of Engineering Mechanics, vol. 115, no. 11, pp. 2406–2426, 1989. View at: Publisher Site  Google Scholar
 J. Lin and R. Betti, “Online identification and damage detection in nonlinear structural systems using a variable forgetting factor approach,” Earthquake Engineering and Structural Dynamics, vol. 33, no. 4, pp. 419–444, 2004. View at: Publisher Site  Google Scholar
 M. J. Perry and C. G. Koh, “Outputonly structural identification in time domain: numerical and experimental studies,” Earthquake Engineering and Structural Dynamics, vol. 37, no. 4, pp. 517–533, 2008. View at: Publisher Site  Google Scholar
 “A review of structural health monitoring literature, 1996–2001,” Los Alamos National Laboratory Report LA13976MS, 2003. View at: Google Scholar
 J. Zhao and J. T. DeWolf, “Sensitivity study for vibrational parameters used in damage detection,” Journal of Structural Engineering, vol. 125, no. 4, pp. 410–416, 1999. View at: Publisher Site  Google Scholar
 H. N. Li, T. H. Yi, L. Ren, D. S. Li, and L. S. Huo, “Reviews on innovations and applications in structural health monitoring for infrastructures,” Structural Monitoring and Maintenance, vol. 1, no. 1, pp. 1–45, 2014. View at: Google Scholar
 H.N. Li, T.H. Yi, M. Gu, and L.S. Huo, “Evaluation of earthquakeinduced structural damages by wavelet transform,” Progress in Natural Science, vol. 19, no. 4, pp. 461–470, 2009. View at: Publisher Site  Google Scholar
 H. N. Li, X. Y. He, and T. H. Yi, “Multicomponent seismic response analysis of offshore platform by wavelet energy principle,” Coastal Engineering, vol. 56, no. 8, pp. 810–830, 2009. View at: Publisher Site  Google Scholar
 T. H. Yi, H. N. Li, and H. M. Sun, “Multistage structural damage diagnosis method based on “energydamage” theory,” Smart Structures and Systems, vol. 12, no. 34, pp. 345–361, 2013. View at: Google Scholar
 A. K. Pandey and M. Biswas, “Damage detection in structures using changes in flexibility,” Journal of Sound and Vibration, vol. 169, no. 1, pp. 3–17, 1994. View at: Publisher Site  Google Scholar
 A. K. Pandey and M. Biswas, “Damage diagnosis of truss structures by estimation of flexibility change,” The International Journal of Analytical and Experimental Modal Analysis, vol. 10, no. 2, pp. 104–117, 1995. View at: Google Scholar
 D. Bernal, “Load vectors for damage localization,” Journal of Engineering Mechanics, vol. 128, no. 1, pp. 7–14, 2002. View at: Publisher Site  Google Scholar
 D. Bernal and B. Gunes, “Flexibility based approach for damage characterization: benchmark application,” Journal of Engineering Mechanics, vol. 130, no. 1, pp. 61–70, 2004. View at: Publisher Site  Google Scholar
 Y. Gao, M. RuizSandoval, and B. F. Spencer, “Flexibilitybased damage localization employing ambient vibration,” in Proceedings of the 15th ASCE Engineering Mechanics Conference, Columbia University, New York, NY, USA, June 2002. View at: Google Scholar
 Y. Gao, B. F. Spencer, and D. Bernal, “Experimental verification of the flexibilitybased damage locating vector method,” Journal of Engineering Mechanics, vol. 133, no. 10, pp. 1043–1049, 2007. View at: Publisher Site  Google Scholar
 D. Huynh, J. He, and D. Tran, “Damage location vector: a nondestructive structural damage detection technique,” Computers & Structures, vol. 83, no. 28–30, pp. 2353–2367, 2005. View at: Publisher Site  Google Scholar
 D. Bernal, “Damage localization from the null space of changes in the transfer matrix,” AIAA Journal, vol. 45, no. 2, pp. 374–381, 2007. View at: Publisher Site  Google Scholar
 D. Bernal, “Damage localization from transmission zeros of delta systems,” Journal of Engineering Mechanics, vol. 135, no. 2, pp. 93–99, 2009. View at: Publisher Site  Google Scholar
 D. Bernal, “Flexibilitybased damage localization from stochastic realization results,” Journal of Engineering Mechanics, vol. 132, no. 6, pp. 651–658, 2006. View at: Publisher Site  Google Scholar
 P. van Overschee and B. de Moor, “N4SID: subspace algorithms for the identification of combined deterministicstochastic systems,” Automatica, vol. 30, no. 1, pp. 75–93, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 P. van Overschee and B. de Moor, Subspace Identification for Linear Systems: Theory Implementation Applications, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1996. View at: MathSciNet
 J. N. Juang, “System realization using information matrix,” Journal of Guidance, Control, and Dynamics, vol. 20, no. 3, pp. 492–500, 1997. View at: Publisher Site  Google Scholar
 H. Akaike, “Markovian representation of stochastic processes by canonical variables,” SIAM Journal on Control and Optimization, vol. 13, no. 1, pp. 162–173, 1975. View at: Publisher Site  Google Scholar  MathSciNet
 P. Faure, “Stochastic realization algorithms,” in System Identification: Advances and Case Studies, R. Mehra and D. Lainiotis, Eds., Academic Press, 1976. View at: Google Scholar
 The Matrix Cookbook, Petersen & Pedersen, 2006.
 Y. T. Lin, Experimental verification of damage localization of outputonly structural systems using seismic response data [M.S. thesis], Department of Civil Engineering, National ChiaoTung University, Hsinchu, Taiwan, 2012.
Copyright
Copyright © 2014 MingChih Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.