ANovelMethod for ConstructingMain-Aftershock Sequences and Its Application in the Global Damage Accumulation Effects Analysis of Gravity Dams

*e traditional linear elastic and Drucker–Prager (DP) models cannot truly reflect the strong nonlinear characteristics of the concrete and rock foundation of the dam under earthquake.*erefore, for comprehensive evaluation of the cumulative damage of the gravity dam structure caused by aftershock, the dynamic damage of the dam body concrete is analyzed by many scholars through the plastic damage mechanics method, but there is little research on rock material at the dam foundation with the method utilized; thus, the simulation of the whole dynamic damage evolution is worthy of investigation of the dam body and dam foundation. According to the randomness of ground motion, the transcendental probability (P) is introduced to express the statistical characteristics of aftershock intensity, and a newmethod for constructing main-aftershock sequences of groundmotion is proposed in this paper. And then, the law of the damage evolution and energy characteristics of the concrete gravity dam under the combined action of the main shock and aftershock sequences is studied. *e results are shown as follows: the smaller aftershocks do not cause further damage to the dam; as the aftershock intensity increases, the energy characteristics of the dam body and foundation have shown different changing rules; when the ratio of peak aftershock acceleration to peak main shock acceleration (∇PGA) approximately equals 0.68, the aftershock will cause larger secondary damage to the dam.


Introduction
e measured data show that strong earthquakes are often followed with lots of aftershocks. After the severe damage of the structure caused by the strong main shock, the damage of the structure will be aggravated by the strong aftershock.
is repeated damage will be undoubtedly worse than the damage caused by the single main shock, which calls for higher requirements on the seismic performance of the structure [1]. Many post earthquake investigations also reveal that the damage to structures caused by aftershocks can not be ignored [2][3][4][5]. However, the current seismic design code [6] only considers the action of the main shock, while the secondary damage of the dam structure caused by the strong aftershocks is ignored. is may lead to insufficient earthquake resistance of the dam structure under the main-aftershock sequence. Once the dam is damaged in the earthquake, it will not only endanger the safety of the project itself but also bring immeasurable loss to people's life and property safety in the downstream. Studying the effects of aftershocks on the dynamic characteristic of the concrete gravity dam is helpful to comprehensively and accurately assess the seismic resistance of large dam projects.
Due to the lack of the record of the main-aftershocks in practice, many scholars mostly use the main-aftershocks of tectonics in their research. erefore, it is vital to developing a reasonable method of seismic ground motion construction based on main-aftershocks sequences. For major projects, according to the probability analysis of seismic hazards, the main seismic parameters of the foundation rocks of the engineering site under different transcendence probabilities can be obtained. So, the core of constructing the main-aftershock sequence is to determine the seismic parameters of the aftershock sequence. It should be noted that the magnitude and location of aftershocks depend on the main shock [7]; in this way, more consistency can be maintained in the probability seismic risk analysis of the main-aftershock. Different scholars have proposed different methods of seismic ground motion construction. Das et al. [8] studied the aftershock and main shock ground motion parameters and gave the relationship between the two parameters; however, they only chose Chi-Chi earthquake records and did not consider the saturation effect in the near-field region. Boore et al. [9] did not consider the relationship between the main shock and aftershock ground motion parameters. Graizer and Kalkan [10] did not tell the difference between the main shock and aftershock when establishing the attenuation relationship of ground motion. Vahedian et al. [11] used the artificial neural network method to generate aftershock ground motion. Actually, due to the randomness of earthquakes and the complexity of their seismogenic mechanism, probabilistic analysis can better describe the actual situation. Combined with the seismic risk analysis method, it should be more practical to analyze the effect of aftershocks with different transcendental probabilities on the cumulative damage of structures. However, in the existing methods of main-aftershock sequences ground motion construction, no scholars have quantitatively described the transcendental probability of aftershocks with different intensities. In addition, the influence of duration and intensity envelope is also rarely considered. is paper considers the above factors comprehensively and proposes a new construction method of the main-aftershock sequences ground motions.
At present, scholars have analyzed and studied the seismic performance of the structure under the action of the main-aftershock. Alliard and Leger [12] studied the effect on the drainage efficiency of gravity dams considering the aftershocks. Wang and Zhang [13] analyzed the effects of the main-aftershock seismic sequences on the dynamic characteristic of gravity dams. Wang et al. [14], respectively, analyzed the damage evolution of the gravity dam under the action of the main shock, aftershock, and main-aftershock. Abdollahzadeh et al. [15] evaluated the seismic performance of reinforced concrete structures under the action of main and residual earthquakes. In the current study, the linear elastic and Drucker-Prager (DP) model is mainly used to simulate the dam foundation. Actually, microfissures are randomly distributed in the dam foundation rock mass, whose tensile strength is lower than that of the dam concrete. erefore, the damage analysis to dam foundation rock is an important part for the nonlinear dynamic analysis of the high-dam multicoupling system [16]. At present, there is no literature about the influence of aftershock on the whole dam foundation damage under different transcendental probabilities by using the damage mechanics model. In this paper, the whole damage model of the gravity dam is established based on an actual project. Using the construction method proposed in this paper, main-aftershock earthquake waves are constructed and synthesized. e cumulative damage of gravity dams subjected to strong aftershocks with different transcendental probability after the main earthquake is analyzed. is study quantifies the effect of aftershocks on the whole dynamic damage to the dam body and foundation under different transcendental probabilities and reveals the evolution law of the whole damage of dam body and dam foundation under different aftershocks.

