
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39261626
72369
10.1038/s41598-024-72369-0
Article
Dynamic fracture mechanism and fragment characteristics of sandstone specimens with asymmetrical conjugate fissures under static pre-compression
Feng Peng fengpengfx@163.com

123
Cao Pan 1
Li Juntao 1
Tang Ran 13
Li Huajin lihuajin@cdu.edu.cn

13
1 https://ror.org/034z67559 grid.411292.d 0000 0004 1798 8975 School of Architecture and Civil Engineering, Chengdu University, Chengdu, Sichuan China
2 https://ror.org/01t001k65 grid.440679.8 0000 0000 9601 4335 Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education, Chongqing Jiaotong University, Chongqing, China
3 https://ror.org/034z67559 grid.411292.d 0000 0004 1798 8975 Sichuan Engineering Research Center for Mechanical Properties and Engineering Technology of Unsaturated Soils (Chengdu University), Chengdu, China
11 9 2024
11 9 2024
2024
14 212681 3 2024
6 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
This paper systematically investigates the dynamic mechanical response and fragment characteristics of sandstone with asymmetrical conjugate fissures subjected to preexisting static stress based on the split Hopkinson pressure bar apparatus. The cross-fissured sandstone exhibits significant variations in mechanical behavior when exposed to higher dynamic strain rates under identical static pre-stress. The cross-fissured sandstone with a higher dynamic strain rate is characterized by a greater coupled strength under the same static pre-stress; for a given dynamic load, the highest coupled strength occurs under the static pre-stress of 60% UCS. The failure mode of the cross-fissured sandstone is predominantly governed by the dynamic strain rate, independent of static pre-stress variations considered in this investigation. At lower dynamic strain rates, the specimens typically exhibit a mixed tensile-shear failure mode, characterized by the dominance of larger fragments in the broken specimens. In contrast, under high dynamic impacts, the sandstone tends to fail in a shear-dominated manner, resulting in smaller fragments with a more uniform size distribution. Furthermore, the study explores how varying dynamic strain rates and static pre-stress influence the fragment characteristics of the cross-fissured sandstone. Higher dynamic strain rates and increased static pre-stress generally lead to smaller mean fragment sizes. This phenomenon is quantitatively described by fitting fragment size distributions using the Generalized Extreme Value (GEV) distribution, revealing a decrease in the location parameter (μ) and an increase in fractal dimension. These metrics indicate that higher dynamic strain rates and static pre-stress result in sandstone specimens breaking into smaller fragments with a more homogeneous size distribution. These findings contribute to a deeper understanding of rock dynamics, with potential implications for engineering applications involving similar geological configuration under dynamic loading conditions.

Keywords

Mechanical response
Failure mechanism
Fragmentation characteristic
Asymmetrical conjugate fissures
Dynamic loading
Subject terms

Petrology
Mechanical properties
The Open Research Subject of Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education, Chongqing Jiaotong UniversitySLK2023A02 Feng Peng Natural Science Foundation of Sichuan Province, China2024NSFSC0983 2023NSFSC0789 Feng Peng Li Huajin Open Fund of Sichuan Engineering Research Center for Mechanical Properties and Engineering Technology of Unsaturated SoilsSC-FBHT2024-06 Feng Peng issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

After experiencing the long-term geological activity, buried field rock masses generally contain myriads of cross fissures, which significantly influence the mechanical response of rock masses when suffering exterior dynamic loads1. For instance, as shown in Fig. 1, buried field rock masses with cross fissures under static tectonic stress can be influenced by the external dynamic disturbance from excavation, and the propagation of cross fissures affects the mechanical properties and fragment size distribution of rock masses2. In this case, static tectonic stress corresponding with static energy accumulation within the rock can be regarded as the intrinsic cause of rock failure after dynamic disturbance. So technically speaking, both the static pre-stress and dynamic disturbance influences the failure of deep rock masses3–5. Thus, it is necessary to investigate the effect of dynamic strain rate on the mechanical properties and fragment characteristics of rocks with cross discontinuities under static pre-stress.Fig. 1 Schematic diagram of the field rock masses with cross fissures under static tectonic stress and external dynamic disturbance from excavation.

For the study of cross-fissured rocks, existing literatures mainly focus on static loading situations6. Liu et al.7,8 studied the mechanical behaviors of molded gypsum specimens with one set of cross fissures under uniaxial compression and biaxial compression, and they concluded nine crack types. Zhang et al. and Zhou et al.9–12 further conducted a series of uniaxial compression tests on rock-like materials containing one or two sets of cross fissures with different fissure configurations and they analyzed the mechanical response and crack coalescence behaviors. The preexisting cross fissure can significantly weaken the strength of rock-like materials, and the cross-fissure geometry influences the rock bridge coalescence patterns13. Zhou et al.14 systematically investigated the difference between two-dimensional and three-dimensional cracking behaviors of PMMA specimens with cross‑embedded fissures, involving crack growth, coalescence and wrapping. In addition, Cao et al.15 experimentally and numerically studied the mechanical behaviors of rock-like materials with multiple cross fissures, including the uniaxial compressive strength and failure patterns of the specimens (stepped path failure pattern, planar failure pattern, and two other shear type failure patterns). Existing literatures on the cross-fissured rock specimens under static loading are summarized in Table 1. However, the research for cross-fissured specimens under dynamic loads is rare compared to static loading conditions. Feng et al.16 studied the effects of individual dynamic strain rates on the mechanical behaviors of sandstones with symmetrical and asymmetrical cross fissures based on the Split Hopkinson Pressure Bar (SHPB) loading system. Table 1 The summary of existing literatures on the cross-fissured specimens under individual static loading.

Author (year)	Brittle materials	Cross fissure types	Loading path	
Zhang et al.9	Rock-like materials	One set of cross fissures with different inclination angles	Uniaxial static compression	
Zhang et al.10	Rock-like materials	Two aligned X-type flaws	Uniaxial static compression	
Zhang et al.11	Rock-like materials	Two sets of cross fissures with different aligned stepped, and collinear configurations	Uniaxial static compression	
Liu et al.7	Gypsum	One set of cross fissures	Uniaxial static compression	
Liu et al.8	Gypsum	One set of cross fissures	Biaxial static compression	
Zhou et al.12	Rock-like materials	Two sets of cross fissures	Uniaxial static compression	
Zhou et al.14	PMMA	Two sets of cross fissures with different inclination angles and horizontal spacings	Uniaxial static compression	
Liang et al.13	Mudstone	Two sets of cross fissures with different cross angles and rock bridging angles	Uniaxial static compression	
Cao et al.15	Rock-like materials	Multiple sets of cross fissures	Uniaxial static compression	

