Conditional simulation of spatially variable seismic ground motions based on evolutionary spectra

This paper proposes a computational procedure for the conditional simulation of spatially variable seismic ground motions for long span bridges with multiple supports. The seismic ground motions, with part of their time histories measured at some supports, are regarded as zero‐mean nonstationary random processes characterized by predefined evolutionary power spectral density. To conditionally simulate unknown seismic ground motion time histories at other supports, the Kriging method is first described briefly for the conditional simulation of a random vector comprised of zero‐mean Gaussian variables. The multivariate oscillatory processes characterized by the evolutionary power spectral density matrix are then introduced, and the Fourier coefficients of the oscillatory processes and their covariance matrix are derived. By applying the Kriging method to the random vector of the Fourier coefficients and using the inverse Fourier transform, unknown nonstationary seismic ground motion time histories can be simulated. A numerical example is selected to demonstrate capabilities of the proposed simulation procedure, and the results show that the procedure can ensure unbiased time‐varying correlation functions, especially the cross correlation between known and unknown time histories. The procedure is finally applied to the Tsing Ma suspension bridge in Hong Kong to generate ground accelerations at its multiple supports using limited seismic records. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of orthogonal horizontal ground motion components for specified earthquake and site characteristics

A method for generating an ensemble of orthogonal horizontal ground motion components with correlated parameters for specified earthquake and site characteristics is presented. The method employs a parameterized stochastic model that is based on a time‐modulated filtered white‐noise process with the filter having time‐varying characteristics. Whereas the input white‐noise excitation describes the stochastic nature of the ground motion, the forms of the modulating function and the filter and their parameters characterize the evolutionary intensity and nonstationary frequency content of the ground motion. The stochastic model is fitted to a database of recorded horizontal ground motion component pairs that are rotated into their principal axes, a set of orthogonal axes along which the components are statistically uncorrelated. Model parameters are identified for each ground motion component in the database. Using these data, predictive equations are developed for the model parameters in terms of earthquake and site characteristics and correlation coefficients between parameters of the two components are estimated. Given a design scenario specified in terms of earthquake and site characteristics, the results of this study allow one to generate realizations of correlated model parameters and use them along with simulated white‐noise processes to generate synthetic pairs of horizontal ground motion components along the principal axes. The proposed simulation method does not require any seed recorded ground motion and is ideal for use in performance‐based earthquake engineering. Copyright © 2011 John Wiley & Sons, Ltd.     

地震波双模态时变修正Kanai-Tajimi非平稳随机模型的改进及参数识别

在Vlachos等提出的双模态时变修正Kanai-Tajimi功率谱模型及其参数识别方法的基础上,利用杜修力等提出的Kanai-Tajimi功率谱滤波方法并引进遗传算法及二次优化识别技术进行改进,建立地震动时变功率谱的参数模型化方法。通过集集地震波的时变功率谱模型参数识别及模拟地震动算例,验证改进后的双模态时变修正Kanai-Tajimi功率谱模型的可行性和有效性,其方法可运用到重大工程结构抗震分析的设计地震动输入中。     

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.     

Stochastic model for simulation of near‐fault ground motions

A parameterized stochastic model of near‐fault ground motion in two orthogonal horizontal directions is developed. The major characteristics of recorded near‐fault ground motions are represented. These include near‐fault effects of directivity and fling step; temporal and spectral non‐stationarity; intensity, duration, and frequency content characteristics; directionality of components; and the natural variability of ground motions. Not all near‐fault ground motions contain a forward directivity pulse, even when the conditions for such a pulse are favorable. The proposed model accounts for both pulse‐like and non‐pulse‐like cases. The model is fitted to recorded near‐fault ground motions by matching important characteristics, thus generating an ‘observed’ set of model parameters for different earthquake source and site characteristics. A method to generate and post‐process synthetic motions for specified model parameters is also presented. Synthetic ground motion time series are generated using fitted parameter values. They are compared with corresponding recorded motions to validate the proposed model and simulation procedure. The use of synthetic motions in addition to or in place of recorded motions is desirable in performance‐based earthquake engineering applications, particularly when recorded motions are scarce or when they are unavailable for a specified design scenario. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of synthetic ground motions for specified earthquake and site characteristics

A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.     

全非平稳地震作用的结构随机反应与可靠度分析

