Stability Analysis and Support Design Methods for Rock Foundation Pit with Combination of Structural Plane and Karst Cave

In this paper, the failure mode of the “deep foundation pit slope with high dip angle structural surface and cave development” is analyzed by numerical simulations, and an innovative numerical calculation method for calculating the stability of foundation pit slope based on the principle of upper limit analysis and the prediction method of rock slope surface is proposed. The lateral earth pressure acting on the support structure is further deduced, which can be used more widely than the traditional earth pressure calculation equations. Based on ﬁnite element analysis of a pole system, we developed an improved numerical calculation method of the force, and deformation of retaining piles is proposed to carry out the analysis of the deformation of structures under force. Our method avoids the need to artiﬁcially place earth pressure on either side of the embedded portions of the piles. Compared with the measured horizontal displacement of the retaining pile, this improved numerical calculation method can reﬂect the force and deformation law of retaining piles at certain construction stages. In this paper, the upper limit analysis method based on the stability analysis of slopes and earth pressure calculation equations is proposed, taking into account the working conditions of the karst cave in diﬀerent locations, and the minimum supporting force required can be calculated more accurately while predicting the range of slope sliding surface. An improved numerical calculation method of the retaining piles is proposed in this paper, and the calculation results can reﬂect the deformation law of the retaining piles during the excavation of the rock foundation pit, and the whole analysis and calculation process has certain advantages compared with the traditional method.


Introduction
In traditional support design concepts, the stability of rock pit slopes is generally considered to be relatively easy to control compared to the stability of soil pit slopes. However, in actuality, the support design of rock foundation pit slopes is often more difficult. e fundamental reason for this is that the stability of rock slopes is mainly controlled by the presence of unfavorable geological structures, such as caves, joints, and weak interlayers. us far, many researchers have conducted in-depth studies on various complex geological rock slopes [1,2]. Most of these studies have focused on natural large-scale landslides, but little attention has been paid to man-made slopes. e stability analysis of slopes is generally based on limit equilibrium methods [3][4][5], limit analysis [6][7][8], and numerical analysis [9][10][11][12]. Although researchers have considered the influence of unfavorable geological conditions in their analyses, these analyses have been limited to the general reduction of the overall stability coefficient due to the presence of unfavorable geological structures [13][14][15][16], without focusing on the effects of specific characteristics, such as the distribution of karst caves, on the overall stability.
is has led to a disconnection between the stability analysis of rock slopes and the actual operating conditions. For the treatment of adverse geological measures such as caves, scholars have proposed various methods [17][18][19], but considering the environmental impact, it is difficult to use traditional methods to construct many urban pits, except for large caves that are reinforced by grouting. And as a result, in rock support design, the elastic fulcrum method is usually adopted for the internal force analysis of foundation support structures [20][21][22][23]. e design process often relies on the designer's engineering experience regarding the laws of force and deformation of a pile, especially the calculation of the earth pressure in the embedded portion of the pile.
However, for complex rock slopes containing both karst caves and structural planes, using engineering experience as a design reference has a high degree of uncertainty in guiding the design. For this reason, some scholars have studied the distribution of lateral earth pressure on piles based on indoor tests [24,25] or field tests [26], and based on this, earth pressure models have been proposed to describe the relationship between earth pressure and displacement [27][28][29], but this part of the research is limited to soil slopes. erefore, to avoid subjectivity in the design process, it is necessary to develop methods for predicting the range of the sliding fractures and to support design calculation methods for different failure modes of sliding fractures for complex rock slopes containing karst caves and structural planes.
is paper proposes a stability and numerical calculation method for lateral earth pressure on piles based on the principle of upper limit analysis for different failure modes caused by karst caves and structural planes, and compares the calculation results with the traditional earth pressure calculation methods to determine its accuracy. At the same time, the traditional pile deformation calculation model will be improved and made more suitable for rock foundation pit deformation analysis, so that it can provide a reference for rock foundation pit design.

