Advances in Civil Engineering

Advances in Civil Engineering / 2020 / Article
Special Issue

Deep Rock Behaviour in Engineering Environments

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 8820158 |

Chao Peng, Qifeng Guo, Zhenxiong Yan, Minglong Wang, Jiliang Pan, "Investigating the Failure Mechanism of Jointed Rock Slopes Based on Discrete Element Method", Advances in Civil Engineering, vol. 2020, Article ID 8820158, 19 pages, 2020.

Investigating the Failure Mechanism of Jointed Rock Slopes Based on Discrete Element Method

Academic Editor: Fengqiang Gong
Received19 Jun 2020
Revised10 Aug 2020
Accepted24 Aug 2020
Published12 Sep 2020


This paper presents a comprehensive engineering method to investigate the failure mechanism of the jointed rock slopes. The field geology survey is first carried out to obtain the slope joint data. A joint network model considering the structural complexity of rock mass is generated in the PFC software. The synthetic rock mass (SRM) approach for simulating the mechanical behavior of jointed rock mass is employed, in which the flat-jointed bonded-particle model (FJM) and smooth joint contact model (SJM) represent intact rock and joints, respectively. Subsequently, the effect of microparameters on macromechanical properties of rock is investigated for parameter calibration. Moreover, the scale effect is analyzed by multiscale numerical tests, and the representative elementary volume (REV) size in the selected research area is found as 16 m × 16  m × 16 m. The microparameters of the SRM model are calibrated to match the mechanical properties of the engineering rock mass. Finally, an engineering case from Shuichang open-pit mine is analyzed and the failure process of the slope during the excavation process from micro- to macroscale is obtained. It has been found that failure occurs at the bottom of the slope and gradually develops upwards. The overall failure of the slope is dominated by the shallow local tension fracture and wedge failure.

1. Introduction

The rock mass is a geological body with fractures complexly distributed. These natural discontinuities such as joints and faults are ubiquitous in rock, which have an important impact on the mechanical behaviors of rock material and lead to the anisotropy, heterogeneity, and scale effect [1]. The mechanical properties and failure mechanism of the rock mass are always a major field of rock mechanics. The accurate description of geometrical characteristics of joint network is critical for the safe design, construction, and stability analysis in rock mass engineering. Moreover, investigating the mechanical properties and failure mechanism of jointed rock from micromechanical viewpoint can both enrich rock mass theory and guide the rock engineering practice.

Considerable research studies have been carried out on the joint network modelling and analysis methods of jointed rock mass [219]. Due to the complexity and uncertainty of joints system in the rock, the stochastic approach based on the theories of probability and statistics is widely used to describe joint geometries from limited measurements such as borehole logging and scanline [2], which is also the mathematical theoretical basis of numerical simulation of joint network. Beacher et al. [3, 4] assumed joints to be planar discs and considered the geometrical properties (e.g., spacing, size, and aperture) as independent random variables obeying certain probability distributions. Hudson and Priest [5], Priest and Hudson [6, 7], and Rives et al. [8] arrived at the conclusion that joint spacing usually followed a negative exponential, lognormal, or normal distribution depending on the degree of joint saturation in the network. Similarly, the properties of frequency [9], size [10, 11], and aperture [12, 13] were analyzed based on geological survey. Einstein and Baecher [14] processed measurement data using a rosette or stereogram to group the joints into different sets with their orientation feature by a uniform, normal, or Fisher distribution. Hammah et al. [1519] presented separately different data analysis methods for optimized and determined dominant orientations, which promoted the progress of statistical method of joint. Three-dimensional modelling and analysis numerical methods have been developed with the computer-aided technology in the past two decades [2022]. Xi et al. [23] investigated the initiation and propagation of precracks in rock. Einstein and Locsin [24] and Ivanova et al. [25] proposed 2D or 3D fracture simulation models, which promoted the application of modelling theory.

Particle-based models, which were originally introduced by Cundall and Strack [26] to simulate granular materials (e.g., soils and sands) and which gradually evolved to commercial code, i.e., particle flow code (PFC) [27], are now widely used for jointed rock mass. Conventional PFC modelling of intact rock considers the bonded-particle model (BPM) [28], which has two limitations, i.e., an unreasonable UCS/TS (compressive to tensile strength) ratio and low friction angle. Some options were suggested to solve these problems: changing the particle size distribution to reduce the porosity [29, 30] or changing particle shape into cluster [28] and clumps [31]. Potyondy et al. [32] created a new model, flat-jointed bonded-particle model (FJM), to overcome these limitations and used it to reproduce a realistic hard rock. To suppress such an artificial roughness, smooth joint contact model (SJM), which was independent of the arrangement of local particles in contact, was proposed to simulate the geometric-mechanical characteristics of joints [33]. Based on the two well-established techniques, BPM and SJM, Ivars et al. [33] and Pierce et al. [34] developed a synthetic rock mass (SRM) to represent the mechanical behaviors of jointed rock mass and obtained the scale effect, anisotropy, and brittleness that cannot be obtained using empirical methods of property estimation. The uses of SRM are mostly related to simulating crack propagation and the slip of rock mass under selected loading conditions, such as direct shear test [3539] and uniaxial and triaxial compression test [4042]. Some researchers employed SRM to study the effect of sample size on rock mass strength and derive the equivalent rock mass properties [4345]. SRM was also used to simulate the mechanical properties and failure mechanism of rock mass in engineering projects [4649]. However, due to complicated geometry of joints, SRM has rarely been used for the slope stability analysis of open-pit mine and the failure mechanisms of jointed rock slope during the ore excavation have not been fully understood.