在平稳地震动过程的Clough-Penzien功率谱基础上,采用林家浩非均匀调制函数建立全非平稳地震动过程的演变功率谱。根据我国现行的《建筑抗震设计规范》进行全非平稳地震动演变功率谱的参数识别研究。应用非平稳随机过程模拟的谱表示-随机函数方法,生成建筑结构抗震设计所用地震动的代表性样本集合。通过代表性样本集合的二阶统计值及地震反应谱与目标值的拟合比较,验证本文方法的有效性。最后结合概率密度演化方法,进行以层间位移角为控制准则的结构随机地震反应分析与抗震可靠度计算。     

地震地面运动模型的研究

地震地面运动被模拟成均值为零的两次过滤Ｇａｕｓｓ白噪声随机过程。第一次过滤削减白噪声的高频含量;第二次过滤削减白噪声的低频含量。根据地震记录的功率谱,使用非线性函数的最小二乘法,确定了两次过滤Ｇａｕｓｓ白噪声随机过程的功率谱密度函数的参数。     

基于强度和能量控制的随机地震动模型

考虑地震动的随机性和频率与强度非平稳性,通过理论分析,提出了一般随机地震动模型,并给出了确定模型参数的原则和方法。该模型以地震动强度、地震动能量以及地震动持时等宏观指标作为控制随机地震动模型参数的指标,而对其内在的频谱组成等指标只要求满足一般地震动的特征。该模型可以用于描述平稳随机过程、强度非平稳随机过程以及强度和频率完全非平稳随机过程。通过与常用功率谱模型的比较,验证了该模型的合理性。     

Wavelet-Hilbert transform-based simulation of pulse-like ground motion

In nonlinear dynamic structural analysis, a suite of pulse-like ground motions is required for the performance-based design of structures near active faults. The dissimilarity in the amplitude and frequency content of the earthquake time series referred to nonstationary properties in temporal and spectral, respectively. An approach is proposed based on the nonstationary properties of the far-field records and the seismological information in an event for simulating pulse-like records. The pulse-like earthquake time history is estimated via the superposition of the residual part of the earthquake with the estimated pulse. The wavelet-based Hilbert transform is utilized to characterize the nonstationary properties, the instantaneous amplitude, and frequencies of far-field records to model residual part. The effects of near-fault and pulse are estimated based on the seismological properties of the region. The validation of the procedure is indicated by comparing simulated time-series, response spectra, and Arias intensity with recorded pulse-like records in two different earthquakes in California; the M_w 6.7 1994 Northridge and the M_w 6.5 1979 Imperial valley.

     

An evolutionary stochastic ground-motion model defined by a seismological scenario and local site conditions

This paper is concerned with modeling earthquake-induced ground accelerations and the simulation of the dynamic response of linear structures through the principles of stochastic dynamics. A fully evolutionary approach, with nonstationarity both in amplitude and in frequency content, is proposed in order to define the seismic action, based on seismological information in the form of a small number of input parameters commonly available in deterministic or probabilistic seismic design situations. The signal is obtained by filtering a Gaussian white-noise. The finite duration and time-varying amplitude properties are obtained by using a suitable envelope function. By utilizing a subset of the records from the PEER-NGA strong-motion database, and time-series analysis tools extended to nonstationary processes, the key transfer-function properties, in terms of circular frequency, damping ratio and spectral intensity factor, are identified. A regression analysis is conducted for practical and flexible application of this model, in order to empirically relate the identified time-varying parameters of the filter to the characteristics defining earthquake scenarios such as magnitude, rupture distance and soil type. A validation study and a parametric investigation using elastic response spectra is also included. Results show that the final seismic model can reproduce, with satisfactory accuracy, the characteristics of acceleration records in a region, over a broad range of response periods.     

A numerical study on the cumulative out‐of‐plane damage to church masonry façades due to a sequence of strong ground motions

Seismic shocks occur sometimes as a sequence, close in space and time, of destructive events of comparable intensity. In these cases, a significant portion of the damage to historical buildings can be related with the cumulated damage on structures that become progressively more vulnerable. This research investigates the specific increase of damage determined by a sequence of strong ground motions, focusing the interest on the out‐of‐plane response of 2 church masonry façades. The dynamic analyses were performed by a specific rigid body and spring model RBSM, which only accounts for out‐of‐plane damage mechanisms. Two idealized models of façade, each made of 2 different masonry bonds, have been studied by applying various sequences of recorded accelerograms. The results highlighted a complex relationship between the spectral content of the seismic shocks and the characteristics of the structures that change in the course of the loading sequence due to the development of damage. The Housner spectral intensity proved to be a reliable scalar measure of the ground motion destructiveness for these façades. Moreover, when considering a design‐consistent accelerogram that causes a relevant damage pattern, ie, with a significant elongation of the effective first period of vibration, the numerical results indicated a possible spectral intensity threshold below which the occurrence of repeated seismic shocks, both before and after the reference design shock, can be considered as irrelevant. On the other hand, a catastrophic increase of damage should be expected when this threshold is overcome.     

