Analytical fragility assessment using unscaled ground motion records

It is desirable that nonlinear dynamic analyses for structural fragility assessment are performed using unscaled ground motions. The widespread use of a simple dynamic analysis procedure known as Cloud Analysis, which uses unscaled records and linear regression, has been impeded by its alleged inaccuracies. This paper investigates fragility assessment based on Cloud Analysis by adopting, as the performance variable, a scalar demand to capacity ratio that is equal to unity at the onset of limit state. It is shown that the Cloud Analysis, performed based on a careful choice of records, leads to reasonable and efficient fragility estimates. There are 2 main rules to keep in mind for record selection: to make sure that a good portion of the records leads to a demand to capacity ratio greater than unity and that the dispersion in records' seismic intensity is considerable. An inevitable consequence of implementing these rules is that one often needs to deal with the so‐called collapse cases. To formally consider the collapse cases, a 5‐parameter fragility model is proposed that mixes the simple regression in the logarithmic scale with logistic regression. The joint distribution of fragility parameters can be obtained by adopting a Markov Chain Monte Carlo simulation scheme leading directly to the fragility and its confidence intervals. The resulting fragility curves compare reasonably with those obtained from the Incremental Dynamic Analysis and Multiple Stripe Analysis with (variable) conditional spectrum–compatible suites of records at different intensity levels for 3 older reinforced concrete frames with shear‐, shear‐flexure‐, and flexure‐dominant behavior.     

Comparing the adequacy of alternative ground motion intensity measures for the estimation of structural responses

Practical methods for the probability‐based seismic assessment of structures make use of estimates of demands produced by earthquakes of different intensities. The uncertainties associated with these estimates are highly dependent on the variable adopted as the intensity measure (IM, e.g., PGA, spectral acceleration, etc.). This generates the need to compare the efficiency of an originally adopted IMwith that of a new candidate. This implies comparing the dispersion of the demand measure (DM, e.g., maximum interstorey drift ratio, ductility demand, etc.) conditional to each of the two IMs. In order to obtain the demand estimates in a conventional way, a full set of dynamic response analyses should be performed for each IM under scrutiny, i.e., multiple records scaled at several fixed values of each IM. The procedure developed here serves to accelerate this comparison avoiding the effort required to evaluate the dynamic responses of the structure for all the ground motion time histories considered every time that a new IM is adopted. For this purpose, use is made of available results of analyses performed for a different (i.e., the original) IM. Two methods are proposed: the direct method involves performing a regression of the results obtained from the original analyses, taking the candidate IM as the independent variable. The indirect method involves rebuilding the probability density function of the DM given a defined value of the candidate IM by means of the total probability theorem, using the results of the original analyses and certain data relating the two IMs. The proposed methods have been tested by application to several SDOF systems with different periods and different cyclic‐response backbone curves. The conditions affecting their approximation are explored, and some criteria to improve them are identified. The procedure can also be used to determine the optimum value of a parameter to be used in a parameter‐based IM. Copyright © 2004 John Wiley & Sons, Ltd.     

Cloud to IDA: Efficient fragility assessment with limited scaling

Incremental dynamic analysis (IDA) leads to curves expressed in terms of structural response versus intensity, commonly known as the IDA curves. It is known that implementation of IDA usually involves significant computational effort and most often significant scaling of the original records to various intensity levels. Employing as the performance variable the critical demand to capacity ratio (DCR) throughout the structure, which is equal to unity at the onset of the limit state, facilitates the identification of the intensity values at the onset of a desired limit state and hence the implementation of the IDA procedure. Employing the structural response to un‐scaled records and the corresponding regression‐based response predictions (a.k.a., the “Cloud Analysis”) helps in identifying the range of intensity values corresponding to demand to capacity ratio values in the vicinity of unity. The Cloud to IDA procedure for structural fragility assessment is proposed on the premise of exploiting the Cloud Analysis results to obtain the IDA curves both with minimum number of analyses and minimum amount of scaling. The transverse frame of a shear‐critical 7‐story older RC building in Van Nuys, CA, which is modeled in Opensees with fiber‐section considering the flexural‐shear‐axial interactions and the bar slip, is employed as a case study. It is demonstrated, by comparing the results to IDA and other state of the art non‐linear dynamic procedures based on no scaling or spectral‐shape‐compatible scaling, that the Cloud to IDA procedure leads to reliable results in terms of structural fragility and risk for the prescribed limit state.     

A web‐based methodology for the prediction of approximate IDA curves

