Intensity measure selection for vulnerability studies of building classes

The selection of a scalar Intensity Measure (IM) for performing analytical vulnerability (loss) assessment across a building class is addressed. We investigate the ability of several IM choices to downgrade the effect of seismological parameters (sufficiency) as well as reduce the record‐to‐record variability (efficiency) for both highrise and lowrise sets of ‘index’ buildings. These characteristics are explored in unprecedented detail, employing comparisons and statistical significance testing at given levels of local engineering demand parameters (story drift ratios and peak floor accelerations) that relate to losses, instead of global variables such as the maximum interstory drift. Thus, a detailed limit‐state‐specific view is offered for the suitability of different scalar IMs for loss assessment. As expected, typical single‐period spectral values are found to introduce unwanted bias at high levels of scaling, both for a single as well as a class of buildings. On the other hand, the geometric mean of the spectral acceleration values estimated at several periods between the class‐average second‐mode and an elongated class‐average first‐mode period offers a practical choice that significantly reduces the spectral‐shape bias without requiring the development of new ground motion prediction equations. Given that record selection remains a site‐ and building‐specific process, such an improved IM can help achieve reliable estimates for building portfolios, as well as single structures, at no additional cost. Copyright © 2015 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.     

Evaluation of a recently proposed record selection and scaling procedure for low‐rise to mid‐rise reinforced concrete buildings and its use for probabilistic risk assessment studies

This paper evaluates a recent record selection and scaling procedure of the authors that can determine the probabilistic structural response of buildings behaving either in the elastic or post‐elastic range. This feature marks a significant strength on the procedure as the probabilistic structural response distribution conveys important information on probability‐based damage assessment. The paper presents case studies that show the utilization of the proposed record selection and scaling procedure as a tool for the estimation of damage states and derivation of site‐specific and region‐specific fragility functions. The method can be used to describe exceedance probabilities of damage limits under a certain target hazard level with known annual exceedance rate (via probabilistic seismic hazard assessment). Thus, the resulting fragility models can relate the seismicity of the region (or a site) with the resulting building performance in a more accurate manner. Under this context, this simple and computationally efficient record selection and scaling procedure can be benefitted significantly by probability‐based risk assessment methods that have started to be considered as indispensable for developing robust earthquake loss models. Copyright © 2013 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.     

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.     

Record selection for aftershock incremental dynamic analysis

The back‐to‐back application of mainshock records as aftershock is often considered in conducting aftershock incremental dynamic analysis. In such an approach, the characteristics of mainshock records are considered to be similar to those of major aftershock records within the same mainshock–aftershock sequences. The underlying assumption is that the characteristics of selected mainshocks, other than those used for record selection, are not significant in the assessment of structural responses. A case study is set up to investigate the effects of aftershock record selection on the collapse vulnerability assessment. The numerical results for a specific wood‐frame structure indicate that the aftershock fragility can be affected by the aftershock record characteristics, particularly response spectral shape. Copyright © 2014 John Wiley & Sons, Ltd.     

Conditional spectrum bidirectional record selection for risk assessment of 3D structures using scalar and vector IMs

Several proposals are explored for the hazard and intensity measure (IM) consistent selection of bidirectional ground motions to assess the performance of 3D structural models. Recent studies have shown the necessity of selecting records that thoroughly represent the seismicity at the site of interest, as well as the usefulness of efficient IMs capable of estimating the response of buildings with low scatter. However, the advances realized are mostly geared towards the structural analysis of 2D models. Few are the combined record, and IM selection approaches suggested expressly for nonlinear dynamic analysis of 3D structural models, especially when plan asymmetry and torsion sensitivity come into play. Conditional spectrum selection is leveraged and expanded here to offer a suite of approaches based on both scalar and vector IMs that convey information from two orthogonal horizontal components of the ground motion. Applications on multiple 3D building models highlight the importance of (a) employing the same IM for both record selection and response assessment and (b) maintaining hazard consistency in both horizontal components, when using either a scalar or a vector IM. All tested approaches that respect these conditions can be viable, yet the one based on the geometric mean of multiple spectral ordinates from both components over a period range seems to hold the most promise for general use.     

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.     

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.     

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.     

Seismic vulnerability assessment to slight damage based on experimental modal parameters

The aim of this paper is to adjust behaviour models for each class of structure for vulnerability assessment by using ambient vibration. A simple model based on frequencies, mode shapes and damping, taken from ambient vibrations, allows computation of the response of the structures and comparison of inter‐storey drifts with the limits found in the literature for the slight damage grade, considered here as the limit of elastic behaviour. Two complete methodologies for building fragility curves are proposed: (1) using a multi‐degree of freedom system including higher modes and full seismic ground‐motion and (2) using a single‐degree of freedom model considering the fundamental mode f₀ of the structure and ground‐motion displacement response spectra S_D(f₀). These two methods were applied to the city of Grenoble, where 60 buildings were studied. Fragility curves for slight damage were derived for the various masonry and reinforced concrete classes of buildings. A site‐specific earthquake scenario, taking into account local site conditions, was considered, corresponding to an M_L = 5.5 earthquake at a distance of 15 km. The results show the benefits of using experimental models to reduce variability of the slight damage fragility curve. Moreover, by introducing the experimental modal model of the buildings, it is possible to improve seismic risk assessment at an overall scale (the city) or a local scale (the building) for the first damage grade (slight damage). This level of damage, of great interest for moderate seismic‐prone regions, may contribute to the seismic loss assessment. Copyright © 2011 John Wiley & Sons, Ltd.     

