Abstract

To calculate the required strength of a cemented backfill with high aspect ratio, the confirmation of lateral pressure is fundamental and needs to be determined first. As for the backfill with a high aspect ratio of height to length, the shape of the slip surface is not straight when in the active state due to the limited space, which is different from the general backfill. For this reason, a formulation of the slip surface with a curved shape and a lateral pressure calculation method based on this curved slip surface were proposed. The proposed equation of the slip surface is affected by the geometry parameters of the backfill, internal friction angle of the backfill, and the friction angle of the backfill-rock interface. Then, by the combination of the minor principal stress trajectory method and the horizontal slice method, an ordinary differential equation of stresses was established and then solved numerically. Finally, the method based on Mitchell’s three-dimensional limit equilibrium model was used to calculate the required strength of the cemented backfill. The calculated results were compared with previous studies and validated with numerical models. The results showed good consistency for the backfills with high aspect ratios.

1. Introduction

Cemented backfill mining is widely used to deal with tailings efficiently and ensure safety of underground stopes [1, 2]. The performance of the cemented backfill is influenced by many factors (such as cement type, addition of reinforcement, loading rate, multiphysical consolidation process [35], chemical shrinkage, and temperature) after its placement. And all these factors can influence the stability of the cemented backfill. But, from the mechanical view, the effects of these factors can all be combined in the strength parameters of cohesion and friction angle. Besides that, the backfill-rock interface is also important. Due to the roughness of the rock face [6, 7], the cohesion and friction developed at the backfill-rock interface function as shear resistance preventing the failure of the backfill [8, 9]. Considering the common two-step mining method of extracting and then backfilling, the cemented backfill body of first-step stope will have a free face when extracting the second-step stope. In a series of methods for designing the required strength of the filling body by the limit equilibrium method, the shear resistance force from lateral rock mass is a major concern in the process of evaluating the stability of the backfill body [1012]. Thus, the lateral pressure force which contributes to the friction resistance has an important influence on the safety factor of the backfill body.

After determining the cohesion and friction angle of the backfill-rock interface, how to determine the lateral pressure force of the filling body remains an important problem. The product of the pressure force and the friction coefficient of the backfill-rock interface is part of shear resistance force preventing the filling body from sliding to the open space. Overestimating this pressure can lead to overestimation of the stability of the filling body, which may result in a lower strength of the backfill. This may pose a potential hazard, such as the collapse of the backfill to the adjacent mining area, the liquidation resulted by the insufficiency of strength of the backfill, and the collapse of barricade. Ting et al. [13] simulated the backfilling body of vertical and inclined stope through laboratory experiments and considered that the earth pressure calculated from the active earth pressure state is closer to the experimental results compared with the static earth pressure state. Sobhi et al. [14] also pointed out that when calculating the lateral pressure, the active earth pressure state is more practical than the static earth pressure state which assumes a very large stiffness of the surrounding rock. The calculation of the lateral pressure of the filling body by taking the active earth pressure state is also adopted in [7, 11, 14]. There are two kinds of slip surface when calculating the design strength of the cemented backfill. One is the slip surface used to determine the pressure on the lateral rock. Another is the slip surface used for calculating the safety factor of the cemented backfill when in the limit equilibrium state. And these two are different.

However, in the abovementioned literature, the lateral pressure of the filling body is calculated based on the linear slip surface. In fact, for the filling body with high aspect ratio (the ratio of the height of free face to the length of free face > tan (45° + φ/2)), Rankine’s slip surface (namely a straight face with an inclination angle of π/4 + φ/2) cannot develop under active earth pressure state, and often a curved slip surface will occur [1519]. In the case where a linear slip surface cannot be developed, how to obtain a usable lateral pressure force calculation method remains a very important problem. Thus, this paper focuses on the rarely studied backfills with high aspect ratios. The core objective of the present study is to develop a strength design method for cemented backfills with high aspect ratios. To achieve this, a suitable calculation method for the lateral pressure of the filling body with high aspect ratio is required first. And then, a strength design formula for filling body in this case is developed.

2. Calculation of Lateral Pressure Force

