Crack Initiation , Propagation , and Failure Characteristics of Jointed Rock or Rock-Like Specimens : A Review

School of Resources and Safety Engineering, Central South University, Changsha, Hunan 410083, China School of Civil, Environmental and Mining Engineering, #e University of Western Australia, Perth 6009, Australia Engineering laboratory of Spatial Information Technology of Highway Geological Disaster Early Warning in Hunan Province, Changsha University of Science & Technology, Changsha 410114, Hunan, China School of Highway, Chang’an University, Xi’an 710064, China School of Resources and Environment Engineering, Wuhan University of Technology, Wuhan 430070, China


Introduction
After hundred millions years of geological movements, there is a large number of discontinuities in natural rock mass.e natural rock mass is clearly a kind of heterogeneity material [1], and the discontinuities have a great influence on the mechanical behavior of rock mass.In most engineering cases, such as roadway excavation, rock slope, and pillars in deep mining activities, there often requires the estimation of the strength and failure characteristics of a rock mass that contains a large number of discontinuities [2][3][4][5].And, many engineering disasters are closely associated with crack initiation and propagation between preexisting joints (Figure 1).For example, the collapse of the Malpasset arch dam (French) results from the crack propagation in rocks of bam foundation.And, there are also many excavation activities that have to face with the crack propagation and coalescence in rock mass.en, it is very necessary to investigate the strength and failure characteristics of jointed rock mass.At the same time, the crack initiation, propagation, and failure characteristics are of great interest to engineers and scientists.
is review concentrates on experimental techniques for mechanical behavior of jointed natural rocks or brittle rocklike materials.All literature available to the authors (in total, 196 references) concerning this topic was extensively reviewed.
is review is structured in the following way.After "Introduction," the synthesis, analysis, and evaluation are discussed in Section 2. In Section 3, the loading techniques are discussed, and for each loading pattern, the mechanical characteristics of the specimen with different joint geometry parameters have also been discussed.Section 4 presents widely used numerical simulation techniques for crack initiation, propagation, and failure process of rock or rock-like materials.Finally, a brief summary and some prospective research are presented.

Experimental Materials and Specimen Preparation
In recent years, the failure behavior of jointed rock mass has been widely studied by scholars, and the experimental test is the most common way.So far, many types of materials have been employed for such tests.As we all know, the ideal material for experiment is the natural rocks.And, there are many kinds of rocks which have been used in previous studies, such as granite, sand stone, and marble (Table 1).As shown in Figure 2, the specimens are made from natural rocks.e fissures in the specimen are open, and the fissure is produced via high-pressure water-jet cutting method.
Apart from the natural rocks, many kinds of rock-like materials have also been widely used by scholars, such as glass [7], Columbia Resin 39 [8], and molded gypsum [9][10][11][12][13][14][15][16]; the most common among them is gypsum and cement mortar (Table 1).e cement mortar is used because the main framework of the specimen consists of sand and cement, and the sand can provide frictional behavior of the modelling material.
en, the cement mortar is suitable for modelling rock mass.Compared to the cement mortar, the frictional behavior on the failure plane is inapparent.
For most rock-like material specimens, the fissure or joint was mainly handcrafted.In Figures 3(a)-3(c), the open fissure in PMMA, gypsum, and cement mortar specimen are shown, respectively.In the specimen preparation, preexisting fissures are created by inserting metal sheets into the fresh mixture at the desired locations of the fissures and removing them after rock-like materials become hardened.
At the same time, in the previous studies, the closed fissure or joint can also be created by scholars through different ways.In Figure 4(a), the ubiquitous joints created by inserting mica sheets into fresh cement mortar are shown [20].e mica sheet has also been used by other scholars to investigate the failure characteristics of jointed rock mass [21].Because of the strength of the mica sheet and relatively low interface between mica sheet and mortar, the mica sheet is suitable for modelling the joints which can be filled up with clay.In Figure 4(b), the closed fissure created by inserting steel sheets into the mortar is shown [22].During preparation process, the steel sheets were removed carefully after about 3 ± 5 min of curing, and then, further expansion of the plaster mixture during curing closed the slots, thereby creating preexisting cracks with closed surfaces.And, to obtain a crack surface with different roughness, three kinds of roughness of the stainless steel sheets were used.Like the fissure shown in Figure 4(c), the plat closed fissure is created by inserting a video tape strip [23].e video tape strips are pulled out of the mold 30 min after vibration, and the gypsum is soft enough and can expand to fill the gap left behind by the strips.Moreover, instead of pull out, some scholars also attempt to create closed joint via the left behind the galvanized sheets in the mixture [24].Actually, the stiff of galvanized sheet is far below the rock-like material.Because the interface is relatively smooth, the slide effect on the interface is also obvious; it is similar to the mica sheetmortar interface.

