Reliability Analysis of High Concrete-Face Rockﬁll Dams and Study of Seismic Performance of Earthquake-Resistant Measures Based on Stochastic Dynamic Analysis

: The randomness of earthquake excitation has a signiﬁcant impact on the seismic performance of high earth-rock dams. In this paper, the seismic performance of geosynthetic-reinforced soil structures (GRSS) of high concrete face rockﬁll dams (CFRDs) is evaluated from the stochastic perspective. Multiple groups of seismic ground motions are generated based on spectral expression-random function non-stationary model. Taking Gushui CFRD as an example, this study calculates the failure probability of each damage level of non-reinforce slopes and reinforce slopes based on generalized probability density evolution method (GPDEM) and reliability analysis is presented though multiple evaluation indicators. The result shows that GRSS can reduce the mild damage of CFRDs during earthquake and restrain the moderate and severe damage. The inﬂuence of vertical spacing and length of GRSS on the seismic performance is obtained, which provides a reference for the seismic design and risk analysis of CFRDs.


Introduction
In recent years, the high earth-rock dams under construction or proposed in China are mainly distributed in the western regions, where earthquakes occur frequently. Once a dam breaks, it could cause immeasurable losses. According to domestic and international experimental results [1,2] and existing earthquake experiences [3,4] with earth-rock dams, the safety of earth-rock dams is closely related to the stability of the dam slope downstream.
Geosynthetic-reinforced soil structures (GRSS) perform well during strong earthquakes, in comparison with other earthquake-resistant measures, in maintaining the stability of slope rockfill [5,6]; therefore, the study of GRSS should receive greater attention. Li et al. [7] used the Newmark sliding block displacement method to evaluate the effect of reinforcement technology on dam crest rockfill, and Noorzad and Omidvar [8] performed a parametric analysis to study the effect of reinforcements on the seismic behavior of reinforced dams. The results showed that reinforcement measures can reduce dam settlement and maximum shear strain, but increase the maximum horizontal peak acceleration. Zhu et al. [9] analyzed the influence of various parameters of GRSS on seismic performance, based on the Fast Lagrangian Analysis of Continuum (FLAC) method. Yang [10] evaluated the seismic stability of reinforced earth-rock dams, based on the upper-bound theorem of limit analysis, and studied the influence of geogrid length on seismic performance. Glovatsky et al. [11] developed the theoretical foundations for the modeling and design of test stands in the

Reliability Analysis Method Based on Stochastic Dynamics
The reliability analysis of engineering has made great progress in the field of earthquakeresistant engineering research in recent years [16,17]. The system of reliability analysis of earthquake-resistant measures has four main steps: (1) multiple groups of random ground motion processes are generated, based on the non-stationary ground motion model; (2) the finite element model and the input the generated groups of random ground motions are utilized for batch computing; (3) the probability density function (PDF) and the cumulative distribution function (CDF) are obtained for dam safety performance indices, combined with the generalized probability density evolution theory; and (4) based on these steps, a probability model is constructed to analyze the reliability of a dam before and after the application of earthquake-resistant measures, so seismic performance may be evaluated.

