New method of generating spectrum compatible accelerograms using neural networks

A new method is proposed for generating artificial earthquake accelerograms from response spectra. This method uses the learning capabilities of neural networks to developed the knowledge of the inverse mapping from the response spectra to earthquake accelerogram. In the proposed method the neural networks learn the inverse mapping directly from the actual recorded earthquake accelerograms and their response spectra. A two-stage approach is used. In the first stage, a replicator neural network is used as a data compression tool. The replicator neural network compresses the vector of the discrete Fourier spectra of the accelerograms to vectors of much smaller dimension. In the second stage, a multi-layer feed-forward neural network learns to relate the response spectrum to the compressed Fourier spectrum. A simple example is presented, in which only 30 accelerograms are used to train the two-stage neural networks. This example demonstrates how the method works and shows its potential. © 1998 John Wiley & Sons, Ltd.     

Wavelet-based generation of spatially correlated accelerograms

For the seismic analysis of complex or nonlinear extended structures, it is useful to generate a set of properly correlated earthquake accelerograms that are consistent with a specified seismic hazard. A new simulation approach is presented in this paper for the generation of ensembles of spatially correlated accelerograms such that the simulated motions are consistent with (i) a parent accelerogram in the sense of temporal variations in frequency content, (ii) a design spectrum in the mean sense, and (iii) with a given instantaneous coherency structure. The formulation is based on the extension of stochastic decomposition technique to wavelet domain via the method of spectral factorization. A complex variant of the modified Littlewood-Paley wavelet function is proposed for the wavelet-based representation of earthquake accelerograms, such that this explicitly brings out the phase information of the signal, besides being able to decompose it into component time-histories having energy in non-overlapping frequency bands. The proposed approach is illustrated by generating ensembles of accelerograms at four stations.     

Enrichment of seismic ground motion databases using Karhunen–Loève expansion

Best estimate seismic analysis are generally based on time‐domain simulations of structural responses. The seismic load is then modeled by a stochastic process representing ground motion. For this purpose, the analyst can use recorded accelerograms or work with synthetically generated ones. The number of ground motion time‐histories available for a given scenario and site condition is limited and generally not sufficient for carrying out more advanced probabilistic structural response analysis. It is then necessary to have at our disposal methods that allow for generating synthetic accelerograms that realistically characterize earthquake ground motions. However, most of the methods proposed in literature for generating synthetic accelerograms do not accurately reproduce the natural variability of ground motion parameters (such as PGA, cumulative absolute velocity, and Arias intensity) observed for recorded time histories. In this paper, we introduce a new method for generating synthetic ground motion, based on Karhunen‐Loève decomposition and a non‐Gaussian stochastic model. The proposed method enables the structural analyst to simulate ground motion time histories featuring the properties mentioned above. To demonstrate its capability, we study the influence of the simulation method on different ground motion parameters and on soil response spectra. We finally compute fragility curves to illustrate the practical application of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of strength reduction factors via normalized pseudo‐acceleration response spectrum

Estimation of design forces in ductility‐based earthquake‐resistant design continues to be carried out with the application of response modification factors on elastic design spectra, and it remains interesting to explore how best to estimate strength reduction factors (SRFs) for a design situation. This paper considers the relatively less explored alternative of modelling SRF spectrum via a given response spectrum. A new model is proposed to estimate the SRF spectrum in terms of a pseudo‐spectral acceleration (PSA) spectrum and ductility demand ratio with the help of two coefficients. The proposed model is illustrated for an elasto‐plastic oscillator, in case of 10 recorded accelerograms and three ductility ratios. The proposed model is convenient and is able to predict SRF spectrum reasonably well, particularly at periods up to 1.0 s. Coefficients of the proposed model may also be determined in case of a given design spectrum when there is uncertainty in SRF spectrum due to uncertainty in temporal characteristics of the ground motion. This is illustrated with the help of 474 accelerograms recorded in western U.S.A. and different scaled PSA spectra. It is shown that probabilistic estimates may be obtained in this situation for SRF spectrum by assuming the error residuals to be log normally distributed with period‐dependent parameters. Copyright © 2006 John Wiley & Sons, Ltd.     

Generation of spectrum compatible accelerograms