The remainder of this paper is organized as follows. Section 2 reviews the main aspects and methodology of the components of SRM. In section 3, based on the statistic characteristics of joints, a numerical model of the jointed rock mass is built by PFC software. In section 4, the microparameters of FJM and SJM are calibrated to reproduce the macroparameters derived from laboratory tests. Then, multiscale SRM models are established to investigate the scale effect and determine the representative elementary volume (REV) size by numerical tests. Next, the microparameters are recalibrated to represent the behaviors of rock mass according to the REV size and Hoek–Brown method. In Section 5, a jointed slope of Shuichang open-pit iron mine is selected as an engineering case, and the failure mechanism and evolution process of the slope are analyzed from a micromechanical viewpoint.

2. Methodology

2.1. Synthetic Rock Mass (SRM) in PFC

The synthetic rock mass (SRM) approach has been implemented in PFC software, where particles are rigid spherical bodies joined by deformable contacts, to solve problems by using the explicit formulation of the distinct element method (DEM). This new technique can be used as a virtual laboratory to conduct numerical experiments to represent in a qualitative and quantitative mechanical behavior of jointed rock mass and has already made it possible to extend the volumes of rock at the scale of 10–100 m containing thousands of nonpersistent joints. In PFC software, the SRM brings together two well-established techniques: bonded-particle model for intact rock deformation and fracture and discrete fracture network (DFN) for joints. That is to say, the SRM model represents the jointed rock mass as an assembly of intact rock and an embedded DFN. The main components of SRM sample are shown in Figure 1. The following section presents the main aspects of each model.

2.2. Representation of the Intact Rock

FJM (flat-jointed bonded-particle model) was proposed based on the distinct element method to simulate the mechanical behavior of intact rock in PFC software by Potyondy [32]. This model can overcome the limitations of the standard BPM model and well explain the mechanical behavior of various aspects of rock, especially for hard rock. In FJM, the intact rock is represented as assemblage of rigid circular (in 2D) or spherical (in 3D) particles and piece of finite-size interface between two notional surfaces. A typical flat joint contact (Figure 2) can be deformed, partially damaged, and broken by the influence of external force. The particle can be considered as a skirted particle, which can provide interlocking and rotational resistance even after the interface breaking. Therefore, a more reasonable unconfined compressive to tensile strength (UCS/TS) ratio for hard rock can be obtained through FJM.

2.3. Representation of the DFN

The smooth joint contact model (SJM) simulates the behavior of a smooth interface, regardless of the local particle contact orientations along the interface [33]. This kind of contact is described as smooth because particle pairs joined by a smooth joint contact may overlap and “slide” past each other, instead of being forced to move around one another (Figure 3).

A smooth joint can be envisioned as a set of elastic springs uniformly distributed over a circular cross section, centered at the contact point and oriented parallel with the joint plane. The area of the smooth-joint cross section is given bywhere R(1) and R(2) are the particle radii and is the radius multiplier.

When the unbonded SJM model is specified, the incremental force calculation is performed depending on the elastic portion of the displacement increment .where kn and ks are the normal and shear stiffness, respectively, and (Fn)0 and (Fs)0 are the smooth joint normal force and shear force at the beginning of the timestep, respectively.

If , then . Otherwise, Fs = μFn, sliding is assumed to occur. While slipping, shear displacements produce an increase in normal force due to dilation:where μ is the friction coefficient and ψ is the dilation angle.

3. Construction of Complex Joint Network Model

3.1. Geometrical Characteristics of Joints

The natural joints often comprise complex network and dominate the mechanical behaviors of jointed rock mass. Therefore, it is important to obtain the geometrical characteristic representation of complex two or three-dimensional joint network for modelling jointed rock mass. However, it is impossible to directly obtain the in situ joint data deep in the rock. The details of 3D structures of joint network are usually inferred from the observation information of lower dimensional limited exposures like borehole logging, outcrop scanline, or window mapping. A large number of field survey data show that the joints in the rock mass mostly appear in sets, and it is necessary to identify the dominant sets of joints which highly affect the mechanical behavior of rock.

In this paper, a field survey was carried out in the north of Shuichang open-pit iron mine by using borehole logging and scanline methods. A total of 311 sets of joint are measured, and the data show that joints in this area are well developed. In addition, some joints have potential adverse effects on slope stability because the strike and dip direction of joints are almost the same with those of the slope. Statistical analyses are conducted on the data by using rose diagram or stereogram so that joints can be grouped into different sets with their orientations (Figure 4).

By using the hierarchical cluster method, the measured data were grouped into four dominant orientations, namely, A, B, C, and D (Table 1). And the number of joints in these four orientations is 248, accounting for 79.7% of the total.


Dip direction (°)40–70100–170210–240280–340
Dip angle (°)13–7519–8625–7719–69

The characteristic of distribution of joints can be described by geometric parameters such as dip direction, dip angle, and spacing. To provide necessary information for the subsequent establishment of joint network model, probability distribution modes of those parameters are studied with orientation characterized by a uniform, normal, or negative exponential distribution, and the histogram and fitting curve are shown in Figure 5.

The basic parameters and distributions for dominant orientations in Shuichang open-pit iron mine are obtained as shown in Table 2.


