Intensity measure conversion of fragility curves

In seismic risk assessment of structures, fragility functions are the probabilistic characterization of vulnerability at the component and/or structural level, expressing the probability of failure as a function of a ground motion intensity measure (IM). Fragility curves, in general, are structure- and site-specific, thus a comparison of fragility curves, then of vulnerability, is not straightforward across multiple structures. Also, it could be the case that hazard at a site of interest is not available for the IM originally considered in the fragility assessment. These situations require to convert fragility curves from an original IM to a target one. The present study addresses a hazard-consistent probabilistic framework for converting spectral acceleration-based IMs from an original IM to a target IM at a given site. In particular, three conversion cases, under different assumptions on the explanatory power of the involved IMs with respect to structural failure, are discussed: (a) a vector-valued IM consisting of the original and target IMs, magnitude, and source-to-site distance; (b) a vector-valued IM consisting of the original and target IMs; and (c) the original (scalar) IM only, assuming that structural response, given the IM, is statistically independent of the other ground motion variables. In this framework, the original fragility functions are characterized using the state-of-the-art methods in performance-based earthquake engineering, then the fragility curves as a function of the target IM are evaluated through applications of the probability calculus rules, ensuring consistency with the seismic hazard at the site of interest. The conversion strategy is illustrated through the applications to three-, six-, and nine-story Italian code-conforming reinforced concrete buildings designed for a high-hazard site in Italy. The study shows that, in most of the cases, the converted fragility curves have agreement with the reference curves directly developed in terms of the target IM. Cases in which least agreement was found are likely due to the models used to obtain the terms required by the conversion equations.     

Vector-intensity measure based seismic vulnerability analysis of bridge structures

This paper presents a method for seismic vulnerability analysis of bridge structures based on vector-valued intensity measure(v IM), which predicts the limit-state capacities efficiently with multi-intensity measures of seismic event. Accounting for the uncertainties of the bridge model, ten single-bent overpass bridge structures are taken as samples statistically using Latin hypercube sampling approach. 200 earthquake records are chosen randomly for the uncertainties of ground motions according to the site condition of the bridges. The uncertainties of structural capacity and seismic demand are evaluated with the ratios of demand to capacity in different damage state. By comparing the relative importance of different intensity measures, Sa(T1) and Sa(T2) are chosen as v IM. Then, the vector-valued fragility functions of different bridge components are developed. Finally, the system-level vulnerability of the bridge based on v IM is studied with DunnettSobel class correlation matrix which can consider the correlation effects of different bridge components. The study indicates that an increment IMs from a scalar IM to v IM results in a significant reduction in the dispersion of fragility functions and in the uncertainties in evaluating earthquake risk. The feasibility and validity of the proposed vulnerability analysis method is validated and the bridge is more vulnerable than any components.     

Impact of seismic excitation direction on the fragility analysis of horizontally curved concrete bridges

The focus of this study is the impact of the seismic excitation direction on the fragility of horizontally curved bridges. Nonlinear time history analyses are performed on a typical, curved concrete bridge in China using a set of real ground motions with different incident angles. To build reliable probabilistic seismic demand models, ten commonly used intensity measures (IMs) are assessed in terms of various metrics to determine the optimal IMs, which account for the influence of the seismic excitation directions. Subsequently, fragility surfaces with respect to both the optimal IM and incident angles are generated to qualify the fragility sensitivity for various components and the bridge system to the seismic excitation directions. Moreover, the rationality and applicability of the methods recommended by the Caltrans, Eurocode 8 and Chinese codes for determining the seismic excitation direction of curved bridges are evaluated. The results indicate that the excitation direction imposes a minor impact on the optimal IM rankings. Compared to structure-independent IMs, structure-dependent IMs are more appropriate for predicting the demands of horizontally curved concrete bridges. However, the seismic excitation direction significantly affects the component fragilities, and the level of the effect intensifies with increasing limit states. If the incident angle occurrence probability is not provided, the Chinese code method for the seismic excitation direction is more suitable for the horizontally curved concrete bridge fragility assessment, which has the advantages of computational efficiency when compared to the Caltrans code and relatively conservative results when compared to Eurocode 8.     

Selection of optimal intensity measures in probabilistic seismic demand models of highway bridge portfolios

