Stability of Reinforced Retaining Wall under Seismic Loads

: Based on the horizontal slice method (HSM) and assuming a log spiral slip surface, a method to analyze the stability of a reinforced retaining wall under seismic loads was established in this study by calculating the tensile force of the reinforcement. A parametric study was conducted on the normalized tensile force of the reinforcement, and it was observed that the normalized tensile force tends to increase with acceleration of the seismic load and the height of the backﬁll. Moreover, it also increases with soil unit weight, while it decreases with increased friction angle of the backﬁll soil, and the inﬂuence of soil cohesion on the normalized tensile force is not signiﬁcant. The HSM method is proved to be suitable for analyzing the tensile force of reinforcement in retaining walls under seismic loads.


Introduction
The reinforced retaining wall is a composite structure widely adopted in seismic regions due to its excellent deformation characteristics under seismic loads, such as earthquakes, owing to the friction between the reinforcement and the soil [1][2][3][4][5][6][7][8][9][10]. Okabe [11] proposed the Mononobe-Okabe method, which has been applied in analyzing the stability of reinforced retaining walls under dynamic loads [12]. Zeng and Steedman [13] studied the propagation of the shear wave in the backfill of a retaining wall by the pseudodynamic method and validated the analysis by centrifuge tests. Based on the pseudostatic method, Ling and Leshchinsky [14] and Kwak [15] investigated the vertical seismic response of reinforced soil structures under earthquake conditions. Shahgoli et al. [16] adopted the horizontal slice method (HSM) to analyze the stability of reinforced soil structures subjected to seismic loads. Choudhury and Singh [17] used the pseudostatic method to calculate the active earth pressure behind a retaining wall under seismic load. The slip surface of the reinforced retaining wall was also detected by vibration table tests [18,19], and it was found that the slip surface was like a log spiral. However, the mechanism of collaboration between reinforcement and soil remains unclear. Calculation of the tensile force of the reinforcement is still far from being fully understood and established [20,21]. Based on HSM and log spiral slip surface [16,22,23], the collaborative work of reinforcement and soils is considered in this study to establish a method of analyzing the stability of a reinforced retaining wall under seismic loads. Calculation of the tensile force of the reinforcement is carried out by considering the influence of soil parameters and loading conditions [24,25]. The slice method is proposed to determine the critical failure angle of the backfill wedge under complex conditions, and an iterative calculation method is presented to determine the tension crack depth of the active earth pressure under seismic loads [26,27]. The axial stress of the anchor is greatly increased by seismic excitation, and the increment is mainly induced within the excitation period of great acceleration amplitude, evaluation of the effect of earthquake frequency content on seismic behavior of cantilever retaining wall including soil-structure interaction [28,29]. By taking the failure surface as a plane, the analytical expression for the critical angle of the failure surface, earthquake-induced displacements of gravity retaining walls and anchor-reinforced slopes are also derived [30,31].
A critical review of literature shows that a very few studies have been reported for the stability analysis of a reinforced retaining wall. The present paper extends the latter work by developing a number of formulations, including a pseudostatic method of analysis. The formulations are described and compared with other published methods used for the seismic analysis of reinforced soil structures. However, no study seems to be available to quantify the effect of seismic acceleration coefficients and soil and reinforcing structure on the stability of reinforced soil walls. In view of this observation, the present study is directed to conduct a seismic stability analysis of the reinforced retaining wall using the HSM and considering the pseudo static forces. The limit equilibrium equation of a reinforced retaining wall under seismic loads is established, and the formula for calculating the tensile force in reinforcement is proposed to address the stability of reinforced retaining wall. In the end, the parametric study has been carried out to investigate the effect of various parameters like reinforcing structure, angle of internal friction and cohesion of soil, horizontal seismic loading.

Basic Assumptions of HSM
In the HSM approach, the sliding wedge is divided into horizontal slices with rigid-plastic behavior. The seismic force is considered as a pseudostatic force acting at the center of gravity of each slice [23]. The geotechnical interaction between the reinforcing structure and the surrounding soil is ignored. The assumptions are as follows: (1) The log spiral slip surface line of the reinforced retaining wall passes through the toe of the wall.
(2) The sliding body of the retaining wall is divided into several horizontal slices, the number of which is equal to that of reinforcement layers, and each slice contains reinforcement [32]. (3) The horizontal shear force between reinforcements is ignored. (4) The safety factor is assumed to be equal for various slices [23]. (5) Both horizontal and vertical seismic actions are applied to the center of gravity of the slice. (6) The analysis is done by the limited equilibrium method. (7) The anchoring strength of the reinforcement and the resistance of the reinforcement against pulling out are not taken into account.

