A mathematical model to predict the inelastic response of a steel frame: Formulation of the model

The purpose of this research is to use data from experiments to formulate a mathematical model that will predict the non-linear response of a single-storey steel frame to an earthquake input. The process used in this formulation is system identification. The form of the model is a second-order non-linear differential equation with linear viscous damping and Ramberg—Osgood type hysteresis. The damping coefficient and the three parameters in the hysteretic model are to be established. An integral weighted mean squared error function is used to evaluate the [goodness of fit] between the model's response and the structure's response when both are subjected to the same excitation. The function includes errors in displacement and acceleration and is integrated from zero to a time T, which may be the full duration of the recorded response or only a portion of it. The parameters are adjusted using a modified Gauss-Newton method until the error function is minimized. The computer program incorporating these steps in the system identification process is verified with simulated data. Results given in the paper show that in every case the program converges in few iterations to the assigned set of parameters.     

Application of singular spectrum analysis to identify the degrading structure using deteriorating distributed element model

Reinforced concrete structure may exhibit significant inelastic hysteretic behavior when subject to strong earthquake excitation. To investigate such an inelastic behavior, in this study, a new system identification technique is applied by using the deteriorating distributed element (DDE) model to simulate the hysteretic behavior of a degrading structure. Through the advanced signal processing technique, the multiple singular spectrum analysis (SSA) and the nonlinear SSA, the recorded inelastic restoring force of a deteriorating structure can be decomposed into several independent additive components in its sequentially degrading order and with decreasing weight. With each decomposed hysteresis loop, the model parameters of the DDE model are identified. The evolutionary properties of the progressive stiffness degradation behavior of reinforced concrete structure can be observed from the identified model parameters. Finally, comparison on the physical properties of the identified DDE model with respect to the seismic response data of the deteriorating structure is also discussed. The result shows that the proposed identification technique can have a good estimation on the seismic behavior of the degrading structure. Copyright © 2012 John Wiley & Sons, Ltd.     

Identification of Fei-Tsui arch dam from both ambient and seismic response data

This purpose of this paper is to study the dynamic characteristics of the Fei-Tsui arch dam using the seismic response data and the ambient vibration data. For the identification of dam properties from seismic response data, the multiple inputs from the abutment of the dam to represent the nonuniform excitations of seismic input motion are considered, and the ARX model is applied using the discrete-time linear filtering approach with least-squares approximation to identify the dynamic characteristics of the dam. The system modal dampings, natural frequencies and frequency response functions are identified. A comparison of the identified modal parameters is made among different seismic events. Post-earthquake safety evaluation of the dam can be made based on the identified model. Finally, the ambient vibration test of the dam is performed to identify the mode shapes along the dam crest.     

Base isolation systems for earthquake protection of multi-storey shear structures

This paper is a study of the effectiveness of a wide range of bilinear hysteretic isolation systems in shielding multistorey 2-D shear structures from earthquake excitations. Important parameters of the isolation system are identified and their effect on structure response noted. It is shown that isolation systems can be constructed which allow the structure proper to remain purely elastic even during very strong ground motions. It is further shown that the shear responses and base displacements of structures on these isolation systems can be accurately estimated from elastic response spectra of the forcing earthquakes. The philosophy of structure isolation is discussed and an introduction given to the physical devices currently available to provide it.     

An identification methodology for a class of hysteretic structures

This paper presents a method of identification for determining non-linear dynamic models for certain hysteretic structures. Particular attention is given to modelling and identifying the hysteretic behaviour of structures from strong-motion earthquake data. In this method, the response is separated into mode-like components which are analogous to those of a linear system. Based on modelling of the generalized restoring force of each mode-like component, both non-hysteretic and hysteretic non-linear models are incorporated into the general methodology. A non-hysteretic model provides an initial estimate for a final hysteretic model. The approach is applicable even when data are available from only a small number of locations in the structure. The structural model identified from this method provides a means to predict the response to future events and, ultimately, to examine the damage to a structure as a result of an earthquake.     

各向异性介质弹性波多参数全波形反演

各向异性介质弹性波方程全波形反演过程中多参数之间的相互耦合,使得弱参数在反演过程中难得到理想的结果.本文以VTI介质为例,在各参数辐射模式分析的基础上,基于改进的散射积分算法实现目标函数梯度的直接求取,进一步构建高斯牛顿方向,实现Hessian矩阵的有效利用,以考虑Hessian矩阵非主对角线元素包含的各参数间的耦合效应,在不使用任何反演策略的情况下实现高精度的VTI介质弹性波方程多参数同步反演.同时,该方法在计算过程中无需存储庞大的核函数矩阵,且无需传统截断牛顿法中额外的正演计算,因此内存占用小,计算效率高.本文数值试验验证了该方法的有效性,为各向异性多参数全波形反演提供了一种新的解决方案.     