Ductility‐dependent intensity measure that accounts for ground‐motion spectral shape and duration

The calculated nonlinear structural responses of a building can vary greatly, even if recorded ground motions are scaled to the same spectral acceleration at a building's fundamental period. To reduce the variation in structural response at a particular ground‐motion intensity, this paper proposes an intensity measure (IM_comb) that accounts for the combined effects of spectral acceleration, ground‐motion duration, and response spectrum shape. The intensity measure includes a new measure of spectral shape that integrates the spectrum over a period range that depends on the structure's ductility. The new IM is efficient, sufficient, scalable, transparent, and versatile. These features make it suitable for evaluating the intensities of measured and simulated ground motions. The efficiency and sufficiency of the new IM is demonstrated for the following: (i) elastic‐perfectly plastic single‐degree‐of‐freedom (SDOF) oscillators with a variety of ductility demands and periods; (ii) ductile and brittle deteriorating SDOF systems with a variety of periods; and (iii) collapse analysis for 30 previously designed frames. The efficiency is attributable to the inclusion of duration and to the ductility dependence of the spectral shape measure. For each of these systems, the transparency of the intensity measure made it possible to identify the sensitivity of structural response to the various characteristics of the ground motion. Spectral shape affected all structures, but in particular, ductile structures. Duration only affected structures with cyclic deterioration. Copyright © 2015 John Wiley & Sons, Ltd.     

Energy‐compatible and spectrum‐compatible (ECSC) ground motion simulation using wavelet packets

A stochastic ground‐motion simulation and modification technique is developed to generate energy‐compatible and spectrum‐compatible (ECSC) synthetic motions through wavelet packet characterization and modification in both frequency and time domains. The ECSC method significantly advances traditional spectral matching approaches, because it generates ground motions that not only match the target spectral accelerations, but also match Arias intensity build‐up and significant durations. The great similarity between the ECSC simulated motions and the actual recorded motions is demonstrated through one‐to‐one comparison of a variety of intensity measures. Extensive numerical simulations were also performed to validate the performance of the ECSC ground motions through nonlinear analyses of elasto‐plastic oscillators. The ECSC method can be easily implemented in the generalized conditional intensity measure framework by directly simulating a set of motions following a targeted distribution of multiple intensity measures. Therefore, the ECSC method has great potential to be used in performance‐based earthquake design and analysis. Copyright © 2017 John Wiley & Sons, Ltd.     

Nonstationary seismic responses of structure with nonlinear stiffness subject to modulated Kanai–Tajimi excitation

A simple structure under earthquake excitation is modeled as a single‐degree‐of‐freedom system with nonlinear stiffness subject to modulated Kanai–Tajimi excitation. The nonstationary responses including the nonstationary probability densities of the system responses and the statistical moments are obtained in semi‐analytical form. By applying the stochastic averaging method based on the generalized harmonic functions, the averaged Fokker–Planck–Kolmogorov(FPK) equation governing the nonstationary probability density of the amplitude is derived. Then, the solution of the FPK equation is approximately expressed by a series expansion in terms of a set of properly selected basis functions with time‐dependent coefficients. According to the Galerkin method, the time‐dependent coefficients are solved from a set of linear first‐order differential equations. Thus, the nonstationary probability densities of the amplitude and the state responses as well as the statistic moments of the amplitude are obtained. Finally, two types of the modulating functions, i.e. constant function and exponential function, are considered to give some semi‐analytical formulae. The proposed procedures are checked against the Monte Carlo simulation. The effects of the structure natural frequency and the intensity of the excitation as well as the ground stiffness on the system responses are discussed. It should be pointed out that the proposed method is good for broadband excitation and light damping. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling spatially correlated spectral accelerations at multiple periods using principal component analysis and geostatistics

