Evaluation of modal pushover analysis using generic frames

An Erratum has been published for this article in Earthquake Engineering and Structural Dynamics 2003; 32:1795. The recently developed modal pushover analysis (MPA) has been shown to be a significant improvement over the pushover analysis procedures currently used in structural engineering practice. None of the current invariant force distributions accounts for the contribution of higher modes—higher than the fundamental mode—to the response or for redistribution of inertial forces because of structural yielding. By including the contributions of a sufficient number of modes of vibration (generally two to three), the height‐wise distribution of responses estimated by MPA is generally similar to the ‘exact’ results from non‐linear response history analysis (RHA). Although the results of the previous research were extremely promising, only a few buildings were evaluated. The results presented below evaluate the accuracy of MPA for a wide range of buildings and ground motion ensembles. The selected structures are idealized frames of six different heights: 3, 6, 9, 12, 15, and 18 stories and five strength levels corresponding to SDF‐system ductility factor of 1, 1.5, 2, 4, and 6; each frame is analysed for 20 ground motions. Comparing the median values of storey‐drift demands determined by MPA to those obtained from non‐linear RHA shows that the MPA predicts reasonably well the changing height‐wise variation of demand with building height and SDF‐system ductility factor. Median and dispersion values of the ratios of storey‐drift demands determined by MPA and non‐linear‐RHA procedures were computed to measure the bias and dispersion of MPA estimates with the following results: (1) the bias and dispersion in the MPA procedure tend to increase for longer‐period frames and larger SDF‐system ductility factors (although these trends are not perfect); (2) the bias and dispersion in MPA estimates of seismic demands for inelastic frames are usually larger than for elastic systems; (3) the well‐known response spectrum analysis (RSA), which is equivalent to the MPA for elastic systems, consistently underestimates the response of elastic structures, e.g. up to 18% in the upper‐storey drifts of 18‐storey frames. Finally, the MPA procedure is simplified to facilitate its implementation in engineering practice—where the earthquake hazard is usually defined in terms of a median (or some other percentile) design spectrum for elastic systems—and the accuracy of this simplified procedure is documented. Copyright © 2003 John Wiley & Sons, Ltd.     

Equating incremental dynamic analysis with static nonlinear analysis at near-field excitation

In an incremental dynamic analysis(IDA) using a set of ground motion records,nonlinear time history analysis needs to be performed on structures.It is well recognized that IDA calls for high computational efforts and the results are highly sensitive to selected ground motions.As a result,alternative static methods are needed.This study aims to introduce a new double-stage(N1- N2) static method to estimate capacity curves of MR frames.The technique is regulated to resemble IDA results with specific emphasis on near-field ground motions.Using an ensemble of 56 near-field earthquake records,required ID As have been carried out for SAC-Los Angeles 3-,9- and 20-story buildings and an additional 15-story building.The results of the proposed static method are compared with those from IDA,displacement-based adaptive procedure(DAP),and multimodal procedure(MMP).The results indicate that in addition to enhanced accuracy,very little time is required in the case of N1-N2 method.Thus,for the 3-story structure,the time required is less than 1 minute.The proposed N1-N2 method shows the best accuracy in terms of lateral mechanisms for the 15-story frame while for the other cases,the first mode load pattern leads to the best accuracy.     

Application of MPA to estimate probability of collapse of structures

This paper develops a modal pushover analysis‐ (MPA) based approximate procedure to quantify the collapse potential of structural systems. The computationally demanding incremental dynamic analysis (IDA) of the structural system is avoided by MPA of the structure in conjunction with empirical equations for the collapse strength ratio for the first‐mode single‐degree‐of‐freedom (SDF) system; higher modes of vibration play essentially no role in estimating the ground motion intensity required to cause collapse of the structure. Presented are collapse fragility curves for 6‐, 9‐, and 20‐story regular special moment‐resisting teel frames computed by the exact and approximate procedures, demonstrating that the MPA‐based approximate procedure requires only a small fraction (1% in one example) of the computational effort inherent in exact IDA and still achieves highly accurate results. Copyright © 2010 John Wiley & Sons, Ltd.     

The extended N2 method taking into account higher mode effects in elevation

