Wavelet‐based non‐stationary response analysis of a friction base‐isolated structure

A wavelet‐based stochastic formulation has been presented in this paper for the seismic analysis of a base‐isolated structural system which is modelled as a two‐degree‐of‐freedom (2‐DOF) system. The ground motion has been modelled as a non‐stationary process (both in amplitude and frequency) by using modified Littlewood–Paley basis wavelets. The proposed formulation is based on replacing the non‐linear system by an equivalent linear system with time‐dependent damping properties. The expressions of the instantaneous damping and the power spectral density function (PSDF) of the superstructure response have been obtained in terms of the functionals of input wavelet coefficients. The proposed formulation has been validated by simulating a ground motion process. The effect of the frequency non‐stationarity on the non‐linear response has also been studied in detail, and it has been clearly shown how ignoring the frequency non‐stationarity in the ground motion leads to inaccurate non‐linear response calculations. Copyright © 2000 John Wiley & Sons, Ltd.     

Dynamic soil–structure interaction effects on the seismic response of suspension bridges

A stochastic approach has been formulated for the linear analysis of suspension bridges subjected to earthquake excitations. The transfer functions of various responses have been formulated while including the effects of dynamic Soil–Structure Interaction (SSI) via the use of the fixed-base modes of the structure. The excitation has been characterized by the ‘equivalent stationary’ processes corresponding to the free-field motions at each support and by an assumed coherency function between these motions. The proposed formulation considers the non-stationarity in the structural response due to sudden application of excitation by considering (i) the time-dependent frequency response functions, and (ii) the order statistics formulation for the peak factors in evolutionary response processes. The formulation has been illustrated by analysing the seismic response of the Golden Gate Bridge at San Francisco for two example excitations conforming to USNRC-specified design spectra. The significance of various governing parameters on the dynamic soil–structure interaction effects on the seismic response of suspension bridges has also been studied. It has been found that the contribution of the vertical component of ground motion to the bridge response increases with increasing soil compliance. Also, the extent to which the spatial variation of ground motion affects the bridge response depends on how significant the SSI effects are. Copyright © 1999 John Wiley & Sons Ltd.     

Stochastic seismic response of multiply-supported secondary systems in flexible-base structures

A formulation has been proposed for the transfer function of a secondary system response while the primary system is supported on a compliant soil and the excitation comprises of translational ground motion at its base. For this purpose, the earlier formulation of the authors for the fixed-base case, which exactly considers the interaction between the two sub-systems and is based on the use of their individual modal properties, has been extended. Also, the concept of modifying the input excitation for the interaction accelerations (associated with the soil–structure interaction) has been used. An example P–S system and three example earthquake excitations have been considered to illustrate the proposed formulation and to estimate the expected response peak amplitudes in the secondary system. This study shows that ‘detuning’ of the tuned systems may occur in case of significant soil–structure interaction. Further, for the reasons of both safety and economy, ignoring the interaction effects in designing the secondary systems may not always be justified. Copyright © 1999 John Wiley & Sons, Ltd.     

A damage-based definition of effective peak acceleration

This paper presents a rational basis for obtaining the Effective Peak Acceleration (EPA) of a given ground motion process. The proposed formulation considers the statistical variability in the ground motion, and is centred on the idea of explicitly linking EPA with expected cumulative damage in the structures due to the inelastic excursions. The structural behaviour has been modelled by a Single-Degree-Of-Freedom (SDOF) bilinear hysteretic oscillator. EPA is considered to be the expected PGA of a scaled ground motion process such that this oscillator undergoes a specified expected damage under the unscaled process if it is linearly designed for the scaled process. For estimation of the damage, the oscillator has been replaced by an equivalent linear oscillator through stochastic averaging. A parametric study has been carried out to investigate the dependence of EPA on several governing parameters, and it has been shown that despite the strong dependence of EPA on oscillator time period, it may be possible to obtain ‘period-averaged’ EPA values for several ground motion processes for engineering applications. © 1998 John Wiley & Sons, Ltd.     

Prediction of the median IDA curve by employing a limited number of ground motion records

