Conditional spectrum‐based ground motion record selection using average spectral acceleration

The use of a seismic intensity measure (IM) is paramount in decoupling seismic hazard and structural response estimation when assessing the performance of structures. For this to be valid, the IM needs to be sufficient;that is, the engineering demand parameter (EDP) response should be independent of other ground motion characteristics when conditioned on the IM. Whenever non‐trivial dependence is found, such as in the case of the IM being the first‐mode spectral acceleration, ground motion selection must be employed to generate sets of ground motion records that are consistent vis‐à‐vis the hazard conditioned on the IM. Conditional spectrum record selection is such a method for choosing records that are consistent with the site‐dependent spectral shape conditioned on the first‐mode spectral acceleration. Based on a single structural period, however the result may be suboptimal, or insufficient, for EDPs influenced by different period values, for example, peak interstory drifts or peak floor accelerations at different floors, potentially requiring different record suites for each. Recently, the log‐average spectral acceleration over a period range, AvgSA, has emerged as an improved scalar IM for building response estimation whose hazard can be evaluated using existing ground motion prediction equations. Herein, we present a recasting of conditional spectrum record selection that is based on AvgSA over a period range as the conditioning IM. This procedure ensures increased efficiency and sufficiency in simultaneously estimating multiple EDPs by means of a single IM. Copyright © 2017 John Wiley & Sons, Ltd.     

A vector-valued intensity measure for near-fault ground motions

Near-fault ground motions containing high energy and large amplitude velocity pulses may cause severe damage to structures. The most widely used intensity measure (IM) is the elastic spectral acceleration at the fundamental period of the structure (Sa(T₁)); however, Sa(T₁) is not a sufficient IM with respect to the effects of the pulse-like ground motions on structural response. For near-fault ground motions, including pulse-like and non–pulse-like time histories, we propose a vector-valued IM consisting of a new IM called instantaneous power (IP(T₁)) and the Sa(T₁). The IP(T₁) is defined as the maximum power of the bandpass-filtered velocity time series over a time interval of 0.5T₁. The IP(T₁) is period-dependent because the velocity time series is filtered over a period range (0.2T₁-3T₁). This allows the IP(T₁) to represent the power of the near-fault ground motions relevant to the response of the structure. Using two-dimensional models of the 2- and 9-story steel-frame buildings, we show that the proposed [Sa(T₁), IP(T₁)] vector IM gives more accurate estimates of the maximum inter-story drift and collapse capacity responses from near-fault ground motions than using the vector IM consisting of the Sa(T₁), the presence of the velocity pulse, and the period of the velocity pulse. Moreover, for the structures considered, for a given Sa(T₁), the IP(T₁) is more strongly correlated with structural damage from near-fault ground motions than the combination of the velocity pulse and pulse period.     

Spectral shape proxies and nonlinear structural response

In this paper, spectral-shape-based intensity measures (IMs) are discussed with respect to ordinary, pulse-like and narrow-band records. First, the analyses address the ability of these IMs to capture the peak and cumulative damage potentials of ground motions. Second, a new vector-valued ground motion IM based on the spectral acceleration at the first mode of the structure, Sa(T₁), and a parameter proxy for the spectral shape, namely N_p, is introduced. The vector 〈Sa,N_p〉 is compared to other state-of-the-art IMs in terms of estimation of the seismic response of nonlinear single degree of freedom systems, reinforced concrete and steel moment resisting frames. Results show that 〈Sa,N_p〉 may be especially useful to represent the ground motion potential in the case of records with peculiar spectral shape. Further, it is shown that 〈Sa,N_p〉 has the properties of efficiency, sufficiency and scaling robustness. Finally, a scalar ground motion IM based on Sa(T₁) and N_p is also discussed, and the possibility to compute the seismic hazard analysis for it is illustrated.     

Evaluation of the damping modification factor for structures subjected to near-fault ground motions