Failure Characteristics of Jointed Rocks or
Rock-Like Materials under Different Loading Methods

Uniaxial Compressive
Loading. e uniaxial compressive loading is the most common load mode in previous works.And, in many engineering cases, the rock mass failure occurs from the crack coalescence under axial loading.In Figure 5, the pillar failure from the coalescence between two inclined fissures is shown.It is clear that the preexisting fissure has a great influence on the strength, deformation, and failure characteristics of rock mass.In order to better understand the failure process of engineering rock masses, experimental studies have focused on the jointed specimens that contain preexisting joints. 2

Advances in Civil Engineering
For the specimen with single joint or fissure, under uniaxial loading, two types of cracks will initiate from the preexisting fissure, namely, wing and secondary cracks [7,28,40,92,101,102].With further loading, new cracks propagate and reach the edge of the specimen and result in overall failure.In Figure 6, the various crack types initiated from the preexisting flaws identified by Yang et al. [28] and Wong and Einstein [43] are shown.Furthermore, substantial experimental efforts have illustrated cracks with similar characteristics in the specimens with single fissure [10,11,34,77,78,99,[103][104][105].
Apart from the single joint, the coalescence pattern between cracks has been investigated via experiment tests as well.For a specimen containing two or three fractures, some studies focused on crack coalescence between parallel fractures [12,14,18,22,23,40,82,101,[106][107][108][109][110][111][112].In Figures 7(b) and 7(c), the experimental results from Wong and Einstein [15] and Sagong and Bobet [14] are shown, respectively.Wong and Chau [22] presented nine types of coalescence pattern between two parallel fissures (Figure 7(c)).Overall, the coalescence pattern mainly concentrates on tensile pattern, shear pattern, and mixed pattern.For the tensile or shear pattern, the coalescence results from the tensile or shear crack propagation and linkage.However, for the mixed pattern, it is the combination of tensile and shear crack.At the same time, there are many kinds of combinations.
For the unparallel fissures, Lee and Jeon [17] investigated the penetration mode between two unparallel fissures by using PMMA, Hwangdeung granite, and Diastone.e unparallel fissures were created through the water-jet system.e unparallel fissure is the combination of a horizontal fissure and an inclined fissure underneath, and the inclined fissure had an inclination angle of 30 °-90 °in 15 °intervals.Under uniaxial compressive loading, coalescence occurred mainly through the tensile cracks or tensile and shear cracks (Figure 8).In addition, the crack initiation and coalescence stresses were analysed.e results indicated that the shielding effect of the horizontal fissure plays an important role in the stresses, and the stresses reached a critical point when the inclination angle is 60 °.
Except for the literatures mentioned above, there are also laboratory tests conducted by other scholars to investigate the failure behavior of the specimen with unparallel fissures.Zhang et al. [69] studied crack   (a) Hwangdeung granite [17].(b) Sand stone [18].
Advances in Civil Engineering coalescence between two nonparallel flaws; five types of linkage were observed between two flaws: tensile crack linkage, tensile crack linkage with shear coalescence at the tip, shear crack linkage, mixed linkage, and indirect crack linkage.Using photographic monitoring and acoustic emissions monitoring techniques, Yang et al. [113] investigated the relationship between the real-time crack coalescence process and axial stress-time behavior for red sandstone containing two unparallel fissures.
In addition to studies on parallel or unparallel fractures, the mechanical behavior of multifissure specimens under uniaxial loading are also well discussed in the literature.As mentioned above, two types of cracks initiated from the tips of fissure.With further loading, new cracks propagate and link with others to form overall failure.In the previous studies, many kinds of joint geometry are also considered, such as joint inclination angel [53,71,80,88], joint distance [53,71,80,88], and overlap distance [88].Cao et al. [53] (a) 4 Advances in Civil Engineering conducted laboratory tests on physical models of rock with nonpersistent joints.e specimen is made of cement mortar, and the dimensions (height × width × thickness) are 200 × 150 × 30 mm. e preexisting fissures in the specimens are created by inserting metal shims into the fresh mortar and removing it after it becomes hardened.ere are four kinds of inclination angles in the specimen: 25 °, 45 °, 75 °, and 90 °.And, for each angle, the fissure number varied from 5 to 20 at 5-unit increments.e experimental test results indicated that the strength of the specimen increases with the increasing of the joint inclination angle and decreases with the increasing of the fissure number.Moreover, the failure modes of specimens can be classified into four categories: mixed failure, stepped path failure, Shearing failure, and intact failure (Figure 9).
For the specimen with different inclination, Yang's [88] result indicates that the jointed rock mass is in significant dependence of the joint orientation, and four kinds of failure modes are identified: tensile failure across the joint plane, shear failure along the joint plane, tensile failure along the joint plane, and intact material failure.For the ubiquitous joint, Cao et al. [20] identified four types of failure pattern through experimentation on rock-like specimens, the failure patterns of the multiple jointed specimen including stepped path failure, planar failure, shear-I failure, and shear-II failure (Figure 10).

Biaxial and Triaxial
Loading.As mentioned above, the failure process of jointed rocks or rock-like materials has been fully investigated by scholars.However, the confining pressure also has a great influence on the mechanical behavior of jointed rock mass.For the failure behavior considering a confining pressure, experimental studies have been conducted via biaxial or triaxial loading.Bobet and Einstein [11] investigated the fracture coalescence in the gypsum specimen with two open or closed preexisting fissures under biaxial compression.e tests showed wing cracks appeared under unconfined and slightly confined compressive loads and disappeared entirely at high confining stresses.Mughieda and   e experimental results showed three basic failure modes: failure through a planar surface, stepped failure, and failure by the rotation of new blocks (Figure 12).Furthermore, the strength and deformation of the jointed specimen are related to the failure mode.To be specific, the specimen with planar failure and stepped failure exhibits higher strength and small failure strain.However, the rotational failure mode is associated with low strength and large deformation.

Advances in Civil Engineering
For the triaxial compression, there are mainly two types of loading methods, through either a pressure or displacement boundary condition [114].As shown in Figure 13, for the conventional triaxial test, the specimen is cylindrical and has been put inside a pressure chamber and isotropically loaded by hydrostatic pressure using various confining fluids (e.g., air, water, and hydraulic oil).
Huang et al. [39] conducted conventional triaxial compression experiments on sandstone specimens with two closed nonoverlapping fissures.As shown in Figure 14, four types of sandstone specimens were tested.Moreover, in order to form closed fissures, the mixture of water, cement, and sand has been used to fill up the open fissures.e test results show that the arrangement of the fissure has a significant influence on the deformation, strength, and failure pattern of the specimen.Compared with the stress-strain curves of the intact specimen (Type A in Figure 14), the postpeak curves obtain from the precracked rock specimens exhibit obvious fluctuation, especially for Type D in Figure 14.For the Type B, the stress-strain curves exhibit Z-shape characteristics and show a double-peak stress.Moreover, the influence of fissure geometry is greater than that of confining pressure.In addition, the stress for crack initiation, dilation, and the peak strength of jointed sandstone specimens is significantly lower than that of intact rock.
Using a combination of experiments and PFC (3D), Hang et al. [115] investigated the internal damage behavior of rock under conventional triaxial compression tests.As shown in Figure 15, there are two unparallel fissures in the cylindrical specimens (50 mm diameter and 100 mm height).And, four kinds of confining pressures have been considered: 0 MPa, 5 MPa, 15 MPa, and 25 MPa.e results show that the crack  Advances in Civil Engineering evolution and failure characteristics are dependent on the fissure geometry and confining pressure.At the same time, when the confining pressure is low, the failure mode of the specimen is mainly affected by fissure geometry.When the confining pressure is high, the effect of confining pressure on the failure mode is greater than that of fissure geometry.
However, for the true triaxial compression, there are few studies focused on the specimen with fissure or joint.Most of the literatures concentrated on the mechanical behavior of intact rocks under true triaxial compression.Chen and Feng [116] carried out true triaxial unloading tests on new granite specimens, and the rock triaxial strength has also been discussed.With unloading, the minor principal stress, Du et al. [117] investigated true triaxial strength and failure modes of cubic rock specimens (Figure 16).e triaxial tests have been carried on granite, sandstone, and cement mortar cubic specimens.Under true triaxial unloading condition, the test results show that both the strength and failure modes of cubic rock specimens are affected by the intermediate principal stress.e end effect of specimen seems not to play an important role in strength and failure modes of cubic specimen under triaxial tests.For the hard rock, when σ 2 increases to a critical value, the fracture angle of the specimen may change from shear to slabbing failure.However, for medium strong and weak rocks, they exhibit shear failure with a large amount of plastic deformation.

Tension Loading.
For the tension loading, there are two types of standard methods: direct tension [118,119] and Figure 8: Sketch of the observed cracks around the bridge area in the Diastone specimens [17].
indirect tension [119,120] tests.And, these two kinds of methods were suggested by the ISRM and ASTM for determining the quasistatic tensile strength of rock materials.
For the crack initiation under tension loading, the Brazilian tensile test is the most common in previous works.Khanlari et al. [121] conducted Brazilian tensile tests on laminated sandstone to investigate the tensile strength and failure pattern of jointed rock mass under tension loading.Li et al.(c) Fissure geometry and loading condition [115].
Advances in Civil Engineering 11 [122] investigated the effects of loading direction on failure load and failure modes for Brazilian tests on jointed coal rock.Using a combination of experiments and numerical simulation, Haeri et al. [62] investigated the crack propagation and coalescence in rock-like disks.ere are two kinds of jointed specimen made by inserting thin metal slim into the fresh mortar.For the specimen with a single crack, the width and thickness of the crack are 30 mm and 1 mm, respectively.And, the inclination of preexisting crack is 0 °, 15 °, 30 °, 45 °, 60 °, and 90 °.For the double-cracked specimens, the width and thickness of the crack are 20 mm and 1 mm, respectively.One crack at the horizontal, and another is oriented at different angles with respect to the horizontal (0 °, 30 °, 60 °, and 90 °).e test result shows that the crack orientation has a great influence on the failure load.In Figure 17, the failure patterns for precracked specimens under the Brazilian tensile test are shown.For the crack initiation, the fissure orientation has a significant influence on the crack initiation angle.And, the coalescence mode between the fissures is mainly tensile mode.
Compared to the indirect tension test, because of the misalignment and stress concentration during the test, the direct tension tests are difficult to perform.But, there are few people who used this method to investigate the failure behavior of jointed rocks or rock-like materials.Yang et al. [123] conducted uniaxial tension experiments on precracked rocklike specimens (Figure 18); the strength and failure behavior of specimen have been analyzed and discussed.e results indicated that crack geometry parameters including crack dip angle, crack spacing, and crack intensity, have significant influence on the strength and failure modes of the samples.

Shear Loading.
As mentioned above, the crack initiation, propagation, and failure characteristics of jointed rock or rock-like specimens have been investigated by many scholars.Compared with the compressive tests, relatively few experiments were done to investigate the pattern of crack coalescence under direct shear loading [67,68,83,84,[124][125][126][127][128][129].Under shear loading, the failure behavior of jointed rock mass is different from that under compression.Lajtai [124,125] conducted direct shear testing on rock-like specimen with nonpersistent fissures, and the test result indicated that the failure mode changes with increasing normal stress.At the same time, he also proposed a composite failure envelope to describe the changes of strength.Savilahti et al. [126] also investigated the failure behavior of jointed rock-like specimen under direct shear loading, and the influence of joint separation and overlapping on failure behavior of the jointed specimen have been analyzed.Sarfarazi et al. [67] has studied the effect of joint overlap on the full failure behavior of a rock bridge in the direct shear test via combination of the experimental test and numerical simulation.Four kinds of specimens have been tested; all of the ligament length are kept at 45 mm and the ligament angle is 0 °, 25 °, 90 °, and 115 °, respectively.During testing, the normal stress was set at 0.1 MPa, and the test results show that the failure stress decrease with increasing of the ligament angle.In Figure 19, the failure modes for specimen with different ligament angles and most of the specimen failure from the propagation of tensile crack are shown.
For the failure characteristics of the specimen with multifissures under direct shear loading, the breakage and shear behavior of intermittent joints (Figure 20) have been investigated by Gehle and Kutter [83].e results show that the failure of specimen can be divided into three phases; both the geometrical parameters and loading conditions have been found to influence the activated shear resistance in each phase.Moreover, the mechanisms which govern the different shear phases could be identified as (1) tensile rupturing, (2) rolling and sliding friction of dilatant joint zones, and (3) sliding within the joint filling composed of brecciated material.
Apart from the direct shear loading, there is another shear test method, namely, restrictive shear test.In this loading, the compressive load and shear load on the specimen increases with the increasing of the load P. It can be used to investigate the failure process of jointed rock mass in compressive and shear environment, especially for the rock mass in rock slope.In Figure 21, the stepped failure in rock slope is shown.
Actually, there are relatively few scholars who conducted experiments to investigate the failure process of jointed rock mass under restrictive shear loading.Zhang et al. [21] has investigated the mechanical behavior of rock-like specimen mixed flaw, and the strength, fragmentation, and fractal properties have been discussed.For the mixed flaws, the length of edge-notched flaw and imbedded flaw is 10 and 30 mm, respectively.e inclination of imbedded flaw changes from 0 °to 90 °with an increment of 15 °.Moreover, with two kinds of shear angle (45 °and 60 °), the jointed specimen are loaded under compressive-shear loading until failure.Based on the experimental results, three different patterns of tensile cracks and shear cracks are observed (Figure 22).At the same time, it has been found that the shear strength is a function of the flaw geometry and the shear-normal stress ratio.

Cycle Loading.
Recently, for the cycle loading, most studies in the literature have focused on the mechanical behavior of intact rocks or rock-like materials under cyclic loadings.Based on the experimental tests, the hysteresis of 12 Advances in Civil Engineering the stress-strain curve has been revealed [131][132][133][134].At the same time, the fatigue strength and deformation of the intact specimen have also been investigated by scholars, and the results indicated that the dynamic strength and elastic modulus decrease exponentially with the increase of cycles [135][136][137][138][139].
As mentioned above, the joint or preexisting fissure has a significant influence on the failure behavior of rock mass.
e dynamic response of jointed rocks will be different from their static properties, and the mechanical behavior of jointed rock mass has also been investigated by scholars.e previous studies [140][141][142] reported that jointed rock mass are very sensitive to the cyclic loading and the joint confirmation has a significant influence on the dynamic strength and deformation behavior.Erarslan and Williams [143] conducted the static and cyclic loading test on inclined cracked chevron notched Brazilian disc (CCNBD) specimens, and the experimental results show that the failure load obtained through cyclic loading decreased between 30 and 45% compared with those in static loading.Liu et al.Advances in Civil Engineering 13 [63] performed the cyclic uniaxial compression test on jointed rock-like specimens, and the influence of joint geometric parameters such as dip angle, persistency, density, and spacing on mechanical properties of intermittent jointed specimen have been investigated.e test result indicated that the stress-strain curve of jointed rock under cyclic loadings is dominated by its curve under monotonic uniaxial loadings.At the same time, under cyclic loading, two types of cracks were observed in the jointed specimen.To be specific, shear cracks mainly occur in the specimen with higher joint inclination angle or higher persistency.

FEM or XFEM.
Based on the finite element method (FEM) and nonlinear dynamics method, Li and Wong [173] has investigated the influence of inclination angle and loading condition on crack initiation and propagation.e numerical results indicated that the loading condition has influence on crack initiation sequences and the overall crack pattern.Under a relatively low loading rate or a small magnitude of maximum loading pressure, tensile crack appears first.However, under a relatively high loading rate and a large magnitude of maximum loading pressure, shear crack would occur first (Figure 23).For the closed fissure, Xie et al. [174] used XFEM to investigate the crack initiation and propagation in rock-like material with closed fissure under uniaxial compression.e numerical simulation results show that (1) for the specimen with inclination angles 30 °and 45 °, minor effects are exerted by crack surface friction on the stress distribution around the fissures, and the effect are much more obvious when inclination angle is 60 °(Figure 24); (2) when the inclination angle is 45 °, it is the most favorable value for crack propagation; (3) the friction seems only to play a minor role on the initiation location and angle of the wing cracks, but the friction has a great influence on the propagation length.
RFPA2D is developed by Northeastern University, People's Republic of China, and it can be used to model the evolution of damage in brittle materials by allowing the linear elastic elements to fail in a brittle manner.In recent years, this method has been used for modelling progressive failure in rocks or jointed rock-like materials.Based on RFPA2D, Tang and Kou [145] presented two particular cases concerning crack propagation and coalescence in brittle materials (the model containing a row of small flaws and several larger flaws and the model containing randomly distributed homogeneities as shown in Figure 25).e numerical results show that, under axial compression, wing crack occurs at the tips of the preexisting flaws, and propagation occurs along with the direction of maximum far-field compression.At the same time, the coalescence between preexisting flaws may be in tensile, shear, and combination of tensile/shear pattern.With a confining pressure, the crack becomes stable.However, lateral tensile stress has a significant influence on the crack growth; even a small value will result in unstable growth.For the specimen containing homogeneities on a grain scale, the numerical results indicated that the failure characteristics strongly depend on the mechanical and geometric properties of the grains and inclusions.For the failure behavior of the jointed specimen under axial compression, Tang et al. also conducted similar numerical studies through RFPA2D, such as crack coalescence in rock-like materials containing three flaws [108]; the numerical results also show qualitatively a reasonably good agreement with reported experimental results.Wong et al. [175] conducted numerical simulation by using RFPA2D to investigate the splitting failure in brittle rocks containing open joints under uniaxial compressive loading.And, the influence of preexisting joint length on wing crack growth has been studied.
Apart from the axial compression, RFPA2D can also be used to investigate the failure process of the jointed specimen under shear loading.Zhang et al. [176] has investigated the shear behavior of rock specimens with several intermittent joints (Figure 26).Based on the numerical simulation, the whole failure process and the failure patterns are observed.For the failure pattern under shear loading, it is mostly affected by joint geometry parameters, and the shear strength of the specimen is related to the failure pattern.Moreover, both the joint separation and azimuth angle have influence on the wing crack propagation, and the wing crack dominates the overall failure of the specimen.Furthermore, the results also show that the macroshear crack is the result from the accumulation of microtensile damage.
Da Huang et al. [31] used the ANSYS AUTODYN-2D and investigated the crack initiation and propagation in three types of precracked sandstone specimens under conventional triaxial compression.In Figure 27, the comparison between numerical and experimental results is shown.e coalescence pattern in numerical results shows a great agreement with that in experimental results.e result also indicated that the crack initiation stress, critical stress of dilation, and peak strength for precracked specimens are far below than those in intact rocks.At the same time, all of them increase with the confining pressure.

DEM.
In recent years, the DEM has been widely used by scholars to model the failure process of the jointed rocks or brittle rock-like materials.DEM has experienced decades of development, and many kinds of DEM software have sprung up in recent years, such as PFC2D/3D, UDEC, and 3DEC.For crack initiation and propagation, the PFC has been widely accepted by scholars, and the numerical results show a great agreement with the experimental results.
Based on the parallel bond model in PFC2D, Zhang et al. [177] has investigated the crack initiation and propagation under uniaxial compressive loading.e numerical simulation result shows that the inclination of joint has a strong influence on the crack initiation and propagation behavior.And, by analyzing the parallel bond forces and displacement fields, the crack initiation location has been identified.Moreover, two types of displacement fields, namely, type I (DF_I) and type II (DF_II), have been proposed to distinguish the tensile crack and shear crack (Figure 28).In Figure 19, two displacement field types associated with different microcracking processes are shown.For the type I (DF_I), if the two displacement trend lines diverged from each other and there Advances in Civil Engineering was very limited relative shear movement, relative tensile displacement occurred dominantly in the region in between.For the type II (DF_I), if the two displacement trend lines exhibited both a relative tensile displacement and a shear displacement in the region in between.Based on the parallel bond model, the crack initiation, propagation, and  16 Advances in Civil Engineering coalescence in between two fissures have been further investigated by Zhang et al. [111].Lee and Jeon [17] have investigated the crack initiation, propagation, and coalescence between unparallel fissures under uniaxial compression.e relationship between crack initiation, coalescence, and failure stress and fissure inclination angle in the single and double crack has been analysed.Apart from the specimen with open or closed fissure, the failure characteristic of multifissure specimen has also been extensively discussed in literatures.Cao et al. [53] used the PFC2D to model the failure behavior of multifissure specimen under uniaxial compression, and the fissure geometry parameters such as inclination (α) and fusser number (N f ) have been taken into consideration.Based on the numerical results, the influence of fissure geometry on peak strength of multifissure specimen has been discussed.To be specific, the peak strength increased with increasing of inclination angle.e material strength was lowest for inclination angle 25 °, and highest for 90 °. e influence of N f on the peak strength depended on α.For α � 25 °and 45 °, N f had a strong effect on the peak strength, while for higher α values, especially for the 90 °sample, there were no obvious changes in peak strength with different N f values.At the same time, the failure modes in the simulated results also agree very well with those in experimental results.
Fan et al. [172] usd the PFC3D to investigate the macromechanical behavior of jointed blocks with multinonpersistent joints under uniaxial loading (Figure 29).e effect of joint inclination, size, and joint mechanical properties on the strength, deformability, stress-strain, and failure modes have been studied.e simulation results indicated that the joint particle stiffness only play a minor to a significant role on peak strength depending on inclination and joint density.And, the joint particle stiffness plays a negligible role on deformation modulus depending on inclination and joint density value.For the stress-strain curves, they can be further classified into four types.Moreover, the failure modes such as (1) splitting failure; (2) plane failure; (3) stepped path; and (4) intact material failure occurring in numerical results also agree well with those in experimental results.
By use of numerical direct shear tests (PFC2D), Sarfarazi et al. [67] has investigated the effect of joint overlap on the failure process of the jointed specimen.In Figure 30, the planar nonpersistent joints and en echelon e ligament angle is 0 °, 25 °, 90 °, 115 °, and 140 °, and the ligament length is kept constant for all specimens (18 mm).Moreover, the lengths of the edge-notched joints are 21, 21.8, 24.2, 33.8, and 36.9 mm for ligament angles of 0 °, 25 °, 50 °, 90 °, 115 °, and 140 °, respectively.e simulation results indicated that the rock bridge lost their capacity when nearly 50 % of total cracks developed within the rock bridge.Furthermore, the simulation results also show that the macroshear failure in rock-bridge results from the cluster of microtensile cracks.
Based on the discrete element software UDEC, Vergara et al. [178] studied the mechanical behavior of rock containing parallel nonpersistent joint (Figure 31).And, the results indicated the large anisotropy in the strength resulting from variation of the joint orientation and lower strength of the specimens was caused by the coalescence of fractures belonging to parallel joint sets.
Jin et al. [179] used 3DEC software to simulate the anisotropic mechanical behaviors of columnar jointed basalt under compression (Figure 32).e results show that, with different primary joint angles, the jointed specimen exhibits different mechanical properties under the presence of confining pressure.e primary joints play an important role on the anisotropy of the jointed specimen.At the same time, it is found that increasing confining pressure can reduce the influence of columnar joints on the anisotropy of mechanical properties obviously.

FEM-DEM.
As a kind of heterogeneity material, natural rocks containing fracture and the fragmentation processes limit the applicability of continuum-based models to model the failure behavior of natural rock mass in rock engineering.In order to resolve or improve these limitations, Munjiza et al. [180] proposed the finitediscrete element method (FEM-DEM).For the FEM-DEM method, the discrete element has been discretised into finite elements.It is indicated that there is a finite element mesh associated with each discrete element.
en, the For example, Lisjak et al. [182] used the finite-discrete element method (FEM/DEM) to model the mechanical behavior of layered materials (Figure 33), and the results also demonstrated the effectiveness of the finite-discrete element method in simulating the short-term mechanical response of the Opalinus clay.

Other Numerical Methods.
e Fast Lagrangian Analysis of Continua 3D (FLAC3D) has also been used by scholars to investigate the crack initiation and propagation in jointed rocks or rock-like materials [98,[185][186][187][188].For example, based on the FLAC3D, Fu et al. [98] proposed a new kind of elastic-brittle model to model the crack initiation and   e uniaxial and biaxial numerical results show an excellent consistency with experiment results.In Figures 34(a) and 34(b), the 3D views of the crack propagation process under uniaxial loading and biaxial loading, respectively, are shown.To be specific, under uniaxial loading, secondary cracks are all wing cracks and formed at both ends of the precrack's major axis.Under biaxial loading, (the lateral stress is 20% of the peak stress), a few antiwing cracks appear.Accordingly, only in high lateral stress, antiwing cracks will appear.
Combined with maximum tensile stress criterion and the Mohr-Coulomb criterion, the extended nonordinary statebased peridynamics has been used for modelling the initiation, propagation, and coalescence of the jointed rocks under compressive loads by Wang et al. [189][190][191][192][193][194][195][196].Different types of cracks includes wing crack, oblique secondary crack, quasi-coplanar secondary crack, and antiwing crack are modeled and distinguished by the proposed method, and the numerical results show a great agreement with the previous experimental ones (Figure 35).
Pramanik and Deb et al. [161] used a kind of methodology which is developed in the SPH framework to investigate the failure behavior of rock material containing multiple discontinuities or joints.In this method, the joint is represented by a set of particles at the location of the joints.At the same time, based on the Drucker-Prager yield criterion, the free-sip, no-sip, and symmetric boundary conditions are also implemented in this method.In Figure 36 the failure mode of the specimen with different joint inclination is shown, and the numerical results are in good agreement with derived theoretical results.e efficacy of the numerical method is successfully demonstrated by two samples under uniaxial and gravitational loading conditions.Moreover, this method also has shown promises to model the failure process of jointed rock mass in three dimensions.For this kind of method, more examples have been discussed in previous works [159][160][161][162].
Yao et al. [197] proposed an extended rigid block spring method (RBSM) to model the damage and failure of anisotropic rock mass.In this method, the tensile failure Advances in Civil Engineering between interface occurs when normal stress reaches the tensile strength.However, the shear failure is described by a nonlinear criterion in terms of the local normal and shear stresses.In Figure 37, the failure modes of rock mass under different confining pressure are shown.e numerical results indicated that the inherent bedding plane and confining pressure have a significant influence on the macroscopic mechanical strength and failure mode of rock mass.In addition, the numerical results agree well with typical experimental data for both elastic properties and mechanical strength.Meng et al. [198][199][200][201] formulated a series of discrete numerical mechanics models as standard second-order cone programs.Advanced optimisation algorithms can be employed to solve the problem.It is necessary to note that the purely static discrete element method can be employed.Furthermore, both the hard-sphere and soft-sphere discrete numerical models can be recovered.Notably, apart from the numerical methods mentioned above, there are also other methods developed by scholars and been used to model the failure process of rocks or rock-like materials.However, they cannot be listed all here one by one.

Conclusion
e discontinuities have a great influence on the mechanical behavior of rock mass, and under loading, the failure of natural rock mass results from the crack propagation and coalescence in rock mass.Compared with the intact rocks, the jointed or fractured rock mass usually exhibits weaker and highly anisotropic mechanical characteristics.
e experimental results indicated that the joint confirmation parameters have a significant influence on the mechanical behavior of rock mass.e strength parameters of jointed rocks decrease with the increasing of joint number, length, and persistent degree.At the same time, the strength of the jointed rocks usually has the lowest value when inclination is around 45 °.For the failure characteristics of jointed rocks, the coalescence between joints includes tensile, shear, and mixed modes.
e numerical simulation has been viewed as a kind of economical and practical method and has also been accepted by many scholars.Compared with FEM or XFEM, the DEM exhibits strong advantage in crack initiation and propagation.In recent years, the DEM has been widely used to model the failure behavior of jointed rocks or rock-like materials, and most of the numerical values are in good agreement with experimental results.Moreover, other kinds of numerical methods such as NMM, smoothed-particle hydrodynamics, peridynamics, and RBSM have also been successfully applied to the analysis of failure process of jointed rocks.Natural Science Foundation of China (2018JJ2500), Open Fund of Engineering Laboratory of Spatial Information Technology of Highway Geological Disaster Early Warning in Hunan Province (Changsha University of Science & Technology) (Nos.kfj70603, kfj170407, and kfj170406), and Natural Science Basic Research Plan in Shaanxi Province of China (2018JQ4015).e authors wish to acknowledge these supports.[197].
quasi-coplanar secondary cracks Type VI: oblique secondary cracks and wing crack Type II: Quasi-coplanar secondary cracks and out of plane tensile secondary crack Type VII: Oblique secondary cracks and out of plane tensile secondary crack Type III: Quasi-coplanar secondary crack and wing crack Type VIII: Oblique secondary cracks Type IV: Wing crack Type IX: Oblique secondary crack and quasi-coplanar secondary crack Type V: Quasi-coplanar secondary crack and out of plane secondary shear crack I

mm 100mm Figure 14 :Figure 15 :
Figure 14: Four types of tested sandstone specimen and geometry of two preexisting nonoverlapping filled flaws: type A is intact specimen; types B-D are specimens with two different preexisting flaws [39].

Figure 17 :Figure 18 :
Figure 17: Experimental results illustrating the coalescence path of rock-like disk specimens containing double cracks [62].

Figure 22 :
Figure 22: Restrictive shear test, failure patterns, and peak shear strength of jointed specimen [21].(a) Layout of the loading system.(b) Patterns of crack propagation and coalescence.

Figure 20 :Figure 21 :
Figure 20: Rotating mechanism of the rock bridges in the second phase of shearing [83].

Figure 23 :Figure 24 :
Figure 23: Initiation and propagation of cracks under different loading rates [173].(a) A relatively low loading rate.(b) A high loading rate.

Figure 25 :Figure 26 :
Figure 25: e numerical models for the simulations.Model containing (a) a row of small flaws and several larger flaws and (b) randomly distributed homogeneities [145].

Figure 37 :
Figure37: Failure modes for specimens with various orientations under compression with different confining pressures[197].

Table 1 :
Major materials for the specimen in testing the failure behavior of jointed rocks.