The N2 method has been extended in order to take into account higher mode effects in elevation. The extension is based on the assumption that the structure remains in the elastic range when vibrating in higher modes. The seismic demand in terms of displacements and storey drifts can be obtained by enveloping the results of basic pushover analysis and the results of standard elastic modal analysis. The approach is consistent with the extended N2 method used for plan‐asymmetric buildings. The proposed procedure was applied to three variants of three steel frame buildings used in the SAC project. The structural response was investigated for two sets of ground motions. Different ground motion intensities were used in order to investigate the influence of the magnitude of plastic deformations. The N2 results were compared with the results of nonlinear response‐history analysis, two other pushover‐based methods (modal pushover analysis (MPA) and modified MPA (MMPA)), and pushover analysis without consideration of higher modes. It was found that a considerable influence of higher modes on storey drifts is present at the upper part of medium‐and high‐rise structures. This effect is the largest in the case of elastic behaviour and decreases with ground motion intensity. The higher mode effects also depend on the spectral shape. The approximate methods (extended N2, MPA and MMPA) are able to provide fair estimates of response in the case of the test examples. Accuracy decreases with the height of the building, and with the intensity of ground motion. The N2 results are generally conservative. Copyright © 2011 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.     

Evaluation of the influence of vertical irregularities on the seismic performance of a nine‐storey steel frame

A methodology based on incremental dynamic analysis (IDA) is presented for the evaluation of structures with vertical irregularities. Four types of storey‐irregularities are considered: stiffness, strength, combined stiffness and strength, and mass irregularities. Using the well‐known nine‐storey LA9 steel frame as a base, the objective is to quantify the effect of irregularities, both for individual and for combinations of stories, on its response. In this context a rational methodology for comparing the seismic performance of different structural configurations is proposed by means of IDA. This entails performing non‐linear time history analyses for a suite of ground motion records scaled to several intensity levels and suitably interpolating the results to calculate capacities for a number of limit‐states, from elasticity to final global instability. By expressing all limit‐state capacities with a common intensity measure, the reference and each modified structure can be naturally compared without needing to have the same period or yield base shear. Using the bootstrap method to construct appropriate confidence intervals, it becomes possible to isolate the effect of irregularities from the record‐to‐record variability. Thus, the proposed methodology enables a full‐range performance evaluation using a highly accurate analysis method that pinpoints the effect of any source of irregularity for each limit‐state. Copyright © 2006 John Wiley & Sons, Ltd.     

Evaluation of three‐dimensional modal pushover analysis for unsymmetric‐plan buildings subjected to two components of ground motion

The accuracy of the three‐dimensional modal pushover analysis (MPA) procedure in estimating seismic demands for unsymmetric‐plan buildings due to two horizontal components of ground motion, simultaneously, is evaluated. Eight low‐and medium‐rise structures were considered. Four intended to represent older buildings were designed according to the 1985 Uniform Building Code, whereas four other designs intended to represent newer buildings were based on the 2006 International Building Code. The median seismic demands for these buildings to 39 two‐component ground motions, scaled to two intensity levels, were computed by MPA and nonlinear response history analysis (RHA), and then compared. Even for these ground motions that deform the buildings significantly into the inelastic range, MPA offers sufficient degree of accuracy. It is demonstrated that PMPA, a variant of the MPA procedure, for nonlinear systems is almost as accurate as the well‐known standard response spectrum analysis procedure is for linear systems. Thus, for practical applications, the PMPA procedure offers an attractive alternative to nonlinear RHA, whereby seismic demands can be estimated directly from the (elastic) design spectrum. In contrast, the nonlinear static procedure specified in the ASCE/SEI 41‐06 Standard is demonstrated to grossly underestimate seismic demands for some of the unsymmetric‐plan buildings considered. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Three‐dimensional modal pushover analysis of buildings subjected to two components of ground motion,including its evaluation for tall buildings

The modal pushover analysis (MPA) procedure, presently restricted to one horizontal component of ground motion, is extended to three‐dimensional analysis of buildings—symmetric or unsymmetric in plan—subjected to two horizontal components of ground motion, simultaneously. Also presented is a variant of this method, called the practical modal pushover analysis (PMPA) procedure, which estimates seismic demands directly from the earthquake response (or design) spectrum. Its accuracy in estimating seismic demands for very tall buildings is evaluated, demonstrating that for nonlinear systems this procedure is almost as accurate as the response spectrum analysis procedure is for linear systems. Thus, for practical applications, the PMPA procedure offers an attractive alternative whereby seismic demands can be estimated directly from the (elastic) design spectrum, thus avoiding the complications of selecting and scaling ground motions for nonlinear response history analysis. Copyright © 2010 John Wiley & Sons, Ltd.     

