Statistical procedures for developing earthquake damage fragility curves

This paper describes statistical procedures for developing earthquake damage fragility functions. Although fragility curves abound in earthquake engineering and risk assessment literature, the focus has generally been on the methods for obtaining the damage data (i.e., the analysis of structures), and little emphasis is placed on the process for fitting fragility curves to this data. This paper provides a synthesis of the most commonly used methods for fitting fragility curves and highlights some of their significant limitations. More novel methods are described for parametric fragility curve development (generalized linear models and cumulative link models) and non‐parametric curves (generalized additive model and Gaussian kernel smoothing). An extensive discussion of the advantages and disadvantages of each method is provided, as well as examples using both empirical and analytical data. The paper further proposes methods for treating the uncertainty in intensity measure, an issue common with empirical data. Finally, the paper describes approaches for choosing among various fragility models, based on an evaluation of prediction error for a user‐defined loss function. Copyright © 2015 John Wiley & Sons, Ltd.     

Development and assessment of damage‐to‐loss models for moment‐frame reinforced concrete buildings

The assessment of earthquake loss often requires the definition of a relation between a measure of damage and a quantity of loss, usually achieved through the employment of a damage‐to‐loss model. These models are frequently characterized by a large variability, which inevitably increases the uncertainty in the vulnerability assessment and earthquake loss estimation. This study provides an insight on the development of damage‐to‐loss functions for moment‐frame reinforced concrete buildings through an analytical methodology. Tri‐dimensional finite element models of existing reinforced concrete buildings were subjected to a number of ground motion records compatible with the seismicity in the region of interest, through nonlinear dynamic analysis. These results were used to assess, for a number of damage states, the probability distribution of loss ratio, taking into consideration member damage and different repair techniques, as well as to derive sets of fragility functions. Then, a vulnerability model (in terms of the ratio of cost of repair to cost of replacement, conditional on the level of ground shaking intensity) was derived and compared with the vulnerability functions obtained through the combination of various damage‐to‐loss models with the set of fragility functions developed herein. In order to provide realistic estimates of economic losses due to seismic action, a comprehensive study on repair costs using current Portuguese market values was also carried out. The results of this study highlight important issues in the derivation of vulnerability functions, which are a fundamental component for an adequate seismic risk assessment. © 2015 The Authors. Earthquake Engineering & Structural Dynamics published by John Wiley & Sons Ltd.     

On the treatment of uncertainty in seismic vulnerability and portfolio risk assessment

In a related study developed by the authors, building fragility is represented by intensity‐specific distributions of damage exceedance probability of various damage states. The contribution of the latter has been demonstrated in the context of loss estimation of building portfolios, where it is shown that the proposed concept of conditional fragility functions provides the link between seismic intensity and the uncertainty in damage exceedance probabilities. In the present study, this methodology is extended to the definition of building vulnerability, whereby vulnerability functions are characterized by hazard‐consistent distributions of damage ratio per level of primary seismic intensity parameter—Sa(T₁). The latter is further included in a loss assessment framework, in which the impact of variability and spatial correlation of damage ratio in the probabilistic evaluation of seismic loss is accounted for, using test‐bed portfolios of 2, 5, and 8‐story precode reinforced concrete buildings located in the district of Lisbon, Portugal. This methodology is evaluated in comparison with current state‐of‐the‐art methods of vulnerability and loss calculation, highlighting the discrepancies that can arise in loss estimates when the variability and spatial distributions of damage ratio, influenced by ground motion properties other than the considered primary intensity measure, are not taken into account.     

On the treatment of uncertainties in the development of fragility functions for earthquake loss estimation of building portfolios