In this paper a new method is presented for generating earthquake accelerograms which have pre-established response spectra. The non-stationary random nature and other salient features of the accelerograms can be taken care of by the procedure developed. The method leads to a sample spectrum which lies above the given spectrum. The generation of records to suit several spectra simultaneously can also be handled by this approach. The method is detailed first. This is followed by several numerical examples.     

根据渐进谱的统计规律生成地震加速度时程

基于美国西部80条基岩上的近场强震记录, 采用Nakayama方法生成记录的渐进谱, 并参照Kameda方式,用统计方法建立了根据震级、距离等地震参数预测渐进谱的统计模型. 提出一种以渐进谱为目标谱的生成幅值和频率非平稳地震加速度时程的迭代方法. 由于考虑了渐进谱幅值和相位的相互影响,所生成的时程的相位也是时频非平稳的,并在相位调整中识别了相位谱增量符号以加速迭代收敛进程. 最后根据统计回归的目标渐进谱模型和本文提出的拟合目标渐进谱的方法,可生成不同震级、距离条件下的幅值和频率均非平稳的地震加速度时程.      

Earthquake accelerogram simulation with statistical law of evolutionary power spectrum

By using the technique for evolutionary power spectrum proposed by Nakayama and with reference to the Kameda formula, an evolutionary spectrum prediction model for given earthquake magnitude and distance is established based on the 80 near-source acceleration records at rock surface with large magnitude from the ground motion database of western U.S.. Then a new iteration method is developed for generation of random accelerograms non-stationary both in amplitude and frequency which are compatible with target evolutionary spectrum. The phase spectra of those simulated accelerograms are also non-stationary in time and frequency domains since the interaction between amplitude and phase angle has been considered during the generation. Furthermore, the sign of the phase spectrum increment is identified to accelerate the iteration. With the proposed statistical model for predicting evolutionary power spectra and the new method for generating compatible time history, the artificial random earthquake accelerograms non-stationary both in amplitude and frequency for certain magnitude and distance can be provided.     

Numerical simulation of ground acceleration spectra and accelerograms for engineering application

A numerical approach to the earthquake ground motion analysis is proposed for regions where no accelerograms are available. Using Haskell matrix techniques, the response spectra of a layered substratum for SV waves were calculated and then multiplied by the spectra corresponding to Brune's type pulses. The ground acceleration spectra were obtained for different angles of pulse incidence at the substratum base. The spectrum shape depends upon the substratum response and the pulse shape, while its level was related to the maximum ground acceleration corresponding to the expected maximum intensity. Transformation of the ground spectra into the time domain produced numerical accelerograms for horizontal and vertical components and for different angles of pulse incidence. Finally, a standard statistical procedure was applied to obtain the design response spectra used in engineering applications.     

Simulation of three-dimensional strong ground motions along principal axes,San Fernando earthquake

Power spectral density which describes frequency content is considered one of the most significant properties to be taken into account when generating ground motions through the use of stochastic processes. Using a smoothed and normalized Fourier amplitude spectrum, frequency content for components of motion along a set of principal axes is estimated. Fourier amplitude spectra obtained by this moving-window technique are presented which show the time dependency of frequency content for motions produced by the San Fernando earthquake of 9 February 1971. A mathematical model to simulate ground motion processes is proposed for which both the intensity and frequency content are non-stationary. Using this mathematical model with parameter characteristics along principal axes similar to those of the motions recorded during the San Fernando earthquake, three-dimensional ground motions are synthetically generated. The properties of the simulated motions show general characteristics similar to the characteristics observed in real accelerograms. The suggested model is considered adequate for engineering purposes.     

Parametric modeling of earthquake response spectra

A mathematical model for the response spectra is determined using statistical analysis. The form of the model is first established using fifty computer simulated accelerograms. The final form is then used on twenty-five accelerograms from fifteen past United States earthquakes. This model smooths out peaks and valleys which are characteristic of the response spectrum of any single earthquake. Thus it serves as a ‘smooth design spectrum’ and can be used to approximate structural response to a future seismic event.     

Eurocode 8-compatible synthetic time-series as input to dynamic analysis