Methodology for the development of analytical fragility curves for retrofitted bridges

Fragility curves for retrofitted bridges indicate the influence of various retrofit measures on the probability of achieving specified levels of damage. This paper presents an analytical methodology for developing fragility curves for classes of retrofitted bridge systems. The approach captures the impact of retrofit on the vulnerability of multiple components, which to date has not been adequately addressed, and results in a comparison of the system fragility before and after the application of different retrofit measures. Details presented include analytical modeling, uncertainty treatment, impact of retrofit on demand models, capacity estimates, and component and system fragility curves. The findings indicate the importance of evaluating the impact of retrofit not only on the targeted response quantity and component vulnerability but also on the overall bridge fragility. As illustrated by the case study of a retrofitted multi‐span continuous (MSC) concrete girder bridge class, a given retrofit measure may have a positive impact on some components, yet no impact or a negative impact on other critical components. Consideration of the fragility based only on individual retrofitted components, without regard for the system, may lead to over‐estimation or under‐estimation of the impact on the bridge fragility. The proposed methodology provides an opportunity to effectively compare the fragility of the MSC concrete bridge retrofit with a range of different retrofit measures. The most effective retrofit in reducing probable damage for a given intensity is a function of the damage state of interest. Copyright © 2008 John Wiley & Sons, Ltd.     

Mainshock-consistent ground motion record selection for aftershock sequences

In recent years, the additional risk posed to the built environment due to aftershock sequences and triggered events has been brought to attention, and several efforts have been directed towards developing fragility functions for structures in damaged conditions. Despite this rise of interest, a rather fundamental component for such tasks, namely that of aftershock ground motion record selection, has remained under-scrutinized. Herein, we propose a pragmatic procedure that can be applied for the selection of mainshock-aftershock ground motion pairs using consistent causal parameters and accounting for the correlation between their spectral accelerations. In addition, a structural analysis strategy that can be employed for the analytical derivation of damage-dependent fragility functions is outlined and presented through a case study. A more conventional back-to-back IDA analysis is also carried out in order to compare the derived damage-dependent fragility functions with the ones obtained with the proposed procedure. The results indicate that record selection remains a crucial factor even when assessing the structural vulnerability of damaged buildings, and should thus be treated cautiously.     

A methodology for defining seismic scenario‐structure‐based limit state criteria for rc high‐rise wall buildings using net drift

As a result of population growth and consequent urbanization, the number of high‐rise buildings is rapidly growing worldwide resulting in increased exposure to multiple‐scenario earthquakes and associated risk. The wide range in frequency contents of possible strong ground motions can have an impact on the seismic response, vulnerability and limit states definitions of RC high‐rise wall structures. Motivated by the pressing need to derive more accurate fragility relations to be used in seismic risk assessment and mitigation of such structures, a methodology is proposed to obtain reliable, Seismic Scenario‐Structure‐Based (SSSB) definitions of limit state criteria. A 30‐story wall building, located in a multi‐seismic scenario study region, is utilized to illustrate the methodology. The building is designed following modern codes and then modeled using nonlinear fiber‐based approach. Uncertainty in ground motions is accounted for by the selection of forty real earthquake records representing two seismic scenarios. Seismic scenario‐based building local response at increasing earthquake intensities is mapped using Multi‐Record Incremental Dynamic Analyses (MRIDAs) with a new scalar intensity measure. Net Inter‐Story Drift (NISD) is selected as a global damage measure based on a parametric study involving seven buildings ranging from 20 to 50 stories. This damage measure is used to link local damage events, including shear, to global response under different seismic scenarios. While the study concludes by proposing SSSB limit state criteria for the sample building, the proposed methodology arrives at a reliable definition of limit state criteria for an inventory of RC high‐rise wall buildings under multiple earthquake scenarios. Copyright © 2016 John Wiley & Sons, Ltd.     

Loss estimation for non‐ductile reinforced concrete building in Victoria,British Columbia,Canada: effects of mega‐thrust Mw9‐class subduction earthquakes and aftershocks