Failure Modes of Complex Rock Slopes
e stability of rock slopes with high dip angles is generally controlled by the structural planes. However, for rock slopes containing both karst caves and structural planes, the stability of the rock slopes is affected by the development of karst caves, and thus, the stability of the rock slopes is influenced by both structural planes and karst caves. From the perspective of the different positions and shapes of the caves and structures, the stability failure characteristics of such type of rock slopes are mainly divided into the three categories that follow Figure 1. (1) Sliding failure occurs along a particular structural plane, and the karst cave does not affect the overall stability of the slope. (2) e karst cave is penetrated, and the upper part of the fracture extends to the top surface of the slope, resulting in sliding failure characterized by a combination of curved and flat surfaces. (3) Karst cave penetration occurs between multiple structural planes, resulting in polyline-like sliding failure.
Of the above three failure modes, stability evaluation and lateral earth pressure calculation equations only exist for sliding failure along a single structural plane.
ere are currently no direct calculation equations for the other two failure modes.

Stability Analysis of a Rock Foundation Pit Based on the Upper Bound Analysis Method
For the failure mode due to penetration between karst caves described in Section 2, a simplified model was developed for foundation pit analysis with a single retaining pile. Based on the virtual power equation, the total internal energy dissipation rate of the system, _ W int , is equal to the internal energy dissipation rate of the unstable rock mass on the velocity discontinuity surface, _ W sint , plus the dissipation rate of the internal energy of the pile resistance force in the kinematically admissible velocity field, _ W pint . e balance equation is where _ W pint is the lateral pressure and the slip surface function of the rock, that is, the dissipation rate of the internal energy by the pile body's resistance force in a kinematically admissible velocity field, and _ W sint is the internal energy dissipation rate of the unstable rock mass itself on the velocity discontinuity surface.

Rock-Soil Lateral Pressure Calculation Method Based on
Complex Rock Slopes. In the upper bound analysis, it is possible to obtain the virtual power equation by setting the internal energy dissipation rate of the system equal to the power exerted by the external force [30].
where the internal energy represents the internal energy dissipation rate of the unstable rock mass and the velocity discontinuity surface. e external power includes the equivalent load generated by the self-weight of the rock mass and the virtual power made in the kinematically admissible velocity field. e final function of the sliding failure surface of the rock mass needs to be considered as a piecewise function. Specifically, the expression for the entire slope's surface is established in the following manner. e structure itself and the rock bridges between the structural planes, and the caves are expressed as linear functions, while the penetrated sections inside a cave and between the caves are expressed as logarithmic spirals.
As shown in Figure 2, the rock bridge between the karst cave and the structural planes is considered as individual segments. MN is defined as the upper structural plane, and GR is defined as the lower structural plane. e upper rock bridge DE is a straight line. e lower rock bridge and the karst cave are approximated by a logarithmic helix through the failure section, denoted by EF. OE is the starting point of the logarithmic spiral. D is the starting point of the rock bridge facture. e horizontal line DK, which crosses vertical excavation face AG at point K, is the constructed virtual working surface. DK is subjected to the gravity of the upper rock mass ACDK. e gravity function of segment DK corresponds to a uniform load with a length L and amount q. erefore, the force analysis of the entire analysis model can be divided into three parts: the logarithmic spiral surface composed of KDEO, OEF, and OFG. H C represents the depth of the lower structural plane on the slope surface. e angle between DB and DC is 45°+ φ/2. c 1 and φ 1 are the shear strength indicators for the rock mass. c 2 and φ 2 are the shear strength indicators for the structural plane. If k is the safety factor of the foundation pit, then the relationship between the strength parameters of the rock mass (c 1 , φ 1 ) and the reduced strength parameters (c m1 , φ m1 ) and the relationship between the strength parameters of the structural plane (c 2 , φ 2 ) and the reduced strength parameters (c m2 , φ m2 ) are as follows: e relationships between v 0 , v 1 , v 2 , and v can be expressed as (4) e total work done by the internal force is equal to the sum of the work done on straight line DE, the logarithmic spiral slip surface EF, and the straight line FG, and can be expressed as where E total is the total internal force of the work, 2E r is the total internal force of the work of the logarithmic spiral, and E EF and E FG are the internal forces of the work on EF and FG, respectively.
T � e θ 1 tanφ m1 , e total external force is the sum of the weight of the sliding soil, the earth pressure, the external force equivalent to the weight of the overlying soil, and the work done by the retaining force.

