System identification of a concrete arch dam and calibration of its finite element model

A magnitude 4.3 earthquake occurred near Pacoima Dam on 13 January 2001. An accelerometer array that had been upgraded after the Northridge earthquake recorded the motion with 17 channels on the dam and the dam–foundation interface. Using this data, properties of the first two modes are found from a system identification study. Modal properties are also determined from a forced vibration experiment performed in 2002 and indicate a significantly stiffer system than is estimated from the 2001 earthquake records. The 2001 earthquake, although small, must have induced temporary nonlinearity. This has implications for structural health monitoring. The source of the nonlinear behaviour is believed to be loss of stiffness in the foundation rock. A finite element model of Pacoima Dam is constructed and calibrated to match modal properties determined from the system identification study. A dynamic simulation of the 2001 earthquake response produces computed motions that agree fairly well with the recorded ones. Copyright © 2006 John Wiley & Sons, Ltd.     

A numerical approach for modeling near-fault ground motion and its application in the 1994 Northridge earthquake

An approach for simulating near-fault ground motion was presented by combining the finite fault model with a numerical algorithm, named investigated lump method presented in this paper for wave propagation. The investigated lumps are constructed from the auxiliary quadrilateral grids. The dynamic equilibrium equations of a typical investigated lump have been derived and obtained by integrating the stresses along the contour of the investigated lump. The stresses are calculated using the constitutive relations and the interpolation techniques. The investigated lump method is then implemented using the equilibrium equations of investigated lumps and the calculations of stresses alternately in time domain. The stability criterion of the algorithm has been given. Comparisons with the discrete wave-number method solutions for predicting the ground motions at the Pacoima Dam during the San Fernando earthquake show the validity of the method presented in this paper for simulating near-field ground motions. A finite fault source model has been implemented in the algorithm here. The source parameters given by Wald et al. (1996) [18] are applied to synthesize the ground motions at three stations during the 1994 Northridge earthquake. The simulating results qualitatively match to the corresponding ground motion records. The studies demonstrated that the approach presented in this paper is an effective tool for the numerical simulation of near-fault ground motion.     

Functional series TARMA modelling and simulation of earthquake ground motion

The non‐stationary Functional Series time‐dependent autoregressive moving average (TARMA) modelling and simulation of earthquake ground motion is considered. Full Functional Series TARMA models, capable of modelling both resonances and antiresonances, are examined for the first time via a novel mixed parametric/non‐parametric estimation scheme, and critical comparisons with pure TAR and recursive ARMA (RARMA)‐recursive maximum likelihood (RML) adaptive filtering type modelling are made. The study is based upon two California ground motion signals: a 1979 El Centro accelerogram and a 1994 Pacoima Dam accelerogram. A systematic analysis, employing various functional subspaces and model orders, leads to two Haar function based models: a TARMA(2,4)₈ model for the El Centro case and a TARMA(6,2)₁₀ model for the Pacoima Dam case. Both models are formally validated and their simulation (synthesis) capabilities are demonstrated via Monte Carlo experiments focusing on important time domain signal characteristics. The Functional Series TAR/TARMA models are shown to achieve parsimony, as well as superior accuracy and simulation capabilities, over their RARMA counterparts. Copyright © 2001 John Wiley & Sons Ltd.     

Response of an arch dam to non-uniform excitation generated by a seismic wave scattering model

Non-uniform ground motions are generated based on a single record available at a site and seismic wave scattering analysis. The Chino Hills 2008 earthquake records at the Pacoima Dam site are used to indicate the accuracy of the method. Dynamic analysis of the Pacoima dam-reservoir-foundation under uniform and non-uniform ground motions is carried out using the EACD-3D2008 software, and the results are compared to recorded responses at different locations on the dam. There is good agreement between computed and recorded displacements of the dam for non-uniform excitation. For uniform excitation, the displacements are underestimated in comparison with those obtained from recorded excitation. Significant intensification of stresses, especially near the foundation, and different patterns of stress distribution are observed for non-uniform excitation in comparison with uniform excitation. For uniform excitation maximum stresses occur in the crown cantilever near the crest, but for non-uniform excitation the maximum stresses occur along the sides and near the foundation.     