Non-Stationary Ground Motion Model
In this paper, the spectral expression-random function non-stationary ground motion model [18] was established, based on the improved Clough-Penzien power spectrum model [19], and used for the generation of multiple groups of random ground motion.
The acceleration time series of the stochastic seismic ground motions are generated based on the spectral representation of the random function method of non-stationary stochastic processes [20]. This ground motion generation method has good applicability in the calculation considering the randomness of ground motions [21]. The random process of non-stationary ground motion acceleration with zero mean can be generated by the following formula: where ω k = k∆ω(∆ω = ω/N). {X k , Y k } (k = 1, 2, . . . , N) are the standard orthogonal random variables with an interval frequency of ∆ω = 0.15 rad/s, N is the number of the truncated items with N = 1600 here [18]. S .. X g is the bilateral evolutionary power spectral density function, and the expression is as follows [18]: Mathematics 2021, 9,3124 3 of 17 where A(t) is the intensity modulation function and calculated as: where c is the average time of peak ground acceleration (PGA) emergence, d is the shape control index of A(t). In this paper, c is taken as 4 s and d is taken as 2, according to [18]. In the evolutionary power spectral density function, the frequency modulation function can be determined by the following parameters: where ω 0 and ξ 0 are the initial angular frequency and the initial damping ratio of the site soil; a and b are parameters determined on the basis of the field classification and seismic design categories; and T is the duration of the ground motion acceleration time history, which differs according to different sites. In this study, the site type used was I 1 , and the parameter values ω 0 , ξ 0 , a, b, and T were 25 (rad/s), 0.45, 3.5, 0.3, and 15 (s), respectively, in accordance with the China Hydraulic Seismic Design Code (NB 35047-2015). The spectral parameters S 0 (t) in Equation (2) reflecting the spectral intensity can be expressed as: where a max is the mean value of PGA with a max = 0.340 g, which is the checking ground motion of Gushui CFRD; and γ is the equivalent peak factor, depended on the sort of seismic site, taken here to be 2.6 according to [18].
When ω = 0, the following equation should be satisfied: In Equation (1), the variables {X k , Y k } (k= 1, 2, . . . , N) are standard orthogonal random variables, which are uniquely determined by the orthogonal basis function constructed based on the idea of random function through mapping, and meet the following basic conditions: where E [•] represents mathematical expectation and δ jk is the Kronecker delta. The method of constructing standard orthogonal random variables is as follows: Suppose X n and Y n (n= 1, 2, . . . , N) are two independent random variables respectively Θ 1 and Θ 2 , then the random function can be recorded as: where cas(x) = cos(x) + sin(x) is the Hartley orthogonal basis function [22], basic random variables Θ 1 and Θ 2 are distributed uniformly and independent in the interval [0, 2π], which can usually be obtained by number-theoretic method. After certain deterministic mapping, they become the standard orthogonal random variables required by the Equation (1). The mean square error of the above non-stationary ground motion acceleration process simulation can be expressed as: where ω u = N∆ω is the truncation frequency, and generally the mean square error of peak acceleration is limited to far less than 1.0 [20].

Generalized Probability Density Evolution Theory
The generalized probability density evolution method (GPDEM) was proposed by Li et al. [23]. Starting from random events, the theory combines the decouple system's physical equations, based on the principle of conservation of probability, to obtain the generalized probability density evolution equation and to establish the relationship between each physical quantity of interest in the project and in the dynamic system through constitutive relations and deformation coordination relations. Thus, the theory is expressed as a function of basic random variables, and then the probability density function of the studied physical quantity is solved by combining the initial conditions and the boundary conditions. In recent years, the theory has made good progress in the analysis of uncertainty reliability and in the application of large-scale nonlinear structures [24,25].
The motion equation of n degree of freedom system can be expressed as:
where M, C are the n × n order mass and damping matrix, G (·) is the linear or nonlinear restoring force vector, ..
X and X are the acceleration, velocity and displacement vectors of the structural response respectively, the Γ is n × r order excitation influence matrix and the F(Θ, t) is r order excitation vector. And n is the number of degrees of freedom of the system and r is the order of external excitation here.
For a general well posed dynamic system, the physical solution of Equation (12) exists, uniquely and continuously depends on the basic parameters, so the solution of Equation (12) can be expressed as: The speed process can be expressed as: The information of other physical quantities Z = (Z 1 , . . . , Z m ) T in practical engineering can also be expressed as a function of basic random variables by establishing relations with .. X and . X through constitutive relations and deformation coordination relations: The time change rate (speed) can be expressed as: .
Since Equation (16) itself can be considered as a random dynamic process, the randomness comes entirely from Θ. The extended random process (Z t , Θ) is described according to the random event of probability conservation in the whole evolution process. The gener-alized probability density evolution equation can be obtained, considering the arbitrariness of Ω Θ [26]: where p ZΘ (z, θ, t) refers to the joint PDF of (Z, Θ), in which the source random factors are completely described by Θ. Z refers to the physical quantity studied. The augmented system composed of (Z, Θ) is a conservative probability system, which follows the law of probability conservation. It is worth pointing that the dimension m of this equation is the number of physical quantities studied. When only a certain response physical quantity is considered, the equation can further degenerate into a one-dimensional partial differential equation: Equation (18) has only partial differential for z and t, while θ is in the form of a parametric equation. Therefore, a series of deterministic values can be obtained, that is, for a given Θ = θ q (q = 1, 2, . . . , n sel ), where q represents different ground motion processes, and n sel = 89 in this paper. Derivative (velocity) • Z j (θ q , t m )(j = 1, 2, . . . , m) of the required physical quantity and θ q can be obtained by solving the physical equation. When the random parameters are determined by selecting points in the probability space, the partial differential equation of the random dynamic system is transformed into a set of deterministic dynamic equations. The engineering structure can be solved by various numerical simulation methods such as finite element method and finite difference method, and Equation (18) becomes the following series of equations: The initial condition of Equation (19) is: The boundary conditions of Equation (19) are: ∂p ZΘ (z, θq, t) z j →±∞ = 0 , j = 1, 2, . . . , m.
The discrete numerical solutions p ZΘ (z, θq, t) can be obtained by bringing in the initial conditions and boundary conditions. The solutions p Z (z, t) can be obtained by accumulating all the above discrete numerical solutions, and the result when m = 1 is as follows:

