Abstract

In recent studies, spectral matching is the most commonly proposed method for selecting earthquake records for time-history analysis of structures. However, until now, there have been no serious investigations of the effects of coordinate values on the scaling of ground motions. This paper investigated the influence of using arithmetic and logarithmic values of response spectra in spectral matching procedures (i.e., ASM and LSM methods) on the results of nonlinear structural time-history analysis. Steel moment resisting frame structures of the 3-, 9-, and 20-stories, which represent low-, medium-, and high-rise buildings, respectively, were used as examples. Structural benchmark responses were determined by calculating the arithmetic mean and median of peak interstory drift ratio (PIDR) demands based on the three record sets developed by the American SAC Steel Project. The three record sets represent seismic hazard levels with 50%, 10%, and 2% probabilities exceeded in 50 years, and their average acceleration spectra were also taken as the target spectrum. Moreover, another 40 record components for selection were scaled both by ASM and LSM methods. The seven components whose spectra were best compatible with the target spectra were selected for the structural time-history analysis. The scale factors obtained by the LSM method are nearly larger than that of the ASM method, and their ranking and selection of records are different. The estimation accuracies of structural mean (median) responses by both methods can be controlled within an engineering acceptable range (±20%), but the LSM method may cause larger structural responses than the ASM method. The LSM method has a better capacity for reducing the variability of structural responses than the ASM method, and this advantage is more significant for longer-period structures (e.g., 20-story structure) with more severe nonlinear responses.

1. Introduction

The nonlinear time-history analysis has been a widely accepted method for seismic design and assessment of structures, also being the most prevalent process used for performance-based design [14]. According to the seismic design codes (e.g., Eurocode 8 [5] and ASCE 7-16 [6]), time-history analyses are typically conducted with 7 or 11 ground motions records, and the arithmetic mean of the engineering demands parameters (EDPs) (e.g., story drifts and rotation capacity of plastic hinge) are taken as the reference results for the design [2, 7, 8]. However, the results of time-history analysis are significantly influenced by the ground motions input. In order to reduce record-to-record variability and obtain a more accurate estimate of structural responses, a few methods for selection and scaling of ground motions have been developed [2, 4, 817]. The nonlinear time-history analysis is routinely adopted for the seismic design of base-isolated buildings, and some more advanced methods have been proposed, such as the weighted scaling methods [1820] or the selection of standard records and synthetic accelerograms scaled by the PGV [21]. Among these methods, more attention has been given to amplitude linear scaling for a real record and its response spectral matching with the target spectrum [813, 15, 1824].

In the ground motion scaling and spectral matching procedure, the spectral compatibility between the record’s response spectrum and the target spectrum is mainly verified by the sum of the squared error (SSE) [8, 11, 1820] or the general root mean square error [9, 10, 15, 24]. Error indexes should be minimized so that the record can realize the optimal spectral matching. Either the arithmetic values or the logarithm values of response spectra were adopted to calculate the error indexes in several different published studies. Arithmetic values were used in references [810, 1315, 1820, 24] since engineers are more familiar with arithmetic mean and design spectra in seismic codes can be taken directly as target spectra as well. Meanwhile, some researchers [11, 12, 22, 23] used the logarithm values of response spectra because logarithmic values are consistent with the assumption that the structural seismic responses are log-normally distributed. Logarithmic values are commonly used in the probabilistic seismic hazard analysis (PSHA) research field. The target spectra derived from the PSHA of a given site, such as the uniform spectrum (UHS) [25, 26] and the conditional mean spectrum (CMS) [11], were also fit for using the logarithm values. Until now, there have been no serious investigations regarding the effects of coordinate values on the scaling of the ground motions. It has not been ascertained whether the scale factors (SFs) for a record are the same under the two coordinates systems. Furthermore, for the same target spectra, it has not been ascertained if the spectral matching procedure is performed using the arithmetic values and logarithmic values and whether the structural EDPs predicted by the two procedures are similar.

In this paper, the 3-, 9-, and 20-story steel moment resisting frame structures, which represent low-, medium-, and high-rise buildings, were taken as examples. Their seismic responses were evaluated by time-history analysis. The arithmetic values and logarithmic values of response spectra were both used in spectral matching with the target spectra and then the ground motions were selected and scaled for structural time-history analysis. The target spectra that matched with the NEHRP 1994 design spectra [27] were defined as the average spectral accelerations of the three ground motion sets. The three record sets were developed in the American SAC Steel Project and represent 50%, 10%, and 2% probabilities of being exceeded in 50 years. The EDPs predicted by the two coordinate systems were compared with each other and influence factors, such as the structural periods and nonlinear response degrees were also discussed.

