Tracking modal parameters from building seismic response data using recursive subspace identification algorithm

Tracking modal parameters and estimating the current structural state of a building from seismic response measurements, particularly during strong earthquake excitations, can provide useful information for building safety assessment and the adaptive control of a structure. Therefore, online or recursive identification techniques need to be developed and implemented for building seismic response monitoring. This paper develops and examines different methods to track modal parameters from building seismic response data. The methods include recursive data‐driven subspace identification (RSI‐DATA) using Givens rotation algorithm, and RSI‐DATA using Bona fide algorithm. The question on how well the results of RSI‐DATA reflect the real condition is investigated and verified with a bilinear SDOF simulation study. Time‐varying modal parameters of a four‐story reinforced concrete school building are identified based on a series of earthquake excitations, including several seismic events, large and small. Discussions on the different methods' ability to track the time‐varying modal parameters are presented. The variation of the identified building modal frequencies and damping ratios from a series of event‐by‐event seismic responses, particularly before and after retrofitting of the building is also discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Identification of base‐excited structures using output‐only parameter estimation

This paper presents a new identification technique for the extraction of modal parameters of structural systems subjected to base excitation. The technique uses output‐only measurements of the structural response. A combined subspace‐maximum likelihood algorithm is developed and applied to a three‐degree‐of‐freedom simulation model. Five ensembles of synthetically generated input signals, representing varying input characteristics, are employed in Monte Carlo simulations to illustrate the applicability of the method. The technique is able to circumvent some of the difficulties arising from short data sets by employing the Expectation Maximization (EM) algorithm to refine the subspace state estimates. This approach is motivated by successful application by previous authors on speech signals. Results indicate that, for certain system characteristics, more accurate pole estimates can be identified using the combined subspace‐EM formulation. In general, the damping ratios of the system are difficult to identify accurately due to limitations on data set length. The applicability of the technique to structural vibration signals is illustrated through the identification of seismic response data from the Vincent Thomas Bridge. Copyright © 2003 John Wiley & Sons, Ltd.     

Two‐stage damage detection algorithms of structure using modal parameters identified from recursive subspace identification

Structural damage assessment under external loading, such as earthquake excitation, is an important issue in structural safety evaluation. In this regard, appropriate data analysis and feature extraction techniques are required to interpret the measured data and to identify the state of the structure and, if possible, to detect the damage. In this study, the recursive subspace identification with Bona‐fide LQ renewing algorithm (RSI‐BonaFide‐Oblique) incorporated with moving window technique is utilized to identify modal parameters such as natural frequencies, damping ratios, and mode shapes at each instant of time during the strong earthquake excitation. From which the least square stiffness method (LSSM) combined with the model updating technique, called efficient model correction method (EMCM), is used to estimate the first‐stage system stiffness matrix using the simplified model from the previously identified modal parameters (nominal model). In the second stage, 2 different damage assessment algorithms related to the nominal system stiffness matrix were derived. First, the model updating technique, called EMCM, is applied to correct the nominal model by the newly identified modal parameters during the strong motion. Second, the element damage index can be calculated using element damage index method (EDIM) to quantify the damage extent in each element. Verification of the proposed methods through the shaking table test data of 2 different types of structures and a building earthquake response data is demonstrated to specify its corresponding damage location, the time of occurrence during the excitation, and the percentage of stiffness reduction.     

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.     

Subspace system identification of support‐excited structures—part I: theory and black‐box system identification

This paper reviews the theoretical principles of subspace system identification as applied to the problem of estimating black‐box state‐space models of support‐excited structures (e.g., structures exposed to earthquakes). The work distinguishes itself from past studies by providing readers with a powerful geometric interpretation of subspace operations that relates directly to theoretical structural dynamics. To validate the performance of subspace system identification, a series of experiments are conducted on a multistory steel frame structure exposed to moderate seismic ground motions; structural response data is used off‐line to estimate black‐box state‐space models. Ground motions and structural response measurements are used by the subspace system identification method to derive a complete input–output state‐space model of the steel frame system. The modal parameters of the structure are extracted from the estimated input–output state‐space model. With the use of only structural response data, output‐only state‐space models of the system are also estimated by subspace system identification. The paper concludes with a comparison study of the modal parameters extracted from the input–output and output‐only state‐space models in order to quantify the uncertainties present in modal parameters extracted from output‐only models. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of subspace identification technique to long‐term seismic response monitoring of structures