Reliability Calculation Based on GPDEM
The dynamic reliability of engineering structures usually includes two kinds of problems: first exceedance probability failure and cumulative damage failure reliability, which can be obtained by constructing a virtual random process and solving the corresponding generalized probability density evolution equation [27].
Taking a random process X(Θ, t) (random ground motions in this study) as an example, its extreme value can be expressed as follows [28]: Mathematics 2021, 9, 3124 6 of 17 The extreme value or cumulative value depends on the source random vector Θ and the random variables of the time period [0, T], so a virtual process can be constructed: Obviously, the condition Equation (25) can be met: For Equation (24), the derivative of τ is: Since the randomness of the virtual random process Z X (τ) comes entirely from Θ, (Z X (τ), Θ) constitutes a probabilistic conservative system. According to the GPDEM, the joint PDF of (Z X (τ), Θ), that is p ZΘ (z, θ, t), satisfies the following generalized probability density evolution equation: The equation is solved to obtain the cumulative distribution curve of the studied physical quantity based on the solution method in Section 2.2, so as to obtain the transcendence probability corresponding to each value.

Calculation Indexes of Reliability Analysis of Concrete Face Rockfill Dam
At present, the finite element dynamic time history analysis method is used for the stability analysis of the dam body, and the indexes include cumulative time of Fs < 1.0 and cumulative slip. In this paper, the two indexes are used to judge the damage grade of the dam body.

Overall Stability Index
(1) Safety factor The safety factor is the ratio of the maximum shear strength provided by the soil on the potential sliding surface to the actual shear stress generated by the external load. The preearthquake stress of the dam and the instantaneous dynamic stress during the earthquake are calculated by the finite element method. The stability of the dam is calculated based to the static and dynamic superposition results of the element based on Newmark method [29,30], and the safety factor is calculated with the following formula: where c i and ϕ i are the cohesion and internal friction angle of the soil in the ith cell respectively; l i is the length of the slip arc through the ith cell; σ i and τ i are the normal stress and tangential stress on the slip arc surface of the ith cell, which are obtained by superposition of pre-earthquake stress and dynamic stress.
(2) Cumulative time of Fs < 1.0 The time cumulative time of Fs < 1.0 in the whole earthquake process is obtained by determining and judging the minimum safety factor at each time and cumulative addition. In this paper, the criteria for classifying the damage level according to cumulative time is taken as: it is considered that slight damage occurs when the cumulative time is greater than 0 s. The time of 0.5 s is the moderate failure limit, 1.5 s is the severe failure limit, and 2 s is the local non dam break limit according to [31].

Local Sliding Failure Index
In this paper, the cumulative slip is used as the index to evaluate the local slip failure degree of the dam body. Based on the existing research [32,33] and the China Hydraulic Seismic Design Code (NB 35047-2015), the evaluation criteria of the cumulative slip index used in this paper are as follows: when the slip begins, it is judged as mild failure; 20 cm is the moderate failure limit, 100 cm is the severe failure limit, and 150 cm is the local non dam break limit.
For any sliding arc, the sliding angular velocity of the slider around the center of the circle can be calculated by Equation (29) [30]: where, I is the moment of inertia of the sliding body; α(t) is the sliding angular velocity of the sliding body after instantaneous instability; M is the rotational moment acting on the sliding body, and R is the slip arc radius.
When an instantaneous slip occurs in a slip arc at a certain time, the slip amount of the slip arc is: Multiple instantaneous sliding may occur in the whole time period, and the cumulative sliding amount is: Maximum slip of dam slope is the maximum cumulative slip of all possible slip arcs:

