Abstract

The hydraulic erosion of the subgrade slope can lead to serious loss of filling materials and landslides, affecting stability. In particular, there is no quantitative consideration of the scouring calculation method for the design codes of High-Speed Railway Subgrade Slope (HRSS). This paper mainly focuses on the anti-erosion stability of HRSS with coarse-grained soil and analyzes the current research shortcomings. The scouring on the subgrade slope was analyzed by establishing an incipient flow velocity formula that optimized the derivation process and parameters. Compared with the two types of cases, the exactitude of the formula derivation is verified. Finally, combined with orthogonal grouping, a single-factor analysis was conducted on the slope gradient and soil shear strength. With the same particle size, the incipient flow velocity decreases with the increase of the slope gradient (29.7° – 45°). When the particle size is upwards of 0.1 mm, the slope gradient has a great impact on the incipient flow velocity. And as the slope gradient decreases, the acceleration of incipient flow velocity with the same particle size increases. The incipient flow velocity increases with the increase of C (Cohesion) or φ (internal friction angle). The value of φ influences more deeply in the particle size range of 0.1 – 20 mm, that is, to advance the resisting erosion capacity, the shear strength of the slope soil is encouraged to improve when designing slope protection measures.

1. Introduction

High-Speed Railway Subgrade Slope (HRSS), shown in Figure 1, perennially exposed to the natural environment, is prone to instability under rainfall. Landslides caused by the rainfall infiltration or erosion failure are the main issues that lead to slope instability and other foundation risks [15]. Some measures are adopted to reduce the possibility of landslides and erosion problems, including ecological protection and masonry protection. The slope instability caused by rainfall infiltration has been studied comprehensively [6]. HRSS is divided into embankment slopes and cutting slopes. The height and slope gradient directly determine the rainfall area, affecting the incipient flow velocity induced by erosion on the slope surface. Referring to the “Code for Design of High-Speed Railways” (TB10621-2014) [7], the height of embankment slopes for high-speed railways generally does not exceed 20 m, the upper height does not exceed 12 m, and the slope ratio is 1 : 1.5. The lower part shall not exceed 8 m in height, with a slope rate of 1 : 1.75, and the slope form shall be a broken line or stepped type. Depending on the engineering performance and grading characteristics, subgrade slope fillers can be divided into groups A, B, C, and D. Fillers of HRSS mainly use groups A and B fillers, coarse-grained soil. Nevertheless, there is room for improvement in the calculating system of erosion failure for subgrade slopes in the design codes. Rainfall erosion often happens on the slope surface, such as sheet erosion and rills. As the continuous and multiple scouring, erosion will further develop into the interior of the subgrade, resulting in landslides [810]. Therefore, the erosion failure of the subgrade slope which may lead to serious security risks should be quantitatively studied.

Closely related to the disciplines such as soil science, soil mechanics, hydraulics, and sediment dynamics, erosion issues are affected by various factors [11, 12]. Soil particles may be scoured out of the slope surface under the direct impact of rainwater and the drag force of the water flow, which should be quantitatively analyzed based on the particle size and the flow velocity. According to the German research results, Figure 2 noted the critical states of erosion for soil particles, which depend on the functional relationship between the soil particle size and flow velocity [13, 14]. Since the second half of the 19th century, people have deeply realized the necessity of quantitative research on soil erosion, and are eager to try to establish a simple and efficient mathematical method for quantitative prediction [1517]. Based on the observations in 10,000 areas over 30 years, the Universal Soil Loss Equation (USEL) [18, 19] was extensively used in predicting the average annual soil erosion in the United States. The design codes in Japan and China only mention some subgrade slope protection measures but do not mention the anti-erosion stability analysis method of subgrade slope caused by rainfall. The slope protection engineering of HRSS includes ecological protection and engineering protection. The engineering protection measures mainly include physical slope protection (wall), skeleton slope protection, anchor rod (cable) frame beam slope protection, and sprayed concrete slope protection. The ecological protection measures mainly include various forms such as spraying grass, planting belt (bag) grass, and guest soil planting.

