Nonlinear FE model updating and reconstruction of the response of an instrumented seismic isolated bridge to the 2010 Maule Chile earthquake

Nonlinear finite element (FE) modeling has been widely used to investigate the effects of seismic isolation on the response of bridges to earthquakes. However, most FE models of seismic isolated bridges (SIB) have used seismic isolator models calibrated from component test data, while the prediction accuracy of nonlinear FE models of SIB is rarely addressed by using data recorded from instrumented bridges. In this paper, the accuracy of a state‐of‐the‐art FE model is studied through nonlinear FE model updating (FEMU) of an existing instrumented SIB, the Marga‐Marga Bridge located in Viña del Mar, Chile. The seismic isolator models are updated in 2 phases: component‐wise and system‐wise FEMU. The isolator model parameters obtained from 23 isolator component tests show large scatter, and poor goodness of fit of the FE‐predicted bridge response to the 2010 Mw 8.8 Maule, Chile Earthquake is obtained when most of those parameter sets are used for the isolator elements of the bridge model. In contrast, good agreement is obtained between the FE‐predicted and measured bridge response when the isolator model parameters are calibrated using the bridge response data recorded during the mega‐earthquake. Nonlinear FEMU is conducted by solving single‐ and multiobjective optimization problems using high‐throughput cloud computing. The updated FE model is then used to reconstruct response quantities not recorded during the earthquake, gaining more insight into the effects of seismic isolation on the response of the bridge during the strong earthquake.     

Efficient seismic response data storage and transmission using ARX model‐based sensor data compression algorithm

This paper presents a linear predictor (LP)‐based lossless sensor data compression algorithm for efficient transmission, storage and retrieval of seismic data. Auto‐Regressive with eXogenous input (ARX) model is selected as the model structure of LP. Since earthquake ground motion is typically measured at the base of monitored structures, the ARX model parameters are calculated in a system identification framework using sensor network data and measured input signals. In this way, sensor data compression takes advantage of structural system information to maximize the sensor data compression performance. Numerical simulation results show that several factors including LP order, measurement noise, input and limited sensor number affect the performance of the proposed lossless sensor data compression algorithm concerned. Generally, the lossless data compression algorithm is capable of reducing the size of raw sensor data while causing no information loss in the sensor data. Copyright © 2005 John Wiley & Sons, Ltd.     

Tangential‐projection algorithm for manifold representation in unidentifiable model updating problems

The problem of updating a structural model and its associated uncertainties by utilizing structural response data is addressed. In an identifiable case, the posterior probability density function (PDF) of the uncertain model parameters for given measured data can be approximated by a weighted sum of Gaussian distributions centered at a number of discrete optimal values of the parameters at which some positive measure‐of‐fit function is minimized. The present paper focuses on the problem of model updating in the general unidentifiable case for which certain simplifying assumptions available for identifiable cases are not valid. In this case, the PDF is distributed in the neighbourhood of an extended and usually highly complex manifold of the parameter space that cannot be calculated explicitly. The computational difficulties associated with calculating the highly complex posterior PDF are discussed and a new adaptive algorithm, referred to as the tangential‐projection (TP) algorithm, allowing for an efficient approximate representation of the above manifold and the posterior PDF is presented. Using this approximation, expressions for calculating the uncertain predictive response are established. A numerical example involving noisy data is presented to demonstrate the proposed method. Copyright © 2002 John Wiley & Sons, Ltd.     

飘浮体系斜拉桥纵向抗震计算的单塔模型

通过全桥模型计算结果与单塔模型计算结果的对比,讨论了集中质量在塔上锚固区处堆聚位置,定义了一个结构动力特性和结构地震反应的修正系数,提出了飘浮体系斜拉桥塔纵桥向地震反应计算的单塔模型。多个工程实例的分析结果表明,采用单塔模型的典型误差在15%左右,精度可以满足工程抗震设计的要求,使桥塔的抗震计算得到简化。     

Parameter identification for on‐line model updating in hybrid simulations using a gradient‐based method