Reliability Analysis
This paper mainly selects two characteristic physical quantities: safety factor cumulative time of Fs < 1.0 and cumulative slip to determine the stability of the dam body, and carries out reliability analysis by calculating the failure probability of different grades of the dam body. The state function of the structure can be expressed as: where U is the structural state function, which is used to represent the safe state of the structure. At that time U = 0, it means that the structure reaches the failure limit state, when U < 0, it means that the structure is damaged, when U > 0, it means that the structure is safe; R indicates the comprehensive resistance of the structure (i.e., the allowable failure state of the structure in this paper), and S indicates the bearing effect of the structure (i.e., the calculated maximum value). If expressed as a function, R and S can be comprehensively expressed as functions of basic random variables (such as strength parameters, seismic load, displacement response, etc.). The reliability analysis of the structures (i.e., the calculation formula of failure probability) is following: where P f is failure probability. In this paper, it represents the probability that the calculated values of indicators exceed the allowable values for the specified failure level.

Model Establishment
This paper has used the Gushui CFRD as an example [34]. The elevation of the dam crest is 2287 m; the elevation of the dam bottom is 2042 m; the dam height is 242 m (the wave wall is not considered); the upstream dam slope is 1:1.5; the downstream dam slope above the path is 1:1.6; and the slope below the path is 1:1.5. The width of the dam crest is 20 m and the length of the dam crest is 437 m. The cushion zone and the transition zone are set under the face slab. The width of the cushion zone is 4 m, the total thickness of the transition zone is 8 m, and the face slab thickness is 0.4 m~1.24 m. The face slab was poured in three phases to 2167 m, 2235 m, and 2285 m, respectively, and the water storage was stored in two phases to 2267 m (the designed normal water level). The model was densified at the upper part of the dam, where the vertical spacing was set at 1m, considering the dam size, the research purpose, and existing studies [35,36], so the number of model nodes was 8432 and the number of elements was 8277 ( Figure 1) before the application of GRSS for the dam slope. The setting of earthquake-resistant measures is shown in Figure 2. The dam element was simulated by the quadrilateral isoparametric element, which is a continuous medium block isoparametric element. The non-thickness Goodman element was used for the indirect contact surface between the panel and the cushion zone, based on the assumption that the two contact surfaces were connected by countless tiny tangential and normal springs. The contact surface and adjacent contact surface elements only had force connected at the node. The above two element types are often used in the field of earth-rock dams [37], and the accuracy meets engineering needs. After filling and water storage, static, dynamic, and two-dimensional stability calculations were carried out. The hydrodynamic pressure on the panel was simulated by the added-mass method [38]. mic load, displacement response, etc.).
The reliability analysis of the structures (i.e., the calculation formula of failure probability) is following: where f P is failure probability. In this paper, it represents the probability that the calculated values of indicators exceed the allowable values for the specified failure level.

Model Establishment
This paper has used the Gushui CFRD as an example [34]. The elevation of the dam crest is 2287 m; the elevation of the dam bottom is 2042 m; the dam height is 242 m (the wave wall is not considered); the upstream dam slope is 1:1.5; the downstream dam slope above the path is 1:1.6; and the slope below the path is 1:1.5. The width of the dam crest is 20 m and the length of the dam crest is 437 m. The cushion zone and the transition zone are set under the face slab. The width of the cushion zone is 4 m, the total thickness of the transition zone is 8 m, and the face slab thickness is 0.4 m~1.24 m. The face slab was poured in three phases to 2167 m, 2235 m, and 2285 m, respectively, and the water storage was stored in two phases to 2267 m (the designed normal water level). The model was densified at the upper part of the dam, where the vertical spacing was set at 1m, considering the dam size, the research purpose, and existing studies [35,36], so the number of model nodes was 8432 and the number of elements was 8277 ( Figure 1) before the application of GRSS for the dam slope. The setting of earthquake-resistant measures is shown in Figure  2. The dam element was simulated by the quadrilateral isoparametric element, which is a continuous medium block isoparametric element. The non-thickness Goodman element was used for the indirect contact surface between the panel and the cushion zone, based on the assumption that the two contact surfaces were connected by countless tiny tangential and normal springs. The contact surface and adjacent contact surface elements only had force connected at the node. The above two element types are often used in the field of earth-rock dams [37], and the accuracy meets engineering needs. After filling and water storage, static, dynamic, and two-dimensional stability calculations were carried out. The hydrodynamic pressure on the panel was simulated by the added-mass method [38].

Constitutive Model and Material Parameters
In this paper, the static calculation of rockfill is simulated by Duncan-Chang E-B model [39]. Based on Duncan Chang E-μ, volume modulus B is used instead of Poisson's ratio t ν as the calculation parameter in Duncan-Chang E-B model, solving the problem which is quite different from the actual situation in the hyperbolic assumption. The dynamic calculation is analyzed by equivalent linear analysis method based on equivalent visco-elastic model [40]. In this paper, the static, dynamic and stability calculation are based on GEODYNA, which has good performance in the calculation of large-scale nonlinear structures [41]. Based on the dam site conditions and specifications, the static and

