Probabilistic estimation of maximum inelastic displacement demands for performance‐based design

A probabilistic approach to estimate maximum inelastic displacement demands of single‐degree‐of‐freedom (SDOF) systems is presented. By making use of the probability of exceedance of maximum inelastic displacement demands for given maximum elastic spectral displacement and the mean annual frequency of exceedance of elastic spectral ordinates, a simplified procedure is proposed to estimate mean annual frequencies of exceedance of maximum inelastic displacement demands. Simplifying assumptions are thoroughly examined and discussed. Using readily available elastic seismic hazard curves the procedure can be used to compute maximum inelastic displacement seismic hazard curves and uniform hazard spectra of maximum inelastic displacement demands. The resulting maximum inelastic displacement demand spectra provide a more rational way of establishing seismic demands for new and existing structures when performance‐based approaches are used. The proposed procedure is illustrated for elastoplastic SDOF systems having known‐lateral strength located in a region of high seismicity in California. Copyright © 2007 John Wiley & Sons, Ltd.     

Uniform fragility spectra for the performance‐based seismic design of structures considering variabilities in structural properties

This paper presents a new methodology based on structural performance to determine uniform fragility design spectra, i.e., spectra with the same probability of exceedance of a performance level for a given seismic intensity. The design spectra calculated with this methodology provide directly the lateral strength, in terms of yield‐ pseudo‐accelerations, associated with the rate of exceedance of a specific ductility characterizing the performance level for which the structures will be designed. This procedure involves the assessment of the seismic hazard using a large enough number of seismic records of several magnitudes; these records are simulated with an improved empirical Green function method. The statistics of the performance of a single degree of freedom system are obtained using Monte Carlo simulation considering the seismic demand, the fundamental period, and the strength of the structure as uncertain variables. With these results, the conditional probability that a structure exceeds a specific performance level is obtained. The authors consider that the proposed procedure is a significant improvement to others considered in the literature and a useful research tool for the further development of uniform fragility spectra that can be used for the performance‐based seismic design and retrofit of structures.     

A procedure on ground motion selection and scaling for nonlinear response of simple structural systems

This study presents a ground-motion selection and scaling methodology that preserves the basic seismological features of the scaled records with reduced scatter in the nonlinear structural response. The methodology modifies each strong-motion recording with known fundamental seismological parameters using the estimations of ground-motion prediction equations for a given target hazard level. It provides robust estimations on target building response through scaled ground motions and calculates the dispersion about this target. This alternative procedure is not only useful for record scaling and selection but, upon its further refinement, can also be advantageous for the probabilistic methods that assess the engineering demand parameters for a given target hazard level. Case studies that compare the performance of the proposed procedure with some other record selection and scaling methods suggest its usefulness for building performance assessment and loss models. Copyright © 2012 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.     

Assessing uncertainty in estimation of seismic response for PBEE

State‐of‐the‐art approaches to probabilistic assessment of seismic structural reliability are based on simulation of structural behavior via nonlinear dynamic analysis of computer models. Simulations are carried out considering samples of ground motions supposedly drawn from specific populations of signals virtually recorded at the site of interest. This serves to produce samples of structural response to evaluate the failure rate, which in turn allows to compute the failure risk (probability) in a time interval of interest. This procedure alone implies that uncertainty of estimation affects the probabilistic results. The latter is seldom quantified in risk analyses, although it may be relevant. This short paper discusses some basic issues and some simple statistical tools, which can aid the analyst towards the assessment of the impact of sample variability on fragility functions and the resulting seismic structural risk. On the statistical inference side, the addressed strategies are based on consolidated results such as the well‐known delta method and on some resampling plans belonging to the bootstrap family. On the structural side, they rely on assumptions and methods typical in performance‐based earthquake engineering applications. Copyright © 2017 John Wiley & Sons, Ltd.     

Modelling considerations in probabilistic performance‐based seismic evaluation: case study of the I‐880 viaduct