To improve the efficiency of model fitting, parameter identification techniques have been actively investigated. Recently, the applications of parameter identification migrated from off‐line model fitting to on‐line model updating. The objective of this study is to develop a gradient‐based method for model updating to advance hybrid simulation also called hybrid test. A novel modification of the proposed method, which can reduce the number of design variables to improve the identification efficiency, is illustrated in detail. To investigate the model updating, simulated hybrid tests were conducted with a 5‐story steel frame equipped with buckling‐restrained braces (BRBs) utilized in the shaking table tests conducted in E‐Defense in Japan in 2009. The calibrated analytical model that was verified with the test results can serve as the reference model. In the simulated hybrid tests, the physical BRB substructure is numerically simulated by utilizing a truss element with the 2‐surface model identical to the part of the reference model. Such numerical verification allows simulation of measurement errors for investigation on the performance of the proposed method. Moreover, the feasibility of sharing the identified parameter values, which were obtained from the physical substructure responses, with the relevant numerical models is also verified with the artificial component responses derived from the physical experiments.     

Dynamic characteristics and seismic response of adjacent buildings linked by discrete dampers

Coupling adjacent buildings using discrete viscoelastic dampers for control of response to low and moderate seismic events is investigated in this paper. The complex modal superposition method is first used to determine dynamic characteristics, mainly modal damping ratio and modal frequency, of damper-linked linear adjacent buildings for practical use. Random seismic response of linear adjacent buildings linked by dampers is then determined by a combination of the complex modal superposition method and the pseudo-excitation method. This combined method can effectively and accurately determine random seismic response of non-classically damped systems in the frequency domain. Parametric studies are finally performed to identify optimal parameters of viscoelastic dampers for achieving the maximum modal damping ratio or the maximum response reduction of adjacent buildings. It is demonstrated that using discrete viscoelastic dampers of proper parameters to link adjacent buildings can reduce random seismic responses significantly. Copyright © 1999 John Wiley & Sons Ltd.     

Theoretical models and recorded response in the estimation of cumulative seismic damage on non‐linear structures

A damage‐function model is proposed for the inelastic response of structures in terms of initial damage and of displacement amplitudes and secant stiffnesses of response cycles. The model is used to obtain an analytical closed‐form solution for the probability distribution of cumulative damage after an earthquake ground motion given the distribution prior to such excitation and information on the inelastic structural response. The formulation is applied to a reinforced concrete frame and the results show the capabilities of the method to yield updated distributions of damage. Copyright © 2000 John Wiley & Sons, Ltd.     

多维地震输入下曲线匝道桥梁伸缩缝碰撞响应

针对曲线匝道桥梁伸缩缝地震碰撞破坏现象,依托某市多层互通式立交体系中一座独立、单支、多联、曲线、箱梁匝道桥,用Kelvin接触单元模拟了伸缩缝处地震和碰撞效应,建立了全桥空间动力分析模型。利用非线性时程分析法,分别输入单维(X方向、Y方向)和多维共8种地震动工况,分析了曲线匝道桥梁伸缩缝的地震碰撞响应差异。结果表明:曲线匝道桥梁由于曲率的影响,不同方向地震动响应存在耦合效应,其伸缩缝碰撞响应的最不利值应考虑多维地震输入工况,并用碰撞力和碰撞位移双指标控制伸缩缝的地震碰撞响应。     

Effects of seismic isolation on the seismic response of a California high‐speed rail prototype bridge with soil‐structure and track‐structure interactions

With the launch of the high‐speed train project in California, the seismic risk is a crucial concern to the stakeholders. To investigate the seismic behavior of future California High‐Speed Rail (CHSR) bridge structures, a 3D nonlinear finite‐element model of a CHSR prototype bridge is developed. Soil‐structure and track‐structure interactions are accounted for in this comprehensive numerical model used to simulate the seismic response of the bridge and track system. This paper focuses on examining potential benefits and possible drawbacks of the a priori promising application of seismic isolation in CHSR bridges. Nonlinear time history analyses are performed for this prototype bridge subjected to two bidirectional horizontal historical earthquake ground motions each scaled to two different seismic hazard levels. The effect of seismic isolation on the seismic performance of the bridge is investigated through a detailed comparison of the seismic response of the bridge with and without seismic isolation. It is found that seismic isolation significantly reduces the deck acceleration and the force demand in the bridge substructure (i.e., piers and foundations), especially for high‐intensity earthquakes. However, seismic isolation increases the deck displacement (relative to the pile cap) and the stresses in the rails. These findings imply that seismic isolation can be promisingly applied to CHSR bridges with due consideration of balancing its beneficial and detrimental effects through using appropriate isolators design. The optimum seismic isolator properties can be sought by solving a performance‐based optimum seismic design problem using the nonlinear finite‐element model presented herein. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of vertical seismic response for a large‐scale beam‐supported aqueduct

