Abstract

Installing geosynthetic reinforcement at the bottom of the embankment will provide positive support for subgrade stability and settlement control, if there is a void or a weak foundation. In routine design work, the geomembrane effect must be well estimated and the tensile strain should be precisely predicted. Conventional analytical methods often adopt the limit state method to calculate the overlying load on the deflected geosynthetic. However, this assumption does not necessarily apply to all conditions, especially when the foundation soil can provide certain resistance. In this study, a semiempirical prediction method for evaluating the geomembrane effect of the basal reinforcement was proposed, and an iterative solution for calculating tensile strains of a deflected geosynthetic was deduced. In derivation, a virtual inclined slip surface and interaction between the geosynthetic and soil were quantitatively evaluated by coupling the arching effect and the geomembrane effect. Moreover, the development of shear stress along the slip surface can be considered, as well as different segments of the basal reinforcement. Then, the proposed method was validated by two large-scale experiments. Comparison of the results of this method with measurements and results of other analytical models confirmed that this analytical method can take good care of the varying process of the localized sinking, regarding the overlying loads on the geosynthetic and the subsequent tensile strains.

1. Introduction

Constructing road embankments on soft soil foundation, or there exists sinkholes or voids underneath, faces the threat of excessive settlement or different settlement. Although lightweight fill, piles, and subgrade enhancing technologies are always adopted to solve the problem, as well as other engineering issues, implementing a geosynthetic is often necessary to reduce settlement and enhance stability [17]. In reinforced structures, the geosynthetic mainly plays a load redistribution role between the subsided zone and the surrounding less-deformed zone. Moreover, the geomembrane effect always occurs along with the soil arching effect. Subjected to local subsidence at the fill base, the overlying fill deflects the basal reinforcement and is also prevented from subsiding by the geosynthetic simultaneously. Aiming to desire a reasonable design of the reinforced structure, an important task is to accurately evaluate the geomembrane effect, i.e., calculate the tensile strains. However, it is not practical to consider the tension of a geosynthetic solely without investigating the overlying loads. Therefore, the arching effect and the geomembrane effect, both of which are key design points, should be quantitatively evaluated. The former is to obtain the load acted on the geosynthetic, while the latter is to determine the maximum tensile strain of the installed geosynthetic.

In general, the portrayals used for depicting soil arches can be divided into two categories. One is to derive the static stress balance of a horizontal thin strip analysis unit in the vertical direction, represented by the Terzaghi [8] method. The other is to deduce the stress expression of the loosening zone of Earth pressure influenced by arching based on different nonvertical curvilinear slip surfaces [9, 10]. Arching behavior is always influenced by the geomembrane effect and the resistance from the compressible subsoil [1114]. For the family of frictional arching models, the shear stress along the slip surface is mainly dependent on the differential settlement at the fill base. Several methods have been found in the literature for evaluating the mechanism of shear friction that is significantly affected by arch behaviors or geosynthetic layers [1518]. Those methods are usually calculated in two steps: the first is to deduce the pressure, and the second is to compute the tensile strain of the geosynthetic subjected to the previous computational pressure.

Based on the fixed geometry assumption (such as circular, parabolic, or their combination), traditional geomembrane evaluation methods usually consider that the geosynthetic layer is fixed at one edge and the tensile strain is uniformly distributed within the full-length range [16, 19, 20]. Moreover, the tensile strain is only dependent on the ratio of the reinforcement’s deflection to the cavity diameter or subsiding width, instead of the overburden load and other material properties. However, nonuniform distribution and peak strain existence within the geosynthetic reinforcement have been gradually confirmed in real scenarios and numerical simulations [2124], indicating that the working condition of the horizontal installed reinforcement should be characterized more realistically. Then, improved evaluation methods are strongly recommended, and more factors (such as the interfacial friction between the fill material and reinforcement) have to be considered in the geomembrane effect [2529].

In a low embankment, as the two-dimensional stress state is considered, the settlement area at the top of the embankment is larger than that of the base. This indicates that the slip surface of the subsiding embankment fill against the motionless part resting on both sides undergoes a process of changing from the vertical direction to inclined directions. Therefore, the reinforcing behavior of the geosynthetic layer is quite different from the conventional ideas, accounting for the soil arching effect by utilizing a vertical [8], a spiral [30], a semicircular [31, 32], or a triangular [33] physical model presented by the aforementioned literature. Using a similar configuration, the aim of this study is to present a calculation method for the geomembrane effect with regard to nonuniform tensile strain of the geosynthetic reinforcement. In derivation, an improved arching mechanism based on mobilized shear stress within the inclined slip surface is evaluated, as well as the identified reinforcing mechanism of the horizontal geosynthetic layer. Finally, the method is validated by two physical model tests.