Study on Construction Method of the Main-Aftershock Sequence
In this paper, the main-aftershock (mainshock-aftershock) sequence is defined as the combination of the largest aftershock after the main shock and the main shock; that is, the sequence of two shocks is superimposed. ere will be dozens or even hundred of aftershocks after the strong main shocks; some scholars have found that smaller aftershocks have less impact on the further damage of the dam. erefore, the main-aftershock sequence relationship in this paper mainly considers the influence of the largest aftershock after the main shock. is paper combines the statistical results of the main shock and aftershock and attenuation relation of aftershock ground motion, introducing transcendental probability function to express statistical characteristics of aftershock intensity. At the same time, the attenuation law of the envelope function and the effect of the duration of the ground motion are considered. e process of the proposed seismic method for the mainaftershock sequence is as follows.
(1) Determine the (max) aftershock magnitude according to the statistical relationship between the main earthquake magnitude and the aftershock magnitude. (2) Choose the appropriate attenuation relationship of aftershock ground motion according to the magnitude, the epicenter distance, and site characteristics of the main shock. e transcendental probability (P) is introduced to express the aftershock intensity. e parameters of aftershock such as PGA and PSA under different transcendental probability are predicted.
(3) According to the magnitude and the epicenter distance of the main shock and aftershock obtained, determine the duration and intensity envelope functions. e main-aftershock sequence can be obtained by combining the time history of the main and aftershock seismic sequences.

Magnitude Statistical Relationship between Main Shock and Aftershock.
Magnitude is one of the most commonly used parameters to describe the magnitude of earthquakes and the characteristics of earthquakes. In [17], based on the definition of strong aftershock, data of 123 earthquakes with magnitude 6 or above satisfying the principle of defining strong aftershock in the world are collected, after the regression analysis by the least square method. e regression relationship between the main shock magnitude (M m ) and the strong aftershock magnitude (M a ) is obtained as follows: (1) Using the statistical formula of variance, the unbiased estimation of variance for linear regression equation (1) is obtained σ 2 � 0. 18. e difference between the maximum aftershock magnitude and the main shock magnitude is mostly in accordance with the bath's law of about 1.2.