2. Spectral Matching with Target Spectra

When the target spectrum is given, the sum of squared errors can be used to measure the compatibility between the spectrum of the selected ground motion and the target spectrum. When using arithmetic values, the error index, SSEA, is calculated by the following equation [8]:where and are the target spectral acceleration and the record’s (unscaled) spectral acceleration at period , respectively, is the scale factor, is the total number of discrete period points and should be greater than 50 over the entire range of the periods (i.e., [, 1.5T1], T1 is the first mode period of the structure) [11], and is the number of vibration modes that can ensure that the cumulative modal mass participation ratio is greater than 90%, which considers the effects of the higher modes for the structural responses [13]. The upper bound is 2T1 according to ASCE 7-16 [6], indicating that the lowest useable frequency of the selected records has to be extremely low if the fundamental period of the structure is longer (e.g., T1 = 4.11 s for the 20-story structure in this study). It will decrease the number of available records that can be used for selecting from the PEER NGA database. Hence, the upper bound is adopted as 1.5T1 in this study according to ASCE 7-10 [28].

To obtain a minimum for each record, the between the record’s scaled spectrum and the target spectrum are minimized using the least squares method; i.e., , and the is determined by the following equation:

When the logarithmic values are considered, the sum of squared errors () and the scale factors () are determined by the following equations [11]:

In addition, 0.02 and 0.04 were both used as the intervals of the period under the arithmetic and logarithmic coordinates systems and the SFs calculated using the both period intervals are all the same. Thus, 0.02 was used as the interval of the period in this study. The dense interval can ensure that the SFs are not influenced by the chosen coordinates systems.

In this paper, the spectral matching methods that use arithmetic values and logarithmic values of response spectra are referred to as the ASM and LSM method, respectively.

3. Definition of the Target Spectrum

Somerville [29] put forward three sets of ground motions in the SAC project that were selected based on the probabilistic seismic hazard analysis (PSHA) and matched with the NEHRP 1994 design response spectra [27]. These records represented 50/50, 10/50, and 2/50 seismic hazard levels in Los Angeles (i.e., 50%, 10%, and 2% probabilities of being exceeded in 50 years at the building site). Each set consisted of 20 records (i.e., 10 ground motions each with two orthogonal components). For simplicity, the three probabilities of exceedance or seismic hazard levels are all denoted as 50/50, 10/50, and 2/50 in the figures and tables in this work. Most of these records are real earthquake records, but there are several simulated ground motions in the record set of the 2/50 hazard level. The average (arithmetic mean) spectrum of each set of records is defined as the target spectrum representing the given hazard level. The target spectra representing the three hazard levels with a damping ratio of are plotted in Figure 1 on the arithmetic scale (Figure 1(a)) and logarithmic scale (Figure 1(b)). The damping ratio is consistent with that of the three steel structures used as analysis models.

4. The Referenced Structures for Analysis

The 3-, 9-, and 20-story steel moment resisting frame structures used as analysis models in this work were initially designed for the SAC Phase II Steel Project [30], assuming that they would be located in Los Angeles. The buildings were not actually constructed, but they represent typical short-, medium-, and high-rise buildings in regions where there is a high risk of earthquakes. Therefore, these buildings were adopted to investigate the effects of the ASM and LSM methods on the results generated by nonlinear structural time-history analysis. More detailed information is provided in references [30, 31].

The finite element models (FEMs) of the buildings were built based on the platform of ABAQUS 6.12 [32] according to the structural details in reference [30]. The moment resisting frames are located at the perimeter of the structure. This design practice resulted in the presence of many interior frames which were designed to suffer from gravity loads only. There is no need to build the interior frames in analysis models. Only the moment resisting frames which resistant the horizontal forces, such as wind loads and earthquake loads, were modeled as two-dimensional frames that represented half of the structure in one direction [30]. B22 elements were used to model the beams and columns. It was assumed that plastic hinges were present only at the end of all structural beams and at the end of the columns at the first floor. The bilinear material constitutive model with the strain-hardening ratio of 0.01 was applied to the steel material. Rayleigh damping with a damping ratio of 2% was used for the first and second vibration modes. The Finite Element Modeling (FEM) took into account the effect of the P-Delta and the shear behavior of the Panel Zone according to reference [31].