Nonlinear dynamic analysis of existing or planned structures often requires the use of accelerograms that match a target design spectrum. Here, our main concern is to generate a set of motions with a good level of fit to the Eurocode 8 design spectra for France. Synthetic time series are generated by means of a non-stationary stochastic method. To calibrate the input parameters in the stochastic approach, we select a reference set of accelerograms for a Eurocode 8 type B site category from the PEER Ground-Motion Database, which are then adjusted to the target spectrum through wavelet addition. Then, we compute nonlinear seismic responses of a soil column, including pore pressure effects, and brittle and ductile structures to the stochastic time-series, the natural accelerograms and time-series generated using stationary stochastic approaches. The results of these calculations reveal considerable variability in response despite the similarities in terms of spectral acceleration.     

An envelope for Mohr's circle in seismically excited three‐dimensional structures

A previously developed response‐spectrum‐based procedure for computing the envelope that bounds the time‐varying realizations of Mohr's circle at any prescribed location within a two‐dimensional structure is extended for use with three‐dimensional structures subjected to as many as three translational components of ground acceleration. The proposed envelope, which is completely defined by quantities that are routinely used and calculated in conventional response spectrum analyses, is developed for the general case in which the principal directions of the earthquake, along which the ground accelerations are uncorrelated, are unknown. The accuracy of the proposed envelope is evaluated by comparing it to the results of an ensemble of time‐history analyses performed on a concrete arch dam using simulated accelerograms. It is found that the proposed envelope has a level of accuracy that is suitable for structural design and analysis. The largest observed difference between the simulated and predicted mean envelopes is less than 5%. Copyright © 2004 John Wiley & Sons, Ltd.     

Scaling of strength reduction factors for degrading elasto‐plastic oscillators

The inelastic (design) spectra characterizing a seismic hazard are generally obtained by the scaling‐down of the elastic (design) spectra via a set of response modification factors. The component of these factors, which accounts for the ductility demand ratio, is known as the strength reduction factor (SRF), and the variation of this factor with initial period of the oscillator is called an SRF spectrum. This study considers scaling of the SRF spectrum in the case of an elasto‐plastic oscillator with strength and stiffness degradation characteristics. Two models are considered: one depending directly on the characterization of source and site parameters and the other depending on the normalized design spectrum characterization of the seismic hazard. The first model is the same as that proposed earlier by the second author, and is given in terms of earthquake magnitude, strong‐motion duration, predominant period, geological site conditions, ductility demand ratio, and ductility supply‐related parameter. The second model is a new model proposed here in terms of the normalized pseudo‐spectral acceleration values (to unit peak ground acceleration), ductility demand ratio and ductility supply‐related parameter. For each of these models, least‐square estimates of the coefficients are obtained through regression analyses of the data for 956 recorded accelerograms in western U.S.A. Parametric studies carried out with the help of these models confirm the dependence of SRFs on strong‐motion duration and earthquake magnitude besides predominant period and site conditions. It is also seen that degradation characteristics make a slight difference for high ductility demands and may lead to lower values of SRFs, unless the oscillators are very flexible. Copyright © 2004 John Wiley & Sons, Ltd.     

地震动随机函数模型研究(Ⅰ)——模型建立

基于随机过程的随机函数解释,发展了一种随机地震动模型.模型包括两个部分,分别是基于理论物理关系的Fourier幅值谱模型和基于经验关系的Fourier相位谱模型.基于拟加速度谱概念建立的幅值谱模型能够反映强震记录的幅值特性,而Fourier相位谱模型则利用累积相位谱的单值连续特性,以线性多项式与三角级数组合的形式给出.模拟记录与实测记录对比结果显示两者不仅在波形结构上相似性良好,而且反应谱层次亦吻合良好.     

Rocking strong earthquake accelerations

In this paper, the method presented by Lee and Trifunac (1985) for generating synthetic torsional accelerograms has been extended to the estimation of synthetic rocking accelerograms and of their response spectra. Results from our previous regression analyses for the characterization of strong shaking in terms of (1) earthquake magnitude and epicentral distance, or (2) Modified Mercalli Intensity at the site are utilized here again. The effects of geologic environment, in terms of site parameters or the representative depth of sediments, which influence amplification, and the dispersive properties of ground motion are also included. The synthetic rocking accelerogram is then constructed from the horizontal and vertical acceleration components.     