As shown in Figure 1, a vertical backfill with geometric sizes of is considered. The physical experiments performed by Mitchell et al. [20] showed that the sliding body is a wedge on the slip surface when the filling body is unstable. And the backfill-rock pressure is the lateral pressure of the filling body in the unstable state. As mentioned above, the two-dimensional active earth pressure model is generally used to calculate the lateral pressure. The monitored data showed that the earth pressure distribution at the midline and the edge is very close [16, 21], so the determined earth pressure is generally considered to be the earth pressure distribution along the midline of lateral face. Moreover, for conservative design and simplicity, it is generally considered that the filling body is cohesionless when calculating the side pressure [12]. In fact, the difficulty in calculating earth pressure of cohesive filling is the same as that of a cohesionless one and can be eliminated by the translation of the normal stress coordinate of the Mohr stress circle.

2.1. Determination of the Shape of the Slip Surface

The formulation of a curved slip surface involves the specific geometry of the filling body, which means that it is more targeted to the filling body and closer to the actual situation. However, the disadvantage is that it needs to be calculated for each project, and the formulation of the curved slip surface is more complex than that of a straight line resulting in difficulty to obtain an analytical solution for earth pressure distribution. But, with the development of numerical calculation methods, it is very easy to numerically solve out the earth pressure. Therefore, using a curved slip surface to calculate the lateral pressure should be considered as a refined design method.

The slip surface shown in Figure 2 is under the active earth pressure state. Inspired by the experiments of Tsagareli [15] and Yang and Tang [19], an exponential expression of the curved slip surface was selected in this paper. The specific constraints of the formula are as follows: the slip surface passes through the bottom point and its diagonal vertex, that is, the points with coordinates of (0, 0) and (Lb, H) in the given coordinate system shown in Figure 2.

The equation of the slip surface is specifically written aswhere φ is the internal friction angle of the filling body; C is a constant associated with the geometry parameters of the given filling body, which can be obtained by the abovementioned constraint conditions of the slip surface. It should be noted that an analytical expression for C does not exist for the complexity of the formula of the slip surface. Thus, it is necessary to substitute certain parameters of the given filling body to obtain a specific equation of slip surface.

As shown in Figure 3, comparing the slip surfaces obtained from the experimental data of Yang and Tang [19], the exponential slip surface proposed in this paper can satisfactorily describe the slip surface morphology of the vertical filling body under active earth pressure.

Then, the angle of the slip surface changing with height can be written as

The width of the sliding body varying with height can be deduced as

2.2. Estimation of Stress State

Handy [22] has long noticed the existence of friction at the wall-backfill interface and that the direction of the principal stress in the filling body will be deflected at a certain angle anyway. Then, the method of calculating the earth pressure with the minor principal stress trajectory was proposed. As shown in Figure 4, the minor principal stress trajectory method fully considers the disturbance of the wall-backfill friction to the stress distribution state of the filling body. So, it has been widely used in the calculation of earth pressure, especially the arc-shaped minor principal stress trajectory [18, 23, 24]. For this reason, the horizontal slice method combined with arc-shaped minor principal stress trajectory is also adopted in this paper for lateral pressure calculation, as shown in Figures 4 and 5.

For a vertical rock face, the deflection angle of principal stress at the backfill-rock interface is [25]

As for the deflection angle of principal stress at the curved slip surface, the Mohr-stress circle at the slip surface drawn by the Mohr-Coulomb criterion shows that the angle between and is always . Meanwhile, if the angle between the slip surface of the horizontal differential layer and the horizontal plane is notated as , that is, the angle between and is , then the following expression of the deflection angle of principal stress at the slip surface can be obtained:

Therefore, for a curved slip surface, equations (6) and (7) can be obtained by using the arc-shaped minor principal stress trajectory hypothesis.

Thus, the following equations are derived:

2.3. Establishment of Equilibrium Equation

Figure 6 is a schematic diagram of the force analysis of the differential strip of the filling body considering the horizontal shear stress. Assuming that the backfill-rock interface and the slip surface are both in the limit equilibrium state defined by the Mohr–Coulomb criterion, the following force analysis can be performed.

Ty and Vy are the shear force and vertical force in the given height y, respectively. and are the shear force and normal force acting on the backfill-rock interface in the given height y, respectively. And ; . Ts and Ns are the shear force and normal force acting on the slip surface in the given height y, respectively. And ; . And the expressions for differential forces in Figure 6 can be listed as follows:

According to the force balance in horizontal direction, equation (10) can be obtained.

Dividing equation (10) by and eliminating the higher-order infinitesimal leads to the following equation:

