Incremental dynamic analysis for estimating seismic performance sensitivity and uncertainty

Incremental dynamic analysis (IDA) is presented as a powerful tool to evaluate the variability in the seismic demand and capacity of non‐deterministic structural models, building upon existing methodologies of Monte Carlo simulation and approximate moment‐estimation. A nine‐story steel moment‐resisting frame is used as a testbed, employing parameterized moment‐rotation relationships with non‐deterministic quadrilinear backbones for the beam plastic‐hinges. The uncertain properties of the backbones include the yield moment, the post‐yield hardening ratio, the end‐of‐hardening rotation, the slope of the descending branch, the residual moment capacity and the ultimate rotation reached. IDA is employed to accurately assess the seismic performance of the model for any combination of the parameters by performing multiple nonlinear timehistory analyses for a suite of ground motion records. Sensitivity analyses on both the IDA and the static pushover level reveal the yield moment and the two rotational‐ductility parameters to be the most influential for the frame behavior. To propagate the parametric uncertainty to the actual seismic performance we employ (a) Monte Carlo simulation with latin hypercube sampling, (b) point‐estimate and (c) first‐order second‐moment techniques, thus offering competing methods that represent different compromises between speed and accuracy. The final results provide firm ground for challenging current assumptions in seismic guidelines on using a median‐parameter model to estimate the median seismic performance and employing the well‐known square‐root‐sum‐of‐squares rule to combine aleatory randomness and epistemic uncertainty. Copyright © 2009 John Wiley & Sons, Ltd.     

Fast performance uncertainty estimation via pushover and approximate IDA

Approximate methods based on the static pushover are introduced to estimate the seismic performance uncertainty of structures having non‐deterministic modeling parameters. At their basis lies the use of static pushover analysis to approximate Incremental Dynamic Analysis (IDA) and estimate the demand and capacity epistemic uncertainty. As a testbed we use a nine‐storey steel frame having beam hinges with uncertain moment–rotation relationships. Their properties are fully described by six, randomly distributed, parameters. Using Monte Carlo simulation with Latin hypercube sampling, a characteristic ensemble of structures is created. The Static Pushover to IDA (SPO2IDA) software is used to approximate the IDA capacity curve from the appropriately post‐processed results of the static pushover. The approximate IDAs allow the evaluation of the seismic demand and capacity for the full range of limit‐states, even close to global dynamic instability. Moment‐estimating techniques such as Rosenblueth's point estimating method and the first‐order, second‐moment (FOSM) method are adopted as simple alternatives to obtain performance statistics with only a few simulations. The pushover is shown to be a tool that combined with SPO2IDA and moment‐estimating techniques can supply the uncertainty in the seismic performance of first‐mode‐dominated buildings for the full range of limit‐states, thus replacing semi‐empirical or code‐tabulated values (e.g. FEMA‐350), often adopted in performance‐based earthquake engineering. Copyright © 2009 John Wiley & Sons, Ltd.     

Analytical modelling of near‐source pulse‐like seismic demand for multi‐linear backbone oscillators

Nonlinear static procedures, which relate the seismic demand of a structure to that of an equivalent single‐degree‐of‐freedom oscillator, are well‐established tools in the performance‐based earthquake engineering paradigm. Initially, such procedures made recourse to inelastic spectra derived for simple elastic–plastic bilinear oscillators, but the request for demand estimates that delve deeper into the inelastic range, motivated investigating the seismic demand of oscillators with more complex backbone curves. Meanwhile, near‐source (NS) pulse‐like ground motions have been receiving increased attention, because they can induce a distinctive type of inelastic demand. Pulse‐like NS ground motions are usually the result of rupture directivity, where seismic waves generated at different points along the rupture front arrive at a site at the same time, leading to a double‐sided velocity pulse, which delivers most of the seismic energy. Recent research has led to a methodology for incorporating this NS effect in the implementation of nonlinear static procedures. Both of the previously mentioned lines of research motivate the present study on the ductility demands imposed by pulse‐like NS ground motions on oscillators that feature pinching hysteretic behaviour with trilinear backbone curves. Incremental dynamic analysis is used considering 130 pulse‐like‐identified ground motions. Median, 16% and 84% fractile incremental dynamic analysis curves are calculated and fitted by an analytical model. Least‐squares estimates are obtained for the model parameters, which importantly include pulse period T_p. The resulting equations effectively constitute an R ? μ ? T ? T_p relation for pulse‐like NS motions. Potential applications of this result towards estimation of NS seismic demand are also briefly discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

应用弹塑性反应谱对IDA方法的改进研究