Formulation of the HSM
HSM was proposed by Lo and Xu [22] and was aimed at overcoming the difficulties in analyzing reinforced soil structures by the vertical slice method (Figure 1). In Figure 1, H is the height of the retaining wall (in m); S y,j denotes the vertical distance between the reinforcement of layer j and layer j+1 (in m), and it is usually treated as a certain value in the design; S t,i and N N , i are the shear force and normal force acting on the ith slice, respectively; l ai is the length of reinforcement in the slip zone behind the wall; and l ei refers to the length of reinforcement in the anchorage zone. Figure 2 zooms in on the ith slice taken from Figure 1, which contains one reinforcement. In the figure, F N,I and F N,i+1 represent the normal force acting on the top and bottom of the ith slice, respectively; T i refers to the tensile force of the ith reinforcement [33]; Q hi and Q vi denote horizontal and vertical seismic force, respectively, acting on the top of the slice; α i is the inclination angle of slip surface of the slice in the horizontal direction; and W i is the weight of the ith slice. It should be noted that the slip surface is assumed to be a log spiral shape passing through the toe of the wall [34]. The sliding wedge is assumed to be in a limited equilibrium state, while both horizontal and vertical seismic loads are applied on the center of gravity of the slice. Then the safety for each slice is: where τ f and τ r refer to shear strength and shear stress on the sliding surface, respectively.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 13 assumed to be in a limited equilibrium state, while both horizontal and vertical seismic loads are applied on the center of gravity of the slice. Then the safety for each slice is: where f  and r  refer to shear strength and shear stress on the sliding surface, respectively.
Nomenclature could be found in Appendix A.

Seismic Loads
Under seismic loads, the shear wave velocity s V and compression wave velocity p V are characterized as: where G ,  , and  are shear modulus, density, and Poisson's ratio of soils, respectively. In this study, we assumed  is 0.3 for soils, therefore: The period for the horizontal seismic load T is: Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 13 assumed to be in a limited equilibrium state, while both horizontal and vertical seismic loads are applied on the center of gravity of the slice. Then the safety for each slice is: where f  and r  refer to shear strength and shear stress on the sliding surface, respectively.
Nomenclature could be found in Appendix A.

Seismic Loads
Under seismic loads, the shear wave velocity s V and compression wave velocity p V are characterized as: where G ,  , and  are shear modulus, density, and Poisson's ratio of soils, respectively. In this study, we assumed  is 0.3 for soils, therefore: The period for the horizontal seismic load T is:

Seismic Loads
Under seismic loads, the shear wave velocity V s and compression wave velocity V p are characterized as: where G, ρ, and υ are shear modulus, density, and Poisson's ratio of soils, respectively. In this study, we assumed υ is 0.3 for soils, therefore: The period for the horizontal seismic load T is: α h and α v are horizontal and vertical seismic acceleration, respectively. We assumed that the direction of propagation of seismic waves is vertical. Then the acceleration of z deep location from the surface of the retaining wall at moment t is: The body mass of the element in z with m depth and dz height m i (z) is: while the self-weight of the ith slice is: where γ is the unit weight of backfill soil, and l i and l i+1 are the length of the top and bottom side of the ith slice, respectively. Horizontal inertia force on the ith element is: Then the total horizontal seismic load on the slice Q h is: The vertical inertia force on the ith element is: Then the total horizontal seismic load on the retaining wall Q v is where λ is the vertical wavelength of the shear wave (λ = T × V S ), and η denotes the vertical wavelength of the compression wave ( . By considering q hi and q vi on the critical direction of the element, the limits of Q h and Q v according to Ling and Leshchinsky [14] are, respectively: where k h = α h /g and k v = α v /gX. When the ith slice is in limited equilibrium state, the vertical and horizontal stress applied on the slice are F y =0 and F x =0 [21,24], then: Appl. Sci. 2019, 9, 2175 5 of 13

Tensile Force of Reinforcement
According to the definition of safety factor [21], the shear force at the base of the slice is calculated by Equation (18): where c and ϕ are cohesion and friction angle of soil, respectively. The following formula can be obtained by compiling Equations (16) and (18): Then Equation (19) can be obtained by compiling Equations (17)-(19): After that, the tensile force of reinforcement is: For each layer of reinforcement, the tensile force is: where S y,j is the vertical distance between reinforcements and h j denotes the vertical distance between the jth reinforcement and the top of the wall. The normalized tensile force K for all reinforcements is [21]: For the stability analysis of a reinforced retaining wall under seismic loads [35], the maximum tensile force of each reinforcement needs to be calculated by Equation (22). This is also to determine the most critical log spiral slip surface, in which α i is obtained. This problem is the optimal solution of the linear programming problem. Matlab was used to get the interactive solution of n i=1 T i and α i . Figure 3 shows the geometry of the reinforced retaining wall as a case study. It can be seen that the wall is 6 m in height (H). The retaining wall is made of 6 pieces of reinforced concrete panel, each of which has a height of 1 m. The toe of the wall is made of plain concrete, with a width of 0.6 m and buried depth of 0.6 m. The reinforcement is made of metal, with a cross-section of 3 × 0.2 cm. Peak seismic acceleration α max is assumed to be 0.2 g. Factors affecting the stability of the retaining wall investigated are seismic acceleration coefficients (k h and k v ) and soil parameters (c, γ, and ϕ); parameters for the study can be seen in Table 1.