Probabilistic seismic demand models are a common and often essential step in generating analytical fragility curves for highway bridges. With these probabilistic models being traditionally conditioned on a single seismic intensity measure (IM), the degree of uncertainty in the models is dependent on the IM used. Selection of an optimal IM for conditioning these demand models is not a trivial matter and has been the focus of numerous studies. Unlike previous studies that consider a single structure for IM selection, this study evaluates optimal IMs for use when generating probabilistic seismic demand models for bridge portfolios such as would be found in HAZUS‐MH. Selection criteria such as efficiency, practicality, sufficiency, and hazard computability are considered in the selection process. A case study is performed considering the multi‐span simply supported steel girder bridge class. Probabilistic seismic demand models are generated considering variability in the geometric configurations and material properties, using two suites of ground motions—one synthetic and one recorded motion suite. Results show that of the 10 IMs considered, peak ground acceleration (PGA) and spectral acceleration at the fundamental period are the most optimal for the synthetic motions, and that cumulative absolute velocity is also a close contender when using recorded motions. However, when hazard computability is considered, PGA is selected as the IM of choice. Previous studies have shown that spectrally based quantities perform better than PGA for a given structure, but the findings of this study indicate that when a portfolio of bridges is considered, PGA should be used. Copyright © 2007 John Wiley & Sons, Ltd.     

Performance-based seismic assessment of a large diameter extended pile shaft in a cohesionless soil

The seismic behavior of a large diameter extended pile shaft founded on a dense sandy site is investigated in this paper. First, a deterministic analysis is conducted including both nonlinear dynamic analysis(NDA) and pushover analysis to gain insights into the behavior of the pile and make sure an appropriate modeling technique is utilized. Then a probabilistic analysis is performed using the results of NDA for various demands. To this end a set of 40 pulse-like ground motions are picked and subsequently 40 nonlinear dynamic and pushover analyses are performed. The data obtained from NDA are used to generate probabilistic seismic demand model(PSDM) plots and consequently the median line and dispersion for each plot are computed. The NDA and pushover data are also plotted against each other to find out to what extent they are correlated. These operations are done for various engineering demand parameters(EDPs). A sensitivity analysis is done to pick the most appropriate intensity measure(IM) which would cause a minimum dispersion in PSDM plots out of 7 different IMs. Peak ground acceleration(PGA) is found to be the most appropriate IM. Pushover coefficient equations as a function of PGA are proposed which can be applied to the pushover analysis data to yield a better outcome with respect to the NDA. At the end, the pacific earthquake engineering research(PEER) center methodology is utilized to generate the fragility curves using the properties obtained from PSDM plots and considering various states of damage ranging from minor to severe. The extended pile shaft shows more vulnerability with a higher probability with respect to minor damage compared to severe damage.     

地面运动强度参数对RC桥梁地震需求的影响

在基于性能的地震工程学(PBEE)中,建立概率地震需求模型(PSDM)时需要对桥梁结构的工程需求参数(EDP)进行概率估计。其中,强地面运动参数(IM)的选择对EDP的概率估计影响很大,因此需要正确选择IM。分别采用目前最广泛使用的结构第一模态周期弹性谱加速度(5%阻尼比)Sa(T1,5%)和峰值地面加速度PGA作为IM,选择实际地震波并进行合理的调值,对一座钢筋混凝土桥墩进行IDA分析,其计算结果表明:对于不同性质EDP的概率估计值,以PGA作为IM计算所得的结果明显偏于非保守,且离散度一般也更大。说明可以针对不同性质的EDP,根据地面运动强度的大小,选择不同的IM,通过合理的调值对EDP进行概率估计,可以更加精确、高效地建立PSDM。     

Prediction of spatially distributed seismic demands in specific structures: Ground motion and structural response

The efficacy of various ground motion intensity measures (IMs) in the prediction of spatially distributed seismic demands (engineering demand parameters, (EDPs)) within a structure is investigated. This has direct implications to building‐specific seismic loss estimation, where the seismic demand on different components is dependent on the location of the component in the structure. Several common IMs are investigated in terms of their ability to predict the spatially distributed demands in a 10‐storey office building, which is measured in terms of maximum interstorey drift ratios and maximum floor accelerations. It is found that the ability of an IM to efficiently predict a specific EDP depends on the similarity between the frequency range of the ground motion that controls the IM and that of the EDP. An IMs predictability has a direct effect on the median response demands for ground motions scaled to a specified probability of exceedance from a ground motion hazard curve. All of the IMs investigated were found to be insufficient with respect to at least one of magnitude, source‐to‐site distance, or epsilon when predicting all peak interstorey drifts and peak floor accelerations in a 10‐storey reinforced concrete frame structure. Careful ground motion selection and/or seismic demand modification is therefore required to predict such a spatially distributed demands without significant bias. Copyright © 2009 John Wiley & Sons, Ltd.     

