Simulation Study on Strength and Failure Characteristics for Granite with a Set of Cross-Joints of Different Lengths

+e strength and failure characteristics for granite specimen with a set of cross-joints of different lengths were studied using PFC2D software. +e results show that when the included angle of α between the main joint and loading direction is 30° or 45°, no matter what the included angle of β betweenmain and secondary joints is, the main joint controls crack propagation and failure of granite specimen, which occurs the shear failure propagating frommain joint tips, and the corresponding uniaxial compressive strength is low. Meanwhile, the secondary joint is the key joint for crack propagation and failure at α of 0° and 90° except when β is 90°. +e granite specimen occurs the shear failure propagating from secondary joint tips. And, the shear failure crossing upper tips of main and secondary joints is found at α of 0° or 90° and β of 90°. +eir uniaxial compressive strengths are large. Also, the combined actions of main and secondary joints determine crack propagation and failure at α of 60° except when β is 90°. +e granite specimen occurs the hybrid failure, including shear failure propagating from main joint tips and tensile failure propagating from main and secondary joints center or secondary joint tips. And, when α is 60° and β is 90°, the granite specimen occurs the shear failure along secondary joint plane direction, and its uniaxial compressive strength is small. Generally, when α or β is a fixed value, the uniaxial compressive strength firstly decreases and then increases with the increase of β or α. Additionally, when α is 60° and β> 45°, the uniaxial compressive strength represents a decreasing trend.+e uniaxial compressive strength at α and β between 30° and 60° is generally small. Finally, the microdisplacement field distributions of granite specimen were discussed.


Introduction
Under the effects of geological tectonic stress field, there are a significant amount of native defects in a natural rock mass, such as joints, fractures, faults, various microinclusions, and pores.ese native defects play an important role in the mechanical properties of the rock mass, further influencing the stability of underground engineering [1][2][3].Among them, the joint is a common native defect.It is significant to study the strength and failure characteristics of jointed rock mass.
Many experimental and numerical simulation studies have been conducted on the strength and failure characteristics of natural jointed rock masses by many scholars, and they have obtained many interesting results.Li et al. performed the conventional triaxial compression tests on the marble specimens containing natural joints, and the effects of joint angle (between the joint plane and major principal stress) varying from 29.3 °to 56.0 °on the strength and failure model were studied [4].rough conventional triaxial compression tests on the shale, three failure modes were found by Niandou et al., which was affected by the included angle between joint plane and major principal stress, and confining pressure [5].Mao and Yang have carried out the uniaxial and triaxial compression tests of slate with single set of structural planes, and they established one-dimensional and three-dimensional intensity variation law of slate with single set of structural planes using single discontinuities theory [6].Prudencio and Jan performed the biaxial compression tests on the physical models of rock with nonpersistent joints, and the effects of joint spacing, joint angle, and joint arrangement on the characteristics of failure mode, strength, and deformation of rock mass were studied [7].By conducting the uniaxial compression tests on gypsum specimens with a set of preexisting open flaws, the influences of joint inclination angle and joint connectivity rate on uniaxial compression strength, elastic modulus, and stress-strain curves of rock mass were well investigated by Chen et al. [8].Using a new shear device, Gehle and Kutter studied the breakage and shear behavior of intermittent rock joints [9].Xiao et al. studied the deformation and strength of columnar jointed rock-like material made of the mixtures of plaster, cement, and water under uniaxial compression and proposed four types of failure modes of columnar jointed rock mass under uniaxial compression [10].Goldstein et al. studied the size effects on the strength of jointed rock mass under uniaxial compression [11].
e true-triaxial compression tests focused on the influences of orientation and the intermediate principal stress of joints on the rock mass strength were carried out by Reik and Zacas [12].Laijtai performed the direct shear test on the intermittent jointed rock mass, and the composite damage theory was put forward [13].Based on the equivalent rock mass (ERM) technique, Zhou et al. performed the triaxial compression simulation on the jointed rock mass using the PFC 3D software [14].And, they found that the deformation and failure behavior were affected by joint degradation, the microcrack formation in the intact rock, the interaction between two joints, and the interactions of microcracks and joints.Using FLAC software, the effects of joint angle on the peak strength, mechanical behavior, and pattern of shearband (SB) of the anisotropic jointed rock specimen [15] were studied.Using PFC 2D software, Yang and Huang carried out the Brazilian splitting test for the jointed rock specimen, and the effects of joint angle and notch angle on the tensile strength and failure mode of jointed rock specimens were detailed analyzed [16].Yang et al. studied the growth mechanism of three-dimensional (3D) cracks in rocks under tension as well as the strength and failure behavior of samples [17].Morgan et al. studied the cracking processes in Barre granite, including the fracture process zones and crack coalescence [18].Also, the effects of joint angle in coal on the failure mechanical behavior of a roof rock-coal combined body were studied by Yin et al. [19].
e aforementioned achievements have great important significances for understanding the strength and failure characteristics of jointed rock mass.However, most of the above studies are focused on the rock mass containing single joint, groups of parallel joints, or nonpersistent joints.And, in the actual rock mass engineering, the joints generally occur in the form of a cross.Due to the difficulty in sampling for the rock specimen with cross-joints, there are few studies on the strength and failure characteristics of rock mass with cross-joints.e strength and failure characteristics of the rock specimen containing two sets of cross-joints with equal lengths were studied by Kulatilake et al. [20] and Huang and Yang [21].At the same time, the cross-joints in rock mass may be of different lengths.
For these, in this paper, the uniaxial compression tests for granite specimens containing a set of cross-joints with different lengths were simulated by PFC 2D software.And, on the basic of variations of the included angles of α between the main joint and loading direction, and β between the main joint and secondary joint, the effects of cross-joints on strength and failure characteristics were analyzed.e above achievements are of great significances for understanding mechanical properties of jointed rock mass.