Influence of Seismic Acceleration Coefficients
Variation of the normalized tensile force K with friction angle was also evaluated for various seismic acceleration coefficients.  is fixed as 20 kN/m 3 , while is c is treated as 0. The horizontal seismic acceleration coefficient kh [36] is calculated by: while kh and kv need to fulfill the relationship of:  The effect of soil's relative density on the friction angle has not been taken into account.

Influence of Seismic Acceleration Coefficients
Variation of the normalized tensile force K with friction angle was also evaluated for various seismic acceleration coefficients. γ is fixed as 20 kN/m 3 , while is c is treated as 0. The horizontal seismic acceleration coefficient k h [36] is calculated by: while k h and k v need to fulfill the relationship of: Figures 4 and 5 show the influence of seismic acceleration coefficients on the normalized tensile force, while Figure 6 illustrates the influence of seismic acceleration coefficients on the inclination angle. It can be identified that normalized tensile force increases with the increase of seismic acceleration coefficient, especially when k h is increased from 0.1 to 0.2. Comparing Figure 4; Figure 5 shows that k v also contributes to K, but it only has limited effect. This implies that the horizontal seismic load is the dominant force for a reinforced retaining wall. It also shows that when the friction angle increases, the normalized tensile force will decrease. This means that the use of a well-graded backfill material with a large friction angle will be beneficial to the stability of the reinforced retaining wall under seismic loads [37]. Figure 6 shows that the slip surface is almost linear from the toe to the middle part of the wall, since the variation of the inclination angle for each slice is very small, while from the middle to the top surface of the wall, the slip surface is a vertical line parallel to the wall. It also shows that the slip surface of the retaining wall moves toward the wall side with the increasing seismic load. the friction angle increases, the normalized tensile force will decrease. This means that the use of a well-graded backfill material with a large friction angle will be beneficial to the stability of the reinforced retaining wall under seismic loads [37]. Figure 6 shows that the slip surface is almost linear from the toe to the middle part of the wall, since the variation of the inclination angle for each slice is very small, while from the middle to the top surface of the wall, the slip surface is a vertical line parallel to the wall. It also shows that the slip surface of the retaining wall moves toward the wall side with the increasing seismic load.    the friction angle increases, the normalized tensile force will decrease. This means that the use of a well-graded backfill material with a large friction angle will be beneficial to the stability of the reinforced retaining wall under seismic loads [37]. Figure 6 shows that the slip surface is almost linear from the toe to the middle part of the wall, since the variation of the inclination angle for each slice is very small, while from the middle to the top surface of the wall, the slip surface is a vertical line parallel to the wall. It also shows that the slip surface of the retaining wall moves toward the wall side with the increasing seismic load.   It shows that with increased soil cohesion, the normalized tensile force tends to increase, and its increment is obvious. It can be identified from Figure 8 that the slip surface of the retaining wall moves toward the wall side with increasing cohesion. To investigate the effect of soil unit weight, soil cohesion and friction angle are fixed at 10 kPa and 33°, respectively, k = 0.5, and k = 0.2. It can be seen in Figure 9 that bigger unit weight leads to higher normalized tensile force. This implies that better stability of the retaining wall will be achieved by using lightweight backfill materials. Figure 10 illustrates that the slip surface of the retaining wall moves toward the wall side with increasing soil unit weight. Figure 11 shows the influence of soil friction angle on normalized tensile force (k = 0.5, k = 0.2,  = 20 kN/m 3 , c = 10 kPa), and that normalized tensile force decreases with the increased soil friction angle. Therefore, increasing the soil friction angle (such as by using some coarse materials) will contribute to the stability of the retaining wall. As demonstrated in Figure 12, the slip surface of the retaining wall moves toward the wall side with increments of friction angle.