This paper addresses the issue of structural system identification using earthquake‐induced structural response. The proposed methodology is based on the subspace identification algorithm to perform identification of structural dynamic characteristics using input–output seismic response data. Incorporated with subspace identification algorithm, a scheme to remove spurious modes is also used to identify real system poles. The efficiency of the proposed method is shown by the analysis of all measurement data from all measurement directly. The recorded seismic response data of three structures (one 7‐story RC building, one midisolation building, and one isolated bridge), under Taiwan Strong Motion Instrumentation Program, are analyzed during the past 15 years. The results present the variation of the identified fundamental modal frequencies and damping ratios from all the recorded seismic events that these three structures had encountered during their service life. Seismic assessment of the structures from the identified system dynamic characteristics during the period of their service is discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

基于经验模式分解的随机子空间识别方法

提出了基于经验模式分解（EMD）的环境激励结构模态参数随机子空间识别（SSI）方法。该方法用设置间断频率的EMD将结构环境振动响应原始信号分解成若干个基本模式分量（IMF）,使每一个基本模式分量仅为结构的某一阶固有模态,进而用随机子空间方法进行模态参数识别。实桥环境振动实验分析结果表明,该方法能有效地避免结构各阶模态之间的相互影响,能够更清晰方便地得到结构的模态参数。     

Response‐only modal identification of structures using strong motion data

Dynamic characteristics of structures — viz. natural frequencies, damping ratios, and mode shapes — are central to earthquake‐resistant design. These values identified from field measurements are useful for model validation and health‐monitoring. Most system identification methods require input excitations motions to be measured and the structural response; however, the true input motions are seldom recordable. For example, when soil–structure interaction effects are non‐negligible, neither the free‐field motions nor the recorded responses of the foundations may be assumed as ‘input’. Even in the absence of soil–structure interaction, in many instances, the foundation responses are not recorded (or are recorded with a low signal‐to‐noise ratio). Unfortunately, existing output‐only methods are limited to free vibration data, or weak stationary ambient excitations. However, it is well‐known that the dynamic characteristics of most civil structures are amplitude‐dependent; thus, parameters identified from low‐amplitude responses do not match well with those from strong excitations, which arguably are more pertinent to seismic design. In this study, we present a new identification method through which a structure's dynamic characteristics can be extracted using only seismic response (output) signals. In this method, first, the response signals’ spatial time‐frequency distributions are used for blindly identifying the classical mode shapes and the modal coordinate signals. Second, cross‐relations among the modal coordinates are employed to determine the system's natural frequencies and damping ratios on the premise of linear behavior for the system. We use simulated (but realistic) data to verify the method, and also apply it to a real‐life data set to demonstrate its utility. Copyright © 2012 John Wiley & Sons, Ltd.     

Discussion of user-defined parameters for recursive subspace identification: Application to seismic response of building structures

Structural damage assessment under external loading, such as earthquake excitation, is an important issue in structural safety evaluation. In this regard, an appropriate data analysis and system identification technique is required to interpret the measured data and to identify the state of the structure. Generally, the recursive system identification algorithm is used. In this study, the recursive subspace identification (RSI) algorithm based on the matrix inversion lemma algorithm with oblique projection technique (RSI-Inversion-Oblique) is applied to investigate the time-varying dynamic characteristics. The user-defined parameters used in the RSI-Inversion-Oblique technique are carefully discussed, which include the size of the data Hankel matrix (i), model order to extract the physical modes, and forgetting factor (FF) to detect the time-varying system modal frequencies. Response data from the Northridge earthquake from the Sherman Oaks building (CSMIP) is used as an example to examine a systematic method to determine the suitable user-defined parameters in RSI. It is concluded that the number of rows in the data Hankel matrix significantly influences the identification of the time-varying fundamental modal frequency of the structure. An algorithmic model order selection method using the eigenvalue distribution of RSI-Inversion can detect the system modal frequencies at each appending data window without causing any abnormality.     

Vibration‐based damage assessment of steel structure using global and local response measurements

Damage assessment of a structure involves acquiring and identifying dynamic characteristics of the structure and using these characteristics to evaluate behavior and performance. In this study, an unsymmetrical three‐story steel structure (fabricated with one weak column in the first floor) was tested on shaking table and subjected to a series of earthquake excitations with increasing level of excitation back to back. Besides, white noise excitation was also applied in between the earthquake excitation to serve as the reference state. Both the traditional sensing system (accelerometer and linear variable differential transformer) and the local optical tracker system were implemented in the structure to collect the vibration‐based responses. For operational modal analysis, structural response from white noise excitation will be used in this study. First, the traditional system identification using global response data is used (multivariate autoregressive (AR)‐model) to extract system natural frequencies and mode shapes from all different set of white noise responses after earthquake excitation. The migration of AR‐coefficient ellipse error from each sensor response was used to identify the damage location. Second, blind source separation technique was used to identify the modal contribution of the structure from each test, which provide information to detect the damage severity. Finally, from the local optical tracker array data, the principal component analysis was applied to quantify the earthquake‐induce local stress of the structural member. Combine the result from damage detection using global measurement and the identified local element stress, one can locate and quantify the damage. Copyright © 2015 John Wiley & Sons, Ltd.     