Limitations associated with deterministic methods to quantify demands and develop rational acceptance criteria have led to the emergence of probabilistic procedures in performance‐based seismic engineering. The Pacific Earthquake Engineering Research performance‐based methodology is one such approach. In this paper, the impact of certain modelling decisions made at different stages of the evaluation process on the performance assessment of a typical multi‐bent viaduct is examined. Modelling, in the context of this paper, covers hazard modelling, structural modelling and loss modelling. The specific application considered in this study is a section of an existing viaduct in California: the I‐880 interstate highway. Several simulation models of the viaduct are developed, a series of nonlinear time‐history analyses are carried out to predict demands, measures of damage are evaluated and the probability of closure of the viaduct is estimated using the specified hazard for the site. It is concluded that the methodology offers several advantages over existing deterministic performance‐based procedures. Results of the investigation indicate that the assessment methodology is particularly sensitive to the reliability of decisions made by bridge inspectors following a seismic event, and to the dispersion in the demand estimation, which in turn is influenced by several factors including soil–structure interaction effects and ground motion scaling procedures. Copyright © 2005 John Wiley & Sons, Ltd.     

Seismic performance of wood‐frame houses in south‐western British Columbia

The seismic performance of conventional wood‐frame structures in south‐western British Columbia is analytically investigated through incremental dynamic analysis by utilizing available UBC‐SAWS models, which were calibrated based on experimental test results. To define an adequate target response spectrum that is consistent with information from national seismic hazard maps, record selection/scaling based on the conditional mean spectrum (CMS) is implemented. Furthermore, to reflect complex seismic hazard contributions from different earthquake sources (i.e. crustal events, interface events, and inslab events), we construct CMS for three earthquake types, and use them to select and scale an adequate set of ground motion records for the seismic performance evaluation. We focus on the impacts of adopting different record selection criteria and of using different shear‐wall types (Houses 1–4; in terms of seismic resistance, House 1>House 2>House 3>House 4) on the nonlinear structural response. The results indicate that the record selection procedures have significant influence on the probabilistic relationship between spectral acceleration at the fundamental vibration period and maximum inter‐story drift ratio, highlighting the importance of taking into account response spectral shapes in selecting and scaling ground motion records. Subjected to ground motions corresponding to the return period of 2500 years, House 1 is expected to experience very limited extent of damage; Houses 2 and 3 may be disturbed by minor damage; whereas House 4 may suffer from major damage occasionally. Finally, we develop statistical models of the maximum inter‐story drift ratio conditioned on a seismic intensity level for wood‐frame houses, which is useful for seismic vulnerability assessment. Copyright © 2010 John Wiley & Sons, Ltd.     

Practice‐oriented estimation of the seismic demand hazard using ground motions at few intensity levels

This paper examines the calculation of the seismic demand hazard in a practice‐oriented manner via the use of seismic response analyses at few intensity levels. The seismic demand hazard is a more robust measure for quantifying seismic performance, when seismic hazard is represented in a probabilistic format, than intensity‐based assessments, which remain prevalent in seismic design codes. It is illustrated that, for a relatively complex bridge–foundation–soil system case study, the seismic demand hazard can be estimated with sufficient accuracy using as little as three intensity measure levels that have exceedance probabilities of 50%, 10% and 2% in 50 years which are already of interest in multi‐objective performance‐based design. Compared with the conventional use of the mean demand from an intensity‐based assessment(s), it is illustrated that, for the same number of seismic response analyses, a practice‐oriented ‘approximate’ seismic demand hazard is a more accurate and precise estimate of the ‘exact’ seismic demand hazard. Direct estimation of the seismic demand hazard also provides information of seismic performance at multiple exceedance rates. Thus, it is advocated that if seismic hazard is considered in a probabilistic format, then seismic performance assessment, and acceptance criteria, should be in terms of the seismic demand hazard and not intensity‐based assessments. Copyright © 2013 John Wiley & Sons, Ltd.     

Alternative non‐linear demand estimation methods for probability‐based seismic assessments

