Synthesis of accelerograms compatible with the Chinese GB 50011-2001 design spectrum via harmonic wavelets： artificial and historic records

A versatile approach is employed to generate artificial accelerograms which satisfy the compatibility criteria prescribed by the Chinese aseismic code provisions GB 50011-2001. In particular, a frequency dependent peak factor derived by means of appropriate Monte Carlo analyses is introduced to relate the GB 50011 -2001 design spectrum to a parametrically defined evolutionary power spectrum (EPS). Special attention is given to the definition of the frequency content of the EPS in order to accommodate the mathematical form of the aforementioned design spectrum. Further, a one-to-one relationship is established between the parameter controlling the time-varying intensity of the EPS and the effective strong ground motion duration. Subsequently, an efficient auto-regressive moving-average (ARMA) filtering technique is utilized to generate ensembles of non-stationary artificial accelerograms whose average response spectrum is in a close agreement with the considered design spectrum. Furthermore, a harmonic wavelet based iterative scheme is adopted to modify these artificial signals so that a close matching of the signals' response spectra with the GB 50011-2001 design spectrum is achieved on an individual basis. This is also done for field recorded accelerograms pertaining to the May, 2008 Wenchuan seismic event. In the process, zero-phase high-pass filtering is performed to accomplish proper baseline correction of the acquired spectrum compatible artificial and field accelerograms. Numerical results are given in a tabulated format to expedite their use in practice.     

基于小波变换的拟合规范反应谱多维地震动模拟

本文提出一种基于小波变换的拟合规范反应谱的多维地震动模拟算法。首先将规范反应谱推广到三维相关设计反应谱,然后将已有的三维地震动加速度时间历程曲线分解为一系列不同频段上的地震动分量,调整每一个地震动分量的幅值使其在相应的频率范围内拟合设计反应谱,最后经过调整后的地震动分量进行重构得到更新的地震动时间历程曲线。将该时间历程曲线的反应谱与目标反应谱进行比较,重复该过程直到误差位于特定的范围内。该方法可以保留原始地震动的局部时-频特性,为多维地震动的模拟提供了一种新的方法。     

Effective linear damping and stiffness coefficients of nonlinear systems for design spectrum based analysis

A stochastic approach for obtaining reliable estimates of the peak response of nonlinear systems to excitations specified via a design seismic spectrum is proposed. This is achieved in an efficient manner without resorting to numerical integration of the governing nonlinear equations of motion. First, a numerical scheme is utilized to derive a power spectrum which is compatible in a stochastic sense with a given design spectrum. This power spectrum is then treated as the excitation spectrum to determine effective damping and stiffness coefficients corresponding to an equivalent linear system (ELS) via a statistical linearization scheme. Further, the obtained coefficients are used in conjunction with the (linear) design spectrum to estimate the peak response of the original nonlinear systems. The cases of systems with piecewise linear stiffness nonlinearity, along with bilinear hysteretic systems are considered. The seismic severity is specified by the elastic design spectrum prescribed by the European aseismic code provisions (EC8). Monte Carlo simulations pertaining to an ensemble of nonstationary EC8 design spectrum compatible accelerograms are conducted to confirm that the average peak response of the nonlinear systems compare reasonably well with that of the ELS, within the known level of accuracy furnished by the statistical linearization method. In this manner, the proposed approach yields ELS which can replace the original nonlinear systems in carrying out computationally efficient analyses in the initial stages of the aseismic design of structures under severe seismic excitations specified in terms of a design spectrum.     

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

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

SMART-1台阵记录的长周期反应谱特性

利用台湾SMART－1台阵的数字强震仪加速度记录分析了同类场地上的地震动长周期分量特性，结果表明：同一次地震震中距基本相同的同类场地加速度反应谱短周期部分有较明显的差别，长周期部分则基本相同；震级和震中距对短周期反应谱影响明显，对长周期反应谱也有一定影响。最后将SMART－1台阳场地的平均谱与新的建筑抗震设计规范修订稿中提出的设计反应谱作了比较。     

BUILDING DUCTILITY DEMAND: INTERPLATE VERSUS INTRAPLATE EARTHQUAKES

The inelastic seismic response behaviour for a range of simplified single-degree-of-freedom models has been analysed using 180 random phase angle synthetic accelerograms with different frequency contents and different durations and 105 real accelerograms collected from different regions worldwide. Results from the analyses have identified that the frequency content of the excitation can greatly influence the ductility demand ratio due to inelastic amplification effects. Consequently, results derived from intraplate earthquake records (typically of higher frequency content) were generally different to those from interplate records. However, the commonly used El Centro accelerogram has significantly lower ductility demand in the low period range than the average of records with similar elastic response spectral shape. Apart from this, there was little evidence to suggest any inherent differences in the inelastic response behaviour of buildings from intraplate and interplate earthquakes which possessed similar frequency content. Thus, the average ductility demand ratios from future earthquakes in an area can be predicted by interpolation of the results presented in this paper assuming the elastic response spectrum has been defined. Ductility demand ratios derived from the synthetic accelerograms and the real accelerograms with similar frequency content have been shown to be consistent. However, results from synthetic records derived for the idealised code design spectra (such as the Uniform Building Code and the Australian Standard AS1170.4) indicate a significantly higher ductility demand in the long period range.     

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.     

Artificial ground motion compatible with specified ground shaking peaks and target response spectrum

1 Introduction Time history analysis is usually needed not only for the seismic design of important structures, such as nuclear plants, large concrete dams and super high-rise buildings, but also for seismic performance based design. Performance-based design requires tha knowledge of the entire process of nonlinear seismic responses of structures be obtained, which depends on the ground motion time history as input for this analysis. However, because the observed strong-motion recordings are l…     

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.     