Dip direction (°)Mean value57131249318
Standard deviation8.3920.38.5820.76
Distribution lawP11P22P11P11
Dip angle (°)Mean value44544142
Standard deviation16.2514.1116.5315.44
Distribution lawP22P22P22P22
Half trace length (m)Mean value1.682.181.771.96
Standard deviation0.870.850.631.01
Distribution lawP33P33P33P33
Spacing (m)Mean value0.260.350.310.43
Standard deviation0.
Distribution lawP33P33P33P33

1P1 means uniform; 2P2 means normal; 3P3 means negative exponential.
3.2. Generation and Testing of Joint Network Model

The shape of the joint in three-dimensional model is regarded as planar discs of radius in PFC. In this case, the probability density function of the disc’s radius is calculated aswhere is the mean value of negative exponential function.

The mean radius of joint disc is given aswhere is the mean radius of joint disc and is the mean trace length.

In a relatively homogeneous area, the volume density is independent of direction and relatively constant and has the following relationship with the average linear density [50]:

Based on formulas (4) and (6), the volume density of the joint with the radius approximately complied with the negative exponential distribution can be obtained aswhere is the average linear density of joints, is the volume density of joints, d is the diameter, and is the average radius of joints.

According to formulas (4) to (7), the parameters of joint disc were calculated, as shown in Table 3.

IDMean radius of disc (m)Linear density of joint (crack/m)Volume density of joint (crack/m3)


The Monte Carlo method is a good tool for joint network modelling, which can approximately solve uncertain problems with a series of random numbers. In this paper, combined with the Monte Carlo method, the built-in language Fish provided by PFC is used to construct a three-dimensional joint network model. The random numbers of joint disc parameters in accordance with corresponding probability distribution are generated in a 10 m × 10  m × 10 m cube space. Based on the OPEN GL technique, the joint network model in PFC is shown in Figure 6.

Then, the reliability of this model is tested by comparing the distribution of main joints in the effective area between the two-dimensional section of the model and the window of field survey. Figure 7 shows the joint distribution from the field survey and numerical model. The three-dimensional joint network model has a good statistical similarity with the real field distribution in attitude, scale, and density. The joints obtained from numerical models can well describe the real distribution characteristics of joints.

4. Calibration of the Microparameters

The input microparameters in PFC software cannot be measured directly through conventional laboratory tests, which also have highly nonlinear relationships with the parameters such as Young’s modulus, uniaxial compression strength, and Poisson’s ratio in continuum numerical simulation. Therefore, parameter calibration is an essential part of simulation process in PFC to ensure that the microparameters can represent the macromechanical properties of rock well.

Before the calibration process, a series of conventional laboratory tests are carried out on Shuichang iron mine granite to obtain the mechanical parameters. The mechanical test equipment is shown in Figure 8. There are two sizes of experimental sample: 50 mm diameter and 100 mm height for compression tests and 50 mm diameter and 25 mm height for Brazilian tests, which are in good agreement with the sample sizes recommended by ISRM. Uniaxial compression tests and Brazilian tests are conducted at loading rates of 0.5–1.0 MPa/s and 0.1–0.3 MPa/s, respectively. Triaxial compression tests were stress path control with a constant loading rate of 0.5 MPa/s, and the confining pressures were 10, 15, 20, and 25 MPa. In addition, direct shear tests were conducted on samples to obtain the joint parameters. The sample size is 10  mm × 10  mm × 10 mm, with a constant horizontal shear rate of 0.01 m/s under different normal stresses of 300, 500, and 700 kPa. The mechanical parameters of the granite obtained from experimental tests are given in Table 4.


RockDensity, ρ (kg·m−3)2640Young’s modulus, E (GPa)55.23
Poisson’s ratio, ν0.183Uniaxial compression strength, σc (MPa)125.71
Cohesion, c (MPa)12.35Internal friction angle, φ (°)41.51
JointTensile strength, σt (MPa)6.65
Cohesion, cj (kPa)27.83Internal friction angle, φj (°)21.06

4.1. Calibration of the FJM

As is known, calibration of the microparameters in PFC is a complicated process because it is impossible to describe the relationship between macro- and microparameters with quantitative mathematical relation. The basic method to determine the microparameters is “trial-and-error.” The specific method is as follows: by keeping other parameters constant, investigate the effect of one single parameter on simulation results, compare the macroparameters from laboratory tests results, and then repeat the above steps time and again until the parameter could be in good agreement with the observed macroscale response.

In this paper, the microparameters of the FJM in PFC, which affect the macromechanical properties of numerical sample, should be determined as shown in Table 5. Those microparameters are calibrated to match Young’s modulus (E), Poisson’s ratio (ν), uniaxial compression strength (σc), tensile strength (σt), cohesion (c), and internal friction angle (φ) of the model with respect to the rock sample. The tensile strength is matched based on the numerical direct tension test, and the others are matched by using numerical uniaxial and triaxial compression tests. In addition, the other microparameters can be determined empirically by using the references.


dminMinimum particle radiusmm
dmax/dminParticle size ratio
Radius multiplier
NrNumber of elements in radial direction
NaNumber of elements in circumferential direction
EcEffective modulus of both particle and bondGPa
kratNormal-to-shear stiffness ratio of both particle and bond
μFriction coefficient
σbBond tensile strengthMPa
cbBond cohesion strength,MPa
ΦbFriction angle°