Wavelet‐based characterization of design ground motions

With the recent emergence of wavelet‐based procedures for stochastic analyses of linear and non‐linear structural systems subjected to earthquake ground motions, it has become necessary that seismic ground motion processes are characterized through statistical functionals of wavelet coefficients. While direct characterization in terms of earthquake and site parameters may have to wait for a few more years due to the complexity of the problem, this study attempts such characterization through commonly available Fourier and response spectra for design earthquake motions. Two approaches have been proposed for obtaining the spectrum‐compatible wavelet functionals, one for input Fourier spectrum and another for input response spectrum, such that the total number of input data points are 30–35% of those required for a time‐history analysis. The proposed methods provide for simulating ‘desired non‐stationary characteristics’ consistent with those in a recorded accelerogram. Numerical studies have been performed to illustrate the proposed approaches. Further, the wavelet functionals compatible with a USNRC spectrum in the case of 35 recorded motions of similar strong motion durations have been used to obtain the strength reduction factor spectra for elasto‐plastic oscillators and to show that about ±20% variation may be assumed from mean to 5 and 95% confidence levels due to uncertainty in the non‐stationary characteristics of the ground motion process. Copyright © 2002 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.     

适用于结构时程分析法的输入地震波

本文对建筑结构抗震设计时程法所需的输入地震波提出一种基于规范反应谱的、与设防烈度、场地条件、设计近震或远震、结构自振特性有关的选择原则和方法。按远、近震和四类场地标定了反应谱、延性谱和积累损伤谱。通过对六幢不同高度的剪切和弯曲型结构模型的弹塑性分析,表明离差很小。建议在进行结构时程法分析时,按地震加速度反应谱特定的分布规律选择4条加速度记录作为输入计算。两个实际的例子表明按上述方法计算的结果与按底部剪力法计算的结果基本相符。     

Wavelet-based generation of spectrum-compatible time-histories

This paper deals with the well-known problem of generating spectrum-compatible synthetic accelerograms for the linear and non-linear time-history analyses of structural systems. A wavelet-based procedure has been used to decompose a recorded accelerogram into a desired number of time-histories with non-overlapping frequency contents, and then each of the time-histories has been suitably scaled for matching of the response spectrum of the revised accelerogram with a specified design spectrum. The key idea behind this iterative procedure is to modify a recorded accelerogram such that the temporal variations in its frequency content are retained in the synthesized accelerogram. The proposed procedure has been illustrated by modifying five recorded accelerograms of widely different characteristics such that those are compatible with the same USNRC design spectrum.     

HORIZONTAL-TO-VERTICAL SPECTRUM RATIO OF EARTHQUAKE GROUND MOTION FOR SITE CHARACTERIZATION

Nakamura's method, which uses a horizontal-to-vertical Fourier spectrum ratio of microtremor, has become popular to determine the predominant period and amplification of a site. In this study, this method is extended for earthquake ground motion recordings using new strong motion data recorded by JMA-87-type accelerometers. From the analysis of these accelerograms, horizontal-to-vertical Fourier spectrum ratios of a site for different earthquakes are also found to be stable irrespective of magnitude, distance and depth. To establish this fact, attenuation relations of velocity response spectra for horizontal and vertical components are derived for three damping ratios (0, 2 and 5 per cent) using the JMA data. Then the horizontal-to-vertical ratios of the velocity response spectra are obtained. The results show that the horizontal and vertical velocity response spectra are dependent on magnitude, distance and depth, but that their ratios are almost independent of magnitude, distance and depth. However, since the current data set consists of mostly intermediate to far field data, this observation should be limited to records of these distance ranges. Introducing station coefficients, representing site amplification, to this relation yields the value comparable to the horizontal-to-vertical Fourier spectrum ratio at a specific site. The stability of the spectrum ratio is explained by the transfer function between the ground surface and stiff-soil outcrop due to S-wave propagation. These results suggest that site amplification characteristics can be evaluated by one-point two-component surface recordings of earthquake ground motion, in a similar manner as proposed by Nakamura for microtremor. © 1997 by John Wiley & Sons, Ltd.