Influence of Soil Parameters
1.4 Figure 6. Inclination angle of slip surface and horizontal direction for various k h . Figure 7 shows the influence of cohesion of soil on normalized tensile force (k = 0.5 k h , γ = 20 kN/m 3 , ϕ = 33 • ), while Figure 8 illustrates the influence of cohesion of soil on the inclination angle. It shows that with increased soil cohesion, the normalized tensile force tends to increase, and its increment is obvious. It can be identified from Figure 8 that the slip surface of the retaining wall moves toward the wall side with increasing cohesion. To investigate the effect of soil unit weight, soil cohesion and friction angle are fixed at 10 kPa and 33 • , respectively, k = 0.5, and k = 0.2. It can be seen in Figure 9 that bigger unit weight leads to higher normalized tensile force. This implies that better stability of the retaining wall will be achieved by using lightweight backfill materials. Figure 10 illustrates that the slip surface of the retaining wall moves toward the wall side with increasing soil unit weight. Figure 11 Appl. Sci. 2019, 9, 2175 8 of 13 shows the influence of soil friction angle on normalized tensile force (k = 0.5, k = 0.2, γ = 20 kN/m 3 , c = 10 kPa), and that normalized tensile force decreases with the increased soil friction angle. Therefore, increasing the soil friction angle (such as by using some coarse materials) will contribute to the stability of the retaining wall. As demonstrated in Figure 12, the slip surface of the retaining wall moves toward the wall side with increments of friction angle.

Influence of Soil Parameters
It shows that with increased soil cohesion, the normalized tensile force tends to increase, and its increment is obvious. It can be identified from Figure 8 that the slip surface of the retaining wall moves toward the wall side with increasing cohesion. To investigate the effect of soil unit weight, soil cohesion and friction angle are fixed at 10 kPa and 33°, respectively, k = 0.5, and k = 0.2. It can be seen in Figure 9 that bigger unit weight leads to higher normalized tensile force. This implies that better stability of the retaining wall will be achieved by using lightweight backfill materials. Figure 10 illustrates that the slip surface of the retaining wall moves toward the wall side with increasing soil unit weight. Figure 11 shows the influence of soil friction angle on normalized tensile force (k = 0.5, k = 0.2,  = 20 kN/m 3 , c = 10 kPa), and that normalized tensile force decreases with the increased soil friction angle. Therefore, increasing the soil friction angle (such as by using some coarse materials) will contribute to the stability of the retaining wall. As demonstrated in Figure 12, the slip surface of the retaining wall moves toward the wall side with increments of friction angle.

Influence of the Reinforcing Structure
Vertical distance between reinforcements (Sy,j) is the main parameter influencing the tensile force of the reinforcement. It can be identified from Equation (21) that tensile force has a positive relationship with the vertical distance between reinforcements. In other words, the tensile force will

Influence of the Reinforcing Structure
Vertical distance between reinforcements (Sy,j) is the main parameter influencing the tensile force of the reinforcement. It can be identified from Equation (21) that tensile force has a positive relationship with the vertical distance between reinforcements. In other words, the tensile force will increase with increasing vertical distance between reinforcements. Thus, reducing the vertical

Influence of the Reinforcing Structure
Vertical distance between reinforcements (S y,j ) is the main parameter influencing the tensile force of the reinforcement. It can be identified from Equation (22) that tensile force has a positive relationship with the vertical distance between reinforcements. In other words, the tensile force will increase with increasing vertical distance between reinforcements. Thus, reducing the vertical distance between reinforcements is beneficial to reduce the tensile force, which is favorable for wall stability. Equation (22) also indicates that the reinforcing structure will not have a direct impact on the slip surface of the retaining wall.

Discussion
Based on a parametric analysis of the tensile force of reinforcement in a retaining wall, it can be seen that the horizontal seismic load, the vertical distance between reinforcements, and the soil friction angle are the three dominant factors influencing the tensile force. The findings are in agreement with Saeed and Ali [38] and Nouri et al. [23,32,39]. Vertical seismic load and soil cohesion have an insignificant impact on the tensile force of reinforcement, which is also close to the conclusions drawn by Nimbalkar et al. [34], Syed et al. [40], and Chandaluri et al. [41]. In summary, the HSM method adopted in this study could lead to similar results from a pseudostatic or experimental approach in previous research, proving that this method is suitable for analyzing the tensile force of reinforcement in retaining walls under seismic loads.

Conclusions
The horizontal slice method was adopted to analyze the stability of the reinforced retaining wall under seismic loads, while the slip surface was assumed to have a log spiral shape. Formulas for calculating the tensile force of reinforcement were established in this study. Parametric studies revealed that normalized tensile force increases with increased seismic acceleration coefficient and soil unit weight, while the tensile force of reinforcement will decrease with the increment of soil friction angle. Moreover, the influence of soil cohesion on normalized tensile force is not significant. Regarding the log spiral slip surface, it looks similar in various cases. It is almost linear from the toe to the middle part of the wall, since the variation of the inclination angle for each slice is very small, while from the middle to the top surface, it is a vertical line parallel to the wall. The slip surface of the retaining wall will move toward the wall side with increasing seismic load, soil cohesion, unit weight, or friction angle. The findings in this paper can provide theoretical basis and guidance for the design of reinforced retaining wall, and they could also strengthen the safety management of a reinforced retaining wall by considering earthquake conditions.