# An improved multidimensional MPA procedure for bidirectional earthquake excitations.

1. IntroductionIn recent years, a challenge for performance based seismic design is to establish an effective and feasible procedure to evaluate structural seismic capacities. The pushover analysis procedure, due to its simplicity and efficiency, is increasingly applied to estimate seismic demands of structures. During past years, the theory and application of pushover analysis procedure have been developed adequately [1-7], and the technique is accepted by more guideline documents or codes such as FEMA-440 and Vision 2000.

A problem for pushover analysis is that the higher mode effects of building structures cannot be considered in the procedure [8]. In order to solve this, Chopra and Goel [9-11] propose a modal pushover analysis (MPA) procedure to deal with higher mode effects, in which the seismic demand of individual terms of the modal expansion is determined by pushover analysis using the inertia force distribution for each mode, and the modal demands are then combined by the SRSS rule to obtain the total seismic demand. On the basis of this, some improved approach is presented [12-16]. Another problem is that the pushover procedure is restricted with an assumption of single-mode response in which the deformation distributions along the height of structures are fixed. For solving the problem, adaptive pushover methods [17-23] are presented which aim to capture the changes that occur in the vibration properties of structures.

The previous studies about pushover analysis are almost based on symmetric building structures and single-directional earthquake excitation. Both theoretical studies and seismic disasters indicate that asymmetric-plan structures with irregular distributions of mass or stiffness are likely to undergo torsional responses coupled with the translational vibrations, and this type of structures is likely to suffer more severe displacement demands at the corner elements under earthquake excitations. In addition, the torsion coupling response will induce the structural space effects that cannot be solved in two-dimensional analysis. So, several research efforts have been made to extend and apply the pushover analysis to asymmetric-plan structures whose inelastic seismic responses are intricate [24-26]. At present, the MPA [27-31] is extended to the asymmetric-plan structures in which the inertia force distribution for each mode includes two lateral forces and torque at each floor level, and the CQC rule is used to combine the modal demands to obtain the total demand of asymmetric-plan structures.

For multidimensional MPA, which is applicable to the analysis of asymmetric-plan structures, the modal equivalent single-degree-of-freedom (ESDOF) system is actually subjected to the superposition of bidirectional earthquake excitations. In the current MPA procedures, the response of ESDOF system is calculated according to single-directional earthquake excitation, the x- and y-axes, respectively, and peak response of two directions is combined through SRSS rule, which is applicable to elastic analysis but not applicable to inelastic analysis for having no regard for inelastic coupling of structure responses caused by earthquake excitations of different components. For solving the problem, an improved multidimensional MPA (IMMPA) idea is presented in the paper: superposed bidirectional earthquake acceleration acting on modal ESDOF system can be regarded as an unidirectional "combined earthquake excitation"; in the procedure, to solve structural responses, seismic demand spectra may be established based on SDOF system subjected to the "combined earthquake excitation."

The main objective of this paper is to propose an MPA procedure for evaluating the seismic capacities of asymmetric-plan structures subjected to bidirectional earthquake excitations. An outline of this paper can be expressed as follows. First, the improved procedure is presented based on the assumption of modal equivalent SDOF system which is subjected to the superposition of bidirectional earthquake excitations. Then, the corresponding strength reduction factor spectra are discussed and compared with traditional spectra. Next, the step-by-step procedure is proposed. At the end, two examples are considered to verify the accuracy of the proposed procedure.

2. Motion Equations of Equivalent Systems

Consider an n-storey building. Its plane is asymmetric about the x- or/and y-axes. Each floor diaphragm is rigid in its own plane and has three degrees of freedom defined at the center of mass (CM). The equation of motion governing the response of the n-storey asymmetric structure subjected to bidirectional earthquake excitations (along the x- and y-axes) is as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where M is diagonal mass matrix of order 3n which includes three submatrices m, m, and [I.sub.0]; m is associated with x and y-lateral degrees of freedom, and [I.sub.0] is associated with torsion degrees of freedom. C is the damping matrix and F(u) is the vector of resisting forces. The displacement vector u(f) is equal to [[[u.sub.x] [u.sub.y] [u.sub.[theta]]].sup.T]; [u.sub.x], [u.sub.y], and [u.sub.[theta]] are the x, y, and torsion-directional displacement subvectors, respectively. The influence vectors, [i.sub.x] and [i.sub.y], are as follows:

[i.sub.x] = [{1} {0} {0}], [i.sub.y] = [{0} {1} {0}]. (2)

In MPA procedure, a major assumption is that the responses of structures can be expressed as superposition of the responses of appropriate SDOF systems just like that in the linear range. The displacement vector of the inelastic system can be expanded in terms of the natural vibration modes of the corresponding linear system:

u (t) = [PHI]q = [3n.summation over (i=1)][[phi].sub.i][q.sub.i] (t), (3)

in which the ith modal vector [[phi].sub.i] of size 3n x 1 includes three n x 1 subvectors [[phi].sub.xi], [[phi].sub.yi], and [y.sub.[theta]i]. Substituting u(f) = [PHI]q into (1) and premultiplying both sides of (1) by [[phi].sup.T.sub.i], the following equations are given:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)

where the ith modal resisting force quantity, [F.sub.si](t), equals [[phi].sup.T.sub.i][K.sub.ep](t)[[phi].sub.i][q.sub.i](t); [K.sub.ep](t) is elastic-plastic instantaneous stiffness matrices; [[phi].sup.T.sub.i]M[[phi].sub.i] is interpreted as the ith modal equivalent mass, [M.sub.i]; Modal participation factors, [[GAMMA].sub.xn] and [[GAMMA].sub.yn], along the x- and y-axes, respectively, are expressed as

[[GAMMA].sub.xi] = [[phi].sup.T.sub.i]M[i.sub.x]/[[phi].sup.T.sub.i]M[[phi].sub.i], [[GAMMA].sub.yi] = [[phi].sup.T.sub.i]M[i.sub.y]/[[phi].sup.T.sub.i]M[[phi].sub.i], i = 1, ..., 3n. (5)

Dividing both sides of (4) by modal equivalent mass, [M.sub.i], and using the orthogonality property of modes, 3n uncoupled equations can be derived as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6)

which are the motion equations of equivalent systems from 1th to 3nth mode of the considered structure. In (6), [[omega].sub.i] is natural vibration frequency of the ith modal equivalent system of the structure and [xi] is damping ratio. The relationship between the resisting force parameter, [F.sub.si]/[M.sub.si], and the modal coordinate, [q.sub.i], of the equivalent systems can be estimated from an idealized base shear and roof displacement pushover curve of the nth mode. The equivalent displacement responses, [q.sub.1](t) ... [q.sub.i](t) ... [q.sub.3n](t), can be calculated by (6), and then the displacement vector, u(t), can be solved by 3). Transposing M[??](t) from the left side of (1) to the right side, the following is given:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

in which [P.sub.ine] can be interpreted as earthquake inertia force vector of the structure. Substituting (3) into (7), [P.sub.ine] can be expressed as a summation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

in which [p.sub.ine,i] is the ith modal inertia force vector of the structure.

3. Traditional Multidimensional MPA

In the traditional multidimensional MPA procedure, the peak modal responses [r.sub.x] and [r.sub.y] to x- and y-components of earthquake motion are determined by nonlinear static analysis, respectively. Associated with components of earthquake motion, the total inertia force vector can be expanded as a summation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)

In the procedure, the contribution of responses for one component (x or y) of ground motion is assumed to be uncoupled, and the uncoupled equation of motion is expressed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

The peak responses of x- and y-directions are determined, respectively, by combining the peak modal responses using the CQC rule as follows:

[r.sub.x] = [([summation over i][summation over j] [[rho].sub.ij][r.sub.x,i][r.sub.x,j]).sup.1/2]; [r.sub.y] = [([summation over i][summation over j] [[rho].sub.ij][r.sub.y,i][r.sub.y,j]).sup.1/2]. (11)

The responses [r.sub.x] and [r.sub.y] are combined to determine the total response r by the SRSS rule:

r = [square root of ([r.sup.2.sub.x] + [r.sup.2.sub.y])]. (12)