State‐of‐the‐art methods for the assessment of building fragility consider the structural capacity and seismic demand variability in the estimation of the probability of exceeding different damage states. However, questions remain regarding the appropriate treatment of such sources of uncertainty from a statistical significance perspective. In this study, material, geometrical and mechanical properties of a number of building classes are simulated by means of a Monte Carlo sampling process in which the statistical distribution of the aforementioned parameters is taken into consideration. Record selection is performed in accordance with hazard‐consistent distributions of a comprehensive set of intensity measures, and issues related with sufficiency, efficiency, predictability and scaling robustness are addressed. Based on the appraised minimum number of ground motion records required to achieve statistically meaningful estimates of response variability conditioned on different levels of seismic intensity, the concept of conditional fragility functions is presented. These functions translate the probability of exceeding a set of damage states as a function of a secondary sufficient intensity measure, when records are selected and scaled for a particular level of primary seismic intensity parameter. It is demonstrated that this process allows a hazard‐consistent and statistically meaningful representation of uncertainty and correlation in the estimation of intensity‐dependent damage exceedance probabilities. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of earthquake ground motions on fragility curves of highway bridge piers based on numerical simulation

Fragility curves express the probability of structural damage due to earthquakes as a function of ground motion indices, e.g., PGA, PGV. Based on the actual damage data of highway bridges from the 1995 Hyogoken‐Nanbu (Kobe) earthquake, a set of empirical fragility curves was constructed. However, the type of structure, structural performance (static and dynamic) and variation of input ground motion were not considered to construct the empirical fragility curves. In this study, an analytical approach was adopted to construct fragility curves for highway bridge piers of specific bridges. A typical bridge structure was considered and its piers were designed according to the seismic design codes in Japan. Using the strong motion records from Japan and the United States, non‐linear dynamic response analyses were performed, and the damage indices for the bridge piers were obtained. Using the damage indices and ground motion indices, fragility curves for the bridge piers were constructed assuming a lognormal distribution. The analytical fragility curves were compared with the empirical ones. The proposed approach may be used in constructing the fragility curves for highway bridge structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Seismic reliability assessment and the nonergodicity in the modelling parameter uncertainties

Modelling uncertainty can significantly affect the structural seismic reliability assessment. However, the limit state excursion due to this type of uncertainty may not be described by a Poisson process as it lacks renewal properties with the occurrence of each earthquake event. Furthermore, considering uncertainties related to ground motion representation by employing recorded ground motions together with modelling uncertainties is not a trivial task. Robust fragility assessment, proposed previously by the authors, employs the structural response to recorded ground motion as data in order to update prescribed seismic fragility models. Robust fragility can be extremely efficient for considering also the structural modelling uncertainties by creating a dataset of one-to-one assignments of structural model realizations and as-recorded ground motions. This can reduce the computational effort by more than 1 order of magnitude. However, it should be kept in mind that the fragility concept itself is based on the underlying assumption of Poisson-type renewal. Using the concept of updated robust reliability, considering both the uncertainty in ground motion representation based on as-recorded ground motion and non ergodic modelling uncertainties, the error introduced through structural reliability assessment by using the robust fragility is quantified. It is shown through specific application to an existing RC frame that this error is quite small when the product of the time interval and the standard deviation of failure rate is small and is on the conservative side.     

Seismic risk assessment of liquid storage tanks via a nonlinear surrogate model

A performance‐based earthquake engineering approach is developed for the seismic risk assessment of fixed‐roof atmospheric steel liquid storage tanks. The proposed method is based on a surrogate single‐mass model that consists of elastic beam‐column elements and nonlinear springs. Appropriate component and system‐level damage states are defined, following the identification of commonly observed modes of failure that may occur during an earthquake. Incremental dynamic analysis and simplified cloud are offered as potential approaches to derive the distribution of response parameters given the seismic intensity. A parametric investigation that engages the aforementioned analysis methods is conducted on 3 tanks of varying geometry, considering both anchored and unanchored support conditions. Special attention is paid to the elephant's foot buckling formation, by offering extensive information on its capacity and demand representation within the seismic risk assessment process. Seismic fragility curves are initially extracted for the component‐level damage states, to compare the effect of each analysis approach on the estimated performance. The subsequent generation of system‐level fragility curves reveals the issue of nonsequential damage states, whereby significant damage may abruptly appear without precursory lighter damage states.     

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

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

Wavelet‐based characterization of design ground motions

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

基于性能的钢筋混凝土框架剪力墙结构地震易损性分析