Seismic response estimation method for earthquake‐damaged RC buildings

This study proposes a procedure for identifying spectral response curves for earthquake‐damaged areas in developing countries without seismic records. An earthquake‐damaged reinforced concrete building located in Padang, Indonesia was selected to illustrate the identification of the maximum seismic response during the 2009 West Sumatra earthquake. This paper summarizes the damage incurred by the building; the majority of the damage was observed in the third story in the span direction. The damage was quantitatively evaluated using the damage index R according to the Japanese guidelines for post‐earthquake damage evaluation. The damage index was also applied to the proposed spectral response identification method. The seismic performance of the building was evaluated by a nonlinear static analysis. The analytical results reproduced a drift concentration in the third story. The R‐index decreased with an increase in the story drift, which provided an estimation of the maximum response of the building during the earthquake. The estimation was verified via an earthquake response analysis of the building using ground acceleration data, which were simulated based on acceleration records of engineering bedrock that considered site amplification. The maximum response estimated by the R‐index was consistent with the maximum response obtained from the earthquake response analysis. Therefore, the proposed method enables the construction of spectral response curves by integrating the identification results for the maximum responses in a number of earthquake‐damaged buildings despite a lack of seismic records. Copyright © 2016 The Authors. Earthquake Engineering & Structural Dynamics published by John Wiley & Sons Ltd.     

System identification applied to long‐span cable‐supported bridges using seismic records

This paper presents the application of system identification (SI) to long‐span cable‐supported bridges using seismic records. The SI method is based on the System Realization using Information Matrix (SRIM) that utilizes correlations between base motions and bridge accelerations to identify coefficient matrices of a state‐space model. Numerical simulations using a benchmark cable‐stayed bridge demonstrate the advantages of this method in dealing with multiple‐input multiple‐output (MIMO) data from relatively short seismic records. Important issues related to the effects of sensor arrangement, measurement noise, input inclusion, and the types of input with respect to identification results are also investigated. The method is applied to identify modal parameters of the Yokohama Bay Bridge, Rainbow Bridge, and Tsurumi Fairway Bridge using the records from the 2004 Chuetsu‐Niigata earthquake. Comparison of modal parameters with the results of ambient vibration tests, forced vibration tests, and analytical models are presented together with discussions regarding the effects of earthquake excitation amplitude on global and local structural modes. Copyright © 2007 John Wiley & Sons, Ltd.     

On‐line identification of non‐linear hysteretic structural systems using a variable trace approach

In this paper, an adaptive on‐line parametric identification algorithm based on the variable trace approach is presented for the identification of non‐linear hysteretic structures. At each time step, this recursive least‐square‐based algorithm upgrades the diagonal elements of the adaptation gain matrix by comparing the values of estimated parameters between two consecutive time steps. Such an approach will enforce a smooth convergence of the parameter values, a fast tracking of the parameter changes and will remain adaptive as time progresses. The effectiveness and efficiency of the proposed algorithm is shown by considering the effects of excitation amplitude, of the measurement units, of larger sampling time interval and of measurement noise. The cases of exact‐, under‐, over‐parameterization of the structural model have been analysed. The proposed algorithm is also quite effective in identifying time‐varying structural parameters to simulate cumulative damage in structural systems. Copyright © 2001 John Wiley & Sons, Ltd.     

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

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

Identification of linear structural systems using earthquake-induced vibration data

This paper addresses the issue of system identification for linear structural systems using earthquake induced time histories of the structural response. The proposed methodology is based on the Eigensystem Realization Algorithm (ERA) and on the Observer/Kalman filter IDentification (OKID) approach to perform identification of structural systems using general input–output data via Markov parameters. The efficiency of the proposed technique is shown by numerical examples for the case of eight-storey building finite element models subjected to earthquake excitation and by the analysis of the data from the dynamic response of the Vincent-Thomas cable suspension bridge (Long Beach, CA) recorded during the Whittier and the Northridge earthquakes. The effects of noise in the measurements and of inadequate instrumentation are investigated. It is shown that the identified models show excellent agreement with the real systems in predicting the structural response time histories when subjected to earthquake-induced ground motion. Copyright © 1999 John Wiley & Sons, Ltd.     

Stiffness‐damping simultaneous identification using limited earthquake records