2. Developing Arching Effect on Nonvertical Slip Surface

Figure 1 illustrates a profile of a 2D analytical model of the basal reinforced embankment subjected to localized sinking. The subsidence width at the base and top are b and Bs, respectively. The height of the fill is H. Terzaghi [8] has pointed out that the angle of the inclined slip surface would vary from 90° to 45 + φ/2 with the increasing ratio of H to b of the fill mass in the famous trapdoor test. Therefore, a vertical slip surface assumption would not be fully conformed to the actual situation in a low fill mass, and the subsidence basin would extend to a degree due to material cohesion or internal friction angle of the fill. Then, boundary trajectories of the slip surfaces can be simplified by two curves, a′b and cd′, with an inclined angle to the vertical direction.

In Figure 2, a soil strip element ABC at the bottom, representing the mechanical behavior of the fill within the stress release region due to arching, is selected to conduct the equilibrium analysis of vertical stresses. The soil strip element can be idealized as a simple beam that is supported at two pivot points. The deflection of the soil element is assumed to be δ, and the deflection angle at point A is θ. The positive stress perpendicular to the slip surface is σns, and the shear stress is τs. Different from the definition of lateral Earth pressure coefficient K proposed by Krynine [34], the coefficient here changes with the position along the slip surface and is also related to the internal friction angle φ of the fill material.

Figure 3 shows a mobilized stress Mohr circle, in which the shear stress along the slip surface varies, to represent the stress state of a position starting from the edge of the basal subsiding area to the fill top surface. In traditional arching theory, the soil stress state along the slip surface is always considered to have reached the limit state. Therefore, the results calculated by the above limit state methods are always conservative, and their assumptions and factors are prone to fly in the face of the real scenario. It is believed that the magnitude of the shear stress along the slip surface is mainly dependent on the relative displacement at the fill base. Typically, when the embankment fill height is not high enough, or if the foundation soil underneath the geosynthetic provides a certain resistance, the Mohr circle will not be tangent to the shear strength envelope. In other words, an incomplete soil arching effect occurs, indicating that the maximum shear stress is lower than the ultimate shear strength.

Based on the above analytical assumptions, a unit length perpendicular to the paper direction is selected for a simplification of the two-dimensional stress equilibrium condition. At a depth z to the top surface within the subsided area, the width of the soil strip unit of EFG is termed as Bz. In the subsequent derivation, the cohesion of the fill material is not considered, and then through the stress equilibrium condition on vertical direction, the following equation can be obtained:where , is the gravity of the fill material, and is the vertical stress.

As shown in Figure 3, the shear stress can be expressed as follows:

According to the geometric position relation in the stress Mohr circle, the following formulations can be obtained:where and are the horizontal/vertical stresses of the soil arch, respectively; and are the major and minor principle stresses, respectively; and is the angle of related to the horizontal direction.

The lateral active Earth pressure coefficient is defined as follows:

When equation (3) is divided by , it becomes

If the lateral Earth pressure coefficient is defined as the ratio of lateral stress to vertical stress, one can get the following:

Here, the shear stress can be formulated as follows:

While is divided by and then substituted into equation (7), the following relationship can be established:

So, the shear stress of equation (7) can be reformulated as follows:

Assume that

The angle is determined by the geometric relation illustrated in Figure 2; when R is defined as the deflection radius, it can be expressed as follows:

Here,

Likewise, the normal stress acting on the idealized inclined slip surface a′b at a depth z to the top surface can be expressed as follows:where is the angle between the inclined slip surface and the vertical direction.

When equation (13) is divided by and the resulting equation is substituted in equation (5), the following relation can be established:

Let

Then, the shear stress can be reorganized to be as follows:

Let

Substituting equations (14) and (15) into equation (1) and dividing both sides by dz, one can get the following:

The solution of the vertical stress can be expressed as follows:

Letwhere is an unknown parameter. As illustrated in Figure 1, when z = H and , then one can get the following:

3. Equivalent Geomembrane Effect