Constitutive Model and Material Parameters
In this paper, the static calculation of rockfill is simulated by Duncan-Chang E-B model [39]. Based on Duncan Chang E-µ, volume modulus B is used instead of Poisson's ratio ν t as the calculation parameter in Duncan-Chang E-B model, solving the problem which is quite different from the actual situation in the hyperbolic assumption. The dy-namic calculation is analyzed by equivalent linear analysis method based on equivalent visco-elastic model [40]. In this paper, the static, dynamic and stability calculation are based on GEODYNA, which has good performance in the calculation of large-scale nonlinear structures [41]. Based on the dam site conditions and specifications, the static and dynamic material parameters of dam rockfill are as show in Tables 1 and 2, comprehensively considering the existing studies [34,[42][43][44][45]. The density of GRSS is taken as 2.6 × 10 3 kg/m 3 , and the ultimate strength is 85 kN/m.

Ground Motions Input
Viscoelastic artificial boundary and equivalent node load are used for wave input [46,47] to realize ground motion input considering the radiation damping effect of infinite foundation and the influence of traveling wave effect. The horizontal PGA is 0.340 g. The vertical PGA is 2/3 of the horizontal acceleration.
This chapter uses the aforementioned spectral expression-random function nonstationary random ground motion model to generate the seismic input of 89 different action processes due to the difference of ground motions in the magnitude of ground motion frequency, the height of peak and valley value and the duration of fluctuation mode. When 89 sample acceleration time histories are generated, the error between the mean value of peak acceleration and the target value (zero) is 8.0%. The results show that the characteristics of the sample set are consistent with the target in the second-order numerical statistical sense [20]. Figure 3 shows the seismic acceleration time series information of the generated ground motion samples. Figure 3b,c show the mean and standard deviation of the generated ground motion samples. Figure 3e shows that the mean of response spectrum of 89 samples coincides with the specification spectrum in China Hydraulic Seismic Design Code (NB 35047-2015). It can be seen that the ground motions generated by the spectral expression-random function non-stationary ground motion model fit well with the target earthquake motion. of the generated ground motion samples. Figure 3b,c show the mean and standard deviation of the generated ground motion samples. Figure 3e shows that the mean of response spectrum of 89 samples coincides with the specification spectrum in China Hydraulic Seismic Design Code (NB 35047-2015). It can be seen that the ground motions generated by the spectral expression-random function non-stationary ground motion model fit well with the target earthquake motion.

Seismic Response Analysis
In order to study the influence of GRSS on the stability of dam slope during earthquakes, based on the existing research [9][10][11] and with reference to the common parameter settings in engineering [48], the vertical spacing of GRSS layers (Sv) was set at 4 m and the length (Lmax) was set at 50 m, for comparison with the original dam slope as a typical work condition. The random dynamic response of the dam slope was analyzed, and the seismic performance of the proposed measure was evaluated from the perspective of reliability.
Eighty-nine groups of safety factor history curves were obtained by calculating the dynamic response of the high CFRD under 89 groups of random ground motions. The safety factor probability density surface (Figures 4b and 5b) was solved by the finite difference method (FDM) in TVD format [26] to reflect the transmission information of the safety factor in space and time, based on the GPDEM. The PDF of the safety factor fluctuated in time and space, and the safety factor gradually concentrated with increases in time.