This paper presents, within the performance‐based earthquake engineering framework, a comprehensive probabilistic seismic loss estimation method that accounts for main sources of uncertainty related to hazard, vulnerability, and loss. The loss assessment rigorously integrates multiple engineering demand parameters (maximum and residual inter‐story drift ratio and peak floor acceleration) with consideration of mainshock–aftershock sequences. A 4‐story non‐ductile reinforced concrete building located in Victoria, British Colombia, Canada, is considered as a case study. For 100 mainshock and mainshock–aftershock earthquake records, incremental dynamic analysis is performed, and the three engineering demand parameters are fitted with a probability distribution and corresponding dependence computed. Finally, with consideration of different demolition limit states, loss assessment is performed. From the results, it can be shown that when seismic vulnerability models are integrated with seismic hazard, the aftershock effects are relatively minor in terms of overall seismic loss (1–4% increase). Moreover, demolition limit state parameters, uncertainties of collapse fragility, and non‐collapse seismic demand prediction models have showed significant contribution to the loss assessment. The seismic loss curves for the reference case and for cases with the varied parameters can differ by as large as about 150%. Copyright © 2015 John Wiley & Sons, Ltd.     

损失分布和非结构损失对房屋建筑地易损性的影响

房屋建筑的地震易损性是地震损失评估和地震巨灾风险模型的基础。作为房屋建筑的重要组成部分,各类非结构构件的损失在现有的易损性模型中并未得到足够重视。本文以一栋典型钢筋混凝土框架结构教学楼为对象,通过将房屋建筑中的各类构件划分为具有不同地震损伤特性和损失后果的易损性组,考察建筑内的损失分布和非结构损失对房屋建筑地震易损性的影响。分析结果表明:由于许多非结构构件在中小地震作用下即可能发生较严重的破坏,房屋建筑在中小地震下的易损性主要受非结构损失控制;随着地震动强度等级的不断提高,结构损伤渐趋严重,结构损失对整体建筑易损性的影响不断增大;在结构进入震后不可修状态之前,建筑不同楼层的损失分布是评估建筑地震损失时不可忽略的因素。     

A hybrid method for the vulnerability assessment of R/C and URM buildings

The methodology followed by the Aristotle University (AUTh) team for the vulnerability assessment of reinforced concrete (R/C) and unreinforced masonry (URM) structures is presented. The paper focuses on the derivation of vulnerability (fragility) curves in terms of peak ground acceleration (PGA), as well as spectral displacement (s _d), and also includes the estimation of capacity curves, for several R/C and URM building types. The vulnerability assessment methodology is based on the hybrid approach developed at AUTh, which combines statistical data with appropriately processed (utilising repair cost models) results from nonlinear dynamic or static analyses, that permit extrapolation of statistical data to PGA’s and/or spectral displacements for which no data are available. The statistical data used herein are from earthquake-damaged greek buildings. An extensive numerical study is carried out, wherein a large number of building types (representing most of the common typologies in S. Europe) are modelled and analysed. Vulnerability curves for several damage states are then derived using the aforementioned hybrid approach. These curves are subsequently used in combination with the mean spectrum of the Microzonation study of Thessaloniki as the basis for the derivation of new vulnerability curves involving spectral quantities. Pushover curves are derived for all building types, then reduced to standard capacity curves, and can easily be used together with the S _d fragility curves as an alternative for developing seismic risk scenarios.     

Post-earthquake estimates of different ground motion intensity measures for the 2012 Emilia earthquake

Estimates of the earthquake ground motion intensity over a geographical area have multiple uses, that is, emergency management, civil protection and seismic fragility assessment. In particular, with reference to fragility assessment, it is of interest to have estimates of the values of different ground-motion intensity measures in order to correlate them with the observed damage. To this purpose, the present paper uses a procedure recently proposed in the literature to estimate the ground-motion intensity for the 2012 Emilia mainshocks, considering different ground motion intensity measures and directionality effects. Ground motion prediction equations based on different site effect models, and spatial correlation models are calibrated for the Emilia earthquakes. The paper discusses the accuracy of the shakemaps obtained using the different soil effect models considered and presents the obtained shakemaps as supplementary material. The procedure presented in the paper is aimed at providing ground motion intensity values for seismic fragility assessment and is not intended as a tool to estimate shakemaps for rapid emergency assessment.     

Use of the conditional spectrum to incorporate record‐to‐record variability in simplified seismic assessment of RC wall buildings

A procedure for incorporating record‐to‐record variability into the simplified seismic assessment of RC wall buildings is presented. The procedure relies on the use of the conditional spectrum to randomly sample spectral ordinates at relevant periods of vibration. For inelastic response, displacement reduction factors are then used to relate inelastic displacement demand to the spectral displacement at the effective period for single‐degree‐of‐freedom systems. Simple equations are used to convert back and forth between multi‐degree‐of‐freedom RC wall buildings and equivalent single‐degree‐of‐systems so that relevant engineering demand parameters can be obtained. Consideration is also given to higher‐mode effects by adapting existing modal combination rules. The proposed method is applied to several case study buildings, showing promising results in the examination of inter‐storey drift ratio and shear forces. The proposed method captures the variation in the distribution of structural response parameters that occurs with variations in structural configuration, intensity, engineering demand parameter of interest and site characteristics. Discussion is provided on possible ways to improve the accuracy of the procedure and suggestions for additional future work. Copyright © 2015 John Wiley & Sons, Ltd.     

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.