Substructural first‐ and second‐order model identification for structural damage assessment

This paper presents two methods to perform system identification at the substructural level, taking advantage of reduction in the number of unknowns and degrees of freedom (DOFs) involved, for damage assessment of fairly large structures. The first method is based on first‐order state space formulation of the substructure where the eigensystem realization algorithm (ERA) and the observer/Kalman filter identification (OKID) are used. Identification at the global level is then performed to obtain the second‐order model parameters. In the second method, identification is performed at the substructural level in both the first‐ and second‐order model identification. Both methods are illustrated using numerical simulation studies where results indicate their significantly better performance than identification using the global structure, in terms of efficiency and accuracy. A 12‐DOF system and a fairly large structural system with 50 DOFs are used where the effects of noisy data are considered. In addition to numerical simulation studies, laboratory experiments involving an eight‐storey frame model are carried out to illustrate the performance of the proposed method. The identification results presented in terms of the stiffness integrity index show that the proposed methodology is able to locate and quantify damage fairly accurately. Copyright © 2005 John Wiley & Sons, Ltd.     

随机子空间识别方法计算效率的改进

参数识别是结构健康监测领域研究中的重点。随机子空间法是近年来发展起来的一种线性系统辩识方法,可以有效地从环境激励的结构响应中获取模态参数。该方法的基本原理是将“将来”数据向“过去”数据进行垂直投影,进而根据该投影计算出可观测矩阵和一个Kalman滤波状态序列。而Kalman滤波序列是“过去”输出信号的线性组合,即“过去”输出和“将来”状态估计建立了关系。而从Hankel矩阵的组成来看,由于要使得该矩阵具有单无限的条件,所需的计算时间也比较长。据此对随机子空间方法进行了改进。改进的基本思想是采用一个测点的信号作为“过去”作为输出信号代替全部测点的信号。这样就减少了计算量。最后用一数值模拟算例进行了验证,结果良好。     

Damage assessment through structural identification of a three‐story large‐scale precast concrete structure

This paper investigates the damage assessment of a three‐story half‐scale precast concrete building resembling a parking garage through structural identification. The structure was tested under earthquake‐type loading on the NEES large high‐performance outdoor shake table at the University of California San Diego in 2008. The tests provide a unique opportunity to capture the dynamic performance of precast concrete structures built under realistic boundary conditions. The effective modal parameters of the structure at different damage states have been identified from white‐noise and scaled earthquake test data with the assumption that the structure responded in a quasi‐linear manner. Modal identification has been performed using the deterministic‐stochastic subspace identification method based on the measured input–output data. The changes in the identified modal parameters are correlated to the observed damage. In general, the natural frequencies decrease, and the damping ratios increase as the structure is exposed to larger base excitations, indicating loss of stiffness, development/propagation of cracks, and failure in joint connections. The analysis of the modal rotations and curvatures allowed the localization of shear and flexural damages respectively and the checking of the effectiveness of repair actions. Copyright © 2013 John Wiley & Sons, Ltd.     

Stochastic subspace identification for output‐only modal analysis: application to super high‐rise tower under abnormal loading condition

The objective of this paper is to develop an online system parameter estimation technique from the response measurements through using the recursive covariance‐driven stochastic subspace identification (SSI‐COV) approach. In developing the recursive SSI‐COV, to avoid time‐consumption of singular value decomposition in recursive SSI, the extended instrumental variable version of the projection approximation subspace tracking method is used in SSI‐COV. Besides, to reduce the effect of noise on the results of identification, the preprocessing of data using recursive singular spectrum analysis technique is also presented to remove the noise contaminant measurements to enhance the stability of data analysis. On the basis of the proposed method, both the ambient vibration and seismic response data of a tower (Canton Tower) are used to observe the time‐varying system natural frequencies of a tower from its operating condition. Results from using off‐line SSI‐COV method under normal operating condition are also presented. Comparison on the identified time‐varying dynamic characteristics of the tower under normal operating condition and earthquake response of distanced earthquake event is discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

Parameter identification of framed structures using an improved finite element model‐updating method—Part I: formulation and verification