A new method of stiffness‐damping simultaneous identification of building structures is proposed using limited earthquake records. It is shown that when horizontal accelerations are recorded at the floors just above and below a specific storey in a shear building model, the storey stiffness and the damping ratio can be identified uniquely. The viscous damping coefficient and the linear hysteretic damping ratio can also be identified simultaneously in a numerical model structure. The accuracy of the present identification method is investigated through the actual limited earthquake records in a base‐isolated building. It is further shown that an advanced identification technique for mechanical properties of a Maxwell‐type model can be developed by combining the present method with a perturbation technique. Copyright © 2000 John Wiley & Sons, Ltd.     

APPLICATION OF OFF-LINE AND ON-LINE IDENTIFICATION TECHNIQUES TO BUILDING SEISMIC RESPONSE DATA

The objectives of this paper are to present a comparison of the dynamic characteristics of a seven-storey reinforced concrete building (Van Nuys–Holiday Inn) identified from four recorded strong-motion response data (Whittier earthquake, Landers earthquake, Big Bear earthquake and Northridge earthquake). In the analysis, time-domain methods for estimating the system parameters and the modal properties of the building are studied. Both off-line and on-line identification algorithms are applied to these seismic response data. Under the assumption of a linear time-invariant system the ARX model and ARMAX model are used. Comparison of the identification results using different models are made. In addition, recursive procedures are adapted as on-line identification and the time-varying modal parameters are estimated. For structural systems under strong earthquake excitation, a recursive identification method, adaptive forgetting through multiple models (AFMM), is introduced to identify systems with rapidly changing parameters. Through the analysis of the seismic response data of the building subjected to four earthquakes the identification algorithm and the identification results are discussed.     

Vertical earthquake response of megawatt‐sized wind turbine with soil‐structure interaction effects

The dynamic response of a wind turbine on monopile is studied under horizontal and vertical earthquake excitations. The analyses are carried out using the finite element program SAP2000. The finite element model of the structure is verified against the results of shake table tests, and the earthquake response of the soil model is verified against analytical solutions of the steady‐state response of homogeneous strata. The focus of the analyses in this paper is the vertical earthquake response of wind turbines including the soil‐structure interaction effects. The analyses are carried out for both a non‐homogeneous stratum and a deep soil using the three‐step method. In addition, a procedure is implemented which allows one to perform coupled soil‐structure interaction analyses by properly tuning the damping in the tower structure. The analyses show amplification of the ground surface acceleration to the top of the tower by a factor of two. These accelerations are capable of causing damage in the turbine and the tower structure, or malfunctioning of the turbine after the earthquake; therefore, vertical earthquake excitation is considered a potential critical loading in design of wind turbines even in low‐to‐moderate seismic areas. Copyright © 2015 John Wiley & Sons, Ltd.     

地震下考虑桩-土-结构相互作用的输电塔-线体系响应分析

本文研究了考虑桩-土-结构相互作用的输电塔-线体系在地震作用下的响应。根据实际工程,采用ABAQUS有限元软件,建立了考虑桩-土-结构相互作用效应的输电塔-线体系有限元模型。选取不同场地类型下的12条天然地震波,研究了不同地震波激励下考虑桩-土-结构相互作用效应输电塔-线体系动力响应。通过与考虑刚性基础的输电塔-线体系动力响应对比,得到了输电塔的薄弱位置,并提出了基于刚性基础的输电塔抗震放大系数,可为输电塔抗震设计提供参考。     

A unified formulation of the piecewise exact method for inelastic seismic demand analysis including the P‐delta effect

The non‐linear analysis of single‐degree‐of‐freedom (SDOF) systems provides the essential background information for both strength‐based design and displacement‐based evaluation/design methodologies through the development of the inelastic response spectra. The recursive solution procedure called the piecewise exact method, which is efficiently used for the response analysis of linear SDOF systems, is re‐formulated in this paper in a unified format to analyse the non‐linear SDOF systems with multi‐linear hysteresis models. The unified formulation is also capable of handling the P‐delta effect, which generally involves the negative post‐yield stiffness of the hysteresis loops. The attractiveness of the method lies in the fact that it provides the exact solution when the loading time history is composed of piecewise linear segments, a condition that is perfectly satisfied for the earthquake excitation. Based on simple recursive relationships given for positive, negative and zero effective stiffnesses, the unified form of the piecewise exact method proves to be an extremely powerful and probably the best tool for the SDOF inelastic time‐history and response spectrum analysis including the P‐delta effect. A number of examples are presented to demonstrate the implementation of the method. Copyright © 2003 John Wiley & Sons, Ltd.