近场强地震动数值模拟的简化计算方法

近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动。在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型。通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性     

Stochastic simulation of regionalized ground motions using wavelet packets and cokriging analysis

Performance‐based earthquake engineering often requires ground‐motion time‐history analyses to be performed, but very often, ground motions are not recorded at the location being analyzed. The present study is among the first attempt to stochastically simulate spatially distributed ground motions over a region using wavelet packets and cokriging analysis. First, we characterize the time and frequency properties of ground motions using the wavelet packet analysis. The spatial cross‐correlations of wavelet packet parameters are determined through geostatistical analysis of regionalized ground‐motion data from the Northridge and Chi‐Chi earthquakes. It is observed that the spatial cross‐correlations of wavelet packet parameters are closely related to regional site conditions. Furthermore, using the developed spatial cross‐correlation model and the cokriging technique, wavelet packet parameters at unmeasured locations can be best estimated, and regionalized ground‐motion time histories can be synthesized. Case studies and blind tests using data from the Northridge and Chi‐Chi earthquakes demonstrate that the simulated ground motions generally agree well with the actual recorded data. The proposed method can be used to stochastically simulate regionalized ground motions for time‐history analyses of distributed infrastructure and has important applications in regional‐scale hazard analysis and loss estimation. Copyright © 2014 John Wiley & Sons, Ltd.     

Strong Ground Motion Estimation During the Kutch, India Earthquake

In the absence of strong motion records, ground motion during the 26^th January, 2001 Kutch, India earthquake, has been estimated by analytical methods. A contour map of peak ground acceleration (PGA) values in the near source region is provided. These results are validated by comparing them with spectral response recorder data and field observations. It is found that very near the epicenter, PGA would have exceeded 0.6 g. A set of three aftershock records have been used as empirical Green's functions to simulate ground acceleration time history and 5% damped response spectrum at Bhuj City. It is found that at Bhuj, PGA would have been 0.31 g–0.37 g. It is demonstrated that source mechanism models can be effectively used to understand spatial variability of large-scale ground movements near urban areas due to the rupture of active faults.     

Simulation of broadband seismic ground motions at dam canyons by using a deterministic numerical approach

As a deterministic numerical approach for simulation of earthquake ground motions, the spectral element method (SEM) is applied to generate a broadband acceleration array for dam-canyons instead of the traditional empirical or stochastic methods. Specifically, the SEM analysis model with an extra fine mesh is used for the Pacoima Canyon to simulate the entire path starting from earthquake source rupture via the propagation medium to the local site. The source and the 3D earth model (velocity structure) are validated through the modeling of the Newhall earthquake on 28 October 2012 at a frequency of up to 8 Hz. Subsequently, the San Fernando earthquake records on 13 January 2001 are further used to study the effects of propagation path in simulation. Finally, the spatially varying ground motions at the Pacoima Canyon are obtained for different source mechanisms. The results show that the source mechanism and the local site topography significantly affect the distribution of the peak accelerations along the canyon.     

Earthquake response analysis of the olive view hospital psychiatric day clinic