Similarly, the equilibrium equation in vertical direction can be obtained as

After introducing and , equation (10) becomes

By , equation (13) is further written as

Similarly, the vertical balance equation can be rewritten as

Then, the ordinary differential equation for average vertical stress in the backfill can be obtained as equation (16) by combining equations (14) and (15).

The use of circular minor stress trajectory may also have other substitutes because the method of minor stress trajectory is just a semi-analytical method. Furthermore, the established ordinary differential equation cannot be solved analytically because several transcendental functions are within this function.

3. Strength Design of the Cemented Backfill

3.1. Calculation of Lateral Pressure Stress

The specific expressions of L (y) and α (y) in equation (16) are determined by the formula of slip surface; the expressions of (y) and (y) are determined by the stress distribution described by the arc-shaped minor principal stress trajectory. The distribution of average vertical stress with height can be obtained by solving equation (16) numerically. Then, the solution can be used to obtain the distribution of the lateral pressure.

3.2. Analysis on the Limit Equilibrium Model

As shown in Figure 1, the weight of the sliding body iswhere is the cohesion of the backfill-rock interface; ; P is the lateral pressure force acting on the filling body; and the sliding angle β is generally taken as , according to Mitchell’s experiments [20].

The safety factor of the sliding body is then calculated aswhere c is the apparent cohesion along the slip surface. If it is considered that FS = 1, , the cohesion c required for any given internal friction angle is

From the abovementioned expression, it can be seen that when other parameters are given and fixed, the larger the , the smaller the required cohesion. To ensure a conservative design, it is necessary to find a value that is accurate enough but little less than the actual value to achieve a safe design.

3.3. Process for Designing the Required Strength

As mentioned above, the calculation of the lateral pressure in the case of curved siding face is specific to the parameters of given filling body, that is, the geometry of filling body and the friction angle of the backfill-rock interface and the internal friction angle of the filling body. These data can be directly selected from typical data, such as studies of Liu et al. [21], Chen et al. [26], Falaknaz et al. [27], Jahanbakhshzadeh et al. [28], and Yang et al. [29]. And it is recommended to obtain approximate parameters in advance according to the results of the proportioning test in a laboratory. Also, one can select parameters within a certain range and then take the mean value of the lateral pressure force. If a conservative calculation is needed, it is suggested to take smaller values of φ and δ.

If the abovementioned lateral pressure calculation method is used to obtain a dimensionless lateral pressure distribution of , the lateral pressure force is calculated as

The abovementioned integrals also need to be numerically integrated, and the lateral pressure force is obtained and then substituted into equation (19) to obtain the required strength design parameters.

4. Verification and Discussion of the Proposed Method

In order to verify the reliability of the new calculation method for lateral pressure and the required strength of the backfill, the calculated lateral pressure and the designed strength of the filling body were verified. Lateral pressure on sides of the backfill was calculated and then submitted into equation (19) to evaluate the required cohesion. After that, the calculated cohesions were compared with those of previous studies and then used as cohesion parameters in FLAC3D models to evaluate the stability of the backfill.

4.1. Computation of Lateral Pressure

The solid data points in Figure 7 are the data from [29] and also [19]. The hollow data points are calculated using the proposed method. It can be seen from the comparison that the proposed method can predict the lateral pressure well. When the buried depth is relatively large, the calculated results are kind of different from the measured data, which can be viewed from the following perspectives.

Firstly, Wu et al. [30] pointed out that the experimental data given by Yang et al. [29] is consistent with similar experimental data in the shallow part, but there is a problem of larger earth pressure when the depth is large, which is not consistent with similar experimental data. Then, Yang [31] replied that this may be due to the fact that the experimental conditions are not fully guaranteed according to theoretical analysis, including the failure to ensure that the wall is completely in translational move and the fact that the friction of wall-backfill is complicatedly distributed along the wall. Thus, the actual experiment does not completely correspond to the experimental design that is theoretically hypothesized and desired. Secondly, as mentioned above, the purpose of this paper is to find a pressure that is close to the actual lateral pressure but less than the actual side pressure. This will not overestimate the lateral friction and antisliding force of the filling body and can guarantee a safe design as well as certain excess safety factor for dynamic disturbance.

Figures 8 and 9 depict the results of lateral pressure force and lateral friction calculated with the proposed method under the condition that δ=φ, respectively.