Under a dynamic loading case, the load increases at a higher rate, and the deformation and fracture propagation occurs in a limited time. Thus, the fissured rock might exhibit more complicated mechanical response and fracture behaviors under dynamic loading than static loading17. Due to more energy absorption, rock materials under dynamic strain rate are generally broken into irregular-distributed fragments in a certain scale range18,19. Essentially, the formation of rock fragments experiences a complex process, which is related to the exterior loading condition and inherent constitutive behavior20,21. After suffering dynamic impact, the recovered rock fragments can be regarded as following a statistical distribution, such as Rosin–Rammler distribution22, Gate-Gaudin-Schuhmann distribution23, Weibull distribution24, and generalized extreme value (GEV) distribution25. Moreover, the fragments with fractal-based structures generally exhibit a self-similarity characteristic, where the structures show a similar pattern with different fragment scales. The fragment size can be described by the fractal behavior, and the distribution of fragment sizes can be quantitatively characterized using the fractal dimension26, a higher fractal dimension generally corresponds with more fragments with smaller volumes and higher crush degrees27.

So far, research focused on the dynamic fragmentation characteristics of fissured rocks is limited, especially for rock specimens under coupled static-dynamic loads. In this study, we systematically investigate the dynamic mechanical behaviors of cross-fissured sandstones under static pre-compression based on the SHPB dynamic loading system. In addition, the fragments of the fractured specimens are carefully collected, and a series of sieving tests are performed to reveal the fragment characteristics of broken rocks based on the fractal theory and statistical principle.

Specimen preparation and experimental equipment

A relatively homogeneous sandstone material from Neijiang, Sichuan province of China, is selected for specimen preparation in this experimental campaign. The sandstone is composed of feldspar (43%), quartz (35%), rock fragments (20%) and clay minerals (2%). All the rock specimens are prismatic in appearance, and their geometric dimensions (width × thickness × height) are 35 mm × 35 mm × 45 mm, as shown in Fig. 2a. To avoid the influence of potential anisotropy on experimental results, all the rock specimens are processed to have consistent loading directions relative to the virgin sandstone block. The preparation process of the fissured specimens is as follows: (a) place the relatively homogeneous sandstone block on a rock-cutting machine, and cut it according to the preset size to obtain the intact specimen; (b) use a high-strength alloy to process a micro-hole with 1.5 mm diameter at the geometric center of the surface where the width and height are located; (c) employ a 1-mm-thick diamond wire saw to cut the fissures with 7 mm length and desired inclination angles; and (d) polish the specimen surfaces to control their roughness to be within 0.02 mm. Three specimens are prepared for each set of experimental parameters. In this manuscript, fissure represents the preexisting through-going fractures of the specimen, and flaw stands for the generated cracks within specimens during dynamic loading process. The configuration of asymmetrical cross fissures can be defined via three geometrical parameters: fissure length (b = 7 mm), fissure inclination angle I (α = 45°) and fissure inclination angle II (β = 75°). Fissure inclination angle of zero parallel the orientation of applied dynamic load. F1 and F2 represent the tips of the preexisting fissure with inclination angle of 45°, while the F3 and F4 stand for the tips of the preexisting fissure with inclination angle of 75°, respectively (Fig. 2a). Static mechanical parameters of the cross-fissured sandstone acquired through laboratory tests include the uniaxial compressive strength of 51.2 MPa and the elastic modulus of 5.18 GPa. In addition, each rock specimen is notated as R-Dn, where R represents the ratio of static pre-stress to static uniaxial compressive strength of the fissured specimen, and Dn denotes the number of the dynamic strain rate. For example, 0.5-D3 represents the cross-fissured sandstone subjected to the static pre-stress of 0.5 times uniaxial compression strength and the dynamic strain rate numbered 3.Fig. 2 (a) Geometrical characteristics of the asymmetrical cross-fissured specimen and (b) schematic diagram of the SHPB rock dynamic loading system.

A series of coupled static-dynamic compressive loading tests are conducted on the SHPB dynamic loading system, which consists of a striker system, a bar system, a data-acquiring system, and a high-speed camera system, as shown in Fig. 2b. Before loading, cross-fissured sandstones coated with Vaseline are clamped between the incident bar and the transmitted bar. Then, the specimen is axially compressed to the designated stress level in a quasi-static manner. The pre-applied static stress is assumed to be constant throughout the dynamic loading phase of the experiment and added to the dynamic stress history derived from the SHPB system. During the dynamic loading process, a compressive wave (incident wave) is generated from the striker and propagates through the incident bar28. Once reaching the incident bar-specimen interface, the incident wave can be transferred into two parts: one is the transmitted stress wave through the fissured specimen and then enter the transmitted bar, and the other one is the reflected stress wave reflected into the incident bar. The three pulse signals can be recorded via the strain gauges, which are mounted on the surfaces of incident and transmitted bars at the prescribed position and connect to the data acquiring system29,30. Figure 2b shows a typical equal arm bridge, where RG and R0 are the resistance of the strain gauge and the precise standard resistance; ΔU and U0 represent the output voltage and the input voltage of the bridge. If the RG is increased with the value of ΔR, the corresponding strain ε(t) measured via the strain gauge and the voltage signal ΔU(t) can be expressed as follows19:1 εt=2ΔUtK1K2U0

here, K1 and K2 stand for the sensitivity coefficient and the amplification factor of the corresponding strain gauge.

Based on the three recorded pulse signals including the incident wave signal (εi), reflected wave signal (εr) and transmitted wave signal (εt), the strain rate ε˙t, the strain ε(t), and the dynamic stress σ(t) of the cross-fissured rocks can be acquired31,32:2 ε˙t=CLεi-εr-εtr=E/ρLεi-εr-εtrεt=∫0tε˙tdt=CL∫0tεi-εr-εtrdtσt=AE2Asεi+εr+εtr

In the formula, C, A, ρ, and E represent the longitudinal wave velocity, cross-sectional area, density, and the elastic modulus of the elastic SHPB, respectively. As and L stand for the cross-sectional area and length of the cross-fissured rock specimen.

Note that the deformation of cross-fissured sandstones occurs in a very limited time during dynamic loading process. The high-speed camera is thus employed herein to capture the fracture and deformation process on the surface of cross-fissured specimens. The recorded pictures are recorded with resolution of 256 × 256 pixel at a frame rate of 180,000 frames per second (fps). Subsequently, the digital image correlation (DIC) technique is utilized to measure the full-field strain fields of rock specimens after covering speckles based on the captured photos.

Methodology

Sieving tests

A series of sieving tests are conducted after collecting the fragments of broken sandstone in the SHPB loading tests. In this manuscript, the sieve with bore diameters of 0.075, 0.25, 0.5, 1, 2, 5, 10, 20 and 40 mm is utilized. Based on the sieving results, the mean fragment size (M¯) is utilized to represent the damage degree of the fissured rocks, which can be expressed as follows27:3 M¯=∑i=1nmidi/Mt

Here, Mt stands for the total mass of the cross-fissured rock specimen; n is the gradient value of the sieves; mi represents the mass of the fragment size between two adjacent sieves, in which the average diameter of the two sieves is represented by di.

In addition, for the statistical perspective, an exponential-like function, generalized extreme value (GEV) distribution is adopted herein to describe the fragment distribution of rock specimens, which can be expressed as follows33,34:4 Fd;ξ,μ,λ=exp-1+ξd-μλ-1/ξ

where represents the shape parameter; μ is the location parameter, indicating the average size of the rock fragments; λ stands for the scale parameter, describing the range of the fragment size distribution of recovered specimens. The deep analysis on the parameters μ and λ can provide rich information on the fragment size distribution and help to illuminate the strain rate effects on the distribution of rock fragmentation.