By the theories of potential flow and structural vibration, the formulae for evaluating the ‘wet’ (with water) frequencies and mode shapes of the beam‐supported aqueduct are derived through a simplified fluid‐structure interaction analysis. The time‐history formulae of structural responses to the vertical seismic excitation are obtained. Applying the response‐spectrum principle, the equivalent vertical earthquake load exerted on the beam and the corresponding effects are also derived. Several illustrative examples are conducted. The analytical results show that: (i) The ‘wet’ frequencies of the structure are lower than the corresponding ‘dry’ (without water) frequencies due to the participating water mass, but the ‘wet’ mode shapes are identical to the corresponding ‘dry’ ones. (ii) The water mass plays an important role in the vertical seismic response, which varies with the different geological sites. For the different seismic inputs, the deeper the water is, the greater are the structural responses. (iii) The vertical seismic effects on the beam are generally not too small to be neglected and should be considered in the structural designs of a beam‐supported aqueduct. Copyright © 2002 John Wiley & Sons, Ltd.     

Improved dynamic correction method in seismic analysis of both classically and non‐classically damped structures

A step‐by‐step approximate procedure taking into consideration high‐frequency modes, usually neglected in the modal analysis of both classically and non‐classically damped structures, is presented. This procedure can be considered as an extension of traditional modal correction methods, like the mode‐acceleration method and the dynamic correction method, which are very effective for structural systems subjected to forcing functions described by analytical laws. The proposed procedure, herein called improved dynamic correction method, requires two steps. In the first step, the number of differential equations of motion are reduced and consequently solved by using the first few undamped mode‐shapes. In the second step, the errors due to modal truncation are reduced by correcting the dynamic response and solving a new set of differential equations, formally similar to the original differential equations of motion. The difference between the two groups of differential equations lies in the forcing vector, which is evaluated in such a way as to correct the effects of modal truncation on applied loads. Copyright © 2001 John Wiley & Sons, Ltd.     

Optimal design of tuned mass dampers for a multi‐storey cross laminated timber building against seismic loads

In this paper, the effectiveness of different design solutions for tuned mass dampers (TMD) applied to high‐rise cross‐laminated (X‐Lam) timber buildings as a means to reduce the seismic accelerations was investigated. A seven‐storey full‐scale structure previously tested on shaking table was used as a reference. The optimal design parameters of the TMDs, i.e. damping and frequency ratios, were determined by using a genetic algorithm on a simplified model of the reference structure, composed by seven masses each representing one storey. The optimal solutions for the TMDs were then applied to a detailed finite element model of the seven‐storey building, where the timber panels were modelled with shell elements and the steel connectors with linear spring. By comparing the numerical results of the building with and without multiple TMDs, the improvement in seismic response was assessed. Dynamic time‐history analyses were carried out for a set of seven natural records, selected in accordance with Eurocode 8, on the simplified model, and for Kobe earthquake ground motion on the detailed model. Results in terms of acceleration reduction for different TMD configurations show that the behaviour of the seven‐storey timber building can be significantly improved, especially at the upper storeys. Copyright © 2016 John Wiley & Sons, Ltd.     

Full‐scale dynamic response of an RC building under weak seismic motions using earthquake recordings,ambient vibrations and modelling