Based on previous research, there are six microparameters that have significant effects on the numerical sample in the FJM, i.e., Ec, krat, μ, σb, cb, and Φb. Therefore, a cylindrical sample model with a 50 mm diameter and a 100 mm height, which is in agreement with the rock sample size recommended by ISRM, is generated to simulate laboratory test (Figure 9). The sample model consisted of 14036 balls with radius varying from 1.00 to 1.66 mm.

In the FJM, the calibration process is time-consuming due to no one-to-one correspondence existing between macroparameters and microparameters. For example, the UCS, Poisson’s ratio, and cohesion occur in varying degrees of change with the change of microparameter cb. To improve the efficiency and avoid adjusting parameters blindly, a L25 (56) orthogonal design test and range analysis were used to investigate the response law of macroparameters and microparameters. The orthogonal test table and calculated macroparameters are given in Appendix (Tables 12 and 13). The results of range analysis are shown in Table 6, where R is the range value, to represent the influence of the microparameter on macroparameters, and a higher value indicates a higher ranking.

MacroparametersRange value and rankingMicroparameters

Rankingcb > σb > μ > krat > Ec > Φb
RankingEc > krat > Φb > μ > cb > σb
Rankingkrat >σb > cb > μ > Ec > Φb
Rankingσb > krat > μ > cb > Ec > Φb
Rankingcb > krat > μ> Ec > σb > Φb
RankingΦb > μ > cb > krat > σb > Ec

On the basis of aforementioned analytical results, a new improved trial-and-error method with high efficiency is used to calibrate the FJM. The final microparameters obtained are given in Table 7.

Microparameterdmin (mm)dmax/dminNrNaEc (GPa)kratμσb (MPa)cb (MPa)Φb (°)


Then, the microparameters were verified by using the numerical test. Figure 10 illustrates the stress-strain curves of uniaxial compressive test and failure modes from experimental result and numerical result. The stress-strain curves, both before and after the peak, are consistent with those of the experiment, and the failure patterns from the experiment and numerical test are quite similar.

The comparison of basic mechanical parameters between the experimental and numerical results is shown in Table 8. The errors of Poisson’s ratio, uniaxial compression strength, and tensile strength are 3.83%, 2.31%, and 3.16%, respectively. Young’s modulus and internal friction angle are just 0.94% and 0.53% larger than the experimental results, respectively. The value of cohesion from simulation yields the largest errors, which is 9.64% higher than the experimental value, but that is still less than 10%. Obviously, the properties obtained from the numerical tests with calibrated parameters are in good agreement with those from experiment.

Mechanical propertyE (GPa)νσc (MPa)c (MPa)φ (°)σt (MPa)

Experimental result55.230.183125.7112.3541.516.65
Numerical result55.750.190128.6113.5441.736.86
Error (%)0.943.832.319.640.533.16

4.2. Calibration of the SJM

The macroparameters of the SJM are usually determined by simulating the direct shear test sample of rock joint. In this paper, an unbonded SJM is selected, that is, the parameter “sj_state” is set to 0. As in the case of FJM calibration, a model with the same size as the experimental sample was generated in cube shear box, and a suitable set of microparameters of SJM was obtained based upon numerical direct shear tests matching the laboratory test response of the rock, as shown in Table 9.


sj_kn (N·m−1)Normal stiffness2.5 × 107
sj_ks (N·m−1)Shear stiffness2.5 × 107
sj_fricFriction coefficient0.5
sj-da (°)Dilation angle5
sj_stateJoint bond state0
sj_ten (MPa)Bonded tensile strength0
sj_bcoh (MPa)Bonded cohesion0
sj_bfa (°)Bonded friction angle0

The shear stress-shear displacement curves under different normal stresses from experimental and numerical direct shear tests are shown in Figure 11. The cohesion and internal friction angle of SJM can be calculated by using strength envelope according to those curves. The values of these two parameters are 30.58 kPa and 20.2°, which are 9.88% higher and 4.08% lower than the experimental value, respectively. Clearly, the calibration parameters can represent the macromechanical properties of rock mass well.

4.3. Scale Effect

Due to the randomness in the distribution of joints, scale effect can substantially influence the rock strength, namely, the strength of a region decreases with increase in region size up to the point at which a representative size is reached. Compared with the engineering rock mass, the laboratory samples are usually small and do not contain systematic joints which affect the rock strength. For this reason, the microparameters of engineering rock mass in the SRM model should be considered for the potential impacts of joints and should be calibrated to match a strength that reflects the size of a typical rock mass, instead of a core sample. Hence, a concept named “representative elementary volume (REV)” is introduced, which is the minimum size with which the rock mechanical properties can be treated as equivalent continuous [51, 52].

In order to investigate the scale effect and REV size of rock mass, seven SRM models with different sizes were established to conduct numerical uniaxial compression tests. These cubical models with side lengths of 1, 2, 6, 10, 14, 16, and 20 m, respectively, are shown in Figure 12. Moreover, the geometric centers of these models were fixed to make calculated results comparable. The calibrated microparameters in Tables 7 and 9 were used to represent macromechanical properties of rock.

Figure 13 shows the calculated stress-strain curve of multiscale models. As expected, the strength parameters decrease with the increase of model size, showing obvious scale effect.

The results are presented in Figure 14. Uniaxial compression strength decreases approximately linearly with the increase of size and begins to fluctuate around a stable value of 30 MPa after the size is larger than 16 m. The same goes for Young’s modulus; the values are close to 23 GPa when the size increases from 14 m to 20 m. Therefore, it can be concluded that the REV size is 16 m × 16  m × 16 m. In other words, the minimum size of a rock sample representing the property of engineering rock in Shuichang open-pit iron mine is16  m × 16 m × 16 m.