针对传统IDA方法中结构的加速度反应远高于SPO分析结果的现象,采用弹塑性反应谱作为IDA的烈度度量,同时给出随地震烈度改变而不断变化的屈服强度系数,提出了改进的IDA方法,更有效地考虑了结构的非线性地震反应特征。实际工程分析表明,改进的IDA方法能够体现结构的非线性受力反应,同时兼顾结构的动力反应特征,与SPO曲线和由规范计算得到的结构承载力有更好的可比较性。     

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

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

Near‐optimal piecewise linear fits of static pushover capacity curves for equivalent SDOF analysis

The piecewise linear (‘multilinear’) approximation of realistic force‐deformation capacity curves is investigated for structural systems incorporating generalized plastic, hardening, and negative stiffness behaviors. This fitting process factually links capacity and demand and lies at the core of nonlinear static assessment procedures. Despite codification, the various fitting rules used can produce highly heterogeneous results for the same capacity curve, especially for the highly‐curved backbones resulting from the gradual plasticization or the progressive failures of structural elements. To achieve an improved fit, the error introduced by the approximation is quantified by studying it at the single‐degree‐of‐freedom level, thus avoiding any issues related to multi‐degree‐of‐freedom versus single‐degree‐of‐freedom realizations. Incremental dynamic analysis is employed to enable a direct comparison of the actual backbones versus their candidate piecewise linear approximations in terms of the spectral acceleration capacity for a continuum of limit‐states. In all cases, current code‐based procedures are found to be highly biased wherever widespread significant stiffness changes occur, generally leading to very conservative estimates of performance. The practical rules determined allow, instead, the definition of standardized low‐bias bilinear, trilinear, or quadrilinear approximations, regardless of the details of the capacity curve shape. Copyright © 2012 John Wiley & Sons, Ltd.     

A static predictor of seismic demand on frames based on a post‐elastic deflected shape

Predictors of seismic structural demands (such as inter‐storey drift angles) that are less time‐consuming than nonlinear dynamic analysis have proven useful for structural performance assessment and for design. Luco and Cornell previously proposed a simple predictor that extends the idea of modal superposition (of the first two modes) with the square‐root‐of‐sum‐of‐squares (SRSS) rule by taking a first‐mode inelastic spectral displacement into account. This predictor achieved a significant improvement over simply using the response of an elastic oscillator; however, it cannot capture well large displacements caused by local yielding. A possible improvement of Luco's predictor is discussed in this paper, where it is proposed to consider three enhancements: (i) a post‐elastic first‐mode shape approximated by the deflected shape from a nonlinear static pushover analysis (NSPA) at the step corresponding to the maximum drift of an equivalent inelastic single‐degree‐of‐freedom (SDOF) system, (ii) a trilinear backbone curve for the SDOF system, and (iii) the elastic third‐mode response for long‐period buildings. Numerical examples demonstrate that the proposed predictor is less biased and results in less dispersion than Luco's original predictor. Copyright © 2006 John Wiley & Sons, Ltd.     

Effect of the spectral shape of ground motion records on the collapse fragility assessment of degrading SDOF systems

The main purpose of this paper is to study the collapse capacity of single degree of freedom(SDOF) systems and to produce fragility curves as well as collapse capacity spectra while considering a broad range of structural parameters, including system degradation, the P-Δ effect, ductility capacity and the post-capping stiffness ratio. The modified IbarraKrawinkler deterioration model was used to consider hysteretic behavior. A comprehensive study was conducted to extract the collapse capacity spectrum of SDOF systems with a wide range of periods, varying from 0.05 to 4 s, to cover short, intermediate and long periods. Incremental dynamic analysis(IDA) was performed for SDOF systems to identify the condition in which the collapse capacity of the system is determined. The IDAs were performed using different sets of seismic ground motions. The ground motion records were categorized into different sets based on three spectral shape parameters, including the epsilon, Sa Ratio and N_p. The collapse fragility curves of SDOF systems with different periods were extracted to illustrate the collapse capacity at different probability levels. The results show that structural degradation and ductility as well as the spectral shape parameters significantly affect the collapse capacity of SDOF systems. On the other hand, the post-capping stiffness ratio and small levels of the P-Δ effect do not remarkably change collapse capacity. Also, the collapse capacity of SDOF systems is more sensitive to the records categorized based on Sa Ratio and N_p than those classified based on epsilon.     

Envelope‐based pushover analysis procedure for the approximate seismic response analysis of buildings