A web‐based methodology for the prediction of approximate IDA curves, which consists of two independent processes, is proposed. The result of the first process is a response database of the SDOF model, whereas the second process involves the prediction of approximate IDA curves from the response database by using n‐dimensional linear interpolation. Such an approach enables user‐friendly prediction of the seismic response parameters with high accuracy. In order to demonstrate the capabilities of the proposed methodology, a web application for the prediction of the approximate 16th, 50th and 84th fractile responses of an RC structure was developed. For the presented case study, the response database was computed for a set of 30 ground motion records and the discrete values of six structural parameters. Very good agreement between the computed and the approximated IDA curves of the four‐storey RC building was observed. Copyright © 2012 John Wiley & Sons, Ltd.     

A ground motion selection procedure for enforcing hazard consistency and estimating seismic demand hazard curves

This paper develops a procedure to select unscaled ground motions for estimating seismic demand hazard curves (SDHCs) in performance‐based earthquake engineering. Currently, SDHCs are estimated from a probabilistic seismic demand analysis, where several ensembles of ground motions are selected and scaled to a user‐specified scalar conditioning intensity measure (IM). In contrast, the procedure developed herein provides a way to select a single ensemble of unscaled ground motions for estimating the SDHC. In the context of unscaled motions, the proposed procedure requires three inputs: (i) database of unscaled ground motions, (ii) I M , the vector of IMs for selecting ground motions, and (iii) sample size, n; in the context of scaled motions, two additional inputs are needed: (i) a maximum acceptable scale factor, SF_max, and (ii) a target fraction of scaled ground motions, γ. Using a recently developed approach for evaluating ground motion selection and modification procedures, the proposed procedure is evaluated for a variety of inputs and is demonstrated to provide accurate estimates of the SDHC when the vector of IMs chosen to select ground motions is sufficient for the response quantity of interest. Copyright © 2015 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.     

The effective number of cycles of earthquake ground motion

The seismic response of any system that accumulates damage under cyclic loading is dependent not only on the maximum amplitude of the motion but also its duration. This is explicitly recognized in methods for estimating the liquefaction potential of soil deposits. Many researchers have proposed that the effective number of cycles of the ground motion is a more robust indicator of the destructive capacity of the shaking than the duration. However, as is the case with strong‐motion duration, there is no universally accepted approach to determining the effective number of cycles of motion, and the different methods that have been proposed can give widely varying results for a particular accelerogram. Definitions of the effective number of cycles of motion are reviewed, classified and compared. Measurements are found to differ particularly for accelerograms with broad‐banded frequency content, which contain a significant number of non‐zero crossing peaks. The key seismological parameters influencing the number of cycles of motion and associated equations for predicting this quantity for future earthquakes are identified. Correlations between cycle counts and different duration measures are explored and found to be rather poor in the absence of additional parameters. Copyright © 2004 John Wiley & Sons, Ltd.     

地震强度指标选择对长周期钢结构余震IDA分析的影响

结构在主震作用下发生损伤,由于刚度和强度退化导致结构周期的延长,再加上高阶周期的影响,在随后的余震对主震损伤结构抗震性能的影响分析中,基于弹性周期的传统谱加速度可能不适用,有必要对地震动强度指标的选择进行进一步研究。本文以9层Benchmark钢框架为研究对象,选取10种地震动强度指标,分别从离散性和相关性两个维度对比分析了其在余震IDA分析中的优劣性。研究结果表明:与最大层间位移角相比,最大残余层间位移角更能够准确地量化主震对结构造成的损伤;对于钢框架结构,基于一段周期的几何平均谱加速度的IM_Boj&Lev指标是最优的地震动强度指标,相对于S_a（T₁）离散性降低了30%左右,相关性提高了6%左右。     

Direct estimation of the seismic demand and capacity of oscillators with multi‐linear static pushovers through IDA

SPO2IDA is introduced, a software tool that is capable of recreating the seismic behaviour of oscillators with complex quadrilinear backbones. It provides a direct connection between the static pushover (SPO) curve and the results of incremental dynamic analysis (IDA), a computer‐intensive procedure that offers thorough demand and capacity prediction capability by using a series of nonlinear dynamic analyses under a suitably scaled suite of ground motion records. To achieve this, the seismic behaviour of numerous single‐degree‐of‐freedom (SDOF) systems is investigated through IDA. The oscillators have a wide range of periods and feature pinching hysteresis with backbones ranging from simple bilinear to complex quadrilinear with an elastic, a hardening and a negative‐stiffness segment plus a final residual plateau that terminates with a drop to zero strength. An efficient method is introduced to treat the backbone shape by summarizing the analysis results into the 16, 50 and 84% fractile IDA curves, reducing them to a few shape parameters and finding simpler backbones that reproduce the IDA curves of complex ones. Thus, vast economies are realized while important intuition is gained on the role of the backbone shape to the seismic performance. The final product is SPO2IDA, an accurate, spreadsheet‐level tool for performance‐based earthquake engineering that can rapidly estimate demands and limit‐state capacities, strength reduction R‐factors and inelastic displacement ratios for any SDOF system with such a quadrilinear SPO curve. Copyright © 2006 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.     