A methodology has been proposed which can be used to reduce the number of ground motion records needed for the reliable prediction of the median seismic response of structures by means of incremental dynamic analysis (IDA). This methodology is presently limited to predictions of the median IDA curve only. The reduction in the number of ground motion records needed to predict the median IDA curve is achieved by introducing a precedence list of ground motion records. The determination of such a list is an optimization problem, which is solved in the paper by means of (1) a genetic algorithm and (2) a proposed simple procedure. The seismic response of a simple, computationally non‐demanding structural model has been used as input data for the optimization problem. The presented example is a three‐storey‐reinforced concrete building, subjected to two sets of ground motion records, one a free‐field set and the other a near‐field set. It is shown that the median IDA curves can be predicted with acceptable accuracy by employing only four ground motion records instead of the 24 or 30 records, which are the total number of ground motion records for the free‐field and near‐field sets, respectively. Copyright © 2007 John Wiley & Sons, Ltd.     

Kriging metamodeling in seismic risk assessment based on stochastic ground motion models

An efficient computational framework is presented for seismic risk assessment within a modeling approach that utilizes stochastic ground motion models to describe the seismic hazard. The framework is based on the use of a kriging surrogate model (metamodel) to provide an approximate relationship between the structural response and the structural and ground motion parameters that are considered as uncertain. The stochastic character of the excitation is addressed by assuming that under the influence of the white noise (used within the ground motion model) the response follows a lognormal distribution. Once the surrogate model is established, a task that involves the formulation of an initial database to inform the metamodel development, it is then directly used for all response evaluations required to estimate seismic risk. The model prediction error stemming from the metamodel is directly incorporated within the seismic risk quantification and assessment, whereas an adaptive approach is developed to refine the database that informs the metamodel development. The ability to efficiently obtain derivative information through the kriging metamodel and its utility for various tasks within the probabilistic seismic risk assessment is also discussed. As an illustrative example, the assessment of seismic risk for a benchmark four‐story concrete office building is presented. The potential that ground motions include near‐fault characteristics is explicitly addressed within the context of this example. The implementation of the framework for the same structure equipped with fluid viscous dampers is also demonstrated. Copyright © 2015 John Wiley & Sons, Ltd.     

逆断层地震作用下三维沉积盆地地震动谱元法模拟

近断层地震动对近场区域会造成严重的地震灾害,尤其区域内存在沉积地形时会进一步加重灾害。目前,针对释放能量更大、破坏力更强的逆冲断层地震作用下沉积盆地的地震动响应研究还未见报道。本文即采用谱元法,研究了动力学逆断层地震作用下的三维沉积盆地的动力响应。文中以椭球形沉积盆地为例,对其逆断层地震动响应特性进行了分析,并探讨了改变内外介质波速比和沉积厚度时沉积盆地内部观测点地震动时程和峰值变化规律。研究表明:①沉积内外介质波速比对沉积盆地的地震动影响显著,当沉积内部介质波速比降低时,盆地内部地表的峰值响应增大,地震动持续时间明显延长,尤其是位于盆地中心加速度峰值放大2.08倍,持时延长1.97倍;②沉积盆地厚度同样对其地震动响应产生影响,当沉积厚度增加时盆地中心位置地震动响应减小,加速度峰值缩小约0.64倍,而盆地边缘区域的地震动响应明显增大,峰值放大约1.35倍。      

三维地震断层动力破裂的显式并行有限元解法

本文在近场波动有限元的基础上,根据断层动力破裂的摩擦模型,详细地推导出了断层面节点在破裂的不同时刻的计算公式,提出了一种模拟断层动力破裂的显式并行有限元方法,利用此方法,可以模拟断层的破裂过程、地表破裂及由此产生的地面运动。     

Making sense of Australian earthquakes: the Ambraseys legacy

Through his students, Professor Nick Ambraseys has had a strong impact on the introduction of earthquake engineering practices in Australia, including historical earthquake studies, strong motion instrumentation and analysis, foundation studies including liquefaction, and building code formulation. In Australia the process of upgrading the third and current edition of the earthquake code and hazard map has begun. There are now about 150 digital strong motion recorders installed in cities, on major structures and at Australian National Seismograph Network sites. Three volumes of an isoseismal atlas have been published totalling more than 300 maps, mainly historical earthquakes. Significant progress has been made in paleoseismological studies across the continent, adding to the complexity of the intraplate seismicity model. With time and more installed accelerographs, the peak ground acceleration recorded in Australia has increased from .25 g in 1984 to .5 g in 1988 and 1 g in 1994, all from earthquakes smaller than magnitude 5, supporting Ambraseys contention that PGA alone is not a suitable parameter for a design ground motion.     