Advances in Civil Engineering
en, the general equation can be written as where A � c cos φ m1 · tan β.
L is a function of r 0 . In order to find the position in the cave between the two structural planes where penetration occurs, it is necessary to determine the minimum weight generated by the overlying rock mass, that is, the minimum value of L, defined as dL/dr 0 � 0.
When the retaining force F p � 0, the safety factor can be expressed as where c m1 and c m2 are the cohesion and internal friction angle of the rock mass, respectively, and φ m1 and φ m2 are the cohesion and internal friction angle of the structural plane, respectively. H C is the depth of the underlying structural plane from the slope's surface. When the rock mass slides along the fracture surface (defined as k ≤ 0), the lateral earth pressure F p in formula (9) can be expressed as e above equations can be programmed into MATLAB to constrain several parameters to find the optimal solution to the equation.
It should be noted that if the location of the cave is placed on the top surface of the foundation pit, the above equations can be used for the case of the slope sliding along a single 4 Advances in Civil Engineering structural plane in Figure 1(a). If the above equations do not take into account the shift in rock mass ratio of the virtual action surface, then they can be used when the slope is damaged by sliding from a curved surface to a plane in Figure 1(b). Applying the equations sequentially as a combination of multiple structural planes and karst caves, the equations can be used in the case of multiple lines of sliding failure between multiple structural planes occurs in the slope.

Comparison of the Safety Factor Results Obtained Using the eoretical Method and Simulations.
To verify the accuracy of the theoretical equations, karst cave penetration failure along parallel structural planes in a foundation pit slope was used as a test case. A calculation was carried out using MATLAB to obtain the safety factor change curve, which was then compared with the numerical simulation calculation results. e comparison and verification were carried out by substituting the numerical simulation parameters (c 1 � 300 kPa, φ 1 � 30°, c 2 � 144 kPa, φ 2 � 30°, c � 23 kN/m 3 , and H C � 25 m) into equation (13), and by setting the dip angle of the structural plane to 60°and 70°(the point labels of the karst cave model are shown in Figure 3). e curves of the relationship between the safety factor and the position of the karst cave's centroid are shown in Figure 4. At a dip angle of 60°, the load of the overlying rock mass at point 3 is 43.55 kN/m, L � 9.63 m. e factor of safety k is 1.24 for the theoretical calculation and 1.3 for the numerical simulation. e result of the theoretical calculation is slightly smaller than that of the numerical simulation in Figure 4(a). e result of the theoretical calculation tends to be more conservative.
e same calculation was carried out at point 1 for a dip angle of 70°. As shown in Figure 4(b), the result of the theoretical calculation is slightly larger than the numerical simulation result, with a maximum difference in the safety factor of 0.08. e calculation results are in general agreement with the simulation results.