Uniaxial Compression Model for Granite
Specimen Containing a Set of Cross-Joints with Different Lengths 2.1.Particle Flow Code.e particle flow code (PFC) is an effective method to study the macromechanics problems of an analytic object (including construction and rock mass) at the microlevel [22][23][24][25][26], which has been widely applied in simulation tests of uniaxial compression, biaxial compression, and triaxial compression for the rock specimen [19,[27][28][29][30]. ere are two bonding models in PFC 2D software, including the contact bond model and parallel bond model [26][27][28][29][30].Among them, the parallel bond can transmit both force and moment between particles, which is well applied to simulate the compact material, such as rock material.In this study, we built the uniaxial compression model for the granite specimen containing a set of crossjoints with different lengths using the parallel bond.

Microparameters of Granite Specimen.
In the parallel bond model, the macromechanical properties of the simulation model are mainly determined by microparameters of particles, such as the contact modulus of particles, particle stiffness ratio, elastic modulus of parallel bonding particles, normal strength of parallel bonding particles, shear strength of parallel bonding particles, and ratio of parallel bonding normal stiffness to parallel bonding shear stiffness.However, these microparameters cannot be directly achieved by the laboratory tests.erefore, before running the uniaxial compression tests, a lot of numerical simulation tests under same conditions as laboratory tests should be conducted until that it can well reproduce the macromechanical properties measured in the laboratory test.And, the determinations of microparameters are a process of minimizing the errors between the numerical results and experimental results.
Due to the limitation of laboratory test condition, the microparameters of granite provided by Zhang et al. [31,32] were used to carry out numerical test.Firstly, Zhang et al. conducted the triaxial compression tests for the standard granite specimen [31].And, the mechanical parameters (shown in Table 1) of triaxial compression of the standard granite specimen were selected as the reference basis for determining the microparameters.en through the method of "trial and error" repeated check comparison, the microparameters of granite were well achieved, as shown in Table 2.

Advances in Civil Engineering
At the same time, the di erence between the test values of cohesive force and internal friction angle is larger than that of the simulated values.is is the result of circular particles used in the simulation model [32].Figures 1 and 2 show the comparison of stress-strain curves and failure modes obtained by experiment and simulation under triaxial test of con ning pressure of 6 MPa, respectively.In Figures 1 and 2, the simulated stress-strain curve and failure mode are in good agreement with the laboratory test of the real granite specimen.

Model Construction for Granite Specimen Containing a Set of Cross-Joints with Di erent Lengths.
A particle ow model of the granite specimen with a set of cross-joints with different lengths was established by radius extension, as shown in Figure 3. e model size is 70 mm × 140 mm.And, a total of 16459 particle samples were generated in the model.e minimum particle radius is 0.3 mm, and the ratio of largest to smallest radius of particle is chosen as 1.66.A set of crossjoints with di erent lengths were created through JSET command.For distinguishing the main joint and secondary joint, particles passing through the main joint plane are represented by red, while particles passing through the secondary joint plane are represented by blue.And, centers of the main joint, secondary joint, and granite specimen are coincident.e lengths of the main joint and secondary joint are 30 mm and 20 mm, respectively.Additionally, the included angle between the main joint and loading direction is α, and β is the included angle between the main joint and secondary joint.
A joint plane contact is de ned in PFC 2D software as a special contact that exists between particles that fall on opposite sides of the joint plane [33].According to the previous studies, the microparameters of the joint plane were weakened and generally set as very small values [19,[26][27][28][29].In this study, the microparameters of the main  joint plane were taken the same as that of the secondary joint plane.Among them, the friction coe cient was set as 0.1, and parallel bond normal strength and parallel bond tangential strength were all set as 0. Now, the uniaxial compression model for the granite specimen with a set of cross-joints with di erent lengths was built.For this work to be executed, the loading was performed by moving the upper and lower walls at a loading rate of 0.1 m/s.In PFC 2D software, the time step is about 4.2 × 10 −8 s/step.us, 0.1 m/s can be converted into 4.2 × 10 −9 m/step.at is, the upper and lower walls moving 1 mm require about 238095 steps.Also, the loading rate of 0.2 m/s is su cient to ensure that the rock samples are under quasistatic loading [34].
erefore, in this study, the loading rate of 0.1 m/s can ensure that uniaxial compression models are under quasistatic loading.

Simulation Test Conditions.
In this study, α and β were all taken as 0 °, 30 °, 45 °, 60 °, and 90 °, respectively.en, α and β were combined with each other, and there are 25 test simulation conditions, as shown in Table 3. e simulation test conditions of 1, 6, 11, 16, and 21 are single jointed granite specimens.

E ects of α and β on the Stress-Strain
Behavior.Figure 4 illustrates the stress-strain curves of granite specimens under a de nite value of α.And, Table 4 gives the simulation results, including the uniaxial compressive strength and elastic modulus.
In Figure 4, the stress-strain curves of granite specimens under di erent simulation test conditions are basically consistent, which can be divided into three stages, ranging from linear elastic deformation stage, nonlinear deformation stage, to postpeak strain softening stage.And, the included angles of α and β a ect the uniaxial compressive strength and axial strain values required for each stage.At the same time, the elastic modulus almost does not change with the included angles of α and β, indicating that the stress-stain curves are coincident in the linear elastic deformation stage.And, these are veri ed in Table 4.In Table 4, the elastic modulus of granite specimens is basically similar. is is principally because that there is no initial compression stage in PFC 2D software, which is di erent form the physical test [35].And, this phenomenon can be found in many simulation studies using PFC 2D software [19,[26][27][28][29].

E ects of α and β on the Strength Characteristics.
Figure 5 shows the variation regularities of the uniaxial compressive strength with α and β, respectively.
In Figure 5(a), when β is a xed value, the uniaxial compressive strength decreases in the rst place and then increases with the increase of α.When β is 0 °, the granite specimens are the single jointed specimens.e variation regularity of uniaxial compressive strength with α is consistent with that of the single jointed rock specimen, verifying the accuracies of the simulation.However, when β is 45 °or 60 °, the uniaxial compressive strength shows a decreasing trend at α of 90 °compared with that at α of 60 °.rough previous studies and macrofailure patterns of granite specimens (shown in Figure 6), the uniaxial compressive strength corresponding to the shear failure along the joint plane direction is relatively low.e granite specimen at α of 90 °and β of 45 °is mainly destroyed by the shear failure basically along the secondary joint plane direction.e hybrid failures occur in the granite specimen at α of 90 °and β of 60 °, including the shear failure basically along the secondary joint plane direction and tensile failure propagating from the lower tip of the secondary joint.us, their uniaxial compressive strengths decrease at a certain extent.
In    Advances in Civil Engineering main joint tips and tensile failure propagating from the joint center or the secondary joint tips.eir shear crack propagation paths are similar.At the same time, tensile crack propagation paths are not the same, especially the crack 1 and crack 2 as shown in Figure 6.In PFC 2D software, only when the stress intensity transmitted between the particles is larger than the parallel bonding forces between the particles, the microcrack can be generated.Propagation path of crack 1 is larger than crack 2, which means more external forces are needed to the overcome parallel bonding forces for microcracks formation.Correspondingly, the uniaxial compressive strength is larger.us, the uniaxial compressive strength at β of 45 °is larger than that at β of 60 °.  e in uences of α and β on the uniaxial compressive strength are summarized in Figure 7. Generally, when α and β are all between 30 °and 60 °, the uniaxial compressive strength of the granite specimen is generally small.

Failure Characteristic Analysis
4.1.Macrofailure Patterns of Granite Specimen.Figure 6 gives the macrofailure patterns under di erent simulation test conditions.e failure characteristics of granite specimens were concretely analyzed as follows.
4.1.1.α of 0 °or 90 °.When β is 0 °, the main joint and secondary joint are coincident.Now the granite specimens are single jointed specimens.And, their failures are the hybrid failure pattern, including the shear failure and splitting failure.When β ranges from 30 °to 60 °, the secondary joint is the key joint for crack propagation and failure of the granite specimen.e granite specimens are mainly destroyed by the shear failure propagating from the secondary joint tips.And, the secondary failure cracks propagating from the main joint tips occur in di erent degrees.However, when α is 90  No matter what the value of β is, the shear failure propagating from the main joint tips causes the failure of the granite specimen, illustrating that the main joint is the key joint for crack propagation and failure of the granite specimen.

α of 60 °.
When β is 0 °, the shear failure propagating from single joint tips is the main failure of the granite specimen, accompanied by tensile failure propagating from the joint center.And, when β are 30 °, 45 °, and 60 °, the failures of granite specimens are the hybrid failure patterns, including the shear failure propagating from the main joint tips and tensile failure propagating from the joint center or secondary joint tips.And, the tensile failure degree is larger than that at β of 0 °, illustrating that the secondary joint enhances the tensile failure.Additionally, the secondary joint is the key joint for crack propagation and failure of the granite specimen at β of 90 °. e granite specimen occurs shear failure along the secondary joint plane direction.

Microdisplacement Field Distributions of Granite Specimen.
e microdisplacement distribution is helpful to reveal micromechanism of crack propagation [15,18].Figure 8 presents the microdisplacement distributions of the macrotensile crack and macroshear crack.And, the icons of AA, BB, CC, DD, and EE correspond to AA, BB, CC, DD, and EE in Figure 6.
In Figure 8, there are three types of mesodisplacement eld distributions for macrotensile crack.Firstly, the particles move in the same direction, and then parts of them change the movement direction (shown in Figure 8(a) AA), or they move apart (shown in Figure 8(a) BB).us, the microcracks are generated.And, these two kinds of mesodisplacement eld distributions for macrotensile crack are basically consistent with research results studied by Zhang and Wong [36].e third one is that particles move in opposite directions for microcrack generation, as shown in Figure 8(a) CC.Moreover, particles located on both sides of macroshear crack move in the parallel and opposite directions.

Conclusions
According to the above analysis, the following conclusions were drawn: (1) e included angles of α and β a ect the uniaxial compressive strength of the granite specimen.Generally, when α or β is a xed value, the uniaxial  8 Advances in Civil Engineering compressive strength firstly decreases and then increases with the increase of β or α.At the same time, when α is 60 °and β > 45 °, the uniaxial compressive strength represents a decreasing trend.And, the uniaxial compressive strength at α and β between 30 °and 60 °is generally small.(2) e failure pattern of the granite specimen depends on the included angles of α and β.When α is 30 °or 45 °, no matter what the value of β is, the main joint controls the crack propagation and failure of the granite specimen.e granite specimen mainly occurs the shear failure propagating from the main joint tips.At the same time, the secondary joint is the key joint for crack propagation and failure of the granite specimen at α of 0 °and 90 °except when β is 90 °. e granite specimen is mainly destroyed the shear failure propagating from secondary joint tips.And, the shear failure crossing upper tips of the main and secondary joints is found in the granite specimen at α of 0 °or 90 °and β of 90 °.Also, the combined actions of the main and secondary joints determine the crack propagation and failure of the granite specimen at α of 60 °except when β is 90 °. e granite specimen failure is the hybrid failure pattern, including the shear failure propagating from the main joint tips and tensile failure propagating from the joint center or secondary joint tips.And, when α is 60 °and β is 90 °, the granite specimen occurs the shear failure along the secondary joint plane direction.(3) ere are three types of mesodisplacement field distributions for macrotensile crack.Firstly, the particles move in the same direction, and then parts of them change the movement direction or, they move apart.e third one is that particles move in opposite directions for microcrack generation.Moreover, particles located on both sides of macroshear crack move in the parallel and opposite directions.

Figure 1 :Figure 2 :
Figure 1: Comparison of stress-strain curves obtained by experiment and simulation under triaxial test.

Figure 5 (
b), the uniaxial compressive strength rstly decreases and then increases with the increase of β under a xed value of α.At the same time, when α is 60 °and β > 45 °, the uniaxial compressive strength represents a decreasing

Figure 3 :
Figure 3: Simulation model for the granite specimen with a set of cross-joints with di erent lengths.

Figure 7 :
Figure 7: 3D graph of uniaxial compressive strengths of granite specimens containing a set of cross-joints with di erent lengths.