An envelope‐based pushover analysis procedure is presented that assumes that the seismic demand for each response parameter is controlled by a predominant system failure mode that may vary according to the ground motion. To be able to simulate the most important system failure modes, several pushover analyses need to be performed, as in a modal pushover analysis procedure, whereas the total seismic demand is determined by enveloping the results associated with each pushover analysis. The demand for the most common system failure mode resulting from the ‘first‐mode’ pushover analysis is obtained by response history analysis for the equivalent ‘modal‐based’ SDOF model, whereas demand for other failure modes is based on the ‘failure‐based’ SDOF models. This makes the envelope‐based pushover analysis procedure equivalent to the N2 method provided that it involves only ‘first‐mode’ pushover analysis and response history analysis of the corresponding ‘modal‐based’ SDOF model. It is shown that the accuracy of the approximate 16th, 50th and 84th percentile response expressed in terms of IDA curves does not decrease with the height of the building or with the intensity of ground motion. This is because the estimates of the roof displacement and the maximum storey drift due to individual ground motions were predicted with a sufficient degree of accuracy for almost all the ground motions from the analysed sets. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of deterioration modelling on the seismic response of steel moment frames designed to Eurocode 8

This paper assesses the influence of cyclic and in‐cycle degradation on seismic drift demands in moment‐resisting steel frames (MRF) designed to Eurocode 8. The structural characteristics, ground motion frequency content, and level of inelasticity are the primary parameters considered. A set of single‐degree‐of‐freedom (SDOF) systems, subjected to varying levels of inelastic demands, is initially investigated followed by an extensive study on multi‐storey frames. The latter comprises a large number of incremental dynamic analyses (IDA) on 12 frames modelled with or without consideration of degradation effects. A suite of 56 far‐field ground motion records, appropriately scaled to simulate 4 levels of inelastic demand, is employed for the IDA. Characteristic results from a detailed parametric investigation show that maximum response in terms of global and inter‐storey drifts is notably affected by degradation phenomena, in addition to the earthquake frequency content and the scaled inelastic demands. Consistently, both SDOF and frame systems with fundamental periods shorter than the mean period of ground motion can experience higher lateral strength demands and seismic drifts than those of non‐degrading counterparts in the same period range. Also, degrading multi‐storey frames can exhibit distinctly different plastic mechanisms with concentration of drifts at lower levels. Importantly, degrading systems might reach a “near‐collapse” limit state at ductility demand levels comparable to or lower than the assumed design behaviour factor, a result with direct consequences on optimised design situations where over‐strength would be minimal. Finally, the implications of the findings with respect to design‐level limit states are discussed.     

Seismic assessment of concrete gravity dams using capacity estimation and damage indexes

A new concept to determine state of the damage in concrete gravity dams is introduced. The Pine Flat concrete gravity dam has been selected for the purpose of the analysis and its structural capacity, assuming no sliding plane and rigid foundation, has been estimated using the two well‐known methods: nonlinear static pushover (SPO) and incremental dynamic analysis (IDA). With the use of these two methods, performance and various limit states of the dam have been determined, and three damage indexes have been proposed on the basis of the comparison of seismic demands and the dam's capacity. It is concluded that the SPO and IDA can be effectively used to develop indexes for seismic performance evaluation and damage assessment of concrete gravity dams. Copyright © 2012 John Wiley & Sons, Ltd.     

Prediction of the median IDA curve by employing a limited number of ground motion records

A methodology has been proposed which can be used to reduce the number of ground motion records needed for the reliable prediction of the median seismic response of structures by means of incremental dynamic analysis (IDA). This methodology is presently limited to predictions of the median IDA curve only. The reduction in the number of ground motion records needed to predict the median IDA curve is achieved by introducing a precedence list of ground motion records. The determination of such a list is an optimization problem, which is solved in the paper by means of (1) a genetic algorithm and (2) a proposed simple procedure. The seismic response of a simple, computationally non‐demanding structural model has been used as input data for the optimization problem. The presented example is a three‐storey‐reinforced concrete building, subjected to two sets of ground motion records, one a free‐field set and the other a near‐field set. It is shown that the median IDA curves can be predicted with acceptable accuracy by employing only four ground motion records instead of the 24 or 30 records, which are the total number of ground motion records for the free‐field and near‐field sets, respectively. Copyright © 2007 John Wiley & Sons, Ltd.     

Seismic response of self‐centring hysteretic SDOF systems