Additionally, as an effective tool to reflect the broken degree of rock fragmentation, fractal theory originally defined by Mandelbrot is employed to quantify the irregular characteristics and self-similarity of fragment size distribution35–37. The fractal dimension D can be calculated based on the mass and equivalent size of the broken fragments with the following formula38:5 D=3-k=3-Δlg(Mr<ri/Mt)Δlg(ri/rmax)

where Mr<ri represents the cumulative fragment mass with the associated fragment size less than or equal to ri; k is the slope of the fitting line of lnMr<ri/M0 and lnri/rmax.

Digital image correlation technique

Based on the captured images of rock specimen during dynamic loading process, the digital image correlation (DIC) technique is employed here to calculate the deformation of cross-fissured rock specimen via tracing the difference of each point in the referenced and deformed images, as shown in Fig. 3. The deformation of the point q’ (x’, y’) can be calculated as follows39:6 x′=x+u+∂u∂xΔx+∂u∂yΔyy′=y+v+∂v∂xΔx+∂v∂yΔy

Fig. 3 Schematic diagram of the calculation principle of DIC algorithm.

In this formula, u, v, ∂u∂x, ∂u∂y, ∂v∂x and ∂v∂y stand for the six deformation parameters. u and v represent x- and y- displacement components; ∂u∂x, ∂u∂y, ∂v∂x and ∂v∂y stand for the first-order displacement gradients of the reference subset. The displacement gradient tensor ∇u can be expressed in the form of matrix as40:7 ∇u=∂u∂x∂u∂y∂v∂x∂v∂y

The strain tensor ε can be described as follows:8 ε=∂u∂x12∂u∂y+∂v∂x12∂u∂y+∂v∂x∂v∂y

Note that the arbitrary orientation of the generated cracks during dynamic loading process. Based on the coordinate transformation, the displacement gradients and strain tensor in any direction can thus be acquired as follows:9 ∇u′=∂u′∂x′∂u′∂y′∂v′∂x′∂v′∂y′=TT·∇u·T,ε′=TT·ε·TT=cosαsinα-sinαcosα

Here, T represents the coordinate transformation matrix, α is the angle between transformed and original axis.

Results and discussion

Mechanical response

In the coupled static-dynamic loading test for cross-fissured rock specimen, the stress equilibrium of the two ends of specimen before failure should be guaranteed to avoid premature failure prior to reaching peak stress in the rock specimen33. The momentum trap is an important consideration in SHPB experiments, which was not employed in the present study. In this study, a shaper made via copper sheet with dimension of Φ10*10 mm is utilized to generated a slowly rising half sine wave to promote the stress equilibrium of the specimen during dynamic loading process. Figure 4 compares the incident stress, reflected stress, transmitted stress, and the sum of the incident stress and reflected stress of the two typical fissured specimens (Nos. 0.2-D1 and 0.8-D4) during loading process. The transmitted stress approach a favorable agreement with the sum of the incident stress and reflected stress, indicating that the valid stress equilibrium condition is achieved and maintained for the cross-fissured sandstones in the present coupled static-dynamic loading tests. Figure 5 depicts the stress–strain curves of the cross-fissured rock specimens subjected to varied coupled static-dynamic loads. All the dynamic stress–strain curves exhibit the elastic–plastic deformation characteristics before reaching the peak stress, which can be concluded into two stages. The one is the linear deformation stage, where the stress increases linearly with the increase of strain. The other one is the nonlinear deformation stage, corresponding with the unstable crack initiation and propagation. After peak stress, the stress–strain curves feature an immediate decrease, associated with the macro failure of the fissured specimen. Figure 6 shows the effects of static and dynamic strain rate on the coupled strength of the cross-fissured rock specimens under varied coupled static-dynamic loads; the details are listed in Table 2. In this study, the direct dynamic strength can be acquired from the peak values of dynamic stress–strain curves. The coupled strength is the sum of the static pre-stress and direct dynamic strength. Under the same static pre-stress, the cross-fissured specimens with higher dynamic strain rate can be characterized via the greater coupled strength, featuring evidently loading-rate dependence. However, for a given dynamic strain rate, the static pre-stress shows dual effects on the coupled strength of the fissured sandstones. The data from tests at 0.2, 0.4, and 0.6 UCS converge on a similar trend line, while the 0.8 UCS shows failure at lower dynamic stresses. When the lower static pre-stress applied on the fissured specimen, a part of preexisting micro discontinues are compressed, and the propagation of stress wave within rock specimen can be enhanced, which greatly inhibits the deterioration of rock material and induce the higher coupled strength. The greater static pre-stress can strengthen the enhancing effect of stress wave propagation, corresponding with the higher coupled strength. However, once the axial static pre-stress exceeds the limitation of elastic deformation of the rock materials, the closed micro discontinues open again, and more cracks generate within rock specimen, which provides abundant reflected surfaces for the stress wave, resulting in the decrease of the coupled strength. The 0.2, 0.4, and 0.6 UCS test are effectively testing the undeformed rock while some amount of inelastic deformation has accrued during the application of the 0.8 UCS which has weakened the rock. This complexity partly might results from the inherently strain hardening nature of microcrack growth under compressive loads before cracks reach a critical length41,42. In other words, the lower UCS tests did not achieve the critical crack length while the higher 0.8 UCS did. Similar behavior was noted in Aben et al.43 and in Braunagel and Griffith44 wherein a certain threshold of pre-existing damage was required in a sample in order to achieve a significant decline in dynamic compressive strength. Therefore, under a given dynamic strain rate, the coupled strength increases with increasing static pre-stress first and then decreases, where the maximum value occurs under the static pre-stress of 60% UCS; the similar results were also observed in the intact rock specimen and rock specimen containing combined fissures and rectangular holes5.Fig. 4 Dynamic stress equilibrium check for the typical cross-fissured specimens (a) (No.0.2- D1) and (b) (No. 0.8- D4).

Fig. 5 Three dimensional stress–strain curves of the cross-fissured rock specimens under coupled static-dynamic loading tests with different pre-stress of (a) 0.2 UCS, (b) 0.4 UCS, (c) 0.6 UCS, and (d) 0.8 UCS.

Fig. 6 (a) Determination of the dynamic strength and coupled dynamic strength; (b) the coupled strengths of the cross-fissured rock specimens under varied coupled static-dynamic load in three-dimensional perspective; (c) and (d) show the effects of dynamic strain rate and static pre-stress on the coupled strength and total strength of the fissured specimens.

Table 2 Coupled strength of the fissured specimens under varied coupled static-dynamic strain rates.