The recent research on subgrade slope erosion mainly focuses on the mechanism of erosion starting, the erosion process and calculating the scouring amount. Based on the conversion relationship between the flow velocity field and shear field and the feature of erosion for granite eluvium cutting slope, a calculating formula model, including the drag force and the flow velocity, was established by Luo et al. [20]. However, the cohesion force, one of the parameters for soil, was not considered. Depending on the soil characteristics, the homogeneity and cohesion of soil on slopes were classified and discussed by Wang [21] to analyze the critical conditions for subgrade slopes notionally. Based on the law of energy conservation, a slope erosion calculation model was put forward by Li et al. [22], who compared the quantity of soil loss calculated by his model with the data of the erosion amount obtained by the on-site measure, which verified that the model had a good prediction effect. Figuring out the influencing factors of rainfall erosion and the hydraulic erosion process, an erosion calculation model and the evaluation method for railway slope proposed by Kuang [23] to optimize slope protection measures. To promote the calculating model, Guo [24] analyzed the mechanism of slope erosion initiation of cohesionless filler and cohesive filler, taking into account the structure of railway subgrade, filler, rainfall, and catchment characteristics, and obtained the variation rules of slope velocity and threshold particle size with rainfall intensity. Model tests are commonly used to study slope instability, due to the frequent use of artificial rainfall erosion in the experiments, which is often limited to the implementation of local rainfall patterns, there is a boundary constraint effect between the rainfall area and the non-rainfall area, manifested as the lateral infiltration of rainwater and the lag of slope deformation in the non-rainfall area on the rainfall area. A deep-buried isolation trench test method has been proposed to address these boundary effects [25]. The problem with the existing rainfall warning system for the prevention of railway disasters is that the determination of rainfall warning values is mainly by statistical data or experience, so lacking deep analysis of the mechanism of slope topsoil disaster induced by rainfall [26]. This situation causes the rainfall warning values to not match actual disaster grades. The main railway subgrade slope disaster is topsoil slip during rainy days.

To sum up, the existing scouring starter theory has not comprehensively pondered the soil properties, and the scouring starter formula does not fully reflect the soil-related parameters. Furthermore, existing formulas for the initiation of slope erosion are often derived indirectly through certain assumptions using the conversion from the horizontal surface to the slope surface, resulting in miscalculation. The coarse-grained soil, with firm water permeability and high-compressive strength, is normally applied to high-speed railway subgrade engineering. This paper quantitatively deduced and verified a formula for critical incipient flow velocity of scouring of HRSS with coarse-grained soil induced by rainfall erosion. On this basis, the erosion mechanism for soil particles has been clarified.

2. Basic Assumptions

To optimize the derivation process, the basic assumptions are as follows:(1)The soil particles are spherical;(2)When the water flow satisfies the Bernoulli principle, and the soil particles leave the slope height relative to three times the particle size, the uplifting force basically disappears;(3)The intergranular action, charge gravity, bite force, cementation force, and cohesion of soil particles and shear strength parameters are uniformly simplified to adhesion;(4)The surface infiltration conditions of the slope surface are good, and the infiltration pressure is consistent with the infiltration direction;(5)The direction of rainfall is the same as the direction of gravity;(6)When the critical start of scouring, the friction resistance flow rate is a fixed value;(7)When the critical start of scouring, the friction coefficient is a constant fixed value;(8)The critical starting of scouring is static;(9)The slope erosion is a strong turbulent effect, and the viscosity of the water flow is not considered;(10)The length of turbulence mixing can be determined using Carmen’s turbulence similarity hypothesis.

3. Derivation of the Critical State of Velocity Scouring Formula on the Slope Surface

3.1. Force Analysis for Soil Particles on the Slope

Under the action of water flow, soil particles on the slope surface of the subgrade are subjected to the floating gravity w′ (N), the drag force FD (N), the lift force FL (N), and the infiltration pressure FS (N) and cohesion C (N).where γs is the unit weight of soil particles (N/m3); γw is the unit weight of water (N/m3); D is the grain size (m).

On the basis of sediment motion mechanics, drag forces and lift forces are exerted by liquid phase water flow on solid particles. During the water scouring, friction will occur when the water scours the surface of soil particles. When the Reynolds number of the water flows is slightly higher, the streamlines at the top of the soil particles will separate, and vortices will exist on the surface of the soil particles, resulting in a pressure differential between the front and rear of the soil particles, forming the form resistance. The drag force FL is the combination of the friction induced by water flow on the surface of soil particles and the shape resistance generated by vortices. In light of Bernoulli’s principle, when water flows on the slope, the velocity at the top of the soil particles is significantly greater than the velocity of the infiltration flow between the soil particles at the bottom, resulting in an uplift force FL.where CD is the drag force coefficient; CL is the lift force coefficient; ρ is the density of water (kg/m3); u0 is the velocity of flow on the soil surface (m/s).