抗震结构最优设防水平决策与全寿命优化设计的简化易损性分析方法

地震破坏状态概率分析是抗震结构最优设防水平决策与全寿命优化设计的关键问题之一,作者们曾根据"三水准"设防原则和地震烈度的危险性分析结果提出了二种实用分析方法.本文建立了结构破坏状态概率分析的地震易损性理论基础,提出了极限状态地震易损性和破坏状态地震易损性的概念.从地震易损性理论的角度,对二种地震破坏状态概率简化分析方法进行了深入分析,指出第1种方法是一种考虑结构宏观抗震能力不确定性的半理论半经验易损性分析方法,而第2种方法则是一种不考虑结构宏观抗震能力不确定性的简化易损性分析方法,根据地震破坏状态概率的计算结果对二种方法进行了对比分析,并提出了相关的建议.     

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.     

Fully probabilistic seismic displacement analysis of spatially distributed slopes using spatially correlated vector intensity measures

Earthquake‐induced slope displacement is an important parameter for safety evaluation and earthquake design of slope systems. Traditional probabilistic seismic hazard analysis usually focuses on evaluating slope displacement at a particular location, and it is not suitable for spatially distributed slopes over a large region. This study proposes a computationally efficient framework for fully probabilistic seismic displacement analysis of spatially distributed slope systems using spatially correlated vector intensity measures (IMs). First, a spatial cross‐correlation model for three key ground motion IMs, that is, peak ground acceleration (PGA), Arias intensity, and peak ground velocity, is developed using 2686 ground motion recordings from 11 recent earthquakes. To reduce the computational cost, Monte Carlo simulation and data reduction techniques are utilized to generate spatially correlated random fields for the vector IMs. The slope displacement hazards over the region are further quantified using empirical predictive equations. Finally, an illustrative example is presented to highlight the importance of the spatial correlation and the advantage of using spatially correlated vector IMs in seismic hazard analysis of spatially distributed slopes. Copyright © 2013 John Wiley & Sons, Ltd.     

Vulnerability of floating tunnel shafts for increasing earthquake loading

Fragility curves constitute the cornerstone in seismic risk evaluations and performance-based earthquake engineering. They describe the probability of a structure to experience a certain damage level for a given earthquake intensity measure, providing a relationship between seismic hazard and vulnerability. In this paper a numerical approach is applied to derive fragility curves for tunnel shafts built in clays, a component that is found in several critical infrastructure such as urban metro networks, airport facilities or water and waste water projects. The seismic response of a representative tunnel shaft is assessed using tridimensional finite difference non-linear analyses carried out with the program FLAC^3D, under increasing levels of seismic intensity. A hysteretic model is used to simulate the soil non-linear behavior during the seismic event. The effect of soil conditions and ground motion characteristics on the soil-structure system response is accounted for in the analyses. The damage is defined based on the exceedance of the concrete wall shaft capacity due to the developed seismic forces. The fragility curves are estimated in terms of peak ground acceleration at a rock or stiff soil outcrop, based on the evolution of damage with increasing earthquake intensity. The proposed fragility models allows the characterization of the seismic risk of a representative tunnel shaft typology and soil conditions considering the associated uncertainties, and partially fill the gap of data required in performing a risk analysis assessment of tunnels shafts.     

Seismic fragility of base-isolated water storage tanks under non-stationary earthquakes

Seismic fragility curves for fixed-base and base-isolated liquid storage tanks are developed under non-stationary earthquakes, and their seismic performance are compared. The correlation between different earthquake intensity measure (IM) parameters and peak response quantities of the base-isolated liquid storage tanks are investigated. The failure criteria are chosen based on (1) the elastic buckling strength of the tank wall, which is defined in terms of critical base shear and critical overturning moment, and (2) in terms of the critical isolation displacement. The uncertainty involved is considered in the earthquake characteristics. Non-stationary earthquake ground motions are generated using Monte Carlo (MC) simulation. Influence of the isolator characteristic parameters and modeling approaches on the seismic fragility of the base-isolated liquid storage tanks is also investigated. Peak ground acceleration is found to be the well correlated IM parameter with the peak response quantities of the base-isolated liquid storage tanks. Substantial decrease in the seismic fragility of the base-isolated liquid storage tanks is observed as compared to the fixed-base tanks. Significant influence of the isolator characteristic parameters on the seismic fragility of the base-isolated liquid storage tanks are reported in the present study.     

Development of empirical and analytical fragility functions using kernel smoothing methods

