The comparative performance of seismic response spectrum combination rules in building analysis

The peak dynamic responses of two mathematical models of a fifteen-storey steel moment resisting frame building subjected to three earthquake excitations are computed by the response spectrum and time history methods. The models examined are: a ‘regular’ building in which the centres of stiffness and mass are coincident resulting in uncoupled modes with well-separated periods in each component direction of response; and an ‘irregular’ building with the mass offset from the stiffness centre of the building causing coupled modes with the translational modes having closely spaced periods. Four response spectrum modal combination rules are discussed and are used to predict the peak responses: (1) the square root of the sum of the squares (SRSS) method; (2) the double sum combination (DSC) method; (3) the complete quadratic combination (CQC) method; and (4) the absolute sum (ABS) method. The response spectrum results are compared to the corresponding peak time history values to evaluate the accuracy of the different combination rules. The DSC and the CQC methods provide good peak response estimates for both the regular and irregular building models. The SRSS method provides good peak response estimates for the regular building, but yields significant errors in the irregular building response estimates. The poor accuracy in the irregular building results is attributable to the effects of coupled modes with closely spaced periods. It is concluded that the DSC and CQC methods produce response estimates of equivalent accuracy. Both methods are recommended for general use. In addition to the DSC and CQC rules, the SRSS method is recommended for systems where coupled modes with closely spaced periods do not dominate the response.     

Envelopes of maximum seismic response for a partially symmetric single storey building model

An analysis is made of the coupled lateral-torsional response of a partially symmetric single-storey building model to horizontal translatory earthquake excitation. Interest centres on the evaluation of realistic estimates for two equivalent static actions (a shear and a torque) which account for the worst dynamic consequences of torsional unbalance. The results substantiate the findings of previous investigations which have given rise to the belief that strong modal coupling and severely coupled lateral and torsional responses are possible even in nominally symmetric buildings. The response of the model is assumed to be linearly elastic and viscously damped. In a preliminary analysis the equations of motion are solved using the modal analysis technique and the conditions necessary for full modal coupling are ascertained. Then by employing the design spectrum concept, together with suitably conservative procedures for combining the modal maxima, dimensionless forms of the equivalent static actions are evaluated as functions of two independent parameters. The final results are furnished by modified square root of the sum of the squares (SRSS) combination functions which take account of the spacing between the translational and torsional frequencies. Examples at the end of the paper illustrate the practical significance of the work.     

Simplified transverse seismic analysis of buried structures

This paper presents a simplified method for the analysis of square cross-section buried structures (tunnel) subjected to seismic motion. Finite element analyses are performed to assess the fundamental modes of vibration of the soil layer with and without the tunnel. The influence of the tunnel on the modes of vibrations is taken into account by comparing the modal deformations in the free-field to those in the presence of the tunnel. From this comparison the zone of influence of the modal displacements due to the presence of the structure is determined. The resulting model is subjected to horizontal and vertical excitation of statistically independent accelerograms compatible with the response spectra of the Regulatory Guide 1.6 of the Nuclear Energy Commission. The free-field displacement is introduced at the boundaries of the zone of influence. The proposed simplified static analysis yields a state of stresses similar to that obtained from a full dynamic analysis of the complete soil–tunnel system. Several examples are solved to corroborate the validity of the method.     

Seismic design response by an alternative SRSS rule

The square root of the sum of the squares (SRSS) procedure and its modified forms are often used to obtain seismic design response. The design inputs for such procedures are usually defined in terms of pseudo velocity or acceleration response spectra. Erroneous results have been obtained with these existing SRSS procedures, especially in the calculation of responses where high frequency effects dominate. Here an alternative SRSS procedure is developed using the so-called mode acceleration approach of structural dynamics. The design input in this procedure is defined in terms of relative acceleration and relative velocity spectra. The relative spectra can be related to pseudo spectra. For a given number of modes to be included in the analysis the new SRSS rule proposed here will reduce the error due to the so-called ‘missing mass’ effect and predict a more accurate response value than the rules which use pseudo spectra as input, for systems either with or without dominant high frequency mode effects.     