In this study, we formulate an improved finite element model‐updating method to address the numerical difficulties associated with ill conditioning and rank deficiency. These complications are frequently encountered model‐updating problems, and occur when the identification of a larger number of physical parameters is attempted than that warranted by the information content of the experimental data. Based on the standard bounded variables least‐squares (BVLS) method, which incorporates the usual upper/lower‐bound constraints, the proposed method (henceforth referred to as BVLSrc) is equipped with novel sensitivity‐based relative constraints. The relative constraints are automatically constructed using the correlation coefficients between the sensitivity vectors of updating parameters. The veracity and effectiveness of BVLSrc is investigated through the simulated, yet realistic, forced‐vibration testing of a simple framed structure using its frequency response function as input data. By comparing the results of BVLSrc with those obtained via (the competing) pure BVLS and regularization methods, we show that BVLSrc and regularization methods yield approximate solutions with similar and sufficiently high accuracy, while pure BVLS method yields physically inadmissible solutions. We further demonstrate that BVLSrc is computationally more efficient, because, unlike regularization methods, it does not require the laborious a priori calculations to determine an optimal penalty parameter, and its results are far less sensitive to the initial estimates of the updating parameters. Copyright © 2006 John Wiley & Sons, Ltd.     

随机子空间方法在桥塔模态参数识别中的应用

基于环境振动的结构模态参数识别方法正逐渐成为国内外研究的一大热点。环境振动方法就是仅仅利用结构测试的输出信号进行结构的模态参数识别,随机子空间方法就是其中的一种。随机子空间法是近年来发展起来的一种线性系统辩识方法,可以有效地从环境激励的结构响应中获取模态参数。它属于时域的方法,该方法不需要进行FFT变换,它不仅可以识别结构的频率,而且可以识别结构的阻尼和振型。文章首先介绍了随机子空间的理论,然后用该方法对正在施工中的南京长江三桥的南塔进行模态参数识别,通过与其他方法的识别结果进行比较,证明随机子空间方法不失为一种有效的模态参数识别方法。     

Output‐only structural identification in time domain: Numerical and experimental studies

By identifying changes in stiffness parameters, structural damage can be detected and monitored. Although considerable progress has been made in this research area, many challenges remain in achieving robust structural identification based on incomplete and noisy measurement signals. The identification task is made even more difficult if measurement of input force is to be eliminated. To this end, an output‐only structural identification strategy is proposed to identify unknown stiffness and damping parameters. A non‐classical approach based on genetic algorithms (GAs) is adopted. The proposed strategy makes use of the recently developed GA‐based method of search space reduction, which has shown to be able to accurately and reliably identify structural parameters from measured input and output signals. By modifying the numerical integration scheme, input can be computed as the parameter identification task is in progress, thereby eliminating the need to measure forces. Numerical and experimental results demonstrate the power of the strategy in accurate and efficient identification of structural parameters and damage using only incomplete acceleration measurements. Copyright © 2007 John Wiley & Sons, Ltd.     

A reassessment of dynamic characteristics of the Quincy Bayview Bridge using output‐only identification techniques

A reassessment of the dynamic characteristics of the 542 m cable‐stayed Bayview Bridge in Quincy, Illinois, is presented using a newly developed output‐only system identification technique. The technique is applied to an extensive set of ambient vibration response data acquired from the bridge in 1987. Vertical, torsional and transverse modal frequencies of the deck are identified, and uncertainty in damping values are estimated using an automated procedure on several redundant measurements at four locations. Important practical implementation issues associated with the implementation of the procedure and selection of algorithm design parameters for stochastic subspace identification techniques are discussed. An overall mean and standard deviation of damping of 1.0±0.8% is estimated considering all identified vertical, torsional and transverse modes in the 0–2 Hz band. The mean damping for the fundamental vertical mode (0.37 Hz) is identified as 1.4±0.5%, and for the first coupled torsion–transverse mode (0.56 Hz) is identified as 1.1±0.8%. Variability in the damping estimates is shown to decrease as estimated modal RMS acceleration levels increase. Standard deviations on estimated damping range from 0.05% to 2%. The results are shown to be a substantial improvement in the evaluation of damping compared to earlier spectral analysis conducted on the same data set. Copyright © 2005 John Wiley & Sons, Ltd.     

Probabilistic approach for modal identification using non‐stationary noisy response measurements only