4.4. Final Calibration

In this paper, combined with the results from the laboratory tests, the mechanical parameters of engineering rock mass were estimated by using RocLab software based on the Hoek–Brown (H-B) method. The failure envelope curve of engineering rock mass in Shuichang open-pit iron mine is shown in Figures 15 and 16.

Subsequently, a series of sample models were established based on the REV size to conduct numerical tests. Making use of the method mentioned in Section 4.1, the microparameters of SRM were calibrated to match the mechanical parameters of engineering rock mass. The final microparameters of engineering rock mass in SRM are presented in Table 10.

Ec (GPa)kratμσb (MPa)cb (MPa)Φb (°)sj_kn (N·m−1)sj_ks (N·m−1)sj_fricsj_da (°)

Value181. × 1072.5 × 1070.55

The comparison of basic mechanical parameters between the Hoek–Brown method and numerical results is shown in Table 11. While the errors of the result are limited within 10%, it is considered to be a suitable set of microparameters to represent the macromechanical properties of engineering rock mass well.

Mechanical propertyE’ (GPa)σc’ (MPa)c’ (MPa)φ’ (°)σt’ (MPa)

H-B method14.099.1533.51341.730.154
Numerical result13.1510.0413.81443.360.151
Error (%)6.6719.7028.5683.9061.948

5. Engineering Case

5.1. Engineering Background

The field joint survey was carried out on the north side of the open-pit mine slope, No. II mining area of Shuichang iron mine. Rock mass quality in this area is declined due to the well-developed discontinuous nature of rock such as fault and joint, causing local instability of slope. Records show multiple rock fall and landslips occurred in this area, especially at elevation 44m–116 m. A 3D numerical model is developed using the PFC software in order to investigate the micromechanism of failure of jointed rock slope during excavation process.

5.2. Construction of Jointed Rock Model

The model takes the direction of perpendicular to the stope as the x-axis, extension direction of slope is y-axis, and vertical direction is z-axis. The lengths along the axis of x, y, and z are 140 m, 20 m, and 90 m, respectively. For obtaining reliable simulation results, the ratio of smallest characteristic length for the model to the median particle radius should be 50–100 [33, 34, 5355]. Considering the computational efficiency and accuracy, the minimum particle radius and particle size ratio are set as 0.3 m and 1.66, respectively. The four groups’ dominant orientations of joints obtained in section 3.1 were added into the model successively. The model before excavation is composed of 81886 particles and 54219 cracks, as shown in Figure 17. The height and slope angle of bench in this model are 15 m and 65°, respectively, which are entirely consistent with the facts. In what concerns the model boundary conditions, the horizontal displacements in the model’s vertical boundaries were fixed, along with all of the displacements in the lower boundary. The slope surface and higher boundary were free. Because the main rock of this area is granite, the strength parameters of rock mass and joint are selected according to Table 10.

5.3. Micromechanical Analysis of Slope Failure

In PFC, particles are rigid spherical bodies with bonded contacts representing intact rock. And those contacts will be broken in shear mode or tension mode when external force exceeds the bonding strength. Therefore, the number, mode, and propagation of crack can be used to analyze the failure evolution. In this paper, the model is excavated by four steps, one bench at a time. Subsequently, the failure mechanism during the excavation process was investigated from micromechanical viewpoint based on PFC software.

Figure 18 illustrates the relationship between crack number and timestep. It should be noticed that generation of cracks is concentrated in a short time (about 4000 timestep) after excavation, which indicates that excavation may cause instantaneous failure of slope. Moreover, with the progress of excavation, the number of cracks is increasing to a total of 14800, and the number of tensile cracks is much more than that of shear cracks at any timestep. It can be inferred that the failure is essentially caused by the dominated tensile microcracks.

Figure 19 shows the evolutionary process of microcracks during excavation, where the tensile crack and shear crack are shown in red and blue, respectively. In the first two steps, a small number of cracks appeared around the toe of the slope. From the third step, cracks begin to increase sharply in number and expand constantly inside the slope. At the end of the excavation process, the distribution of cracks reaches a depth of about 50 mm, and the tensile crack almost propagates through from the middle of the third bench to the toe of the fourth bench, forming a macro fracture zone.

Diagrams of velocity vector (Figure 20) and displacement vector (Figure 21) clearly show the distribution of potential unstable rock mass when the slope is unstable. The distribution is the same as that of microcrack. The velocity and displacement around the third and top of the fourth bench are larger than other regions.

Figure 22(a) shows the horizontal displacement of particles at the end of the simulation. Particles with large displacement are mainly concentrated in the middle and lower part of the slope, some of which have been separated from parent rocks and form a rock fall. And due to the cross cuttings of joints, the main failure mode of slope is wedge failure, being relatively complete in the upper and seriously damaged in the lower of the whole slope, which is consistent with the actual situation, as shown in Figure 22(b).

6. Conclusions