Modal decomposition method for stationary response of non-classically damped systems

The stationary response of multi-degree-of-freedom non-classically damped linear systems subjected to stationary input excitation is studied. A modal decomposition procedure based on the complex eigenvectors and eigenvalues of the system is used to derive general expressions for the spectral moments of response. These expressions are in terms of cross-modal spectral moments and explicitly account for the correlation between modal responses; thus, they are applicable to structures characterized with significant non-classical damping as well as structures with closely spaced frequencies. Closed form solutions are presented for the important case of response to white-noise input. Various quantities of response of general engineering interest can be obtained in terms of these spectral moments. These include mean zero-crossing rate and mean, variance and distribution of peak response over a specified duration. Examples point out several instances where non-classical damping effects become significant and illustrate the marked improvement of the results of this study over conventional analysis based on classical damping approximations.     

Mode‐acceleration approach to seismic response of multiply‐supported secondary systems

A mode‐acceleration approach has been proposed for estimating the seismic response of a linear, classically‐damped, multiply‐supported secondary system within the framework of a power spectral density function (PSDF)‐based stochastic approach, while the primary system is linear and classically‐damped. Response transfer functions have been formulated in terms of chosen numbers of fixed‐base modes of the primary and secondary systems. The proposed approach does not involve the determination of combined system properties, and is applicable to the secondary systems with high mass ratios also. Through a few example primary–secondary systems and an example band‐limited white noise excitation, it has been shown that this approach leads to reasonably accurate results when only a few primary and secondary modes are to be considered. The proposed formulation has been used to obtain input data for a decoupled response spectrum analysis of secondary systems. This data accurately accounts for the effects of interaction between the primary and secondary systems. It is shown to lead to substantial reductions in the errors associated with the envelope spectrum method in the case of moderately heavy to heavy secondary systems and when the spatial coupling does not play a major role. Copyright © 2002 John Wiley & Sons, Ltd.     

Seismic response of torsionally coupled base isolated structures

An analytical study of the seismic response of typical base isolated structures mounted on rubber bearings is presented. Isolated buildings are liable to have closely spaced lower modes of vibration with small eccentricity between centres of mass and rigidity. The isolated structure is modelled as a rigid deck with lumped masses supported on axially inextensible elastomeric rubber bearings. This simplified system has three degrees of freedom (dof), two translations and one rotation in the horizontal plane. The Green's functions for the displacement response of the 3 dof system are derived for both undamped and damped cases with small and large eccentricities. The small eccentricity case is taken from a specific isolated building, while the large eccentricity case arises from the 5 per cent accidental eccentricity which is required by various seismic codes. An interaction equation for normalized displacements is established for an idealized flat velocity spectrum or hyperbolic acceleration spectrum. An isolated building on rubber bearings would have its fundamental period fall into this range of a design spectrum. Numerical results for the specific building subjected to the El Centro earthquake of 1940 are presented. Both the time history and the response spectrum modal superposition analysis were performed. In the response spectrum analysis, the Complete Quadratic Combination (CQC) showed superiority over the Square Root of the Sum of Squares (SRSS) in estimating maximum responses. It is concluded that the effect of torsional coupling on the transient response of base isolated structures is insignificant, due to the combined effect of the time lag between the maximum translational and torsional responses and the influence of damping in the isolation system which for elastomeric bearings can be as high as 8 to 10 per cent.     

A response spectrum method for random vibration analysis of mdf systems

A response spectrum method for stationary random vibration analysis of linear, multi-degree-of-freedom systems is developed. The method is based on the assumption that the input excitation is a wide-band, stationary Gaussian process and the response is stationary. However, it can also be used as a good approximation for the response to a transient stationary Gaussian input with a duration several times longer than the fundamental period of the system. Various response quantities, including the mean-squares of the response and its time derivative, the response mean frequency, and the cumulative distribution and the mean and variance of the peak response are obtained in terms of the ordinates of the mean response spectrum of the input excitation and the modal properties of the system. The formulation includes the cross-correlation between modal responses, which is shown to be significant for modes with closely spaced natural frequencies. The proposed procedure is demonstrated for an example structure that is subjected to an ensemble of earthquake-induced base excitations. Computed results based on the response spectrum method are in close agreement with simulation results obtained from time-history dynamic analysis. The significance of closely spaced modes and the error associated with a conventional method that neglects the modal correlations are also demonstrated.     