Fragility functions that define the probabilistic relationship between structural damage and ground motion intensity are an integral part of performance‐based earthquake engineering or seismic risk analysis. This paper introduces three approaches based on kernel smoothing methods for developing analytical and empirical fragility functions. A kernel assigns a weight to each data that is inversely related to the distance between the data value and the input of the fragility function of interest. The kernel smoothing methods are, therefore, non‐parametric forms of data interpolation. These methods enable the implicit treatment of uncertainty in either or both of ground motion intensity and structural damage without making any assumption about the shape of the resulting fragility functions. They are particularly beneficial for sparse, noisy, or non‐homogeneous data sets. For illustration purposes, two types of data are considered. The first is a set of numerically simulated responses for a four‐story steel moment‐resisting frame, and the second is a set of field observations collected after the 2010 Haiti earthquake. The results demonstrate that these methods can develop continuous representations of fragility functions without specifying their functional forms and treat sparse data sets more efficiently than conventional data binning and parametric curve fitting methods. Moreover, various uncertainty analyses are conducted to address the issues of over‐fitting, bias, and confidence intervals. Copyright © 2014 John Wiley & Sons, Ltd.     

Correlation properties of integral ground-motion intensity measures from Italian strong-motion records

This study investigates the correlation properties of integral ground-motion intensity measures (IMs) from Italian strong-motion records. The considered integral IMs include 5–95% significant duration, Housner intensity, cumulative absolute velocity, and Arias intensity. Both IM spatial correlation and the correlation between different integral and amplitude-based IMs (i.e., cross-IM correlation) are addressed in this study. To this aim, a new Italian ground-motion model (GMM) with spatial correlation for integral IMs is first introduced. Based on the newly developed GMM, the empirical correlation coefficients from interevent and intraevent residuals are investigated and various analytical correlation models between integral IMs and amplitude-based IMs are proposed. The effective range parameter representing spatial correlation properties and the trend in the cross-IM correlations are compared with existing models in the literature. The variability of the effective range parameters with respect to event-specific features is also discussed. Modeling ground-motion spatial and cross-IM correlations is an important step in seismic hazard and risk assessment of spatially distributed systems. Investigating region-specific correlation properties based on Italian strong-motion records is of special interest as several correlation models have been developed based on global datasets, often lacking earthquakes in extensional regions such as Italy.     

Intensity measures for the assessment of the seismic response of buried steel pipelines

To estimate the demand of structures, investigating the correlation between engineering demand parameters and intensity measures (IMs) is of prime importance in performance-based earthquake engineering. In the present paper, the efficiency and sufficiency of some IMs for evaluating the seismic response of buried steel pipelines are investigated. Six buried pipe models with different diameter to thickness and burial depth to diameter ratios, and different soil properties are subjected to an ensemble of 30 far-field earthquake ground motion records. The records are scaled to several intensity levels and a number of incremental dynamic analyses are performed. The approach used in the analyses is finite element modeling. Pipes are modeled using shell elements while equivalent springs and dashpots are used for modeling the soil. Several ground motion intensity measures are used to investigate their efficiency and sufficiency in assessing the seismic demand and capacity of the buried steel pipelines in terms of engineering demand parameter measured by the peak axial compressive strain at the critical section of the pipe. Using the regression analysis, efficient and sufficient IMs are proposed for two groups of buried pipelines separately. The first one is a group of pipes buried in soils with low stiffness and the second one is those buried in soils with higher stiffness. It is concluded that for the first group of pipes, \(\sqrt {{\text{VSI}}[\upomega_{1} ({\text{PGD}} + {\text{RMS}}_{\text{d}} )]}\) followed by root mean square of displacement (RMS_d) are the optimal IMs based on both efficiency and sufficiency; and for the second group, the only optimal IM is PGD²/RMS_d.     

Intensity measures for response prediction of a torsional building subjected to bi-directional earthquake ground motion

In the framework of performance-based assessment of existing RC buildings, an important step is the definition of an intensity measure (IM) of the seismic action that be accurate and efficient in correlating with the engineering demand parameters (EDP) of interest and in reducing the variability in their prediction. Objective of this work is to assess five different IMs, when they are applied to an in-plan irregular three-dimensional building subjected to bi-directional earthquake motion. The objective is pursued through multiple regression of the results obtained from nonlinear dynamic analyses. The study shows that certain IMs are more suitable to evaluate EDPs when dealing with 3D structures subjected to 2D earthquakes.     

椭圆空洞对盾构隧道地震易损性曲线影响研究