A critical examination of seismic response uncertainty analysis in earthquake engineering

The last decade of performance‐based earthquake engineering (PBEE) research has seen a rapidly increasing emphasis placed on the explicit quantification of uncertainties. This paper examines uncertainty consideration in input ground‐motion and numerical seismic response analyses as part of PBEE, with particular attention given to the physical consistency and completeness of uncertainty consideration. It is argued that the use of the commonly adopted incremental dynamic analysis leads to a biased representation of the seismic intensity and that when considering the number of ground motions to be used in seismic response analyses, attention should be given to both reducing parameter estimation uncertainty and also limiting ground‐motion selection bias. Research into uncertainties in system‐specific numerical seismic response analysis models to date has been largely restricted to the consideration of ‘low‐level’ constitutive model parameter uncertainties. However, ‘high‐level’ constitutive model and model methodology uncertainties are likely significant and therefore represent a key research area in the coming years. It is also argued that the common omission of high‐level seismic response analysis modelling uncertainties leads to a fallacy that ground‐motion uncertainty is more significant than numerical modelling uncertainty. The author's opinion of the role of uncertainty analysis in PBEE is also presented. Copyright © 2013 John Wiley & Sons, Ltd.     

Approximate incremental dynamic analysis using the modal pushover analysis procedure

Incremental dynamic analysis (IDA)—a procedure developed for accurate estimation of seismic demand and capacity of structures—requires non‐linear response history analysis of the structure for an ensemble of ground motions, each scaled to many intensity levels, selected to cover the entire range of structural response—all the way from elastic behaviour to global dynamic instability. Recognizing that IDA of practical structures is computationally extremely demanding, an approximate procedure based on the modal pushover analysis procedure is developed. Presented are the IDA curves and limit state capacities for the SAC‐Los Angeles 3‐, 9‐, and 20‐storey buildings computed by the exact and approximate procedures for an ensemble of 20 ground motions. These results demonstrate that the MPA‐based approximate procedure reduces the computational effort by a factor of 30 (for the 9‐storey building), at the same time providing results to a useful degree of accuracy over the entire range of responses—all the way from elastic behaviour to global dynamic instability—provided a proper hysteretic model is selected for modal SDF systems. The accuracy of the approximate procedure does not deteriorate for 9‐ and 20‐storey buildings, although their dynamics is more complex, involving several ‘modes’ of vibration. For all three buildings, the accuracy of the MPA‐based approximate procedure is also satisfactory for estimating the structural capacities for the limit states of immediate occupancy, collapse prevention, and global dynamic instability. Copyright © 2006 John Wiley & Sons, Ltd.     

Seismic response analysis using characteristic ground motion records for risk‐based decision‐making (3R method)

The concept of intensity‐based assessment for risk‐based decision‐making is introduced. It is realized by means of the so‐called 3R method (response analysis, record selection and risk‐based decision‐making), which can be used to check the adequacy of design of a new building or of the strengthening of an existing building by performing conventional pushover analysis and dynamic analysis for only a few ground motions, which are termed characteristic ground motions. Because the objective of the method is not a precise assessment of the seismic risk, a simple decision model for risk acceptability can be introduced. The engineer can decide that the reliability of a no‐collapse requirement is sufficient when collapse is observed in the case of less than half of, for example, seven characteristic ground motions. From the theoretical point of view, it is shown that the accuracy of the method is acceptable if the non‐linear response history analyses are performed at a low percentile of limit‐state intensity, which is also proven by means of several examples of multi‐storey reinforced concrete frame buildings. The 3R method represents a compromise between the exclusive use of either pushover analysis or dynamic analysis and can be easily introduced into building codes provided that its applicability is further investigated (e.g. asymmetric structures and other performance objectives) and that the procedure for the selection of characteristic ground motions is automated and readily available to engineers (www.smartengineering.si).     

Seismic intensity measures for above‐ground liquid storage tanks

A series of scalar and vector intensity measures is examined to determine their suitability within the seismic risk assessment of liquid storage tanks. Using a surrogate modelling approach on a squat tank that is examined under both anchored and unanchored support conditions, incremental dynamic analysis is adopted to generate the distributions of response parameters conditioned on each of the candidate intensity measures. Efficiency and sufficiency metrics are used in order to perform the intensity measure evaluation for individual failure modes, while a comparison in terms of mean annual frequency of exceedance is performed with respect to a damage state that is mutually governed by the impulsive and convective modes of the tank. The results reveal combinations of spectral acceleration ordinates as adequate predictors, among which the average spectral acceleration is singled out as the optimal solution. The sole exception is found for the sloshing‐controlled modes of failure, where mainly the convective period spectral acceleration is deemed adequate to represent the associated response due to their underlying linear relationship. A computationally efficient method in terms of site hazard analysis is finally proposed to serve in place of the vector‐valued intensity measures, providing a good match for the unanchored tank considered and a more conservative one for the corresponding anchored system.     