Response of tall cantilever wall buildings to strong pulse type seismic excitation

This paper investigates the seismic response of tall cantilever wall buildings subjected to pulse type ground motion, with special focus on the relation between the characteristics of ground motion and the higher‐modes of response. Buildings 10, 20, and 40 stories high were designed such that inelastic deformation was concentrated at a single flexural plastic hinge at their base. Using nonlinear response history analysis, the buildings were subjected to near‐fault seismic ground motions and simple closed‐form pulses, which represented distinct pulses within the ground motions. Euler–Bernoulli beam models with lumped mass and lumped plasticity were used to model the buildings. The response of the buildings to the closed‐form pulses fairly matched that of the near‐fault records. Subsequently, a parametric study was conducted for the buildings subjected to three types of closed‐form pulses with a broad range of periods and amplitudes. The results of the parametric study demonstrate the importance of the ratio of the fundamental period of the structure to the period of the pulse to the excitation of higher modes. The study shows that if the modal response spectrum analysis approach is used — considering the first four modes with a uniform yield reduction factor for all modes, and with the square root of sum of squares modal combination rule — it significantly underestimates bending moment and shear force responses. A response spectrum analysis method that uses different yield reduction factors for the first and the higher modes is presented. Copyright © 2011 John Wiley & Sons, Ltd.     

Probabilistic analysis of peak response of MDOF systems with uncertain PSD function

The use of uniform hazard spectra which have the same probability of exceedance at different frequencies has been proposed for the future version of the National Building Code of Canada. Commonly used combination rules to estimate the peak responses of multi‐degree‐of‐freedom (MDOF) systems are the square root of sum of squares rule and the complete quadratic combination rule. However, the probability that the peak response of a MDOF system exceeds the one estimated by using these rules with the peak modal responses from the uniform hazard spectra cannot be inferred directly. The assessment of the probability of exceedance of the peak response of MDOF systems is presented by considering that the uncertainty in seismic excitation due to all potential earthquakes can be lumped in the power spectral density function of the ground acceleration with uncertain model parameters. This probability is evaluated based on the random vibration of linear systems and the first‐order reliability method. It is found that the under‐ or over‐estimations are less than about 5 or 10% if the modal contributions are not within 10–90% of, or not within 20–80% of, the absolute sum of the effective modal peak responses, respectively. Otherwise, severe under‐ or over‐estimation could result. Copyright © 2002 John Wiley & Sons, Ltd.     

Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading

A method for parametric system identification of classically damped linear system in frequency domain is adopted and extended for non‐classically damped linear systems subjected up to six components of earthquake ground motions. This method is able to work in multi‐input/multi‐output (MIMO) case. The response of a two‐degree‐of‐freedom model with non‐classical damping, excited by one‐component earthquake ground motion, is simulated and used to verify the proposed system identification method in the single‐input/multi‐output case. Also, the records of a 10 storey real building during the Northridge earthquake is used to verify the proposed system identification method in the MIMO case. In this case, at first, a single‐input/multi‐output assumption is considered for the system and modal parameters are identified, then other components of earthquake ground motions are added, respectively, and the modal parameters are identified again. This procedure is repeated until all four components of earthquake ground motions which are measured at the base level of the building are included in the identification process. The results of identification of real building show that consideration of non‐classical damping and inclusion of the multi‐components effect of earthquake ground motions can improve the least‐squares match between the finite Fourier transforms of recorded and calculated acceleration responses. Copyright © 2006 John Wiley & Sons, Ltd.     

Dynamic characteristics of non-classically damped structures