The seismic response of single‐degree‐of‐freedom (SDOF) systems incorporating flag‐shaped hysteretic structural behaviour, with self‐centring capability, is investigated numerically. For a SDOF system with a given initial period and strength level, the flag‐shaped hysteretic behaviour is fully defined by a post‐yielding stiffness parameter and an energy‐dissipation parameter. A comprehensive parametric study was conducted to determine the influence of these parameters on SDOF structural response, in terms of displacement ductility, absolute acceleration and absorbed energy. This parametric study was conducted using an ensemble of 20 historical earthquake records corresponding to ordinary ground motions having a probability of exceedence of 10% in 50 years, in California. The responses of the flag‐shaped hysteretic SDOF systems are compared against the responses of similar bilinear elasto‐plastic hysteretic SDOF systems. In this study the elasto‐plastic hysteretic SDOF systems are assigned parameters representative of steel moment resisting frames (MRFs) with post‐Northridge welded beam‐to‐column connections. In turn, the flag‐shaped hysteretic SDOF systems are representative of steel MRFs with newly proposed post‐tensioned energy‐dissipating connections. Building structures with initial periods ranging from 0.1 to 2.0s and having various strength levels are considered. It is shown that a flag‐shaped hysteretic SDOF system of equal or lesser strength can always be found to match or better the response of an elasto‐plastic hysteretic SDOF system in terms of displacement ductility and without incurring any residual drift from the seismic event. Copyright © 2002 John Wiley & Sons, Ltd.     

Force reduction factor for building structures equipped with added viscous dampers

This research work focuses on the analysis of the hysteretic seismic behaviour of inelastic SDOF systems equipped with viscous dampers. In detail, it is aimed at obtaining a practical tool useful for the seismic design of building structures with added dampers, within the framework of the traditional seismic design based on ductility. The objective is to evaluate the appropriate force reduction factor for highly damped (i.e. damping ratio greater than 5 %) SDOF systems able to guarantee a prescribed level of structural safety.     

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.     

Approximate incremental dynamic analysis using the modal pushover analysis procedure

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

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.     

Incremental dynamic analysis with consideration of modeling uncertainties

Incremental dynamic analysis (IDA) has been extended by introducing a set of structural models in addition to the set of ground motion records which is employed in IDA analysis in order to capture record‐to‐record variability. The set of structural models reflects epistemic (modeling) uncertainties, and is determined by utilizing the latin hypercube sampling (LHS) method. The effects of both aleatory and epistemic uncertainty on seismic response parameters are therefore considered in extended IDA analysis. The proposed method has been applied to an example of the four‐storey‐reinforced concrete frame, for which pseudo‐dynamic tests were performed at the ELSA Laboratory, Ispra. The influence of epistemic uncertainty on the seismic response parameters is presented in terms of summarized IDA curves and dispersion measures. The results of extended IDA analysis are compared with the results of IDA analysis, and the sensitivity of the seismic response parameters to the input random variable using the LHS method is discussed. It is shown that epistemic uncertainty does not have significant influence on the seismic response parameters in the range far from collapse, but could have a significant influence on collapse capacity. Copyright © 2008 John Wiley & Sons, Ltd.     

Non‐linear seismic analysis and vulnerability evaluation of a masonry building by means of the SAP2000 V.10 code

The aim of the paper is to explore the possibilities offered by SAP2000^® v.10, a software package with user‐friendly interface widely used by practising engineers, for seismic analyses of masonry buildings. The reliability of the code was first investigated by carrying out static push‐over (SPO) analyses of two walls, already analysed by other researchers using advanced programs. The equivalent frame modelling was employed in all analyses carried out. The code was then used to investigate the seismic performance of an existing two‐storey building typical of the north‐east of Italy, with the walls being made of roughly squared stones. An SPO analysis was performed first on the most significant wall, followed by a number of time‐history analyses aimed to evaluate the dynamic push‐over curves. Finally, the seismic fragility curves were derived, considering the seismic input as a random variable. Copyright © 2007 John Wiley & Sons, Ltd.     

The effects of repeated earthquake ground motions on the non‐linear response of SDOF systems

In many parts of the world, the repetition of medium–strong intensity earthquake ground motions at brief intervals of time has been observed. The new design philosophies for buildings in seismic areas are based on multi‐level design approaches, which take into account more than a single damageability limit state. According to these approaches, a sequence of seismic actions may produce important consequences on the structural safety. In this paper, the effects of repeated earthquake ground motions on the response of single‐degree‐of‐freedom systems (SDOF) with non‐linear behaviour are analysed. A comparison is performed with the effect of a single seismic event on the originally non‐damaged system for different hysteretic models in terms of pseudo‐acceleration response spectra, behaviour factor q and damage parameters. The elastic–perfect plastic system is the most vulnerable one under repeated earthquake ground motions and is characterized by a strong reduction of the q‐factor. A moment resisting steel frame is analysed as well, showing a reduction of the q‐factor under repeated earthquake ground motions even larger than that of an equivalent SDOF system. Copyright © 2002 John Wiley & Sons, Ltd.