In the procedure, the responses of different horizontal components of earthquake excitation are uncoupled. The superposition rule is utilized here to solve the problem which is accurate in linear response stage of the structures under bidirectional ground motions but it is not applicable in nonlinear stage.

4. Improved Multidimensional MPA (IMMPA)

4.1. Equivalent SDOF System under Superposition Excitation. If abs [absolute value of ([[GAMMA].sub.x,i])] [greater than or equal to] abs [absolute value of ([[GAMMA].sub.y,i])], extract [[GAMMA].sub.x,i] from the right side of equation and define [[GAMMA].sub.x,i] = [[GAMMA].sub.i]. The motion equation of equivalent systems is then transformed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

In addition, if abs [absolute value of ([[GAMMA].sub.y,i])] [greater than or equal to] abs [absolute value of ([[GAMMA].sub.x,i])], a similar equation can be established, and it is omitted here because of theoretical consistency. In (13), the parameter [[gamma].sub.i] equals [[GAMMA].sub.y,i]/[[GAMMA].sub.x,i] which is in the range from -1 to 1. The superposition of bicomponent ground acceleration time history, [[??].sub.gx](t) + [gamma][[??].sub.gy](t), in the right side of (13) can be interpreted as a single-component acceleration time history, [[??].sub.gi](t), and the corresponding motion equations are reduced to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

Defining [d.sub.i](t) = [q.sub.i](t)/[[GAMMA].sub.i], (14) is transformed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (15)

where [d.sub.i](t) is governed by the equation of motion for the ith modal equivalent SDOF system subjected to earthquake excitation, [[??].sub.gi](t). The resisting force term, [f.sub.i](t), equals [F.sub.si](t)/[M.sub.i][[GAMMA].sub.i]. Equation (15) shows that, due to the aforementioned assumptions, the nonlinear response of a structure system with 3n degrees of freedom subjected to bidirectional earthquake excitations can be expressed as the sum of the responses of 3n equivalent SDOF systems under single-directional excitation, [[??].sub.gx](t) + [[gamma].sub.i][[??].sub.gy](t), each one corresponding to a vibration "mode"

According to the derivation above, the expression of inertia force vector, Pine, can be expanded as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (16)

where the modal inertia force vector, [p.sub.ine, i], can be reduced as

[p.sub.ine] = [s.sub.i]A (17)

in which [s.sub.i] is distribution factor vector of inertia force, and [s.sub.i] equals

[s.sub.i] = [[GAMMA].sub.i][[m[[phi].sub.xi] m[[phi].sub.yi] [I.sub.0][[phi].sub.[theta]i]].sup.T]. (18)

The factor vector, [s.sub.i], can be interpreted as the distribution

factor of modal pushover analysis.

4.2. Inelastic Response Demands for Equivalent Systems under Superposition Excitation. Generally, the response demands of inelastic systems can be obtained by inelastic response spectra, such as strength reduction factor spectra (i.e., R-[mu]-T relation [32]).

The yield resisting force expression [f.sub.i,y] = [k.sub.x][d.sub.i,y] is defined in which [k.sub.x] is linear stiffness and [d.sub.i,y] is yield displacement, and the relationship for [[mu].sub.i](t) = [d.sub.i](t)/[d.sub.y] is defined, and (15) can be then normalized as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (19)

where [beta]([[omega].sub.i], [xi]) is amplification coefficient spectrum corresponding to fixed damping ratio, [xi], natural vibration frequency, [[omega].sub.i], and excitation; [R.sub.bi] is strength reduction factor; [R.sub.bi] equals [f.sub.e,i](t)/[f.sub.i,y](t) in which [f.sub.e,i](t) is elastic resisting force of corresponding elastic equivalent system. By solving (19), the maximal value of [[mu].sub.i](t) is ductility factor, [[mu].sub.i], that is defined as the ratio between maximal displacement, [d.sub.i,max], and yield displacement, [d.sub.y]. Based on 19), the relationship for [R.sub.bi]-[[mu].sub.i]-[T.sub.i] ([T.sub.i] = 2[pi]/[[omega].sub.i]) can be established through iterative calculation.