The collapse of the Olive View Hospital Psychiatric Day Clinic is studied using three biaxial force-deflection models to represent the columns of the building. These models are: shear collapse, elastic and inelastic. The biaxial models for shear and inelastic behaviour are new developments and are useful for non-linear structural dynamic studies. In the present study, the shear collapse model is intended to represent the actual prototype behaviour. The inelastic model, which is based on a hardening rule of plasticity, is used to study the performance of a hypothetical structure with the same storey shear capacity as the prototype but which exhibits ductile behaviour. The prototype structure had a base storey shear capacity of 25 per cent, and actually failed by shearing of all of the first floor columns. In the present study, the shear collapse model predicted this behaviour even with the El Centro accelerogram as input. This result may have far-reaching significance because many low-rise reinforced concrete buildings which were designed according to recent codes have similar storey shear capacity coefficients and column properties. According to this study, such buildings may collapse even in a moderate earthquake. In the inelastic representation, the structure was found to have a base storey shear capacity of 80 per cent when moment hinging was assumed to occur at the top and bottom of the columns. Even with this high strength capacity, the permanent offset computed from the inelastic model corresponded to a ductility factor of 5 when the Pacoima Dam accelerogram was used as input. On the basis of damage to other structures observed on the site, it seems likely that ground motion of about the Pacoima Dam intensity occurred at Olive View. From this it is concluded that a low-rise ductile frame concrete building, even with this high shear force capacity, may not prove satisfactory for hospital use when subjected to strong ground motion.     

基于小田原试验场地数据的不同基岩输入模型对预测地震动特征的影响研究

随着强震台网的密布及观测记录的增加,为研究各类局部场地地震反应预测模型的合理性提供了有效的参考依据,也使利用强震记录及场地条件研究地震动特征成为可能。选取场地地质参数资料和地震记录数据齐全的日本小田原（Ashigara Valley）盲测试验场地,通过对比不同地震动输入方式及场地反应分析模型,研究地震动特征,分析现有模型的优劣。基于1990年8月5日M5.1强震事件的地表基岩记录和地下基岩地震记录,采用地下台强震记录直接输入、地表基岩台强震记录减半为基底地震动输入、地表基岩台强震记录反演为基底地震动输入作为3种基岩地震动输入。基于局部场地条件分别建立一维等效线性模型、二维黏弹性模型及二维时域等效线性化模型等工程中常用的场地数值分析模型,进行局部场地地震反应分析,预测该盲测场地的地表地震动特征,并与对应的实测强震记录结果进行对比,分析不同基岩地震动输入方式对预测地震动特征及地表土层反应谱特征的影响,重点分析地震动输入、土体非线性、场地横向不均匀性及几何与非线性特征共同作用等因素对地表地震动特征的影响,以期为地表地震动的合理预测提供参考。     

长周期地震动强度指标与SDOF体系残余变形相关性分析

为分析两类长周期地震动作用下如何选取合理的强度指标作为结构抗震设计的输入,从2011年东日本9.0级地震、2003年十胜冲8.0级地震和2016年熊本7.3级地震数据库中选取90条可靠的远场长周期地震动。从1994年美国北岭6.7级地震和1999年台湾集集7.6级地震数据库中选取60条近场长周期地震动,以SDOF体系为研究对象,讨论了阻尼比、屈服刚度折减系数和强度折减系数对残余变形与18个地震动强度指标的相关系数的影响,对比分析了两类长周期地震作用下相关系数之间的异同。研究结果表明:在考虑近场长周期地震动作用时,建议应根据结构周期的大小来选择合适的强度指标作为地震动的输入。远场长周期地震动作用下,以PGA、PGV和PGD为代表的强度指标与残余变形的相关程度均较高,PGV稳定性略好于PGA和PGD,建议PGV作为地震动输入的控制指标。残余变形相关系数受阻尼比、屈服刚度折减系数以及强度折减系数影响不大。     

A new solution to estimate the time delay on the topographic site using time domain 3D boundary element method

This study focuses on investigating spatial variation of ground motion that has great influence on the dynamic behavior of the large structures located on the surface topography. One of the most effective parameters on the spatial variation of ground motion is the difference between the arrival time of seismic waves in different points located on the abutments. In this research, a three-dimensional model of the Pacoima Dam site is prepared. The time domain 3 D boundary element method is used to apply non-uniform excitation at the dam supports. This model is subjected to vertically propagating incident SH and P waves. The time delay can be characterized by calculating the value of the time delay for which the cross-correlation between two records is maximized. Finally, to obtain the time delay in a topographic site, a function considering effective parameters such as the height from the canyon base, wave velocity and predominant frequency, is presented. Furthermore, a code was developed for generating the spatially variation of seismic ground motions. The results show that the proposed functions have an acceptable accuracy in estimating the time delay to generate non-uniform ground motion.     