This paper addresses the problem of identification of the modal parameters for a structural system using measured non‐stationary response time histories only. A Bayesian time‐domain approach is presented which is based on an approximation of the probability distribution of the response to a non‐stationary stochastic excitation. It allows one to obtain not only the most probable values of the updated modal parameters and stochastic excitation parameters but also their associated uncertainties using only one set of response data. It is found that the updated probability distribution can be well approximated by a Gaussian distribution centred at the most probable values of the parameters. Examples using simulated data are presented to illustrate the proposed method. Copyright © 2002 John Wiley & Sons, Ltd.     

Discrete wavelet‐based trend identification in hydrologic time series

Trend identification is a substantial issue in hydrologic series analysis, but it is also a difficult task in practice due to the confusing concept of trend and disadvantages of methods. In this article, an improved definition of trend was given as follows: ‘a trend is the deterministic component in the analysed data and corresponds to the biggest temporal scale on the condition of giving the concerned temporal scale’. It emphasizes the intrinsic and deterministic properties of trend, can clearly distinguish trend from periodicities and points out the prerequisite of the concerned temporal scale only by giving which the trend has its specific meaning. Correspondingly, the discrete wavelet‐based method for trend identification was improved. Differing from those methods used presently, the improved method is to identify trend by comparing the energy difference between hydrologic data and noise, and it can simultaneously separate periodicities and noise. Furthermore, the improved method can quantitatively estimate the statistical significance of the identified trend by using proper confidence interval. Analyses of both synthetic and observed series indicated the identical power of the improved method as the Mann–Kendall test in assessing the statistical significance of the trend in hydrologic data, and by using the former, the identified trend can adaptively reflect the nonlinear and nonstationary variability of hydrologic data. Besides, the results also showed the influences of three key factors (wavelet choice, decomposition level choice and noise content) on discrete wavelet‐based trend identification; hence, they should be carefully considered in practice. Copyright © 2012 John Wiley & Sons, Ltd.     

Blind modal identification of output‐only structures in time‐domain based on complexity pursuit

Output‐only modal identification is needed when only structural responses are available. As a powerful unsupervised learning algorithm, blind source separation (BSS) technique is able to recover the hidden sources and the unknown mixing process using only the observed mixtures. This paper proposes a new time‐domain output‐only modal identification method based on a novel BSS learning algorithm, complexity pursuit (CP). The proposed concept—independent ‘physical systems’ living on the modal coordinates—connects the targeted constituent sources (and their mixing process) targeted by the CP learning rule and the modal responses (and the mode matrix), which can then be directly extracted by the CP algorithm from the measured free or ambient system responses. Numerical simulation results show that the CP method realizes accurate and robust modal identification even in the closely spaced mode and the highly damped mode cases subject to non‐stationary ambient excitation and provides excellent approximation to the non‐diagonalizable highly damped (complex) modes. Experimental and real‐world seismic‐excited structure examples are also presented to demonstrate its capability of blindly extracting modal information from system responses. The proposed CP is shown to yield clear physical interpretation in modal identification; it is computational efficient, user‐friendly, and automatic, requiring little expertise interactions for implementations. Copyright © 2013 John Wiley & Sons, Ltd.     

1D System identification of a 54‐story steel frame building by seismic interferometry

A 54‐story steel, perimeter‐frame building in downtown Los Angeles, California, is identified by a wave method using records of the Northridge earthquake of 1994 (M_L = 6.4, R = 32 km). The building is represented as a layered shear beam and a torsional shaft, characterized by the corresponding velocities of vertically propagating waves through the structure. The previously introduced waveform inversion algorithm is applied, which fits in the least squares sense pulses in low‐pass filtered impulse response functions computed at different stories. This paper demonstrates that layered shear beam and torsional shaft models are valid for this building, within bands that include the first five modes of vibration for each of the North–South (NS), East–West (EW), and torsional responses (0–1.7 Hz for NS and EW, and 0–3.5 Hz for the torsional response). The observed pulse travel time from ground floor to penthouse level is τ ≈1.5 s for NS and EW and τ ≈ 0.9 s for the torsional responses. The identified equivalent uniform shear beam wave velocities are β_eq ≈ 140 m/s for NS and EW responses, and 260 m/s for torsion, and the apparent Q ≈ 25 for the NS and torsional, and ≈14 for the EW response. Across the layers, the wave velocity varied 90–170 m/s for the NS, 80–180 m/s for the EW, and 170–350 m/s for the torsional responses. The identification method is intended for use in structural health monitoring. Copyright © 2013 John Wiley & Sons, Ltd.