Table 1 lists the lower three or four vibration periods and the modal-mass participation factors for the three structures. For the 3-, 9-, and 20-story structures, the lower two, four, and three modes should be considered in the selection and scaling of ground motions, respectively.

5. Selection and Scaling of Ground Motions

5.1. Ground Motions for Selection

We only used real records for selection and scaling in this study. Based on preliminary selection criteria [33],including the earthquake magnitude (MW ≥ 6) and the fault rupture distance (20 km ∼ 40 km), eligible real records were retrieved by the PEER-NAG Database. The records came from six earthquake events and were comprised of 40 components taken from the 20 stations with two orthogonal components (see Table 2). The near-fault ground motions were not included in these records.

5.2. Scale Factors of Ground Motions

For the 3-, 9-, and 20-story structures at the three seismic hazard levels, the relative errors of SFs for the 40 record components are shown in Figure 2, which were calculated by the ASM method using equation (2) and by the LSM method using equation (4). The relative errors are defined as the differences between the SFs calculated by the LSM method and by the ASM method divided by the SFs calculated by the ASM method. For the 3-story structure, nearly all of the relative errors are positive. For the 9- and 20-story structures, the negative relative errors appear at the 50/50 seismic hazard level; however, more than half of the record components have positive relative errors. In general, for most of the record components, the SFs calculated by the LSM method are significantly larger than those calculated by the ASM method, and the relative errors are as high as 50% to 200%. Additionally, except for the 2/50 set for the 20-story structure, the relative errors of the SFs increase as the seismic hazard levels increase.

Figure 2 also shows that the relative errors of SFs for the 9-story structure are larger than those of the other two buildings. It should be noted that the larger values (absolute values) in the response spectra would contribute more to the error index SSEA (SSEL) than the lower values in the response spectra, when SSEA and SSEL are calculated using equations (1) and (3), respectively.

If both the target spectrum and the response spectrum of a record use arithmetic values, the value of the SSEA will mostly depend on the spectral values in the acceleration and velocity sensitive regions (corresponding to the short and medium-short periods, respectively) in the response spectrum; the SFs of the record will also be influenced. The degree of influence on SFs depends on the period range for spectral matching, which is determined by the structural periods.

While logarithmic values were used in the spectral matching procedure, the values of the SSEL and SFs mostly depend on the different period ranges from the ASM method. Because the values of the target spectral acceleration and the record’s spectral acceleration will both tend to zero at the very long periods, then their logarithmic values will go toward minus infinity. Because logarithmic functions have monotonicity, the spectral values in the displacement sensitive region (corresponding to long and very long periods) in the response spectrum will contribute more to the values of SF. The character of the LSM will benefit the seismic analysis of more flexible structures (e.g., high-rise buildings and long-span bridges) that have long first-mode periods.

For example, for the component I-ELC180 of the 9-story structure at the 2/50 hazard level, Figure 3(a) shows the relationship between the arithmetic values of the target spectral acceleration and the scaled record’s spectral acceleration (i.e., “” and “SF × ”); their corresponding logarithm values (i.e., “ln ” and “ln (SF × )”) are illustrated in Figure 3(b). The main control period ranges, in which the SSEA (SSEL) mostly depend on the spectral values, are presented in Figures 3(a) and 3(b). For clarity, the absolute values of “ln ” and “ln (SF × )” are illustrated by the dash lines (see Figure 3(b)). The period matching range ([0.29 s, 3.23 s]) determined by the fundamental period of the 9-story structure (T1 = 2.15 s) is shown in Figures 3(a) and 3(b).

5.3. Selection of Ground Motions