After the soil particles leave the slope and move in a push mode, on the one hand, they begin to be subjected to upper water with faster flow velocity, resulting in increasing the drag force, and on the other hand, the lift force tends to disappear due to a decrease in the pressure differential between the upper and bottom sides of the particles.

When rainfall accumulates on the slope, the infiltration pressure FS generated by the infiltration of the slope surface is consistent with the infiltration direction.where Js is the hydraulic gradient of infiltration flow in soil. For unsaturated soil slopes, when the slope gradient is relatively gentle, the seepage pressure enhances the stability of soil particles. Nevertheless, when the slope gradient is steep, the component of seepage pressure along the flow direction will intensify slope erosion.

According to the soil mechanics theory, the positive charge of the water molecule, attracted by the surface of soil particles with negative charges, is aligned to form binding water, generating cohesion C. Besides, C may exist from the cementation of soil particles by substances such as silicon, iron, and carbonate contained in the soil.

3.2. The Critical State for the Starting

Figure 3 demonstrates the force analysis for the soil particles on the slope. The combined force of infiltration pressure FS and the floating gravity w′ is, among the slope direction and perpendicular to the slope direction, (FS + w′) sin θ and (FS + w′) cosθ separately. Compared with the force mode under horizontal plane conditions, the floating weight, with the component force in the slope direction, exacerbates the erosion of the slope. Besides, contrary to the direction of FD, the angle between C and the slope axis is β.

The FD for soil particles in the direction β = 90° can be determined as follows:where φ is the interior friction angle.

Substituting Equations (1)–(4) into Equation (5) yields:where u0 = UR·fi(UR·D/v); UR is the frictional velocity; and UR = (τc/ρ)0.5; fi(UR·D/v) is the function of Reynolds number. The critical state of drag force can be obtained as follows:where

Substituting Equation (8) into Equation (7) yields:

From Equation (9), the critical state of drag force for soil particles is primarily determined by the soil particle size, the slope gradient, the internal friction angle, and cohesion.

For the laminar flow, the relationship between the flow velocity and the flow shear force can be simply determined as follows:where μ is the coefficient of viscosity for water; u is the flow velocity (m/s); τ is the flow shear force (Pa); y is the distance from the flow section to the slope surface (m); and v is the dynamic viscosity coefficient.

For the turbulent flow, the shear force of the flow is much more complex. Based on the mixing length theory for turbulent flow, the fluctuating velocity among the flowing direction and depth direction can be obtained, respectively:where l is the mixing length of turbulent flow (m).

Based on the fluid mechanics, the turbulent shear formula can be determined as follows:

Substituting Equations (11) and (12) into Equation (13) yields:

When the turbulence of water flow is weak, the water flow shear force should be the sum of the viscous shear force and the turbulent shear force as follows:

However, it is unnecessary to consider the viscous force of water flow in the analysis of slope erosion owing to its strong turbulence.

The turbulent mixing length can be determined by Kármán’s turbulent similarity hypothesis. Assuming that the flow velocities at the distances of y1 and y2 from the slope surface are u1 and u2, respectively, the variation of which can be determined as follows:

According to Kármán’s turbulent similarity hypothesis, there is the following relationship:and

The turbulent mixing length can be obtained as follows:where k is the proportional coefficient, that is, the Kármán constant.

The velocity field can be obtained as follows:

The linear distribution of shear force along a vertical line in a two-dimensional flow can be determined as follows:where τy is the water flow shear force at the distance of y from the slope surface (Pa); τ0 is the water flow shear force at the slope surface (Pa).

Substituting Equation (20) into Equation (21) yields:where uh is the flow velocity at the water flow surface (m/s); uy is the flow velocity at the distance of y from the slope surface (m/s).

When UR remains constant, (uh − uy) can be considered only varied with y. Besides, near the slope surface where erosion occurs, the flow shear force, varied insignificantly, can be approximated as equal to the slope shear force τ0.

The vertical velocity distribution formula for slope flow can be obtained as follows:where y0 is the distance from the slope surface where the velocity of water flow is equal to zero (m).

Integrating Equation (24), the average velocity of the water flow section, divided into two situations, can be obtained:(1)Smooth slope surface(2)Rough slope surfacewhere U is the average velocity of the water flow section (m/s); R is the hydraulic radius (m); ks is the surface roughness (m).

Considering the roughness of the slope surface, Equation (26) is transformed as follows:where χ is the correction coefficient.

Substituting Equation (9) into Equation (27) yields:

4. Case Study and Equation Verification

To verify the accuracy of Equation (28), a comparison is conducted with some case studies, including horizontal channel cases and a subgrade case.