On the other hand, for a SDOF system, according to the definitions of strength reduction factor, [R.sub.b], and ductility factor, [mu], the relationship is given as

[d.sub.p,max] = [[omega].sup.2] x [A.sub.e] x [mu]/[R.sub.b], (20)

where [A.sub.e] is elastic maximum acceleration response (or elastic acceleration spectra); [d.sub.p,max] is inelastic maximum displacement response (or inelastic displacement spectra). Substituting [A.sub.e] and "[R.sub.b]-[mu]-T" model into (19), then the inelastic maximum displacement, [d.sub.p,max], of the SDOF system can be solved.

During the past years, some researchers have presented convenient [R.sub.s] spectra models [32-34]. To combine with current [R.sub.s] spectra, it is recommended in this paper to establish [R.sub.b] spectra by analyzing the relationship between [R.sub.s] spectra and [R.sub.b] spectra. The approach is as follows: establish [R.sub.b]/[R.sub.s] ratio spectra based on different soil sites, analyze the influences of soil site classification, ductility factor, and participation coefficient ratio of vibration mode of two spindle axes, and establish simplified relationship of [R.sub.b]/[R.sub.s]. On this basis, establish [R.sub.b] spectra with [R.sub.s] spectra and simplified relation of [R.sub.b]/[R.sub.s]. The [R.sub.s] spectra for three kinds of soil sites are shown in Figure 1.

The [R.sub.b]/[R.sub.s] spectra ([mu] = 2 and 4, three kinds of soil sites) are shown in Figure 2, in which [mu] is ductility factor and [gamma] is participation coefficient ratio of vibration mode of two spindle axes. The following can be drawn from [R.sub.b]/[R.sub.s] spectra analysis. (1) The values of [R.sub.b]/[R.sub.s] fluctuate around 1 and have no significant regularity as the period increases. (2) The regularity of [R.sub.b]/[R.sub.s] for [gamma] > 0 and that for [gamma] < 0 are different from each other, so it is suggested to analyze them, respectively.

Through the above analysis, the regularity of [R.sub.b]/[R.sub.s] along with period change is not very significant, so it is suggested to take the same [R.sub.b]/[R.sub.s] spectra value within the entire period. In addition, for the reason of disadvantage in [R.sub.b]/[R.sub.s] < 1, "mean value - 1.65 x standard deviation" of [R.sub.b]/[R.sub.s] within 0.05~5.0 seconds is defined as modification coefficient, [alpha], of traditional [R.sub.s]. The distribution curves of modification coefficient, [alpha], along with ductility factor changes are shown in Figure 3. The following can be drawn from Figure 3. (1) For all kinds of soil sites, the minimum value of a is 0.910 and the maximum value is 0.985, which shows that the influence of the analysis variables on [R.sub.s] spectra is little. (2) For [gamma] > 0, a decreases while [gamma] increases; for [gamma] < 0, [alpha] decreases while [gamma] decreases (it is not so significant for [gamma] < 0 at hard soil site and for [gamma] < 0 at intermediate soil site). (3) For hard soil site, [alpha] is maximized while [mu] = 1.5; [alpha] decreases and does not change significantly along with increase of ductility factor while [mu] [greater than or equal to] 2; for intermediate soil site, a increases along with increase of the ductility factor while [mu] [greater than or equal to] 2; for soft soil site, [alpha] decreases along with increase of ductility factor.

4.3. Step-by-Step Procedure. A step-by-step summary of the IMMPA procedure is presented as follows.

(1) The natural frequencies, [[omega].sub.i], and modes, [[phi].sub.i], for elastic vibration of the multidimensional building are calculated and the corresponding mode shapes are normalized.

(2) For the ith mode, determine the base shear-roof displacement, [V.sub.bi]-[u.sub.ri], pushover curve by multidimensional nonlinear static analysis applying the force distribution (17). [V.sub.bi] and [u.sub.ri] are chosen to correspond to the direction of the x- and y-components, respectively. The pushover curve is in the direction of the dominant excitation in the mode being considered.