Attenuation Relation of Aftershock Seismic Ground
Motion Parameters. Zhai and Xie [18] selected 1353 main shock and aftershock ground motions from 19 main-aftershocks sequences to analyze the characteristics that the ratio of aftershock ground motion parameters to main shock ground motion parameters changing with magnitude, the fault distance, and site. Based on the existing attenuation relations, the following attenuation relation of the aftershock ground motion parameters is proposed: (2) In this formula, ∇Y represents the ratio of the aftershock parameter (Y as ) to the main shock parameter (Y ms ); ∇M is the ratio of aftershock magnitude (M as ) to main shock magnitude (M ms ); ∇D is the ratio of aftershock fault distance (D as ) to the main shock fault distance (D ms ); V 30 is the average shear-wave velocity of 30 m underground; and ε is described as the residual between the predicted and observed values. Assuming that it obeys the normal distribution of mean value 0 and standard deviation σ, b 1 ∼b 5 are the fitting coefficients of the formula. e first and second terms on the right side of the equation are used to describe the variation of ∇Y with the magnitude of the main shock and aftershock. In the third term, ln∇D denotes the geometric decay of ∇Y with the ratio of fault spacing. ∇M/D ms is used to simulate the near-field saturation effect. And, the fourth item is used to simulate the effect of site type on ∇Y. Equation (2) is relevant to the main shock parameters; thus, the model can conveniently predict the parameters of aftershock ground motion, as well as reveal the relationship among the parameters of the main shock and aftershock.
In fact, because of the randomness of earthquakes and the complexity of their seismogenic mechanism, when using statistical laws to predict the parameters of aftershock ground motions, the probability analysis method can better describe the actual situation.
us, the statistical law of aftershock intensity expressed by the peak acceleration transcendence probability P of aftershock is proposed in this paper when constructing the parameters of aftershock ground motion. In the calculation, it is assumed that the epicenter location of the main shock is close to that of the aftershock.
Relative parameter ∇Y of aftershock ground motion in equation (2) is expressed as relative peak acceleration of the main-aftershock ground motion ∇PGA (the ratio of aftershock peak acceleration to main shock peak acceleration). Given that ε∼N(0, σ 2 ), z + ε∼N(z, σ 2 ), joint equations (2) and (3) can be obtained ∇PGA � e z + ε , where ∇PGA obeys lognormal distribution. e probability density function equation (4) and cumulative distribution function equation (5) of relative peak acceleration ∇PGA of the main shock and aftershock ground motions can be obtained by calculation.
Cumulative distribution function en, ∇PGA's transcendence probability P is derived out, as the following equation: When calculating, different transcendental probabilities P can be selected to determine the corresponding ground motion strength parameter ∇PGA for analysis.

Law of Attenuation of Earthquake Intensity.
Earthquake duration can be divided into three phases, the initial ground motion increases rapidly from small to large, followed by the continuous stage of strong earthquake and then the slow attenuation stage. Intensity nonstationary is the characteristic of seismic wave changing with time. If the peak points of acceleration time history are connected, the envelope function f(t) can be obtained. e time course of earthquake motion a(t) can be expressed as follows: where x is expressed as a stationary random process. Envelope function f(t) adopts the segmented model proposed by Amin and Ang [19]: where t 1 and t 2 represent the start and end time of strong earthquake duration, respectively; c represents the attenuation coefficient; the peak stationary time t s � t 2 − t 1 ; t d is the total holding time; and attenuation time of seismic e envelope parameters t 1 , t 2 , and c are related to the earthquake mechanism, the propagation path of ground motion, and the site parameters. Huo [20] selected 457 strong earthquake records for statistical analysis and obtained the regression reationship among envelope parameters on bedrock site, magnitude, and epicenter distance.
where M and R are, respectively, represented as magnitude and epicenter distance.