Analysis of Support Design Based on the Upper Bound
Method. To better guide the design of engineering support, the safety evaluation equations for a foundation pit need to reflect the retaining forces required by the foundation pits for conditions where the factor of safety is lower than the specification. Figure 5 shows the relationship between the dip angle of the structural plane and the retaining force (the negative sign indicates that a retaining force is required), where the parameters are c 1 � 300 kPa, φ 1 � 30°, c 2 � 144 kPa, φ 2 � 30°, c � 23 kN/m 3 , h � 2.25 m, and k � 2.0. As the dip angle of the structural plane changes from 60°to 80°, the retaining forces required by the foundation pit change from large to small and then to large. For a constant length L, the retaining force required is maximum when the dip angle of the structure plane reaches 80°. In addition, the greater the distance from the cave to the slope of the foundation pit, that is, the closer the cave to the ground surface in the same structural plane, the greater the retaining force required. Similarly, in Figure 6, the curves of the retaining force versus the strength reduction factor also indicate that the higher the reduction factor, the greater the additional retaining force required, which is consistent with the actual conditions. Table 1 lists the common methods used to calculate the lateral earth pressure on a pile. Since a vertical slope is used, and the top layer of the slope is not overloaded, the calculation result of Coulomb's earth pressure (without considering the friction angle between the earth and the back of the wall) is the same as the calculation result of Rankine's earth pressure. e parameters used in the calculation are as follows: rock weight of 23 kN/m , rock mass friction angle of 45°, cohesion of 150 kPa, Poisson's ratio of 0.17, additional load on the top of the slope of 0, angle between the fracture plane and the horizontal plane of 60°, internal friction angle of the outward sloping structural plane of 15°, and the surface cohesion of 0 kPa at the structural surface of the outer slope. e method of calculating the earth pressure for sliding along outwardly sloping rigid structural planes is specific to rock slope engineering. e calculation result that considers the structural plane is closer to the actual range of damage. e calculation result obtained using our proposed upper bound method is smaller than that obtained using the method for an outward sloping rigid structural plane. is shows that when the operating conditions of karst caves are not considered, the calculation method along the outward sloping structural plane provides a slightly more conservative result. e reason for this is that the upper bound analysis is used to satisfy the equilibrium of the virtual powers based on static and dynamic conditions. e calculation result is the minimum remaining sliding force of the rock mass that satisfies the equilibrium equation, which is more stringent than the traditional static equilibrium equation. erefore, the calculated result for the slip along the outward sloping rigid structural plane satisfies the calculation equations of the upper bound method.

Advances in Civil Engineering
Rankine and Coulomb's earth pressure equations are generally used for loose structures. When the retaining wall is short, the calculations are conservative, while when the retaining wall is high, the calculations tend to be aggressive. e calculation result of the wedge analysis is closer to that of the upper bound analysis, but the disadvantage of this method is that it depends on the empirical value of the internal friction angle of the soil-rock interface.

Method of Calculating Pile
Deformation Control e pile displacement calculation method based on numerical analysis is a method for calculating the forces on retaining piles based on the upper bound method. In the analysis, the rock model only restricts the portion of the pile embedded in the rock. e range of the forces on the piles is determined by the actual excavation depth. e incremental rock-soil lateral pressure generated by the excavation can be  6 8 10 12 14 16 18 20 22 24 26 28

Numerical Method of Calculating the Pile Displacement.
e pile deformation calculations based on the finite element method for a pole system can be divided into two methods. e first method uses the total load method and does not take into account the construction process of the foundation pit. e second method adopts the incremental method to consider each phase of the pit construction. e calculations for both methods depend on the forces on the pile and the restraint forces of the soil around the pile that is below the excavation face. A rectangular parallelepiped rock mass of 120 m × 120 m × 100 m was modelled in the calculations using 3DEC distinct-element calculation software to simulate the actual stress state of the foundation pit by applying loads to the rock mass to replace the calculation of the earth reaction force in the embedded portion (Figures 7 and 8). e interaction between the pile and the rock mass is simulated by setting up two sets of nonlinear coupling springs in the normal and tangential directions, obeying the Mohr-Coulomb theory. e cohesive force and angle of internal friction on the interface are set at 0.8 times as high as the adjacent rock formation, and the normal and shear stiffnesses of the interface are set at 10 times the stiffness of the rock mass [31].
When the incremental method is used to consider the excavation process of a foundation pit, the rectangular parallelepiped model was divided into vertical layers, and each layer represented the thickness of each excavation until the designed excavation depth of the foundation pit was reached. In addition, a retaining pile was set at the center point of the model, with the top of the pile being flush with the top surface of the model. An example is shown in Figure 8. A retaining pile with a diameter of 1 m and a length