Parameter identification of hysteretic model of rubber-bearing based on sequential nonlinear least-square estimation

In order to evaluate the nonlinear performance and the possible damage to rubber-bearings (RBs) during their normal operation or under strong earthquakes, a simplified Bouc-Wen model is used to describe the nonlinear hysteretic behavior of RBs in this paper, which has the advantages of being smooth-varying and physically motivated. Further, based on the results from experimental tests performed by using a particular type of RB (GZN110) under different excitation scenarios, including white noise and several earthquakes, a new system identification method, referred to as the sequential nonlinear least-square estimation (SNLSE), is introduced to identify the model parameters. It is shown that the proposed simplified Bouc-Wen model is capable of describing the nonlinear hysteretic behavior of RBs, and that the SNLSE approach is very effective in identifying the model parameters of RBs.     

On the estimation of hysteretic energy demands for SDOF systems

Based on a statistical study of the dynamic response of single degree of freedom (SDOF) systems subjected to earthquake ground motions, a rule to estimate hysteretic energy demands is proposed. Expressions for elastic–perfectly plastic, stiffness degrading and pinching SDOF systems were obtained. The proposed rule does not explicitly depend on the period of the system; instead, it depends on the elastic pseudo‐acceleration and elastic velocity spectra. It is shown that the proposed rule can be applied to compute hysteretic energy demands for systems located at different soil conditions. In addition, information about scatter and bias of the proposed rule is also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Condition assessment of structures under unknown support excitation

A new method is proposed to assess the condition of structures under unknown support excitation by simultaneously detecting local damage and identifying the support excitation from several structural dynamic responses. The support excitation acting on a structure is modeled by orthogonal polynomial approximations, and the sensitivities of structural dynamic response with respect to its physical parameters and orthogonal coeffi cients are derived. The identifi cation equation is based on Taylor’s fi rst orde...     

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.     

A two‐step direct method for estimating the seismic response of nonlinear structures equipped with nonlinear viscous dampers

The insertion of fluid viscous dampers in building structures is an innovative technology that can improve significantly the seismic response. These devices could be very useful also in the retrofit of existing buildings. The effect of this typology of damping system is usually identified with an equivalent supplemental damping ratio, which depends on the maximum displacement of the structure, so that iterative procedures are required. In this paper, a simplified direct assessment method for nonlinear structures equipped with nonlinear fluid viscous dampers is proposed. The method proposed in this study is composed by two steps. The first one yields the direct estimate of the supplemental damping ratio provided by nonlinear viscous dampers in presence of a linear elastic structural response. The second step extends the procedure to structures with nonlinear behavior. Both graphical and analytical approaches have been developed. The proposed method has then been verified through several applications and comparisons with nonlinear dynamic analyses. Moreover, an investigation has been performed with regard to the influence of the relations that define the damping reduction factor and the hysteretic damping. Copyright © 2014 John Wiley & Sons, Ltd.     

基于正交分解的非线性结构外荷载识别方法研究

提出一种基于切比雪夫正交分解的非线性结构外荷载识别方法及分解阶数确定办法.在识别过程中建立非线性结构体系状态空间方程,并将切比雪夫正交多项式展开系数扩展于状态量,对状态量进行递推估计.通过结构反应频域分析筛选频率范围并确定正交多项式项数.文中将通过6层隔震结构、波形钢腹板PC 组合梁桥的数值仿真和3层隔震框架的振动台试...     

Bivariate Gaussian mixture–based equivalent linearization method for stochastic seismic analysis of nonlinear structures

To address challenges in stochastic seismic analysis of nonlinear structures, this paper further develops a recently proposed Gaussian mixture–based equivalent linearization method (GM‐ELM). The GM‐ELM uses a Gaussian mixture distribution model to approximate the probabilistic distribution of a nonlinear system response. Using properties of the Gaussian mixture model, GM‐ELM can decompose the non‐Gaussian response of a nonlinear system into multiple Gaussian responses of linear single–degree of freedom oscillators. With the set of the equivalent linear systems identified by GM‐ELM, response statistics as crossing rate and first‐passage probability can be computed conveniently using theories of linear random vibration analysis. However, the original version of GM‐ELM may lead to an inaccurate estimate because of the heuristic parameters of the linear system introduced to supplement insufficient information. To overcome this limitation and define unique equivalent linear systems, this paper proposes a further developed version of GM‐ELM, which uses a mixture of bivariate Gaussian densities instead of univariate models. Moreover, to facilitate the use of elastic response spectra for estimating the mean peak responses of a nonlinear structure, a new response spectrum combination rule is proposed for GM‐ELM. Two numerical examples of hysteretic structural systems are presented in this paper to illustrate the application of the bivariate GM‐ELM to nonlinear stochastic seismic analysis. The analysis results obtained by the bivariate GM‐ELM are compared with those obtained by the univariate GM‐ELM, the conventional equivalent linearization method, the tail equivalent linearization method, and Monte Carlo simulation. The supporting source code and data are available for download at https://github.com/yisangri/GitHub‐bGM‐ELM‐code.git     