Seismic Response Analysis
In order to study the influence of GRSS on the stability of dam slope during earthquakes, based on the existing research [9][10][11] and with reference to the common parameter settings in engineering [48], the vertical spacing of GRSS layers (S v ) was set at 4 m and the length (L max ) was set at 50 m, for comparison with the original dam slope as a typical work condition. The random dynamic response of the dam slope was analyzed, and the seismic performance of the proposed measure was evaluated from the perspective of reliability.
Eighty-nine groups of safety factor history curves were obtained by calculating the dynamic response of the high CFRD under 89 groups of random ground motions. The safety factor probability density surface (Figures 4b and 5b) was solved by the finite difference method (FDM) in TVD format [26] to reflect the transmission information of the safety factor in space and time, based on the GPDEM. The PDF of the safety factor fluctuated in time and space, and the safety factor gradually concentrated with increases in time. The comparison between the two figures showed that GRSS have no significant impact on the shape of PDF as a dam safety factor. Mathematics 2021, 9, x FOR PEER REVIEW 11 of 18 The comparison between the two figures showed that GRSS have no significant impact on the shape of PDF as a dam safety factor.
(a) (b) (c)  Figure 6 shows the probability density information of the minimum safety factor before and after the application of earthquake-resistant measures, obtained by solving the equivalent extreme-value event. Before GRSS measures were taken, the average value of the minimum safety factor was 0.944 and the standard deviation was 0.131. After the measures were taken, the average value of the minimum safety factor increased to 1.003 and the standard deviation decreased to 0.105 (Figure 6a). The calculation results of the reinforced slopes showed that the minimum safety factor not only improved in value, but also reduced in dispersion, which indicated that GRSS can not only maintain the stability of a dam body when encountering earthquakes, but also improve the stability of the overall response of the dam body when encountering different ground motions. Similar conclusions were also reached in Figure 6b; that is, GRSS are conducive to improving the reliability of a dam under random ground motions. The comparison between the two figures showed that GRSS have no significant impact on the shape of PDF as a dam safety factor.
(a) (b) (c)  Figure 6 shows the probability density information of the minimum safety factor before and after the application of earthquake-resistant measures, obtained by solving the equivalent extreme-value event. Before GRSS measures were taken, the average value of the minimum safety factor was 0.944 and the standard deviation was 0.131. After the measures were taken, the average value of the minimum safety factor increased to 1.003 and the standard deviation decreased to 0.105 (Figure 6a). The calculation results of the reinforced slopes showed that the minimum safety factor not only improved in value, but also reduced in dispersion, which indicated that GRSS can not only maintain the stability of a dam body when encountering earthquakes, but also improve the stability of the overall response of the dam body when encountering different ground motions. Similar conclusions were also reached in Figure 6b; that is, GRSS are conducive to improving the reliability of a dam under random ground motions.  Figure 6 shows the probability density information of the minimum safety factor before and after the application of earthquake-resistant measures, obtained by solving the equivalent extreme-value event. Before GRSS measures were taken, the average value of the minimum safety factor was 0.944 and the standard deviation was 0.131. After the measures were taken, the average value of the minimum safety factor increased to 1.003 and the standard deviation decreased to 0.105 (Figure 6a). The calculation results of the reinforced slopes showed that the minimum safety factor not only improved in value, but also reduced in dispersion, which indicated that GRSS can not only maintain the stability of a dam body when encountering earthquakes, but also improve the stability of the overall response of the dam body when encountering different ground motions. Similar conclusions were also reached in Figure 6b; that is, GRSS are conducive to improving the reliability of a dam under random ground motions.

Overall Stability Analysis of Dam Body
In this section, the probability density function of the cumulative time of Fs < 1.0 was solved by constructing a virtual random process to analyze the stability of non-reinforced slopes and reinforced slopes, based on the 89 groups of random ground motions referred to in Section 3.3. Under different seismic ground motions, the maximum value of the cumulative time of the original dam body was 1.43 s, which was quite different from the minimum value of 0 s (Figure 7a); the seismic performance of the proposed measure in reducing cumulative time was inconsistent, indicating the need to analyze the seismic performance of earthquake-resistant measures from the stochastic perspective.

Overall Stability Analysis of Dam Body
In this section, the probability density function of the cumulative time of Fs < 1.0 was solved by constructing a virtual random process to analyze the stability of non-reinforced slopes and reinforced slopes, based on the 89 groups of random ground motions referred to in Section 3.3. Under different seismic ground motions, the maximum value of the cumulative time of the original dam body was 1.43 s, which was quite different from the minimum value of 0 s (Figure 7a); the seismic performance of the proposed measure in reducing cumulative time was inconsistent, indicating the need to analyze the seismic performance of earthquake-resistant measures from the stochastic perspective.
According to the aforementioned standard, the probability of a mild failure of the original dam slopes was large, up to more than 50%, while it was not easy to cause moderate or severe damage, as GRSS have an obvious effect in maintaining the stability of the dam slope (Table 3).   According to the aforementioned standard, the probability of a mild failure of the original dam slopes was large, up to more than 50%, while it was not easy to cause moderate or severe damage, as GRSS have an obvious effect in maintaining the stability of the dam slope (Table 3).  Figure 8 and Table 4 show the probability information of cumulative slippage of dam slopes under random ground motions. According to the calculated results, smaller parts of ground motions will generate local sliding of a CFRD under proposed seismic ground motions, and the failure probability of moderate or severe sliding of the dam slope is small. From the perspective of reliability, the proposed measure has an obvious inhibitory effect on the local sliding of the dam slope.  Table 4 show the probability information of cumulative slippage of dam slopes under random ground motions. According to the calculated results, smaller parts of ground motions will generate local sliding of a CFRD under proposed seismic ground motions, and the failure probability of moderate or severe sliding of the dam slope is small. From the perspective of reliability, the proposed measure has an obvious inhibitory effect on the local sliding of the dam slope.