The damping modification factor (DMF) has been extensively used in earthquake engineering to describe the variation of structural responses due to varied damping ratios. It is known that DMFs are dependent not only on structural dynamic properties but also on characteristics of ground motions. DMFs regulated in current seismic codes are generally developed based on far-fault ground motions and are inappropriately used in structural design where pulse-like near-fault ground motions are involved. In this paper, statistical investigation of the DMF is performed based on 50 carefully selected pulse-like near-fault ground motions. It is observed that DMFs for pulse-like ground motions exhibit significant dependence on the pulse period T _p in a specific period range. If the period of the structure in response is close to the pulse period, the DMF attains the same level as that derived from far-fault ground motions; as the period of the structure is considerably larger or smaller than the pulse period T _p , the response reduction effect by the increased damping ratio is generally small, except for large earthquakes with long pulse periods, which exhibit significant reduction of response for structures with periods smaller than T _p . Based on the statistical results of DMFs, the empirical formulas for estimating DMFs for displacement, velocity and acceleration spectra are proposed, the effect of structural period, pulse period and damping ratio are considered in the formulas, and the formulas are designed to satisfy the specific reliability requirement in the period range of 0.1 < T/T _p  < 1, which is of engineering interest.     

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.

     

Seismic vulnerability of offshore wind turbines to pulse and non-pulse records

The increasing number of wind turbines in active tectonic regions has attracted scientific interest to evaluate the seismic vulnerability of offshore wind turbines (OWTs). This study aims at assessing the deformation and collapse susceptibility of 2MW and 5MW OWTs subjected to shallow-crustal pulse-like ground motions, which has not been particularly addressed to date. A cloud-based fragility assessment is performed to quantify the seismic response for a given intensity measure and to assess the failure probabilities for pulse-like and non-pulse-like ground motions. The first-mode spectral acceleration S_a(T₁) is found to be an efficient response predictor for OWTs, exhibiting prominent higher-mode behavior, at the serviceability and ultimate conditions. Regardless of earthquake type, it is shown that records with strong vertical components may induce nonlinearity in the supporting tower, leading to potential failure by buckling in three different patterns: (i) at tower base near platform level, (ii) close to tower top, and (iii) between the upper half of the main tower and its top. Type and extent of the damage are related to the coupled excitation of vertical and lateral higher modes, for which tower top acceleration response spectra S_a,i(Top) is an effective identifier. It is also observed that tower's slenderness ratio (l/d), the diameter-to-thickness ratio (d/t), and the rotor-nacelle-assembly mass (m_RNA) are precursors for evaluating the damage mode and vulnerability of OWTs under both pulse-like and non-pulse-like ground motion records.     

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.     

Spectral shape,epsilon and record selection

Selection of earthquake ground motions is considered with the goal of accurately estimating the response of a structure at a specified ground motion intensity, as measured by spectral acceleration at the first‐mode period of the structure, Sa(T₁). Consideration is given to the magnitude, distance and epsilon (ε) values of ground motions. First, it is seen that selecting records based on their ε values is more effective than selecting records based on magnitude and distance. Second, a method is discussed for finding the conditional response spectrum of a ground motion, given a level of Sa(T₁) and its associated mean (disaggregation‐based) causal magnitude, distance and ε value. Records can then be selected to match the mean of this target spectrum, and the same benefits are achieved as when records are selected based on ε. This mean target spectrum differs from a Uniform Hazard Spectrum, and it is argued that this new spectrum is a more appropriate target for record selection. When properly selecting records based on either spectral shape or ε, the reductions in bias and variance of resulting structural response estimates are comparable to the reductions achieved by using a vector‐valued measure of earthquake intensity. Copyright © 2006 John Wiley & Sons, Ltd.     

Structural performance assessment under near‐source pulse‐like ground motions using advanced ground motion intensity measures

This paper demonstrates the effectiveness of utilizing advanced ground motion intensity measures (IMs) to evaluate the seismic performance of a structure subject to near‐source ground motions. Ordinary records are, in addition, utilized to demonstrate the robustness of the advanced IM with respect to record selection and scaling. To perform nonlinear dynamic analyses (NDAs), ground motions need to be selected; as a result, choosing records that are not representative of the site hazard can alter the seismic performance of structures. The median collapse capacity (in terms of IM), for example, can be systematically dictated by including a few aggressive or benign pulse‐like records into the record set used for analyses. In this paper, the elastic‐based IM such as the pseudo‐spectral acceleration (S_a) or a vector of S_a and epsilon has been demonstrated to be deficient to assess the structural responses subject to pulse‐like motions. Using advanced IMs can be, however, more accurate in terms of probabilistic response prediction. Scaling earthquake records using advanced IMs (e.g. inelastic spectral displacement, S_di, and IM _1I&2E; the latter is for the significant higher‐mode contribution structures) subject to ordinary and/or pulse‐like records is efficient, sufficient, and robust relative to record selection and scaling. As a result, detailed record selection is not necessary, and records with virtually any magnitude, distance, epsilon and pulse period can be selected for NDAs. Copyright © 2008 John Wiley & Sons, Ltd.     