The number of records required for dynamic analyses has a significant effect on estimating the mean structural responses. Bommer and Acevedo [1] suggested that seven records might be adequate to produce an acceptable low dispersion of the structural response, while other researchers suggested that 10 records were required at least to obtain a stable mean of the structural dynamic responses. However, Reyes and Kalkan [12] considered that seven records were sufficient because the estimation accuracy of inelastic deformations of the SDOF system did not significantly increase when the number of records increased to 10. Current seismic design codes generally recommend seven (e.g., Eurocode 8 [5] and IBC-2012 [34]) or eleven (ASCE 7-16 [6]) as the required number of records. It should be mentioned that the minimum number of records required for dynamic analyses is affected by many situations, such as the purpose (i.e., the estimated mean or distribution of the structural response), the accuracy of predicted structural responses, the probability of structural collapse, and the degree of structural nonlinearity [35]. If the statistical distribution of the structural responses or the probability of structural collapse is considered in the analysis, the required number of records would be more, i.e., 30 and 60 [22, 23, 36]. This study aims to investigate the structural mean response demands (e.g., interstory drift ratios) calculated by the ASM and LSM methods. Therefore, the adopted number of records is seven herein. Actually, 11 records are also used and the obtained structural response results are similar to that obtained with seven records.

Additionally, the SFs of the records should also be limited. Generally, SFs exceeding 6-7 (or even less) are not adopted for scaling records in time-history analysis. These heavily scaled motions have no longer anything to do with real earthquake episodes that induced strong ground motions [37, 38]. Therefore, the limitations of SFs in this study are four at the 50/50 and 10/50 hazard levels and six at the 2/50 hazard level.

The procedure for records selection by the ASM and LSM methods is summarized as follows: (1) only the components with less than 4.0 or 6.0 are selected; (2) the selected components are preliminarily ranked by the values of (or ); (3) the optimal seven components with minimum (or ) are finally selected for the structural time-history analysis. Moreover, in order to avoid the influence of the interdependency between the two-components recorded by the same station on structural responses, only one horizontal component of a station is selected.

Using the above procedure, seven components were selected for the three structures by the ASM and LSM methods, which are listed in Tables 35. For the 3- and 9-story structures at the same seismic hazard level, four to six components are selected simultaneously by the ASM and LSM methods. However, there are only two to four of these types of components for the 20-story structure. Therefore, the ranking results generated with the ASM and LSM methods are different, especially for longer-period structures (i.e., 20-story structure). The relative errors of the SFs for the components selected simultaneously by both methods are also listed in the rightmost column of Tables 35, and nearly all absolute values of their relative errors are less than 20%.

6. Structural Response Analysis

6.1. Structural Benchmark Demands

In code-based seismic design checks as well as performance-based engineering assessments, the structural responses of interstory drift ratios are often taken as the EDPs. It is commonly assumed that EDPs are log-normally distributed [39]; therefore, the median value and the geometric mean of EDPs are consistent and many researchers used the median to represent structural seismic demands [11, 12, 22, 23]. However, some seismic design codes (e.g., Eurocode 8 and ASCE 7-16) use the average value (arithmetic mean) of EDPs if at least seven ground motions are required for time-history analysis. Therefore, the “true” structural responses—i.e., “benchmark seismic demands”—are defined as both the arithmetic mean and median (geometric mean) of the EDPs under the earthquake excitations of the 50/50, 10/50, and 2/50 sets of ground motions proposed by the SAC project (shown in Figure 4). Since the arithmetic mean of interstory drift ratios is a little larger than the geometric mean for all three structures at the three seismic hazard levels, it will not significantly affect the conclusions of the structural response analysis.

Pushover analysis is conducted to confirm the degree of the structural nonlinear response. As shown in Figure 5, the three structures are in a nonlinear stage when the roof drift ratios increase to 0.07 and must be in a nonlinear stage at the 10/50 and 2/50 hazard levels.

6.2. Structural Response Analysis

The relative errors are defined as the difference between the mean (median) of the structural responses due to the selected seven records and their corresponding benchmark demands divided by the benchmark demands. The relative errors of the peak interstory drift ratio (PIDR) demands over the heights of the three structures based on benchmark mean and median demands using the ASM and LSM methods at the three hazard levels are plotted in Figure 6.

According to either the mean or median of benchmark structural responses, the LSM method generally has a somewhat larger relative error (algebraic value) than the ASM method, if the same structures at the same seismic hazard are considered. Therefore, the LSM method may generate more “conservative” estimates of the structural response as compared to the ASM method. This result is consistent with the trend of the SFs which is that the SFs by the LSM method are usually larger than that by the ASM method. It is important to point out that nearly all of the relative errors of all groups are within the range of ±20%. In other words, both the LSM and ASM methods can be used to predict the mean or median of structural responses. In addition, the component TCU045-W is excluded from the statistical analysis because it could induce a very large displacement response of the 20-story structure at the 2/50 seismic hazard level by the ASM method (Figures 6(e) and 6(f)).