In this paper, the SRM model was established in PFC software to represent the mechanical properties of rock mass from Shuichang open-pit iron mine. Based on the prepared model, the failure mechanism during excavation process was investigated from the micromechanical viewpoint. The conclusions of this study are as follows:(1)statistical analysis method can act as a useful tool to group the dominant orientations of joints and determine the geometric parameters such as dip angle, size, and trace length. The joint distributions of Shuichang slope were obtained by the field survey and quantitatively described in the numerical model.(2)The effects of different microparameters on macromechanical properties of rock are different. Through the orthogonal experiment and variance analysis, the calibration and optimization of the parameters were accomplished. The mechanical properties including stress-strain curve and failure pattern from numerical results were verified with those from experimental results.(3)The existence of joints results in scale effect and anisotropic behavior of rock mass, and these properties tend to gradually weaken with the increasing rock block size. By carrying out numerical tests on multiscale SRM models, the representative elementary volume (REV) in the selected research area was obtained with the size 16 m × 16  m × 16 m. Then, the microparameters of the SRM model were calibrated to match the mechanical parameters of the engineering rock mass.(4)SRM is an effective method to analyze the failure evolution of jointed rock slope. The failure of the slope was dominated by the tensile microcracks between bonded particles. The microcracks primarily occurred at the bottom of the slope and gradually developed upwards. In addition, microcracks were mainly distributed on the shallow part of the slope. After excavation, the wedge occurred in the middle and bottom part of the slope (Appendix) (Tables 12 and 13).


Orthogonal Design Table and Results




Data Availability

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

Conflicts of Interest

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