(3) Idealize the [V.sub.bi]-[u.sub.ri] pushover curve as a bilinear A-D curve by utilizing the well-known procedure [35].

(4) Determine [A.sub.y]-D earthquake demand spectra curve based on the nonlinear dynamic analysis of modal equivalent SDOF system subjected to the superposition excitation.

(5) Determine the roof object displacement of ith mode and then calculate the other ith modal responses (such as drifts, plastic rotations, etc.) by conducting pushover analysis up to the already calculated object displacement.

(6) Repeat steps 1-5 for as many modes as required for sufficient accuracy.

(7) Determine the total extreme responses, r, (including [r.sub.x]- and [r.sub.y]-components) for the excitation of the combined earthquake motion by combining gravity response and the peak modal responses by utilizing the CQC rule.

5. Numerical Example for IMMPA Procedure

5.1. Description of the Example Buildings. In order to clarify how the proposed methodology should be applied, two simple analytical examples are presented. The structures considered are a 10-storey and a 15-storey eccentric reinforced concrete frame buildings, as illustrated in Figure 4. It is considered that the ground motion is acting simultaneously along the two horizontal axes. Each floor has three degrees of freedom (DOF) defined at the center of mass (CM). Each floor diaphragm is rigid in its own plane and the center of stiffness (CS) deviates from the CM. The dimensions of the two buildings in height are 3.6 m for all stories, and the dimensions in plane are shown as Figure 4. For the 10-storey building, the sectional sizes of beams are 300 x 700 [mm.sup.2] for 2th~11th floor; the sectional sizes of columns are 700 x 700 [mm.sup.2] for 1th~3th storey, 600 x 600 [mm.sup.2] for 4th~7th storey, and 500 x 500 [mm.sup.2] for 8th~10th storey. For the 15-storey building, the sectional sizes of beams are 300 x 700 [mm.sup.2] for 2th~5th floor, 300 x 650 [mm.sup.2] for 6th~11th floor, and 300 x 600 [mm.sup.2] for 12th~16th floor; the sectional sizes of columns are 700 x 700 [mm.sup.2] for 1th~4th storey, 650 x 650 [mm.sup.2] for 5th~10th storey, and 500 x 500 [mm.sup.2] for Thh~15th storey. Steel ratios are approximately 15% for beam sections and 2% for column sections of the two buildings. Concrete compression strength is selected as 30 MPa for all columns and beams of the structure. The design dead load and live load are, respectively, 6.6kN/[m.sup.2] (4.7kN/[m.sup.2]) and 1.0kN/[m.sup.2] (2.0kN/[m.sup.2]) for each floor (roof). The damping of the building is modeled by the Rayleigh damping, and damping ratio [xi] equals 5%.

The bidirectional "Taft" earthquake records (Kern County 1952,1095 Taft Lincoln School) are selected to verify the accuracy of the presented method. The acceleration peak values of the records are adjusted to 3.1m/[sec.sup.2] and 4.0 m/[sec.sup.2], respectively. The nonlinear response time history analyses (NL-RHA) are implemented for the two example building subjected to the selected earthquake record and the time history responses of roof displacements are illustrated in Figure 5, in which 8 peak displacement points are defined as 8 objective control points for variance analysis.

5.2. Comparative Evaluation Method. The method to compare estimation results of IMMPA with exact results of NLRHA for verifying the procedure is as follows.

(1) The restoring force model of modal equivalent SDOF system is established by modal pushover analysis.

(2) The dynamic response [Q.sub.is] of the ith mode equivalent SDOF system is solved by NL-RHA, in which the restoring force model obtained in step 1 is used. The dynamic response [Q.sub.is] is converted into the ith modal responses [Q.sub.i], and total structural dynamic responses Q can be obtained by superposition of modal responses; that is, Q = [summation] [Q.sub.i].

(3) Corresponding to the eight objective control points given in Section 5.1, the modal roof objective displacements, [u.sub.roi], and structural roof objective displacements, [u.sub.ro], are solved, respectively, according to step 2.