4.1. Verification by Field Tests and Equations of the Horizontal Slope (θ = 0°)

The same parameters and filed data, conducted by Dou [27], are used for the comparison. Moreover, some empirical formulas for incipient flow velocity on the horizontal slope were proposed by Dou [27] and Ma [28], respectively. The unit weight of soil particles γs is 26.5 kN/m3. The interior friction angles φ and cohesion C equal 38° and 0 kN, respectively. The correction coefficient χ is taken for 1, and the hydraulic radius R is 0.15 m. The surface roughness ks and the function of Reynolds number f(UR·D/v) equal 0.005 and 0.04, respectively. For the coarse-grained soil, size between 0.075 and 20 mm, the curve of the present study is located at the lower edge of the test results or other empirical formulas, shown in Figure 4, which is in good agreement with the test results.

4.2. Verification by the Equation of the Subgrade Slope (θ = 33.7°)

The formula in this paper is directly derived from the soil particles on the subgrade slope, instead of assuming the transformation relationship from the horizontal surface to the slope surface. The same parameters in the calculation formula for the subgrade slope, conducted by Guo [24], are used for comparison, shown in Figure 5. The unit weight of soil particles γs is 26.5 kN/m3. The interior friction angles φ and cohesion C equal 38° and 20 kN, respectively. The correction coefficient χ is taken for 1, and the hydraulic radius R is 0.015 m. The surface roughness ks and the function of Reynolds number f(UR·D/v) equal 0.005 and 0.04, respectively. For the whole size between 0.075 and 20 mm, the curve of the present study is in good agreement with Guo’s empirical equation. In the vast majority of the curve, the maximum difference between Guo’s and the present study is within 5%.

To sum up, the incipient flow velocity equation for soil particles in this paper is not only suited for the horizontal slope but also has a high accuracy for the subgrade slope.

5. Analysis of Impact Factors of the Incipient Flow Velocity

5.1. The Incipient Flow Velocity under Different Reynolds Number

The incipient flow velocity varies with different Reynolds numbers, shown in Figure 6. The unit weight of soil particles γs equals 26.5 kN/m3. The interior friction angles φ and cohesion C equal 38° and 20 kN, respectively. The correction coefficient χ is taken for 1, and the hydraulic radius R is 0.015 m. The surface roughness ks and the slope gradient θ equal 0.005° and 33.7°, respectively.

For coarse-grained soil, the function of Reynolds numbers is usually taken as 0.04–0.045. Thus, in this range, with the same particle size, the incipient flow velocity increases with the increase of the f(UR·D/v) but impacts less. Thus, the f(UR·D/v) is taken as 0.04 in the following analysis.

5.2. The Incipient Flow Velocity under Different Slope Gradients

To analyze the affection of slope gradients to the incipient flow velocity, four common subgrade slope gradients (tan θ) are shown in Figure 7. The incipient flow velocity under different slope gradients is shown in Figure 7. The unit weight of soil particles γs equals 26.5 kN/m3. The interior friction angles φ and cohesion C equal 38° and 20 kN, respectively. The correction coefficient χ is taken for 1, and the hydraulic radius R is 0.015 m. The surface roughness ks and the function of Reynolds number f(UR·D/v) equal to 0.005 and 0.04, respectively.

With the same particle size, the incipient flow velocity decreases with the increase of the slope gradient. When the particle size is less than 0.1 mm, the incipient flow velocity is slightly affected by the slope gradient, conversely, the slope gradient has a great impact on it. And as the slope gradient decreases, the acceleration of the incipient flow velocity with the same particle size increases.

5.3. The Incipient Flow Velocity under Different C or φ

The orthogonal grouping of C and φ was conducted to analyze the affection of soil shear strength on the incipient flow velocity on the slope, shown in Table 1 and Figure 8. The unit weight of soil particles γs and the slope gradient θ equal 26.5 kN/m3 and 33.7°, respectively. The correction coefficient χ is taken for 1, and the hydraulic radius R is 0.015 m. The surface roughness ks and the function of Reynolds number f(UR·D/v) equal 0.005 and 0.04, respectively.

The incipient flow velocity increases with the increase of C or φ. The value of C impacts greater in the particle size range of 0.075–0.1 mm, while the value of φ influences more deeply in the particle size range of 0.1–20 mm. With the φ decreasing by 20.8%, the incipient flow velocity of Group 3 is much larger than that of Group 5 or Group 8, although the C increases to 200%.

5.4. The Incipient Flow Velocity under Different Slope Height

