Experimental and Numerical Study of the Failure Behavior of Intermittent Rock Joints Subjected to Direct Shear Load

Two series of intermittent rock joints containing three joints arranged along the central shear axis were considered in this study. The failure behavior under direct shear loads was investigated by means of both physical tests and numerical simulations. The cracking behavior was found to be distinctly associated with the joint arrangement. Several types of main and secondary cracks were identified. The variation trends of the crack initiation stress ratio with inclination angle were analyzed and found to be partly different for the two series of intermittent joints. The whole fracturing process was characterized by three phases. Not all samples have to experience all three phases. The second phase is alternative and can be reflected by the shearing curve. Hence, two types of shearing curves, including single and double peaks, were identified. The double peak is due to the extrusion or sawteeth cutting in the second phase. Moreover, the numerical micromechanical analysis was performed to explain the shear behavior using the contact force and microcrack within the specimen. Based on the numerically measured local stresses, maximum and minimum principal stresses around the middle joint at crack initiation stress and peak shear stress were analyzed.


Introduction
Discontinuous joints are widely known to exist within natural rock masses.Types of joints play a dominant role in the mechanical behavior of jointed rock such as the shear or compressive behavior.In the past decades, many approaches including the laboratory test, numerical simulation, and theoretical analysis, have been conducted to get a better understanding of the failure behavior of jointed rocks subjected to various loadings [1][2][3][4][5][6][7].Deep knowledge of rock mass failure characteristics benefits the engineering construction, for instance, mining roadway, tunnel, hydropower dam, and so on, so the relevant investigations still attract intense interest.
Recently, for the discontinuous joint (intermittent joint) in the solid rock under a compressive loading, the failure behavior has attracted much attention in the literature [8][9][10][11][12][13][14].e laboratory test is known to be the most widely used method to investigate cracking behavior.In general, open or closed joints were precast in rock-like materials by inserting thin sheets [8][9][10][11][12].In some studies, the hydraulic cutting technology was adopted to produce open joints in the real rock [15].For a rock block containing intermittent joints, the geometrical parameters, joint orientation, and length are considered as factors directly influencing cracking.Abundant observations show that the crack often initiates at the joint tip, and the crack initiation, propagation, and coalescence are dramatically dependent on the joint arrangement [8,9,[15][16][17][18][19][20].Tensile and shear cracks are two basic crack initiation modes, but the propagation and coalescence trajectory is influenced by the properties of neighboring joints such as joint spacing and the persistent ratio.Many kinds of combinations of influential factors incorporating joint length, orientation, spacing, and number were further considered.erefore, the majority of crack patterns and failure modes were derived from different cases.In addition to the summary of cracking types, the relevant mechanism was also theoretically or numerically revealed [21,22].Obviously, shear and compression are two quite different ways of loading.e shear plane mainly bears the shear load in shear, while the whole block bears the compressive load in compression.For the rock block with the same joint arrangement, the observed failure behavior under shear loading may be different from the failure behavior under compression.During shearing, the failure characteristic is dependent mainly on the joints close to the shear plane.Moreover, cracking features are related to the joint geometrical parameters.Using a specially designed shear device, Gehle and Kutter [23] conducted a series of direct shear tests on intermittent joints, and three phases, including crack initiation, friction, and dilation sliding, were identified.e shear failure behavior of the intermittent joint depends significantly on the joint orientation.Gerolymatou and Triantafyllidis [24] prepared gypsum specimens with different joint orientations and carried out a largeshear-displacement direct shear test, they found that the shear strength as well as the crack pattern is influenced by both joint spacing and orientation.As stated for the shear failure of the intermittent joint, the research focuses mainly on the cracking types and the relationship with joint orientation.Even though the visual cracks were identified, the internal stress distribution surrounding a joint was not well exhibited under shearing.erefore, a more comprehensive understanding of the intermittent joint under shearing is necessary.But it is quite difficult to derive the internal stress for a physical test.However, a numerical approach covers the shortage.
Numerical approaches are often used to further explain the physical failure phenomenon.To date, continuum methods, discontinuum methods, and hybrid continuum/discontinuum methods have mostly been used to simulate the fracture behavior of rock masses.e differences between these methods have been discussed by researchers [4,25,26].Continuum methods cannot present explicitly the crack initiation, propagation, and nucleation because mesh grids in the problem domain should remain in the same neighborhood during the whole deformation or transport process and thus block rotation so that opening of the mesh will not happen.In contrast, discontinuum methods are able to model zone detachment, large displacements, and block rotation, and therefore, the fracturing process in rock or rock-like materials can be presented explicitly by microcracks formed in models based on the discrete element method (DEM) [5,[27][28][29][30][31][32].Additionally, based on the rock failure analysis process, the numerical direct shear test was performed to investigate the shear behavior of intermittent joints with different geometrical parameters.Meanwhile, the whole shear failure process was visually recorded, which is the greatest advantage relative to the physical test [33].Among the developed numerical approaches, particle flow code (PFC) is useful.Using the particle flow code approach, the internal stress evolution could be well exhibited, especially the internal stress evolution surrounding an open joint.Moreover, the magnitude of the internal stress can also be derived and can be used to account for the cracking mechanism.According to the published reports regarding the shear failure of intermittent joints, particle flow code has not yet been used to compute the principal stresses around an open joint in shear.Additionally, macrocrack propagation is composed of microcracks.However, the microcracks are impossible to identify in a physical test.
As stated above, the shear failure behavior of the intermittent rock joint has not been investigated currently relative to the compressive failure behavior.To further explore the shear failure mechanism of intermittent rock joints, physical and numerical tests both were conducted in the study.Two series of intermittent joint were molded using gypsum material by changing the joint arrangement along the shear axis.With the aid of the bonded particle model, the increase in microcracks was related to macrocrack propagation, and the internal stress distribution around the joint was computed to explain the shear failure mechanism.

Physical Test and Numerical Simulation
2.1.Joint Arrangement.In the present work, two series of joint configurations were considered.e dimensions of the prismatic specimen were 150 mm × 150 mm × 30 mm (height × length × width).ree joints were arranged within each sample along the shear plane.To better describe samples hereafter, the two series were named A and B. In Series A illustrated in Figure 1(a), two joints were located at the two lateral edges and kept the inclination angle of 0 °fixed while the middle one was located at the sample centroid.e middle joint inclination angle, β A , rotating in the clockwise direction, was set from 0 °to 150 °with an interval of 30 °.In Series B illustrated in Figure 1(b), all three joints rotated along the joint center in the counterclockwise direction.e inclination angle, β B , was also set from 0 °to 150 °with an interval of 30 °.Each joint length (L) was set to 20 mm. e shear directions for Series A and B are shown in Figures 1(a) and 1(b), respectively.e ligament connecting the two tips of each neighboring joint is called the rock bridge.e detailed joint configuration can be found in Table 1.

Sample Preparation.
Gypsum is a brittle material like rock material, so it is usually used as a rock-like material to model the behavior of rock specimens [9,10,23,24].First, gypsum powder and water were mixed at a ratio of 1 : 0.35 in weight and then stirred well.Next, the liquid gypsum was quickly cast into steel molds and vibrated.e liquid gypsum often consolidated in less than twenty minutes.For jointed specimens, the mica sheets should be inserted according to the designed location and inclination angle.Since gypsum is the air setting material, the samples do not need a water spray cure.

Physical Test.
A direct shear test was performed to investigate the failure characteristics of the jointed rock.e test was conducted on the shear-rheology instrument RYL-600 shown in Figure 2. is instrument consists mainly of the control system, shear loading system, normal loading system, and input/output devices including a computer and printer.
e jointed sample was placed in the shear box first, and then the normal stress was applied at 2.5 MPa with a speed of 50 N/s.en, the shear loading was applied at a speed of 2 mm/min.
2 Advances in Civil Engineering 2.4.Numerical Simulation.Among the numerical approaches mentioned in Section 1, DEM simulations o er another merit of allowing microscale behavior such as contact forces, particle displacement vectors, and internal stresses to be obtained.Considering these inherent advantages of simulating the fracturing process in rock or rock-like materials, Particle Flow Code (PFC) was used to explore the rupture mechanism of rock-like specimens with joints in shearing.e parallel-bonded particle model (PBPM) was employed in this study, which has been extensively used and validated in simulating the fracturing process of rock or rock-like materials in a ood of literature reviews [5,16,22,28,29].Since the microproperties of the real physical material are not determined directly in the laboratory, calibrating the PFC synthetic specimen requires an iterative trial-and-error procedure.e microparameters for the PFC model were determined by conducting uniaxial compression and Brazilian numerical simulations on synthetic intact specimens without any preexisting aws or joints.As the speci c calibration procedure is rather mature, it will not be introduced in this paper.Please refer to the PFC 2D manual for more details concerning microparameter calibration [34].e model porosity is not equivalent to the porosity of real intact rock attributed to the utilization of circle disks to represent the real grains [35].Moreover, the mode of tensile or shear failure of a contact is dependent on the applied local stress exceeding either the tensile or shear strength of the contact.After an iterative procedure, all the microparameters were determined as listed in Table 2. e PFC synthetic intact specimen is composed of 43,298 balls and      Advances in Civil Engineering 97,872 parallel-bonded bonds.e macroproperties of the PFC intact specimen are comparable to those tested by the physical experiment as listed in Table 3. e numerical direct shear setup (Figure 3(a)) is a simplification of the experimental direct shear setup in the laboratory.e gap between the upper and lower parts of the shear box is 2 mm.All of the contacts between the wall and particle are frictionless.First, an initial velocity was exerted on the top wall of the shear box (Wall 2) until the normal stress reached the desired value.e normal stress was obtained through averaging the reaction forces acting along both the top and bottom walls of the shear box (Wall 2 and Wall 1) divided by the length of the synthetic specimen.
en, the shear displacement was applied to the upper part of the shear box (Wall 3 and Wall 4) by a fixed shearing velocity when the normal stress remained constant.To ensure that the whole shearing process was in a quasistatic loading condition, a shearing velocity of 0.006 m/s was adopted.e loading rate in the simulation is different from the loading rate in the physical experiment.For instance, the loading rate of 0.006 m/s in this situation means approximately 1.5 × 10 −9 m/step because the timestep is 2.47 × 10 −7 s.To be consistent with the experiment, the lower part of the shear box (Walls 1, 5, and 6) was kept fixed throughout the shearing.e average shear stress was calculated by the reaction forces exerted along Wall 6 divided by the specimen length.Contact forces denoted by force chains exist between particles.Compressive and tensile contact forces are marked black and green, respectively, as shown in Figures 3(b

Category of Crack Patterns.
From the published reports [1,20,23,28,30], we know that in the condition of compressive loading, the crack often initiates at the joint tips, and the crack initiation type is related to the joint orientation.e loading modes of shear and compression are quite different.Under shear loading, the shear plane is subjected mainly to the shear stress.For coplanar and non-coplanar joints on the shear plane, internal stresses around joints resulting from the shear loading are different, so that cracking behaviors are different.In the study, a main crack (MC) and secondary crack (SC) were observed.e main crack is the category of crack connecting the tips of two neighboring joints.e SC is the kind of crack initiating at the joint tips before or after the peak shear stress (τ p ). e SC seldom connects neighboring joints except in B120 and B150 cases.Both Samples A0 and B0 can be regarded as the coplanar intermittent joints.e joint orientation is consistent with the shear direction as well as the shear plane.erefore, cracks are first initiated at the tips and then easily propagated along the shear plane.
is MC is named MC-I, as shown in Table 4.In this case, since the direct shear failure leads to the rock bridge rupture, the generated cracks are shear cracks, and no SC was observed.From 30 °to 150 °, one (Series A) or three joints (Series B) are inclined to the shear direction, and the rock bridge is longer relative to the rock bridge of 0 °samples.e cracking is quite different from the cracking in 0 °cases.In Series A, since only the middle joint orientation varied, cracks initiated from both tips of the middle joint and one tip of the edge joints, and then propagated oppositely and finally coalesced.In the case of the joint arrangement, as shown in Figure 5, cracks always propagated through A (solid line), not B (dashed line), that is, the cracking trajectory between two neighboring joints followed the shortest way. is MC type is named MC-II.In Series B from 30 °to 150 °, MCs always also initiated at the joint tips, and cracks obliquely propagated towards another tip of adjacent joints.For another two joints close to the sample edge, cracks initiated at the joint tip approaching the sample  4).If the shear loading continued after MCs completely coalesced, in some cases of Series A, an SC may be generated at the lower middle joint tip (SC-I) such as that shown in Figure 6(a).For B30 and B60, SCs initiated at the upper and/or lower joint tips and propagated in the direction approximately perpendicular to the joint orientation (Figure 6(b)).ey are also regarded as SC-I.SC-I generation is due to mutual squeezing between upper and lower blocks.For B90, B120, and B150, SCs (SC-II) also initiated at the upper and/or lower joint tips, but the propagating direction of SC-II is approximately consistent with the shear direction that is quite di erent from the shear direction of SC-I.SC-II generation is due to the tensile internal stress acting at the joint tips before the peak shear stress, which will be analyzed in Section 3.6.Once the MCs have completely coalesced, regular interfaces with one sawtooth in Series A or several in Series B were created.e inclination angle of the sawtooth varies with the joint orientation.

Advances in Civil Engineering
A detailed description of the failure of the sawtooth interface will be given in the following context.e MC cannot be simply regarded as a pure tensile or shear crack.Actually, the MC is a mixed mode, as shown in Figure 6(a).
e shear dust and the tensile crack could be observed on one MC trajectory.e MC mode is related to the joint orientation.Certainly, two types of SCs are both generated at the sawtooth roots of the upper or lower block under the tensile stress, so all SCs are regarded as tensile.e mechanisms of SC-I and SC-II are different, which will be numerically explained in Section 3.4.

Curve Types of the Shear Stress-Shear Displacement.
e shear stress-shear displacement curve (shearing curve) reveals the variation tendency of shear resistance.In addition, the cracking process can be indicated partly by the shearing curve.For instance, the SC propagation may create a second peak shear stress.It will be explained in detail in the following context.Among all curves of shear stress versus shear displacement, two typical curve types are identified.As shown in Figure 7(a), the shear stress monotonically increases, with the shear displacement increasing before the peak value.After one certain shear displacement, the shear stress reaches the peak value.Actually, the shear displacement at the peak shear stress is dependent on the joint configuration.After the peak value, the shear stress drops rapidly and then enters into a long residual stage.Due to the continuous friction and sliding behavior on the rupture interface, within the residual stage, the shear curve is not smooth but slightly fluctuating.is described curve type has one single peak only, so it is named a single-peak curve (Type I).
e single-peak curve was mostly found in the shearing process of Samples A0, A90, A120, A150, and B0.
e other curve type (Type II) is shown in Figure 7(b).Like the single-peak curve prior to the peak, the shear stress also  6 Advances in Civil Engineering purely increased with the shear displacement increasing.After the rst peak, the shear stress dropped approximately vertically to some degree.With new cracks generated or the extrusion on the sawtooth interface after MC coalescence, the shear stress rose up again after experiencing a short shear stress valley.Due to the second increase in the shear resistance, a second shear stress peak was generated.e second peak shear stress is as a rule lower than the rst peak shear stress.After the second peak, the shear stress gradually comes into the uctuated residual stage.According to the shear curve features of Type II, it is named a double-peak curve.e type II curve can be found during the shearing process of Samples A30, A60, B30, B60, B90, B120, and B150. is curve type has two peaks, which is the biggest di erence between the two curve types.e di erence between the two curve types denotes the fracturing di erence.According to the crack patterns described above, after the MC coalescence between neighboring joints, the SC propagation or extrusion between the upper and lower blocks on the sawtooth face nally results in the second peak.e curve shape indirectly reveals the fracturing process during the shear test, which will be described in detail in the following section.A shear displacement that is long enough can bring about a complete shear curve.

Crack Initiation Stress Ratio.
Fracturing may be revealed by the shear stress-shear displacement curve to some extent, as stated above.In the numerical simulation, crack initiation and coalescence can easily be recorded.According to published reports, the crack initiation stress is associated with the joint con guration.In particle ow code, the crack initiation stress is de ned as the stress when the number of microcracks reaches 1% of the peak stress, so the crack initiation stress is associated with the fracturing process.e crack initiation stress ratio is de ned as the ratio of the crack initiation stress to the peak shear stress.e real crack initiation stress is very di cult to be accurately obtained.Instead, it is easy to be recorded using the numerical approach.Figure 8 shows that the crack initiation stress ratio varies with the inclination angle.e crack initiation stress ratio for 30 °samples in both series is the lowest ratio, indicating that the most further shear stress is needed from the rst crack appearance to sample failure when the joint inclination equals 30 °. e trend of the crack initiation stress ratio for Series A and Series B is similar from 0 °to 90 °, decreasing rst and then increasing.e trend is distinguished from the trend obtained in compression tests, which increase from 0 °to 75 ° [34].
us, it demonstrates further that the fracturing behavior of jointed rocks is related to the loading condition.However, the variation trend is di erent between Series A and Series B from 90 °to 150 °. e total length of sheared rock bridges from the crack initiation stress to the peak stress may be contributed by this di erence.For instance, in Series A, the MC propagation and coalescense are along the shortest rock bridge (rotating compared with that of samples with lower joint inclinations), and the total length of the sheared rock bridges decreases with the joint inclination increasing from 90 °to 150 °.In contrast, the total length of the sheared rock bridges for Series B increases with the joint inclination increasing from 90 °to 150 °, as shown in Table 4.
erefore, the crack initiation stress ratio for Series A rises while the crack initiation stress ratio for Series B is reduced from 90 °to 150 °.

Phase Division of the Cracking Process.
On the basis of the function of cracking recording in PFC, the cracking process is described based on numerical results using four typical fracturing behaviors.To clearly state the shear failure process, according to the fracturing features at di erent stages, four typical fracturing moments, that is, visual crack initiation (Point A), one MC coalescence (Point B), complete crack coalescence (Point C), and the residual stage (Point D), are presented to stand for the whole fracturing process.As mentioned above, particles in the numerical model are cemented by the parallel bond, which may break in the shear or normal direction when the local stress acting on the parallel bond exceeds the predetermined shear or tensile strength of the bond.Each broken bond creates one e number of microcracks can be recorded in the shearing process.
According to the definition of crack initiation stress, cracks are not obviously distinguishable at the initiation stress.e cracks are assemblies of a few microcracks, and the cracking trend is hard to identify.erefore, to more distinctly evaluate the initial crack, visual cracks are presented and the corresponding shear stress is marked on the shearing curve.A combination of Figures 9 and 10 shows the fracturing process of A30 at different shear stresses.When the shear stress was approaching 0.9 τ p , it was observed that cracks first initiated at the tips of the middle joint and the other two edge joints (Point A, Figure 10 Under the shear load, the sawtooth may be sheared off under a high normal stress or the upper block climbs up along the contact face under a low normal stress.In the present work, the strength of the modeling material, gypsum, is quite low.For A30, under a high normal stress, despite climbing is limited, when cracks coalesced completely, the upper block climbed a little, so the normal displacement increased as shown in Figure 9.During climbing, the sawtooth face is extruded, which creates a second peak shear stress.Due to extrusion, the SC (Figure 10(f)) may be generated at the lower tip of a middle joint like A60 shown in Figure 6(a).Before Point A, the number of microcracks increased a little, but from Point A to Point C, the number of microcracks increased rapidly.After Point C, the number of microcracks increased slowly.Failure of A60 is analogous to the failure of A30.
As shown in Figures 11 and 12, cracks first initiated at approximately 0.8 τ p (Point A, Figure 12(a)).
e crack initiation locations are identical to the crack initiation locations of A30. e MC first coalesced at the left rock bridge from Point A to Point B (Figure 12(b)), within which the shear stress reached the peak value.Cracks completely coalesced from Point B to Point C, that is, the rapid shear stress declining stage (Figure 12(c)).Comparing Figures 12(c) and 12(d), the numerical crack coalescence is in good agreement with the experimental crack coalescence.As Figures 12(c) and 12(d) show, the crack propagation still followed the shortest trajectory.
e MC coalescence separated the sample into upper and lower blocks and created an interface with one sawtooth (Figure 12(e)).Since the upper block moved away from the middle joint, the joint gap got wider and wider with the shear displacement increasing.Compared with Figure 10, since no climbing took place, the normal displacement did not keep increasing like Sample A30. e peak normal displacement may be resulted from that the upper block slide of the irregular interface.Namely, the interface between the upper and lower blocks is composed of particles, the interface is just a little fluctuate because of the arrangement of unequal radii particles.erefore, no second peak and extrusion on the sawtooth face were observed.e variation trend of Sample A150 is like the variation trend of A30.Before the visual crack initiation, the crack number increased slowly.Cracks sharply increased from initiation to crack coalescence.During the residual stage, the microcrack number increased slightly.Failure of A120 and A90 is analogous to the failure of A150.
A combination of Figures 13 and 14 shows the failure process of Sample B0.When the shear stress is close to the peak value (Point A, Figure 13), cracks initiated at the joint tips (Figure 14(a)).Once cracks initiated, cracks grew sharply, and completely coalesced along the shear plane (Figure 14(b)).Certainly, the peak shear stress was obtained from Point A to Point B, namely, during the crack coalescence.Comparing Figures 14(b) and 14(c), a good agreement between numerical and experimental failures can be found.For 0 °cases, the shear plane was directly subjected to the shear stress, so the failure occurred along the shear direction on the shear plane.For the cases of A0 and B0, which are coplanar joints, the joint orientations are consistent with the shear plane.Once the cracks coalesced, A0 and B0 can be regarded as smooth joints.In this case, the upper block just slides along the ruptured interface.e normal displacement is also quite small.e reason is the same like the above-stated A150.If the shear load continued to apply, the friction on the shear plane contributed to a stable residual strength (Figure 14(d)).Since no SC was generated, the single-peak shear stress was observed.ough the joint arrangement in Series B is slightly different from the joint arrangement in Series A, the failure of A0 is like the failure of B0. e variation trend of Sample B0 is also like the variation trend above, but microcracks are basically observed on the shear plane.
For Sample B90, the joints are perpendicular to the shear plane as well as the shear direction.e shear plane traverses the midpoints of three joints.At a shear stress of approximately 0.9 τ p (Figure 15), visual tensile cracks, namely, SCs, first initiated at most joint tips (Figure 16(a)).For one identical joint, two tensile cracks initiating at the upper and lower tips propagated towards opposite directions.As shown in Figure 16(b), even though MCs have not completely coalesced yet, the shear stress has reached the peak and then starts decreasing.During the shear stress decreasing stage, MCs completely coalesced obliquely (Figure 16(c)).Meanwhile, the shear stress was still dropping.Comparing Figures 16(c) and 16(d), the main and SCs are basically identical for numerical and experimental failures.Complete MCs separated the original sample into two blocks with a sawtooth interface (Figure 16(e)).e sawteeth mutually interlocked and prevented the upper block from climbing along the interface.With the continuous shear load, some sawteeth could be sheared off along the initial tensile cracks (Figure 16(f)), so that a second shear stress Advances in Civil Engineering peak was created (Figure 15).e normal displacement trend with shear displacement for B90 is different from B0.
e normal displacement increased with increasing shear displacement.It is because the upper block climbed slightly along the sawteeth interface.Due to cracking along the first initial SC, the second peak is lower than the first peak.e variation trend of the microcrack number is different from the aforementioned variation trends.e variation trend will be described in Section 3.7 in detail.
Figures 15 and 16 represent typical ruptures.en, samples withJ different joint arrangements are known to experience different failure processes.e failure process is closely related to the joint arrangement.According to the whole fracturing process of all cases, three phases can be summarized as in the following sections.

First Phase: Cracking.
e cracking phase includes the crack initiation, propagation, and coalescence as well as the SC initiation.All cases in Series A and Series B have to experience the cracking phase.Certainly, cracking types are distinguished.In Series A, the MC coalescence types are the same, but the SC initiation is influenced by the combination of shear direction and joint orientation.at is why the SC 10 Advances in Civil Engineering initiation can be observed for A30 and A60 but not for A90, A120, and A150.In Series B, MCs coalesced obliquely between two neighboring joints.SCs were generated in all cases except B0.However, the cracking mechanism of SCs for B30 and B60 is different from the cracking mechanism of B90, B120, and B150.Cracking is also influenced by the combination of shear direction and joint orientation.

Second Phase: Extrusion and Sawtooth Cutting.
According to the abovementioned fracturing process, the second failure phase is obviously not necessary for all samples.As stated above, parts of samples such as A30 and A60, B30 and B60 experienced the extrusion at joint faces as shown in Figure 17(a) and Figure 17(b).A sawtooth interface was created due to MC coalescence.Sawteeth on the upper and lower blocks interlocked, under continuous shear loading.e extrusion effect took place on the contact face between the sawteeth.For B90, B120, and B150, after MC coalescences, some sawteeth were cut off along the first initial tensile cracks like that shown in Figure 17(c).

3.4.3.
ird Phase: Friction and Sliding.When the shear stress experienced a first and second peak, it came into the residual stage.On the interface, the friction and sliding effect played the main role.For A0 and B0, once cracks completely coalesced, blocks oppositely slid on the interface, so they did not experience the second phase and directly came into the third phase.For some cases, they came into the third phase after sawtooth cutting or extrusion.e sliding and friction effect on the interface contributes to the residual shear stress.For the smooth interface such as A0 and B0, the residual shear stress is low relative to that for the case in which the interface is composed of small fragments.Fragments result in undulate residual shear stress.Figure 19 shows the contact force evolution of Sample B90.At the beginning of loading, within the whole sample, compressive contact forces have not yet distinctly presented at joint tips, but tensile contact forces started intensively locating along joints (Figures 19(a1) and 19(a2)).
en, compressive contact forces concentrated at the joint tips.Meanwhile, tensile contact forces also presented at the joint tips.Cracks initiated at the joint tips, and tensile contact force concentration can account for the SC initiation (Figure 16(b)), which will be analyzed further using the principal stress around the middle joint in Section 3.6.en, compressive contact forces mainly concentrated at rock bridges (Figures 19(c1  e location of the contact force concentration may be where cracks initiated, that is, where the local stress is high.
e contact force concentration qualitatively represents the location of high internal stress and also indicates where cracks may initiate and propagate.Even though fracturing processes are various for samples with different joint orientations, the contact force evolution has to experience the following phases.e compressive contact force evolution should experience at least three stages, including the first uniform distribution, the concentration at joint tips and rock bridges, and the slight concentration on the failure plane.e tensile contact force evolution is partly different from the compressive contact force evolution, experiencing mainly a scattered force distribution close to the sample edges, concentration on joint sides and rock bridges, and concentration at the tensile crack tips.erefore, the compressive contact force evolution is good to explain the MC propagation and coalescence, while the tensile contact force evolution is for the SC initiation and propagation.

Analysis of Principal Stresses around the Middle Joint.
As shown in Figures 10,12 and 16, cracks often initiated at the joint tips.To explain crack initiation at the joint tips and understand the stress evolution around joints, particularly the middle joint, the internal stress evolution of A30 and B90 was investigated as examples.Owing to the discreteness of the PFC model, stress, which is a continuum variable, actually does not exist anywhere in the particle assembly.To overcome this problem, an averaging procedure is adopted, which also bridges the gap between the microscale and a continuum.erefore, the stress at one point is the average value in a measurement region of volume V, computed as the following formula: where N c is the number of contacts within the measurement region, F (c) and L (c) are the contact vector and the branch vector joining the centroids of the two contact ends, respectively, and ⊗ represents the outer product.F (c) includes both normal and shear force components at contact but neglects the parallel-bond moment [34].
On the periphery of the middle joint, 16 measurement circles are arranged.e schematic diagrams of the measurement circle arrangement are shown in Figure 20.e diameter of each measurement circle is 3 mm.Each circle contains approximately 14 particles.Within one measurement circle, σ xx , σ yy , and τ xy can be obtained, according to the following equation: en, the maximum principal stress (σ max ) and minimum principal stress (σ min ) for each measurement circle are e positive and negative signs represent the tensile and compression stresses, respectively.e principal stress distribution around the middle joint is illustrated using two typical samples, that is, A30 and B90.As seen in Figures 21(a) and 21(b), σ max and σ min are found to be varied at di erent locations in A30, and the highest stresses take place at circles 1 and 9, which are located at the joint tips.At circles 1 and 9, both the maximum and minimum principal stresses are always positive ones, which means that they are under the biaxial compression and tend to be shear failure.In other circles, the magnitude of computed principal stresses is much lower than the magnitude of computed principal stresses at circles 1 and 9. Furthermore, the magnitude of minimum principal stresses at circles 8 and 16 is the highest, exceeding the tensile strength of the synthetic sample.Combining the maximum principal stresses at these two locations (Figure 21(a)), circle 8 is under axial tension, and the stress state at circle 16 resembles that in the Brazilian tests.
ese well explain the appearance of tensile microcracks at circles 8 and 16 (Figure 10(a)).Circles 1 and 16 and circles 8 and 9 are located at the left-top and right-bottom of the middle joint, respectively.Because the shearing direction for the upper part is right to left, the location of circles 1 and 16 relative to the joint is along the shearing direction.is rule is also suitable for the location of circles 8 and 9. Taken together, all of these factors drive the crack initiation at the joint tips and guide the crack propagation direction.
Figure 22 shows the maximum and minimum principal stresses at di erent circle locations in B90 for the di erent shearing stages.Like sample A30, at circles 1 and 16 and 8 and 9 in sample B90, both the maximum and minimum principal stresses remain positive and their magnitude is much larger than the magnitude at other measurement circles.ey are under biaxial compression and inclined to fail by shearing.Circles 1 and 16 and 8 and 9 are located at the right-top and left-bottom of the middle joint, respectively, which, relative to the joint, are along the shearing direction.It can be used to explain cracks formed at the joint tips and guide the crack propagation direction.In addition, the stress states at circles 2, 10, and 11 (Figure 22(a)) are analogous to the Mode I fracture, and the magnitude of minimum principal stresses at these locations exceeds the tensile strength of the synthetic sample.
us, these two factors result in the formation of tensile microcracks at these locations (Figure 16(a)).Namely, it well accounts for the formation of SC-II cracks.
Overall, the stress concentration at the joint tips is associated with the crack initiation, and the location relative to the joint is along the shearing direction, which guides the crack propagation direction.Moreover, the formation of SC-II cracks is attributed to the tensile stress.

Variation Trend of the Number of Microcracks.
In the numerical modeling process, the increase in the microcrack number is directly related to the cracking behavior.Cracking behaviors are various for di erent joint-arranged samples.
e visual crack during the shear failure is an assembly of microcracks.Since the macrocrack propagation is associated with the shear displacement as shown in Figures 10,12,14,and 16, the increase in the microcrack number is also associated with the shear displacement as well as the fracturing behavior.Figure 23 shows the varying relationship between microcrack number and shear displacement.At di erent stages of shear displacement, growth rates for the microcrack number are di erent, which can be evaluated by the slope of the microcrack number versus the shear displacement curve.A big slope means a rapid growth rate.At the beginning of shearing, since the internal stress has not reached one value that can break parallel bonds, when the shear displacement increased from 0 to 1 mm, even though the shear stress dramatically increased with the shear displacement, the microcrack number obviously increased a little.is increase can be named a slow increase phase.
is phase corresponds to the crack initiation.When the shear displacement ranged within 1 to 2 mm marked with two vertical red dashed lines, the microcrack number sharply increased approximately vertically, which corresponds to the MC propagation and coalescence like that shown in Figures 12  and 14. e approximately vertical phase is called the rst rapid increase.For Series A, the rst rapid increase in Advances in Civil Engineering microcrack number was within 1 to 2 mm, but for Series B, the rst rapid increase in microcrack number nished in less than 2 mm, such as the dashed curves of B30 and B60. e rst rapid increase mainly corresponds to MC propagation and coalescence.In some cases of Series B such as the curves for B30 and B60, due to the SC propagation or extrusion, the second rapid increase could be observed but hardly took place in Series A. As the shear displacement continued after MC coalescence, due to the friction e ect on the rupture interface, the microcrack number still continued to increase, but its growth rate became slow relative to the growth rate of the rst rapid increase phase, which can be evaluated upon the curve trend.Generally, four phases of microcrack number increases, that is, slow increase, rst rapid increase, second rapid increase, and continuous increase, are summarized.e second rapid increase can be observed in a few cases in Series B.

Conclusions
In the present work, two series of intermittent joints with di erent joint orientations were used to carry out direct shear tests.Experimental and numerical approaches were used to explain the macro-and microfracture mechanisms.Some interesting conclusions are drawn as follows: (1) Main and SCs are identi ed according to mechanism and crack pattern.ey both initiate at the joint tips.MCs connect the tips of two adjacent joints.In Series A, MCs always propagate following the shortest trajectory.In Series B, they obliquely go through rock bridges.SCs may initiate prior to or after the peak shear stress and propagate in the direction that is approximately perpendicular to the joint orientation or parallel to the shear direction.SCs are all tensile cracks, while the MC is usually a mixed mode of tensile and shear cracks.(2) e fracturing process is signi cantly related to the joint orientation.e whole fracturing process is identi ed as three phases, which are characterized by cracking, extrusion or sawtooth cutting, and friction and sliding, respectively.But three phases are not undergone for all cases.e rst and third phases are necessary for all cases, but the second phase is necessary only for A30, A60, B30, B60, B90, B120, and B150.(3) Shearing curves are characterized by the number of peaks, so two types, single-peak and double-peak curves, are identi ed.A single-peak curve is derived from A0, A90, A120, A150, and B0, while the doublepeak curve is derived from A30, A60, B30, B60, B90, B120, and B150.(4) e variation trend of the crack initiation stress ratio is associated with the loading condition and the joint inclination.From the joint inclination of 0 °to 90 °, the crack initiation stress ratio under the shearing decreases rst reaching the valley at 30 °and then increases, which is di erent from the crack initiation stress ratio under the compression loading.From the joint inclination of 90 °to 150 °, the crack initiation stress ratio keeps increasing for Series A while decreasing for Series B. In essence, this di erence originates from the variation of crack propagation and coalescense paths.Advances in Civil Engineering (5) e compressive contact force evolution may undergo initial uniform distribution, concentration at the joint tips, concentration at rock bridges, and an extrusion surface.e tensile contact force evolution may experience initial scattered distribution close to the sample edges, concentration on both joint sides, with concentration at the rock bridges and SC tips.e compressive contact force evolution is good to explain MC propagation and coalescence, while the tensile contact force evolution is good for SC initiation and propagation.(6) Principal stresses around joints are associated with crack initiation and propagation.e highest maximum and minimum principal stresses take place at the joint tips and guide crack propagation, accounting well for crack initiation at the joint tips and propagation trajectory.(7) At the beginning of shearing, the number of microcracks increases slowly within a certain shear displacement mainly ranging from 0 to 1 mm.e increase in microcracks can be classi ed into four phases, namely, slow increase, rst rapid increase, second rapid increase, and continuous increase.e four increase phases correspond to crack initiation, MC propagation, extrusion or SC propagation, and sliding on failure interface.e second rapid increase is observed only in the shearing of several samples that may undergo the second phase of shearing.
As analyzed above, the failure of discontinuous joints under a shear load mainly takes place on the shear plane.e failure characteristic is dependent signi cantly on the joint orientation and joint con guration.e presented results can be used in other rock mechanics applications, for instance, the sliding analysis of rock slope with discontinuous joints and the stability of an underground cavern with jointed rock mass.

Figure 1 :
Figure 1: Schematic diagram of the joint con guration for (a) Series A and (b) Series B.

2 . 0
Particle friction coe cient, μ 0.5 Parallel-bond modulus, E c (GPa) 0.6 Ratio of parallel-bond normal to shear sti ness, k n /k s 2.0 Parallel-bond normal strength, mean σ m (MPa) 9.0 Parallel-bond normal strength, standard deviation σ sd (MPa) 2.0 Parallel-bond shear strength, mean τ m (MPa) 5.4 Parallel-bond shear strength, standard deviation τ sd (MPa) 1.0 ) and 3(c).Before loading, the compressive contact force exists basically uniformly inside samples, but the tensile contact force exists mainly near the sample edges.For the jointed specimens, the open joints were created by deleting particles located in the corresponding regions shown in Figure4(a).e arrangement of joints is the same in the physical flawed specimens.Even though open joints were created, compressive and tensile contact force distributions (Figures 4(b) and 4(c)) are still identical to the compressive and tensile contact force distributions for the intact model, respectively.No contact force concentration is observed around joints.

Figure 3 :
Figure 3: Intact numerical model: (a) schematic diagram of direct shear setup; (b) compressive contact force chain; (c) tensile contact force chain.

Figure 5 :
Figure 5: Schematic diagram of the main cracking trajectory in Series A.

Figure 8 :Figure 9 :
Figure 8: Crack initiation stress ratio of Series A and Series B samples.

Figure 17 :
Figure 17: Schematic diagram of extrusion and sawtooth cutting.

Figure 18 :
Figure 18: Contact force evolution for Sample A30 (shear direction: right to left): (a) compressive contact force; (b) tensile contact force.

Figure 19 :
Figure 19: Contact force evolution of B90 (shear direction: left to right): (a) compressive contact force; (b) tensile contact force.

Figure 20 :
Figure 20: Schematic diagram of the measurement circle arrangement.

Figure 21 :
Figure 21: Principal stresses around middle joint for A30 at (a) crack initiation stress and (b) peak shear stress.

Table 1 :
Joint inclination angle and sample number.

Table 2 :
Microparameters in bonded particle model.Ratio of particle normal to shear sti ness, k n /k s

Table 3 :
Comparison of experimental and numerical results.

Table 4 :
Category of crack patterns.