采用增量动力分析方法,探究水平向地震下地层空洞对盾构隧道地震响应特征的影响规律。针对管片损伤及周边地层应力,选择弯矩比作为性能评价指标,峰值加速度(PGA)及峰值速度(PGV)作为衡量地震强度指标(IM),阐明椭圆形空洞对管片抗震性能的影响,得到隧道结构的地震易损性曲线。研究表明：椭圆空洞加大了浅埋盾构隧道的地震破坏概率;PGA与PGV均可作为IM并获得相应的隧道易损性曲线;使用弯矩比作为破坏指标,PGV作为地震动指标,其对应的易损性曲线对地层变异性更敏感。研究结论可为潜在空洞发育区防震方案的制定提供参考。     

Ground-motion intensity measure correlations observed in Italian strong-motion records

Ground-motion models (GMMs) are widely used in probabilistic seismic hazard analysis (PSHA) to estimate the probability distributions of earthquake-induced ground-motion intensity measures (IMs) at a site, given an earthquake of a certain magnitude occurring at a nearby location. Accounting for spatial and cross-IM correlations in earthquake-induced ground motions has important implications on probabilistic seismic hazard and loss estimates. This study first develops a new Italian GMM with spatial correlation for 31 amplitude-related IMs, including peak ground acceleration (PGA), peak ground velocity (PGV), and 5%-damped elastic pseudo-spectral accelerations (PSAs) at 29 periods ranging from 0.01 to 4 seconds. The model estimation is performed through a recently developed one-stage nonlinear regression algorithm proposed by the authors, known as the Scoring estimation approach. In fact, current state-of-practice approaches estimate spatial correlation separately from the GMM estimation, resulting in inconsistent and statistically inefficient estimators of interevent and intraevent variances and parameters in the spatial correlation model. We test whether this affects the subsequent cross-IM correlation analysis. To this aim, based on the newly developed GMM, the empirical correlation coefficients from interevent and intraevent residuals are investigated. Finally, a set of analytical correlation models between the selected IMs are proposed. This is of special interest as several correlation models between different IMs have been calibrated and validated based on advanced GMMs and global datasets, lacking earthquakes in extensional regions; however, modeling the correlation between different IM types has not been adequately addressed by current, state-of-the-art GMMs and recent ground-motion records for Italy.     

Incremental dynamic analysis of wood-frame houses in Canada: Effects of dominant earthquake scenarios on seismic fragility

Seismic fragility can be assessed by conducting incremental dynamic analysis (IDA). This study extends the current conditional mean spectrum (CMS)-based record selection approach for IDA by taking into account detailed seismic hazard information. The proposed method is applied to conventional wood-frame houses in Canada, across which dominant earthquake scenarios and associated hazard levels vary significantly. Effects due to different seismic environments, site conditions, CMS-based record selection methods, and house models are investigated by comparing various seismic fragility models. Moreover, relative impact of the key characteristics is evaluated in terms of seismic loss curve for a group of wood-frame houses. Importantly, a close examination of regional seismic hazard characteristics using seismic hazard curve and seismic deaggregation facilitates the deeper understanding of the impact of ground motion characteristics on seismic fragility. A comprehensive and systematic assessment of key uncertainties associated with seismic fragility is provided.     

Deriving fragility functions from bilinearized capacity curves for earthquake scenario modelling using the conditional spectrum

The development of fragility curves to perform seismic scenario-based risk assessment requires a fully probabilistic procedure in order to account for uncertainties at each step of the computation. This is especially true when developing fragility curves conditional on an Intensity Measure that is directly available from a ground-motion prediction equation. In this study, we propose a new derivation method that uses realistic spectra instead of design spectral shapes or uniform hazard spectra and allows one to easily account for the features of the site-specific hazard that influences the fragility, without using non-linear dynamic analysis. The proposed method has been applied to typical school building types in the city of Basel (Switzerland) and the results have been compared to the standard practice in Europe. The results confirm that fragility curves are scenario dependent and are particularly sensitive to the magnitude of the earthquake scenario. The same background theory used for the derivation of the fragility curves has allowed an innovative method to be proposed for the conversion of fragility curves to a common IM (i.e. spectral acceleration or PGA). This conversion is the only way direct comparisons of fragility curves can be made and is useful when inter-period correlation cannot be used in scenario loss assessment. Moreover, such conversion is necessary to compare and verify newly developed curves against those from previous studies. Conversion to macroseismic intensity is also relevant for the comparison between mechanical-based and empirical fragility curves, in order to detect possible biases.