The tension of the geosynthetic is dependent on the deflection, and the maximum tensile strain is recognized to occur at the edges of the subsiding basin [21, 22]. It is acceptable that the geomembrane effect would be more predominant due to the increasing differential settlement. In Figure 4, the reinforcement layer consists of three segments, the anchorage section on the supporting soil, the transition section due to tensile elongation, and the deflection section within the subsided area. To simplify the analysis, the shape of the deflected geosynthetic layer is idealized to be a fixed form, accompanied with the assumption that the vertical stresses on soil and reinforcement are uniformly distributed (see Figure 5). Theoretically, the load acting on the geosynthetic layer (within the deflected span) equals to the arch zone weight deducting the friction from the adjacent stationary fill mass and also equals to the conjunction of the geomembrane lift capacity and resistance from the compressible subsoil.

Refining from the analytical models proposed by Villard and Briançon [22] and Feng et al. [29], a similar evaluation method for the geomembrane effect of reinforcement is also presented in this study. In Figure 4, the geosynthetic layer deflects due to a uniform load qg acting on its upper surface, but generates a maximum tension Tc at the edge of the subsided area. The tangent angle at point C (in Figure 4) of the reinforcement is . At a random point M within the deflection, its tensile force is TM, with an oblique angle of . Adopting the Coulomb friction law (Figure 6) to define the interface friction between the soil and geosynthetic, one obtains the following expression:where and are the relative displacement and critical relative displacement between the soil and geosynthetic, respectively; is the friction angle between the soil and geosynthetic; and is the normal stress.

Checking a microelement of (see Figure 7), the equilibrium relationships can be established aswhere , , and are the tensile, vertical, and horizontal tensile component forces, respectively, of the geosynthetic at point M; is the inclination angle to the horizontal direction; and and are the deformed deflection of the geosynthetic and its derivative form at point M with an abscissa position x to the midpoint.

In order to make a further simplification, an assumed relationship can be established as follows:where and are the upper and lower friction strength of the interfaces between the soil and geosynthetic, respectively; is the reduction coefficient of the lower interface to the upper one.

Assuming that , where is the equivalent geomembrane stress upward due to deflection, and are the upper and lower vertical stresses, respectively, acting on the geosynthetic element. The following relationship can be established:

Then, one can obtain

and

Deriving the third term in equation (23), it can be reformulated as

Substituting equations (26) and (28) in equation (27) yields

Originally, there is also a relationshipwhere is the tensile strain of the selected analysis element and is the tensile stiffness of the geosynthetic. Then, the geomembrane effect of equation (29) can be reformulated as

Also, the tensile strain at point M can be expressed aswhere is the derivative form of the geosynthetic displacement at point M. Substituting equation (32) in equation (31) yields

The overlying load of the geosynthetic layer under soil arching is averaged in the range of the subsided area and yieldswhere ks is the reaction coefficient of foundation soil.

In theory, equation (19) is equal to equation (34). In this study, the parabolic shape of the geosynthetic is adopted to evaluate the geomembrane effect due to localized sinking. The only uncertainty is how to determine the deflection of the geosynthetic with a particular tensile stiffness.

As the transition zone between the deflection section and the anchorage section, there is an angle change of tensile forces, which is accompanied by the elongation of the installed geosynthetic. Combining the above effects yieldswhere and are the tensile force at point C of the deflected section and the tensile force of anchorage section, respectively; is the inclination angle of the tangent direction of the deflected geosynthetic; and is considered to be a comprehensive friction angle of the geosynthetic and soil. Villard and Briançon [22] suggested that this angle equals to the lower friction angle of the geosynthetic and soil. Occasionally, there exists a relationship that , which connects the deflection section and the anchorage section.

As illustrated in Figure 8, the anchoring behavior of the geosynthetic should be considered in two cases. One is that the relative displacement between the geosynthetic and soil does not exceed the critical displacement as defined in Figure 6, and the other is that the relative displacement between the geosynthetic and soil within the anchorage section has reached the critical state.

Selecting an element dr at a distance r to the edge point C of the subsided area, the force balance relation in the horizontal direction can be formulated as follows:where and are the upper and lower friction stresses, respectively, of the interfaces between the geosynthetic and soil within the anchorage section. To simplify the derivation, cohesion of interface is not considered, and the Coulomb friction law is adopted. Under a specific relative displacement U between the geosynthetic and soil, the friction can be formulated as follows:where is the normal stress acting on the anchorage length of the geosynthetic; and are the upper and lower friction angles, respectively, of interfaces within the anchorage section.