It is showed that the lateral pressure increases at an increasing rate with the increase in aspect ratio , and the smaller the friction angle, the larger the side pressure.

The reason why a smaller friction angle will lead to a larger side pressure under the same geometric conditions is that the friction angles of the filling body and the backfill-rock interface is synchronously reduced, and when a lower value is taken, the backfill is closer to the hydrostatic pressure state, and thus the surrounding rock is subjected to the maximum pressure in the hydrostatic pressure state. Although the lateral pressure P contributing to frictional resistance becomes larger under smaller values of φ and δ, the overall instability resistance decreases. This can be seen from the increase in the cohesion required for the backfill in the following section, which fully demonstrates that the stability of the backfill is not simply linear. If the δ used in the calculation of lateral pressure is larger than φ, a complex value with an imaginary part of P may occur, which means that δ should be no more than φ in the physical sense. Similarly, if the friction angle used to calculate the minimal cohesion is greater than the friction angle used to calculate the side pressure, then the resulting cohesion will also be wrong and there may exist a negative value, which indicates that the friction angle has an important influence on the stability of the filling body and need to be identified carefully.

4.2. Computation of the Required Cohesion

The calculated results of the required cohesion under the condition that δ = φ, FS = 1 and are shown in Figure 10.

The minimal cohesion increases with the increase in aspect ratio. Under the same geometric parameters, the larger the friction angle is, the smaller the required cohesion is. Similarly, under the same cohesion, a larger friction angle would allow for a larger aspect ratio. After the aspect ratio exceeds a certain threshold, the cohesion required seems to decrease slightly. However, the abovementioned results are based on the fact that the damage of the backfill is close to the wedge-shaped failure; for the filling body with a larger aspect ratio, whether the failure mode is still the wedge sliding or it is toppling over or of some other forms, need to be identified with the experiment, and then it is decided whether to use the calculated cohesion value or not.

4.3. Comparison of the Required Cohesion

In the case that the safety factor is 1 and the cohesion and friction angle of the rock-backfill interface are equal to those of the backfill, the calculated cohesion were compared with results from previous methods by Mitchell [20], Li and Aubertin [32], Li [9], and Li and Aubertin [33]. The comparison results are shown in Figure 11.

The curves show that the cohesions calculated with the same given parameters are lower than the results given by Mitchell et al. [20], Li and Aubertin [32], and Li [9]. And when the aspect ratio is less than 4, the calculated results are larger than those of Li and Aubertin, [33], but when the aspect ratio is greater than 4, the calculated results are slightly smaller than those given by Li and Aubertin [33]. These comparisons show that the design of the backfill strength proposed in this paper is effective and slightly conservative, but it has been greatly improved compared with the traditional method.

4.4. Verification with Numerical Modeling

The effective stress in the cemented backfill is evolutive [6, 34], which means that the strength parameters of the cemented backfill will develop from zero to a peak and then decline gradually. Furthermore, the stress boundary condition of the backfill is changing too [7]. These dynamic changes may come from many factors, such as disturbance from blasting or other dynamic loads, multiphysical consolidation process, and the attenuation effect from corrosive underground water and sulphate. However, the research on the influence of these factors is still ongoing [3, 4, 35] and cannot be quantified easily in laboratory modelling and numerical modelling. And to assess the validation of the proposed method, numerical modelling was proved to be an effective way [2, 11, 12, 21, 33, 3641]. To assess the stability of the backfill using 3D numerical simulations was proved to be effective [2, 2832]. Geometric and mechanical parameters used in numerical modeling were taken as typical values from previous studies [7, 8, 11, 12, 20, 36, 38]. The width and length of the backfill are fixed to 10 m; the heights of the backfill were taken as 16 m, 20 m, 24 m, 30 m, 40 m, and 50 m, a total of 6 different values. The size of the whole numerical model is varying from 330 m × 350 m × 400 m to 330 m × 350  m × 500 m according to the height of the backfill to avoid the boundary effect. Four different friction angle values of 25°, 30°, 35°, and 40° were used for each aspect ratio. Other parameters of the filling body are unit weight γ = 18 kN/m³, Young’s modulus E = 300 MPa, and Poisson’s ratio  = 0.3. Mechanical parameters of rock were taken as typical values from previous studies [7, 8, 11, 12, 36, 38]: unit weight γ = 27 kN/m³, Young’s modulus E = 30 GPa, and Poisson’s ratio  = 0.25. The four vertical boundaries of meshes were constricted using roller support. The bottom of meshes was fixed, and the top of the meshes was set as free face.