The coefficient of variation (COV), defined as the ratio between the sample standard deviation and the sample mean or median, was used to assess the effectiveness of the ASM and LSM methods to reduce dispersion in the PIDR demands. Figure 7 shows the COVs of the PIDR demands over the height of the three structures calculated by both the LSM and ASM methods at the three hazard levels. For the benchmark mean and median demands, the COVs are very similar for the same structure and at the same hazard level.

For long-period structures (i.e., 9-story and 20-story structures), the LSM method has lower COVs compared to the ASM method at all three hazard levels. This trend is even more pronounced for the 20-story structure with a more severe nonlinear response, and the COVs at the weakest floor of the 20-story structure by the LSM method are only 30% of that by the ASM method at the 2/50 hazard level. Moreover, if we only compare the COVs by the LSM method, the values of COV are very similar at all three hazard levels. Therefore, the structural nonlinear responses have little effect on the above trend as well as on the values of COV when using the LSM method. However, the COVs using the ASM method are significantly affected by the structural responses, and the values are larger for a more severe nonlinear response. Therefore, for longer-period structures, such as the 9- and 20-story buildings analyzed in this work, using logarithmic values of response spectra is more applicable for spectral matching than using arithmetic values.

The COVs of the 3-story structure are a little different with that of the 9- and 20-story structures, and the COVs are dependent on the structural responses. While the 2/50 ground motion set was used as input and the 3-story structure suffered the most severe nonlinear responses (shown in Figure 5), the COVs for the ASM method are larger than that of the LSM method. In contrast, the COVs by the ASM method are smaller than those of the LSM method, when the structure is subjected to the excitations of the other two ground motion sets selected at the 50/50 and 10/50 hazard levels.

7. Conclusions

The least square method is generally used to verify the compatibility between the spectrum of the selected ground motion and the target spectrum and to calculate the scale factors of the selected records. However, until now, there have not been any serious investigations performed regarding the effects of the coordinate values on the scaling of ground motions. This study aimed to investigate the influence of using arithmetic and logarithmic values in a spectral matching procedure (i.e., the ASM and LSM method) on the results of nonlinear structural time-history analysis. The 3-, 9-, and 20-story steel moment resisting frame structures, which represent low-, medium-, and high-rise buildings in the Los Angeles area in the SAC Steel Project, were used as examples. Based on the records that were selected and scaled by the ASM and LSM methods, the structural dynamic responses were compared at the three seismic hazard levels (i.e., 50%, 10%, and 2% probabilities exceeded in 50 years). The main conclusions are summarized below:(1)The SFs by the LSM method are significantly larger than those derived from the ASM method for low-, medium-, and high-rise buildings (i.e., 3-, 9-, and 20-story structures) at the three seismic hazard levels. This is mainly because the values of the SSEA for the ASM method mostly depends on the spectral values in the acceleration and velocity sensitive regions. Nevertheless, the values of the SSEL for the LSM method mostly depends on the spectral values in the displacement sensitive region. The trends will affect the values of SF. This indicates that the ASM method may be more suitable for short-period structures and the LSM method is more suitable for long-period structures. Meanwhile, the ranking and selection of the records are also different. There are only approximately four (for 3- and 9-story structures) or two (for 20-story structure) components simultaneously selected by both methods among the optimal seven components selected for nonlinear structural time-history analysis.(2)The relative errors of the PIDR demands can all be controlled within an acceptable range (±20%) by both methods, when the demands obtained by both methods are compared with the benchmark mean (or median) demands. Thus, the LSM and ASM methods both can be used to predict the mean or median of structural responses. However, the relative errors (algebraic value) calculated by the LSM method are larger than that of the AMS method; thus, the structural mean (or median) responses predicted by the LSM method will be slightly larger.(3)For long-period structures (i.e., 9-story and 20-story structures), the LSM method performed better at reducing variability in the structural responses results as compared to the ASM method. This advantage is more significant for longer-period structures (e.g., 20-story structure) with more severe nonlinear responses (e.g., 2/50 seismic hazard level).

Data Availability

Data used in this study can be found in the Pacific Earthquake Engineering Research Center Strong Motion Database at ngawest2.berkeley.edu (last accessed January 2019). All data, models, or code generated or used in this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 1778206).