An Example of the Main-Aftershock Sequence
Construction. Taking a concrete gravity dam project as an example, the main-aftershock sequence ground motion is constructed. In the project, the design earthquake is based on the ground motion parameters with a probability of 2% over the base period of 100 years, and the peak acceleration of the site design bedrock is 0.316 g. e corresponding earthquake magnitude M m is 7.3. e epicenter distance D ms is 25 km, and the average underground shear-wave velocity V 30 of 30 m in bedrock is 1300 m/s. e aftershock magnitude can be obtained from equation (1), and the maximum aftershock magnitude is 6.09 (∇M � 0.843). In reference [21], homologous earthquakes are selected for main shock and aftershock, corresponding ∇D � 1. According to equation (3), Z � − 0.735. e standard deviation σ � 0.569 is known from reference [18]. Substitute the data into equations (4) and (6). e probability density function and cumulative distribution function are derived out, as shown in the following equations: Figures 1 and 2 are probability density function curves and transcendental probability function curves, respectively. is paper takes the aftershock PGA transcendence probability P � 50% as an example, corresponding relative acceleration of the main shock and aftershock is ∇PGA � 0.479. Let ∇Y � ∇Sa, and take the site bedrock acceleration response spectrum (PSA) as the main shock acceleration response spectrum, substituting it into equation (2) to obtain the acceleration response spectrum of aftershock site, as shown in Figure 3. According to the attenuation relation of strength envelope equation (5), the duration of the main shock and aftershock ground motions is calculated to be 33s and 22s, respectively. e main-aftershock earthquake waves with nonstationary characteristics are synthesized (with a damping ratio of 5%). e interval between the main shock and aftershock ground motions is 10 s, and then the seismic waves are combined. Taking ∇PGA � 0.479 as an example, the acceleration time history curve of horizontal principal residual sequence seismic wave is shown in Figure 4. T1, T2, and T3 in the figure represent the duration of the main shock, the time interval between the main-aftershock, and the duration of the aftershock, respectively. Peak acceleration of vertical seismic wave takes 2/3 of horizontal.

Plastic Damage Model of Concrete.
Model constitutive relationship is as follows [22]: where σ and d are expressed as stress and damage variables, respectively; ε and ε pl are the elastic and plastic strains; and D el 0 is the initial elastic stiffness. e elastic modulus after damage is calculated by the following formula: e yield equation is expressed as follows: where α and c are material constants; p � − (1/3)trσ; σ is the effective stress tensor parameters; σ c and σ t are effective compressive stress tensor and tensile stress tensor, respectively; σ ⌢ max is the maximum eigenvalue of σ; and qis the effective Mises equivalent stress.

Model Validation.
A simple uniaxial tension compression and tension compression cycle problem is numerically verified and compared with the existing experiments [22,23]. e model adopts the following material parameters: modulus of elasticity E 0 � 3.1 × 104 MPa, tensile strength F t � 3.48 MPa, and fracture energy G f � 40 N/m. e unit feature length L t � 82.6 mm. Figure 5 shows the comparison between the uniaxial cyclic load test data and the numerical simulation results. e test data are the local area data measured by the wide strain gauge; if the average stress-strain relationship in the small local area is considered, the absolute value of the post peak curve slope of the final test data will be reduced. Since the stress curve is assumed to be the difference of two exponential functions, the slope of the numerical simulation curve after the peak in the tensile stress-strain curve is more gentle than the test data. It can be seen that the numerical results can well simulate the stress-strain curve of concrete under uniaxial tension cyclic load. e stiffness degradation of the model material is simulated in each unloading and loading cycle. It can be seen from Figure 6 that the stiffness decreases after tensile damage under cyclic load; the "unilateral effect" of concrete can be well reflected by the continuous reverse pressure and the recovery of elastic modulus in the compression zone.   Shock and Vibration method of seismic ground motion construction for the main-aftershock sequence proposed in Section 1, the simulated earthquake wave of the main-aftershock sequence is constructed and synthesized. e nonlinear dynamic response of concrete gravity dams with different aftershock transcendental probabilities is studied.

