The joints have great influence on the strength of jointed rock mass and lead to the multiscale, nonhomogeneous, and anisotropic characteristics. In order to consider these effects, a new model based on a genetic algorithm is proposed for locating the critical slip surface (CSS) in jointed rock mass slope (JRMS) from its stress field. A finite element method (FEM) was employed to analyze the stress field. A method of calculating the mechanical persistence ratio (MPR) was used. The calculated multiscale and anisotropic characteristics of the MPR were used in the fitness function of genetic algorithm (GA) to calculate the factor of safety. The GA was used to solve optimization problems of JRMS stability. Some numerical examples were given. The results show that the multiscale and anisotropic characteristics of the MPR played an important role in locating the CSS in JRMS. The proposed model calculated the CSS and the factor of safety of the slope with satisfactory precision.
1. Introduction
Slope stability analysis is important for assessing the safety of earth dams and natural slopes. However, stability analysis of JRMS is one of the important problems for geological engineering. Locating the CSS of JRMS is the chief concern. Because the rock masses consist of a large number of joints, they have the multiscale, nonhomogeneous, and anisotropic material characteristics. These reasons lead to the complex of deformation and stress distributions in JRMS and make the analysis of JRMS difficult. As a powerful tool, GAs are used widely in the homogeneous slope stability analysis. Goh [1] employed a GA to search for the minimum factor of safety of homogeneous slope. Mccombie and Wilkinson [2] presented a simple GA to locate the CSS of a homogeneous layered slope. Zolfaghari et al. [3] proposed a simple GA to locate the critical noncircular failure surface of a homogeneous layered slope. Sun et al. [4] applied a spline curve in conjunction with a GA to locate the CSS in slope stability analysis. Sengupta and Upadhyay [5] presented a GA to search for the minimum factor of safety of a homogeneous slope. Li et al. [6] used a real coded GA to search for the noncircular CSS of a homogeneous layered slope. However, these methods cannot solve such problems in nonhomogeneous and anisotropic material. FEM, on the other hand, has been used for the analysis of slope stability. In contrast to the simplified techniques, FEM can solve the problem of the growth of inelastic zones with time and a complex loading sequence. Some tests have to be carried out according to the concept and application. Kim and Lee [7] used FEM to locate the CSS of a slope. Cho and Lee [8] proposed a FE flow-deformation coupled analysis program for solving unsaturated slope stability problems. Wang et al. [9] employed GA and FEM to locate the critical circular slip surface. Zheng et al. [10] used FEM and strength reduction techniques to calculate the safety factors for slopes. Others of researchers focused on the influence of joints on the material strength of jointed rock mass. Du et al. [11] analyzed the mutual influences of joints and rock bridges. Chen et al. [12] studied the influence of joints on the shape of the CSS in JRMS and described the shape of CSS in JRMS as noncircular. Zhang et al. [13] used the minimum joints frequency to locate the CSS of JRMS.
In this paper, based on the recent studies, a new model based on GA was proposed to search for the CSS in JRMS from its stress field. The studies on the influences of joints on multiscale and nonhomogeneous characteristics of jointed rock masses were made. The methods for calculating the factor of safety of slope based on the stress field, as calculated by FEM, were proposed. Then the GA was used to locate the CSS of JRMS. And then the numerical examples were analyzed to verify the proposed model.
2. The Algorithm for Locating the Multiscale CSS2.1. The Factor of Safety in the Method of Slices
The factor of safety F in the method of slices is expressed as
(1)F=∫ΓτfdΓ∫ΓτndΓ,
where τn is the mobilized shear stress at any point along the slip surface, τf is the shear strength at the corresponding point, and Γ is the length of the slip surface.
The mobilized shear stress τn and the shear strength τf can be calculated from the stress field which has been obtained by the FEM:
(2)τf=ci+σntanφi,ci=(1-ki)cr+kicj,tanφi=(1-ki)tanφr+kitanφj,σn=σx+σy2+σx-σy2cos2α-τxysin2α,τn=σx-σy2sin2α+τxycos2α,
where ci and φi are the cohesion of a material of ith segment of a slip surface and the friction angle of a material of ith segment of a slip surface, respectively, cr and φr are the cohesion and friction angle of rock bridges, respectively, ci and φi are the cohesion and friction angle of joints, respectively, σn is the stress acting normal to a specified slip surface, α is the angle of normal direction of slip surface to the horizontal plane, σx, σy, and τxy are the stresses on the slip surface, and ci and φi are the function of the MPR of ith segment of slip surface ki.
2.2. The MPR k
Figure 1 shows rock bridges in jointed rock masses which are present because of the nonpersistent nature of the joints. In order to calculate the decrease in strength of jointed rock masses in different directions, Wang et al. [14] define the MPR of a rock mass as the ratio of the joint network on the shear failure path when the jointed rock mass is sheared to the damaged state in a certain direction. Figure 2 shows that the MPR k is calculated by
(3)kβ0=∑JL∑JL+∑RBR,
where JL and RBR are the projection length of joints and rock bridges in the shear failure path, respectively. As shown in Figure 2, B and L are the projection width and projection length, respectively. β0 is the visual angle, which can be used to express the direction of the joints.
A natural joint system.
The failure path of joints and rock bridges.
From (3) and Figure 2, we can observe that the MPR k is related to the projection width B and projection length L. Thus the MPR k has the multiscale characteristics.
In the 2D joint network, there are many combinations of forms of the joints and rock bridges. In order to simplify the calculation, according to the research of Du et al. [11], it is assumed that there are three basic combinations of the joints and rock bridges, for which the shear resistances are calculated as Figure 3.
The three basic combination forms of the joints and rock bridges.
The two adjacent joints overlap but do not intersect as shown in Figure 3(a), and the shear resistant force R is expressed as
(4)R=τj1l1+σth+τj2l2,
where τj1 and τj2 are the sheared strength of the two joints, respectively. σt is tensile strength of the rock bridges. l1 and l2 are the length of the two joints along the shear plane, respectively. h is the length of the rock bridge along the vertical component of the shear plane and d is the length of the rock bridge along the shear plane.
The two adjacent joints are neither overlapping nor intersecting, as shown in Figure 3(b). The shear resistant force R is expressed as [15]
(5)R={τj1l1+σth+τj2l2,β≥θ,τj1l1+τbh+τj2l2,β<θ
in which
(6)τb=[σt(σt+σn)]1/2,θ=arctan(2τb/σn)2,
where τb is the shear strength of the rock bridge along the shear direction, θ is the angle between the failure plane of the rock bridge and the shear direction, and σn is the average value of the normal stress.
The two adjacent joints are intersecting, as shown in Figure 3(c). The shear resistant force R is expressed as
(7)R=τj1l1+τj2l2.
The dynamic programming method is used to search for the failure path of the joints and rock bridges included in the shear bandwidth and then calculate k. The coded program flow diagram is shown in Figure 4.
The program flow diagram for calculating k.
From the aforementioned analysis, we can observe that k has multiscale characteristics. When the bandwidth changes, k also changes.
2.3. The Stress Integration for a Slip Line within an Element
As represented in (1), the line integral of stress along a slip surface is required to calculate the factor of safety. The stress at any point within an element is
(8)σ=∑b=1nnNbσbnode,
where nn is the number of nodes in the element, Nb is a shape function for a nodal point b, and σbnode is the nodal stress.
With simpliication of the theoretical formulations, as shown in Figure 5, only one slip line segment is assumed within an element. And the stress integration for a slip line within an element can be calculated by
(9)∫P1P2σdΓ=L2∑i=1nintσ(P(si))Wi
in which
(10)P(si)=P1+P22+P2-P12si,
where nint is the number of integration points in the local coordinate system and si and Wi are the local coordinates of the ith integration point and the weight value of the ith integration point, respectively. P1 and P2 are intersection point vectors between the boundary of an element and a slip surface in the global coordinate system, and L is the length between P1 and P2.
Slip surface within an element.
A slip line segment
Local coordinate system
When the angle of the slip surface within an element is identified, the ki within an element can be also identified, and (1) can be calculated.
2.4. GA Applied in Locating the CSS
Basically, GA [16] is based on Darwin’s theory of survival of the fittest, which is based on the principle that solutions to a problem can be obtained through evolution. The algorithm starts with a set of possible solutions. The set of possible solutions is called the “population.” Each possible solution within the population is called a “chromosome.” Associated with each chromosome is a fitness value, which is found by evaluating the chromosome with respect to the objective function. Solutions from one population are taken and used to construct a new population so that the new population (offspring) will be fitter than the old one. This process is repeated until the best chromosome representing the optimum solution is produced or some termination criterion, such as a set number of generations, is reached. The fundamentals of GA are described in the following sections.
2.4.1. Fitness Function
The search for the CSS can be mathematically expressed as the minimization of the factor of safety F. And the problem can be written as
(11)minFsubjecttoα1≤α2≤⋯≤αi≤⋯≤αN,1≤i≤N,
where αi is the angle of ith segment of the slip surface, which is shown in Figure 6. N is the number of segments making up the slip surface.
The segments of the slip surface.
The constraint condition can be converted into a penalty function G(α1,…,αN) and substituted into the objective function, which can be rewritten as
(12)min(F+G(α1,…,αN))
in which
(13)G(α1,…,αN)=fpunFwhenαi>αi+1,
where fpun is the penalty factor. And fpun is adopted as a big value, which is far more than 1.
2.4.2. Initialization
The GA starts with a population of N possible solutions. Each solution is created by n binary encoding.
2.4.3. Selection
The selection operator is used to determine which chromosomes are chosen as parents that will create offspring for the next generation. In this study, the roulette wheel selection is used in conjunction with elitism. This usually involves retaining the best chromosome at each generation to ensure that the best chromosome is not lost if it is not selected in reproduction or if it is destroyed by crossover or mutation.
2.4.4. Crossover and Mutation
The single-point crossover operator is employed in this study and the randomly flipping mutation procedure is also used.
2.5. The Main Idea of Proposed Model
The program flow diagram for the proposed model is shown in Figure 7.
The program flow diagram for the proposed model.
3. Numerical Results and Discussions3.1. The Numerical Example for Calculating Multiscale Characters of k
The distribution of joints was shown in Figure 8. The parameters used in the numerical example were shown in Table 1. The length and depth of model were both 100 m. The average normal stress and the tensile strength of the rock bridges were 2.5 MPa and 2.0 MPa, respectively. The sheared strengths of the joints were all 1.95 MPa. The influences of the projection length and projection width on the MPR k were shown in Figures 9 and 10.
Occurrence geometry parameters of the joints.
Distribution form
Dip (°)
Dip angle (°)
Spacing (m)
Trace length (m)
Normal
Uniform
Normal
Negative exponential
Average value
129.25
28.62
2.71
3.7
Standard deviation
36.88
11.29
6.03
0.58
The distribution of joints.
The influence of the projection length on the MPR k (m).
The influence of the projection width on the MPR k (m).
Example 1 was used to describe the multiscale characteristics of the MPR k. Through observing the results in Figure 8 to Figure 10, they showed that the MPR k was related to the projection width B and projection length L. Thus the MPR k had the multiscale characteristics. Thus when calculating the CSSs of JRMS, it was necessary to consider the relationship between MPR k and the length of each segment of the slip surface.
3.2. Numerical Example: A Homogeneous Slope
This example was a homogeneous slope. The geometry and FE meshes of the slope were shown in Figures 11 and 12, respectively. The coefficient of cohesion c and friction angle φ of the slope were 9.8 KPa and 10°, respectively. The elastic modulus E and Poisson ratio ν of the slope were 200 MPa and 0.3, respectively. The density of the slope ρ was 1764 kg/m3. The distribution of von Mises stress in the slope was shown in Figure 13. The CSS of slope was shown in Figure 14. The factor of safety of slope F=1.3228. And other results of the CSS and factor of safety of slope were shown in Figure 14 and Table 2.
Statistics of the results of Example 2.
F calculated by Kang et al. [17]
F calculated by Sun et al. [4]
F calculated by proposed model
Min.
Max.
Mean.
SD.
1.3206
1.3210
1.3228
1.3206
1.3228
1.3215
1.172×10-3
Min.: minimum value; Max.: maximum value; Mean.: mean value; SD.: standard deviation.
The geometry of the slope.
FE meshes of the slope.
The distribution of von Mises stress of slope.
The CSSs for the slope [4, 17].
Example 2 was given for verifying the proposed model. The slope was homogeneous. And the results calculated by proposed model were compared with those calculated by other models in related literatures. In Figure 14, Kang et al. used an artificial bee colony algorithm (ABC) to locate the CSS of slopes and the factor of safety F was calculated by the Spencer method. Sun et al. used a Spline-based GA to locate the CSS of slopes.
Through observing the results shown in Figure 11 to Figure 14 and Table 2, they showed that the mean value, minimum value, and maximum value of F are 1.3215, 1.3206 and 1.3228, respectively. The standard deviation of the factor of safety F of the slope calculated by the proposed model and those used elsewhere was 0.1172%. And the value of F calculated by proposed model based on GA was greater than that calculated by the model based on ABC and the model based on Spline-based GA. The reason was mainly that the methods to calculate the factor of safety F and to locate the CSS of slopes are different. The proposed model used the FEM to analyze the stress field near the CSS, which considered the effects of the deformation on the stress in the slope. However, the values of F calculated by different models were close. These results showed that the proposed model could calculate the factor of safety F and the CSS of homogeneous slope well.
3.3. Numerical Example: A Homogeneous Layered Slope
This example was a homogeneous layered slope. The geometry and FE meshes of the slope were shown in Figures 15 and 16, respectively. The physical parameters for the layered slope were shown in Table 3. The distribution of von Mises stress in the slope was shown in Figure 17. The CSS of the slope was shown in Figure 18. The factor of safety of slope F=1.3418. And other results of the CSS and factor of safety of slope were shown in Figure 18 and Table 4.
Physical parameters for layered slope.
Layers
c (KPa)
φ (°)
ρ (kg/m3)
E (MPa)
ν
1
15.0
20.0
1900
200
0.3
2
17.0
21.0
1900
200
0.3
3
5.0
10.0
1900
200
0.3
4
35.0
28.0
1900
200
0.3
Statistics of the results of Example 3.
F calculated by Kang et al. [17]
F calculated by Li et al. [6]
F calculated by proposed model
Min.
Max.
Mean.
SD.
1.2917
1.3360
1.3418
1.2917
1.3418
1.3232
2.740×10-2
The geometry of the layered slope.
FE mesh of the layered slope.
The distribution of von Mises stress of layered slope.
The CSSs for a layered slope [6, 17].
Example 3 was also given for verifying the proposed model. This example was a homogeneous layered slope. And the results calculated by the proposed model were compared with those calculated by other models in related literatures. In Figure 19, Kang et al. used an artificial bee colony algorithm (ABC) to locate the CSS of slopes and the factor of safety F was calculated by the Spencer method. Li et al. used a real-coded genetic algorithm (RCGA) to locate the CSS of slopes.
The geometry of the earth rockfill dam.
Through observing the results of Figure 15 to Figure 18 and Table 4, they showed that the mean value, minimum value, and maximum value of F were 1.3232, 1.2917, and 1.3418, respectively. The standard deviation of the factor of safety F of the homogeneous layered slope calculated by the proposed model and those used elsewhere was 2.740%. And the value of F calculated by proposed model was also greater than that calculated by the model based on ABC and the model based on RCGA. The reason was mainly that the methods to calculate the factor of safety F of the CSS of the slope were different. The proposed model used the FEM to analyze the stress in the layered slope, which considered the effects of the deformation on the stress in the layered slope. However, the values of F calculated by different models were close. These results showed the proposed model could calculate the factor of safety F and the CSS of layered slope well.
3.4. Numerical Example: Homogeneous Layered Earth Rockfill Dam
This example was a homogeneous layered earth rockfill dam. The geometry and FE meshes of the earth rockfill dam were shown in Figures 19 and 20, respectively. The physical parameters for the layered earth rockfill dam were shown in Table 5. The distribution of von Mises stress in the earth rockfill dam was shown in Figure 21. The CSS for the earth rockfill dam was shown in Figure 22. The factor of safety of the earth rockfill dam F=1.3868. And other results of the CSS and factor of safety of earth rockfill dam were shown in Figure 22 and Table 6.
Physical parameters for layered slope.
Layers
c (KPa)
φ (°)
ρ (kg/m3)
E (MPa)
ν
1
49.0
29.0
2038
200
0.3
2
0.0
30.0
1764
200
0.3
3
7.84
20.0
2038
200
0.3
4
0.0
30.0
1764
200
0.3
Statistics of the results of Example 4.
F calculated by Kang et al. [17]
F calculated by Sun et al. [4]
F calculated by proposed model
Min.
Max.
Mean.
SD.
1.3776
1.3950
1.3868
1.3776
1.3950
1.3865
8.705×10-3
FE mesh of the earth rockfill dam.
The distribution of von Mises stress in the earth rockfill dam.
The CSSs for the earth rockfill dam [4, 17].
Example 4 was also given for verifying the proposed model. This example was a homogeneous layered slope. And the results calculated by proposed model were compared with those calculated by other models in related literatures. In Figure 22, Kang et al. used an artificial bee colony algorithm (ABC) to locate the CSS of slopes and the factor of safety F was calculated by the Spencer method. Sun et al. used a Spline-based GA to locate the CSS of slopes.
Through observing the results of Figure 19 to Figure 22 and Table 6, they showed that the mean value, minimum value, and maximum value of F were 1.3865, 1.3776, and 1.3950, respectively. The standard deviation of the factor of safety F of the earth rockfill dam calculated by the proposed model and those used elsewhere was 0.8705%. And the value of F calculated by proposed model was more close to the mean value. The reason was mainly that the methods to calculate the factor of safety F of the CSS of the earth rockfill dam were different. The proposed model used the FEM to analyze the stress in the earth rockfill dam, which considered the effects of the deformation on the stress in the earth rockfill dam. However, the values of F calculated by different models were close. These results showed that the proposed model could calculate the factor of safety F and the CSS of the earth rockfill dam well.
3.5. Numerical Example: JRMS
This example was a JRMS. The geometry and FE meshes of the slope were the same as those used in Example 2. The coefficient of cohesion c and friction angle φ of the rock bridges in the slope were 9.8 KPa and 10°, respectively, which were the same as those used in Example 2. The coefficients of cohesion c and friction angle φ of the joints of slope were 3.0 KPa and 7°, respectively. The elastic modulus E, Poisson ratio ν, and density ρ of the slope were the same as those used in Example 2. Thus, the distribution of von Mises stress of slope was the same as that in Example 2. The MPR k using different projection lengths was shown in Figure 23. The CSS of JRMS was shown in Figure 24. And the results of the CSS and factor of safety of JRMS were shown in Figure 24 and Table 7.
Statistics of the results of Example 5.
F for homogeneous slope
F for JRMS in projection length = 100 m
F for JRMS in multiscale projection length
1.3228
1.255
1.302
The rose diagram of k using a multiscale projection length.
The CSS of the JRMS.
Example 5 was used to analyze a JRMS. This example firstly used the model, which did not consider the multiscale characteristics of the MPR k. In other words, the MPR k did not change when the length of any segment of the slip surface was different. Then the material parameters of joints and rock bridges such as cohesion and friction angle remained constant. Then this example used the model, which considered the multiscale characteristics of the MPR k. In other words, the MPR k changed when the length of segment of slip surface was different. Then the material parameters of joints and rock bridges such as cohesion and friction angle also changed. And through observing the results of Figure 24 and Table 7, they showed that the multiscale characteristics of the MPR k had great influence on the CSS and factor of safety of JRMS. The factors of safety of JRMS calculated by different models were both less than that of homogeneous slope. And if the distribution of the MPR k changed much in different projection length, the CSS and factor of safety of JRMS would vary significantly. The factor of safety of the JRMS as calculated by the proposed model, which considered the multiscale characteristics of the MPR, was 3.75% more than that calculated by the proposed model, which did not consider the multiscale characteristics of the MPR. It showed that the multiscale characteristics of the MPR played an important role in the stability of the JRMS.
4. Conclusion
A model based on a GA for locating the multiscale CSSs in a statistical JRMS from the stress field was proposed. The proposed model could consider the influence of different sets of joints on the material strength of a jointed rock mass and the multiscale characteristics of the MPR k. The results showed that the MPR was related to the projection width and projection length. Thus the MPR had the multiscale characteristics and these played an important role in locating the CSS in a JRMS. The results also showed that the proposed model could calculate CSSs and factors of safety of slopes accurately.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (Grant nos. 51109029, 51178081, 51138001, and 51009020) and the State Key Development Program for Basic Research of China (no. 2013CB035905).
GohA. T. C.Search for critical slip circle using genetic algorithms20001731812112-s2.0-003370887110.1080/02630250008970282MccombieP.WilkinsonP.The use of the simple genetic algorithm in finding the critical factor of safety in slope stability analysis20022986997142-s2.0-003637781710.1016/S0266-352X(02)00027-7ZolfaghariA. R.HeathA. C.McCombieP. F.Simple genetic algorithm search for critical non-circular failure surface in slope stability analysis20053231391522-s2.0-1744437345710.1016/j.compgeo.2005.02.001SunJ.LiJ.LiuQ.Search for critical slip surface in slope stability analysis by spline-based GA method200813422522562-s2.0-3834912726510.1061/(ASCE)1090-0241(2008)134:2(252)SenguptaA.UpadhyayA.Locating the critical failure surface in a slope stability analysis by genetic algorithm2009913873922-s2.0-5374909491210.1016/j.asoc.2008.04.015LiY.-C.ChenY.-M.ZhanT. L. T.LingD.-S.CleallP. J.An efficient approach for locating the critical slip surface in slope stability analyses using a realcoded genetic algorithm20104778068202-s2.0-7925158474110.1139/T09-124KimJ. Y.LeeS. R.An improved search strategy for the critical slip surface using finite element stress fields19972142953132-s2.0-003142900110.1016/S0266-352X(97)00027-XChoS. E.LeeS. R.Instability of unsaturated soil slopes due to infiltration20012831852082-s2.0-003505244710.1016/S0266-352X(00)00027-6WangC. H.XiaX. Y.LiG.Genetic algorithm for searching for critical slip surface in soil slopes based on stress fields20044434254282-s2.0-2942546052ZhengH.SunG. H.LiuD.A practical procedure for searching critical slip surfaces of slopes based on the strength reduction technique2009361-2152-s2.0-5694908667010.1016/j.compgeo.2008.06.002DuJ. C.ChenZ. Y.MiH. L.WangX. G.ZhouJ. C.Determination of comprehensive shear strength for jointed rock masses in 3D condition using genetic algorithm and Monte-Carlo method20042313215721632-s2.0-4344678783ChenJ. P.LuB.GuX. M.FanJ. H.ChengY. H.Numerical simulation study on 3D comprehensive shear strength of fractured rockmass2006257146314682-s2.0-33748465317ZhangW.ChenJ. P.WangQ.Determination of the most dangerous sliding surface for fractured rock mass slope based on fracture frequency201220692696WangX. G.ChenZ. Y.LiuW. S.Determination of joint persistence and shear strength parameters of rock masses by Monte-Carlo method1992114345355LajtaiL. E.Strength of discontinuous rocks in direct shear196919221822-s2.0-001452012010.1680/geot.1969.19.2.218GohA. T. C.Genetic algorithm search for critical slip surface in multiple-wedge stability analysis19993623823912-s2.0-003323711210.1139/t98-110KangF.LiJ. J.MaZ. Y.An artificial bee colony algorithm for locating the critical slip surface in slope stability analysis201345220722310.1080/0305215X.2012.665451