In countries with a moderate seismic hazard, the classical methods developed for strong motion prone countries to estimate the seismic behaviour and subsequent vulnerability of existing buildings are often inadequate and not financially realistic. The main goals of this paper are to show how the modal analysis can contribute to the understanding of the seismic building response and the good relevancy of a modal model based on ambient vibrations for estimating the structural deformation under weak earthquakes. We describe the application of an enhanced modal analysis technique (frequency domain decomposition) to process ambient vibration recordings taken at the Grenoble City Hall building (France). The frequencies of ambient vibrations are compared with those of weak earthquakes recorded by the French permanent accelerometric network (RAP) that was installed to monitor the building. The frequency variations of the building under weak earthquakes are shown to be less (∼2%) and therefore ambient vibration frequencies are relevant over the elastic domain of the building. The modal parameters extracted from ambient vibrations are then used to determine the 1D lumped‐mass model in order to reproduce the inter‐storey drift under weak earthquakes and to fix a 3D numerical model that could be used for strong earthquakes. The correlation coefficients between data and synthetic motion are close to 80 and 90% in horizontal directions, for the 1D and 3D modelling, respectively. Copyright © 2009 John Wiley & Sons, Ltd.     

Modal response analysis of multi‐support structures using a random vibration approach

A formulation is developed for modal response analysis of multi‐support structures using a random vibration approach. The spectral moments of the structural response are rigorously decomposed into contributions from spectral moments of uncoupled modal responses. An advantage of the proposed formulation is that the total dynamic response can be obtained on the basis of mode by mode uncoupled analyses. The contributions to the total response from modal responses under individual support ground motions and under cross‐correlated pairs of support ground motions can be recognized explicitly. The application and performance of the formulation is illustrated by means of an example using a well‐established coherency spectrum model and widely known power spectra models, such as white noise and Kanai–Tajimi. The first three spectral moments of displacement, shear, and bending moment responses are computed, showing that the formulation produces the same results as the exact solution. Copyright © 2015 John Wiley & Sons, Ltd.     

Envelope‐based pushover analysis procedure for the approximate seismic response analysis of buildings

An envelope‐based pushover analysis procedure is presented that assumes that the seismic demand for each response parameter is controlled by a predominant system failure mode that may vary according to the ground motion. To be able to simulate the most important system failure modes, several pushover analyses need to be performed, as in a modal pushover analysis procedure, whereas the total seismic demand is determined by enveloping the results associated with each pushover analysis. The demand for the most common system failure mode resulting from the ‘first‐mode’ pushover analysis is obtained by response history analysis for the equivalent ‘modal‐based’ SDOF model, whereas demand for other failure modes is based on the ‘failure‐based’ SDOF models. This makes the envelope‐based pushover analysis procedure equivalent to the N2 method provided that it involves only ‘first‐mode’ pushover analysis and response history analysis of the corresponding ‘modal‐based’ SDOF model. It is shown that the accuracy of the approximate 16th, 50th and 84th percentile response expressed in terms of IDA curves does not decrease with the height of the building or with the intensity of ground motion. This is because the estimates of the roof displacement and the maximum storey drift due to individual ground motions were predicted with a sufficient degree of accuracy for almost all the ground motions from the analysed sets. Copyright © 2013 John Wiley & Sons, Ltd.     

Development and validation of p‐y modeling approach for seismic response predictions of highway bridges

This study aims to realistically simulate the seismic responses of typical highway bridges in California with considerations of soil–structure interaction effects. The p‐y modeling approaches are developed and validated for embankments and pile foundations of bridges. The p‐y approach models the lateral and vertical foundation flexibility with distributed p‐y springs and associated t‐z and q‐z springs. Building upon the existing p‐y models for pile foundations, the study develops the nonlinear p‐y springs for embankments based on nonlinear 2D and 3D continuum finite element analysis under passive loading condition along both longitudinal and transverse directions. Closed‐form expressions are developed for two key parameters, the ultimate resistant force p_ult and the displacement y₅₀, where 0.5p_ult is reached, of embankment p‐y models as functions of abutment geometry (wall width and height, embankment fill height, etc.) and soil material properties (wall‐soil friction angle, soil friction angle, and cohesion). In order to account for the kinematic and site responses, depth‐varying ground motions are derived and applied at the free‐end of p‐y springs, which reflects the amplified embankment crest motion. The modeling approach is applied to simulate the seismic responses of the Painter Street Bridge and validated through comparisons with the recorded responses during the 1992 Petrolia earthquake. It is demonstrated that the flexibility and motion amplification at end abutments are the most crucial modeling aspects. The developed p‐y models and the modeling approach can effectively predict the seismic responses of highway bridges. Copyright © 2016 John Wiley & Sons, Ltd.     

