Comparison and Selection of Ground Motion Intensity Measures for Seismic Design of Super HighRise Buildings Xiao Lu^{1}, Xinzheng Lu^{1,} [*] , Hong Guan^{2} and Lieping Ye^{1}^{} ^{1} Key Laboratory of Civil Engineering Safety and Durability of China Education Ministry, Department of Civil Engineering, Tsinghua University, Beijing, P.R. China, 100084; ^{2}Griffith School of Engineering, Griffith University Gold Coast Campus, Queensland 4222, Australia Advances in Structural Engineering, 2013, 16(7): 12491262. Abstract: Ground motion intensity measures (IMs) are an important basis of structural seismic design. Extensive research has been conducted on the selection of IMs for building structures with fundamental periods shorter than 6 s. However, minimal work has been performed for super highrise buildings whose fundamental periods are much longer than 6 s. To fill the gap in this research area, this paper aims to develop a simplified analytical model for super highrise buildings based on the flexuralshear coupled beam model. The variation in correlation and dispersion between different IMs is evaluated and the structural seismic response demand measures (DMs) are analyzed for different structural fundamental periods. Subsequently, rational IMs for the seismic design of super highrise buildings are suggested. In addition, the influence of different flexural and shear stiffness ratios, a_{0}, on the selection of an IM for super highrise buildings is also discussed. Finally, a series of incremental dynamic analyses (IDA) of the Shanghai Tower, with a total height of approximately 632 m, is performed to verify the rationality of using different IMs for super highrise buildings. The numerical results indicate that with a minimum dispersion, the peak ground velocity (PGV) has a better correlation to the story drift ratio than any other IMs. While considering structural nonlinearity, PGV still yields the minimum coefficient of variation in the collapse analysis of the actual super highrise building. It is therefore recommended that PGV be used as an IM for the seismic design of super highrise buildings. Key words: super highrise building; collapse analysis; ground motion intensity measure; simplified model. DOI:10.1260/13694332.16.7.1249 If you need the PDF version of this paper, please email to luxinzheng@sina.com