Then, equation (36) can be reformulated as

Considering the relation , when the relative displacement between the geosynthetic and soil does not exceed the critical condition, as illustrated in Figure 8(a), equations (36) to (38) can be deduced to bewhere . The general solution of equation (39) iswhere and are the constant coefficients depending on boundary conditions. When r = L (the anchorage length) and U = 0; when r = 0 and , and one can get

Now, the displacement of point C can be solved to be

When the relative displacement between the geosynthetic and soil within the anchorage section partially reaches the critical state, the anchorage section AC should be divided into two parts, the AD segment and the DC segment. The AD segment shares the same deducing procedure as derived above. For the DC segment, one can getwhere is the tensile strain of the geosynthetic at distance r to point C; and .

Then, the solution can be deduced to be

For the dividing point D, when , there exists relations, and . Therefore, the dividing point position can be determined bywhere .

Then, the relationship between the displacement and tensile forces can be established as

According to a fixed analytical model, such as definite mechanical properties of reinforcement materials and soils, as well as given geometries, the only unknown is the deflection of the geosynthetic. Obviously, the iterative method should be utilized to solve the above question, and some simplified conditions need to be assumed. For example, the stress diffusion in the horizontal direction is not considered, and load transfer adjustment only occurs in the vertical direction. Moreover, the separation between the geosynthetic and soil is not allowed, and the stress and displacement along the geosynthetic layer are continuously coordinated. Integrating the entire derivation process, the iterative calculation procedure can be illustrated as in Figure 9. To begin with an arbitrary original setting of the (e.g., a linear expression), the calculation starts with a fixed geosynthetic trajectory. Once the deflection of the geosynthetic is determined, the tensile strain can be deduced, as well as the equivalent applying stress. Meanwhile, the overlying load that is averaged within the subsided area can be calculated, and a comparison is made to the previous geomembrane effect. If the deviation is within the acceptable range (such as one thousandth), the iterative solution ends. Then, the anchorage section can be evaluated in two cases.

4. Model Validation and Discussion

This part compares the results of the present method and a full-scale experiment test, as well as the results of the analytical methods proposed by Villard and Briançon [22] and Feng et al. [29]. The test was carried out by a sudden deflation of an air balloon to investigate the reinforcing behavior of the basal geosynthetic (see Figure 10). Detailed information about the physical model is as follows: the cohesionless fill height is 0.5 m with gravity 17 kN/m3, the width and height of the balloon are 2.0 m and 0.25 m, respectively. A single layer of the geosynthetic is installed on the balloon, with more than 2.0 m of the anchorage length. The ultimate tensile force is 125 kN/m with regard to the tensile strain of 12%. The model achieves the effect of basement hollowing by means of airbag deflation. Therefore, there is no support from the underneath foundation soil, indicating that the geosynthetic layer fully accommodates the upper embankment load. The relevant calculation parameters are shown in Table 1, and the tensile strains are compared, as illustrated in Figure 11.

In Figure 11, the strains of the deflected geosynthetic calculated by the present method, the Villard and Briançon [22] method, and the Feng et al.’s [29] method are all close to the test results. Due to the lack of settlement data on the embankment surface, the calculations are all founded on the vertical slip surface assumption. Obviously, the present method yields the maximum tensile strains within the subsided area, and they are closer to the results observed five months after the end of the experiment. The Villard and Briançon [22] method does not consider the interface friction of geosynthetic and soil within the subsided area, so the peak tensile strain is the smallest and the length of anchorage is the smallest, too. Feng et al. [29] argued that ignoring the upper friction results in undervaluation of maximum geosynthetic strain, and their method obtains a relatively higher geomembrane effect. However, the overlying vertical load and its distribution significantly influence the geomembrane effect of the geosynthetic. Feng et al. [29] adopted the Terzaghi’s [8] fully mobilized arching theory to calculate the vertical stress on the geosynthetic layer, so the load obtained is less than that of the present method, as well as the corresponding geosynthetic tensile strains.