以某典型的20层钢筋混凝土框架剪力墙结构作为研究对象,研究基于性能的RC框架剪力墙结构易损性分析方法。首先选择合适的地震动记录,以0.2g为步长进行调幅后,建立300个结构-地震动样本空间,并确定结构损伤指标和性能参数;然后应用增量动力分析方法计算结构的地震动力响应,选择基本周期加速度反应谱为地震动参数,以研究结构反应的不确定性,并深入分析地震动参数与结构地震需求参数的关系;在此基础上,建立该结构基于加速度反应谱的易损性曲线进行结构易损性分析与评估。结果表明：随着地震动强度的增大,IDA曲线由单调增加变为非单调增加,分位曲线（16%,50%和84%）可以准确地衡量结构的性能;框剪结构在地震作用下的抗震性能表现良好,随着地震强度的增长,各性能超越概率大小的增长速度是不同的。     

Development of seismic fragility surfaces for reinforced concrete buildings by means of nonlinear time‐history analysis

Fragility curves are generally developed using a single parameter to relate the level of shaking to the expected structural damage. The main goal of this work is to use several parameters to characterize the earthquake ground motion. The fragility curves will, therefore, become surfaces when the ground motion is represented by two parameters. To this end, the roles of various strong‐motion parameters on the induced damage in the structure are compared through nonlinear time‐history numerical calculations. A robust structural model that can be used to perform numerous nonlinear dynamic calculations, with an acceptable cost, is adopted. The developed model is based on the use of structural elements with concentrated nonlinear damage mechanics and plasticity‐type behavior. The relations between numerous ground‐motion parameters, characterizing different aspects of the shaking, and the computed damage are analyzed and discussed. Natural and synthetic accelerograms were chosen/computed based on a consideration of the magnitude‐distance ranges of design earthquakes. A complete methodology for building fragility surfaces based on the damage calculation through nonlinear numerical analysis of multi‐degree‐of‐freedom systems is proposed. The fragility surfaces are built to represent the probability that a given damage level is reached (or exceeded) for any given level of ground motion characterized by the two chosen parameters. The results show that an increase from one to two ground‐motion parameters leads to a significant reduction in the scatter in the fragility analysis and allows the uncertainties related to the effect of the second ground‐motion parameter to be accounted for within risk assessments. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of peak ground velocity on deformation demands for SDOF systems

The effect of peak ground velocity (PGV) on single‐degree‐of‐freedom (SDOF) deformation demands and for certain ground‐motion features is described by using a total of 60 soil site records with source‐to‐site distances less than 23 km and moment magnitudes between 5.5 and 7.6. The observations based on these records indicate that PGV correlates well with the earthquake magnitude and provides useful information about the ground‐motion frequency content and strong‐motion duration that can play a role on the seismic demand of structures. The statistical results computed from non‐linear response history analyses of different hysteretic models highlight that PGV correlates better with the deformation demands with respect to other ground motion intensity measures. The choice of PGV as ground motion intensity decreases the dispersion due to record‐to‐record variability of SDOF deformation demands, particularly in the short period range. The central tendencies of deformation demands are sensitive to PGV and they may vary considerably as a function of the hysteretic model and structural period. The results provided in this study suggest a consideration of PGV as a stable candidate for ground motion intensity measure in simplified seismic assessment methods that are used to estimate structural performance for earthquake hazard analysis. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐structure‐specific intensity measure parameters and characteristic period of near‐fault ground motions