(4) The IMMPA is implemented and the modal static responses [Q.sup.s.sub.i], corresponding to [u.sub.roi], are obtained. The total static responses, [Q.sup.s], are determined by linear combining the response values of each modal pushover analysis; that is, [Q.sup.s] = [sumamtion] [Q.sup.s.sub.i].

(5) The demand responses are compared between IMMPA and NL-RHA.

5.3. Comparison of Responses. The comparisons of story drift ratios and floor displacements between NL-RHA method and IMMPA method are illustrated in Figure 6. The values gotten with NL-RHA method are considered to be the exact one, which will be utilized to analyze the estimation accuracy of structural responses with IMMPA method. It can be drawn from Figure 6 that (1) in general IMMPA method can be used to estimate deformation distribution of structures under bidirectional earthquake excitations, rationally. (2) With peak value of earthquake acceleration and structural plastic deformation increase, the deviation of structural response estimated with MMPA may increase. The reason is that the concept of elastic modal decomposition is not applicable to inelastic response analysis in theory. In addition, it is not rational to estimate dynamic response of structures with static analysis method, which may also result in deviation. (3) By comparing 10-storey structure and 15-storey structure in the case, it is found that the response deviation of the 15storey structure, estimated with IMMPA, is larger than that of the 10-storey structure, so it can be drawn that structural deviation estimated with IMMPA will increase as the number of total floors of structure increases. The reason is that, with the increase of the number of total floors, the vibration modes are gradually complex and the effect of the higher mode on structural response increases gradually.

6. Conclusion

A new IMMPA idea is presented in the paper based on traditional MPA method: superposed bidirectional earthquake excitation acting on modal ESDOF system can be regarded as an unidirectional "combined earthquake excitation"; in pushover procedure, the static force replacing "combined earthquake excitation" is assigned to three components of each floor of a structure based on model of vibration.

In accordance with this idea, in solution of structural displacement response, it is required to input "combined earthquake acceleration" to modal ESDOF system in order to establish [R.sub.b] spectra and convert it to seismic demand spectra. It is recommended in this paper to get modification coefficient of [R.sub.s] spectra by analyzing the relation of [R.sub.b]/[R.sub.s]. Define modification coefficient, [alpha], as "mean value - 1.65 x standard deviation" of [R.sub.b]/[R.sub.s] within 0.05~5 seconds. Based on selected earthquake records, analyze the influence of soil site classification, ductility factor, and participation coefficient ratio of vibration mode on [alpha].

In order to check accuracy of the IMMPA procedure presented in this paper, two structures with mass eccentricity are designed, a bidirectional earthquake ground motion, Taft, is taken as earthquake excitation of the structural system, and deformation distribution of the structure is calculated, respectively, with NL-RHA and IMMPA, and the results are concluded as follows after comparative analysis.

(1) In general, IMMPA method can be used to accurately estimate the distribution feature of plastic deformation of the structure along its floors.

(2) Along with increase of peak of earthquake acceleration, structural response deviation estimated with IMMPA method may also increase.

(3) Along with increase of the number of total floors of structures, structural response deviation estimated with IMMPA method may also increase.

http://dx.doi.org/10.1155/2014/320756

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research work was jointly supported by the National Natural Science Foundation of China (Grant no. 51108067), the "973" Project (Grand no. 2011CB013605), and Dalian Nationalities University project (Grant no. DC120101095).

References

[1] P. Fajfar and P. Gaspersic, "The N2 method for the seismic damage analysis of RC buildings," Earthquake Engineering and Structural Dynamics, vol. 25, no. 1, pp. 31-46, 1996.

[2] P. Fajfar, D. Marusic, and I. Perus, "Torsional effects in the pushover-based seismic analysis of buildings," Journal of Earthquake Engineering, vol. 9, no. 6, pp. 831-854, 2005.

[3] E. Kalkan and S. K. Kunnath, "Assessment of current nonlinear static procedures for seismic evaluation of buildings," Engineering Structures, vol. 29, no. 3, pp. 305-316, 2007

[4] E. Erduran, "Assessment of current nonlinear static procedures on the estimation of torsional effects in low-rise frame buildings," Engineering Structures, vol. 30, no. 9, pp. 2548-2558, 2008.