Evaluation of the exact conditional spectrum and generalized conditional intensity measure methods for ground motion selection

Two existing, contemporary ground motion selection and modification procedures – (i) exact conditional spectrum (CS‐exact) and (ii) generalized conditional intensity measure (GCIM) – are evaluated in their ability to accurately estimate seismic demand hazard curves (SDHCs) of a given structure at a specified site. The amount of effort involved in implementing these procedures to compute a single SDHC is studied, and a case study is chosen where rigorous benchmark SDHCs can be determined for evaluation purposes. By comparing estimates from ground motion selection and modification procedures with the benchmark, we conclude that estimates from CS‐exact are unbiased in many of the cases considered. The estimates from GCIM are even more accurate, as they are unbiased for most – but not all – of the cases where estimates from CS‐exact are biased. We find that it is possible to obtain biased SDHCs from GCIM, even after employing a very diverse collection of intensity measures to select ground motions and implementing its bias‐checking feature, because it is usually difficult to identify intensity measures that are truly ‘sufficient’ for the response of a complex, multi‐degree‐of‐freedom system. Copyright © 2015 John Wiley & Sons, Ltd.     

Earthquake ground motion modeling of induced seismicity in the Groningen gas field

A physics‐based numerical approach is used to characterize earthquake ground motion due to induced seismicity in the Groningen gas field and to improve empirical ground motion models for seismic hazard and risk assessment. To this end, a large‐scale (20 km × 20 km) heterogeneous 3D seismic wave propagation model for the Groningen area is constructed, based on the significant bulk of available geological, geophysical, geotechnical, and seismological data. Results of physics‐based numerical simulations are validated against the ground motion recordings of the January 8, 2018, M_L 3.4 Zeerijp earthquake. Taking advantage of suitable models of slip time functions at the seismic source and of the detailed geophysical model, the numerical simulations are found to reproduce accurately the observed features of ground motions at epicentral distances less than 10 km, in a broad frequency range, up to about 8 Hz. A sensitivity analysis is also addressed to discuss the impact of 3D underground geological features, the stochastic variability of seismic velocities and the frequency dependence of the quality factor. Amongst others, results point out some key features related to 3D seismic wave propagation, such as the magnitude and distance dependence of site amplification functions, that may be relevant to the improvement of the empirical models for earthquake ground motion prediction.     

Conditional spectrum‐based ground motion selection. Part II: Intensity‐based assessments and evaluation of alternative target spectra

In a companion paper, an overview and problem definition was presented for ground motion selection on the basis of the conditional spectrum (CS), to perform risk‐based assessments (which estimate the annual rate of exceeding a specified structural response amplitude) for a 20‐story reinforced concrete frame structure. Here, the methodology is repeated for intensity‐based assessments (which estimate structural response for ground motions with a specified intensity level) to determine the effect of conditioning period. Additionally, intensity‐based and risk‐based assessments are evaluated for two other possible target spectra, specifically the uniform hazard spectrum (UHS) and the conditional mean spectrum (CMS, without variability).It is demonstrated for the structure considered that the choice of conditioning period in the CS can substantially impact structural response estimates in an intensity‐based assessment. When used for intensity‐based assessments, the UHS typically results in equal or higher median estimates of structural response than the CS; the CMS results in similar median estimates of structural response compared with the CS but exhibits lower dispersion because of the omission of variability. The choice of target spectrum is then evaluated for risk‐based assessments, showing that the UHS results in overestimation of structural response hazard, whereas the CMS results in underestimation. Additional analyses are completed for other structures to confirm the generality of the conclusions here. These findings have potentially important implications both for the intensity‐based seismic assessments using the CS in future building codes and the risk‐based seismic assessments typically used in performance‐based earthquake engineering applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Preliminary research on the method for determination of ground motion parameters of diffuse earthquakes

This paper proposes that on the basis of near-source saturation characteristic of ground motion, the peak ground motion acceleration and the acceleration response spectra of diffuse earthquakes affecting the nuclear facility site can be determined by the statistic method. In view of the problem of insufficient statistic samples due to fewer qualified strong earthquake records on the bedrock site, it also proposes that according to the features of ground motion response of the soil, the acceleration records ...     

保定市设计地震动参数的研究

本文根据保定市及周围地区的地震地质环境，在地震危险性分析的基础上，采用等效线性一维波动方程进行土层的地震反应分析。给出50年超越概率63％、10％、2％基岩和地面的水平向峰值加速度、反应谱（场址基本烈度Ⅶ度）和地震影响系数最大值。该结果为抗震设计提供了可靠依据，具有应用价值。