Alternative non‐linear dynamic analysis procedures, using real ground motion records, can be used to make probability‐based seismic assessments. These procedures can be used both to obtain parameter estimates for specific probabilistic assessment criteria such as demand and capacity factored design and also to make direct probabilistic performance assessments using numerical methods. Multiple‐stripe analysis is a non‐linear dynamic analysis method that can be used for performance‐based assessments for a wide range of ground motion intensities and multiple performance objectives from onset of damage through global collapse. Alternatively, the amount of analysis effort needed in the performance assessments can be reduced by performing the structural analyses and estimating the main parameters in the region of ground motion intensity levels of interest. In particular, single‐stripe and double‐stripe analysis can provide local probabilistic demand assessments using minimal number of structural analyses (around 20 to 40). As a case study, the displacement‐based seismic performance of an older reinforced concrete frame structure, which is known to have suffered shear failure in its columns during the 1994 Northridge Earthquake, is evaluated. Copyright © 2008 John Wiley & Sons, Ltd.     

A comparison of intensity‐based demand distributions and the seismic demand hazard for seismic performance assessment

This paper compares the seismic demands obtained from an intensity‐based assessment, as conventionally considered in seismic design guidelines, with the seismic demand hazard. Intensity‐based assessments utilize the distribution of seismic demand from ground motions that have a specific value of some conditioning intensity measure, and the mean of this distribution is conventionally used in design verification. The seismic demand hazard provides the rate of exceedance of various seismic demand values and is obtained by integrating the distribution of seismic demand at multiple intensity levels with the seismic hazard curve. The seismic demand hazard is a more robust metric for quantifying seismic performance, because seismic demands from an intensity‐based assessment: (i) are not unique, with different values obtained using different conditioning intensity measures; and (ii) do not consider the possibility that demand values could be exceeded from different intensity ground motions. Empirical results, for a bridge‐foundation‐soil system, illustrate that the mean seismic demand from an intensity‐based assessment almost always underestimates the demand hazard value for the exceedance rate considered, on average by 17% and with a large variability. Furthermore, modification factors based on approximate theory are found to be unreliable. Adopting the maximum of the mean values from multiple intensity‐based assessments, with different conditional intensity measures, provides a less biased prediction of the seismic demand hazard value, but with still a large variability, and a proportional increase the required number of analyses. For an equivalent number of analyses, direct computation of the seismic demand hazard is a more logical choice and provides additional performance insight. Copyright © 2013 John Wiley & Sons, Ltd.     

Near‐source seismic hazard and design scenarios

Earthquakes damage engineering structures near, relatively to the rupture's size, to the source. In this region, the fault's dynamics affect ground motion propagation differently from site to site, resulting in systematic spatial variability known as directivity. Although a number of researches recommend that records with directivity‐related velocity pulses should be explicitly taken into account when defining design seismic action on structures, probabilistic seismic hazard analysis (PSHA), in its standard version, seems inadequate for the scope. In the study, it is critically reviewed why, from the structural engineering point of view, hazard assessment should account for near‐source effects (i.e., pulse‐like ground motions), and how this can be carried out adjusting PSHA analytically via introduction of specific terms and empirically calibrated models. Disaggregation analysis and design scenarios for near‐source PSHA are also formulated. The analytical procedures are then applied to develop examples of hazard estimates for sites close to strike–slip or dip–slip faults and to address differences with respect to the ordinary case, that is, when pulse‐like effects are not explicitly accounted for. Significant increase of hazard for selected spectral ordinates is found in all investigated cases; increments depend on the fault‐site configuration. Moreover, to address design scenarios for seismic actions on structures, disaggregation results are also discussed, along with limitations of current design spectra to highlight the pulse‐like effects of structural interest. Finally, an attempt to overcome these, by means of disaggregation‐based scenarios specific for the pulse occurrence case, is presented. Copyright © 2012 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.     

Application of a Monte‐Carlo simulation approach for the probabilistic assessment of seismic hazard for geographically distributed portfolio