Influence of time‐varying frequency content in earthquake ground motions on seismic response of linear elastic systems

Earthquake ground motion records are nonstationary in both amplitude and frequency content. However, the latter nonstationarity is typically neglected mainly for the sake of mathematical simplicity. To study the stochastic effects of the time‐varying frequency content of earthquake ground motions on the seismic response of structural systems, a pair of closely related stochastic ground motion models is adopted here. The first model (referred to as ground motion model I) corresponds to a fully nonstationary stochastic earthquake ground motion model previously developed by the authors. The second model (referred to as ground motion model II) is nonstationary in amplitude only and is derived from the first model. Ground motion models I and II have the same mean‐square function and global frequency content but different features of time variation in the frequency content, in that no time variation of the frequency content exists in ground motion model II. New explicit closed‐form solutions are derived for the response of linear elastic SDOF and MDOF systems subjected to stochastic ground motion model II. New analytical solutions for the evolutionary cross‐correlation and cross‐PSD functions between the ground motion input and the structural response are also derived for linear systems subjected to ground motion model I. Comparative analytical results are presented to quantify the effects of the time‐varying frequency content of earthquake ground motions on the structural response of linear elastic systems. It is found that the time‐varying frequency content in the seismic input can have significant effects on the stochastic properties of system response. Copyright © 2016 John Wiley & Sons, Ltd.     

基于物理的随机地震动模型研究

基于物理联系研究地震动随机性,建立了随机地震动与基底输入傅氏谱、场地固有圆频率和场地等价阻尼比之间的物理关系,从随机傅氏谱函数角度描述了地震动随机过程的随机性本质。结合Ⅳ类工程场地的实测地震动记录资料,由数值方法识别了给出基本随机变量的概率分布参数。与实测记录对比表明,本文建立的随机地震动模型具有明确的物理概念,可充分反映地震动的变异性特征。     

Effect of earthquake ground motions on fragility curves of highway bridge piers based on numerical simulation

Fragility curves express the probability of structural damage due to earthquakes as a function of ground motion indices, e.g., PGA, PGV. Based on the actual damage data of highway bridges from the 1995 Hyogoken‐Nanbu (Kobe) earthquake, a set of empirical fragility curves was constructed. However, the type of structure, structural performance (static and dynamic) and variation of input ground motion were not considered to construct the empirical fragility curves. In this study, an analytical approach was adopted to construct fragility curves for highway bridge piers of specific bridges. A typical bridge structure was considered and its piers were designed according to the seismic design codes in Japan. Using the strong motion records from Japan and the United States, non‐linear dynamic response analyses were performed, and the damage indices for the bridge piers were obtained. Using the damage indices and ground motion indices, fragility curves for the bridge piers were constructed assuming a lognormal distribution. The analytical fragility curves were compared with the empirical ones. The proposed approach may be used in constructing the fragility curves for highway bridge structures. Copyright © 2001 John Wiley & Sons, Ltd.     

SEISMIC GROUND MOTION SYNTHETICS OF THE 1989 LOMA PRIETA EARTHQUAKE

A discrete wave number approach in conjunction with a propagator-based formalism is used to synthesize the Loma Prieta earthquake ground motion at both the near and the far field, taking into account all kinds of seismic waves (body and surface). A bilaterally propagating shear slip over a rectangular fault is used to describe the seismic source mechanism, while the earth model is based on geological profiles of the Santa Cruz mountain area and consists of three layers overlaying a half-space. The synthesized ground motion is first compared with actual records from the Loma Prieta earthquake and the agreement between the two is found to be satisfactory, as far as magnitude, duration and essential wave form characteristics are concerned. Then, ground motions are synthesized and plotted at a dense grid of observer locations over a large area around the epicenter, at different time instants. Using such plots, it is possible to study the generation ana propagation of different kinds of seismic waves, the spatial variability of ground motion, as well as the development of the permanent gound deformation.     

Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading

A method for parametric system identification of classically damped linear system in frequency domain is adopted and extended for non‐classically damped linear systems subjected up to six components of earthquake ground motions. This method is able to work in multi‐input/multi‐output (MIMO) case. The response of a two‐degree‐of‐freedom model with non‐classical damping, excited by one‐component earthquake ground motion, is simulated and used to verify the proposed system identification method in the single‐input/multi‐output case. Also, the records of a 10 storey real building during the Northridge earthquake is used to verify the proposed system identification method in the MIMO case. In this case, at first, a single‐input/multi‐output assumption is considered for the system and modal parameters are identified, then other components of earthquake ground motions are added, respectively, and the modal parameters are identified again. This procedure is repeated until all four components of earthquake ground motions which are measured at the base level of the building are included in the identification process. The results of identification of real building show that consideration of non‐classical damping and inclusion of the multi‐components effect of earthquake ground motions can improve the least‐squares match between the finite Fourier transforms of recorded and calculated acceleration responses. Copyright © 2006 John Wiley & Sons, Ltd.     

A predictive model for Arias intensity at multiple sites and consideration of spatial correlations

Arias intensity, I_a, has been identified as an efficient intensity measure for the estimation of earthquake‐induced losses. In this paper, a new model for the prediction of Arias intensity, which incorporates nonlinear site response through the use of the average shear‐wave velocity and a heteroskedastic variance structure, is proposed. In order to estimate the effects of ground motions on spatially‐distributed systems, it is important to take into account the spatial correlation of the intensity measure. However, existing loss‐estimation models, which use I_a as an input, do not take this aspect of the ground motion into account. Therefore, the potential to model the spatial correlation of Arias intensity is also investigated. The empirical predictive model is developed using recordings from the Pacific Earthquake Engineering Research Center Next Generation of Attenuation database whereas the model for spatial correlation makes use of the well‐recorded events from this database, that is the Northridge and Chi‐Chi earthquakes. Copyright © 2011 John Wiley & Sons, Ltd.     

Stochastic response of suspension bridges to earthquake forces

The purpose of the paper is to study the applicability of stochastic methods for determining the response in the vertical plane of long-span bridges to earthquakes. The bridges are of the modern type, with flexible towers, box-decks and inclined hangers, Humber and Bosporus being used as the main examples, although some useful information is obtained by studying a small suspension footbridge. The input ground acceleration is that of the Pacoima Dam record of the 1971 San Fernando event, but use is also made of filtered white noise. Because the stochastic approach is essentially an attempt to abstract usable information which could otherwise be obtained from a lengthy time-history approach, the procedure adopted here is to compare the same parameters obtained in these two ways, as far as possible. In particular, the maximum values of displacements and bending moments are considered. The essential question to be answered is whether the short length of earthquake records, coupled with the large natural periods of vibration of long-span bridges, allows sufficient response information to be generated to make statistical parameters meaningful. The short-span (50 m) footbridge gives no cause for concern here, but the detailed comparison with time-history solutions shows that the stochastic approach for both Bosporus and Humber has to be cautiously assessed, particularly if some trailing zeros are not added to the earthquake record. However, general conclusions are not made, because only one earthquake record has been used.     

Linear analysis of concrete arch dams including dam–water–foundation rock interaction considering spatially varying ground motions

The available substructure method and computer program for earthquake response analysis of arch dams, including the effects of dam–water–foundation rock interaction and recognizing the semi‐unbounded size of the foundation rock and fluid domains, are extended to consider spatial variations in ground motions around the canyon. The response of Mauvoisin Dam in Switzerland to spatially varying ground motion recorded during a small earthquake is analyzed to illustrate the results from this analysis procedure. Copyright © 2009 John Wiley & Sons, Ltd.