Viscous damping formulation and high frequency motion propagation in non-linear site response analysis

Non-linear time domain site response analysis is widely used in evaluating local soil effects on propagated ground motion. This approach has generally provided good estimates of field behavior at longer periods but has shortcomings at relatively shorter periods. Viscous damping is commonly employed in the equation of motion to capture damping at very small strains and employs an approximation of Rayleigh damping using the first natural mode only. This paper introduces a new formulation for the viscous damping using the full Rayleigh damping. The new formulation represents more accurately wave propagation for soil columns greater than 50 m thick and improves non-linear site response analysis at shorter periods. The proposed formulation allows the use of frequency dependent viscous damping. Several examples, including a field case history at Treasure Island, California, demonstrate the significant improvement in computed surface response using the new formulation.     

Ground motion selection for simulation‐based seismic hazard and structural reliability assessment

This paper examines four methods by which ground motions can be selected for dynamic seismic response analyses of engineered systems when the underlying seismic hazard is quantified via ground motion simulation rather than empirical ground motion prediction equations. Even with simulation‐based seismic hazard, a ground motion selection process is still required in order to extract a small number of time series from the much larger set developed as part of the hazard calculation. Four specific methods are presented for ground motion selection from simulation‐based seismic hazard analyses, and pros and cons of each are discussed via a simple and reproducible illustrative example. One of the four methods (method 1 ‘direct analysis’) provides a ‘benchmark’ result (i.e., using all simulated ground motions), enabling the consistency of the other three more efficient selection methods to be addressed. Method 2 (‘stratified sampling’) is a relatively simple way to achieve a significant reduction in the number of ground motions required through selecting subsets of ground motions binned based on an intensity measure, IM. Method 3 (‘simple multiple stripes’) has the benefit of being consistent with conventional seismic assessment practice using as‐recorded ground motions, but both methods 2 and 3 are strongly dependent on the efficiency of the conditioning IM to predict the seismic responses of interest. Method 4 (‘generalized conditional intensity measure‐based selection’) is consistent with ‘advanced’ selection methods used for as‐recorded ground motions and selects subsets of ground motions based on multiple IMs, thus overcoming this limitation in methods 2 and 3. Copyright © 2015 John Wiley & Sons, Ltd.     

基于非平稳随机地震动模型的重力坝抗震可靠度分析

基于正交展开的非平稳随机地震动模型,并考虑混凝土材料的非线性和坝体与库水之间的流固耦合,对印度Koyna重力坝进行有限元分析,得到坝顶水平位移和坝颈拉应力,结合概率密度演化方法和等价极值事件的思想,获得丰富的概率信息。这为坝体结构的随机地震反应分析和可靠度研究提供新的途径。     

Evaluation of dynamic eccentricities obtained using a probabilistic approach,response spectrum methods and modern building codes

In this paper a probabilistic approach has been adopted to study both the effects of uncertainty in earthquake frequency content and the correlation between earthquake frequency content and ground motion intensity on the response of a single-storey torsionally coupled elastic structure. The earthquake ground motion has been assumed to be a Gaussian, zero mean, stationary random process which is fully characterized by a power spectrum. The ground acceleration power spectrum is idealized as a probabilistic normalized power spectrum computed from actual earthquake records. The advantage of such an idealization is that it enables the effect of the natural frequency as a controlling structural parameter in torsional coupling to be assessed. Comparisons of the dynamic amplifications of eccentricity with those obtained from modern codes of practice and conventional response spectrum analyses have been made. The results of this study have shown that the variation in the frequency content has a significant effect on the response of low frequency structures, while the correlation between the frequency content and the intensity of seismic ground motion is insignificant for the wide range of structures considered. The structure natural frequency has been shown to be an important controlling parameter in the torsionally coupled response of structures subject to seismic loading. The frequency dependence of the dynamic amplification of eccentricity was found not to be reflected in the response spectrum analysis and the torsional provisions of modern building codes.     