Considering the impact of slope height, the issue has been simplified by comparing the incipient flow velocity for the same soil size on different parts of the slope surface. In this paper, the derivation of the incipient flow velocity is based on soil particles at any position on the slope as the research object. Therefore, the incipient flow velocity for soil particles with the same particle size will not change with the change in slope height. Nevertheless, under the same rainfall intensity, affected by gravity, rainwater flow speed accelerates along the slope. The higher the slope height, the faster the flow velocity at the toe of the slope is, resulting in more severe erosion. Thus, in further research, establishing the relationship between the incipient flow velocity and rainfall intensity, the influence of slope height can be quantitatively calculated.

6. Discussion

In this paper, the incipient flow velocity of HRSS with coarse-grained soil is studied induced by rainfall erosion, based on theoretical derivation and numerical analysis. In terms of applicability, the formula applies to situations where cohesive soil is present on slope surfaces or horizontal river channels. The masonry skeleton, often used on the subgrade slope, is a common method to reduce the rainfall erosion. Thus, with the incipient flow velocity equation established in this paper, we can quantitatively calculate the erosion condition. The spacing of masonry skeleton can be optimized by combining past measured results to reduce costs. The value of C and φ mainly comes from practical engineering experience. An indoor rainfall erosion model test, which operates the essential process, will be conducted and summarized to optimize parameters and validate the formula. In addition, under the same slope height, a change in slope gradient will inevitably change the rainfall area on the slope surface. Although the incipient flow velocity of filler particles decreases with the increase of slope gradient, at the same time, the rainfall area reduces on the slope, which leads to the decrease of the slope runoff. Thus, the impact of slope gradient on the stability of slope erosion resistance needs further research. The flow velocity and raindrop impact are closely related to the rainfall intensity. The raindrop impact and the sheet flow resistance caused by raindrop impact also influence the starting of soil particles [29, 30]. By establishing the quantitative relationship between the flow velocity, raindrop impact and rainfall intensity, the impact of rainfall intensity on the incipient particle size can be figured out. The surface soil on the slope is composed of soil particles with different particle sizes. When analyzing the slope instability under erosion, it is necessary to consider determining the grading criteria or combining random theory for research. Many factors affect the stability of HRSS, such as rainfall characteristics, slope morphology characteristics, slope rock and soil properties, and slope protection measures. The stability evaluation of slope erosion is complex, involving significant uncertainty. Recently, many research methods concentrate on using fuzzy mathematics and neural networks to conduct stability analysis. Subsequently, after finishing model tests, our team will establish the database for HRSS instability induced by erosion with the evaluation indicators, including quantitative, semiquantitative, and qualitative evaluation methods.

7. Conclusion

There are few studies on the incipient flow velocity of HRSS. An optimized formula, considering the cohesion of coarse-grained soil, is established by deriving directly on the subgrade slope. The incipient flow velocity equation for soil particles in this paper is not only suited for the horizontal slope but also has a high accuracy for the subgrade slope. The main conclusions are as follows:(1)This article studies the erosion initiation law of HRSS and establishes the incipient flow velocity equation suitable for multiple working conditions, including coarse-grained soil subgrade slopes with different slope gradients.(2)Different flow velocities induced by rainfall result in different spacing of masonry skeleton on the subgrade slope. The equation incipient flow velocity can be adopted to optimize the spacing of the masonry skeleton to reduce costs.(3)Compared with the existing empirical formulas for subgrade slope, the accuracy in the present study has been verified.(4)With the same particle size, the incipient flow velocity decreases with the increase of the slope gradient. When the particle size is less than 0.1 mm, the incipient flow velocity is slightly affected by the slope gradient, conversely, the slope gradient has a great impact on it. And as the slope gradient decreases, the acceleration of incipient flow velocity with the same particle size increases.(5)The incipient flow velocity increases with the increase of C or φ. The value of C impacts greater in the particle size range of 0.075–0.1 mm, while the value of φ influences more deeply in the particle size range of 0.1–20 mm. With the φ decreasing by 20.8%, the incipient flow velocity of Group 3 (C = 10 kN, φ = 48) is much larger than that of Group 5 (C = 20 kN, φ = 38°) or group 8 (C = 40 kN, φ = 38°), although the C increases to 400%.(6)The research results on the starting flow rate in this article can provide a theoretical basis for the protection of HRSS. The shear strength of the slope soil is encouraged to improve when designing slope protection measures.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This work was supported by the Scientific Research Project of China Academy of Railway Sciences Co., LTD (2021YJ070).