1. INTRODUCTION With the advancement of structural materials and construction technology, super highrise building construction has entered into a new period of vigorous development since the completion of the first super highrise building taller than 500 mTaipei 101 in 2004. The 2012 statistical data from the Council on Tall Buildings and Urban Habitats (CTBUH) (http://buildingdb.ctbuh.org/) indicates that by the end of 2012, a total of thirteen super highrise buildings of 500 m or taller have been constructed or under construction worldwide. This necessitates further studies to understand the seismic performance of such super highrise buildings. One of the most important issues in seismic design is the ground motion intensity measure, which relates the seismic hazards to the structural responses. Selection of a sufficient and efficient IM can effectively reduce the deviations of the structural responses. In view of this, extensive research has been performed worldwide through which a number of IMs have been proposed (Tothong and Luco 2007; Cordova et al. 2001; Luco and Cornell 2007; Baker and Cornell 2005; Vamvatsikos and Cornell 2005; Ye et al. 2013). However, the application of these proposed IMs in seismic analysis and design has largely been based on normal structures whose fundamental periods are shorter than 6 s and are mostly in the range of 04 s. On the other hand, little work has been reported with regard to IMs for super highrise buildings. One of the most distinguished features of super highrise buildings is that their fundamental periods are much longer than those of the low to mediumrise frame or shear wall buildings. For example, the fundamental periods of Shanghai Tower and PingAn Finance Center, two super highrise buildings still under construction in China, will be larger than 9.0 s (Lu et al, 2011, Yang et al, 2011). More importantly, the structural displacement responses of super highrise buildings are not governed by the fundamental vibration mode. Instead, the higherorder vibration modes contribute more significantly to the displacement responses. As such, a systematic evaluation of the rationality of the existing IMs for seismic design of super highrise buildings is both necessary and significant for engineering practice. In this study, a simplified model combining a flexural and a shear beam is developed based on the mechanical properties of an actual super highrise buildingthe Shanghai Tower. Timehistory analysis is also implemented to discuss the correlations between some existing IMs and structural seismic response demand measures (DMs), as well as their corresponding dispersions. A total of 22 pairs of farfield ground motion records suggested by FEMA P695 (FEMA 2009) are adopted for the timehistory analysis. Subsequently, rational IMs for the seismic design of super highrise buildings are suggested. Finally, the rationality of the recommended IMs when considering structural nonlinearity is verified via a series of incremental dynamic collapse analyses of the Shanghai Tower. 
2. EXISTING INTENSITY MEASURES Published literature (Tothong and Luco 2007; Cordova et al. 2001; Luco and Cornell 2007; Baker and Cornell 2005; Vamvatsikos and Cornell 2005; Ye et al. 2013) indicates that the existing IMs can be classified into two types: scalar and vectorvalued. Scalar IMs can be further categorized into singleparameter and multiparameter IMs according to the number of parameters required. The traditional peak values of ground motion records such as PGA, PGV and PGD are typical singleparameter scalar IMs. These IMs are still widely used in many national design codes because of their simplicity (MOC 2010; The Building Centre of Japan, 2000). Some recent literature indicates that PGA or PGV IMs, which only consider the characteristics of ground motions, do not perform well in predicting structural responses within short to medium periods. Therefore, the 5%damped elastic spectral acceleration S_{a}(T_{1}) at fundamental period T_{1 }(Shome et al. 1998) has been widely used for these periods. S_{a}(T_{1}) covers not only the characteristics of ground motions but also the dynamic features of the structures. Compared with PGA, the dispersion of the predicted structural responses based on S_{a}(T_{1}) is remarkably reduced (Ye et al. 2013), particularly for firstmodedominated, moderate period structures. Despite these advancements, S_{a}(T_{1}) only deals with the characteristics of the elastic fundamental period of the structure. When progressing into the nonlinear stage, the fundamental period of a structure elongates gradually, leading to changes in the structural dynamic characteristics. Furthermore, for highrise buildings, higherorder vibration modes greatly influence the structural responses. All these factors however were not considered in the S_{a}(T_{1}) formulation (Lucchini et al. 2011). To overcome the aforementioned limitations, various attempts have been made to improve S_{a}(T_{1}) by considering the effect of period elongation and the contribution of the higher modes (Tothong and Luco 2007; Cordova et al. 2001; Vamvatsikos and Cornell 2005). For example, Cordova et al. (2001) proposed an IM called S * a, which is based on the fundamental period T_{1 }and the longer period T_{f} representing the inelastic (damaged) structure S * a, also taking into consideration the period elongation in the nonlinear stage, as expressed by Eqn 1 below: _{ } , (1) where a and c are two parameters to be calibrated. Note that the dynamic behavior of buildings upon entering into the nonlinear stage is not fully understood, particularly at the early stage of structural design because of the difficulties in accurately predicting the nonlinear behavior of the building. Consequently, IMs that consider period elongation are limited in scope and practical application. With respect to the higher vibration modes, the IM proposed by Vamvatsikos and Cornell (2005) also considers the participation of higherorder vibration modes. In the expression_{ } , the values of t_{a} and t_{b} represent the higher modal periods of concern. Similarly, Luco and Cornell^{ }(2007)^{ }proposed an IM called IM_{1E&2E} based on the first two vibration modes and the squarerootofsumofsquares (SRSS) rule of modal combination. This IM considers the effect of the contribution of higher modes caused by the nearfault pulse, as expressed in Eqn 2. _{ } , (2) where_{ } and S_{d}(T_{1},z_{1}) and S_{d}(T_{2},z_{2}) are the elastic spectral displacements, respectively, at damping ratios z_{1} and z_{2} with corresponding periods T_{1} and T_{2}. Eqn 2 also defines the j^{th }¨C mode participation factor for the story drift ratio of the i^{th} story, q_{i}, as _{ } . Note that PF[2] 1 denotes the first mode participation factor for the maximum peak interstory drift angle, q_{max}, obtained using the SRSS of the first two modes. Similarly, PF[1] 1 represents the firstmode participation factor for q_{max} when only the first mode is considered. The coefficient _{ } reflects the effect of the second mode on q_{max} and the relevant spectral shape. The coefficient _{ } mainly indicates the influence of the firstmode and firsttwomodes estimations on q_{max}, which may occur on different stories. The earlier proposed IM_{1E&2E} has also been improved by the same researchers (Luco and Cornell 2007), where the inelastic spectral displacement SI d(T_{1},z_{1},d_{y}) was adopted as a replacement for the elastic spectral displacement S_{d}(T_{1},z_{1}). The new IM is called IM_{1I&2E.} Despite the advancement of IM_{1E&2E} and IM_{1I&2E}, their complex expressions render the various parameters difficult to be determined, especially at the early design stage, thus limiting practical application of these IMs. In addition to the aforementioned scalar IMs, some vectorvalued IMs have also been proposed in recent years. For example, Baker and Cornell (2005) proposed a twoparameter vectorvalued IM called <S_{a}(T_{1}), e>. The definition of parameter e is originated from seismology engineering. This parameter reflects not only the distinction of the spectral shapes but also, to some extent, the effect of highermode contribution and period elongation. However, due to the absence of the basic data required to obtain e in many countries, more work must be performed before the vectorvalued IM <S_{a}(T_{1}), e> can be applied worldwide. Ye et al. (2013) compared the applicability of 18 commonly used IMs in an elastoplastic analysis of different buildings. The findings indicated that S_{a}(T_{1}) has the best correlation with the structural demand measures (DMs) for short and moderateperiod structures. In addition, the correlation between PGV and the DMs is also acceptable. Furthermore, PGV performs better as an IM when the fundamental period of the structure increases. Based on their study, S_{a}(T_{1}) and PGV have been recommended for use as the IMs in actual structural design. Existing research on IMs as outlined above has mainly targeted structures whose fundamental periods are shorter than 6 s and mostly in the range of 0 to 4 s. It is evident that the rationality of the existing IMs requires further validation for super highrise buildings with fundamental periods as long as approximately 10 s. This validation is the focus of the present study, with an emphasis on the actual requirements of structural design. A total of six different IMs (PGA, PGV, PGD, S_{a}(T_{1}), S_{d}(T_{1}), and IM_{1E&2E}) are evaluated and their rationalities are discussed. PGV and S_{a}(T_{1}) were selected based on the recommendations of Ye et al. (2013). To consider the contribution of the higherorder vibration modes of super highrise buildings, IM_{1E&2E }was also selected. In addition, PGA, PGD and S_{d}(T_{1}) were considered because they are widely used in actual design. The correlation between the IMs and DMs with different vibration periods is analyzed, and rational IMs for super highrise building design are recommended. 3. ANALYTICAL MODEL AND GROUND MOTION SELECTION 3.1 Simplified analytical model This study mainly focuses on super highrise buildings with very long fundamental vibration periods. In general, the commonly used lateral seismicresistant super highrise systems include framed tube, bundled tube, tubeintube, diagonalized (trussed tubes, diagrids or braced frames), core and outrigger (central lateral system linking to the perimeter system through outriggers) and hybrid systems (Lu et al. 2011). The deformation modes of these structural systems are typically a combination of flexure and shear, which differs greatly from the traditional frame or shear wall buildings. This implies that the traditional multistoried shear model and cantilever flexural beam model are not suitable for super highrise buildings. Consequently, the dynamic properties and seismic responses of super highrise buildings are approximated herein by using an equivalent continuum model consisting of a shear and a flexural cantilever beam, which is proposed by Miranda and Taghavi (2005) (Figure 1). The shear and flexural cantilever beams are connected by an infinite number of rigid links to achieve horizontal deformation compatibility at the same height. The response of this model subjected to ground acceleration can be expressed in the following partial differential equation: _{ } , (5) where u(x,t) is the horizontal displacement at nondimensional height x at time t; m(x) is the mass per unit length of the model; x varies between 0 at the base of the building and 1 at the roof level; u_{g}(t) represents ground displacement at time t; H denotes the total height of the model; and EI(x) and GA(x) are, respectively, the flexural stiffness of the flexural beam and the shear stiffness of the shear beam at nondimensional height x. Miranda and Taghavi (2005) assumed that the variations in the shear and flexural stiffness along the height are identical, i.e., S(x) or
in which, S(x) is the dimensionless function that defines the variation of stiffness along the height of the simplified model, and EI_{0} and GA_{0 }are, respectively, the flexural and shear stiffness at the base of the model. The dimensionless parameter a_{0}=H(GA_{0}/EI_{0})^{1/2 }controls the proportions of the flexural and shear deformations in the simplified model and thus the lateral deformation shape of the building. Theoretically, when a_{0}=0, the model degenerates into a pure bending model, and when a_{0}=¡Þ, it degenerates into a pure shear model. Miranda and Taghavi (2005) indicated that a_{0} is normally between 0 and 1.5 for shear wall buildings; between 1.5 to 5 for buildings with dual structural systems consisting of a combination of momentresisting frames and shear walls; and between 5 and 20 for momentresisting frame buildings. In addition, when a_{0}¡30, the deformation mode is dominated by shear deformations. A number of different buildings were analyzed by Miranda and Taghavi (2005) to verify the appropriateness of the model in simulating the elastic vibration of different buildings. To ensure the appropriateness of the simplified model, the 632 m super highrise building Shanghai Tower (see Figure 2) is selected herein to calibrate the key parameters used in the model. The Shanghai Tower has a hybrid lateral seismicresistant system called the ¡°megacolumn/coretube/outriggertruss¡±. According to architecture and functional requirements, eight reinforced stories were built for every 13 to 16 stories, which divide the entire structure into eight zones. The plane layouts of these eight stories are also briefly illustrated in Figure 2. The fundamental period of the Shanghai Tower is nearly 9.83 s (Lu et al. 2011), which is far beyond the range of 6 s specified in the design response spectrum of the Chinese Code for the Seismic Design of Buildings (MOC 2010). The parameters of the simplified model are calibrated such that the fundamental dynamic properties of the simplified model approach those of the Shanghai Tower in the XZ plane. The total gravitational load of the simplified model, which is determined from a combination of 1.0 D + 0.5 L, is identical to that of the Shanghai Tower, in which D is the dead loads, including the effects of self weight and permanently attached equipment and fixtures; L is the live loads. Note that the dead and live loads of the Shanghai Tower have little variation along its height; the simplified model thus has a constant mass distribution along the height. According to the stiffness distribution of the Shanghai Tower, a quadratic curve, i.e., S(x)=1(10.55)x^{2}, is adopted to model the structural stiffness along the height (Eqns 6 and 7). The flexural stiffness at the base, EI_{0}, and the dimensionless flexural and shear ratio, a_{0}, are calibrated according to the dynamic properties of the Shanghai Tower. The values of the main parameters used in the simplified model are summarized in Table 1. The first five free translational vibration periods and the relative deviation D between the simplified model and the fine FE model of Shanghai Tower are compared in Table 2. The relative deviation D, defined by Eqn 8, is adopted to describe the distinction between the two models. In the equation, T_{simplify} and T_{fine} represent, respectively, the free translational vibration periods of the simplified model and the fine FE model of Shanghai Tower.
It is evident that some discrepancies are present in the comparison. Note that a substantial reduction in computational costs is achieved by employing the simplified model compared with a very fine FE model. Note also that the focus of this study is to offer a general rule for super highrise buildings instead of investigating the characteristics of a particular structure. Therefore, the resulting discrepancies are considered acceptable. 
3.2 Ground motion selection Due to significant randomness in the ground motion records, the structural responses may vary significantly in a timehistory analysis. To minimize such a variation while maintaining inherent randomness of the ground motions, 22 pairs of farfield ground motion records recommended by the FEMA P695 (FEMA 2009) are adopted herein as the basis ground motion record set. Note that two horizontal components of each ground motion are considered, thereby resulting in 44 natural records. Their corresponding elastic acceleration spectrums at a 5% damping ratio are presented in Figure 3. This ground motion record set not only maintains the randomness of the ground motions themselves but also leads to a certain degree of consistency in the structural responses. 4. EVALUATION OF THE INTENSITY MEASURES 4.1 Evaluation method The existing literature (Baker and Cornell 2005, 2008; Padgett et al. 2008) indicates that the relationship between the structural seismic response demand measures (DMs) and the ground motion intensity measures (IMs) can be estimated using the powerfunction model presented in Eqn 9: _{ } , (9) in which a and b are the regression coefficients. The equation can be rearranged in a linear regression form of logarithms of DM and IM; i.e., _{ } (10) Eqn 10 clearly satisfies the classical linear regression model. The correlation coefficient r between ln(DM) and ln(IM) and the corresponding dispersion b can then be obtained by regression analysis of n discrete data points (DM_{i}, IM_{i}). This analysis can be performed using the least square method in timehistory analysis. The value of dispersion b is calculated from Eqn 11 below: _{ } (11) The correlation coefficient r ranges between 1 and 1. The closer to 1 r is, the better the correlation between the calculated DMs and IMs. In general, r¡0.8 implies a good correlation. In addition, the smaller the dispersion b, the better the IMs are. Note that all modernday super highrise buildings are rigorously designed with sufficient safety margins to meet the seismic requirement under the Maximum Considered Earthquake (MCE) specified in the design codes. Thus, the major lateral force resisting components remain elastic or slight yielding, instead of exhibiting obvious nonlinear behavior, when subjected to MCE. Therefore, elastic time history analysis is considered appropriate in evaluating the rationality of the existing IMs for super highrise buildings. The specific research procedures are as follows: (1) For a simplified model with a given fundamental period and corresponding parameters, calculate the maximum structural response DM_{i} subjected to the i^{th} ground motion record using the timehistory analysis method; (2) Calculate the IM_{i} of the corresponding i^{th} ground motion record; (3) Calculate n ground motion records and obtain n discrete data points (IM_{i}£¬DM_{i}) and plot them in the ln(IM)ln(DM) coordinate system. The correlation coefficient r between ln(IM) and ln(DM) for the given fundamental period and the corresponding dispersion b are obtained by linear regression. The schematic diagram of correlation between the IM and DM is shown in Figure 4; (4) Adjust the model parameters (EI_{0} or r(x)) to obtain a series of structures with different fundamental periods. Next, repeat steps 1 to 3 to calculate the variations in the correlation coefficient between DM_{i} and IM_{i} and in the corresponding dispersion for different structural fundamental periods. Note that although only EI_{0} or r(x) is modified in this step, the influence of adjusting other parameters will be discussed in later sections of this work. While many DMs are available, such as the maximum top displacement (d_{max}), maximum acceleration (a_{max}), maximum story drift ratio (q_{max}), maximum base shear force (F_{max}), and total input energy (E_{input}), the most commonly used DMs in structural seismic design and seismic response analysis are q_{max}, d_{max} and a_{max}. Therefore, this paper aims to discuss the variations in the correlation coefficients between these DM_{i} (q_{max}, d_{max}, a_{max}) and IM_{i} and their corresponding dispersion with different structural fundamental periods T_{1}. 4.2 Result Analysis Based on the simplified model with a_{0}=4 proposed in the above section in Table 1, five different simplified models with the fundamental periods of 1, 3, 6, 9 and 10 s respectively are established by adjusting the values of EI_{0} and r(x). In these models, classical Rayleigh damping is adopted with a value of 5%. Figures 5 to 10 demonstrate the variations in correlation between the DMs of interest (q_{max}, d_{max}, a_{max}) and the IMs (PGA, PGV, PGD, S_{a}(T_{1}), S_{d}(T_{1}), and IM_{1E&2E}) and their corresponding dispersion for different fundamental periods. Figure 5 demonstrates that the correlation between PGA, S_{a}(T_{1}), S_{d}(T_{1}), and IM_{1E&2E }and q_{max}_{ }decrease gradually with an increase at the fundamental period, particularly for PGA. When the fundamental period approaches 10 s, the correlation coefficient of PGA is less than 0.4, which indicates that PGA has little correlation with q_{max}_{ }for the structures with super long periods. In addition, the correlation coefficient for IM_{1E&2E} is approximately 0.7, which indicates that IM_{1E&2E} has some correlation with q_{max}; however, it does not accurately represent the characteristic of the seismic response of super highrise buildings because only the first two modes are considered. The correlation between PGV and PGD and q_{max} increase gradually with increasing fundamental period, and after 6 s, the correlation between PGV and PGD and q_{max} decreases slightly. In addition, for a fundamental period of 9 s, the correlation coefficient r of PGD increases to 0.6, while the correlation coefficient r of PGV approaches 0.84. Therefore, PGV has the best correlation with q_{max},_{ }and PGA has the least correlation with q_{max }at a period of 10 s. Similarly, Figure 6 indicates that in the range of 10 s, the dispersion of PGV is the smallest, while the dispersion of PGA is the largest. These results indicate that for super highrise buildings, PGV has the best correlation with q_{max }and the minimum dispersion. The existing literature also suggests that q_{max} has a significant correlation with structural damage (Bozorgnia and Bertero 2001; Miranda and Akkar 2006). This further confirmed that PGV can be suggested as an IM for structural seismic design and collapse analysis for super highrise buildings. Figure 7 indicates that S_{a}(T_{1}), S_{d}(T_{1}) and IM_{1E&2E} have the greatest correlation with d_{max}_{,} and the correlation coefficients always remain at approximately 0.94, even when the fundamental period increases from 1 to 10 s. The correlation between PGD and d_{max} increases gradually from less than 0.25 (at a period of 1 s) to 0.94 (at a period of 10 s). The correlation between PGV and d_{max} decreases gradually, and the correlation coefficient r is approximately 0.48 at a period of 10 s. Finally, the correlation between PGA and d_{max }decreases rapidly with increasing period and even becomes negative when the fundamental period is equal to 10 s. Figure 8 also indicates similar results, i.e., S_{a}(T_{1}), S_{d}(T_{1}) and IM_{1E&2E} have the minimum dispersion among these IMs. Figures 9 and 10 demonstrate that PGA has the greatest correlation with the structural maximum acceleration (a_{max}) and the smallest dispersion b; and the value of the correlation coefficient r remains at approximately 0.85 when the fundamental period increases from 1 to 9 s. PGV also shows some correlation with a_{max}_{,}, and the correlation coefficient r remains near 0.6 and increases slightly with an increasing period. In contrast, S_{a}(T_{1}), S_{d}(T_{1}) and IM_{1E&2E} display poor correlations with a_{max }with large dispersions b. The correlation coefficients decrease rapidly with an increase in the fundamental period. When the fundamental period becomes very long, the correlation coefficients become negative. This implies that the two variations have very little correlation, in other words, the more the IMs (S_{a}(T_{1}), S_{d}(T_{1}) and IM_{1E&2E}) are, the smaller the a_{max} is. In addition, PGD also shows little or no correlation with a_{max }and the correlation coefficients is close to 0 in the period range concerned. In performancebased design, acceleration is usually associated with comfort and property loss; therefore, PGA can be selected as the IM for structural acceleration prediction. 4.3 Influence of a_{0} To make the aforementioned results universal and representative, the influence of a_{0 }on the simplified model is discussed. Miranda and Taghavi (2005) reported that when a_{0}=0, the simplified model degenerates into a pure bending model, whereas when a_{0}=¡Þ, this model can be approximately considered as a pure shear model. Although pure bending or pure shear behavior have little possibility of occurrence in super highrise buildings with a vibration period as long as 9 s, for investigation purpose, a series of simplified models with a_{0}=0 and 30 and T_{1}=1, 3, 6, 9, and 10 s (totally 10 simplified models) are created to examine the relationship between a_{0 }and the correlation coefficient between DMs and IMs. The variation in the correlation coefficients and the corresponding dispersion between q_{max} and each IM when a_{0}=0 are illustrated in Figures 11 and 12, respectively. In addition, similar variations in the correlation and dispersion when a_{0}=30 are illustrated in Figures 13 and 14, respectively. For different values of a_{0, }Figures 15 and 16 demonstrate, respectively, the correlation coefficient r and the corresponding dispersion b between PGV and q_{max}. When a_{0}=0, the simplified model represents a pure bending structure, whose higherorder period decreases rapidly. The second and third vibration periods (i.e., T_{2} and T_{3}) are approximately 1/6 and 1/16 of the fundamental period T_{1}, respectively. For the flexural/shear structures with a_{0}=4, on the other hand, T_{2} and T_{3} are approximately 1/3 and 1/6 of T_{1}, respectively. It is evident that the influence of higherorder modes is smaller for pure bending structures. However, IM_{1E&2E}, which considers the first two vibration modes, has a significant correlation with q_{max}. Figures 11 and 12 indicate that the correlation coefficients decrease slightly with increasing T_{1}, and the values of the correlation coefficients remain above 0.9 for different fundamental periods. The dispersion of IM_{1E&2E} is also the smallest. The correlation coefficients between PGV and q_{max} at different fundamental periods when a_{0}=0 are always above 0.7. The correlation between S_{a}(T_{1}) (or S_{d}(T_{1})) and q_{max} when a_{0}=0 decreases gradually when T_{1} increases. When T_{1} is shorter than 6 s, the correlation between S_{a}(T_{1}) (or S_{d}(T_{1})) and q_{max} is slightly greater than that of PGV; however, when T_{1} is longer than 6 s, the correlations are slightly worse than those of PGV. The correlation between PGD and q_{max }increases rapidly with increasing T_{1}, and when T_{1} is longer than 6 s, the correlation decreases slightly, which is similar to the variation in PGV. PGA exhibits the least correlation with q_{max}, and the correlation coefficients decrease rapidly with increasing T_{1}. When T_{1} approaches 10 s, the correlation coefficient increases slightly, but the value is still very low at approximately 0.27. When a_{0}= 30, the variations in the correlation coefficient r and the corresponding dispersion b between q_{max} and each IM at different fundamental periods approximate the values when a_{0}=4. That is, by increasing the structural fundamental period, the correlations between the IMs (PGA, S_{a}(T_{1}), S_{d}(T_{1}), and IM_{1E&2E}) and q_{max}_{ }decrease gradually, while the correlations between PGV (or PGD) and q_{max}_{ }increase steadily. For a super longperiod structure with a fundamental period longer than 9 s, PGV exhibits the best correlation with q_{max}_{ }and minimum dispersion. In contrast, PGA exhibits the least correlation with q_{max }and maximum dispersion. Furthermore, Figures 15 and 16 demonstrate that for a pure bending structure (a_{0}=0), the correlation coefficients between PGV and q_{max} at different fundamental periods are smaller than the values of a flexuralshear deformation combined structure (a_{0}=4) or a shear deformation dominant structure (a_{0}=30). For the bendingshear deformation combined structure (a_{0}=4), the correlation coefficients are slightly smaller than the values for the shear deformation dominant structure (a_{0}=30) if the fundamental period T_{1} is shorter than 5 s. However, when T_{1} is longer than 5 s, the correlation coefficients are slightly larger than the values for the shear deformation dominant structure (a_{0}=30). Moreover, when T_{1} approaches 10 s, these values are very similar. Overall, for different values of a_{0}, PGV always shows a good correlation with q_{max}. Considering that in practical application, most super highrise buildings are flexuralshear deformation combined systems, it is still feasible to select PGV as the IM for the structural seismic design and collapse analysis of super highrise buildings. Similarly, for different values of a_{0}, Figures 17 and 18 demonstrate, respectively, the variations of the correlation coefficient r and the corresponding dispersion b between S_{a}(T_{1}) and d_{max}; and such variations between PGA and a_{max} are presented in Figures 19 and 20, respectively. Figure 17 illustrates that a_{0} has little influence on r between S_{a}(T_{1}) and d_{max}, where r is constantly above 0.95 for different values of a_{0} in the period range concerned. The corresponding b remains below 0.2 (see Figure 18). Figure 19 illustrates that best correlations are achieved between PGA and a_{max} for shear deformation dominant structures (a_{0}=30), and this is followed by flexuralshear deformation combined structures (a_{0}=4). Whereas smallest correlation is found for pure bending structures (a_{0}=0). Such findings suggest that the larger the value of a_{0 }is, the better the correlation is between PGA and a_{max}. It is necessary to note that the smallest correlation coefficient for pure bending structures (a_{0}=0) is still larger than 0.78. This implies that PGA correlates very well with a_{max} when a_{0} varies between 0 and 30; and the influence of a_{0} on the correlation coefficients is insignificant. In general, a_{0 }has little influence on the correlation between IMs and DMs. Although the discussions given in Section 4.2 are based on the analysis of the Shanghai Tower, the findings remains universal and are applicable to other similar super highrise buildings. In summary, for different values of a_{0}, PGV, S_{a}(T_{1}) and PGA always demonstrate good correlations with q_{max}, d_{max} and a_{max}, respectively. 
5. VERIFICATION THROUGH COLLAPSE ANALYSIS OF THE SHANGHAI TOWER It should be noted that the maximum story drift ratio q_{max} is always the controlling factor for seismic safety measure of building structures. Due to a very good correlation between PGV and q_{max}, as indicated in the above analysis, PGV could be selected as an IM for seismic design and response prediction of super highrise buildings. This conclusion is reached based on the simplified model illustrated in Figure 1 with elastic dynamic analysis; however, a super highrise building may enter the nonlinear stage or even collapse under extreme earthquake conditions. Hence, the rationality of this conclusion must be further validated when the super highrise building enters the nonlinear region. For this reason, the following section presents a verification of the above conclusion using an incremental dynamic collapse analysis of the Shanghai Tower presented in Figure 2. Lu et al. (2011, 2013a, 2013b) proposed a nonlinear finite element model based on the general purpose finite element code of MSC.Marc to simulate the collapse of the highrise and super highrise buildings, including the Shanghai Tower. In the FE model of Shanghai Tower, the fiber model is adopted to simulate the columns and beams in the external frame and the components of the outriggers. The multilayer shell model is used to simulate the shear walls, the coupling beams and the mega columns; and elemental deactivation technology is adopted to simulate the failure of the components upon reaching their ultimate deformation stage. The complete FE model of the Shanghai Tower is illustrated in Figure 2. Further details of the structural dimensions, the material constitutive relationships and the collapse simulation procedures can be found in Lu et al. (2011). Incremental dynamic analysis (IDA) is a powerful method in earthquake engineering for performing a comprehensive assessment of the structural behavior under seismic loads. During IDA, several ground motions are considered with increased levels of seismic intensity thereby making the structure undergo an entire range of behavior, from elastic to inelastic and finally to collapse. In this paper, the larger components of each of the 22 pairs of farfield ground motion records suggested by FEMA P695 (FEMA 2009) are adopted as the ground motion set. In addition, data from the ElCentro EW 1940 ground motion record and the Shanghai artificial ground motion record are added into this ground motion set. This makes a total of 23 natural records and 1 artificial record. During the analysis, these ground motion records are input into the fine FE model along the X direction and their intensities are gradually increased up to collapse of the building. One of the typical potential collapse modes is shown in Figure 21. The values of different IMs (PGA, PGV, PGD, S_{a}(T_{1}), S_{d}(T_{1}), and IM_{1E&2E}) that cause the building to collapse are presented in Table 3. Since this building is very important, its safety margin against earthquakeinduced collapse is very large. The coefficient of variation (COV) of each IM is presented in Figure 22. It can be seen that when the building begins to collapse under the 24 ground motion records, the COVs of S_{d}(T_{1}) and S_{a}(T_{1}) are maximum, reaching nearly 0.729 and 0.690. This is followed by the PGA, PGD and IM_{1E&2E}. In addition, PGV has the minimum COV, i.e. 0.182, which is approximately 2/5 of the COV of PGA and 1/4 of the COV of S_{a}(T_{1}). Even when the super highrise building enters the nonlinear stage, PGV still shows a very good correlation with the structural response and collapse possibility. PGV can therefore be selected as the IM for seismic response prediction and collapse analysis. 
6. CONCLUSIONS This paper briefly reviewed some existing IMs and constructed a simplified flexuralshear beamcoupling model based on an actual super highrise building, the Shanghai Tower. The correlations between the fundamental DMs and some existing IMs of this super highrise building are analyzed, and the appropriate scopes for different IMs are also discussed. Consequently, the following preliminary conclusions are obtained: (1) With increasing structural fundamental period, the correlations between PGA, S_{a}(T_{1}), S_{d}(T_{1}) and the maximum story drift ratio (q_{max}) decrease rapidly, and their corresponding dispersions increase gradually. Hence, PGA, S_{a}(T_{1}) and S_{d}(T_{1}) are not suitable IMs for super highrise buildings for seismic design and time history analysis. (2) In the range of super long periods, the correlation coefficient between PGV and q_{max} is more than 0.8 with a smallest dispersion. These results indicate that PGV has a significant correlation with q_{max}. Therefore, PGV can be adopted as the IM for super highrise buildings with super long periods. (3) Over the range of short to super long periods, S_{a}(T_{1}) and S_{d}(T_{1}) always exhibit significant correlations with the structural maximum top displacement (d_{max}); therefore, S_{a}(T_{1}) and S_{d}(T_{1}) can be selected as the IMs for structural global displacement prediction. (4) PGA exhibits a significant correlation with the structural acceleration response; therefore, PGA can be selected as the IM for structural acceleration prediction and comfort design. (5) From the collapse analysis of the Shanghai Tower, it can be concluded that under strong or extreme earthquakes, even when the super highrise building enters the nonlinear region, PGV can still be selected as the IM for seismic response prediction and collapse analysis. It is worth mentioning that the above findings are reached based on the investigation of a single building ¨C the Shanghai Tower. To ascertain these research findings to be used for seismic design guidelines of super highrise buildings, more case studies on other tall buildings are desirable in further work. ACKNOWLEDGMENT The authors are grateful for the financial support received from the National Nature Science Foundation of China (No. 51222804, 51261120377, 51178249), the Tsinghua University Initiative Scientific Research Program (No. 2012THZ022, 2011THZ03) and the Fok Ying Dong Education Foundation (No. 131071). References Baker, J.W. and Cornell, C.A. (2005). ¡°A vectoredvalued ground motion intensity measure consisting of spectral acceleration and epsilon¡±, Earthquake Engineering and Structural Dynamics, Vol. 34, No. 10, pp. 11931217. Baker, J.W. and Cornell, C.A. (2008). ¡°Vectorvalued intensity measures for pulselike nearfault ground motions¡±, Engineering Structures, Vol. 30, No. 4, pp. 10481057. Bozorgnia, Y. and Bertero, V.V. (2001). ¡°Improved shaking and damage parameters for postearthquake applications¡±, Proc., SMIP01 Seminar on Utilization of StrongMotion Data, California Division of Mines and Geology, Los Angeles, pp. 122. Cordova, P.P., Deierlein, G.G., Mehanny, S.S.F. and Cornell, C.A. (2001). ¡°Development of a twoparameter seismic intensity measure and probabilistic assessment procedure¡±, The Second U.S.Japan Workshop on Performancebased Earthquake Engineering Methodology for Reinforced Concrete Building Structures, Sapporo, Japan, pp. 187206. FEMA (2009). FEMA P695: Quantification of building seismic performance factors, Applied Technology Council, Redwood City, C.A. Lu, X., Lu, X.Z., Zhang, W.K. and Ye, L.P. (2011). ¡°Collapse simulation of a super highrise building subjected to extremely strong earthquakes¡±, Science ChinaTechnological Sciences, Vol. 54, No. 10, pp. 25492560. Lu, X., Lu, X.Z., Guan, H. and Ye, L.P. (2013a). ¡°Collapse simulation of reinforced concrete highrise building induced by extreme earthquakes¡±, Earthquake Engineering and Structural Dynamics, Vol. 42, No.5, pp. 705723. Lu, X.Z., Lu, X., Guan, H., Zhang, W.K. and Ye, L.P. (2013). ¡°Earthquakeinduced collapse simulation of a supertall megabraced framecore tube building¡±, Journal of Constructional Steel Research, Vol. 82, pp. 5971. Lucchini, A., Mollaioli, F. and Monti, G. (2011). ¡°Intensity measures for response prediction of a torsional building subjected to bidirectional earthquake ground motion¡±, Bulletin of Earthquake Engineering, Vol. 9, No.5, pp. 14991518. Luco, N. and Cornell, C.A. (2007). ¡°Structurespecific scalar intensity measures for nearsource and ordinary earthquake ground motions¡±, Earthquake Spectra, Vol. 23, No. 2, pp. 357391. MOC (2010). Code for seismic design of buildings GB500112010, Ministry of Construction of the People¡¯s Republic of China, Beijing, China. Miranda, E. and Akkar, S.D. (2006). ¡°Generalized interstory drift spectrum¡±, Journal of Structural Engineering, ASCE, Vol. 132, No. 6, pp. 840852. Miranda, E. and Taghavi, S. (2005). ¡°Approximate floor acceleration demands in multistory buildings I: Formulation¡±, Journal of Structural Engineering, ASCE, Vol. 131, No. 2, pp. 203211. Padgett, J.E., Nielson, B.G. and DesRoches, R. (2008). ¡°Selection of optimal intensity measures in probabilistic seismic demand models of highway bridge portfolios¡±, Earthquake Engineering and Structural Dynamics, Vol. 37, No. 5, pp. 711725. Shome, N., Cornell, C.A., Bazzurro, P. and Carballo, J.E. (1998). ¡°Earthquakes, records and nonlinear responses¡±, Earthquake Spectra, Vol. 14, No. 3, pp. 469500. The Building Centre of Japan. (2000). ¡°The building standard law of Japan¡±, Japan. Tothong, P. and Luco, N. (2007). ¡°Probabilistic seismic demand analysis using advanced ground motion intensity measures¡±, Earthquake Engineering and Structural Dynamics, Vol. 36, No. 13, pp. 18371860. Vamvatsikos, D. and Cornell, C.A. (2005). ¡°Developing efficient scalar and vector intensity measures for IDA capacity estimation by incorporating elastic spectral shape information¡±, Earthquake Engineering and Structural Dynamics, Vol. 34, No.13, pp. 15731600. Yang, X.Q., Fu ,X.Y. and Huang, Y.J. (2011). ¡°Dynamic elastoplastic analysis of the Shenzhen Ping¡¯an Financial Center Tower¡± (in Chinese). Journal of Building Structures, Vol. 32, No.7, pp.4049. Ye, L.P., Ma, Q.L., Miao, Z.W., Guan, H. and Zhuge, Y. (2013). ¡°Numerical and comparative study of earthquake intensity indices in seismic analysis¡±, The Structural Design of Tall and Special Buildings, Vol. 22, No. 4, pp.362381. List of Tables Table 1. Main parameters used in the simplified model Table 2. Comparison of the first five periods between the simplified model and the fine FE model Table 3. The values of different IMs at the critical collapse status of the Shanghai Tower List of figures Figure 1. The simplified model consisting a flexural and a shear beam Figure 2. The FE model of Shanghai Tower Figure 3. The elastic acceleration spectrum of the 44 farfield ground motions at a 5% damping ratio Figure 4. The schematic diagram of correlation between the IM and DM Figure 5. The variations of correlation between q_{max }and_{ }IMs at different fundamental periods Figure 6. The variations of dispersion between q_{max }and_{ }IMs at different fundamental periods Figure 7. The variations of correlation between d_{max} and IMs at different fundamental periods Figure 8. The variations of dispersion between d_{max} and IMs at different fundamental periods Figure 9. The variations of correlation between a_{max} and IMs at different fundamental periods Figure 10. The variations of dispersion between a_{max} and IMs at different fundamental periods Figure 11. The variations of correlation between q_{max }and_{ }IMs at different fundamental periods when a_{0}=0 Figure 12. The variations of dispersion between q_{max }and_{ }IMs at different fundamental periods when a_{0}=0 Figure 13. The variations of correlation between q_{max }and_{ }IMs at different fundamental periods when a_{0}=30 Figure 14. The variations of dispersion between q_{max }and_{ }IMs at different fundamental periods when a_{0}=30 Figure 15. The variation of correlation coefficients between PGV and q_{max} with different value of a_{0} Figure 16. The variation of dispersion between PGV and q_{max} with different value of a_{0} Figure 18. The variation of dispersion between S_{a}(T_{1}) and d_{max} with different value of a_{0} Figure 19. The variation of correlation coefficients between PGA and a_{max} with different value of a_{0} Figure 20. The variation of dispersion between PGA and a_{max} with different value of a_{0} Figure 21. The potential collapse mode of the Shanghai Tower (The failed elements are removed from the model) Figure 22. The distribution of COV for different IMs at the stage of collapse Appendix. Notation The following symbols are used in this paper: IM ground motion intensity measure DM seismic response demand measure PGA the peak ground motion acceleration PGV the peak ground motion velocity PGD the peak ground motion displacement S_{a}(T_{1}) the spectral acceleration at fundamental period S_{d}(T_{1}) the spectral displacement at fundamental period d_{max} the maximum top displacement a_{max } the maximum acceleration q_{max} the maximum story drift ratio IDA incremental dynamic analyses THA timehistory analysis a_{0} the dimensionless parameter controls the proportions of the flexural and shears deformations r the correlation coefficient between DM and IM b the dispersion coefficient between DM and IM COV the coefficient of variation Table 1. Main parameters used in the simplified model
Table 2. Comparison of the first five periods between the simplified model and the fine FE model
Table 3. The values of different IMs at the critical collapse status of the Shanghai Tower
Figure 1. The simplified model consisting a flexural and a shear beam
Figure 2. The FE model of Shanghai Tower
Figure 3. The elastic acceleration spectrum of the 44 farfield ground motions at a 5% damping ratio
Figure 4. The schematic diagram of correlation between the IM and DM
Figure 5. The variations of correlation between q_{max }and_{ }IMs at different fundamental periods
Figure 6. The variations of dispersion between q_{max }and_{ }IMs at different fundamental periods
Figure 7. The variations of correlation between d_{max} and IMs at different fundamental periods
Figure 8. The variations of dispersion between d_{max} and IMs at different fundamental periods
Figure 9. The variations of correlation between a_{max} and IMs at different fundamental periods
Figure 10. The variations of dispersion between a_{max} and IMs at different fundamental periods
Figure 11. The variations of correlation between q_{max }and_{ }IMs at different fundamental periods when a_{0}=0
Figure 12. The variations of dispersion between q_{max }and_{ }IMs at different fundamental periods when a_{0}=0
Figure 13. The variations of correlation between q_{max }and_{ }IMs at different fundamental periods when a_{0}=30
Figure 14. The variations of dispersion between q_{max }and_{ }IMs at different fundamental periods when a_{0}=30
Figure 15. The variation of correlation coefficients between PGV and q_{max} with different value of a_{0} _{ }_{} Figure 16. The variation of dispersion between PGV and q_{max} with different value of a_{0} _{ }
Figure 17. The variation of correlation coefficients between S_{a}(T_{1}) and d_{max} with different value of a_{0} _{ }
Figure 18. The variation of dispersion between S_{a}(T_{1}) and d_{max} with different value of a_{0} _{ }
Figure 19. The variation of correlation coefficients between PGA and a_{max} with different value of a_{0} _{ }
Figure 20. The variation of dispersion between PGA and a_{max} with different value of a_{0} _{ }_{} Figure 21. The potential collapse mode of Shanghai Tower (The failed elements are removed from the model)
Figure 22. The distribution of COV for different IMs at the stage of collapse [*] Corresponding author. Email address: luxz@tsinghua.edu.cn. Tel: +861062795364; fax: +861062795364 