非平稳人造地震动反应谱拟合的震记录应用研究

利用天然地震记录数据提取其非平稳相位数据的方法、论证人造地震动方法的可行性和合理性,并以呼和浩特市基岩场地非平稳人造地震动预测拟合应用为例,进行人造地震动反应谱拟合.本研究方法明显加强地震监测数字化记录与结构地震动输入预测模拟的密切联系,使地震动预测因充分利用近场地震观测资料而得到质的飞跃.     

A simplified fragility analysis of fan type cable stayed bridges

A simplified fragility analysis of fan type cable stayed bridges using Probabilistic Risk Analysis (PRA) procedure is presented for determining their failure probability under random ground motion. Seismic input to the bridge support is considered to be a risk consistent response spectrum which is obtained from a separate analysis. For the response analysis, the bridge deck is modeled as a beam supported on springs at different points. The stiffnesses of the springs are determined by a separate 2D static analysis of cable-tower-deck system. The analysis provides a coupled stiffness matrix for the spring system. A continuum method of analysis using dynamic stiffness is used to determine the dynamic properties of the bridges .The response of the bridge deck is obtained by the response spectrum method of analysis as applied to multidegree of freedom system which duly takes into account the quasi - static component of bridge deck vibration. The fragility analysis includes uncertainties arising due to the variation in ground motion, material property, modeling, method of analysis, ductility factor and damage concentration effect. Probability of failure of the bridge deck is determined by the First Order Second Moment (FOSM) method of reliability. A three span double plane symmetrical fan type cable stayed bridge of total span 689 m, is used as an illustrative example. The fragility curves for the bridge deck failure are obtained under a number of parametric variations. Some of the important conclusions of the study indicate that (i) not only vertical component but also the horizontal component of ground motion has considerable effect on the probability of failure; (ii) ground motion with no time lag between support excitations provides a smaller probability of failure as compared to ground motion with very large time lag between support excitation; and (iii) probability of failure may considerably increase for soft soil condition.     

Variability of tall building response to earthquakes with changing epicentre direction

Using a 40-storey building as an example, the structural response under the excitation of horizontal earthquake ground motion is shown to vary greatly with the direction from the building site to the epicentre. Two mathematical models are used to represent the ground motion: the stationary random process model and the evolutionary-type non-stationary random process model. The former model is suitable if the duration of the most intense portion of an earthquake is much longer than the fundamental natural period of the structure. The maximum standard deviations of selected structural response quantities are computed and plotted for all possible epicentre directions, and the information is believed to be useful for design purposes.     

考虑场地类别和断层距的地震动及结构响应参数相关性分析

为了研究场地条件和断层距对地震动参数及其与结构响应参数相关性的影响,本文基于NGA数据库中5 266条水平分量地震动记录,采用皮尔森相关系数,分析了不同场地条件和断层距下地震动参数之间的相关性及其与结构响应参数的关系。结果表明,场地条件和断层距对绝大部分地震动参数相关性会产生明显影响,但存在一些参数的相关性几乎不受场地条件和断层距的影响;以四层钢筋混凝土框架结构为例,发现在不同场地条件和断层距下,结构响应参数与地震动参数的相关性变化明显,这表明在研究结构响应与地震动参数相关性时需要考虑场地条件与断层距的影响。     

Rotation-invariant measures of earthquake response spectra

A new procedure for combining the response spectra of two horizontal components of recorded ground motion is presented. The presented formulation accounts for different orientations of accelerometer sensors and derives the maximum and the expected (mean) horizontal response spectra at a site, both of which are invariant to rotation of sensor axes. The maximum response spectrum is derived as the peak resultant response of single degree of freedom oscillators subjected to the as-recorded ground acceleration. The expected spectrum is derived by projecting the displacement response (due to as-recorded motion) along two orthogonal axes to a principal axes in which the displacement responses are uncorrelated. This property is used to formulate an approximation for the expected response spectrum over all possible sensor orientations. A large set of accelerometric data from Europe and the Middle East is used to demonstrate the applicability of the proposed response spectral measures.     

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ)

Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and El Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.