Housner谱烈度的一种快速近似算法

Housner谱烈度及修正谱烈度作为基于加速度记录时程直接得到的地震动强度参数,与建筑结构破坏及地震宏观烈度存在较高的相关性,是可靠的地震仪器烈度物理参数指标。然而,相对于地面加速度峰值、地面速度峰值等地震动峰值参数,三分量加速度记录对应的谱烈度计算过程较为复杂,耗时相对较长,影响了利用谱烈度确定地震仪器烈度的时效性。基于对强震动加速度记录的统计分析,本文提出了谱烈度的快速近似算法,仅计算4个方向上的谱烈度值,采用其中3点作圆即可获得水平面内谱烈度迹线的近似最大值,使计算速度提高了45倍,且保持了谱烈度作为地震仪器烈度物理指标的精度。利用在汶川M_S 8.0地震等386次M_S > 3.0地震中获取的2701组强震动加速度记录,经可靠性检验,结果表明所提出的Housner谱烈度快速近似算法的计算误差在±4.5%以内,可以同时满足地震仪器烈度速报的可靠性和时效性需求。     

Near‐fault effects on residual displacements of RC structures

Residual displacements of single‐degree‐of‐freedom systems due to ground motions with velocity pulses or fling step displacements are presented as a function of period T and of its ratio to the pulse period T_p. Four hysteretic behaviors are considered: bilinear elastoplastic, stiffness‐degrading with cycling, stiffness‐cum‐strength degrading, with or without pinching. When expressed in terms of T/T_p, peak inelastic and residual displacements due to motions with a pulse or fling appear similar to those due to far‐fault motions, if the response to far‐field records are expressed in terms of the ratio of T to the record's characteristic period. However, as the latter is usually much shorter than the pulse period of motions with fling, the range of periods of interest for common structures becomes a short‐period range under fling motions and exhibits very large amplification of residual and peak inelastic displacements. Similar, but less acute, are the effects of motions with a velocity pulse. Wavelets of different complexity are studied as approximations to near‐fault records. Simple two‐parameter wavelets for fling motions overestimate peak inelastic displacements; those for pulse‐type motions overestimate residual displacements. A more complex four‐parameter wavelet for motions with a velocity pulse predicts overall well residual and peak displacements due to either pulse‐ or fling‐type motions; a hard‐to‐identify parameter of the wavelet impacts little computed residual displacements; another significantly affects them and should be carefully estimated from the record. Even this most successful of wavelets overpredicts residual displacements for the periods of engineering interest. Copyright © 2016 John Wiley & Sons, Ltd.     

Average spectral acceleration as an intensity measure for collapse risk assessment

This paper investigates the performance of spectral acceleration averaged over a period range (Sa_avg) as an intensity measure (IM) for estimating the collapse risk of structures subjected to earthquake loading. The performance of Sa_avg is evaluated using the following criteria: efficiency, sufficiency, the availability or ease of developing probabilistic seismic hazard information in terms of the IM and the variability of collapse risk estimates produced by the IM. Comparisons are also made between Sa_avg and the more traditional IM: spectral acceleration at the first‐mode period of the structure (Sa(T₁)). Though most previous studies have evaluated IMs using a relatively limited set of structures, this paper considers nearly 700 moment‐resisting frame and shear wall structures of various heights to compare the efficiency and sufficiency of the IMs. The collapse risk estimates produced by Sa_avg and Sa(T₁) are also compared, and the variability of the risk estimates is evaluated when different ground motion sets are used to assess the structural response. The results of this paper suggest that Sa_avg, when computed using an appropriate period range, is generally more efficient, more likely to be sufficient and provides more stable collapse risk estimates than Sa(T₁). Copyright © 2015 John Wiley & Sons, Ltd.     

Seismic response estimation method for earthquake‐damaged RC buildings