Verification and validation of a seismic response analysis code for framed structures using the ASI‐Gauss technique

In this report, the capabilities of the adaptively shifted integration (ASI)‐Gauss code in the analysis of the seismic responses of framed structures are verified and validated by comparing the results with detailed numerical simulations performed by the parallel finite element analysis code, E‐Simulator, and with experimental results obtained by E‐Defense. The numerical results obtained by both codes showed good agreement with the experimental results obtained by E‐Defense. Furthermore, seismic waves with unnaturally large magnitudes are applied to a high‐rise building model to demonstrate the ability of the ASI‐Gauss code to analyze the collapse behaviors of building frames. Copyright © 2013 John Wiley & Sons, Ltd.     

Three‐dimensional models of reservoir sediment and effects on the seismic response of arch dams

The important effects of bottom sediments on the seismic response of arch dams are studied in this paper. To do so, a three‐dimensional boundary element model is used. It includes the water reservoir as a compressible fluid, the dam and unbounded foundation rock as viscoelastic solids, and the bottom sediment as a two‐phase poroelastic domain with dynamic behaviour described by Biot's equations. Dynamic interaction among all those regions, local topography and travelling wave effects are taken into account. The results obtained show the important influence of sediment compressibility and permeability on the seismic response. The former is associated with a general change of the system response whereas the permeability has a significant influence on damping at resonance peaks. The analysis is carried out in the frequency domain considering time harmonic excitation due to P and S plane waves. The time‐domain results obtained by using the Fourier transform for a given earthquake accelerogram are also shown. The possibility of using simplified models to represent the bottom sediment effects is discussed in the paper. Two alternative models for porous sediment are tested. Simplified models are shown to be able to reproduce the effects of porous sediments except for very high permeability values. Copyright © 2004 John Wiley & Sons, Ltd.     

Near‐field seismic effects in a homogeneous medium and their removal in vertical seismic profile attenuation estimates

To better understand (and correct for) the factors affecting the estimation of attenuation (Q), we simulate subsurface wave propagation with the Weyl/Sommerfeld integral. The complete spherical wavefield emanating from a P‐wave point source surrounded by a homogeneous, isotropic and attenuative medium is thus computed. In a resulting synthetic vertical seismic profile, we observe near‐field and far‐field responses and a 90° phase rotation between them. Depth dependence of the magnitude spectra in these two depth regions is distinctly different. The logarithm of the magnitude spectra shows a linear dependence on frequency in the far‐field but not in those depth regions where the near‐field becomes significant. Near‐field effects are one possible explanation for large positive and even negative Q‐factors in the shallow section that may be estimated from real vertical seismic profile data when applying the spectral ratio method. We outline a near‐field compensation technique that can reduce errors in the resultant Q estimates.     

Dynamic characteristics and seismic behavior of prefabricated steel stairs in a full‐scale five‐story building shake table test program

This paper investigates the dynamic characteristics and seismic behavior of prefabricated steel stairs in a full‐scale five‐story building shake table test program. The test building was subjected to a suite of earthquake input motions and low‐amplitude white noise base excitations first, while the building was isolated at its base, and subsequently while it was fixed to the shake table platen. This paper presents the modal characteristics of the stairs identified using the data recorded from white noise base excitation tests as well as the physical and measured responses of the stairs from the earthquake tests. The observed damage to the stairs is categorized into three distinct damage states and is correlated with the interstory drift demands of the building. These shake table tests highlight the seismic vulnerability of modern designed stair systems and in particular identifies as a key research need the importance of improving the deformability of flight‐to‐building connections. Copyright © 2015 John Wiley & Sons, Ltd.