The equations of motion of a structure in undamped modal coordinates may have non-zero off-diagonal terms in the damping matrix. Although these terms are commonly neglected, studies have shown that they may have a significant influence on the response to dynamic loads. In this paper, two independent criteria are developed to determine when these damping terms will affect the structure's modal properties and response. It is found that even small off-diagonal damping values can be significant if the structure has closely spaced natural frequencies. To quantify and understand the influence of these damping terms, closed-form analytical expressions are derived for the modal properties and harmonic and stochastic response of structures with closely spaced natural frequencies. One conclusion is that off-diagonal damping terms will decrease a modal damping ratio for each pair of closely spaced modes. This is significant, since a response analysis performed by neglecting these off-diagonal terms will underestimate the true response.     

Mode-superposition methods in dynamic analysis of classically and non-classically damped linear systems

Mode-superposition analysis is an efficient tool for the evaluation of the response of linear systems subjected to dynamic agencies. Two well-known mode-superposition methods are available in the literature, the mode-displacement method and the mode-acceleration method. Within this frame a method is proposed called a dynamic correction method which evaluates the structural response as the sum of a pseudostatic response, which is the particular solution of the differential equations, and a dynamic correction evaluated using a reduced number of natural modes. The greater accuracy of the proposed method with respect to the other methods is evidenced through extensive numerical tests, for classically and non-classically damped systems.     

Rosenblueth's modal combination rule for systems with non-classical damping

A simple rule is derived to combine, within the framework of a complex mode superposition, the maximum modal responses of systems such as soil-structure and structure-equipment systems, for which closely spaced natural frequencies are likely, and for which, because of the large difference in the damping values of their various components, the assumption of an orthogonal damping matrix may lead to significant errors. The rule constitutes the generalization of Rosenblueth's rule for systems with closely spaced natural frequencies and classical modes, and is expressed in terms of their complex mode shapes and natural frequencies. Its derivation is based on the theory of a complex modal analysis for systems with non-classical modes of vibration and on Rosenblueth's original derivation. As in this original derivation, earthquake ground motions are modelled as a stationary white noise process, but the formulae obtained under this assumption are modified later on to account for the transient nature of actual earthquakes. A numerical example is presented to illustrate the application of the rule, and a comparative study with numerical integration solutions is performed to assess its accuracy. In this comparative study, it predicts the numerical integration solutions with an average error of 0.3 per cent.     

Simplifications of CQC method and CCQC method

The response-spectrum mode superposition method is widely used for seismic response analyses of linear systems. In using this method, the complete quadratic combination （CQC） is adopted for classically damped linear systems and the complex complete quadratic combination （CCQC） formula is adopted for non-classically damped linear systems. However, in both cases, the calculation of seismic response analyses is very time consuming. In this paper, the variation of the modal correlation coefficients of displacement, velocity and displacement-velocity with frequency and damping ratios of two modes of interest are studied, Moreover, the calculation errors generated by using CQC and square-root-of-the-sum-of-thesquares （SRSS） methods （or CCQC and CSRSS methods） for different damping combinations are compared. In these analyses, some boundary lines for classically and non-classically damped systems are plotted to distinguish the allowed minimum frequency ratio at given geometric mean of the damping ratios of both modes if their relativity is neglected. Furthermore, the simplified method, which is a special mode quadratic combination method considering only relativity of adjacent modes in CQC method and named simplified CQC or partial quadratic combination （PQC） method for classically damped linear system, is proposed to improve computational efficiency, and the criterion for determination of how many correlated modes should be adopted is proposed. Similarly, the simplified CCQC or complex partial quadratic combination （CPQC） method for the non-classically damped linear system and the corresponding criterion are also deduced. Finally, a numerical example is given to illustrate the applicability, computational accuracy and efficiency of the PQC and CPQC methods.     

Efficient analysis of dynamic response of linear systems