This study proposes a procedure for identifying spectral response curves for earthquake‐damaged areas in developing countries without seismic records. An earthquake‐damaged reinforced concrete building located in Padang, Indonesia was selected to illustrate the identification of the maximum seismic response during the 2009 West Sumatra earthquake. This paper summarizes the damage incurred by the building; the majority of the damage was observed in the third story in the span direction. The damage was quantitatively evaluated using the damage index R according to the Japanese guidelines for post‐earthquake damage evaluation. The damage index was also applied to the proposed spectral response identification method. The seismic performance of the building was evaluated by a nonlinear static analysis. The analytical results reproduced a drift concentration in the third story. The R‐index decreased with an increase in the story drift, which provided an estimation of the maximum response of the building during the earthquake. The estimation was verified via an earthquake response analysis of the building using ground acceleration data, which were simulated based on acceleration records of engineering bedrock that considered site amplification. The maximum response estimated by the R‐index was consistent with the maximum response obtained from the earthquake response analysis. Therefore, the proposed method enables the construction of spectral response curves by integrating the identification results for the maximum responses in a number of earthquake‐damaged buildings despite a lack of seismic records. Copyright © 2016 The Authors. Earthquake Engineering & Structural Dynamics published by John Wiley & Sons Ltd.     

Ground-motion models for average spectral acceleration in a period range: direct and indirect methods

Average spectral acceleration, AvgSA, is defined as the geometric mean of spectral acceleration values over a range of periods and it is a ground motion intensity measure used for structural response prediction. One of its advantages stands on the assumption that its distribution is computable from the available GMPEs for spectral acceleration, GMPE-SA, (called here indirect method) without the need for deriving new specific GMPEs for AvgSA, GMPE-AvgSA, (called here direct method). To what extent this assumption is valid, however, has never been verified. As such, we derived an empirical GMPE-AvgSA based on RESORCE ground motion dataset and we compared its predicted values with those from a GMPE-SA via the indirect approach. As expected, the results show that the indirect approach yields median AvgSA estimates that are identical to those of the direct approach. However, the estimates of AvgSA variance of the two methods are identical only if both the GMPE-SA and their empirical correlation coefficients among different SA ordinates are derived from the same record dataset.     

2021年青海玛多7.4级地震强震动记录及特征分析

2021年5月22日,青海省果洛州玛多县发生7.4级地震,中国强震动观测网络在主震中捕获16组强震动数据。对48条三分向加速度记录进行基线校正、滤波等常规处理,计算相应的地震动参数,发现位于断层破裂前向位置的63DAW台NS向记录的地震动速度波形具有长周期分量丰富的特征。分析6个典型台站的单自由度加速度反应谱,并与我国建筑抗震设计谱比较,分析此次地震的频谱特性。将实际观测到的PGA、PGV和S_a(T=0.1 s、T=1.0 s、T=2.0 s、T=5.0 s)与国内广泛使用的几种地震动预测模型对比,研究此次地震的影响场。通过分析S_a-S_d曲线,探讨此次地震靠近断层区域地面运动大位移与桥梁落梁震害间可能存在的关系。     

Pounding force assessment in performance‐based design of bridges

Bridges with deck supported on either sliding or elastomeric bearings are very common in mid‐seismicity regions. Their main seismic vulnerabilities are related to the pounding of the deck against abutments or between the different deck elements. A simplified model of the longitudinal behavior of those bridges will allow to characterize the reaction forces developed during pounding using the Pacific Earthquake Engineering Research Center framework formula. In order to ensure the general applicability of the results obtained, a large number of system parameter combinations will be considered. The heart of the formula is the identification of suitable intermediate variables. First, the pseudo acceleration spectral value for the fundamental period of the system (Sa(T_s)) will be used as an intensity measure (IM). This IM will result in a very large non‐explained variability of the engineering demand parameter. A portion of this variability will be proved to be related to the relative content of high‐frequency energy in the input motion. Two vector‐valued IMs including a second parameter taking this energy content into account will then be considered. For both of them, a suitable form for the conditional intensity dependence of the response will be obtained. The question of which one to choose will also be analyzed. Finally, additional issues related to the IM will be studied: its applicability to pulse‐type records, the validity of scaling records and the sufficiency of the IM. Copyright © 2009 John Wiley & Sons, Ltd.     

Near-fault acceleration pulses and non-acceleration pulses: Effects on the inelastic displacement ratio