Calculation of Pile Deformation Based on Numerical
Analysis. Based on the above calculation principles, Table 2 shows the incremental earth pressure per 3 m layer for a 30 m excavation depth using Rankine's earth pressure calculation equation.
Substitution of the values in Table 2 into the model shown in Figure 8 yields the results shown in Figure 9. Prior to the excavation of the pit to a depth of 9 m, the displacement of the pile was very small, and the portion of the retaining pile below the excavation face was barely deformed. After the installation of steel supports, as the excavation depth increased, the upper half of the pile tended to deform toward the outside of the foundation pit, while the lower half of the pile appears to convexly deform toward the inside of the foundation pit, with the maximum displacement point exhibiting a downward trend as the excavation depth of the pit increased. When the excavation depth reached 30 m, the maximum displacement point of the deformation curve occurred at about two-thirds of the pile length. Figure 10 shows the pile displacement curves calculated using the traditional finite element method for a pole system. In the calculation, the rock mass model below the unexcavated surface was not defined, and the active earth pressure was used on the outside of the foundation pit, and the earth reaction force on the embedded portion is used for the inner side of the foundation [32]. It can be seen from Figure 10 that the calculated soil reaction forces produced excessive pile deformation prior to the installation of the steel supports, with the maximum deformation occurring at an excavation depth of 3 m, which is far from consistent with the actual situation. However, when the excavation of the foundation pit reached a depth of 21 m, the trend of the deformation was consistent with the calculation conditions and the rock mass constraints. In addition, when the excavation reached the bottom of the foundation pit (30 m), the maximum displacements reached by the two were basically the same, and the maximum displacement occurred at about two-thirds of the pile length.

Method of Support Design for Rock Slopes
Containing Karst Caves e rock slope stability analysis based on the upper bound analysis principle can predict the range of the slip surfaces. Provided that the approximate location of the cave's centroid can be determined via geophysical field testing methods, equation (12) and MATLAB can be used to calculate the minimum overburden load distance for cave penetration failure L min . When the horizontal distance between the actual position of the cave and the slope of the foundation pit is L < L min , penetration failure will not occur. en, the position of the karst cave in equation (13) is placed on the top surface of the foundation pit, which becomes a safety assessment of the overall slip along the structural plane. By using this method to assess each karst cave where penetration can occur, a final extent of sliding fracture of the rock mass can be calculated. Figure 11 presents a flowchart of this process.
Substituting the calculated safety factor k and rock mass parameters into equation (14), the resultant earth pressure within the range of the slip surface can be obtained. If the excavation process of the foundation pit is considered, and a virtual structural plane is set at the corresponding depth of each excavation layer, the increment of the earth pressure generated by each excavation layer can be calculated. By combining this with a numerical simulation based calculation method for the force and the deformation of the retaining pile, the retaining effectiveness of the pile can be evaluated.