The conventional integral approach is very well established in probabilistic seismic hazard assessment (PSHA). However, Monte‐Carlo (MC) simulations can become an efficient and flexible alternative against conventional PSHA when more complicated factors (e.g. spatial correlation of ground shaking) are involved. This study aims at showing the implementation of MC simulation techniques for computing the annual exceedance rates of dynamic ground‐motion intensity measures (GMIMs) (e.g. peak ground acceleration and spectral acceleration). We use multi‐scale random field technique to incorporate spatial correlation and near‐fault directivity while generating MC simulations to assess the probabilistic seismic hazard of dynamic GMIMs. Our approach is capable of producing conditional hazard curves as well. We show various examples to illustrate the potential use of the proposed procedures in the hazard and risk assessment of geographically distributed structural systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Accounting for directionality as a function of structural typology in performance‐based earthquake engineering design

For the seismic design of a structure, horizontal ground shaking is usually considered in two perpendicular directions, even though real horizontal ground motions are complex two‐dimensional phenomena that impose different demands at different orientations. While the issue of ground motion dependence on the orientation of the recording devices has been the focus of many significant developments during the last decade, the effects of directionality on the characteristics of the structure have received less attention. This work presents a proposal to calculate the probability of exceedance of elastic spectral displacements accounting for structural typology and illustrates its relevance by means of its application to two case‐study buildings. In order to ease its implementation in seismic design codes, a simplification is developed by means of a detailed statistical analysis of the results obtained using four sets of real hazard curves. The framework presented herein is considered to represent an important contribution to the field of performance‐based earthquake engineering, permitting improved treatment of directionality effects within seismic risk design and assessment. Copyright © 2016 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.     

A probabilistic framework for correlated seismic downtime and repair cost estimation of geo‐structures

An important component of probabilistic risk assessment methods is the development of models to quantify the direct consequences of damage to geo‐structural components for a given intensity of the hazard. This paper presents a general probabilistic framework for correlated repair cost and downtime estimation of geo‐structures exposed to seismic hazards. The framework uses as input the results of nonlinear time‐history analysis of geo‐structures for the set of earthquake records that are representative of the seismic hazard models for the region of interest. The repair cost and downtime are estimated for individual earthquakes probabilistically considering the uncertainties associated with damage states. In addition, the formulation of the repair cost and downtime accounts for the reduction in the repair requirements as the number of damaged components in the given damage state increases. An analytical linear and two bilinear regression models are proposed for conditional correlated seismic repair cost and downtime estimation of geo‐structures given the intensity measure. The proposed framework is demonstrated by developing seismic repair models of a typical pile‐supported wharf structure on the west coast of the United States. The presented framework is general and can be applied to other types of geo‐structures and hazards and can include other decision variables such as loss of life as well. Copyright © 2013 John Wiley & Sons, Ltd.     

Probabilistic seismic hazard assessment in SE-Spain based on macroseismic site histories

In many countries such as Spain earthquake databases still mainly comprise macroseismic data from felt effects. The full exploit of this information is of basic importance for seismic risk assessment and emergency planning, given the strict link between macroseismic intensity and damage. A probabilistic procedure specifically developed to handle macroseismic data, mostly relying on site information and seismogenic-source free, has been applied to evaluate seismic hazard in SE-Spain (Alicante-Murcia region). Present seismicity is moderate-low with largest magnitudes slightly over Mw5.0. The historical record includes very destructive earthquakes, maximum EMS98 intensities reaching IX–X and X in the nineteenth century (e.g., Torrevieja 1829 earthquake). Very recently, two events in the area on 11 May 2011 (Mw4.5, Mw5.2) killed nine people, injured 300, and produced important damage in the city of Lorca. Regional hazard maps for the area together with specific hazard curves at selected localities are obtained. Results are compared with the maximum observed intensities in the period 1300–2012, and with the values in the seismic hazard map from the Spanish Building Code in force. In general, the maximum felt intensity values are closer to the hazard values calculated for 2 % probability of exceedance in 50 years, using felt and expected intensity. The intensity-based probabilistic hazard maps obtained through the applied approach reduce the inherent smoothing of those based on standard probabilistic seismic hazard assessment approaches for the region, allowing identifying possible over- or sub-estimates of site hazard values, providing very valuable information for risk reduction strategies or for future updates of the building code hazard maps.     

