The constant-strength ductility spectrum is a nonlinear response spectrum that is commonly used to establish the demand curve of the seismic response during performance-based seismic design. It is affected by many factors. In this paper, to evaluate the effect of the major influencing factors, including the structural parameters and seismic factors, the constant-strength ductility spectrum is calculated under different conditions based on 5535 ground motion records. Conclusions are drawn based on the mean constant-strength ductility spectra. (1) With respect to the effects of structural factors, the variation trend of the ductility spectra is highly consistent with increasing T, and the ductility spectra are usually larger for larger ζ and smaller Cy and k2. (2) With respect to the seismic factors, the ductility spectra show obvious differences in different periods; however, some parameters, such as PGA, have no influence on the ductility spectra. The results of this study can provide a theoretical basis for the calculation of ductility demand in the seismic design of structures.
Institute of Engineering Mechanics, China Earthquake Administration2019A01National Natural Science Foundation of ChinaU19392101. Introduction
Nonlinear response spectra are used to establish the demand curves of seismic responses. It is generally accepted that the most accurate method for the seismic performance evaluation and seismic response analysis of structures is the nonlinear time history analysis. In the 1970s, Newmark and Hall [1] proposed the equal energy principle and the equal displacement principle to study the relationship between the yield strength coefficient (Cy) and the displacement ductility demand (u). Subsequently, many experts have performed considerable research on the relationship between Cy and u [2–4]. In recent years, many studies have focused on the effects of structural parameters on constant-ductility seismic resistance spectra by a large number of ground motion records [5–9]. Xie and Zhai studied the influence of the restoring force model on constant-ductility seismic resistance spectra by using a large number of ground motion records; they mainly studied the variation in the seismic resistance coefficient under different ductility values, postyield stiffness ratio values of the bilinear model (k2), and damping ratios (ζ) [6, 7]. Similarly, some researchers have studied the effects of structural parameters and seismic factors on the constant-strength response spectra and constant-ductility response spectra [10, 11]. Miranda [12, 13] explained the influencing factors of magnitude, epicentre distance, and other factors on the constant-ductility response spectrum in detail. Yi and Zhang [14] compared the advantages and disadvantages of constant-ductility strength spectra and constant-strength ductility spectra; it was concluded that the method for calculating constant-strength ductility spectra is simpler and does not require iterative calculations. Additionally, only one u value can be obtained in a given period, while more than one value of Cy may be applied when the u value is fixed. In addition, Some people made improvements in structural parameters to research inelastic response demand. Michel et al. [15] proposed an improved method to the original EC 8 for displacement demand prediction. And the new proposals significantly improve the reliability of displacement demand prediction for values of strength reduction factors greater than 2 compared to the original EC 8 procedure. Castaldo et al. [16, 17] aimed to propose seismic reliability‐based relationships between the behavior factors and the displacement demand for nonlinear hardening and softening structures isolated by friction pendulum system devices and double concave sliding bearings. There are some researchers also studied the effects of different types of strong motions on constant-strength response spectra such as mainshock-aftershock sequence ground motion and pulse-like ground motions [18–20].
The above studies discussed the influence of some factors on constant-strength ductility spectra. But very few research studies have analyzed amplitude, duration, and frequency spectra parameters comprehensively, which is the most important parameters reflecting the characteristics of ground motion. Zhai and Li [21] studied the effects on the nonlinear displacement ratio spectra under various ground parameters such as PGA and PGV. But the effects of these factors on the constant-strength ductility spectrum were not analyzed.
Therefore, this paper conducted a comprehensive analysis of the factors that affect constant-strength ductility spectra. These factors include structural parameters and seismic parameters. The modified Clough model [22] is selected for calculating the constant-strength ductility spectra. The variation rule of constant-strength ductility spectra can be obtained under different structural parameters and seismic factors. By changing the factors, the constant-strength ductility spectra will change, and the variation rule can be analyzed.
2. Ground Motion and SDOF Structural Systems2.1. Ground Motions
The ground motion records are selected from the NGA-West2 database released by the Pacific Earthquake Engineering Research Center (PEER), which includes 21,540 groups of three-component ground motion records from 1935 to 2011 [23]. This paper mainly focuses on strong ground motions that are destructive to structures. Therefore, horizontal component ground motion records with peak ground accelerations (PGAs) greater than 50 gal are selected. The selected dataset includes 5,535 ground motion records. The space distribution of ground motion records is shown in Figure 1, and it can be concluded that the selected ground motion records are mainly distributed in the Mediterranean and the Pacific Coast regions. The magnitude-epicentre distance distribution is shown in Figure 2.
Data distribution of the selected ground motion records.
Magnitude-epicentre distance distribution of the selected ground motion records.
2.2. Structural Systems
There are many kinds of structural hysteretic model such as elastic perfectly plastic (EPP) model, modified Clough (MC) model, and pinching model [18, 24]. The modified Clough model as shown in Figure 3 is selected with reference to Clough’s study [25] for calculating the constant-strength ductility spectra in this paper. And it was used to simulate the flexural behavior that exhibits stiffness degradation at reloading. After the nonlinear response of the structure, the loading stiffness of the structure will gradually decrease with the reloading process. As shown in Figure 3, K0 is the initial stiffness and other parameters will be described in detail later.
The modified Clough model.
The seismic displacement ductility was chosen in this study as the parameter representative of damage due to large deformations. The definition of the displacement ductility u is relatively uniform. Most researches adopted the following definition [18–20]:(1)u=umuy,where um is the maximum displacement response of the SDOF system under ground motions and uy is the yield displacement of the structure.
The second structural parameter is the yield strength coefficient (Cy), which is defined as(2)Cy=fyfe,where fe is the strength demand on an infinitely elastic SDOF system, and fy is the yield strength of the corresponding inelastic SDOF system with the same mass and initial stiffness. Four values of Cy (Cy = 0.3, 0.5, 0.7, and 0.9) are considered to examine the influence of the structural strength.
Structural parameters also include the postyield stiffness ratio (k2) and damping ratio (ζ) and k2 is the ratio of the postyield stiffness to the initial stiffness in the hysteretic model. And ζ is defined as(3)ξ=c2mw,where c is damping coefficient and w and m are the natural frequency and mass of the system, respectively. In practical engineering, the damping ratio is usually 0.05.
3. Structural and Seismic Factors
To better understand the influence of strong motion records on the constant-strength ductility spectra, it is necessary to identify the various effects of the most important factors. The uncertainty of structural response under strong earthquake mainly comes from two aspects: (a) structural uncertainty and (b) uncertainty of ground motions. Therefore, this paper mainly analyzes these effects in terms of structural parameters and seismic factors. The detailed influencing factors are shown in Table 1.
Influencing factors and parameters.
Different aspects
Parameters
Structural parameters
Period (T)
The yield strength coefficient (Cy)
Damping ratio (ζ)
The postyield stiffness ratio (k2)
Seismic factors
Earthquake magnitude (M)
Epicentre distance (R)
Peak ground acceleration (PGA)
Effective peak velocity (EPV)
Significant duration (Ds)
4. Effect of Structural Parameters on the Constant-Strength Ductility Spectra
The structural parameters of a SDOF system have certain influence on the structural response during an earthquake. The structural parameters mainly include T, Cy, ζ, and k2. Some researchers have studied the influence of structural parameters on the constant-strength ductility spectra [14, 20, 26]. When the structural parameters change, the constant-strength ductility spectra change accordingly. Therefore, according to the existing research results, this paper explains the effect of these parameters on the constant-strength ductility spectra.
The constant-strength ductility spectra are calculated as shown in Figure 4 based on 5535 ground motion records under specific conditions: Cy = 0.3, ζ = 0.05, k2 = 0. It can be seen from Figure 4 that the trend of the mean constant-strength ductility spectra is consistent with that of all constant-strength ductility spectra. Therefore, when explaining the influence of different factors on the constant-strength ductility spectra, we mainly study the effect of different factors on the mean constant-strength ductility spectra.
Constant-strength ductility spectra of selected ground motion records.
4.1. Effect of the Period and the Yield Strength Coefficient
To study the influence of T and Cy on the constant-strength ductility spectra, the mean constant-strength ductility spectra under different Cy values is calculated when ζ is 0.05 and k2 is 0, as shown in Figure 5. Four values of Cy are considered (Cy = 0.3, 0.5, 0.7, and 0.9), and the following conclusions are drawn from the analysis.
Effect of different yield strength coefficients (Cy = 0.3, 0.5, 0.7, and 0.9) on the variation of the mean constant-strength ductility spectra.
For constant period, the smaller the value of Cy, the larger the values of u. For constant Cy, the values of u decrease fast with the increase of the T in short-period section. The values of u do not change significantly with the period in the medium- and long-period ranges, particularly for bigger Cy values, and thus the constant-strength ductility spectra can be characterized as period independent. This variation rule is the same as the result of Li and Yi’s research [14, 20].
4.2. Effect of the Damping Ratio
To analyse the influence of ζ on the constant-strength ductility spectra, the variation in the mean constant-strength ductility spectra under different ζ values is calculated when Cy is 0.3 and k2 is 0. The variation is shown in Figure 6. Four values of ζ (ζ = 0.01, 0.05, 0.1, and 0.15) are considered. The following conclusions are drawn from the analysis.
Effect of ζ (ζ = 0.01, 0.05, 0.1, and 0.15) on the variation of the mean constant-strength ductility spectra.
In short-period range, the variation rule of the constant-strength ductility spectra is complex. However, in medium- and long-period ranges. When ζ changes, the displacement ductility value changes obviously, and the smaller the value of ζ is, the slower the value of u decreases. Additionally, the rate of decrease of u varies based on the damping ratio. Therefore, the effect of ζ on the constant-strength ductility spectra cannot be neglected.
4.3. Effect of the Postyield Stiffness Ratio of the Bilinear Model
To quantitatively analyse the effect of k2 on the constant-strength ductility spectra, the variation in the constant-strength ductility spectra under different k2 conditions is calculated when Cy is 0.3 and ζ is 0.05. Four values of k2 (k2 = 0, 0.05, 0.1, and 0.15) are considered. The variation is shown in Figure 7. The following conclusions are drawn from the analysis.
Effect of k2 (k2 = 0, 0.05, 0.1, and 0.15) on the variation of the mean constant-strength ductility spectra.
The variation rule of the constant-strength ductility spectra is obviously consistent at all periods. For different k2 values, with the increase in k2, the constant-strength ductility spectra values gradually decrease for a constant period.
5. Effect of Seismic Factors on the Constant-Ductility Strength Spectra
Seismic factors mainly include the earthquake magnitude, epicentre distance, and ground motion parameters. These factors also influence the constant-strength ductility spectra value. Before analyzing the influence of seismic factors on the constant-strength ductility spectra, the ground motion records are classified according to the site type, and the structure parameters are defined as follows: (1) Cy = 0.3; (2) ζ = 0.05; (3) k2 = 0.
Using the USGS site classification method, 5535 ground motion records are classified into five site types, and the classification result is shown in Table 2. According to the classification results, the number of ground motion records varies greatly by site type, and some sites have very few ground motion records. To eliminate the influence of the number of ground motion records on the constant-strength ductility spectra, site A type, site B type, and site C type are combined as stiff soil site, and site D type and site E type are combined as soft soil site. The earthquake magnitude and epicentre distance distribution of records in different site classifications are shown in Figure 8.
Classification results of ground motion records based on USGS site classification.
Site type
Classification standard Vs30 (m/s)
Number 1
Site type
Number 2
Site A type
>1500
7
Stiff soil
2850
Site B type
760∼1500
101
Site C type
360∼760
2742
Site D type
180∼360
2533
Soft soil
2685
Site E type
<180
152
(a) Earthquake magnitude-Vs30 distribution and (b) epicentre distance-Vs30 distribution in different site classifications.
5.1. Effect of Earthquake Magnitude
To analyze the effect of earthquake magnitude M on constant-strength ductility spectra, the ground motion records are divided into four types according to earthquake magnitude, that is, M < 5, 5 ≤ M < 6, 6 ≤ M < 7, and M ≥ 7, and the classification results are shown in Table 3. The mean constant-strength ductility spectra are calculated in each type, and the variation rule of the mean constant-strength ductility spectra with different magnitudes is analyzed as shown in Figure 9. The following conclusions can be drawn from the analysis.
Classification results of ground motion records based on earthquake magnitude.
Site type
M < 5
5 ≤ M < 6
6 ≤ M < 7
M ≥ 7
Stiff soil site
301
556
1441
552
Soft soil site
226
589
1323
547
Effect of earthquake magnitude on the variation of mean constant-strength ductility spectra: (a) stiff soil site; (b) soft soil site.
The mean constant-strength ductility spectra show obvious differences for different earthquake magnitudes. For the results in stiff soil site, in short and medium ranges, the greater the magnitude, the larger the ductility spectra value. In long-period range, the ductility spectra values of different magnitudes overlap, and there is no obvious variation. For the results in soft soil sites, in short and medium ranges, there is a variation similar to that of the ductility spectra values in stiff soil sites. However, in long-period range, the larger the magnitude, the smaller the ductility spectra value.
5.2. Effect of Epicentre Distance
To analyze the influence of the epicentre distance on the constant-strength ductility spectra, the ground motion records of every site type are divided into two groups according to epicentre distance: (1) near-field; (2) far-field. The classification result is shown in Table 4. The variation rule of the mean constant-strength ductility spectra with different epicentre distances is analyzed in stiff and soft soil sites, as shown in Figure 10. The following conclusions can be drawn from the analysis.
Classification results of ground motion records based on epicentre distance.
Site type
Type
Epicentre distance (km)
Numbers
Stiff soil
Near-field
<30
1227
Far-field
>30
1623
Soft soil
Near-field
<30
936
Far-field
>30
1749
Effect of epicentre distance on the variation of mean constant-strength ductility spectra: (a) stiff soil site; (b) soft soil site.
The mean constant-strength ductility spectra of the different types show obvious differences. For ground motions in stiff soil, in short- and medium-period ranges, the farther the epicentre distance, the larger the ductility spectra value. In contrast, in long-period ranges, the ductility spectra value is usually smaller in the near-field than in far-field. And this phenomenon is more obvious in soft soil sites.
5.3. Effect of Ground Motion Parameters on Constant-Ductility Strength Spectra
The nonlinear responses of structures are closely related to the ground motion parameters. To explain the effects of different ground motion parameters on constant-strength ductility spectra, the representative parameters of amplitude, frequency spectra, and duration are chosen to analyse the variation of constant-strength ductility spectra in this paper. Therefore, in this paper, the representative parameters of amplitude, frequency spectra, and duration are chosen to analyse the variation of constant-strength ductility spectra. The definitions of the structural parameters are the same as those in the previous section.
The selected representative ground motion parameters are the peak ground acceleration (PGA), the effective peak velocity (EPV), and the significant duration (Ds5−95). For seismic structural analysis purposes, the most simple and widely used intensity index is PGA, and it is directly obtained from the ground motion record itself. The calculation of EPV is shown as follows:(4)EPV=Sv2.5,where Sv is the mean velocity response spectrum with a damping ratio of 0.05 near the period of 1 second, and the period range is usually from 0.8 to 1.2 seconds.
For the Ds5−95, Husid [27] proposed an intensity parameter that indicates the increase of ground motion energy over time. Its standardized form is shown as follows:(5)It=∫0ta2tdt∫0Ta2tdt,where the T is the total duration of ground motion and I (t) is a function from 0∼1. Trifunac and Brady [28] defined a significant duration of ground motion as follows:(6)IT1=0.05,IT2=0.95,Ds5−95=T2−T1,where T1 and T2 are, respectively, the times where 5% and 95% of ground motion energy are reached. The significant duration Ds5−95 as defined in equation (6) is most widely used for evaluating the effective duration of earthquakes.
According to these parameters, 5535 ground motion records are classified. The classification results are shown in Table 5, and the PGA-EPV and PGA-Ds5-90 distribution of ground motion records are shown in Figure 11.
Classification results of ground motion records based on ground motion parameters.
Ground motion parameters
Range of values
Numbers
PGA
≤200 gal
4296
>200 gal
1239
EPV
≤10 cm/s
3728
>10 cm/s
1807
Ds5−95
≤25 s
4456
>25 s
1079
(a) PGA-EPV distribution and (b) PGA-Ds5−90 distribution of selected ground motion records.
5.3.1. Effect of PGA
To study the effect of PGA on the constant-strength ductility spectra, the ground motions are divided into two categories: (1) PGA > 200 gal and (2) PGA ≤ 200 gal. The variation of the mean constant-strength ductility spectra is analyzed, and the results are shown in Figure 12. It can be concluded that PGA has little effect on the mean constant-strength ductility spectra values at all periods.
Effect of PGA on the variation of mean constant-strength ductility spectra.
5.3.2. Effect of EPV
To study the effect of EPV on the constant-strength ductility spectra, the ground motion records are divided into two categories: (1) EPV ≤ 10 cm/s and (2) EPV > 10 cm/s. The variation of the mean constant-strength ductility spectra is analyzed, and the results are shown in Figure 13. It can be concluded that EPV has a significant effect on the mean constant-strength ductility spectra values in short-period range; the larger the EPV of the ground motion record, the larger the ductility spectra. But in the medium- and long-period ranges, EPV has little effect on the ductility spectra values.
Effect of EPV on the variation of mean constant-strength ductility spectra.
5.3.3. Effect of Ds5−95
To analyze the influence of Ds5−95 on the constant-strength ductility spectra, the ground motions are divided into two categories: (1) Ds5−95 ≤ 25 s and (2) Ds5−95 > 25 s. The variation of the mean constant-strength ductility spectra is analyzed, and the results are shown in Figure 14. In the long-period range, long-duration ground motion records usually have larger spectra values than short-duration ground motion records. But in short- and medium-period ranges, the longer the duration of ground motion, the larger the ductility spectra.
Effect of Ds5−95 on the variation of the mean constant-strength ductility spectra.
6. Conclusion and Discussion
This study compares the influences of different factors on the constant-strength ductility spectra based on the 5535 ground motion records in NGA-West2 database, and the effects of structural parameters and seismic factors on the spectra are analyzed. The following conclusions are inferred from the results of this study:
With respect to the effects of structural parameters, except for ζ at short-period ranges, the variation of the ductility spectra value is highly consistent with the increasing T. And the ductility spectra values are usually larger for larger ζ and smaller Cy and k2.
With respect to the effects of seismic factors, the ductility spectra values show obvious differences in different period ranges. In short-period range, when the epicentre distance and magnitude of the ground motion records are larger, the corresponding ductility demand is higher. However, for EPV and Ds5−95, when the EPV and Ds5−95 of the ground motion records are larger, the corresponding ductility demand is small. In the medium- and long-period ranges, the ductility spectra value varies under different conditions; however, some parameters, such as EPV and PGA, have little effect on the ductility spectra values.
As one of the nonlinear response spectra, the constant-strength ductility demand spectra are affected by many factors; in this paper, we studied as many influencing factors as possible from the structural and seismic aspects. The effect of a portion of factors studied in this paper are compared with the results of the former studies, while the effect of some other factors were not compared due to lack of analysis of these factors in the previous studies [14, 20]. So it is still necessary to compare the results with the others in the future. In conclusion, the results of this study provide a statistical and quantitative basis for the calculation of ductility demand in the seismic design of structures.
Data Availability
Strong motion records were provided by the Pacific Earthquake Engineering Research Center (PEER).
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Scientific Research Fund of the Institute of Engineering Mechanics, China Earthquake Administration (Grant no. 2019A01). and the National Natural Science Foundation of China (Grant no. U1939210). The authors also appreciate the assistance of the Pacific Earthquake Engineering Research Center (PEER) in providing the strong motion records.
NewmarkN. M.HallW. J.Seismic design criteria and nuclear reactor facilities197346Gaithersburg, MD, USANational Bureau of Standard, U.S. Department of Commerce209236NassarA. A.KrawinklerH.Seismic demands for SDOF and MDOF systems. Report no. 951991954Stanford, CA, USAStanford University1245PalS.DasakaS. S.JainA. K.Inelastic response spectra198725333534410.1016/0045-7949(87)90126-x2-s2.0-0023241088VidicT.FajfarP.FischingerM.Consistent inelastic design spectra: strength and displacement199423550752110.1002/eqe.42902305042-s2.0-0028184459PengB.-F.ConteJ. P.Statistical insight into constant-ductility design using a non-stationary earthquake ground motion model199726989591610.1002/(sici)1096-9845(199709)26:9<895::aid-eqe683>3.0.co;2-gXieL.-l.ZhaiC.-h.Study on the severest real ground motion for seismic design and analysis200316326027110.1007/s11589-003-0030-92-s2.0-0038443283ZhaiC. H.XieL. L.ZhangM. H.Influence of damping on the spectra of constant-ductility seismic resistance of engineering structures2006381017051708in ChineseXueS. L.YuZ.XuG.Elastoplastic dynamic time history analysis of Yinchuan De Feng building20174753742in ChineseHeH. jXuM. H.SuL.Study on the number of input ground motion in elastoplastic time history analysis of structures2018382150156in ChineseZhouD. S.LuX. L.Ducitility demand spectra and inelastic displacement spectra considering soil conditions and design characteristic periods20042413948in ChineseRupakhetyR.SigbjörnssonR.Ground-Motion Prediction Equations (GMPEs) for inelastic displacement and ductility demands of constant-strength SDOF systems20097366167910.1007/s10518-009-9117-62-s2.0-67651036551MirandaE.BerteroV. V.Evaluation of strength reduction factors for earthquake-resistant design199410235737910.1193/1.1585778MirandaE.Inelastic Displacement ratios for structures on firm sites2000126101150115910.1061/(asce)0733-9445(2000)126:10(1150)2-s2.0-0034298759YiW.-J.ZhangH.-Y.KunnathS. K.Probabilistic constant-strength ductility demand spectra2007133456757510.1061/(asce)0733-9445(2007)133:4(567)2-s2.0-33947314620MichelC.LestuzziP.LacaveC.Simplified non-linear seismic displacement demand prediction for low period structures20141241563158110.1007/s10518-014-9585-12-s2.0-84904321667CastaldoP.PalazzoB.AlfanoG.PalumboM. F.Seismic reliability-based ductility demand for hardening and softening structures isolated by friction pendulum bearings20182511e225610.1002/stc.22562-s2.0-85053438003CastaldoP.AlfanoG.Seismic reliability-based design of hardening and softening structures isolated by double concave sliding devices202012910593010.1016/j.soildyn.2019.105930ZhaiC.-H.WenW.-P.LiS.ChenZ.ChangZ.XieL.-L.The damage investigation of inelastic SDOF structure under the mainshock-aftershock sequence-type ground motions201459304110.1016/j.soildyn.2014.01.0032-s2.0-84893499979ZhaiC.JiD.WenW.LeiW.XieL.Constant ductility energy factors for the near-fault pulse-like ground motions201721234335810.1080/13632469.2016.11575292-s2.0-84966696483LiC.KunnathS.ZhaiC.Influence of early-arriving pulse-like ground motions on ductility demands of single-degree-of-freedom systems201822912410.1080/13632469.2018.14667442-s2.0-85046141168ZhaiC.LiS.SunY.Study on inelastic displacement ratio spectra for near-fault pulse-type ground motions20076435135510.1007/s11803-007-0755-x2-s2.0-39349111246ChopraA. K.20175thHoboken, NJ, USAPearson Education Inc.TimothyD. A.RobertB. D.JonathanP. S.NGA-West2 Database20143039891005MirandaE.Ruiz-GarciaJ.Influence of stiffness degradation on strength demands of structures built on soft soil sites200224101271128110.1016/s0141-0296(02)00052-42-s2.0-0036788001CloughR. W.PenzienJ.1975New York, NY, USAMcGraw-HillShiQ. X.MenJ. J.Ductility spectra method with constant ζ_y inperformance-based seismic evaluation of building structures20072715458in ChineseHusidR. L.Analisis de terremotos:analisis general1969812142TrifunacM. D.BradyA. G.A study on the duration of strong earthquake ground motion1975653581626