Analysis on Influenceing Factors of GRSS Performance
In this section, failure probability was calculated by adjusting the length and vertical spacing of reinforcement based on the 89 groups of random ground motions as described in Section 3.3. The influence of reinforcement-setting on the seismic performance of GRSS was considered, to find the optimal GRSS setting from the perspective of reliability, which can be used to provide a reference for the design of earthquake-resistant measures, considering both project costs and effects.

Impact on Maintaining Overall Stability
In order to analyze the influence of the length and vertical spacing of reinforcement on the stability of CFRD slopes, this section kept other factors unchanged, while different characteristic lengths (Lmax = 40 m, 50 m, and 60 m) and characteristic vertical spacings (Sv = 2 m, 4 m, and 6 m) of GRSS were selected according to previous engineering experience and the research of previous literature [11,12]. The equivalent extreme values were calculated based on the GPDEM for the cumulative time of Fs < 1.0 (as shown in Figure 9). The results showed that with other factors unchanged, the reliability of CFRD increased with the increase in reinforcement length and decreased with the increase in reinforcement vertical spacing. The probability of mild failure of dam slope followed an order of 40.10%, 37.96%, and 37.80%, with the increases in the length of reinforcement. It can be seen that

Analysis on Influenceing Factors of GRSS Performance
In this section, failure probability was calculated by adjusting the length and vertical spacing of reinforcement based on the 89 groups of random ground motions as described in Section 3.3. The influence of reinforcement-setting on the seismic performance of GRSS was considered, to find the optimal GRSS setting from the perspective of reliability, which can be used to provide a reference for the design of earthquake-resistant measures, considering both project costs and effects.

Impact on Maintaining Overall Stability
In order to analyze the influence of the length and vertical spacing of reinforcement on the stability of CFRD slopes, this section kept other factors unchanged, while different characteristic lengths (L max = 40 m, 50 m, and 60 m) and characteristic vertical spacings (S v = 2 m, 4 m, and 6 m) of GRSS were selected according to previous engineering experience and the research of previous literature [11,12]. The equivalent extreme values were calculated based on the GPDEM for the cumulative time of Fs < 1.0 (as shown in Figure 9). The results showed that with other factors unchanged, the reliability of CFRD increased with the increase in reinforcement length and decreased with the increase in reinforcement vertical spacing. The probability of mild failure of dam slope followed an order of 40.10%, 37.96%, and 37.80%, with the increases in the length of reinforcement. It can be seen that when the length increased to more than 50 m, the stability reliability of dam slope did not increase significantly. With the increase in vertical spacing, the probability of mild failure of dam slope was 30.73%, 37.96%, and 39.47%, respectively. Other failure conditions were similar, proving that the stability of dam slope decreases significantly when reinforcement spacing is expanded from 2 m to 4 m ( Table 5). Table 5 compares and analyzes the consistency risk of various working conditions based on performance. Taking working condition 1 (with GRSS vertical spacing of 2 m and length of 50 m) as an example, under random ground motion with a GPA of 0.340 g, the probability of mild failure, moderate failure, and severe failure of the dam slope followed an order of 30.73%, 2.69%, and 0.06%, respectively. The probability basically did not reach the situation of local dam-break, which implied that the overall stability reliability of the dam slope was high when the vertical spacing of the geogrid reinforcement was 2 m and the length was 50 m. Even if damage occurred, it was mild and easy to repair, and had little economic or social impact under condition 1. The analysis of other working conditions is similar. These findings can provide a reference for the reliability analysis of GRSS in engineering.
increase significantly. With the increase in vertical spacing, the probability of mild failure of dam slope was 30.73%, 37.96%, and 39.47%, respectively. Other failure conditions were similar, proving that the stability of dam slope decreases significantly when reinforcement spacing is expanded from 2 m to 4 m ( Table 5). Table 5 compares and analyzes the consistency risk of various working conditions based on performance. Taking working condition 1 (with GRSS vertical spacing of 2 m and length of 50 m) as an example, under random ground motion with a GPA of 0.340 g, the probability of mild failure, moderate failure, and severe failure of the dam slope followed an order of 30.73%, 2.69%, and 0.06%, respectively. The probability basically did not reach the situation of local dam-break, which implied that the overall stability reliability of the dam slope was high when the vertical spacing of the geogrid reinforcement was 2 m and the length was 50 m. Even if damage occurred, it was mild and easy to repair, and had little economic or social impact under condition 1. The analysis of other working conditions is similar. These findings can provide a reference for the reliability analysis of GRSS in engineering.  Similarly, in order to study the influence of adjusting the length and vertical spacing of GRSS on the local sliding of concrete-face rockfill dams subjected to earthquakes, different lengths and vertical spacings of geogrid reinforcement were adopted in this section, and the cumulative sliding index was then calculated. Taking mild failure as an example, with an increase in reinforcement length, the probability of mild failure of dam slope followed an order of 25.41%, 23.63%, and 23.07% under the random ground motions with  Similarly, in order to study the influence of adjusting the length and vertical spacing of GRSS on the local sliding of concrete-face rockfill dams subjected to earthquakes, different lengths and vertical spacings of geogrid reinforcement were adopted in this section, and the cumulative sliding index was then calculated. Taking mild failure as an example, with an increase in reinforcement length, the probability of mild failure of dam slope followed an order of 25.41%, 23.63%, and 23.07% under the random ground motions with GPA of 0.340 g. With increases in reinforcement vertical spacing, the probability was 18.05%, 23.63%, and 26.52%, respectively. Moderate failure was similar, indicating that the local sliding failure probability of CFRD decreases with an increase in reinforcement length and increases with an increase in reinforcement spacing ( Figure 10 and Table 6). GPA of 0.340 g. With increases in reinforcement vertical spacing, the probability was 18.05%, 23.63%, and 26.52%, respectively. Moderate failure was similar, indicating that the local sliding failure probability of CFRD decreases with an increase in reinforcement length and increases with an increase in reinforcement spacing ( Figure 10 and Table 6).