Error estimation of closed‐form solution for annual rate of structural collapse

With the increasing emphasis of performance‐based earthquake engineering in the engineering community, several investigations have been presented outlining simplified approaches suitable for performance‐based seismic design (PBSD). Central to most of these PBSD approaches is the use of closed‐form analytical solutions to the probabilistic integral equations representing the rate of exceedance of key performance measures. Situations where such closed‐form solutions are not appropriate primarily relate to the problem of extrapolation outside of the region in which parameters of the closed‐form solution are fit. This study presents a critical review of the closed‐form solution for the annual rate of structural collapse. The closed‐form solution requires the assumptions of lognormality of the collapse fragility and power model form of the ground motion hazard, of which the latter is more significant regarding the error of the closed‐form solution. Via a parametric study, the key variables contributing to the error between the closed‐form solution and solution via numerical integration are illustrated. As these key variables cannot be easily measured, it casts doubt on the use of such closed‐form solutions in future PBSD, especially considering the simple and efficient nature of using direct numerical integration to obtain the solution. Copyright © 2008 John Wiley & Sons, Ltd.     

Tall building performance‐based seismic design using SCEC broadband platform site‐specific ground motion simulations

The scarcity of strong ground motion records presents a challenge for making reliable performance assessments of tall buildings whose seismic design is controlled by large‐magnitude and close‐distance earthquakes. This challenge can be addressed using broadband ground‐motion simulation methods to generate records with site‐specific characteristics of large‐magnitude events. In this paper, simulated site‐specific earthquake seismograms, developed through a related project that was organized through the Southern California Earthquake Center (SCEC) Ground Motion Simulation Validation (GMSV) Technical Activity Group, are used for nonlinear response history analyses of two archetype tall buildings for sites in San Francisco, Los Angeles, and San Bernardino. The SCEC GMSV team created the seismograms using the Broadband Platform (BBP) simulations for five site‐specific earthquake scenarios. The two buildings are evaluated using nonlinear dynamic analyses under comparable record suites selected from the simulated BBP catalog and recorded motions from the NGA‐West database. The collapse risks and structural response demands (maximum story drift ratio, peak floor acceleration, and maximum story shear) under the BBP and NGA suites are compared. In general, this study finds that use of the BBP simulations resolves concerns about estimation biases in structural response analysis which are caused by ground motion scaling, unrealistic spectral shapes, and overconservative spectral variations. While there are remaining concerns that strong coherence in some kinematic fault rupture models may lead to an overestimation of velocity pulse effects in the BBP simulations, the simulations are shown to generally yield realistic pulse‐like features of near‐fault ground motion records.     

Site dependence and record selection schemes for building fragility and regional loss assessment

When performing loss assessment of a geographically dispersed building portfolio, the response or loss (fragility or vulnerability) function of any given archetype building is typically considered to be a consistent property of the building itself. On the other hand, recent advances in record selection have shown that the seismic response of a structure is, in general, dependent on the nature of the hazard at the site of interest. This apparent contradiction begs the question: Are building fragility and vulnerability functions independent of site, and if not, what can be done to avoid having to reassess them for each site of interest? In the following, we show that there is a non‐negligible influence of the site, the degree of which depends on the intensity measure adopted for assessment. Employing a single‐period (e.g., first‐mode), spectral acceleration would require careful record selection at each site and result to significant site‐to‐site variability of the fragility or vulnerability function. On the other hand, an intensity measure comprising the geometric mean of multiple spectral accelerations considerably reduces such variability. In tandem with a conditional spectrum record selection that accounts for multiple sites, it can offer a viable approach for incorporating the effect of site dependence into fragility and vulnerability estimates. Copyright © 2017 John Wiley & Sons, Ltd.