Some representative results of these 24 cases were selected and shown in Figures 12 and 13. Figure 12 depicts the contours of displacement to reveal the potential slip surface. Figure 13 illustrates the contours of strength to stress ratio to assess the safety factor of the backfill.

White-colored parts in Figure 12 stand for the backfill with relatively large displacement towards the void space. From its outline, the general shape of the sliding body may be obtained.

Figure 12 shows that the shape of potential sliding body is close to the wedge used in limit equilibrium analysis, which is similar to results of physical experiments [20]. In Figure 13, zones with a color closer to white indicate a larger safety factor. And if matched with the analytical method strictly, the safety factor of the backfill should be one. Figure 13 shows that the backfills are in the stable state in the cases calculated since the minimum of strength-stress ratio is one, a totally black region. In most regions of the backfill, the equivalent safety factor (strength-stress ratio) is greater than 1 but less than 3. The safety factor of a small part of the backfill is relatively large, and the maximum value can reach about 10. This shows that the safety factor of the major area is enough, but the margin is not large, which effectively reduces the cost of the cement on the basis of ensuring safety.

However, the stress state and boundary condition of the backfill are evolutive. And such evolution is resulted from many factors, including the development of the hydration process, the evolution of pore water pressure, dynamic loading on the backfill, and the evolutive boundary conditions of the backfill (flow of groundwater, attenuation of cohesion and friction between rock-backfill interface, and changing temperature). Furthermore, the strength parameters of the backfill can be weakened by corrosion of hydration products under the influence of corrosive underground water or sulphate. But, the quantification of these influences still needs further study. So, a safety factor usually larger than 1.5 was adopted in practice, which means the adopted strength data will not be in direct correspondence with the designed parameters. So, for a more precise and economical design, a more complicated analysis (taking time factor, dynamic factor, and multiphysical coupling effects into consideration) is about to be developed.

5. Conclusions

This paper focuses on the calculation method intended for the lateral pressure and strength design of the backfilling body with high aspect ratio (the ratio of the height of free face to the length of free face > tan (45° + φ/2)) for mines. The core objective of the present study is to develop a strength design method for cemented backfills with high aspect ratios.

A formulation of the slip surface and a lateral pressure calculation method based on the curved slip surface was proposed. The calculation method of the lateral pressure of the filling body with high aspect ratio was obtained. The lateral pressure can be calculated from given parameters of the filling body, including the height and length of the backfill, the friction angles of the backfill, and the backfill-rock interface. Then, the method based on Mitchell’s three-dimensional limit equilibrium model was used to calculate the required strength of the cemented backfill. The calculated results were compared with previous studies and validated with the numerical models. The results showed good consistency for the backfills with high aspect ratios.

Abbreviations

:Constant used in the exponential slip surface
C:Cohesion of the backfill
:Cohesion of the backfill-rock interface
:Gravity
:Total height of the backfill
:Buried depth from the top of the backfill
:Ratio of to
:Ratio of to
:Length of the filling body
: = width of sliding body at a given height y
:Normal stress at the rock-backfill interface
:Normal stress at the slip surface
P:Lateral pressure force of the backfill
R:Radius of circle of minor principal stress
:Tangential stress at the rock-backfill interface
:Tangential stress at the slip surface
:Vertical force on the horizontal backfill strip at a given height of y
:Width of the backfill
:Gravity of the sliding body
:Height of the horizontal backfill strip
:Inclination of the horizontal backfill strip
:Sliding angle when the backfill is about to fail
:Unit weight of the backfill
:Friction angle of the rock-backfill interface
:Rotation angle of major principal stress at the slip surface
:Rotation angle of major principal stress at the rock-backfill interface
:Horizontal stress at the slip surface
:Horizontal stress at the rock-backfill interface
:Normal stress at the slip surface
:Normal stress at the rock-backfill interface
:Average vertical stress on the horizontal backfill at a given height of y
:Major principal stress
:Minor principal stress
:Average horizontal shear stress on the backfill layer at a given height of y
:Internal friction angle of the backfill.

Data Availability

The data will be available once required.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was funded by the National Key R&D Program of China, grant no. 2017YFC0602903 and the State Key Program of the National Natural Science Foundation of China, grant no. 51834001.