Inelastic seismic response of stiffening systems and development of demand spectrum: Application to MSSS bridges

A stiffening system is a system that increases its stiffness as it goes under large displacements. Such behavioural characteristic can result from constitutive behaviour or at the structural level often from closure of gaps between various components (sub‐systems) of the structure. An example of the latter situation is multi‐span simply supported (MSSS) bridges under horizontal earthquake ground motion. Unlike softening systems, stiffening systems have not been studied. In addition to the need for more understanding of the seismic response of stiffening systems, there is a need to develop response spectrum that can be used in design. Several parameters including gap size and ratios of sub‐systems stiffness, strength, and mass control the behaviour of a stiffening system. In this study, a simplified stiffening model is developed and over 367 000 cases are analysed to investigate the nonlinear stiffening behaviour and pounding. Parameters considered also include ground motion characteristic. Results are evaluated and compared in terms of displacement and dissipated hysteretic energy. Parameter study results show that, on average, the displacement response is lower for stiffening systems, however, they dissipates higher hysteretic energy, due to higher yield cycles and yield excursions, and can possibly sustain more damage than a bilinear, elastic–plastic system. Using parameter study database, design response spectrum for stiffening systems is also proposed and its practical application is demonstrated through its application to an MSSS bridge. Results of this study goes beyond MSSS bridges and will have application for many structural systems where response is characterized by a stiffening behaviour. Copyright © 2007 John Wiley & Sons, Ltd.     

System identification and damage evaluation of degrading hysteresis of reinforced concrete frames

In this study, signal processing approaches and nonlinear identification are used to measure seismic responses of reinforced concrete (RC) structures using the shaking table test. To analyze structural nonlinearity, an equivalent linear system with time‐varying model parameters, singular spectrum analysis to elucidate residual deformation, and wavelet packet transformation analysis to yield the energy distribution among components are adopted to detect the nonlinearity. Then, damage feature extraction is conducted using both the Holder exponent and the Level‐1 detail of the discrete wavelet component. Finally, the modified Bouc‐Wen hysteretic model and the system identification process are employed to the shaking table test data to evaluate the physical parameters, including the stiffness degradation, the strength deterioration and the pinching hysteresis. Finally, the identified stiffness and strength degradation functions from the test data of RC frames in relation to the degree of ground shaking, damage index and the identified nonlinear features are discussed. Based on the proposed method, both signal‐based and model‐based identifications, the relationship between the damage occurrence and severity of structural damage can be identified. Copyright © 2010 John Wiley & Sons, Ltd.     

Seismic response of self‐centring hysteretic SDOF systems

The seismic response of single‐degree‐of‐freedom (SDOF) systems incorporating flag‐shaped hysteretic structural behaviour, with self‐centring capability, is investigated numerically. For a SDOF system with a given initial period and strength level, the flag‐shaped hysteretic behaviour is fully defined by a post‐yielding stiffness parameter and an energy‐dissipation parameter. A comprehensive parametric study was conducted to determine the influence of these parameters on SDOF structural response, in terms of displacement ductility, absolute acceleration and absorbed energy. This parametric study was conducted using an ensemble of 20 historical earthquake records corresponding to ordinary ground motions having a probability of exceedence of 10% in 50 years, in California. The responses of the flag‐shaped hysteretic SDOF systems are compared against the responses of similar bilinear elasto‐plastic hysteretic SDOF systems. In this study the elasto‐plastic hysteretic SDOF systems are assigned parameters representative of steel moment resisting frames (MRFs) with post‐Northridge welded beam‐to‐column connections. In turn, the flag‐shaped hysteretic SDOF systems are representative of steel MRFs with newly proposed post‐tensioned energy‐dissipating connections. Building structures with initial periods ranging from 0.1 to 2.0s and having various strength levels are considered. It is shown that a flag‐shaped hysteretic SDOF system of equal or lesser strength can always be found to match or better the response of an elasto‐plastic hysteretic SDOF system in terms of displacement ductility and without incurring any residual drift from the seismic event. Copyright © 2002 John Wiley & Sons, Ltd.     