Regional seismic risk assessments and quantification of portfolio losses often require simulation of spatially distributed ground motions at multiple intensity measures. For a given earthquake, distributed ground motions are characterized by spatial correlation and correlation between different intensity measures, known as cross‐correlation. This study proposes a new spatial cross‐correlation model for within‐event spectral acceleration residuals that uses a combination of principal component analysis (PCA) and geostatistics. Records from 45 earthquakes are used to investigate earthquake‐to‐earthquake trends in application of PCA to spectral acceleration residuals. Based on the findings, PCA is used to determine coefficients that linearly transform cross‐correlated residuals to independent principal components. Nested semivariogram models are then fit to empirical semivariograms to quantify the spatial correlation of principal components. The resultant PCA spatial cross‐correlation model is shown to be accurate and computationally efficient. A step‐by‐step procedure and an example are presented to illustrate the use of the predictive model for rapid simulation of spatially cross‐correlated spectral accelerations at multiple periods.     

A spatial cross‐correlation model of spectral accelerations at multiple periods

Many seismic loss problems (such as disruption of distributed infrastructure and losses to portfolios of structures) are dependent upon the regional distribution of ground‐motion intensity, rather than intensity at only a single site. Quantifying ground‐motion over a spatially‐distributed region therefore requires information on the correlation between the ground‐motion intensities at different sites during a single event. The focus of the present study is to assess the spatial correlation between ground‐motion spectral accelerations at different periods. Ground motions from eight well‐recorded earthquakes were used to study the spatial correlations. On the basis of obtained empirical correlation estimates, we propose a geostatistics‐based method to formulate a predictive model that is suitable for simulation of spectral accelerations at multiple sites and multiple periods, in the case of crustal earthquakes in active seismic regions. While the calibration of this model and investigation of its implications were somewhat complex, the model itself is very simple to use for making correlation predictions. A user only needs to evaluate a simple equation relying on three sets of coefficients provided here to compute a correlation coefficient for spectral values at two periods and at a specified separation distance. These results may then be used in evaluating the seismic risk of portfolios of structures with differing fundamental periods. Copyright © 2012 John Wiley & Sons, Ltd.     

Developing efficient scalar and vector intensity measures for IDA capacity estimation by incorporating elastic spectral shape information

Scalar and vector intensity measures are developed for the efficient estimation of limit‐state capacities through incremental dynamic analysis (IDA) by exploiting the elastic spectral shape of individual records. IDA is a powerful analysis method that involves subjecting a structural model to several ground motion records, each scaled to multiple levels of intensity (measured by the intensity measure or IM), thus producing curves of structural response parameterized by the IM on top of which limit‐states can be defined and corresponding capacities can be calculated. When traditional IMs are used, such as the peak ground acceleration or the first‐mode spectral acceleration, the IM‐values of the capacities can display large record‐to‐record variability, forcing the use of many records to achieve reliable results. By using single optimal spectral values as well as vectors and scalar combinations of them on three multistorey buildings significant dispersion reductions are realized. Furthermore, IDA is extended to vector IMs, resulting in intricate fractile IDA surfaces. The results reveal the most influential spectral regions/periods for each limit‐state and building, illustrating the evolution of such periods as the seismic intensity and the structural response increase towards global collapse. The ordinates of the elastic spectrum and the spectral shape of each individual record are found to significantly influence the seismic performance and they are shown to provide promising candidates for highly efficient IMs. Copyright © 2005 John Wiley & Sons, Ltd.     

Modification of stochastic ground motion models for matching target intensity measures

Stochastic ground motion models produce synthetic time‐histories by modulating a white noise sequence through functions that address spectral and temporal properties of the excitation. The resultant ground motions can be then used in simulation‐based seismic risk assessment applications. This is established by relating the parameters of the aforementioned functions to earthquake and site characteristics through predictive relationships. An important concern related to the use of these models is the fact that through current approaches in selecting these predictive relationships, compatibility to the seismic hazard is not guaranteed. This work offers a computationally efficient framework for the modification of stochastic ground motion models to match target intensity measures (IMs) for a specific site and structure of interest. This is set as an optimization problem with a dual objective. The first objective minimizes the discrepancy between the target IMs and the predictions established through the stochastic ground motion model for a chosen earthquake scenario. The second objective constraints the deviation from the model characteristics suggested by existing predictive relationships, guaranteeing that the resultant ground motions not only match the target IMs but are also compatible with regional trends. A framework leveraging kriging surrogate modeling is formulated for performing the resultant multi‐objective optimization, and different computational aspects related to this optimization are discussed in detail. The illustrative implementation shows that the proposed framework can provide ground motions with high compatibility to target IMs with small only deviation from existing predictive relationships and discusses approaches for selecting a final compromise between these two competing objectives.