[5] G. Manoukas, A. Athanatopoulou, and I. Avramidis, "Static pushover analysis based on an energy-equivalent SDOF system," Earthquake Spectra, vol. 27, no. 1, pp. 89-105, 2011.

[6] A. Fiore, F. Porco, D. Raffaele, and G. Uva, "About the influence of the infill panels over the collapse mechanisms actived under pushover analyses: two case studies," Soil Dynamics and Earthquake Engineering, vol. 39, pp. 11-22, 2012.

[7] T.-H. Yi, H.-N. Li, and X.-Y. Zhao, "Noise smoothing for structural vibration test signals using an improved wavelethresholding technique," Sensors, vol. 12, no. 8, pp. 11205-11220, 2012.

[8] M. Kreslin and P Fajfar, "The extended N2 method taking into account higher mode effects in elevation," Earthquake Engineering and Structural Dynamics, vol. 40, no. 14, pp. 1571-1589, 2011.

[9] A. K. Chopra and R. K. Goel, "A modal pushover analysis procedure for estimating seismic demands for buildings," Earthquake Engineering and Structural Dynamics, vol. 31, no. 3, pp. 561-582, 2002.

[10] A. K. Chopra, R. K. Goel, and C. Chintanapakdee, "Evaluation ofa modified MPA procedure assuming higher modes as elastic to estimate seismic demands," Earthquake Spectra, vol. 20, no. 3, pp. 757-778, 2004.

[11] R. K. Goel and A. K. Chopra, "Evaluation of modal and FEMA pushover analyses: SAC buildings," Earthquake Spectra, vol. 20, no. 1, pp. 225-254, 2004.

[12] J. Mao, C. Zhai, and L. Xie, "An improved modal pushover analysis procedure for estimating seismic demands of structures," Earthquake Engineering and Engineering Vibration, vol. 7, no. 1, pp. 25-31, 2008.

[13] M. Poursha, F. Khoshnoudian, and A. S. Moghadam, "A consecutive modal pushover procedure for estimating the seismic demands of tall buildings," Engineering Structures, vol. 31, no. 2, pp. 591-599, 2009.

[14] H.-N. Li, T.-H. Yi, M. Gu, and L.-S. Huo, "Evaluation of earthquake-induced structural damages by wavelet transform," Progress in Natural Science, vol. 19, no. 4, pp. 461-470, 2009.

[15] Y. Jiang, G. Li Gang, and D. Yang, "A modified approach of energy balance concept based multimode pushover analysis to estimate seismic demands for buildings," Engineering Structures, vol. 32, no. 5, pp. 1272-1283, 2010.

[16] S. W. Han, K. H. Moon, and A. K. Chopra, "Application of MPA to estimate probability of collapse of structures," Earthquake Engineering and Structural Dynamics, vol. 39, no. 11, pp. 1259-1278, 2010.

[17] S. Antoniou and R. Pinho, "Development and verification of a displacement-based adaptive pushover procedure," Journal of Earthquake Engineering, vol. 8, no. 5, pp. 643-661, 2004.

[18] E. Kalkan and S. K. Kunnath, "Adaptive modal combination procedure for nonlinear static analysis of building structures," Journal of Structural Engineering, vol. 132, no. 11, Article ID 002611QST, pp. 1721-1731, 2006.

[19] C. Casarotti and R. Pinho, "An adaptive capacity spectrum method for assessment of bridges subjected to earthquake action," Bulletin of Earthquake Engineering, vol. 5, no. 3, pp. 377-390, 2007.

[20] S. Kim and Y. C. Kurama, "An alternative pushover analysis procedure to estimate seismic displacement demands," Engineering Structures, vol. 30, no. 12, pp. 3793-3807, 2008.

[21] K. Shakeri, M. A. Shayanfar, and T. Kabeyasawa, "A story shear-based adaptive pushover procedure for estimating seismic demands of buildings," Engineering Structures, vol. 32, no. 1, pp. 174-183, 2010.

[22] R. Bento, C. Bhatt, and R. Pinho, "Using nonlinear static procedures for seismic assessment of the 3D irregular SPEAR building," Earthquake and Structures, vol. 1, no. 2, pp. 177-195, 2010.