Numbers of scaled and matched accelerograms required for inelastic dynamic analyses

Selecting, scaling and matching accelerograms are critically important to engineering design and assessment, enabling structural response to be determined with greater confidence and through fewer analyses than if unscaled accelerograms are employed. This paper considers the response of an 8‐storey multiple‐degree‐of‐freedom reinforced concrete structure to accelerograms selected, linearly scaled or spectrally matched using five different techniques. The first method consists of selecting real records on the basis of seismological characteristics, while the remaining methods make an initial selection on the basis of magnitude and spectral shape before (1) scaling to the target spectral acceleration at the initial period; (2) scaling to the target spectrum over a range of periods; (3) using wavelet adjustments to match the target spectrum and (4) using wavelet adjustments to match multiple target spectra for multiple damping ratios. The analyses indicate that the number of records required to obtain a stable estimate of the response decreases drastically as one moves through these methods. The exact number varies among damage measures and is related to the predictability of the damage measure. For measures such as peak roof and inter‐storey drift, member end rotation and the Park and Ang damage index, as few as one or two records are required to estimate the response to within ±5% (for a 64% confidence level) if matching to multiple damping ratios is conducted. Bias checks are made using predictive equations of the expected response derived from the results of 1656 nonlinear time‐domain analyses of the structure under the action of unscaled accelerograms. Copyright © 2008 John Wiley & Sons, Ltd.     

Prediction of elastic displacement response spectra in Europe and the Middle East

Empirical equations are presented for the prediction of displacement response ordinates for damping ratios of 2, 5, 10, 20 and 30% of critical and for response periods up to 4s, using 532 accelerograms from the strong‐motion databank from Europe and the Middle East. The records were all re‐processed and only employed for regressions at periods within the usable range, defined as a fraction of the filter cut‐off and depending on the instrument type (digital or analogue), earthquake magnitude and site class. The equations can be applied to predict the geometric mean displacement and pseudo‐acceleration spectra for earthquakes with moment magnitudes ( M ) between 5 and 7.6, and for distances up to 100km. The equations also include style‐of‐faulting and site class as explanatory variables. The predictions obtained from these new equations suggest that earlier European equations for spectral displacements underestimate the ordinates at longer periods as a result of severe filtering and the use of the spectral ordinates at periods too close to the filter cut‐off. The results also confirm that the period defining the start of the constant displacement plateau in the Eurocode 8 (EC8) spectrum is excessively short at 2s. The results not only show that the scaling factor defined in EC8 for estimating the spectral ordinates at damping ratios different from 5% of critical are a good general approximation, but also that this scaling varies with magnitude and distance (reflecting the influence of duration) and also displays a mild dependence on response period. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Selection of spectrum- and seismo-compatible accelerograms for the Tuscany region in Central Italy

This article illustrates the results of a study aimed at developing a methodology for the automatic identification of the seismic input at outcropping rock sites and flat topographic conditions necessary to carry out non-linear dynamic analysis of structures and geotechnical systems. The seismic input is provided in terms of a set of 7 natural accelerograms recorded on outcropping rock and satisfying the average spectral compatibility requirements prescribed by the Italian seismic code (NTC08).The study focuses on the territory encompassing Tuscany region in Central Italy and it has been carried out for six return periods, which are 50, 75, 101, 475, 712 and 949 years. The procedure involved four main steps: (1) grouping of the response spectra with similar features; (2) definition of the reference response spectrum for each group; (3) selection of spectrum-compatible accelerograms using the reference response spectrum of each group; and (4) linear scaling of the accelerograms to satisfy the compatibility requirement with respect to other response spectra of the group. The last step is implemented through an interactive, user-friendly program named SCALCONA 2.0, which provides the seismic input in agreement with the site location and return period specified by the user. The program is freely available at the following web site: http://www.rete.toscana.it/sett/pta/sismica/01informazione/banchedati/input_sismici/progettazione/index.htm.     

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

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

Generating multiple spectrum compatible accelerograms using stochastic neural networks

A new neural‐network‐based methodology for generating artificial earthquake spectrum compatible accelerograms from response spectra was proposed in 1997, in which, the learning capabilities of neural networks were used to develop the knowledge of the inverse mapping from the response spectra to earthquake accelerograms. Recently, this methodology has been further extended and enhanced. This paper presents a new stochastic neural network that is capable of generating multiple earthquake accelerograms from a single‐response spectrum. A new stochastic feature to the neural network has been combined with a new scheme for data compression using the replicator neural networks developed in the original method. A benefit of this extended methodology is gaining efficiency in compressing the earthquake accelerograms and extracting their characteristics. The proposed method produces a stochastic ensemble of earthquake accelerograms from any response spectra or design spectra. An example is presented that used 100 recorded accelerograms to train the neural network and several design spectra and response spectra to test this improved methodology. Copyright © 2001 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.     

A procedure for simulating synthetic accelerograms compatible with correlated and conditional probabilistic response spectra

The modeling of seismic load is a major topic that has to be addressed thoroughly in the framework of performance based seismic analysis and design. In this paper, a simple procedure for simulating artificial earthquake accelerograms matching the statistical distribution of response spectra, as given by median ground motion prediction equations, the standard deviation and correlation coefficients, is proposed. The approach follows the general ideas of the (natural) ground motion selection algorithms proposed by Baker [4] and Wang [43] but using simulated (artificial) “spectrum-compatible” accelerograms. This allows to simulate spectrum-compatible accelerograms featuring variability similar to the one of recorded accelerograms when the match of median and ±1 standard-deviation response spectra is imposed by the regulator. The procedure is illustrated by an application to the NGA ground motion data and models.