Notation	Static pre-stress (MPa)	Strain rate (s−1)	Coupled strength (MPa)	
0.2-D1	10.4	63.2	89.1	
0.2-D2	10.4	73.4	91.7	
0.2-D3	10.4	77.1	95.4	
0.2-D4	10.4	91.8	98.9	
0.2-D5	10.4	132.9	114.8	
0.4-D1	20.9	27.0	76.6	
0.4-D2	20.9	63.0	96.5	
0.4-D3	20.9	98.8	107.0	
0.4-D4	20.9	103.8	110.5	
0.4-D5	20.9	126.6	117.2	
0.6-D1	31.3	28.9	81.1	
0.6-D2	31.3	71.8	106.4	
0.6-D3	31.3	101.3	112.0	
0.6-D4	31.3	113.5	117.1	
0.6-D5	31.3	120.7	122.7	
0.8-D1	41.8	60.2	93.1	
0.8-D2	41.8	77.9	97.8	
0.8-D3	41.8	100.5	107.3	
0.8-D4	41.8	116.8	112.2	
0.8-D5	41.8	128.8	117.9	

Failure mechanism

The typical failure modes of the fissured sandstones under coupled static-dynamic loads are concluded in two categories, i.e., the similar dynamic strain rate (Fig. 7a–c) and the same static pre-stress (Fig. 7d–f). It can be observed that the failure mode is mainly influenced via the dynamic strain rate ignoring the static pre-stress in this investigation. In addition, based on the DIC technology, the full-field strain and displacement field of the fissured rock specimen are analyzed to further the failure mechanism. In this study, referenced from the categories of displacement fields proposed by Zhang and Wong45, two displacement field types (DF_I and DF_II) are observed, as shown in Fig. 8a. In the region of two displacement trend lines, DF_I corresponds with relative tensile displacement, while DF_II is associated with both a relative tensile displacement and a shear displacement. In addition, in terms of macroscopic view, four cracks are observed on the specimen surface during loading process, as shown in Fig. 8b. Wing crack initiates from the fissure tips and grows along the axial loading direction, while the anti-wing crack propagates opposite to the wing crack. The shear crack generally extends coplanar or oblique with the preexisting fissure. The mixed tensile-shear cracks can be regarded as the connection of tensile and shear cracks.Fig. 7 Failure modes of the cross-fissured specimens under different coupled static-dynamic loads.

Fig. 8 (a) Two microscopic displacement field types and (b) macro crack types observed in this study.

According to the calculated results of DIC technology, all the failure modes can be categorized into two groups including the mixed tensile-shear failure mode and the shear failure mode, shown in Figs. 9 and 10, respectively. Field rocks generally contains myriads of micro discontinues with irregular distribution due to the geotechnical activities. When the rock specimen is experienced the low dynamic strain rate with insufficient incident energy, the generated cracks generally extend and propagate along or oblique to the axial loading direction, featuring the mixed tensile-shear failure mode. In this loading case, the fissured rock specimen generally broken into large fragments. When the dynamic strain rate increases, the compressed stress wave fleetly propagates to the fissure tips, and the cross-fissured rocks absorb more energy in a short time, resulting in the initiation of more shear cracks for releasing exterior input energy and appearing in the shear failure mode. As the dynamic strain rate increases, the number of broken fragments and the evolution of shear cracks gradually increase, while the fragment size gradually decreases. Generally, a specimen with low aspect ratio may bear more confinement, due to the end friction effect on the surfaces between the specimen end and the bar46,47.Fig. 9 Progressive failure process of the cross-fissured specimens with mixed tensile-shear failure subjected to dynamic load.

Fig. 10 Progressive failure process of the cross-fissured specimens with shear failure subjected to dynamic load.

Figure 9a–f shows the evolution of full-field strain and displacement field of the fissured rock specimen 0.6-D1 with mixed tensile-shear failure subjected to coupled static-dynamic load; Fig. 9g show the corresponding displacement field of the ultimate failure and the recovered fragments. At the initial loading stage, the main stress concentration occurs around the tips of the fissure with inclination angle of 75°, inferred from some scattered tiny high strain zones around the fissure tips F3 and F4, as shown in Fig. 9b. Then, two anti-wing cracks (A1 and A2) coalescence with the two fissure tips at the same side, forming a pair of tensile band with butterfly shape, as shown in Fig. 9c. In addition, a shear crack (S1) initiates from the fissure tip F3 and propagate to the transmitted bar-specimen surface, and two shear cracks (S2 and S3) emerge from fissure tips F1 and F2, and grow oblique to the fissure plane, as shown in Fig. 9d and e. After that, another tensile cracks (T1) initiates from the fissure tip F4, and extends to the incident bar-specimen surface, as shown in Fig. 9e. Eventually, the development of these generated cracks forms an X-shaped fracture networks, as shown in Fig. 9f.

The typically progressive failure process of the cross-fissured specimen 0.6-D4 with shear failure subjected to coupled static-dynamic load is shown in Fig. 10a–f, and the Fig. 10g show the corresponding displacement field of the ultimate failure and the recovered fragments. At the initial loading stage, similar to the specimen 0.6-D1, the relative high strain zone also distributes around the preexisting fissure with inclination angle of 75°. However, new cracks begin to initiate under the action of greater stress concentration. With further loading, strain localizations at the fissure tip propagate further, a visible fracture appears at the right side wall of the opening. During dynamic loading stage, micro cracks initiate from the tip of two pre-existing fissures and extend perpendicularly to the inclination angle; they are in the fashion of anti-wing cracks (A1 and A2). At 165.0 μs, four typical shear cracks (S1, S2, S3 and S4) initiate respectively from the fissure tips are formed as V-shaped. Then, all the cracks propagates towards the loading end and penetrate the whole specimen, gradually losing the integrity of the specimen.

Fragmentation characteristics

During coupled static-dynamic loading process, numerous fractures can be generated and propagate within the fissured rocks, resulting in broken fragments with different size distributions, accompanied with energy dissipation. Thus, the fragment size distribution of the broken specimens can reflect the deformation degree and the fracture behaviors of fissured rocks. The dynamic fragmentation characteristics of the fissured sandstone are affected by the static pre-stress and dynamic strain rate. Under low dynamic strain rate, large fragments play a dominated role in the broken specimen, and feature with columns or lamellar structures. When the fissured sandstone subjected to high dynamic impacting, the specimens are broken into small fragment with homogenous size distribution, and the relatively great cone-shape fragment can be frequently encountered in this loading condition. All the fragment size distributions of the fissured specimens after suffering coupled static-dynamic loads are shown in Fig. 11, which plots the cumulative mass percentage of rock fragment in each interval, as well as the mean fragment size for each specimen. It can be observed that the main fragment size ranges from 10 to 40 mm, and the increasing dynamic strain rate decrease the proportion of large fragments under a given static pre-stress, resulting in the smaller mean fragment size. Under a given dynamic strain rate, the smaller mean fragment size of fissured sandstone can be characterized via the higher static pre-stress.Fig. 11 (a) Fragment size distribution and (b) mean fragment size of the cross-fissured rock specimens under different coupled static-dynamic loads.