This paper focuses on the effects of long‐period pulse of near‐fault ground motions on the structural damage potential. Two sets of near‐fault ground motion records from Chi‐Chi, Taiwan earthquake and Northridge earthquake with and without distinct pulse are selected as the input, and the correlation analysis between 30 non‐structure‐specific intensity measure parameters and maximum inelastic displacements and energy responses (input energy and hysteretic energy) of bilinear single degree of freedom systems are conducted. Based on the frequency characteristic of near‐fault ground motions with remarkable long‐period components, two intensity indices are proposed, namely, the improved effective peak acceleration (IEPA) and improved effective peak velocity (IEPV). In addition a new characteristic period of these ground motions is defined based on IEPA and IEPV. Numerical results illustrate that the intensity measure parameters related to ground acceleration present the best correlation with the seismic responses for rigid systems; the velocity‐related and displacement‐related parameters are better for medium‐frequency systems and flexible systems, respectively. The correlation curves of near‐fault ground motions with velocity pulse differ from those of ground motions without pulse. Moreover, the improved parameters IEPA and IEPV of near‐fault impulsive ground motions enhance the performance of intensity measure of corresponding conventional parameters, i.e. EPA and EPV. The new characteristic period based on IEPA and IEPV can better reflect the frequency content of near‐fault ground motions. Copyright © 2009 John Wiley & Sons, Ltd.     

Aftershock collapse fragility curves for non‐ductile RC buildings: a scenario‐based assessment

Recent studies have addressed the computation of fragility curves for mainshock (MS)‐damaged buildings. However, aftershock (AS) fragilities are generally conditioned on a range of potential post‐MS damage states that are simulated via static or dynamic analyses performed on an intact building. Moreover, there are very few cases where the behavior of non‐ductile reinforced concrete buildings is analyzed. This paper presents an evaluation of AS collapse fragility conditioned on various return periods of MSs, allowing for the rapid assessment of post‐earthquake safety variations based solely on the intensity of the damaging earthquake event. A refined multi‐degree‐of‐freedom model of a seven‐storey non‐ductile building, which includes brittle failure simulations and the evaluation of a system level collapse, is adopted. Aftershock fragilities are obtained by performing an incremental dynamic analysis for a number of MS–AS ground motion sequences and a variety of MS intensities. The AS fragilities show that the probability of collapse significantly increases for higher return periods for the MS. However, this result is mainly ascribable to collapses occurred during MSs. When collapse cases that occur during a MS are not considered in the assessment of AS collapse probability, a smaller shift in the fragility curves is observed as the MS intensity increases. This result is justified considering the type of model and collapse modes introduced, which strongly depend on the brittle behavior of columns failing in shear or due to axial loads. The analysis of damage that is due to MSs when varying the return period confirms this observation. Copyright © 2017 John Wiley & Sons, Ltd.     

Estimation of strong motion distribution in the 1995 Kobe earthquake based on building damage data

The 1995 Kobe earthquake caused unprecedented damage to buildings and civil infrastructures in the city of Kobe and its surrounding areas. In order to evaluate the structural damage in this area due to the earthquake, it is important to estimate the distribution of earthquake ground motion. However, since the number of strong ground motion records is not enough in the heavily damaged areas, it is necessary to estimate the distribution using other data sources. In this paper, the fragility curves for low‐rise residential buildings were constructed using the recorded motions and the building damage data from the intensive field survey by the AIJ and CPIJ group. The fragility curves obtained were then employed to estimate the strong motion distribution in the district level for Kobe and the surrounding areas during the earthquake. The results may be useful to investigate the various damages caused by the earthquake. Copyright © 2001 John Wiley & Sons, Ltd.     

Development of fragility functions as a damage classification/prediction method for steel moment‐resisting frames using a wavelet‐based damage sensitive feature

Fragility functions are commonly used in performance‐based earthquake engineering for predicting the damage state of a structure subjected to an earthquake. This process often involves estimating the structural damage as a function of structural response, such as the story drift ratio and the peak floor absolute acceleration. In this paper, a new framework is proposed to develop fragility functions to be used as a damage classification/prediction method for steel structures based on a wavelet‐based damage sensitive feature (DSF). DSFs are often used in structural health monitoring as an indicator of the damage state of the structure, and they are easily estimated from recorded structural responses. The proposed framework for damage classification of steel structures subjected to earthquakes is demonstrated and validated with a set of numerically simulated data for a four‐story steel moment‐resisting frame designed based on current seismic provisions. It is shown that the damage state of the frame is predicted with less variance using the fragility functions derived from the wavelet‐based DSF than it is with fragility functions derived from an alternate acceleration‐based measure, the spectral acceleration at the first mode period of the structure. Therefore, the fragility functions derived from the wavelet‐based DSF can be used as a probabilistic damage classification model in the field of structural health monitoring and an alternative damage prediction model in the field of performance‐based earthquake engineering. Copyright © 2011 John Wiley & Sons, Ltd.     