The above validation is based on the vertical slip surface condition. However, for the localized sinking of the base, the subsidence on top is enlarged and the width is greater than the base. In this case, the overlying load of the basal reinforcement is different from that of the vertical slip surface as previously assumed. Under a fixed subsidence at the fill base, when the inclined angle of slip surface decreases from almost 90° to ( is the internal friction angle of the fill), the overlying stress acted on the basal geosynthetic is variable as is shown in Figure 12.

This section compares the results of the analytical calculations by this study and a rubber water bag trapdoor test conducted by Miao et al. [35]. Different from the Villard and Briançon’s [22] test, Miao et al. [35] adopted a graded-drainage subsidence manner to evaluate the reinforcing mechanism of the geosynthetic. The advantage is that the changes of displacement and stress of embankment and the response of reinforcing behavior can be recorded stepwise. Although the cohesion of the embankment fill may have influence for the ease of comparison, it is also worth trying to analyze the change of overlying load. Detailed information about the experimental model is as follows: the fill height is 0.6 m with the gravity of 18.8 kN/m3 and the width and height of the rubber water bag are 1.0 m and 0.30 m, respectively. A single layer of the geosynthetic is installed at 5.0 cm on the rubber water bag, with more than 1.0 m of anchorage length. The ultimate tensile force is 42.3 kN/m with regard to the tensile strain of 10.9 %. The subsidence of the rubber bag is controlled by the water discharge.

Due to the relatively low embankment height, frictional arching models proposed by Terzaghi [8] and Handy [36] are selected to make a comparison. As shown in Figure 13, ranking the stresses with the increasing local subsidence, it is clear that the Terzaghi method obtains moderate stresses while the Handy method yields slightly larger ones, but all remain unchanged. With increasing subsiding width from the observation of embankment surface settlement, i.e., increasing inclined angles of the slip surfaces, the load acting on the geosynthetic decreases. Although there are some deviations from the measured data, the trend of the change coincides very well. It is very difficult to obtain the true reduction coefficient, so the reduction coefficient (Table 1) here is replaced by the ratio of geosynthetic deflection to the maximum settlement of the rubber water bag.

Based on three different local subsidence (see Figure 14), the tensile strains of the geosynthetic vary with different overlying loads, accompanied with the different subsidence on the fill surface. Through comparison, it can be seen that this method can predict the magnitude and distribution of tensile strains of the geosynthetic well according to different deflections. The analytical model is based on a symmetric condition, while the validation test is asymmetric, so the difference exists but is acceptable.

5. Conclusions

In this study, a semiempirical analyzing method, applying the developing arching effect based on mobilized shear stress along inclined slip surface, is presented to evaluate the geomembrane effect of the basal geosynthetic due to localized subsidence. Through comparison and discussion, the following conclusions are obtained:(1)Based on traditional definition of lateral Earth pressure coefficient, an improved definition method considering deflection of subsided soil within the subsided area is proposed. In derivation, the influence of mobilized shear stress along slip surface on arching and interface friction between a geosynthetic and fill can be quantitatively evaluated, which can be solved by an iterative method.(2)The present method is validated reasonable by two physical model tests, as well as results calculated by some current methods. For the overlying load on the geosynthetic, the present method obtains variable stresses according to variable deflections with a good coincidence. Considering the interface friction between the geosynthetic and soil and the incomplete arching effect, the present method obtains maximum tensile strains against the other two methods.(3)Influence of the loading pattern along the geosynthetic layer and cohesion of the fill material on the reinforcing mechanism of the geosynthetic are not discussed in this study. Further studies, including numerical simulations, simplifications, and more experimental tests, are strongly recommended to gain better understanding of the validity and application of the present method.

Data Availability

The experimental data used to support the findings of this study are included within the article. And other previously reported data supporting this article are from the outcomes of publicly published papers. Those prior studies (and datasets) have been cited at relevant places within the text as references. The authors agree to share the data of this paper and allow other researchers to verify the results of this article, replicate the analysis, and conduct secondary analyses.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge and appreciate the financial support provided by the National Natural Science Foundation of China (Project nos. 51508279; 51478054; 51508278) and the Jiangsu Provincial Natural Science Fund Projects (BK20150885). Moreover, the work was also supported by the Open Fund of Key Laboratory of Road Structure and Material of Ministry of Transport (Changsha University of Science & Technology, kfj180303), the Systematic Project of Guangxi Key Laboratory of Disaster Prevention and Structural Safety (2016ZDK005), and the Ministry of Housing and Urban-Rural Development of the People’s Republic of China (2018-K4-008).