[23] K. Shakeri, K. Tarbali, and M. Mohebbi, "An adaptive modal pushover procedure for asymmetric-plan buildings," Engineering Structures, vol. 36, pp. 160-172, 2012.

[24] D. Marusic and P Fajfar, "On the inelastic seismic response of asymmetric buildings under bi-axial excitation," Earthquake Engineering and Structural Dynamics, vol. 34, no. 8, pp. 943-963, 2005.

[25] G. Adhikari and R. Pinho, "Development and application of nonlinear static procedures for plan-asymmetric buildings," Tech. Rep. ROSE-2010/01, IUSS Press, Pavia, Italy, 2010.

[26] T.-H. Yi, H.-N. Li, and H.-M. Sun, "Multi-stage structural damage diagnosis method based on "energy-damage" theory," Smart Structures and Systems, vol. 12, no. 3-4, pp. 345-361, 2013.

[27] A. K. Chopra and R. K. Goel, "A modal pushover analysis procedure to estimate seismic demands for unsymmetric-plan buildings," Earthquake Engineering and Structural Dynamics, vol. 33, no. 8, pp. 903-927, 2004.

[28] J. C. Reyes and A. K. Chopra, "Three-dimensional modal pushover analysis of buildings subjected to two components of ground motion, including its evaluation for tall buildings," Earthquake Engineering and Structural Dynamics, vol. 40, no. 7, pp. 789-806, 2011.

[29] J. C. Reyes and A. K. Chopra, "Evaluation of three-dimensional modal pushover analysis for unsymmetric-plan buildings subjected to two components of ground motion," Earthquake Engineering and Structural Dynamics, vol. 40, no. 13, pp. 1475-1494, 2011.

[30] M. Poursha, F. Khoshnoudian, and A. S. Moghadam, "A consecutive modal pushover procedure for nonlinear static analysis of one-way unsymmetric-plan tall building structures," Engineering Structures, vol. 33, no. 9, pp. 2417-2434, 2011.

[31] G. Manoukas, A. Athanatopoulou, and I. Avramidis, "Multimode pushover analysis for asymmetric buildings under biaxial seismic excitation based on a new concept of the equivalent single degree of freedom system," Soil Dynamics and Earthquake Engineering, vol. 38, pp. 88-96, 2012.

[32] T. Vidic, P Fajfar, and M. Fischinger, "Consistent inelastic design spectra: strength and displacement," Earthquake Engineering & Structural Dynamics, vol. 23, no. 5, pp. 507-521, 1994.

[33] B. Borzi, G. M. Calvi, A. S. Elnashai, E. Faccioli, and J. J. Bommer, "Inelastic spectra for displacement-based seismic design," Soil Dynamics and Earthquake Engineering, vol. 21, no. 1, pp. 47-61, 2001.

[34] R. Riddell, J. E. Garcia, and E. Garces, "Inelastic deformation response of SDOF systems subjected to earthquakes," Earthquake Engineering and Structural Dynamics, vol. 31, no. 3, pp. 515-538, 2002.

[35] P. Fajfar, "Capacity spectrum method based on inelastic demand spectra," Earthquake Engineering and Structural Dynamics, vol. 28, pp. 979-993, 1999.

Feng Wang, (1) Jian-Gang Sun, (1) and Ning Zhang (2)

(1) College of Architecture & Civil Engineering, Dalian Nationalities University, Dalian 116600, China

(2) Dalian Polytechnic University, Dalian 116003, China

Correspondence should be addressed to Feng Wang; win_0803@163.com

Received 30 April 2014; Accepted 26 June 2014; Published 20 July 2014

Academic Editor: Xiao-Wei Ye

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Wang, Feng; Sun, Jian-Gang; Zhang, Ning |

Publication: | The Scientific World Journal |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 5365 |

Previous Article: | Morphology and magnetic properties of ferriferous two-phase sodium borosilicate glasses. |

Next Article: | The burden of hypertension and kidney disease in northeast India: the institute for Indian mother and child noncommunicable diseases project. |

Topics: |