Investigation of the elastic characteristics of a three-storey steel structure using system identification

In this report, three different models in increasing order of complexity have been used to identify the seismic behaviour of a three-storey steel structure subjected to arbitrary forcing functions, all of which excite responses within the elastic range. All of the models are constructed using system identification. In the first model, five parameters have been used to identify the frame. Treating the system as a shear building, we assign one stiffness coefficient to each floor and introduce Rayleigh-type damping with two additional parameters. The mass, assumed to be concentrated at a floor level, is kept constant throughout the study. The parameters are established using a modified Gauss-Newton algorithm. The match between measured and predicted quantities is satisfactory when these quantities are restricted to floor accelerations or displacements. To remove the constraint imposed by assuming that the frame deforms as a shear building, a second model with eight parameters is introduced, allowing rotations of the joints as independent degrees of freedom. Six of the eight parameters are related to the stiffness characteristics of the structural members while the remaining two are related to damping as before. In constructing the eight-parameter model, we learned that it is the effective lengths of the members that change during optimization. We also found that the independent response quantities, floor accelerations and joint rotations, must be used in the cost function for the optimization algorithm to converge. The match between measured and predicted quantities for the eight-parameter model is excellent. The set of parameters derived from the minimum squared error gives a model that shows very good correlation using information on the full duration of the pulse or only a portion of it. Also the same correlation exists between the coefficients obtained from different excitations. In an effort to explain the values of the parameters associated with the girders, an additional degree of freedom, namely, the pitching motion of the shaking table, is introduced as an additional degree of freedom. The paper presents, therefore, a five-, an eight- and, finally, a nine-parameter model.     

System identification of linear structures based on Hilbert–Huang spectral analysis. Part 1: normal modes

Based on the Hilbert–Huang spectral analysis, a method is proposed to identify multi‐degree‐of‐freedom (MDOF) linear systems using measured free vibration time histories. For MDOF systems, the normal modes have been assumed to exist. In this method, the measured response data, which are polluted by noises, are first decomposed into modal responses using the empirical mode decomposition (EMD) approach with intermittency criteria. Then, the Hilbert transform is applied to each modal response to obtain the instantaneous amplitude and phase angle time histories. A linear least‐square fit procedure is proposed to identify the natural frequency and damping ratio from the instantaneous amplitude and phase angle for each modal response. Based on a single measurement of the free vibration time history at one appropriate location, natural frequencies and damping ratios can be identified. When the responses at all degrees of freedom are measured, the mode shapes and the physical mass, damping and stiffness matrices of the structure can be determined. The applications of the proposed method are illustrated using three linear systems with different dynamic characteristics. Numerical simulation results demonstrate that the proposed system identification method yields quite accurate results, and it offers a new and effective tool for the system identification of linear structures in which normal modes exist. Copyright © 2003 John Wiley & Sons, Ltd.     

Mechanical formulations for bilinear and trilinear hysteretic models used in base isolators

The best known model for numerically simulating the hysteretic behavior of various structural components is the bilinear hysteretic system. There are two possible mechanical formulations that correspond to the same bilinear model from a mathematical viewpoint. The first one consists of a linear elastic spring connected in series with a parallel system comprising a plastic slider and a linear elastic spring, while the second one comprises a linear elastic spring connected in parallel with an elastic-perfectly plastic system. However, the bilinear hysteretic model is unable to describe either softening or hardening effects in these components. In order to account for this, the bilinear model is extended to a trilinear one. Thus, two trilinear hysteretic models are developed and numerically tested, and the analysis shows that both exhibit three plastic phases. More specifically, the first system exhibits one elastic phase, while the second one exhibits two elastic phases according to the level of strain amplitude. Next, the change of slope between the plastic phases in unloading does not occur at the same displacement level in the two models. Furthermore, the dissipated energy per cycle in the first trilinear model, as proven mathematically and explained physically, decreases in the case of hardening and increases in the case of softening, while in the second trilinear model the dissipated energy per cycle remains unchanged, as is the case with the bilinear model. Numerical examples are presented to quantify the aforementioned observations made in reference to the mechanical behavior of the two trilinear hysteretic models. Finally, a set of cyclic shear tests over a wide range of strain amplitudes on a high damping rubber bearing is used in the parameter identification of the two different systems, namely (a) trilinear hysteretic models of the first type connected in parallel, and (b) trilinear hysteretic models of the second type also connected in parallel. The results show that the complex nonlinear shear behavior of high damping rubber bearings can be correctly simulated by a parallel system which consists of only one component, namely the trilinear hysteretic system of the first type. The second parallel system was not able to describe the enlargement of the dissipated hysteresis area for large strain amplitudes.