Verification of an Engineering Case: The Foundation Pit of the Nansanhuanlu Station
e Nansanhuanlu Station of the Xuzhou Metro Line 3 is located at the intersection of the Nansanhuan Road and the Beijing Road. e main foundation pit of the station was constructed using the open-cut and bottom-up construction method. Retaining piles and horizontal internal supports (anchor cables) were used as the retaining structures. e excavation depth of the foundation pit was about 30 m, and the stratigraphic structure is dominated by limestone (Table 3). e section we studied is located on the west side of the foundation pit and belongs to the slope of the outward sloping structural plane (at a dip angle of 70°). Bored piles of Ø1000 mm × 1500 mm were used for this section. e first support is made of concrete, and the second and third supports are made of steel. An inclinometer tube was installed on the west side of the deepest intermediate mileage section. e monitoring positions are presented in Figure 12.
e inclinometer tube was damaged during the construction of the W2 pile and is therefore not discussed here.
Detection data from the seismic tomography survey (provided by Professor Song Lei's research group) show the caves (Figure 13). Karst cave 1 exhibits a longitudinal extension of the foundation pit, which is a massive cave with a length of 15 m.
In the calculations, a 5 m deep layer was used to simulate the construction. e first concrete support was constructed at a depth of −5 m. When the excavation reached −15 m, the first steel support was installed at −10 m. When the excavation went to −25 m, the second steel support was installed at −20 m, giving a final excavation depth of −30 m.
In the theoretical calculations, Cave 1 ( Figure 13) was transformed into a plane model, that is, without regard to the longitudinal shape and length. Cave 2 was not analyzed due 8 Advances in Civil Engineering to its small size. e rupture angles of the rock mass were taken as 60°(45°+ φ) and 70°. e parameters presented in Table 3 and the overlying rock layer of the karst cave were converted into an equivalent pressure of 122 kN/m. By substituting the measured distance from the center of the cave to the vertical excavation boundary of the foundation pit (5.8 m) and the safety factor k of 1.2 into equation (14), the remaining sliding force at an excavation depth of 15 m was calculated to be 185.64 kN, which is greater than zero. When the excavation depth is 25 m, the calculated remaining sliding force is 865.77 kN, which is also greater than zero. erefore, the upper and lower structural planes of the cave satisfy the conditions of overall sliding.
Substituting the actual horizontal distance between the cave and the excavation boundary of the foundation pit (L � 4.3 m) into equation (12) yields L min � 10.22 m, which is greater than the actual horizontal distance L, thus indicating that cave penetration failure may occur. Based on the above analysis results, the eventual slope of the foundation pit will suffer overall sliding failure. Figure 14 shows the slope displacement failure heat map calculated using the discrete element method. Karst cave penetration has occurred, but there is also relative sliding of the upper rock mass of the karst cave, which is consistent with the theoretical analysis result. Figures 15(a)-15(d) show a comparison of the theoretical value and the measured value for each excavation stage. Figure 15(a) shows that the displacement curve sliding along the 70°structural plane is similar to the deformation of piles W1 and W3 on the measured-slope curves above the excavation face. In Figures 15(b) and 15(c), the displacement curve sliding along the 60°structural plane reflects the deformation characteristics of piles W1 and W3 on the measured-slope curves above the excavation face. e parts below the excavation face deviate from the theoretical calculation results due to the rock mass structural planes and measurement errors but basically reflect the deformation characteristics of the part above the excavation face.  When the pit was excavated to 30 m, the sliding displacement curve along the 60°structural plane was consistent with the deformation of W4 in Figure 15(d). is is because the displacement of retaining pile W4 occurs at the inflection point of the edge of the foundation pit, and the early excavation is strongly influenced by the spatial effect.
In addition, it can be seen that, for a slope with a 70°dip angle, the slip area is dominated by the surface of the high dip structure until the excavation reaches 15 m, while after

Conclusions
In this paper, calculation equations for stability analysis of slopes based on the upper limit analysis and improved methods for calculating the deformation of retaining piles based on numerical analysis are proposed, respectively, for rock slopes developed with a combination of structural planes and caves, and the major conclusions can be summarized as follows: (1) Calculation equations for stability analysis of slopes based on the upper limit analysis can calculate the stability of sliding slopes caused by caves in three different failure modes separately, so that the range of possible sliding slopes in the foundation pit can be predicted, and the analysis covers a wider range. (2) Calculation equations for the remaining sliding force of a foundation pit under ultimate conditions are derived from rock slope stability analysis, which enables a safety factor to be set according to the importance of the project and the calculation of the minimum supporting force required to be provided by the supporting structure. (3) e results of the remaining sliding force calculation equations proposed in this paper are slightly smaller compared to the traditional earth pressure calculation equations, and the conditions satisfied are more stringent, which can be used as a reference for the value of the structure design.
e improved method of calculating the deformation of the retaining pile based on numerical analysis avoids the process of calculating the lateral Earth pressure in the embedded portion of the pile compared with the traditional elastic fulcrum method, and the results of the pile deformation calculation are more in line with reality.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request. Yansen Wang TB18220020B0@cumt.edu.cn.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
Conceptualization, formal analysis, data curation, and writing were performed by J. X; supervision and funding acquisition were done by Y. W. e above authors have read and agreed to the published version of the manuscript.

12
Advances in Civil Engineering