Computational Model and Method.
Select a dam section of the above concrete gravity dam as the analysis object. e height of the dam is 97m, and the water depth of the normal pool level is 83m. C20 concrete is used below 50m of the dam body, and C15 concrete is used above 50m of the dam body. Geological survey shows that 0-10 below the dam site is class III rock, and the others belong to class II rock. e finite element model of the gravity dam was established based on the ABAQUS finite element software (Figure 7). Both the dam body and the dam foundation use CPE4R elements. e model is divided into 8913 nodes and 8658 elements. Material parameters of class II rocks are f � 1.25 and c � 1.6 MPa (f � tanφ, same as below). Material parameters of class III rocks are f � 1.10 and c � 1.15 MPa. Tensile strength of dam foundation rock mass is deduced by using f t � 2c·cosφ/(1 + sinφ) in M-C criterion.
Considering both vertical and horizontal seismic waves, the mass damping coefficient and stiffness damping coefficient are calculated with the damping ratio of 5%. Viscoelastic artificial boundaries are set on both sides and bottom of the model foundation to simulate the influence of wave energy dissipation to infinite ground under the action of earthquakes. Assuming the incompressibility of the reservoir water, the Westergaard additional mass method is used to simulate the hydrodynamic pressure. Reference [24] obtains the damage parameters of the dam concrete material and obtains the damage parameters of the bedrock material by corresponding reduction of the concrete damage curve relative to the tension strength of rock [16]. After considering the effect of seismic energy escaping to far-field, damage response of dam body-foundation integral structure under earthquake action of the main and aftershock is analyzed. e concrete and bedrock material parameters used in the calculation are shown in Table 1.

Effects of the Main-Aftershock on Damage Evolution of Gravity Dams.
In order to study the damage evolution law of gravity dams under different main-aftershock sequences, based on the sequence of the main-aftershock constructed in Section 2.4, this paper increases the value of transcendence probability P for the aftershock intensity. Given that aftershock ground motion intensity is generally smaller than main shock ground motion; the values of P in this paper are 60%, 40%, 30%, 25%, 20%, 15%, and 10%, respectively. By equation (2), the corresponding ∇PGA is 0.415, 0.555, 0.645, 0.705, 0.775, 0.865, and 0.990, respectively. Eight kinds of conditions including single action of main shock and combined action of main and aftershocks are selected for calculation and research. e calculation results are shown in Figure 8. e blue area indicates that no damage has occurred, and the red area indicates that damage cracks have occurred at this location. Under different conditions, the damage locations of concrete gravity dams are roughly the same. e damage cracks of the concrete material of the dam body mainly appear in the downstream slope change, and the bedrock cracks at the dam heel extend downward along the depth direction. Under the action of single main shock, the damage crack at the downstream slope change extends to upstream and downward. Finally, the concrete crack at the downstream slope change extends about 5.2 m upstream and about 7.5 m downward, and the bedrock cracks at the dam heel extend downward 37.2 m. When the transcendence probability of aftershock PGA is high (the aftershock intensity is small), under the main-aftershock (condition a∼c), the damage areas of the dam body and dam foundation have no obvious change compared with those under a single main shock. With the decrease in aftershock PGA transcendence probability (condition d∼g), that is, the aftershock intensity increases, the damage areas of the dam body and foundation have increased. Take condition g as an example, when the ground motion ends, the concrete crack at the downstream slope change extends about 10.2 m upstream and about 8.5 m downward and the bedrock cracks at the dam heel extend downward 46.5 m.

Dissipated Energy Characteristics of Gravity Dam under Main-Aftershock
Under the action of earthquakes, damage to concrete materials of dams and rock materials of dam foundations will accumulate in the form of dissipated energy. Damage evolution process of dam structure is closely related to the ·· Class III rock Class III rock Class II rock C15 C20 Figure 7: Finite element analysis model of gravity dam.  Shock and Vibration process of dissipated energy [25,26].
erefore, understanding the energy changes of the structure in the whole time history of main-aftershock sequence is of great significance to the damage of the cognitive structure under the action of main-aftershock sequence. Based on the plastic damage model, this paper analyzes and compares the dissipated energy characteristics of concrete gravity dams with different aftershock transcendence probability levels.