Ground motion scaling methods for different site conditions and structure characteristics

Non‐linear dynamic time‐history analyses conducted as part of a performance‐based seismic design approach often require that the ground motion records are scaled to a specified level of seismic intensity. Recent research has demonstrated that certain ground motion scaling methods can introduce a large scatter in the estimated seismic demands. The resulting demand estimates may be biased, leading to designs with significant uncertainty and unknown margins of safety, unless a relatively large ensemble of ground motion records is used. This paper investigates the effectiveness of seven ground motion scaling methods in reducing the scatter in estimated peak lateral displacement demands. Non‐linear single‐degree‐of‐freedom systems and non‐linear multi‐degree‐of‐freedom systems are considered with different site conditions (site soil profile and epicentral distance) and structural characteristics (yield strength, period, and hysteretic behavior). It is shown that scaling methods that work well for ground motions representative of stiff soil and far‐field conditions lose their effectiveness for soft soil and near‐field conditions for a wide range of structural characteristics. Copyright © 2003 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.     

Vector-valued fragility functions for seismic risk evaluation

This article presents a method for the development of vector-valued fragility functions, which are a function of more than one intensity measure (IM, also known as ground-motion parameters) for use within seismic risk evaluation of buildings. As an example, a simple unreinforced masonry structure is modelled using state-of-the-art software and hundreds of nonlinear time-history analyses are conducted to compute the response of this structure to earthquake loading. Dozens of different IMs (e.g. peak ground acceleration and velocity, response spectral accelerations at various periods, Arias intensity and various duration and number of cycle measures) are considered to characterize the earthquake shaking. It is demonstrated through various statistical techniques (including Receiver Operating Characteristic analysis) that the use of more than one IM leads to a better prediction of the damage state of the building than just a single IM, which is the current practice. In addition, it is shown that the assumption of the lognormal distribution for the derivation of fragility functions leads to more robust functions than logistic, log-logistic or kernel regression. Finally, actual fragility surfaces using two pairs of IMs (one pair are uncorrelated while the other are correlated) are derived and compared to scalar-based fragility curves using only a single IM and a significant reduction in the uncertainty of the predicted damage level is observed. This type of fragility surface would be a key component of future risk evaluations that take account of recent developments in seismic hazard assessment, such as vector-valued probabilistic seismic hazard assessments.     

Seismic fragility of traction elevators

Second‐generation performance‐based earthquake engineering (PBEE‐2) requires a library of component fragility functions to estimate probabilistic damage to a wide variety of building components. The present work draws on a large body of (mostly) post‐earthquake reconnaissance and (some) post‐earthquake survey observations of traction elevators to create fragility functions useful in PBEE‐2. Two surveys provide detailed observations of 115 representative elevators at 12 hospitals shaken in the 1989 Loma Prieta and 1994 Northridge earthquakes and selected without regard to or foreknowledge of damage. Of these, 55 failed and 60 did not. Approximately half were installed after an important code change of 1972, so one can distinguish the performance of pre‐1973 and post‐1973 elevator construction. They experienced a range of strong motion: 22 with peak ground acceleration (PGA) < 0.25 g, 93 with 0.25 g < PGA < 0.85 g. The hospitals had elevator failure rates as low as 0% and as high as 100%. A third survey describes damage qualitatively for six sites with PGA ≤ 0.25 and per‐site failure rates of 0% to perhaps 30%. Fragility functions are offered where the damage state is the loss of functionality of the elevator. The elevators in these surveys exhibit a median capacity of PGA ≈ 0.35 g with a logarithmic standard deviation of 0.40. Capacity is modestly sensitive to whether the elevator was installed before or after 1973. Using building‐specific intensity measures such as S_a(T₁) does not improve the fragility functions. Copyright © 2015 John Wiley & Sons, Ltd.