After reviewing briefly a recently proposed procedure for evaluating the dynamic transient response of a classically damped linear system from its corresponding steady-state response, a modified procedure is presented which also appears to be highly efficient for non-classically damped systems of the type encountered in studies of soil-structure interaction. The concepts involved are developed by reference to viscously damped single-degree-of-freedom systems, and numerical solutions are included to illustrate the accuracy and efficiency of the proposed procedure and its superiority over the classical Discrete Fourier Transform approach.     

Higher-mode effect on the seismic responses of buildings with viscoelastic dampers

In conventional modal analysis procedures, usually only a few dominant modes are required to describe the dynamic behavior of multi-degrees-of-freedom buildings. The number of modes needed in the dynamic analysis depends on the higher-mode contribution to the structural response, which is called the higher-mode effect. The modal analysis approach, however, may not be directly applied to the dynamic analysis of viscoelastically damped buildings. This is because the dynamic properties of the viscoelastic dampers depend on their vibration frequency. Therefore, the structural stiffness and damping contributed from those dampers would be different for each mode. In this study, the higher-mode effect is referred to as the response difference induced by the frequency-dependent property of viscoelastic dampers at higher modes. Modal analysis procedures for buildings with viscoelastic dampers distributed proportionally and non-proportionally to the stiffness of the buildings are developed to consider the higher-mode effect. Numerical studies on shear-type viscoelastically damped building models are conducted to examine the accuracy of the proposed procedures and to investigate the significance of the higher-mode effect on their seismic response. Two damper models are used to estimate the peak damper forces in the proposed procedures. Study results reveal that the higher-mode effect is significant for long-period viscoelastically damped buildings. The higher-mode effect on base shear is less significant than on story acceleration response. Maximum difference of the seismic response usually occurs at the top story. Also, the higher-mode effect may not be reduced by decreasing the damping ratio provided by the viscoelastic dampers. For practical application, it is realized that the linear viscous damping model without considering the higher-mode effect may predict larger damper forces and hence, is on the conservative side. Supported by: Science Council, Chinese Taipei, grant no. 88-2625-2-002-006     

Blind modal identification of output‐only structures in time‐domain based on complexity pursuit

Output‐only modal identification is needed when only structural responses are available. As a powerful unsupervised learning algorithm, blind source separation (BSS) technique is able to recover the hidden sources and the unknown mixing process using only the observed mixtures. This paper proposes a new time‐domain output‐only modal identification method based on a novel BSS learning algorithm, complexity pursuit (CP). The proposed concept—independent ‘physical systems’ living on the modal coordinates—connects the targeted constituent sources (and their mixing process) targeted by the CP learning rule and the modal responses (and the mode matrix), which can then be directly extracted by the CP algorithm from the measured free or ambient system responses. Numerical simulation results show that the CP method realizes accurate and robust modal identification even in the closely spaced mode and the highly damped mode cases subject to non‐stationary ambient excitation and provides excellent approximation to the non‐diagonalizable highly damped (complex) modes. Experimental and real‐world seismic‐excited structure examples are also presented to demonstrate its capability of blindly extracting modal information from system responses. The proposed CP is shown to yield clear physical interpretation in modal identification; it is computational efficient, user‐friendly, and automatic, requiring little expertise interactions for implementations. Copyright © 2013 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.     

EUROCODE 8-Part 2, Bridges: time history presentation of seismic action for base isolated bridges

The purpose of the investigations presented herein has been to show a developed and verified method for presentation of seismic action by a time history. These investigations are related to Eurocode 8, Part 2, Bridges, and have involved the use of recorded accelerogrammes obtained from occurred earthquakes. The method introduces optimization by which, for a period equal to the natural period of vibration of a structure, the ratio between the square root of the sum of squares of the spectra computed from both horizontal components of a single record—SRSS and the code spectrum is computed. While performing this task, pairs of horizontal components for which this ratio is equal or greater than 1.0 (up to a tolerance value) are selected for nonlinear dynamic analysis, provided that the scaling factor is limited. Verification of the method is given for a base isolated bridge on E-75 motorway, with a natural period of 1.4 s. Also, some important recommendations and discussions related to the presentation of the code spectrum by a time history are given.