Dissipated Energy Equation.
e dynamic equilibrium equation under earthquake action is shown in the following equation [26]: where M is the mass matrix, C denotes the damping matrix, r denotes the vector of restoring forces, p denotes the vector of preseismic applied loads, and the vector € u t represents the absolute acceleration, which is the sum of the relative acceleration, € u, and the ground acceleration, € u g . Integrating the relative displacement u on both sides of equation (16), the dynamic equilibrium between energy components can be expressed as e three terms on the left side of equation (17) successively represent as the absolute kinetic energy (E K , which can be evaluated in a time step using the absolute velocity vector), the dissipated energy due to viscous damping (E C ), and the internal work done by nonlinear restoring forces e two terms on the right side of equation (17) successively represent as the seismic input energy (E WQ ) and the work done by preseismic applied loads (E WP ). e work done by the nonlinear restoring force can be expressed as follows: where d T is the damage value; σ r represents the nonlinear recovery stress; ε pl , respectively, represents the elastic and plastic strain; E E is expressed as recoverable strain energy; and E D and E P , respectively, represent the damage and plastic dissipation energy.

Dam Body and Foundation Whole Dissipated Energy
Characteristics. For the calculation results of seismic damage in Section 3, this article analyzes and compares two dissipated energy indicators: damage dissipated energy (ALLDMD) and plasticity dissipated energy (ALLPD), as shown in Figures 9 and 10. After the single main shock, the whole damage and plastic dissipated energy of concrete gravity dams are 16.76 and 41.19 kN·m, respectively. Under the main-aftershock, when the aftershock is small (condition a), the energy consumption of the concrete gravity dam has no obvious change compared with those under a single main shock. is shows that after the main shock, smaller aftershock cannot cause the secondary dissipated    Figure 11 shows the whole dissipated energy of the concrete gravity dam after the earthquake. e whole plastic dissipated energy is bigger than the damage dissipated energy of the dam. Figure 12 shows the fitting curve of dam dissipated energy growth rate varying with ∇PGA, where the transverse axis is the ratio (∇PGA) of aftershock peak acceleration to main shock peak acceleration and the vertical axis is the dissipated energy growth rate (the same below). e corresponding ∇PGA of the condition a∼g is 0.415, 0.555, 0.645, 0.705, 0.775, 0.865, and 0.990, respectively. It can be seen from Figure 12, with the increase in ∇PGA, the whole damage and plastic dissipated energy of the concrete gravity dam show different changes. When ∇PGA ≤ 0.44 (the transcendence probability of the aftershock PGA is 61.5%), the aftershock action can not cause further damage to the dam. When 0.44 < ∇PGA < 0.74, the whole damage Shock and Vibration dissipated energy growth rate of the dam is greater than the plastic dissipated energy growth rate. However, when 0.74 ≤ ∇PGA (the transcendence probability of the aftershock PGA is 22.3%), the whole plastic dissipated energy growth rate of the dam is greater than the damage dissipated energy growth rate. On the whole, the fitting curve of the dam damage dissipated energy growth rate caused by the aftershock is flatter than the plastic dissipated energy growth rate. In order to quantify the cumulative damage effect of aftershock on the dam, this paper takes the dissipated energy growth rate of 10% as the critical point that the aftershock cause large secondary damage to the dam. At this time, the ∇PGA corresponding to the fitted curve of the whole damage and plastic dissipated energy growth rate is 0.68 and 0.71, respectively.