This research was funded by the National Natural Science Foundation of China (grant no. 51034001).


  1. J. A. Hudson and J. P. Harrison, Engineering Rock Mechanics: An Introduction to the Principles, Elsevier, Oxford, UK, 2000.
  2. W. S. Dershowitz and H. H. Einstein, “Characterizing rock joint geometry with joint system models,” Rock Mechanics and Rock Engineering, vol. 21, no. 1, pp. 21–51, 1988. View at: Publisher Site | Google Scholar
  3. G. B. Baecher, N. A. Lanney, and H. H. Einstein, “Statistical description of rock properties and sampling,” in Proceedings of the 18th US Symposium on Rock Mechanics, pp. 1–8, American Rock Mechanics Association, New York, NY, USA, June 1977. View at: Google Scholar
  4. G. B. Baecher, “Statistical analysis of rock mass fracturing,” Journal of the International Association for Mathematical Geology, vol. 15, no. 2, pp. 329–348, 1983. View at: Publisher Site | Google Scholar
  5. J. A. Hudson and S. D. Priest, “Discontinuity frequency in rock masses,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 20, no. 2, pp. 73–89, 1983. View at: Publisher Site | Google Scholar
  6. S. D. Priest and J. A. Hudson, “Discontinuity spacings in rock,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 13, no. 5, pp. 135–148, 1976. View at: Publisher Site | Google Scholar
  7. S. D. Priest and J. A. Hudson, “Estimation of discontinuity spacing and trace length using scanline surveys,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 18, no. 3, pp. 183–197, 1981. View at: Publisher Site | Google Scholar
  8. T. Rives, M. Razack, J.-P. Petit, and K. D. Rawnsley, “Joint spacing: analogue and numerical simulations,” Journal of Structural Geology, vol. 14, no. 8-9, pp. 925–937, 1992. View at: Publisher Site | Google Scholar
  9. W. S. Dershowitz and H. H. Herda, “Interpretation of fracture spacing and intensity,” in Proceedings of the 33rd U.S. Symposium on Rock Mechanics (USRMS), pp. 757–766, Santa Fe, NM, USA, June 1992. View at: Google Scholar
  10. P. Davy, “On the frequency-length distribution of the san andreas fault system,” Journal of Geophysical Research, vol. 98, pp. 925–937, 1993. View at: Publisher Site | Google Scholar
  11. E. Bonnet, O. Bour, N. E. Odling et al., “Scaling of fracture systems in geological media,” Reviews of Geophysics, vol. 39, no. 3, pp. 347–383, 2001. View at: Publisher Site | Google Scholar
  12. C. A. Barton and M. D. Zoback, “Self-similar distribution and properties of macroscopic fractures at depth in crystalline rock in the cajon pass scientific drill hole,” Journal of Geophysical Research, vol. 97, no. B4, pp. 5181–5200, 1992. View at: Publisher Site | Google Scholar
  13. J. N. Hooker, J. F. W. Gale, L. A. Gomez, S. E. Laubach, R. Marrett, and R. M. Reed, “Aperture-size scaling variations in a low-strain opening-mode fracture set, Cozzette Sandstone, Colorado,” Journal of Structural Geology, vol. 31, no. 7, pp. 707–718, 2009. View at: Publisher Site | Google Scholar
  14. H. H. Einstein and G. B. Baecher, “Probabilistic and statistical methods in engineering geology,” Rock Mechanics and Rock Engineering, vol. 16, no. 1, pp. 39–72, 1983. View at: Publisher Site | Google Scholar
  15. R. E. Hammah and J. H. Curran, “Fuzzy cluster algorithm for the automatic identification of joint sets,” International Journal of Rock Mechanics and Mining Sciences, vol. 35, no. 7, pp. 889–905, 1998. View at: Publisher Site | Google Scholar
  16. B. Lu, X.-L. Ding, and A.-Q. Wu, “Study on method of orientation data partitioning of randomly distributed discontinuities of rocks,” Chinese Journal of Rock Mechanics and Engineering, vol. 26, no. 9, pp. 1809–1816, 2007, in Chinese. View at: Google Scholar
  17. B. Tokhmechi, H. Memarian, B. Moshiri et al., “Investigating the validity of conventional joint set clustering methods,” Engineering Geology, vol. 118, no. 3-4, pp. 75–81, 2011. View at: Publisher Site | Google Scholar
  18. L.-M. Xu, J.-P. Chen, Q. Wang et al., “Fuzzy c-means cluster analysis based on mutative scale chaos optimization algorithm for the grouping of discontinuity sets,” Rock Mechanics and Rock Engineering, vol. 46, pp. 189–198, 2013. View at: Google Scholar
  19. Y. Li, Q. Wang, J. Chen, L. Xu, and S. Song, “K-means algorithm based on particle swarm optimization for the identification of rock discontinuity sets,” Rock Mechanics and Rock Engineering, vol. 48, no. 1, pp. 375–385, 2014. View at: Publisher Site | Google Scholar
  20. C. Wei, Y. Li, W. Zhu, S. Li, S. Wang, and H. Wang, “Experimental observation and numerical investigation on propagation and coalescence process of multiple flaws in rock-like materials subjected to hydraulic pressure and far-field stress,” Theoretical and Applied Fracture Mechanics, vol. 108, p. 102603, 2020. View at: Publisher Site | Google Scholar
  21. J. Zheng, J. Deng, X. Yang, J. Wei, H. Zheng, and Y. Cui, “An improved Monte Carlo simulation method for discontinuity orientations based on fisher distribution and its program implementation,” Computers and Geotechnics, vol. 61, pp. 266–276, 2014. View at: Publisher Site | Google Scholar
  22. Y. Li, H.-P. Wang, W.-B. Cai et al., “Stability monitoring of surrounding rock mass on a forked tunnel using both strain gauges and FBG sensors,” Measurement, vol. 153, Article ID 107449, 2019. View at: Google Scholar
  23. X. Xi, X. Wu, Q. Guo, and M. Cai, “Experimental investigation and numerical simulation on the crack initiation and propagation of rock with pre-existing cracks,” IEEE Access, vol. 8, pp. 129636–129644, 2020. View at: Publisher Site | Google Scholar
  24. H. H. Einstein and J.-L. Z. Locsin, “Modeling rock fracture intersections and application to the Boston area,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 138, no. 11, pp. 1415–1421, 2012. View at: Publisher Site | Google Scholar
  25. V. M. Ivanova, R. Sousa, B. Murrihy, and H. H. Einstein, “Mathematical algorithm development and parametric studies with the GEOFRAC three-dimensional stochastic model of natural rock fracture systems,” Computers & Geosciences, vol. 67, pp. 100–109, 2014. View at: Publisher Site | Google Scholar
  26. P. A. Cundall and O. D. L. Strack, “Discussion: a discrete numerical model for granular assemblies,” Géotechnique, vol. 30, no. 3, pp. 331–336, 1980. View at: Publisher Site | Google Scholar
  27. Itasca Consulting Group Inc, PFC (Version 5.0) manual, Itasca Consulting Group Inc, Minneapolis, MN, USA, 2014.
  28. D. O. Potyondy and P. A. Cundall, “A bonded-particle model for rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 8, pp. 1329–1364, 2004. View at: Publisher Site | Google Scholar
  29. M. P. J. Schöpfer, S. Abe, C. Childs, and J. J. Walsh, “The impact of porosity and crack density on the elasticity, strength and friction of cohesive granular materials: insights from DEM modelling,” International Journal of Rock Mechanics and Mining Sciences, vol. 46, no. 2, pp. 250–261, 2009. View at: Publisher Site | Google Scholar
  30. X. Ding and L. Zhang, “A new contact model to improve the simulated ratio of unconfined compressive strength to tensile strength in bonded particle models,” International Journal of Rock Mechanics and Mining Sciences, vol. 69, no. 3, pp. 111–119, 2014. View at: Publisher Site | Google Scholar
  31. N. Cho, C. D. Martin, and D. C. Sego, “A clumped particle model for rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 7, pp. 997–1010, 2007. View at: Publisher Site | Google Scholar
  32. D. O. Potyondy, “A flat-jointed bonded-particle material for hard rock,” in Proceedings of the 46th US rock mechanics/geomechanics symposium, pp. 1510–1519, Chicago, IL, USA, June 2012. View at: Google Scholar
  33. D. M. Ivars, M. E. Pierce, C. Darcel et al., “The synthetic rock mass approach for jointed rock mass modelling,” International Journal of Rock Mechanics and Mining Sciences, vol. 48, no. 2, pp. 219–244, 2011. View at: Google Scholar
  34. M. E. Pierce, D. M. Ivars, D. O. Potyondy et al., “A synthetic rock mass model for jointed rock,” in Proceedings of the 1st Canada-US Rock Mechanics Symposium-Rock Mechanics Meeting Society’s Challenges and Demands, pp. 341–349, Vancouver, Canada, June 2007. View at: Google Scholar
  35. J. A. Vallejos, K. Suzuki, A. Brzovic, and D. M. Ivars, “Application of Synthetic Rock Mass modeling to veined core-size samples,” International Journal of Rock Mechanics and Mining Sciences, vol. 81, pp. 47–61, 2016. View at: Publisher Site | Google Scholar
  36. Y. Zhou, A. Misra, S.-C. Wu et al., “Macro- and meso-analyses of rock joint direct shear test using particle flow theory,” Chinese Journal of Rock Mechanics and Engineering, vol. 31, no. 6, pp. 1245–1256, 2012, in Chinese. View at: Google Scholar
  37. C.-C. Xia, S.-Y. Long, Z.-C. Tang et al., “Particle flow numerical simulation for shear behavior of rough joints,” Chinese Journal of Rock Mechanics and Engineering, vol. 31, no. 8, pp. 1545–1552, 2012, in Chinese. View at: Google Scholar
  38. M. Bahaaddini, G. Sharrock, and B. K. Hebblewhite, “Numerical direct shear tests to model the shear behaviour of rock joints,” Computers and Geotechnics, vol. 51, pp. 101–115, 2013. View at: Publisher Site | Google Scholar
  39. X.-R. Liu, Z.-Y. Deng, Y.-Q. Liu et al., “Macroscopic and microscopic analysis of particle flow in pre-peak cyclic direct shear test of rock joint,” Journal of China Coal Society, vol. 44, no. 7, pp. 2103–2115, 2019, in Chinese. View at: Google Scholar
  40. C.-C. Chiu, T.-T. Wang, M.-C. Weng, and T.-H. Huang, “Modeling the anisotropic behavior of jointed rock mass using a modified smooth-joint model,” International Journal of Rock Mechanics and Mining Sciences, vol. 62, pp. 14–22, 2013. View at: Publisher Site | Google Scholar
  41. Y. Zhou, S.-C. Wu, Y.-T. Gao, and A. Misra, “Macro and meso analysis of jointed rock mass triaxial compression test by using equivalent rock mass (ERM) technique,” Journal of Central South University, vol. 21, no. 3, pp. 1125–1135, 2014. View at: Publisher Site | Google Scholar
  42. M. Bahaaddini, G. Sharrock, and B. K. Hebblewhite, “Numerical investigation of the effect of joint geometrical parameters on the mechanical properties of a non-persistent jointed rock mass under uniaxial compression,” Computers and Geotechnics, vol. 49, pp. 206–225, 2013. View at: Publisher Site | Google Scholar
  43. K. Esmaieli, J. Hadjigeorgiou, and M. Grenon, “Estimating geometrical and mechanical REV based on synthetic rock mass models at brunswick, mine,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 6, pp. 915–926, 2010. View at: Publisher Site | Google Scholar
  44. Y. Zhou, L. Wang, J.-F. Ding et al., “Particle flow code analysis of multi-scale jointed rock mass based upon equivalent rock mass technique,” Rock and Soil Mechanics, vol. 37, no. 7, pp. 2085–2095, 2016, in Chinese. View at: Google Scholar
  45. J. A. Vallejos, K. Suzuki, A. Brzovic et al., “Characterization and synthetic simulations to determine rock mass behavior at the El Teniente mine, Chile. part II,” in Proceedings of Caving, pp. 1–9, Santiago, Chile, June 2014. View at: Google Scholar
  46. C.-C. Ma, T.-B. Li, H.-L. Xing et al., “Brittle rock modeling approach and its validation using excavation-induced micro-seismicity,” Rock Mechanics and Rock Engineering, vol. 49, no. 8, pp. 3175–3188, 2016. View at: Publisher Site | Google Scholar
  47. M. E. Pierce and C. Fairhurst, “Synthetic rock mass applications in mass mining,” in Harmonising Rock Engineering and the Environment - Proceedings of the 12th ISRM International Congress on Rock Mechanics, pp. 109–114, London, UK, October 2011. View at: Google Scholar
  48. J. Hadjigeorgiou, K. Esmaieli, and M. Grenon, “Stability analysis of vertical excavations in hard rock by integrating a fracture system into a PFC model,” Tunnelling and Underground Space Technology, vol. 24, no. 3, pp. 296–308, 2009. View at: Publisher Site | Google Scholar
  49. Q.-H. Li, J. P. Latham, and C.-F. Tsang, “The use of discrete fracture networks for modelling coupled geomechanical and hydrological behavior of fractured rocks,” Computers and Geotechnics, vol. 85, pp. 151–176, 2017. View at: Google Scholar
  50. F.-Q. Wu, S.-J. Wang, S.-W. Song et al., “Statistical methods and theories in rock mechanics,” Chinese Science Bulletin, vol. 38, no. 15, pp. 1345–1354, 1993, in Chinese. View at: Google Scholar
  51. Y. Zhou, G. Han, S.-C. Wu et al., “Meso failure mechanism of rock mass and slope with intermittent joints,” Chinese Journal of Rock Mechanics and Engineering, vol. 35, no. 2, pp. 3878–3889, 2016. View at: Google Scholar
  52. S.-C. Wu, Y. Zhou, L.-L. Gao et al., “Application of equivalent rock mass technique to rock mass engineering,” Chinese Journal of Rock Mechanics and Engineering, vol. 29, no. 7, pp. 1435–1441, 2010, in Chinese. View at: Google Scholar
  53. X.-T. Yin, Y.-N. Zheng, and S.-K. Ma, “Study of inner scale ratio of rock and soil material based on numerical tests of particle flow code,” Rock and Soil Mechanics, vol. 32, no. 4, pp. 1211–1215, 2011, in Chinese. View at: Google Scholar
  54. X.-B. Ding, L.-Y. Zhang, H.-H. Zhu et al., “Effect of model scale and particle size distribution on pfc3d simulation results,” Rock Mechanics and Rock Engineering, vol. 47, no. 6, pp. 1–18, 2013. View at: Publisher Site | Google Scholar
  55. S. Wu and X. Xu, “A study of three intrinsic problems of the classic discrete element method using flat-joint model,” Rock Mechanics and Rock Engineering, vol. 49, no. 5, pp. 1813–1830, 2016. View at: Publisher Site | Google Scholar

Copyright © 2020 Chao Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.