Figure 12 depicts the influence of static pre-compression and dynamic strain rate on the μ and λ of the fragment size distribution with GEV fitting; the details are tabled in Table 3. The fragment size distributions of fissured sandstone after coupled static-dynamic loading can be well characterized by the GEV distribution. In this manuscript, the location parameter μ and the scale parameter λ are employed to describe the mean fragment size of broken specimen and the range of the fragment size distribution, respectively. Higher dynamic strain rate and the greater static pre-stress can induce the smaller location parameter μ, associated with smaller fragments size distribution. On the whole, the effect of static pre-stress and dynamic strain rate on the scale parameter λ is not significant. Figure 13 depicts the influence of dynamic strain rate and static pre-stress ratio on the fractal dimension of the fissured specimens under coupled static-dynamic loads in three-dimensional and two-dimensional views. In this investigation, the fractal dimension varies from 2.26 to 2.49, with relative deviation of 10.2%. Fractal dimension supplies a different viewpoint from traditional Euclidian geometry, where the dimensions of plane and body are represented via two and three48. The cross-fissured specimen under higher static pre-stress and greater dynamic strain rate can be characterized via the greater fractal dimension, corresponding with the smaller fragment size with more homogeneous distribution.Fig. 12 (a) and (b) depict the location parameter μ of the fragment size distribution of the cross-fissured specimens with GEV fitting in three-dimensional and two-dimensional perspective; (c) and (d) show the scale parameter λ of the fragment size distribution of the cross-fissured specimens with GEV fitting in three-dimensional and two-dimensional perspective.

Table 3 GEV of rock specimens under different coupled static-dynamic loads.

Notation	Static pre-stress (MPa)	Strain rate (s−1)	GEV fitting	Mean fragment size	Fractal dimension	
ζ	μ	λ	
0.2-D1	10.4	63.2	–0.71	21.11	13.51	15.38	2.26	
0.2-D2	10.4	73.4	–	19.68	13.44	14.88	2.28	
0.2-D3	10.4	77.1	–	17.21	13.61	13.88	2.34	
0.2-D4	10.4	91.8	–	16.10	13.52	13.41	2.36	
0.2-D5	10.4	132.9	–	12.03	13.76	11.57	2.43	
0.4-D1	20.9	27.0	–	20.29	13.49	15.14	2.29	
0.4-D2	20.9	63.0	–	17.78	14.17	14.05	2.34	
0.4-D3	20.9	98.8	–	15.01	14.88	12.89	2.38	
0.4-D4	20.9	103.8	–	13.91	14.53	12.44	2.42	
0.4-D5	20.9	126.6	–	9.80	11.62	10.76	2.43	
0.6-D1	31.3	28.9	–	16.24	13.61	13.46	2.34	
0.6-D2	31.3	71.8	–	16.12	14.05	13.38	2.36	
0.6-D3	31.3	101.3	–	12.95	14.19	12.02	2.44	
0.6-D4	31.3	113.5	–	9.96	11.61	10.60	2.44	
0.6-D5	31.3	120.7	–	9.58	11.93	10.48	2.47	
0.8-D1	41.8	60.2	–	15.59	14.32	13.15	2.37	
0.8-D2	41.8	77.9	–	15.02	14.27	12.89	2.40	
0.8-D3	41.8	100.5	–	10.59	12.47	11.34	2.45	
0.8-D4	41.8	116.8	–	8.59	10.80	10.15	2.47	
0.8-D5	41.8	128.8	–	8.37	11.35	9.98	2.49	

Fig. 13 Influence of dynamic strain rate and static pre-stress ratio on the fractal dimension of the fissured specimens under coupled static-dynamic loads in three-dimensional and two-dimensional perspective.

Conclusions

Field rock masses containing preexisting cross fissures at great deep generally experience both the static pre-stress and dynamic disturbances simultaneously. Understanding the mechanical behaviors of these cross-fissured rock specimens is benefit for improving the efficacy of rock breaking. This study systematically investigated the mechanical response and fragment characteristics of cross-fissured sandstones under coupled static-dynamic loads. Based on the experimental and numerical results, the following conclusions can be drawn.Under the same static pre-load, the coupled strength of cross-fissured sandstones shows evident loading-rate dependence with the increase of dynamic strain rate. Whereas, the coupled strength of fissured specimens increases first and then decreases with increasing static pre-stress under a given dynamic strain rate. The highest coupled strength is observed under the static pre-stress of 60% UCS.

According to the calculated results of DIC technology, the failure mode is mainly influenced via the dynamic strain rate ignoring the static pre-stress in this investigation. In this investigation, all the failure modes can be categorized into two groups including the mixed tensile-shear failure mode under low dynamic strain rate and the shear failure mode under high dynamic strain rate. Regardless of the loading condition, the main stress concentration occurs around the tips of the fissure with greater inclination angle.

The dynamic fragmentation characteristics of the fissured sandstone are affected by the static pre-stress and dynamic strain rate. The higher static pre-stress and greater dynamic strain rate can induce the smaller mean fragment size, the smaller location parameter μ in GEV fitting, and the greater fractal dimension of fissured sandstones, corresponding with smaller fragment size with more homogeneous distribution. Whereas, the effect of static pre-stress and dynamic strain rate on the scale parameter λ is not significant.

Acknowledgements

This research was supported by the Natural Science Foundation of Sichuan Province, China (Grant: No. 2024NSFSC0983), Open Research Subject of Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education, Chongqing Jiaotong University (SLK2023A02), Open Fund of Sichuan Engineering Research Center for Mechanical Properties and Engineering Technology of Unsaturated Soils (No. SC-FBHT2024-06), and the National Natural Science Foundation of China (No. 42407239).

Author contributions