Estimating the seismic displacement of friction pendulum isolators based on non‐linear response history analysis

A procedure for developing equations that estimate the isolator displacement due to strong ground motion is applied to buildings isolated with the friction pendulum system. The resulting design equations, based on rigorous non‐linear analysis, offer an alternative to the iterative equivalent‐linear methods used by current U.S. building codes. The governing equations of the system are reduced to a form such that the median normalized displacement of the system due to an ensemble of ground motions is found to depend on only the isolation period—a function of the curvature of the isolator—and the friction force at incipient slip normalized by peak ground velocity. The normalization is effective in minimizing the dispersion of the normalized displacement for an ensemble of ground motions, implying that the median normalized displacement is a reliable estimate of response. The design equations reflect the significant (20 to 38%) increase in displacement when the excitation includes two lateral components of ground motion instead of just one component. Equivalent‐linear methods are shown to underestimate by up to 30% the exact median displacement determined by non‐linear response history analysis for one component of ground motion, and building codes include at most a 4.4% increase for a second component. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic response analysis using characteristic ground motion records for risk‐based decision‐making (3R method)

The concept of intensity‐based assessment for risk‐based decision‐making is introduced. It is realized by means of the so‐called 3R method (response analysis, record selection and risk‐based decision‐making), which can be used to check the adequacy of design of a new building or of the strengthening of an existing building by performing conventional pushover analysis and dynamic analysis for only a few ground motions, which are termed characteristic ground motions. Because the objective of the method is not a precise assessment of the seismic risk, a simple decision model for risk acceptability can be introduced. The engineer can decide that the reliability of a no‐collapse requirement is sufficient when collapse is observed in the case of less than half of, for example, seven characteristic ground motions. From the theoretical point of view, it is shown that the accuracy of the method is acceptable if the non‐linear response history analyses are performed at a low percentile of limit‐state intensity, which is also proven by means of several examples of multi‐storey reinforced concrete frame buildings. The 3R method represents a compromise between the exclusive use of either pushover analysis or dynamic analysis and can be easily introduced into building codes provided that its applicability is further investigated (e.g. asymmetric structures and other performance objectives) and that the procedure for the selection of characteristic ground motions is automated and readily available to engineers (www.smartengineering.si).     

Optimal design of supplemental viscous dampers for linear framed structures

A methodology for the optimal design of supplemental viscous dampers for framed structures is presented. It addresses the problem of minimizing the added damping subject to a constraint on the maximal interstorey angular drift for an ensemble of realistic ground motion records while assuming linear behaviour of the damped structure. The solution is achieved by actually solving an equivalent optimization problem of minimizing the added damping subject to a constraint on a maximal weighted integral on the squared angular drift. The computational effort is appreciably reduced by first using one ‘active’ ground motion record. If the resulting optimal design fails to satisfy the constraints for other ground motions from the original ensemble, additional ground motions (loading conditions) are added one by one to the ‘active’ set until the optimum is reached. An efficient selecting process which is presented herein will usually require one or two records to attain an optimum design. Examples of optimal designs of supplemental dampers are presented for a 2‐storey shear frame and a 10‐storey industrial frame. The 2‐storey shear frame is required to withstand one given ground motion whereas the 10‐storey frame is required to withstand an ensemble of twenty ground motions. The resulting viscously damped structures have envelope values of interstorey drifts equal or less than the target drifts. Copyright © 2005 John Wiley & Sons, Ltd.     

A macro‐element model for inelastic building analysis