Near-fault ground motions can impose particularly high seismic demands on the structures due to the pulses that are typically observed in the velocity time-histories. The velocity pulses can be further categorized into either a distinct acceleration pulse (acc-pulse) or a succession of high-frequency, one-sided acceleration spikes (non-acc-pulse). The different characteristics of velocity pulses imply different frequency content of the ground motions, potentially causing different seismic effects on the structures. This study aims to investigate the characteristics of the two types of velocity pulses and their impacts on the inelastic displacement ratio (C_R) of single-degree-of-freedom systems. First, a new method that enables an automated classification of velocity pulses is used to compile a ground motion dataset which consists of 74 acc-pulses and 45 non-acc-pulses. Several intensity measures characterizing different seismological features are then compared using the two groups of records. Finally, the influences of acc-pulses and non-acc-pulses on the C_R spectra are studied; the effects of pulse period and hysteretic behavior are also considered. Results indicate that the characteristics of the two types of velocity pulses differ significantly, resulting in clearly distinct C_R spectral properties between acc-pulses and non-acc-pulses. Interestingly, mixing acc-pulses and non-acc-pulses can lead to local “bumps” that were found in the C_R spectral shape by previous studies. The findings of this study highlight the importance of distinguishing velocity pulses of different types when selecting near-fault ground motions for assessing the nonlinear dynamic response of structures.     

Comparing response of SDF systems to near‐fault and far‐fault earthquake motions in the context of spectral regions

In spite of important differences in structural response to near‐fault and far‐fault ground motions, this paper aims at extending well‐known concepts and results, based on elastic and inelastic response spectra for far‐fault motions, to near‐fault motions. Compared are certain aspects of the response of elastic and inelastic SDF systems to the two types of motions in the context of the acceleration‐, velocity‐, and displacement‐sensitive regions of the response spectrum, leading to the following conclusions. (1) The velocity‐sensitive region for near‐fault motions is much narrower, and the acceleration‐sensitive and displacement‐sensitive regions are much wider, compared to far‐fault motions; the narrower velocity‐sensitive region is shifted to longer periods. (2) Although, for the same ductility factor, near‐fault ground motions impose a larger strength demand than far‐fault motions—both demands expressed as a fraction of their respective elastic demands—the strength reduction factors R_y for the two types of motions are similar over corresponding spectral regions. (3) Similarly, the ratio u_m/u₀ of deformations of inelastic and elastic systems are similar for the two types of motions over corresponding spectral regions. (4) Design equations for R_y (and for u_m/u₀) should explicitly recognize spectral regions so that the same equations apply to various classes of ground motions as long as the appropriate values of T_a, T_b and T_c are used. (5) The Veletsos–Newmark design equations with T_a=0.04 s, T_b=0.35 s, and T_c=0.79 s are equally valid for the fault‐normal component of near‐fault ground motions. Copyright © 2001 John Wiley & Sons, Ltd.     

云南漾濞6.4级地震强地面运动三要素基本特征分析

2021年5月21日云南漾濞发生6.4级地震,中国数字化强震动观测网络在主震中捕获28组加速度数据。对84条三分向加速度记录进行常规处理,计算出对应的地震动参数,并绘制震中附近水平向PGA和PGV等值线图,结果显示其走向和形态与震后烈度调查结果相近。将实际观测的相对持时进行曲线回归,求得相关系数,并据此计算出T_d(5%～95%)持时对应的震源持续时间在10.22～13.06 s间,与震源破裂过程反演结果基本一致。再通过分析近场6个台站记录的加速度反应谱,并与抗震设计谱比较,总结此次地震的频谱特性,为该区域抗震设防提供参考依据。     

Evaluation of bias on the probability of collapse from amplitude scaling using spectral-shape-matched records

Although for many years it was thought that amplitude scaling of acceleration time series to reach a target intensity did not introduce any bias in the results of nonlinear response history analyses, recent studies have showed that scaling can lead to an overestimation of deformation demands with increasing scale factors. Some studies have suggested that the bias can be explained by differences in spectral shape between the response spectra of unscaled and scaled records. On the basis of these studies, some record selection procedures assume that if records are selected using spectral-shape-matching procedures, amplitude scaling does not induce any bias on the structural response. This study evaluates if bias is introduced on lateral displacement demands and seismic collapse risk estimates even when spectral shape is carefully taken into consideration when selecting ground motions. Several single-degree-of-freedom and multiple-degree-of-freedom systems are analyzed when subjected to unscaled and scaled ground motions selected to approximately match the mean and the variance of the conditional spectrum at the target level of intensity. Results show that an explicit consideration of spectral shape is not enough to avoid a systematic overestimation of lateral displacement demands and collapse probabilities as the scale factor increases. Moreover, the bias is observed in practically all cases for systems with strength degradation and it increases with decreasing period and decreasing lateral strength relative to the strength required to remain elastic. Key reasons behind the bias are presented by evaluating input energy, causal parameters, and damaging pulse distributions in unscaled and scaled ground motion sets.