P.F.: conceptualization; formal analysis; investigation; methodology; resources; supervision; visualization; roles/writing—original draft. P.C.: supervision; investigation; validation; visualization, roles/writing—original draft. J.L.: resources; investigation, validation; visualization, roles/writing—review and editing. R.T., H.L.: investigation, supervision, roles/writing—review and editing. All authors read and approved the final manuscript.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Ortlepp WD The behaviour of tunnels at great depth under large static and dynamic pressures Tunn. Undergr. Space. Technol. 2001 16 41 48 10.1016/S0886-7798(01)00029-3
Ortlepp, W. D. The behaviour of tunnels at great depth under large static and dynamic pressures. Tunn. Undergr. Space. Technol. 16, 41–48 (2001).10.1016/S0886-7798(01)00029-3
2. Li A Dai F Liu Y Du H Jiang R Dynamic stability evaluation of underground cavern sidewalls against flexural toppling considering excavation-induced damage Tunn. Undergr. Space Technol. 2021 112 103903 10.1016/j.tust.2021.103903
Li, A., Dai, F., Liu, Y., Du, H. & Jiang, R. Dynamic stability evaluation of underground cavern sidewalls against flexural toppling considering excavation-induced damage. Tunn. Undergr. Space Technol. 112, 103903 (2021).10.1016/j.tust.2021.103903
3. Jain M Hefny AM Chen S Rock dynamics research related to cavern development for ammunition storage Tunn. Undergr. Space. Technol. 2002 14 513 526
Jain, M. et al. Rock dynamics research related to cavern development for ammunition storage. Tunn. Undergr. Space. Technol. 14, 513–526 (2002).
4. Feng P Wei M Dai F Tang R Qiu H Gong J Investigation on the mechanical behaviors of flawed specimens subjected to coupled static-dynamic loads Int. J. Soil. Dyn. Earthq. Eng. 2020 135 106 220 10.1016/j.soildyn.2020.106220
Feng, P. et al. Investigation on the mechanical behaviors of flawed specimens subjected to coupled static-dynamic loads. Int. J. Soil. Dyn. Earthq. Eng. 135, 106–220 (2020).10.1016/j.soildyn.2020.106220
5. Zhou Z Cai X Li X Cao W Du X Dynamic response and energy evolution of sandstone under coupled static-dynamic compression: Insights from experimental study into deep rock engineering applications Rock Mech. Rock Eng. 2020 53 1305 1331 10.1007/s00603-019-01980-9
Zhou, Z., Cai, X., Li, X., Cao, W. & Du, X. Dynamic response and energy evolution of sandstone under coupled static-dynamic compression: Insights from experimental study into deep rock engineering applications. Rock Mech. Rock Eng. 53, 1305–1331 (2020).10.1007/s00603-019-01980-9
6. Feng P Zhao JC Dai F Wei MD Liu B Mechanical behaviors of conjugate-flawed rocks subjected to coupled static-dynamic compression Acta Geotech. 2022 17 1765 1784 10.1007/s11440-021-01322-6
Feng, P., Zhao, J. C., Dai, F., Wei, M. D. & Liu, B. Mechanical behaviors of conjugate-flawed rocks subjected to coupled static-dynamic compression. Acta Geotech. 17, 1765–1784 (2022).10.1007/s11440-021-01322-6
7. Liu X Liu Q Wei L Huang X Improved strength criterion and numerical manifold method for fracture initiation and propagation Int. J. Geomech. 2017 10.1061/(ASCE)GM.1943-5622.0000676
Liu, X., Liu, Q., Wei, L. & Huang, X. Improved strength criterion and numerical manifold method for fracture initiation and propagation. Int. J. Geomech.10.1061/(ASCE)GM.1943-5622.0000676 (2017).10.1061/(ASCE)GM.1943-5622.0000676
8. Liu X Liu Q Liu B Zhu Y Zhang P Failure behavior for rocklike material with cross crack under biaxial compression J. Mater. Civil Eng. 2019 10.1061/(ASCE)MT.1943-5533.0002540
Liu, X., Liu, Q., Liu, B., Zhu, Y. & Zhang, P. Failure behavior for rocklike material with cross crack under biaxial compression. J. Mater. Civil Eng.10.1061/(ASCE)MT.1943-5533.0002540 (2019).10.1061/(ASCE)MT.1943-5533.0002540
9. Zhang B Li S Xia K Yang X Zhang D Wang S Reinforcement of rock mass with cross-flaws using rock bolt Tunnel. Undergr. Space Technol. 2016 51 346 353 10.1016/j.tust.2015.10.007
Zhang, B. et al. Reinforcement of rock mass with cross-flaws using rock bolt. Tunnel. Undergr. Space Technol. 51, 346–353 (2016).10.1016/j.tust.2015.10.007
10. Zhang B Li S Yang X Xia K J. Liu Y Guo S Wang S The coalescence and strength of rock-like materials containing two aligned X-type flaws under uniaxial compression Geomech. Eng. 2019 17 47 56
Zhang, B. et al. The coalescence and strength of rock-like materials containing two aligned X-type flaws under uniaxial compression. Geomech. Eng. 17, 47–56 (2019).
11. Zhang B Li Y Yang X Li S Liu B Pei Y Influence of two cross-flaws geometry on the strength and crack coalescence of rock-like material specimens under uniaxial compression Int. J. Geomech. 2020 10.1061/(ASCE)GM.1943-5622.0001757
Zhang, B. et al. Influence of two cross-flaws geometry on the strength and crack coalescence of rock-like material specimens under uniaxial compression. Int. J. Geomech.10.1061/(ASCE)GM.1943-5622.0001757 (2020).10.1061/(ASCE)GM.1943-5622.0001757
12. Zhou X Zhang J Wong L Experimental study on the growth, coalescence and wrapping behaviors of 3D cross-embedded flaws under uniaxial compression Rock Mech. Rock Eng. 2018 51 1379 1400 10.1007/s00603-018-1406-4
Zhou, X., Zhang, J. & Wong, L. Experimental study on the growth, coalescence and wrapping behaviors of 3D cross-embedded flaws under uniaxial compression. Rock Mech. Rock Eng. 51, 1379–1400 (2018).10.1007/s00603-018-1406-4
13. Liang D Zhang N Rong H Xiang Z Experimental and numerical studies on crack initiation and coalescence in sandy mudstone with prefabricated cross-flaws under uniaxial compression Shock Vib. 2021 2021 1 17
Liang, D., Zhang, N., Rong, H. & Xiang, Z. Experimental and numerical studies on crack initiation and coalescence in sandy mudstone with prefabricated cross-flaws under uniaxial compression. Shock Vib. 2021, 1–17 (2021).
14. Zhou X Li L Berto F Cracking behaviors of rock-like specimens containing two sets of preexisting cross flaws under uniaxial compression J. Test. Eval. 2019 47 838 867 10.1520/JTE20170358
Zhou, X., Li, L. & Berto, F. Cracking behaviors of rock-like specimens containing two sets of preexisting cross flaws under uniaxial compression. J. Test. Eval. 47, 838–867 (2019).10.1520/JTE20170358
15. Cao R Cao P Fan X Xiong X Lin H An experimental and numerical study on mechanical behavior of ubiquitous-joint brittle rock-like specimens under uniaxial compression Rock Mech. Rock Eng. 2016 49 4319 4338 10.1007/s00603-016-1029-6
Cao, R., Cao, P., Fan, X., Xiong, X. & Lin, H. An experimental and numerical study on mechanical behavior of ubiquitous-joint brittle rock-like specimens under uniaxial compression. Rock Mech. Rock Eng. 49, 4319–4338 (2016).10.1007/s00603-016-1029-6
16. Feng P Xu Y Dai F Effects of dynamic strain rate on the energy dissipation and fragment characteristics of cross-fissured rocks Int. J. Rock Mech. Min. Sci. 2021 138 104600 10.1016/j.ijrmms.2020.104600
Feng, P., Xu, Y. & Dai, F. Effects of dynamic strain rate on the energy dissipation and fragment characteristics of cross-fissured rocks. Int. J. Rock Mech. Min. Sci. 138, 104600 (2021).10.1016/j.ijrmms.2020.104600
17. Du H Dai F Xu Y Yan Z Wei M Mechanical responses and failure mechanism of hydrostatically pressurized rocks under combined compression-shear impacting Int. J. Mech. Sci. 2020 165 105219 10.1016/j.ijmecsci.2019.105219
Du, H., Dai, F., Xu, Y., Yan, Z. & Wei, M. Mechanical responses and failure mechanism of hydrostatically pressurized rocks under combined compression-shear impacting. Int. J. Mech. Sci. 165, 105219 (2020).10.1016/j.ijmecsci.2019.105219
18. Grady DE Fragment size distributions from the dynamic fragmentation of brittle solids Int. J. Impact Eng. 2008 35 1557 1562 10.1016/j.ijimpeng.2008.07.042
Grady, D. E. Fragment size distributions from the dynamic fragmentation of brittle solids. Int. J. Impact Eng. 35, 1557–1562 (2008).10.1016/j.ijimpeng.2008.07.042
19. Li C Xu Y Chen P Li H Lou P Dynamic mechanical properties and fragment fractal characteristics of fractured coal–rock-like combined bodies in split hopkinson pressure bar tests Nat. Resour. Res. 2020 29 3179 3195 10.1007/s11053-020-09656-w
Li, C., Xu, Y., Chen, P., Li, H. & Lou, P. Dynamic mechanical properties and fragment fractal characteristics of fractured coal–rock-like combined bodies in split hopkinson pressure bar tests. Nat. Resour. Res. 29, 3179–3195 (2020).10.1007/s11053-020-09656-w
20. Grady DE Winfree NA Impact fragmentation of high-velocity compact projectiles on thin plates: A physical and statistical characterization of fragment debris Int. J. Impact Eng. 2001 26 249 262 10.1016/S0734-743X(01)00085-9
Grady, D. E. & Winfree, N. A. Impact fragmentation of high-velocity compact projectiles on thin plates: A physical and statistical characterization of fragment debris. Int. J. Impact Eng. 26, 249–262 (2001).10.1016/S0734-743X(01)00085-9
21. Braunagel MJ Griffith WA Microstructural controls on mixed mode dynamic fracture propagation in crystalline and porous granular rocks J. Geophys. Res-sol. Ea. 2023 127 22528
Braunagel, M. J. & Griffith, W. A. Microstructural controls on mixed mode dynamic fracture propagation in crystalline and porous granular rocks. J. Geophys. Res-sol. Ea. 127, 22528 (2023).
22. Rosin P Rammler E The laws governing the fineness of powdered coal J. Inst. Fuel. 1933 7 29 36
Rosin, P. & Rammler, E. The laws governing the fineness of powdered coal. J. Inst. Fuel. 7, 29–36 (1933).
23. Macías-García A Cuerda-Correa EM Díaz-Díez MA Application of the Rosin-Rammler and Gates-Gaudin-Schuhmann models to the particle size distribution analysis of agglomerated cork Mater. Charact. 2004 52 159 164 10.1016/j.matchar.2004.04.007
Macías-García, A., Cuerda-Correa, E. M. & Díaz-Díez, M. A. Application of the Rosin-Rammler and Gates-Gaudin-Schuhmann models to the particle size distribution analysis of agglomerated cork. Mater. Charact. 52, 159–164 (2004).10.1016/j.matchar.2004.04.007
24. Weibull, W. A. Statistical theory of the strength of materials. Ingeniors Vetenskaps Akademien. Proceedings of Royal Swedish InstituteEngineeringResearching. Stockholm:[sn], Vol. 1, 50. (1939).
25. Du H Feng D Yuan X Yi L Xu H Numerical investigation on the dynamic strength and failure behavior of rocks under hydrostatic confinement in shpb testing Int. J. Rock Mech. Min. Sci. 2018 108 43 57 10.1016/j.ijrmms.2018.05.008
Du, H., Feng, D., Yuan, X., Yi, L. & Xu, H. Numerical investigation on the dynamic strength and failure behavior of rocks under hydrostatic confinement in shpb testing. Int. J. Rock Mech. Min. Sci. 108, 43–57 (2018).10.1016/j.ijrmms.2018.05.008
26. Carpinteri A Lacidogna G Pugno N Scaling of energy dissipation in crushing and fragmentation: A fractal and statistical analysis based on particle size distribution Int. J. Fract. 2004 129 131 139 10.1023/B:FRAC.0000045713.22994.f2
Carpinteri, A., Lacidogna, G. & Pugno, N. Scaling of energy dissipation in crushing and fragmentation: A fractal and statistical analysis based on particle size distribution. Int. J. Fract. 129, 131–139 (2004).10.1023/B:FRAC.0000045713.22994.f2
27. Ma Q Cao Z Experimental study on fractal characteristics and energy dissipation of stabilized soil based on shpb test J. Mater. Civ. Eng. 2019 10.1061/(ASCE)MT.1943-5533.0002916
Ma, Q. & Cao, Z. Experimental study on fractal characteristics and energy dissipation of stabilized soil based on shpb test. J. Mater. Civ. Eng.10.1061/(ASCE)MT.1943-5533.0002916 (2019).10.1061/(ASCE)MT.1943-5533.0002916
28. Li XB Zhou ZL Lok TS Hong L Yin TB Innovative testing technique of rock subjected to coupled static and dynamic loads Int. J. Rock Mech. Min. Sci. 2008 45 739 748 10.1016/j.ijrmms.2007.08.013
Li, X. B., Zhou, Z. L., Lok, T. S., Hong, L. & Yin, T. B. Innovative testing technique of rock subjected to coupled static and dynamic loads. Int. J. Rock Mech. Min. Sci. 45, 739–748 (2008).10.1016/j.ijrmms.2007.08.013
29. Frew DJ Akers SA Chen W Green ML Development of a dynamic triaxial Kolsky bar Meas. Sci. Technol. 2010 21 105704 10.1088/0957-0233/21/10/105704
Frew, D. J., Akers, S. A., Chen, W. & Green, M. L. Development of a dynamic triaxial Kolsky bar. Meas. Sci. Technol. 21, 105704 (2010).10.1088/0957-0233/21/10/105704
30. Xia K Yao W Dynamic rock tests using split Hopkinson (Kolsky) bar system—A review J. Rock Mech. Geotech. Eng. 2015 7 27 59 10.1016/j.jrmge.2014.07.008
Xia, K. & Yao, W. Dynamic rock tests using split Hopkinson (Kolsky) bar system—A review. J. Rock Mech. Geotech. Eng. 7, 27–59 (2015).10.1016/j.jrmge.2014.07.008
31. Zhang Z Kou S Jiang L Lindqvist PA Effects of loading rate on rock fracture: fracture characteristics and energy partitioning Int. J. Rock Mech. Min. Sci. 2000 37 745 762 10.1016/S1365-1609(00)00008-3
Zhang, Z., Kou, S., Jiang, L. & Lindqvist, P. A. Effects of loading rate on rock fracture: fracture characteristics and energy partitioning. Int. J. Rock Mech. Min. Sci. 37, 745–762 (2000).10.1016/S1365-1609(00)00008-3
32. Dai F Huang S Xia K Tan Z Some fundamental issues in dynamic compression and tension tests of rocks using split Hopkinson pressure bar Rock Mech. Rock Eng. 2010 43 657 666 10.1007/s00603-010-0091-8
Dai, F., Huang, S., Xia, K. & Tan, Z. Some fundamental issues in dynamic compression and tension tests of rocks using split Hopkinson pressure bar. Rock Mech. Rock Eng. 43, 657–666 (2010).10.1007/s00603-010-0091-8
33. Hou T Xu Q Yang X Lu P Zhou J Experimental study of the fragmentation characteristics of brittle rocks by the effect of a freefall round hammer Int. J. Fract. 2015 194 169 185 10.1007/s10704-015-0046-x
Hou, T., Xu, Q., Yang, X., Lu, P. & Zhou, J. Experimental study of the fragmentation characteristics of brittle rocks by the effect of a freefall round hammer. Int. J. Fract. 194, 169–185 (2015).10.1007/s10704-015-0046-x
34. Shen W Zhao T Crosta GB Dai F Analysis of impact-induced rock fragmentation using a discrete element approach Int. J. Rock Mech. Min. Sci. 2017 98 33 38 10.1016/j.ijrmms.2017.07.014
Shen, W., Zhao, T., Crosta, G. B. & Dai, F. Analysis of impact-induced rock fragmentation using a discrete element approach. Int. J. Rock Mech. Min. Sci. 98, 33–38 (2017).10.1016/j.ijrmms.2017.07.014
35. Mandelbrot BB Wheeler JA The fractal geometry of nature Am. Math. Mon. 1984 91 594
Mandelbrot, B. B. & Wheeler, J. A. The fractal geometry of nature. Am. Math. Mon. 91, 594 (1984).
36. Wang F Yang K You J Lei X Analysis of pore size distribution and fractal dimension in tight sandstone with mercury intrusion porosimetry Results Phys. 2019 13 102283 10.1016/j.rinp.2019.102283
Wang, F., Yang, K., You, J. & Lei, X. Analysis of pore size distribution and fractal dimension in tight sandstone with mercury intrusion porosimetry. Results Phys. 13, 102283 (2019).10.1016/j.rinp.2019.102283
37. Xu X He M Zhu C Lin Y Cao C A new calculation model of blasting damage degree-Based on fractal and tie rod damage theory Eng. Fract. Mech. 2019 220 106619 10.1016/j.engfracmech.2019.106619
Xu, X., He, M., Zhu, C., Lin, Y. & Cao, C. A new calculation model of blasting damage degree-Based on fractal and tie rod damage theory. Eng. Fract. Mech. 220, 106619 (2019).10.1016/j.engfracmech.2019.106619
38. Deng Y Chen M Jin Y Zou D Theoretical analysis and experimental research on the energy dissipation of rock crushing based on fractal theory J. Nat. Gas Sci. Eng. 2016 33 231 239 10.1016/j.jngse.2016.05.020
Deng, Y., Chen, M., Jin, Y. & Zou, D. Theoretical analysis and experimental research on the energy dissipation of rock crushing based on fractal theory. J. Nat. Gas Sci. Eng. 33, 231–239 (2016).10.1016/j.jngse.2016.05.020
39. Pan B Asundi A Xie H Gao J Digital image correlation using iterative least squares and pointwise least squares for displacement field and strain field measurements Opt. Laser Eng. 2009 47 65 74 10.1016/j.optlaseng.2008.10.014
Pan, B., Asundi, A., Xie, H. & Gao, J. Digital image correlation using iterative least squares and pointwise least squares for displacement field and strain field measurements. Opt. Laser Eng. 47, 65–74 (2009).10.1016/j.optlaseng.2008.10.014
40. Zhu Q Li D Han Z Xiao P Li B Failure characteristics of brittle rock containing two rectangular holes under uniaxial compression and coupled static-dynamic loads Acta Geotech. 2021 10.1007/s11440-021-01196-8
Zhu, Q., Li, D., Han, Z., Xiao, P. & Li, B. Failure characteristics of brittle rock containing two rectangular holes under uniaxial compression and coupled static-dynamic loads. Acta Geotech.10.1007/s11440-021-01196-8 (2021).10.1007/s11440-021-01196-8
41. Ashby MF Sammis CG The damage mechanics of brittle solids in compression Pure Appl. Geophys. 1990 133 489 521 10.1007/BF00878002
Ashby, M. F. & Sammis, C. G. The damage mechanics of brittle solids in compression. Pure Appl. Geophys. 133, 489–521 (1990).10.1007/BF00878002
42. Bhat HS Rosakis AJ Sammis CG A micromechanical based constitutive model for brittle failure at high strain rates J. Appl. Mech. 2012 79 031016 10.1115/1.4005897
Bhat, H. S., Rosakis, A. J. & Sammis, C. G. A micromechanical based constitutive model for brittle failure at high strain rates. J. Appl. Mech. 79, 031016 (2012).10.1115/1.4005897
43. Aben FM Doan ML Mitchell TM Toussaint R Reuschle T Fondreist M Gratier J-P Renard F Dynamic fracturing by successive coseismic loadings leads to pulverization in active fault zones J. Geophys. Res-sol. Ea. 2016 121 2338 2360 10.1002/2015JB012542
Aben, F. M. et al. Dynamic fracturing by successive coseismic loadings leads to pulverization in active fault zones. J. Geophys. Res-sol. Ea. 121, 2338–2360 (2016).10.1002/2015JB012542
44. Braunagel MJ Griffith WA The effect of dynamic stress cycling on the compressive strength of rocks Geophys. Res. Lett. 2019 46 6479 6486 10.1029/2019GL082723
Braunagel, M. J. & Griffith, W. A. The effect of dynamic stress cycling on the compressive strength of rocks. Geophys. Res. Lett. 46, 6479–6486 (2019).10.1029/2019GL082723
45. Zhang X Wong L Crack initiation, propagation and coalescence in rock-like material containing two flaws: A numerical study based on bonded-particle model approach Rock Mech. Rock Eng. 2012 46 1001 1021 10.1007/s00603-012-0323-1
Zhang, X. & Wong, L. Crack initiation, propagation and coalescence in rock-like material containing two flaws: A numerical study based on bonded-particle model approach. Rock Mech. Rock Eng. 46, 1001–1021 (2012).10.1007/s00603-012-0323-1
46. Guneyli H Rusen T Effect of length-to-diameter ratio on the unconfined compressive strength of cohesive soil specimens Bull. Eng. Geol. Environ. 2016 75 793 806 10.1007/s10064-015-0835-5
Guneyli, H. & Rusen, T. Effect of length-to-diameter ratio on the unconfined compressive strength of cohesive soil specimens. Bull. Eng. Geol. Environ. 75, 793–806 (2016).10.1007/s10064-015-0835-5
47. Munoz H Taheri A Specimen aspect ratio and progressive feldstrain pattern development of sandstone under uniaxial compression by 3D Digital Image Correlation J. Rock Mech. Geotech. Eng. 2017 9 599 610 10.1016/j.jrmge.2017.01.005
Munoz, H. & Taheri, A. Specimen aspect ratio and progressive feldstrain pattern development of sandstone under uniaxial compression by 3D Digital Image Correlation. J. Rock Mech. Geotech. Eng. 9, 599–610 (2017).10.1016/j.jrmge.2017.01.005
48. Zhang R Dai F Gao MZ Xu NW Zhang CP Fractal analysis of acoustic emission during uniaxial and triaxial loading of rock Int. J. Rock Mech. Min. Sci. 2015 79 241 249 10.1016/j.ijrmms.2015.08.020
Zhang, R., Dai, F., Gao, M. Z., Xu, N. W. & Zhang, C. P. Fractal analysis of acoustic emission during uniaxial and triaxial loading of rock. Int. J. Rock Mech. Min. Sci. 79, 241–249 (2015).10.1016/j.ijrmms.2015.08.020