Dissipated Energy Characteristics of Dam Body and
Foundation. Under different conditions, the dissipated energy of the dam body is shown Figure 13. Similar to the whole dissipated energy of the dam, when the aftershock intensity is small, the dissipated energy of the dam body has no obvious change compared with those under a single main shock. And, the plastic dissipated energy of the dam body under different conditions is greater than that of the damage dissipated energy. Figure 14 shows the fitting curve of the dissipated energy growth rate of the dam body varying with∇PGA. When ∇PGA ≥ 0.46, the whole damage and plastic dissipated energy of the dam body increase with the increase in ∇PGA, but the variation trend is not the same. When 0.46 < ∇PGA < 0.74, there is little difference between the two variations. When 0.74 ≤ ∇PGA, with the increase in ∇PGA, the increasing rate of plastic dissipated energy of dam body is obviously higher than that of damage dissipated energy of dam body, which indicates that strong aftershock can cause greater plastic strain of the dam body.
Under different conditions, the dissipated energy of the dam body is shown Figure 15. Similar to the dam body, the plastic dissipated energy of the dam foundation under different conditions is greater than the damage dissipated energy. Fitting curve of dissipated energy growth rate of dam foundation varying with ∇PGA is shown in Figure 16. Similar to the whole energy characteristics of the concrete gravity dam, when ∇PGA ≤ 0.44, the aftershock action did not cause further damage to the dam foundation. When 0.44 < ∇PGA < 0.74, the growth rate of the dam foundation plastic dissipated energy is greater than the damage dissipated energy growth rate. However, when 0.74 ≤ ∇PGA, the plastic dissipated energy growth rate of the dam foundation is greater than the damage dissipated energy growth rate. On the whole, the fitting curve of plastic dissipated energy growth rate of the dam foundation is steeper than that of damage dissipated energy growth rate.
Compare Figures 13 and 15. Under different conditions, the cumulative effect of aftershocks on the damage of dam foundation is more significant than that of dam body. For example, when ∇PGA � 0.9, the dam body damage and plastic dissipated energy increased by 10.3% and 22.7%, respectively, while the dam foundation damage and plastic dissipated energy increased by 35.2% and 50.4%, respectively.

Conclusion
On the basis of the previous research results and with the randomness of ground motion taken into consideration, a new method for main-aftershock sequence construction based on probability is proposed in this paper, which is intended for the seismic analysis under different main-aftershock ground motion sequences with different aftershock transcendence probabilities. e integral damage model of the dam system is established, with the theory of plastic damage mechanics applied for the analysis of the rock material. And then, the research on nonlinear dynamic response of concrete gravity dams with different aftershock transcendental probabilities is carried out. e main research results of this paper are as follows: (1) In this paper, the appropriate statistical relationship of the main shock and aftershock magnitude is selected, the attenuation relationship of aftershock ground motion parameters considering the nearfield saturation effect is introduced, the transcendental probability (P) is utilized to represent the statistical characteristics of the aftershock intensity, and a method for determining the ground motion input parameters of the main-aftershock sequence is eventually proposed. (2) Different from merely considering the damage characteristics of the dam body materials, this paper considers the damage of the rock material at the dam foundation under earthquake and it turns out that the cracks at the dam heel appear to expand downward in the depth direction, which may cause severe irreversible deformation of the dam foundation. erefore, it should be taken seriously in the seismic analysis of concrete gravity dams. (3) Under the action of the main-aftershock, the damage and plastic dissipated energy of the dam have shown different growth laws with the decrease in the transcendence probability (P) of the aftershock ∇PGA. It turns out that the cumulative effects of aftershock on the damage and plastic strain of the dam foundation are more significant than those of the dam body. In addition, when ∇PGA ≤ 0.44 (the aftershock transcendence probability is about 61.5%), the aftershock action will not cause further damage to the dam, whereas when ∇PGA > 0.68 (the aftershock transcendence probability is about 27.1%), the aftershock action will cause large secondary damage to the dam structure.
It should be noted that this article only considers bedrock materials as the initial state. Actually, after consolidation and grouting of the dam foundation, its elastic modulus and tensile strength will be improved upon. In addition, though a series of values are selected for the transcendental probability of the relative peak acceleration of the main-aftershock in this paper, how to choose appropriate transcendental probability in the design still needs to be further studied on.

Data Availability
e data used to support the research of this article are listed within the article or marked with references in the corresponding position.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper. Shock and Vibration 11