Conclusions
In this paper, the reliability of the seismic performance of GRSS was evaluated from the perspective of probability and the influence of factors on performance. The conclusions were as follows: (1) A set of reliability analysis systems of high CFRDs was established by combining a non-stationary ground motion model, a generalized probability density evolution theory, and multi-evaluation indicators, providing a method for the evaluation of seismic performance of earthquake-resistant measures. (2) GRSS can not only maintain the stability of a dam body during an earthquake; it can also improve the stability of the overall response of a dam under different ground motions. (3) The dam slope of the Gushui CFRD is more susceptible to slight damage under random ground motions with a PGA of 0.340 g, with a certain probability of moderate damage and a lower probability of severe damage. The calculation results showed that the seismic performance of GRSS can inhibit mild damage of 20-30% of ground motion samples and avoid severe damage with a high probability. (4) With an increase in geogrid reinforcement length and a decrease in vertical spacing, the seismic performance in maintaining dam slope stability increases due to strong constraints. However, the increase in reinforcement length has little influence on the

Conclusions
In this paper, the reliability of the seismic performance of GRSS was evaluated from the perspective of probability and the influence of factors on performance. The conclusions were as follows: (1) A set of reliability analysis systems of high CFRDs was established by combining a non-stationary ground motion model, a generalized probability density evolution theory, and multi-evaluation indicators, providing a method for the evaluation of seismic performance of earthquake-resistant measures. (2) GRSS can not only maintain the stability of a dam body during an earthquake; it can also improve the stability of the overall response of a dam under different ground motions. (3) The dam slope of the Gushui CFRD is more susceptible to slight damage under random ground motions with a PGA of 0.340 g, with a certain probability of moderate damage and a lower probability of severe damage. The calculation results showed that the seismic performance of GRSS can inhibit mild damage of 20-30% of ground motion samples and avoid severe damage with a high probability. (4) With an increase in geogrid reinforcement length and a decrease in vertical spacing, the seismic performance in maintaining dam slope stability increases due to strong constraints. However, the increase in reinforcement length has little influence on the performance of GRSS. When reinforcement vertical spacing is reduced from 4 m to 2 m, the reduction in the earthquake-resistant effect of GRSS is relatively obvious. This is conducive to actual engineering design, considering economics and seismic performance.