A three‐dimensional model for approximate inelastic analysis of buildings is presented herein. The model is based on a single macro‐element per building storey. The inelastic properties of the model are characterized by the so‐called ultimate storey shear and torque (USST) surfaces. Different algorithms for the construction of these surfaces, as well as their applications in building modelling, are presented and discussed. Two alternative procedures are developed to integrate the force‐deformation constitutive relationship of the macro‐elements. The first one follows the exact trajectory of the load path of the structure on the USST, and the second uses linear programming without ever forming the USST surface. The accuracy of the model and integration procedure is evaluated by means of the earthquake response of single‐storey systems. The model and integration procedure developed is finally used to compute the inelastic response of a seven‐storey R/C building. The results of this investigation show that the model proposed, although approximate, can be effective in estimating the inelastic deformation demand of a building. It also enables the engineer to capture and interpret important features of the three‐dimensional inelastic response of a structure even before performing any inelastic dynamic analysis. Copyright © 2000 John Wiley & Sons, Ltd.     

Simplified non‐linear seismic analysis of infilled reinforced concrete frames

The N2 method for simplified non‐linear seismic analysis has been extended in order to make it applicable to infilled reinforced concrete frames. Compared to the simple basic variant of the N2 method, two important differences apply. A multi‐linear idealization of the pushover curve, which takes into account the strength degradation which occurs after the infill fails, has to be made, and specific reduction factors, developed in a companion paper, have to be used for the determination of inelastic spectra. It is shown that the N2 method can also be used for the determination of approximate summarized IDA curves. The proposed method was applied to two test buildings. The results were compared with the results obtained by non‐linear dynamic analyses for three sets of ground motions, and a reasonable accuracy was demonstrated. A similar extension of the N2 method can be made to any structural system, provided that an appropriate specific R–µ–T relation is available. Copyright © 2004 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.     

Evaluation of residual drift demands in regular multi‐storey frames for performance‐based seismic assessment

This paper summarizes results of a comprehensive analytical study aimed at evaluating the amplitude and heightwise distribution of residual drift demands in multi‐storey moment‐resisting frames after earthquake excitation. For that purpose, a family of 12 one‐bay two‐dimensional generic frame models was subjected to an ensemble of 40 ground motions scaled to different intensities. In this investigation, an inelastic ground motion intensity measure was employed to scale each record, which allowed reducing the record‐to‐record variability in the estimation of residual drift demands. The results were statistically processed in order to evaluate the influence of ground motion intensity, number of stories, period of vibration, frame mechanism, system overstrength, and hysteretic behaviour on central tendency of residual drift demands. In addition, a special emphasis was given to evaluate the uncertainty in the estimation of residual drift demands. Results of incremental dynamic analyses indicate that the amplitude and heightwise distribution of residual drift demands strongly depends on the frame mechanism, the heightwise system structural overstrength and the component hysteretic behaviour. An important conclusion for performance‐based assessment is that the evaluation of residual drift demands involves significantly larger levels of uncertainty (i.e. record‐to‐record variability) than that of maximum drift demands, which suggests that this variability and corresponding uncertainty should be explicitly taken into account when estimating residual drift demands during performance‐based seismic assessment of frame buildings. Copyright © 2006 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.     

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.     

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.     

A modal pushover analysis procedure to estimate seismic demands for unsymmetric‐plan buildings

An Erratum has been published for this article in Earthquake Engng. Struct. Dyn. 2004; 33:1429. Based on structural dynamics theory, the modal pushover analysis (MPA) procedure retains the conceptual simplicity of current procedures with invariant force distribution, now common in structural engineering practice. The MPA procedure for estimating seismic demands is extended to unsymmetric‐plan buildings. In the MPA procedure, the seismic demand due to individual terms in the modal expansion of the effective earthquake forces is determined by non‐linear static analysis using the inertia force distribution for each mode, which for unsymmetric buildings includes two lateral forces and torque at each floor level. These ‘modal’ demands due to the first few terms of the modal expansion are then combined by the CQC rule to obtain an estimate of the total seismic demand for inelastic systems. When applied to elastic systems, the MPA procedure is equivalent to standard response spectrum analysis (RSA). The MPA estimates of seismic demand for torsionally‐stiff and torsionally‐flexible unsymmetric systems are shown to be similarly accurate as they are for the symmetric building; however, the results deteriorate for a torsionally‐similarly‐stiff unsymmetric‐plan system and the ground motion considered because (a) elastic modes are strongly coupled, and (b) roof displacement is underestimated by the CQC modal combination rule (which would also limit accuracy of RSA for linearly elastic systems). Copyright © 2004 John Wiley & Sons, Ltd.