1. Introduction
Because the facilities installed during a construction project are usually temporary, when applying an earth retaining method for the ground excavation in an urban setting, a stability analysis is generally conducted using the earth pressure and without considering the stability under a cyclic load or earthquake event. However, even when such a retaining method is applied, there may be special circumstances in which the structure is maintained for longer than the normal two-month period (the standard in the Korea [
1]). However, investigations into the stability under cyclic loading or earthquake events have typically been simplified or even omitted.
The occurrence of earthquakes has recently increased worldwide with the increased activity of the seismic belt. Consequently, the development of stable earth retaining walls is becoming increasingly important. In addition, unlike static loads, the behavior of the ground differs under dynamic load conditions such as cyclic loading. In recent years, the change of earth pressure or earth pressure coefficient have been studied mainly through laboratory model tests or numerical analysis of the grounds with retaining walls or earth retaining walls [
2,
3,
4,
5,
6]. A few studies have been conducted to assess the deformation of structures caused by dynamic loads such as railway loads. Therefore, dynamic behavioral analyses are required for cyclic railway loading.
Figure 1 shows a diagram outlining an earth retaining wall used to improve the stability through the installation of stabilizing piles (hereinafter referred to as the “back support”) at an angle and connecting them to the front wall (hereinafter referred to as the “front support”). This is called an inclined earth retaining (IER) method.
An IER structure increases the stability by preventing the earth pressure from acting on the front wall as a result of the connection between the front and back supports. The stability of an IER structure was confirmed through laboratory model tests and field experiments. First, prior to the laboratory model experiments, an IER model was installed in a plane strain model box placed in ground consisting of sand or clay. The lateral displacement according to the installed back support was then analyzed. The decrease in lateral displacement for the case in which the back support was placed in sandy ground was only 8% compared to the case without a back support [
7]. When placed in clay however, the lateral displacement decreased by 68% when the back support was installed [
8]. Therefore, the lateral displacement is suppressed from the combination of inclined front and back supports, which is a feature of an IER structure. Based on an analysis of the stability and field applicability, as well as field experiments and laboratory tests of the soil deformation in clay ground, it was confirmed that the spacing between the front and back supports is important [
9]. However, such studies did not analyze the dynamic stability required in the field where structures are maintained for a considerably long period, including embankment construction for a double line high-speed railway or rockslide prevention during the installation of a retaining wall (
Figure 2).
In this study, the displacement restraining effect of the IER method applied to embankment construction was analyzed by considering the ground and analysis conditions used by Seo et al. [
10] under dynamic loading. In addition, a two-dimensional (2D) dynamic numerical analysis was conducted. By analyzing the lateral and vertical displacements, the effect of the presence or absence of the back support was compared and the displacement restraining effect was investigated for the same H-pile per length.
3. Two-Dimension Dynamic Numerical Analysis
PLAXIS (ver. 2018, PLAXIS, Delft, Netherlands) [
11] was used for the 2D dynamic numerical analysis. For the analysis site, the earth retaining method was applied for a high-speed double-line railway line through an embankment. An embankment generally uses a retaining wall structure, for which many static and dynamic studies are currently being conducted. Although a retaining wall is structurally different from an IER structure, many studies have been carried out on a dynamic numerical analysis of an embankment ground. In particular, among the retaining wall structures, a reinforced earth retaining wall reduces the load acting on the front wall by increasing the strength of the embankment ground using a grid. To take advantage of the structural characteristics of the back support of an IER structure, studies applying a reinforced earth retaining wall have been further referenced. Ling et al. [
12] analyzed a reinforced earth retaining wall applied to an embankment ground and used special three-node elements for each block and the boundary surface of the retaining wall front plate. In general, a constraint method has been adopted to achieve a more stable numerical analysis of the interface than that of a Goodman element used in FLAC. Leshchinsky and Vulova [
13] applied an analysis method and included a stiffener between the blocks. In the analysis described in [
13], the stiffener was treated with a beam element.
Jiang et al. [
14] presented a 2D numerical analysis of a well-monitored pile–slab-supported embankment for the Beijing–Tianjin high-speed railway in China. The piles, which are arranged in a square shape, are a three-dimensional (3D) problem. To convert a 3D problem into a 2D problem, the equivalent elastic modulus of the pile wall was calculated as follows:
Eeq =
Epas +
Es (1 −
as), where
Ep and
Es are the moduli of the pile and soil, respectively; and
as is the area replacement ratio, defined as the ratio of the pile cross-sectional area to the total area. In this study, an analytical approach that considers the wall area was used, and a retaining wall without a reinforcement was applied as an elastic model. For the 2D analysis, the wall was modeled by calculating the area of the H-beam constituting the retaining wall.
The backfill can be simulated by assuming an elastoplastic material or an elastic material similar to that used in a general geotechnical structure problem. In general, a backfill is modeled using an elastoplastic model that includes a Mohr–Coulomb failure. Rowe and Ho [
15] assumed that the elastic modulus of the backfill is constant regardless of the height and set the dilation angle for various friction angles. This study analyzed the forces acting on the wall for various friction angles and found that the friction angle of the backfill material is an important parameter determining the magnitude of the wall acting force. Bathurst and Hatami [
16] estimated the elastic modulus by considering the increase of the stiffness of backfill along the depth using the Mohr–Coulomb failure model. However, they did not consider the change in the elastic modulus according to the change in stress during the analysis. Hatami and Bathurst [
17] modeled the backfill soil as a nonlinear elastic response prior to the peak strength, which was simulated using the stress-dependent hyperbolic model proposed by Duncan et al. [
18]. In this study, the Mohr–Coulomb model was applied, and results obtained from the field and laboratory tests were applied.
Simulations of the boundary surface vary from researcher to researcher. The method for simulating the interface can be classified as follows: 1) using thin continuum elements [
16], 2) specifying different nodes at the same boundary surface and adjusting the relative displacement between the nodes [
15], 3) designating specific adhesion stiffnesses and strengths for the shear force at the boundary surface to specific elements [
13], 4) applying a separate element to the boundary surface using a constraint method [
12], and 5) applying a connection element comprising a spring plastic section at the same node position between the different elements [
13]. Different researchers have utilized different methods and elements to simulate the boundary surface. However, the term “interface” is used in the same way. Such studies do not seem to satisfy the same interface condition, although the results of a comparative analysis of an interface simulation method have yet to be presented. In this study, the interface is modeled by applying an interface method in which the friction coefficients of the wall and backfill interact with each other.
The boundary condition of the ground is fixed in the horizontal direction along the x-axis and in both the horizontal and vertical directions along the y-axis. Load A has a uniform distribution, and the railway sleepers, rails, and lower ground are applied as fixed loads. Therefore, the applied surcharge loading was 50 kN/m
2, which is the load amount of a high-speed railway [
19]. Load B is a dynamic load acting on the railway and is applied to both sides of the rail in a concentrated manner. In a static analysis [
10], a surcharge load is applied with a uniform distribution. However, in a dynamic analysis, the dynamic sinusoidal loading is applied when considering that the high-speed railway is acted upon repetitively and at rapid speeds.
The railway target used for the analysis moved at a speed of 300 km/h and had a length of 380.14 m, and the railway loads were applied during 46 cycles to simulate the ground points of the railway wheels. The railway load was 170 kN, and a load of 85 kN was applied to each side [
20,
21,
22,
23,
24,
25]. In addition, the time taken for the railway to pass through a single point was 4.56 s. During this interval, a frequency of 10.1 Hz was applied for a load of 46 cycles. Finally, the analysis time for the application of the dynamic load was 4.56 s. To analyze the residual behavior after the railway passed by, an analysis time of 5 s was added without a dynamic load.
Figure 4 shows the analysis section, and
Figure 5 illustrates the sinusoidal input wave with a frequency of 10.1 Hz that was applied for the analysis.
3.1. Analysis Conditions
The analysis conditions are based on the IER structure used in the field test, and to conduct a comparative analysis according to the number of H-piles per 10 m of retaining wall, six analysis cases were applied, as shown in
Table 3.
Figure 6 shows the installation view and number of H-piles per 10 m. Analysis cases 1–4 represent the installation effects of the back support when the number of H-piles per 10 m was 10 or 5. Here, the installation condition of the back support differs. For cases 5 and 6, the effects of the H-pile installation arrangement on the front and back supports were analyzed when the same number of H-piles were utilized for the IER construction, and the installation interval between the front and back supports was reversed.
3.2. Design Parameters
3.2.1. Soil Parameters
The soil parameters used are the same as those applied in the analysis by Seo et al. [
10] and are shown in
Table 4. The material constitutive model used was the Mohr–Coulomb model. The soil parameters of soft rocks and weathered soil were used as the values of the ground survey report [
26]. In the case of the muck layer, a ground survey was conducted, and the test results were utilized. The parameters for all layers were applied to the results of the study by Seo et al. [
10]. These parameters are shown in
Table 4.
3.2.2. H-pile Parameters
The dimensions of the H-pile were 300 × 300 × 10 × 15 mm, which are the same as those used in the field test. During the analysis, the support was applied as a linear elastic model, and the axial stiffness (EA) and bending stiffness (EI) were calculated using Equations (1) and (2), where the elastic modulus of the H-pile (
Es) was 210,000 MPa. These equations were derived with reference to the Plaxis manual [
26].
Table 5 shows the H-pile parameters according to the installation interval (
c.t.c.).
where
Es is the elastic modulus of the H-pile (kPa),
A is the area of the H-pile (m
2),
IX is the moment of the second order of the H-pile (m
4), and
c.t.c. is the installation interval of the H-pile (m).
4. Results of Dynamic Behavior Analysis
A numerical analysis was conducted for types A–F, where type A is D-F2-B2-10, type B is S-F1-10, type C is D-F4-B4-5, type D is S-F2-5, type E is D-F4-B2-7.5, and type F is D- F2-B4-7.5. The results of the dynamic behavior analysis for each type are shown in
Figure 7,
Figure 8,
Figure 9 and
Figure 10.
Figure 7 and
Figure 8 show the lateral and vertical displacements at the top of the retaining wall, respectively, during dynamic loading.
Figure 9 and
Figure 10 show the lateral and vertical displacements at the bottom of the retaining wall, respectively, under the same conditions. In the case of the lateral displacement results, a (+) value indicates movement inside the retaining wall (right direction in
Figure 4), and the (-) value indicates movement outside the retaining wall (left direction in
Figure 4).
Table 6 shows the maximum displacement according to type. The dynamic load of a passing railway at a speed of 300 km/h was utilized, and a residual time of 5 s was used to examine the residual behavior after dynamic loading. As a result, vibrations occurred in the retaining wall during dynamic loading, which gradually decreased during the residual time.
Figure 7 shows the results of the lateral displacement at the top of the retaining wall. As shown, types B and D, which are single supports, have a much larger lateral displacement than a dual support. This result indicates that a dual support has a greater load resistance than a single support. This is because the effects of the front and back supports are combined to transfer the support effect, increasing the efficiency. Comparing the results for the single supports, it is evident that type D, with a smaller number of H-beams, has a larger displacement than type B. Comparing the results for the dual supports, the largest lateral displacement occurred in type C, and the smallest lateral displacement occurred in type A. The larger the number of H-beams, the greater the load resistance. Types E and F were compared with the same number of H-beams but for different numbers in the front and back. It was determined that the lateral displacement of type F was approximately 2 mm larger than that of type E, but no other significant differences were shown. The results indicate that if the support has the same number of H-beams, the distribution of beams, when relatively large in number at the front or back, did not have a significant effect.
Combined with the results for the lateral displacement, it can be determined that the most significant influence on the lateral displacement is the support type, namely, single or dual support, during dynamic loading. The next most influential factor is the number of H-beams, whereas the least influential factor is their arrangement, i.e., their distribution between the front and back for the same number of H-beams.
Figure 8 shows the vertical displacements at the top of the retaining wall. According to the results of the support type, a greater amount of vertical displacement occurred with the single support than with the dual support. This result indicates that a dual support has greater resistance to a vertical displacement under a dynamic load than a single support. To analyze the effect of the number of H-beams, the results for dual support types A and C were compared along with single support types B and D. For the single support, the vertical displacement was larger for type D with fewer H-beams than that of type B. For the dual support, type C, with fewer H-beams, had a larger vertical displacement than type A. When the difference in the number of actual H-beams doubled, the difference in the vertical displacement was 1.7 times for the single support and 1.5 times for the dual support. According to these results, the structure had a greater effect on the dual support. Moreover, even if the number of H-beams is reduced by a factor of 2, the decrease in resistance against the vertical displacement is smaller than that of a single support. To analyze the effect of placing many H-beams on the back support, types E and F were compared. It was determined that type E, with more H-beams on the back support, had smaller vertical displacements than type F, with more H-beams on the front support. As such, it can be concluded that the resistance to a vertical settlement improves when more H-beams are installed on the back support than on the front support.
Based on the results for vertical displacement, the most significant influence is the support type for a single support or dual support during dynamic loading. In a dual support, the next most influencing factor is the H-beam arrangement of the front and back supports, whereas the least influential factor is the number of H-beams.
Figure 9 shows the results of the lateral displacement at the bottom of the front support. The single support resulted in displacement in the ground, whereas the dual support showed the opposite behavior with displacement occurring outside the ground. However, the lateral displacement of the single support was approximately equal to that of the dual support. In this case, the effect of the difference in the number of H-beams was not significant, and the overall amounts of displacement were extremely small and were considered negligible. However, the direction of lateral displacement for a single support occurred outside the ground (left), whereas for a dual support it occurred inside the ground (right). This result shows that the support type affects the direction of displacement.
Figure 10 shows the results for a vertical displacement at the bottom of the front support. These results were very small and similar to those of the lateral displacement. However, the displacement of the single support was approximately 4 to 5 times larger than that of the dual support. As with the lateral displacement shown in
Figure 9, the vertical displacement did not significantly affect the number of H-beams or their arrangement on the front and back supports.
Based on these results, the type of support is the most important factor for the bottom layer of the retaining wall, whereas other factors including the number of H-beams and the arrangement of H-beam supports are considered to have little effect.
6. Conclusions
A dynamic numerical analysis was conducted to investigate the stability of an inclined earth retaining wall used during an embankment construction. The main conclusions are summarized as follows:
1. A dual support is more effective for suppressing lateral and vertical displacements than a single support. The resistance to the dynamic loading of a single support was shown to be smaller than that of a dual support.
2. The most significant factor influencing a lateral displacement at the top of a retaining wall is the support type. The next most influential factor is the number of H-beams, whereas the least influential factor for the same number of H-beams is their arrangement.
3. As with a lateral displacement, the most significant influencing factor on a vertical displacement at the top of the retaining wall is the support type. For a dual support, the next most influential factor was the H-beam arrangement, unlike a lateral displacement, and the least influential factor was the number of H-beams.
4. The type of support is an important factor affecting both the lateral and vertical displacements at the front of the retaining wall. Other factors, such as the number of H-beams and the arrangement of the H-beam of the support, did not have much of an effect. The support type affected the direction of the lateral displacement behavior at the bottom of the retaining wall.
5. The lateral displacement of the IER structure was shown to decrease by 81.05% at maximum compared to a one-row retaining wall with the same H-pile number per 10 m. An IER structure is therefore considered to be more stable.
6. The lateral displacement decreased by 89.24% when the structure was installed under the same conditions with or without a back support. It is thought that the back support applied to an IER structure demonstrates a suppression pile effect. The installation of the back support was shown to increase the structural stability.
7. For the same IER structure, the lateral displacement decreased by approximately 30% when the installation spacing of the back support was reduced by twofold.
8. By using the operating load of the railway, the dynamic numerical analysis was conducted. The results show the substantial displacement suppression effect by using the IER method. Using IER method may significantly affect the stability of the embankment construction so that the use of auxiliary method such as anchor installation may be eliminated. In addition, the construction period can be shortened and economic efficiency is expected to be increased. Therefore, it is recommended to use the IER method in the future construction of embankments.