Owned by: China Association for Science and Technology
Sponsored by: China Coal Society
Published by: Springer Nature
About issue

The International Journal of Coal Science & Technology is a peer-reviewed open access journal. It focuses on key topics of coal scientific research and mining development, serving as a forum for scientists to present research findings and discuss challenging issues.


Coverage includes original research articles, new developments, case studies and critical reviews in all aspects of scientific and engineering research on coal, coal utilizations and coal mining. Among the broad topics receiving attention are coal geology, geochemistry, geophysics, mineralogy, and petrology; coal mining theory, technology and engineering; coal processing, utilization and conversion; coal mining environment and reclamation and related aspects.


The International Journal of Coal Science & Technology is published with China Coal Society, who also cover the publication costs so authors do not need to pay an article-processing charge.


The journal operates a single-blind peer-review system, where the reviewers are aware of the names and affiliations of the authors, but the reviewer reports provided to authors are anonymous.


  • A forum for new research findings, case studies and discussion of important challenges in coal science and mining development

  • Offers an international perspective on coal geology, coal mining, technology and engineering, coal processing, utilization and conversion, coal mining environment and reclamation and more

  • Published with the China Coal Society

Show More
Editors-in-Chief
Suping Peng, Shimin Liu
Managing Editor
Wanjie Wang
Associate Editors
Bo Hyun Kim, Dongjie Xue, Pedram Roghanchi, Wu Xiao, Zhiqiang Wu
Publishing model
Open Access. Learn about publishing OA with us
Home > Volumes and issues > Volume 10, issue 6

Mechanical responses and acoustic emission behaviors of coal under compressive differential cyclic loading (DCL): a numerical study via 3D heterogeneous particle model

Research Article

Open Access

Published: 09 June 2023

0 Accesses

International Journal of Coal Science & Technology Volume 10, article number 31, (2023)

Abstract

The stability of coal walls (pillars) can be seriously undermined by diverse in-situ dynamic disturbances. Based on a 3D particle model, this work strives to numerically replicate the major mechanical responses and acoustic emission (AE) behaviors of coal samples under multi-stage compressive cyclic loading with different loading and unloading rates, which is termed differential cyclic loading (DCL). A Weibull-distribution-based model with heterogeneous bond strengths is constructed by both considering the stress–strain relations and AE parameters. Six previously loaded samples were respectively grouped to indicate two DCL regimes, the damage mechanisms for the two groups are explicitly characterized via the time-stress-dependent variation of bond size multiplier, and it is found the two regimes correlate with distinct damage patterns, which involves the competition between stiffness hardening and softening. The numerical b-value is calculated based on the magnitudes of AE energy, the results show that both stress level and bond radius multiplier can impact the numerical b-value. The proposed numerical model succeeds in replicating the stress–strain relations of lab data as well as the elastic-after effect in DCL tests. The effect of damping on energy dissipation and phase shift in numerical model is summarized.

1.Introduction

Coal is a typical sedimentary deposit principally composed of carbon and hydrocarbon that is easily combustible. So far, coal is still the most plentiful fossil fuel in the world. The extraction of coal via underground mining (opposite to surface mining) currently accounts for an obviously larger share of global coal production than open pit. For the past several decades, China has dominated the coal mining and steadily maintained its position as the largest coal-production country, and accounted for approximately half of the global coal production in 2021 (Garside 2022). Due to the exhaustion and depletion of shallow coal in east China, more and more mega mines with thick coal seams (Wang et al. 2021) are constructed and shifted to the central and western parts of China and the coal production in these regions now can account for over 70% of total domestic coal production (Chi et al. 2019). However, the precipitation volume is unexpectedly low in these regions, e.g., in provinces of Shaanxi, Inner Mongolia, and Ningxia, where are almost arid and semiarid climates. The precipitation volume of these regions merely accounts for approximately 4% of China’s total precipitation volume (Chi et al. 2022). The low precipitation volume results in an ultra-fragile natural and ecological environment, the coal mining in these areas can easily incur permanent damage to natural circumstances and exacerbate the environmental vulnerability, such as desertification, surface water loss, vegetation degradation. To tackle such a serious issue, an engineering conception and framework: the coal mine-based underground water reservoir (CMUWR) was proposed by Gu (2015) to utilize the coal pillars (here refer to the pillars which are intentionally left intact to protect against roof collapse, see the coal pillars indicated in Fig. 1a) in the gob as the dams with the connection of reinforced concrete artificial dam (Yao et al. 2020a, b), a schematic diagram for CMUWR is shown in Fig. 1a. It is noticed that the dam composed of coal pillars is subject to the load of overlying strata, which can be either static or dynamic: e.g., strata weighting, seismic-events, drilling, and mining operations in adjunct mining stopes. With development of seismic technology, the monitoring of seismicity at the in-situ scale is critical to improving mining efficiency and operation safety (Cai et al. 2020). Unravelling the fundamental pattern of seismic waves is a prerequisite for data analysis and numerical model construction. A real seismic wave in the form of vibration velocity vs. time is shown in Fig. 1d. Two areas (A and B) are zoomed-in as shown in Fig. 1e, f, where the distinct rates (slopes) can be noticed at the rising and falling phases. This represents diverse loading and unloading rates and this kind of cyclic stress is termed “Differential Cyclic Loading (DCL)” in our previous research (Song et al. 2022c, b). Although the coal’s response under compressive cyclic loading has a long record of research (Li et al. 2019b; Zhang et al. 2020; Duan et al. 2021), the attempts to understand the fatigue behaviors by incorporating the diverse loading and unloading rates still have massive unresolved issues and remains unclarified. The coal dams of CMUWR suffer from oscillation and dynamic load due to the effect both natural or human-induced activities: e.g., earthquakes (Brown and Hudson 1974), seismic events (Dang et al. 2018), hydraulic drilling (Zhuang et al. 2019; Yu et al. 2021), shearer cutting (Wang et al. 2015), underground blasting (Yu et al. 2021). The previous studies from Duan et al. (2021) and Nikolenko et al. (2021) have revealed that the compressive cyclic loading can obviously deteriorate the strength and modulus of coal samples, a low failure stress is always observed in cyclic loading. Whether the DCL can exacerbate or mitigate the fatigue damage of coal is still not sufficiently revealed.

Fig. 1
figure 1

a An overview of coal mine-based underground water reservoir (CMUWR), after Zhang et al. (2021a), the mechanical model of the coal dam can be thought of as a coal sample exposed to compressive cyclic loading due to the imposed vertical load from strata. b Passive seismic velocity tomography presented by with the aid of seismic technology, from Cai et al. (2016). c A layout of the sensor station to monitor the seismicity and record the seismic wave, after Cao et al. (2016). d A real specific seismic wave. e A zoomed-in view of section A in d, a rapid loading and slow unloading pattern is noted. d A zoomed-in view of section B in d, a slow loading and rapid unloading pattern is noted

Except for in-situ monitoring and lab test, numerical characterization of coal’s response under DCL is also of fundamental relevance to evaluate the stability of CMUWR especially when the coal dams are frequently subject to dynamic disturbances. The particle based numerical modelling can deliver a more intricate and comprehensive analysis on damage mechanisms and fracturing of porous materials from crack initiation, coalescence, and final rupture (Hazzard and Young 2002b; Potyondy and Cundall 2004; Hoek and Martin 2014). This can facilitate the calibration of model not only from the view of its constitutive relations but also the associated mechanical responses, e.g., AE behaviors and crack patterns. The numerical characterization of the time-stress-dependent damage evolution incurred by constant, monotonic or cyclic loading via particle model has been studied by a plethora of academics (Chang et al. 2002; Hossain et al. 2007; Potyondy 2007; Li et al. 2019a; Song et al. 2019a, 2022a). The recent investigations on numerical cyclic loading are highlighted as follows: Liu et al. (2018) used the particle model numerically to reproduce the results of cyclic-Brazilian test based the flattened Brazilian disks, the development of microcracks is consistent with the observed experimental result, and a macrocrack lies vertically along the central line. The damage in numerical model is based on the degradation of elastic modulus and bond strength with ongoing of cycles. Li et al. (2022a, b) numerically investigated the behavior of mono/hybrid pile foundation under horizontal cyclic stress, the results indicate that the rise of pile stiffness and application of hybrid pile foundation could enhance the secant stiffness and restrict cumulative lateral displacement. Bai et al. (2022a) conducted a particle model-based numerical investigation to reveal the effect of cyclic normal stress on frictional behaviors of jointed rocks, the phase shift incurred by dynamic load was interpreted and was found related to the changes of the local slipping and frozen contacts. Our former work (Song et al. 2019a, b, 2022a, 2023a) developed various numerical models to replicate the mechanical and micro-seismic behaviors of concrete and coal in constant and variable amplitude cyclic loading. The damage is characterized via the stress-induced bond property degradation. These abovementioned articles numerically examined the applicability of PFC in simulating mechanical responses of different materials under various loads. However, these studies are based on identical loading/unloading rate. Up to now, the numerical simulation on DCL tests is still rarely reported due to the lack of robust lab data concerning such a novel loading regime.

This work is motivated by the formerly conducted DCL tests on coal and strives to more in-depth understand its underlying damage mechanisms induced by distinct loading and unloading rates (Song et al. 2022d). An explicit modelling is conducted based on 3D particle model to simulate the mechanical behaviors. Then, from the microscale, a specific attention is paid to how to numerically replicate the stress–strain curves, micro-fracturing, acoustic emission behaviors observed in lab test. The study can facilitate to deepen the understanding of the mechanism how the DCL impacts the coal’s mechanical responses and can provide practical guidance to the stability assessment of coal pillars in underground coal mines and infrastructures.

2.Experimental dataset and brief introduction

In this section, we briefly present the major results of previous laboratory tests (Song et al. 2022d). The compressive DCL tests were conducted based on an advanced closed-loop servo-controlled hydraulic machine and the geometry of cylindrical anthracite coal sample is [50 mm diameter × 100 mm length]. The characterization of section morphology of a coal sample used in lab testing is presented in Fig. 2 in three scales (Song et al. 2022d). Figure 2c shows that the intersectant cracks are densely distributed and the size of coal matrix is close to the order of millimeters. If we characterize the coal matrix and fracture with disc in a 2D view (Fig. 2d), it has a similar contact with the 2D particle model. The particles intersected by cracks are subject to restricted movements in normal and shear directions. The micro-structural characteristics indicate that the coal is suitable to be modelled via discrete method due to densely distributed cracks and matrixes, where normal contact and frictional behaviors can be well characterized via assigned normal and shear stiffness and frictional parameters. To avoid repeating published content, we only present the relevant data and compare with the numerical results in this study. The calibrations in numerical model with the experimental data are mainly from two parts: monotonic loading (see C16, C18 in Fig. 3a) and DCL tests (see Regime 1: C5, C11, C19; Regime 2:C6, C22, C21 in Fig. 3a). The physical properties for these samples are presented in Table 1. Figure 3 shows the distribution of microcracks and microvoids via thin section of anthracite. The monotonic tests provide a reasonable reference UCS, and the stiffness parameters are recommended to be separately calibrated with each dataset of cyclic loaded sample. The AE characteristics (AE count, b-value, magnitudes of AE events), stress–strain relations, fracture evolution are considered and calibrated with both monotonic and DCL tests.

Fig. 2
figure 2

a An overview of section morphology of a coal sample in lab test (Song et al. 2022d). b A zoomed-in view with the plotting scale of 1 mm. c A zoomed-in view with the plotting scale of 200 μm, the cracks are indicted with the red solid lines. d A characterization of coal matrix with particle models and the order of the size is indicated

Fig. 3
figure 3

a An overview of the cylindrical coal samples (Song et al. 2022d), the data of samples used in this study is indicated. b A graph of the micro structure of coal thin section, the microvoids and microcracks are densely distributed

Table 1 Physically properties of coal specimens, data from Song et al. (2022d)

No.

Mass (g)

Length (mm)

Diameter (mm)

Density (kg/m3)

Regimes

C16

254.33

100.13

49.45

1322

Monotonic

C18

247.70

99.75

49.55

1287

Monotonic

C5

250.41

100.13

49.48

1300

DCL-1

C11

260.18

100.15

49.51

1349

DCL-1

C19

251.38

100.17

49.50

1304

DCL-1

C6

246.59

99.95

49.47

1283

DCL-2

C22

249.55

100.34

49.45

1295

DCL-2

C21

249.71

100.10

49.43

1300

DCL-2

3.Results of numerical simulation

The numerical simulation is based on PFC3D, the micro contact model used in this work is the linear parallel bonded model (LPBM). The readers want to understand fundamentals of LPBM can visit the issued PFC official manual (Itasca Consulting Group Inc. 2014) or the article from Potyondy and Cundall (2004).

3.1 Monotonic loading and calibration

The coal samples were subject to wetting–drying (WD) treatments preceding cyclic testing. We choose six samples exposed to two DCL regimes (Regime 1: C5, C11, C19; Regime 2: C6, C22, C21) to calibrate and replicate the mechanical responses. Prior to the DCL simulation, two monotonically loaded samples (C16, C18) were selected to calibrate the quasi-static monotonic strength. As documented in our previous work (Song et al. 2023b), the model with heterogeneous microparameters exhibits a more realistic AE behavior as well as a reasonable development process of microcracks. In the calibration of monotonic loading, we compare the results obtained in homogeneous and heterogeneous models. Due to the lack of AE monitoring in cyclic test (Song et al. 2022d), the calibration of AE parameters is conducted based on a dataset from Zhang et al. (2020), where the coal sample holds a similar monotonic strength with the samples in our work (Song et al. 2022d). The AE data were captured via 8 nano sensors. Only AE events detected by at least four sensors are used for source location. As recommended in Song et al. (2023b), a Weibull distribution with the shape parameter m = 2.5 is implemented in both bond tensile strength (τten) and normal and shear parallel bond stiffness (kn’ and ks’,the ratio of ks’/ kn’ is fixed to 2). The homogeneous model with identical (τten) and (kn’ and ks’) was established as a control group. To assign Weibull distribution based-micro parameters, we used a method introduced by Liu et al. (2004) and Wang et al. (2019) and was also used in our former work (Song et al. 2023b). First of all, the parameter in the homogeneous model, here we call it Q0, should be determined with calibration of lab data. Then the parameter under Weibull distribution is termed Qi. The Q0 and Qi conform to such a relation: Qi = Q0*x, the x is called the random variable. The cumulative distribution function (CDF) and probability density function (PDF) for x are expressed in Eqs. (1) and (2), where xu is threshold and is recommended to be 0 (Vales et al. 1999; Liu et al. 2004; Song et al. 2023a). m is shape parameter, and x0 is an expectation related scale parameter (Wang et al. 2019). The Eqs. (1) and (2) are turned to Eqs. (3) and (4) when xu = 0. The mean E(x) and variance Var(x) are given in Eqs. (5) and (6). The generation of x is realized via the solution of inverse function of which is controlled by m and x0 (Wang et al. 2019). F(x) varies between [0–1], therefore can be denoted by a random number u between [0–1] The random number can be obtained from a universal random number generator. The x and x0 can determined from Eqs. (7) and (8):

$$F\left( x \right) = 1 - \exp \left[ { - \left( {\frac{{x - x_{u} }}{{x_{0} }}} \right)^{m} } \right]$$
(1)
$$f\left( x \right) = \frac{\text{d}F}{{\text{d}x}} = \frac{m}{{x_{0} }}\left( {\frac{{x - x_{u} }}{{x_{0} }}} \right)^{m - 1} \exp \left[ { - \left( {\frac{{x - x_{u} }}{{x_{0} }}} \right)^{m} } \right]$$
(2)
$$F\left( x \right) = 1 - \exp \left[ { - \left( {\frac{x}{{x_{0} }}} \right)^{m} } \right]$$
(3)
$$f\left( x \right) = \frac{\text{d}F}{{\text{d}x}} = \frac{m}{{x_{0} }}\left( {\frac{x}{{x_{0} }}} \right)^{m - 1} \exp \left[ { - \left( {\frac{x}{{x_{0} }}} \right)^{m} } \right]$$
(4)
$$E\left( x \right) = \mathop \smallint \limits_{0}^{\infty } x\,f\left( x \right)\,\text{d}x = x_{0} {\Gamma }\left( {1 + \frac{1}{m}} \right)$$
(5)
$$Var\left( x \right) = E\left( {x^{2} } \right) - (E\left( x \right))^{2} = x_{0}^{2} \left[ {{\Gamma }\left( {1 + \frac{2}{m}} \right) - {\Gamma }^{2} \left( {1 + \frac{1}{m}} \right)} \right]$$
(6)
$$x = x_{0} \left[ { - \ln (1 - u} \right)]^{\frac{1}{m}}$$
(7)
$$x_{0} = \frac{E\left( x \right)}{{{\Gamma }\left( {1 + \frac{1}{m}} \right)}}$$
(8)

Figure 4 details the calibration of monotonic strength. Figure 4a, b respectively show the statistical distributions of τten and kn’ under Weibull distribution (m = 2.5) in numerical model via Python. In Fig. 4c, it is seen that the UCS for C16 and C18 are respectively close to 10 and 16 MPa. To guarantee an identical sample strength, a comprise is made on strength diversity and the reference UCS is set to 16 MPa, the higher strength allows for a more reasonable range for the variation of bond radius multiplier in the following cyclic loading. The homogeneous and heterogeneous models both show a similar elastic modulus and UCS with respect to the results of C18 in Fig. 4c, which validates the rationality of parameters. The corresponding discrete (AE counts recorded in the fixed time steps) and cumulative AE counts in two models are plotted in Fig. 4d, where the first AE event occurs at approximately 5% UCS in heterogeneous model whereas the value is much higher in homogeneous model and almost reaches 60% UCS. Due to the lack of AE monitoring in samples C16 and C18, we adopt an AE dataset (H4) from the same series of coal (Zhang et al. 2020) which has a similar UCS. The stress–strain curves and the AE results are shown in Fig. 4e, f. The first located AE event occurs around 3% UCS which is close to the level in heterogeneous numerical model (5%). The pattern of cumulative AE counts for H4 is highly consistent with the heterogeneous model in Fig. 4c. Figure 4g, h respectively present locations of AE sources in heterogeneous model and lab test. Three levels (small, medium, and large) of AE events are classified according to the magnitudes of AE energy. The amount of AE energy is proportional to the released energy at failure of the bond (Song et al. 2019b, 2023b). The evolutions of three levels of cumulative AE events are respectively indicated by black, red and blue solid lines. The general pattern and proportions of three levels of cumulative AE events in heterogeneous model highly coincide with lab data which further validates the reasonability of heterogeneous model.

Fig. 4
figure 4

Illustration of Weibull distribution in a tensile strength of bond and b parallel bond stiffness in heterogeneous model, the shape parameter of m is set to 2.5. c Stress–strain curves of C18, C16, homogeneous and heterogeneous models in calibration of uniaxial compressive loading. d AE events and cumulative AE events in homogeneous and heterogeneous models. e Stress–strain curves of H4, where the coal samples is from the same series (Zhang et al. 2020). f AE events and cumulative AE events in lab test of H4. g Evolution of three levels of AE events in heterogeneous model. h Evolution of three levels of AE events in lab test of H4

We run 5 sets of numerical models with different random seeds and the stress–strain curves are shown in Fig. 5a. The comparisons between numerical models (18.02 ± 1.59 MPa for UCS; 1.125 ± 0.182 GPa for secant modulus) and experimental results (17.87 MPa for UCS; 1.13 GPa for secant modulus) in terms of secant moduli and UCS are close to each other. The scatters in numerical models indicate that the change of random seeds can cause the moderate variation of the UCS and secant modulus. The random features can directly influence the evolutions of the bond radius multiplier in the following cyclic simulation. Usually, a higher strength is associated with a smaller bond radius multiplier at failure and vice versa. The uncertainties incurred by particle fabric should be considered in simulation, which is also another heterogeneous characteristic in PFC model. In LPBM, a part of stored energy will be released at failure of a bond which shares a similar physical mechanism with a realistic AE event (Hazzard and Young 2002a; Bai et al. 2022b). Some parameters (e.g., stored bond energy (Song et al. 2023b), moment tensor (Bai et al. 2022b)) were used to characterize the magnitudes of AE energy. The Gutenberg-Richter relationship (Gutenberg and Richter 1944) has been used as a basic characterization of the size distribution of earthquake events and rock failure. It is mathematically expressed as below:

$$Log_{10}\left( {{N}} \right) = a - b{{M}}$$
(9)

where N is the number of earthquakes with magnitude larger than M, a and b are fitting parameters. b describes the size distribution scaling, namely, b-value. Due to the nature of M is related to emitted energy. The M in PFC can be characterized by a scaled released energy at bond failure. M in lab tests is normalized via a scaled absolute energy, which is consistent with M calculated in numerical simulation. The scaling factors for normalization in PFC and lab tests are respectively, 1 × 103 and 1 × 10−7. The sample sizes for calculation in PFC and lab test are 400 and 40 respectively. 4 time points are selected to calculate the b-value in numerical model and lab test. Due to different scaling factors, the for b-value may exhibit significant disparity. However, the declining of b-value with the rise of stress in numerical model is consistent with lab data ahead of failure. The details of calculation of b-value of the four timepoints are presented in Fig. 5d, e. The calibrated micro-parameters subject to monotonic loading are plotted in Table 2.

Fig. 5
figure 5

Comparisons of b-value in heterogeneous model and H4 data based on 4 timepoints (1–4): a Stress–strain curves for heterogeneous model under 5 random seeds and lab results for H4. b b-value of the four timepoints in heterogeneous model, the results is based on the case of seed 1. c b-value of the four timepoints in lab test of H4. d The calculation of b-value in heterogeneous model, the magnitude is normalized by a factor of 1e3 based on the AE energy captured in the code. e The calculation of b-value in lab test of H4, the magnitude is normalized by a factor of 1e-7 based on the absolute AE energy captured in lab test

Table 2 Calibrated Weibull microparameters in monotonic loading

Parameter

kn

ks

kn'

ks'

σt

c

ϕ

Porosity

Density

 

(MN/m)

(MN/m)

(Pa/m)

(Pa/m)

(Pa)

(Pa)

(°)

(%)

(g/cm3)

Value

30

15

m = 2.5,

P0 = 1.5×1010

m = 2.5,

P0 = 7.5×109

m = 2.5,

P0 = 2.5×106

3.5×106

26.5

32

1.3

3.2 DCL test and calibration

3.2.1 Stress propagation

The frequency in DCL test is 0.1 Hz, the low stress rate makes the DCL tests fall in the range of quasi-static loading (Zhang and Zhao 2014), where the inertial forces are negligible in this case according to Cerfontaine and Collin (2018). The checking of stress propagation in numerical model is essential to evaluate the accuracy of loading state. The stress-controlled loading regime is similar to the one used in our former study concerning concrete fatigue (Song et al. 2019b, 2023b), where the vertical stress is imposed via the clump composed of densely packed (with overlap) pebbles set on the top surface of assembly and the rotation and horizontal movements are completely restricted. The vertical stress monitoring is conducted at five parts: namely, the input signal, contact stress at two ends of sample, stress in three measuring spheres, the positions of three measuring spheres are shown in Fig. 6c. Figure 6a, b present the stress paths of two modes in previous DCL tests: (1) rapid loading and slow unloading (2) slow loading and rapid unloading (Song et al. 2022d). In Fig. 6c, d, the monitored vertical stresses for C5 and C6 are shown, almost no inertial forces are generated during cyclic loading, which means that the vertical stress recorded at different parts are consistent with the pattern of input stress wave, and meets the requirement of stress propagation under quasi-static loading. This substantiates the accuracy of stress-controlled loading regime constructed in PFC.

Fig. 6
figure 6

a Stress paths for the samples (C5, C11, C19) exposed to DCL of rapid loading and slow unloading, after Song et al. (2022d). b Stress paths for the samples (C6, C22, C21) exposed to DCL of slow loading and rapid unloading, after Song et al. (2022d). c Stress monitoring in simulation in the case (C5) of rapid loading and slow unloading. d Stress monitoring in simulation in the case (C11) of slow loading and rapid unloading, the stresses at five parts (two sample ends, input signal, measuring spheres 1, 2, 3) are monitored in simulation

3.2.2 Evolution of bond radius multiplier and axial strain

The micro-scale damage of particle assembly characterized by the deterioration of bond size, termed “parallel-bonded stress corrosion (PSC) model”, was first proposed by Potyondy (2007) to simulate a stress-time-dependent deformation due to the corrosion reaction of rocks in presence of water. Since then, the following research (Song et al. 2019a, 2022a; Shen et al. 2022; Xu et al. 2022) implemented the time-dependent variation of bond radius multiplier to numerically replicate mechanical responses of rocks under fixed/variable load amplitudes or environmental actions. The original PSC model idealizes the rock as a cemented granular material and then removes the connection based on a stress-level dependent rate. The corrosion rate is expressed as shown in Eq. 10 (Potyondy 2007), where σc is the bond strength, σa is the crack-initiation stress, and β1 and β2 are two constants to determine the corrosion rate. The original PSC model assumes that the rate of bond size reduction remains constant when the load level is fixed and exceeds the σa. Our previous work on fatigue behavior of concrete indicates that the variation of bond size can be in various manners, the multiplier can be even larger than 1 to mimic an enhancement of elasticity. The reduction in bond diameter is specified by the scaled ratio β, the diameter is reduced to βD, both the effective stiffness and the bond force were reduced (Potyondy 2007). The strain rate of the assembly is determined by bond radius, stress level and crack number. The accurate evolution of bond radius multiplier should be determined with calibration of stress–strain relations.

$$v = \left\{ {\begin{array}{*{20}l} {0, \;\;\sigma < \sigma_\text{a} } \hfill \\ {\beta_{1} e^{{\beta_{2} \left( {\frac{\sigma }{{\sigma_{\text{c}} }}} \right)}} ,\;\;\sigma_{\text{a}} < \sigma < \sigma_{\text{c}} } \hfill \\ {\infty ,\;\;\sigma > \sigma_{\text{c}} } \hfill \\ \end{array} } \right.$$
(10)

To reach a balance between computation cost and model complexity, the radii of particles (1.8–2.3 mm) in DCL simulation are slightly larger than the values in monotonic simulation (1.25–2 mm). The strength-related micro parameters in DCL simulation is consistent with the calibrated parameters in monotonic loading. Figure 7 shows the path of bond radius multiplier after calibration with stress–strain relations of the six selected samples. The calibrated numerical stress–strain curves are indicated by orange dash lines and shown in Fig. 8 together with its corresponding lab test results. Here we want to clarify a phase shift is imposed to each sample to better compare the stress–strain curves during cyclic loading, the magnitudes of phase shift are indicted at the top left corner of Fig. 8. In numerical simulation, a user-defined phase shift is usually adopted to offset the cycle gaps when only partial cycles were used to calibrate simulation or to indicate other parameters: e.g., the strain rate and secant modulus (Song et al. 2019b, 2023b). As known from Fig. 8, the stress–strain curves exposed to rapid loading and slow unloading display a more evident “loose” characteristic which is a manifest indication of larger energy dissipation and a higher growth rate of axial strain. This well correlates with the low fatigue life occurred in this regime and implies that the violent energy dissipation and drastic strain growth can occur much earlier in the first regime than the other one. The evolution path of bond radius multiplier in Fig. 7 between two loading regimes are completely distinct. For the first regime (C5, C11, C19), the radius multiplier remains constant or is monotonically reduced with ongoing of cycles. This coincides with the principle of PSC model, which means that the corrosion gradually deteriorates with rise of stress and increase of cycles. The secant moduli are indicated via colorful bubble diagram (Song et al. 2022d), a high consistency is noted at the moment between drastic drop of moduli and bond multiplier reduction for the samples exposed to rapid loading and slow unloading. For the second loading regime (slow loading and rapid unloading), it is noted that the multiplier first exhibits a step rise with elevation of stress, then the multiplier is reduced and the failure is reached. It can be interpreted from the view of two competing mechanisms: stiffness hardening and softening. Normally, the growth rate of axial strain is increased when the stress level is elevated and the bond radius multiplier remains unchanged. As known from the stress–strain results shown in Fig. 8, the growth rate of the axial strain doesn’t increase significantly after the elevation of stress level, which means that the bond radius multiplier should be raised to offset or restrain the effect of stress elevation. For Figs. 8d, e, two phases respectively correlating stiffness hardening and softening are clear, the peak of the secant modulus is the connection point of the two phases, the variation of radius multiplier varies distinctively before and after the point. For the samples in the second regimes, the failure all occurred in monotonic loading phase after 62 cycles. Therefore, a linear reduction of bond radius multiplier is implemented in this phase to characterize an intensive development of damage prior to failure. Figure 9 presents the comparisons of zeroed axial strains (values of the axial strain of the first cycle is zeroed) at maximum (P) and minimum (R) stresses between lab tests and numerical simulation. The results clearly indicate the applied path of bond radius multiplier in numerical simulation delivers an excellent match with lab testing results and the accuracy of the model is validated.

Fig. 7
figure 7

Evolution path of applied bond radius multiplier versus timestep and the evolutions of secant loading moduli (colorful bubble diagram) and stress levels are indicated a C5 b C11 c C19 d C6 e C22 e C21

Fig. 8
figure 8

Comparisons of stress–strain curves between lab test and numerical simulation, a phase shift is added between lab test and simulation results to offset the strain gap, the phase shift is indicated at the left top corner of each figure a C5 b C11 c C12 d C6 e C22 f C21, the lab data are from Song et al. (2022d)

Fig. 9
figure 9

Comparisons of the zeroed axial peak and residual strains between lab test and numerical simulation, for the convenience of comparison, the axial strains were zeroed from the first cycle, the lab data are from Song et al. (2022d)

3.2.3 AE behaviors in DCL tests
3.2.3.1 AE energy evolution

Numerical AE system can reduce the labor intensity, especially when AE source location is used. A 3D real-time AE system is implemented in PFC3D during both monotonic and DCL tests. For DCL tests, Fig. 10 plots the evolution of AE parameters from beginning to failure based on the samples exposed to DCL. Seven representative moments (a-g) are indicated in Fig. 10 which are used thereafter to show the spatial distribution of AE event magnitude and to calculate the numerical b-value. The AE evolutions for the samples under the first DCL regime (rapid loading and slow unloading) obviously exhibit several “active” phases particularly when the cyclic stress is raised or bond radius multiplier is sharply reduced. For the results in the second regime (slow loading and rapid unloading), the bond radius multiplier is raised since the second cyclic loading stage (CLS), the enhanced multiplier evidently boost the active growth of AE energy in the phase with lower stress level. A transient spike of AE energy is noticed in Fig. 10d, this occurs without reduction of bond radius multiplier and is possibly induced by the elevation of localized contact force in the heterogeneous model.

Fig. 10
figure 10

Evolution of AE energy in numerical model a C5 b C11 c C19 d C6 e C22 f C21

3.2.3.2 Spatial distribution of AE events magnitude

We particularly select seven moments (indicated in Fig. 10) correlating different stress levels and cycle numbers to present the spatial distribution of AE energy magnitudes, see Fig. 11. The size of a cubic in Fig. 11 is proportional to the amount of released energy at bond failure. The results clearly demonstrate that the sample under the loading regime with slow loading and rapid unloading holds a significantly larger magnitude of AE energy shortly ahead of rupture than the samples subject to the other regime. This is attributed to the higher stress that the sample experienced under the regime 2, the higher stress will result in more stored energy at contact of two particles. Therefore, the amount of released energy is increased accordingly. Apart from the effect of stress level, the value of bond radius multiplier is another factor to impact the released energy. For regime 2, the stiffness hardening is observed for all 3 samples at earlier phase, the increased bond radius multiplier (larger than 1) is used instead of remaining constant or gradually reduced. The increase in bond radius multiplier can make the assembly survive under a higher stress via enhanced elastic moduli (Potyondy 2007), the stiffened assembly shows a good performance in energy storage and displays a better elasticity.

Fig. 11
figure 11

Spatial distribution of AE events with indication of the magnitudes of scaled AE energy a C5 b C11 c C19 d C6 e C22 f C21

3.2.3.3 Evolution of numerical b-value

We calculated the numerical b-value based on the released energy of AE events of the last five moments shown in Fig. 11, the scaling factor of the released energy is set to 1×103 and is consistent with the value used in monotonic loading in Sect. 3.1. As documented by Bai et al. (2022b), the b-value can be computed based on moment tensor of AE events, the events that fall within a specific space window and time window (selected timestep) are considered to affect the size of event magnitudes and need to be carefully calibrated. In our study, the sample size for calculation of b-value is 400 AE events. Figure 12 shows the relations of normalized magnitude and the distribution probability. Figure 13a summarized the interrelations between the b-value and the stress level as well as the bond radius multiplier in a 3D diagram. Figure 13b, c respectively plot the impact of stress level and bond radius multiplier on b-value. The results in Fig. 13b indicate that the stress increase ahead of failure (particularly for regime 2) will result in a reduction in b-value, the samples in regime 1 obviously holds a much higher b-value due to the moments selected are far from final rupture. The relation between the bond radius multiplier and b-value revealed in Fig. 13c shows that the reduction of bond multiplier can also incur a smaller b-value. To sum up, the elevation of stress level and reduction of bond radius multiplier can both intensify the damage of particle assembly from insights of b-value. The quantification of damage in numerical model should simultaneously consider the effect of both stress and bond radius multiplier.

Fig. 12
figure 12

Numerical b-values of the last five moments shown Fig. 10

Fig. 13
figure 13

Distribution of b-values versus stress level and bond radius multiplier a 3D distribution b b-value versus stress level c b-value versus bond radius multiplier. The solid lines in Fig. 13b, c are just plotted for the indication of the general distribution pattern of b-value and are not fitting lines

4.Discussion and outlook

4.1 Discussion on energy dissipation and damping

As revealed in previous DCL tests, the anthracite samples exposed to rapid loading and slow unloading can highly likely incur a lower fatigue life than the samples subjected to a contrary loading pattern (Song et al. 2022d). In addition, more energy was dissipated when the loading rate is larger than the counterpart of unloading. To examine whether this characteristic is well replicated, we present and compare the evolution of discrete energy pillars under two loading regimes, the calculation of the energy of a cycle is based on the summation of discrete energy pillars. Each energy pillar quantifies the energy (positive value for shortened specimen in loading phase and negative value for lengthened specimen in unloading phase) between two consecutive data points (Song et al. 2022d). The results of 15th cycle of C11 and C22 are shown as the example in Fig. 14. The data of the lab test is shown in Fig. 14a (Song et al. 2022d), it is noticed that a 0.36 s phase shift is seen between the peak of stress and cumulative energy pillars. The phase shift is obvious for C11 in lab test, and the counterpart is marginal in C22 which is exposed to an opposite loading pattern. The results plotted in Fig. 14b is from the numerical simulation in this work, the data for the same cycle (15th) is used to make a comparison. The numerical result for C11 shows a 470 timestep (TS) phase shift at the peak stress, and the magnitude of phase shift for C22 is lower. This coincides with the lab tests documented in Fig. 14. Moreover, the significant phase shift in lab tests between stress and cumulative energy pillars can incur excessive discrete positive energy after reaching the peak stress, which is the root reason for the more dissipated energy observed in the first loading regime at the end of the cycle, see the area III. This specific behavior is well replicated in our simulation, see the area III of Fig. 14b. Even the differences of absolute values of energy dissipation between two regimes in lab tests and numerical simulation are not close to each other, the performed simulation still succeeds in replicating the most dominant characteristics observed in lab tests: e.g., phase shift, excessive discrete positive energy, and difference in energy dissipation between two regimes.

Fig. 14
figure 14

Evolution of energy pillars and phase shift between two DCL regimes, the 15th cycle of C11 and C22 is used as an example a Lab test results for C11 and C22, after Song et al. (2022d) b Numerical simulation results for C11 and C22

The damping in LPBM is belong to the dashpot component which can provide viscous behaviors and transmits a dashpot force. The normal Fnd and shear Fsd dashpot forces are generated via normal and shear critical-damping ratios, βn and βs, see Eqs. (1112). Where the \(\dot{{\delta }_{\text{n}}}\) and \(\dot{{\delta }_{\text{s}}}\) are relative normal and shear translational velocities. For compact assemblies, non-zero damping is recommended to use to reach the equilibrium and to conduct quasi-static deformation simulations (Itasca Consulting Group Inc. 2014). The low frequency DCL tests makes the loading lie in the range of quasi-static state. Here we discussed the effect of normal and shear damping parameters on energy dissipation under the two DCL regimes. We assign four levels of damping ratios (both βn and βs): 0.1, 0.25, 0.5, 0.75 in the first cycle. Figure 15a indicates the results for rapid loading and slow unloading and vice versa for Fig. 15b. It is observed that with increase of damping ratios: the fluctuations of the discrete energy pillars are obviously mitigated in both loading and unloading phases under the two regimes. This indicates that a highly-damped system can smooth the evolution of axial strain and can result in a less scattered discrete energy pillars. Here we count the phase shift between the peak stress and the peak of cumulative energy, it is seen that a greater phase shift is observed when the larger damping ratio is used in both two loading regimes. The effect of damping on dissipated energy is not clearly observed from the values indicated in the zoomed-in view of area II, it seems that the dissipated energy during the first cycle is independent on the damping ratios for the values used in this work.

$$F_{\text{n}}^{d} \left( x \right) = 2\beta_{\text{n}} \sqrt {m_{\text{c}} k_{\text{n}} } \dot{\delta }_{\text{n}}$$
(11)
$$F_{\text{s}}^{d} \left( x \right) = 2\beta_{\text{s}} \sqrt {m_{\text{c}} k_{\text{s}} } \dot{\delta }_{\text{s}}$$
(12)
Fig. 15
figure 15

Evolution of energy pillars and phase shift exposed to different damping ratios (0.1, 0.25, 0.5, 0.75) a Rapid loading and slow unloading b Slow loading and rapid unloading

4.2 Outlooks

The recent numerical simulations on fatigue behaviors of brittle geomaterials mainly focus on the macroscopic mechanical response, such as the replication of hysteresis loops and failure patterns. The characterization and replication of progressive damage in meso-micro-scales under cyclic loading are still rarely considered and reported. To better present the innovative points and major advantages of the work, a comparison is made in a tabular form between our work and previous studies in literatures, see Table 3. Furthermore, we also discuss the limitations of current model and the outlooks and future research directions are highlighted to provide a research framework for readers, see Table 4.

Table 3 Comparisons of previous literature with our work

Previous work in literature

This work

(Qi et al. 2022) Cycle number is limited which is not robust to show a fatigue process

(Liu et al. 2018) No AE characteristic was presented and its correlation with stress level in numerical modelling was not discussed

(Song et al. 2019b) The active AE events at initial loading phase and lower cyclic stress were not fully presented due to the assigned identical bond strength

(Song et al. 2019b) The active AE events at initial loading phase and lower cyclic stress were not fully presented due to the assigned identical bond strength, only constant-amplitude cyclic loading was simulated

(Zhang et al. 2021b) The homogenous bond strength was assigned according to different mineral contents; the heterogeneity was not well characterized

(Chen et al. 2022) An assumption of linear strength degradation model was proposed, its applicability under multi-level cyclic loading was not examined

(Xu et al. 2021) A 2D meso-scale elasto-plastic damage model to simulate the fatigue behavior of rock during cyclic loading was proposed, its applicability under multi-level cyclic loading was not examined

Different stress rates (Differential cyclic loading)were simulated in loading and unloading phases, this is a novel regime

The implementation of the heterogeneous bond strengths has well replicated the AE characteistics at the lower stress level observed in lab testing

The evolution of b-value in our model is consistent with the classic b-value evolution preceding to failure reported in literature and provides a new approach in PFC based-simulation to show the AE features

The behaviors of coal under multi-level cyclic loading were well replicated which can vadidate the applicatbility of the model under both fixed and variable stress-amplitude

The number of cycles is robust and the whole cyclic loading process was simulated

The parameters to determine the b-value and the magnitude of the AE event was revealed

Table 4 Limitations and outlooks

Shortcoming of current model

Outlooks

The completely sphere particles used in 3D LPBM failed to reproduce the realistic interlocking between mineral grains which usually results in an underrated resistance when rotation and sliding occur between mineral grains. The simulated uniaxial compressive strength (UCS)/uniaxial tensile strength (UTS) ratios in the framework of LPBM always deviate from a reasonable value, which can incur an excessive rotation of sphere compared with irregular particles like the actual rock minerals

The explicit characterizations of different mineral contents were not considered in this work, the heterogeneity in stiffness, elastic modulus, friction, cohesion is not modelled. The complex interactions under load between different minerals can cause the specific crack paths and failure patterns which need to be well considered

Inherent cracks or complex fracture networks were not considered. The inherent cracks have significant effects on rock mechanical properties: e.g., the UCS, UTS, Young's modulus and stress–strain response

The algorithm incorporating the effect of current cracks on the stress field of the adjacent areas and properties degradation is not considered, which is not able to well reproduce the fracture process zone, which is nucleated at the crack tips and is often observed in the laboratory test

The new version of the particle flow code has made it possible to use convex polygons and polyhedral. This allows modellling of non-spherical objects without using clumps, which are extremely computational expensive. The interactions between rigid blocks are possible using regular contact models

Mineral characterization and flowing behaviors of particles can be efficiently realized via the image processing, the classification and identification of different minerals can be based on the gray-scale imaging and deep learning (Liao et al. 2021; Li et al. 2022b)

The inherent complex fracture network can be precisely modelled via CT technology (Chen et al. 2020; Zhao et al. 2021), the CT scanning can help to accurately visualize the orientations and scales of the most dominant cracks, beddings and voids, this can help to reproduce the micro-structural characteristics of the samples

The damage characterization in numerical model should considering the interaction of the current microcracks and the stress field. The recent work from Han et al. (2022) has proposed a new approach to characterize the damage after considering the degradation of the property after crack appearance. The bond-based damage model should also cover some specific behaviors, e.g., self-healing for rock salt. Where the damage and healing can be interconvertible

5.Concluding remarks

A Weibull-distribution-based heterogeneous 3D particle model is employed to simulate the mechanical and micro-seismic responses of coal samples exposed to multi-level differential cyclic loading (DCL). A microscale damage model is implemented by considering the time-stress-dependent evolution of radius multiplier, which incorporates the characterization of stiffness hardening and softening. The proposed model strives to replicate the major results obtained in previous DCL test (Song et al. 2022d) and numerically show its AE behaviors after model validation and calibration with the monotonic loading and AE data. The following conclusions can be drawn:

  1. (1)

    The Weibull heterogeneous model exhibits more realistic behaviors in terms of AE count as well as the ratio of stress level to UCS when the first located AE count occurs. The shape parameter mis set to 2.5. The employed numerical model highly replicates the axial strain evolution under the DCL regime. The damage is impacted by bond radius multiplier and stress level.

  2. (2)

    The numerical b-value is calculated via the released bond energy in the numerical model, the b-value is reduced with stress elevation and reduction of bond radius multiplier shortly preceding the failure. Two different mechanisms: stiffness hardening and softening coexist in regime of slow loading and rapid unloading. Two loading regimes have totally different damage mechanisms, this also results in the distinct evolution of AE energy in numerical models.

  3. (3)

    The characteristics of energy dissipation within a cycle under two regimes are simulated, the phase shift and larger energy dissipated in regime 1 is well replicated. The effect of damping on energy dissipation and phase shift between stress and cumulation of energy is discussed. The highly damped system holds a larger phase shift as well as the less scattered energy variation.

References

[1] Bai Q, Konietzky H, Dang W (2022a) Microscopic modeling of frictional response of smooth joint under normal cyclic loading. Rock Mech Rock Eng 55:169–186. https://doi.org/10.1007/s00603-021-02652-3
[2] Bai Q, Zhang C, Paul Young R (2022b) Using true-triaxial stress path to simulate excavation-induced rock damage: a case study. Int J Coal Sci Technol 9:49. https://doi.org/10.1007/s40789-022-00522-z
[3] Brown ET, Hudson JA (1974) Fatigue failure characteristics of some models of jointed rock. Earthq Eng Struct Dyn 2:379–386. https://doi.org/10.1002/eqe.4290020407
[4] Cai W, Dou L, Si G et al (2016) A principal component analysis/fuzzy comprehensive evaluation model for coal burst liability assessment. Int J Rock Mech Min Sci 81:62–69. https://doi.org/10.1016/j.ijrmms.2015.09.028
[5] Cai W, Bai X, Si G et al (2020) A monitoring investigation into rock burst mechanism based on the coupled theory of static and dynamic stresses. Rock Mech Rock Eng 53:5451–5471. https://doi.org/10.1007/s00603-020-02237-6
[6] Cao A, Dou L, Cai W et al (2016) Tomographic imaging of high seismic activities in underground island longwall face. Arab J Geosci 9:232. https://doi.org/10.1007/s12517-015-2087-x
[7] Cerfontaine B, Collin F (2018) Cyclic and fatigue behaviour of rock materials: review, interpretation and research perspectives. Rock Mech Rock Eng 51:391–414. https://doi.org/10.1007/s00603-017-1337-5
[8] Chang SH, Yun KJ, Lee CI (2002) Modeling of fracture and damage in rock by the bonded-particle model. Geosyst Eng 5:113–120. https://doi.org/10.1080/12269328.2002.10541196
[9] Chen M, Yang S-Q, Ranjith PG, Zhang Y-C (2020) Cracking behavior of rock containing non-persistent joints with various joints inclinations. Theor Appl Fract Mech 109:102701. https://doi.org/10.1016/j.tafmec.2020.102701
[10] Chen X, Tang M, Zeng P et al (2022) Modeling progressive damage and failure of rock under continuous cyclic loading considering strength degradation of mesoscopic finite elements. Environ Earth Sci 81:335. https://doi.org/10.1007/s12665-022-10449-y
[11] Chi M, Zhang D, Honglin L et al (2019) Simulation analysis of water resource damage feature and development degree of mining-induced fracture at ecologically fragile mining area. Environ Earth Sci 78:88. https://doi.org/10.1007/s12665-018-8039-5
[12] Chi M, Cao Z, Wu B et al (2022) Prediction of water source and water volume of underground reservoir in coal mine under multiple aquifers. Water Supply 22:4067–4081. https://doi.org/10.2166/ws.2022.035
[13] Dang W, Konietzky H, Chang L, Frühwirt T (2018) Velocity-frequency-amplitude-dependent frictional resistance of planar joints under dynamic normal load (DNL) conditions. Tunn Undergr Sp Technol 79:27–34. https://doi.org/10.1016/j.tust.2018.04.038
[14] Duan M, Jiang C, Yin W et al (2021) Experimental study on mechanical and damage characteristics of coal under true triaxial cyclic disturbance. Eng Geol 295:106445. https://doi.org/10.1016/j.enggeo.2021.106445
[15] Garside M (2022) Leading hard coal producing countries worldwide in 2020. https://www.statista.com/statistics/264775/top-10-countries-based-on-hard-coal-production/
[16] Gu D (2015) Theory framework and technological system of coal mine underground reservoir. J China Coal Soc 40:239–246
[17] Gutenberg B, Richter CF (1944) Frequency of earthquakes in California*. Bull Seismol Soc Am 34:185–188. https://doi.org/10.1785/BSSA0340040185
[18] Han Y, Liu X, Zhu W (2022) Numerical simulation on rock burst prevention by cutting stress-relief slot. China Min Mag 31:101–109
[19] Hazzard J, Young R (2002a) Simulating acoustic emissions in bonded-particle models of rock. Int J Rock Mech Min Sci. https://doi.org/10.1016/s1365-1609(00)00017-4
[20] Hazzard JF, Young RP (2002b) Moment tensors and micromechanical models. Tectonophysics 356:181–197. https://doi.org/10.1016/S0040-1951(02)00384-0
[21] Hoek E, Martin CD (2014) Fracture initiation and propagation in intact rock—a review. J Rock Mech Geotech Eng 6:287–300
[22] Hossain Z, Indraratna B, Darve F, Thakur PK (2007) DEM analysis of angular ballast breakage under cyclic loading. Geomech Geoengin 2:175–181. https://doi.org/10.1080/17486020701474962
[23] Itasca Consulting Group Inc. (2014) PFC 5.0 (Particle Flow Code) documentation. Minneapolis
[24] Li H, Yang C, Ding X et al (2019a) Weibull linear parallel bond model (WLPBM) for simulating micro-mechanical characteristics of heterogeneous rocks. Eng Anal Bound Elem 108:82–94. https://doi.org/10.1016/j.enganabound.2019.07.018
[25] Li Z, He X, Dou L et al (2019b) Investigating the mechanism and prevention of coal mine dynamic disasters by using dynamic cyclic loading tests. Saf Sci 115:215–228. https://doi.org/10.1016/j.ssci.2019.02.011
[26] Li L, Zheng M, Liu X et al (2022a) Numerical analysis of the cyclic loading behavior of monopile and hybrid pile foundation. Comput Geotech 144:104635. https://doi.org/10.1016/j.compgeo.2022.104635
[27] Li M, Wu Z, Weng L et al (2022b) Quantitative relationships between the mineral composition and macro mechanical behaviors of granite under different temperatures: Insights from mesostructure-based DEM investigations. Comput Geotech 148:104838. https://doi.org/10.1016/j.compgeo.2022.104838
[28] Liao Z, Yang Y, Sun C et al (2021) Image-based prediction of granular flow behaviors in a wedge-shaped hopper by combing DEM and deep learning methods. Powder Technol 383:159–166. https://doi.org/10.1016/j.powtec.2021.01.041
[29] Liu HY, Roquete M, Kou SQ, Lindqvist P-A (2004) Characterization of rock heterogeneity and numerical verification. Eng Geol 72:89–119. https://doi.org/10.1016/j.enggeo.2003.06.004
[30] Liu Y, Dai F, Xu N et al (2018) Experimental and numerical investigation on the tensile fatigue properties of rocks using the cyclic flattened Brazilian disc method. Soil Dyn Earthq Eng 105:68–82. https://doi.org/10.1016/j.soildyn.2017.11.025
[31] Nikolenko PV, Epshtein SA, Shkuratnik VL, Anufrenkova PS (2021) Experimental study of coal fracture dynamics under the influence of cyclic freezing–thawing using shear elastic waves. Int J Coal Sci Technol 8:562–574. https://doi.org/10.1007/s40789-020-00352-x
[32] Potyondy DO (2007) Simulating stress corrosion with a bonded-particle model for rock. Int J Rock Mech Min Sci 44:677–691. https://doi.org/10.1016/j.ijrmms.2006.10.002
[33] Potyondy DO, Cundall PA (2004) A bonded-particle model for rock. Int J Rock Mech Min Sci 41:1329–1364. https://doi.org/10.1016/j.ijrmms.2004.09.011
[34] Qi Z, Xiao J, Lu J et al (2022) Research on rock triaxial cyclic loading and unloading characteristics based on PFC3D. MATEC Web Conf 358:1012. https://doi.org/10.1051/matecconf/202235801012
[35] Shen Z, Huang D, Wang G et al (2022) A mesoscale bond model for discrete element modeling of irregular cemented granular materials. Comput Geotech 152:105051. https://doi.org/10.1016/j.compgeo.2022.105051
[36] Song Z, Konietzky H, Herbst M (2019a) Bonded-particle model-based simulation of artificial rock subjected to cyclic loading. Acta Geotech 14:955–971. https://doi.org/10.1007/s11440-018-0723-9
[37] Song Z, Konietzky H, Herbst M (2019b) Three-dimensional particle model based numerical simulation on multi-level compressive cyclic loading of concrete. Constr Build Mater 225:661–677. https://doi.org/10.1016/j.conbuildmat.2019.07.260
[38] Song Z, Konietzky H, Herbst M et al (2022a) Fatigue and micro-seismic behaviors of concrete disks exposed to cyclic Brazilian testing: a numerical investigation based on a 3D particle-based model. Int J Fatigue 155:106629
[39] Song Z, Konietzky H, Wu Y et al (2022b) Mechanical behaviour of medium-grained sandstones exposed to differential cyclic loading (DCL) with distinct loading and unloading rates. J Rock Mech Geotech Eng 14:1849–1871
[40] Song Z, Konietzky H, Wu Y, Cai X (2022c) Mechanical and microseismic characteristics of sandstones subject to moderate low-frequency differential cyclic loading (DCL) followed by monotonic loading up to failure. Acta Geotech. https://doi.org/10.1007/s11440-022-01570-0
[41] Song Z, Wu Y, Konietzky H et al (2022d) Mechanical behaviors of anthracite coal subject to low-cycle compressive differential cyclic loading (DCL) after wetting–drying (WD) treatment: an experimental study. Geomech Geophys Geo-Energy Geo-Resour 8:116. https://doi.org/10.1007/s40948-022-00423-0
[42] Song Z, Yang Z, Zhang M et al (2023a) Effect of heterogeneity on mechanical and micro-seismic behaviors of sandstone subjected to multi-level cyclic loading: a discrete element method investigation. J Rock Mech Geotech Eng. https://doi.org/10.1016/j.jrmge.2022.11.020
[43] Song Z, Yang Z, Zhang M et al (2023b) Numerical characterization of sandstone’s mechanical responses under multi-level compressive differential cyclic loading (DCL): a 3D particle-based numerical investigation. Comput Geotech 154:105109. https://doi.org/10.1016/j.compgeo.2022.105109
[44] Vales F, Rezakhanlou R, Olagnon C (1999) Determination of the fracture mechanical parameters of porous ceramics from microstructure parameters measured by quantitative image analysis. J Mater Sci 34:4081–4088. https://doi.org/10.1023/A:1004632503249
[45] Wang J, Yu B, Kang H et al (2015) Key technologies and equipment for a fully mechanized top-coal caving operation with a large mining height at ultra-thick coal seams. Int J Coal Sci Technol 2:97–161. https://doi.org/10.1007/s40789-015-0071-4
[46] Wang F, Konietzky H, Herbst M (2019) Influence of heterogeneity on thermo-mechanical behaviour of rocks. Comput Geotech 116:103184. https://doi.org/10.1016/j.compgeo.2019.103184
[47] Wang J, Yang S, Wei W et al (2021) Drawing mechanisms for top coal in longwall top coal caving (LTCC): a review of two decades of literature. Int J Coal Sci Technol. https://doi.org/10.1007/s40789-021-00453-1
[48] Xu T, Fu M, Yang S et al (2021) A numerical meso-scale elasto-plastic damage model for modeling the deformation and fracturing of sandstone under cyclic loading. Rock Mech Rock Eng 54:4569–4591. https://doi.org/10.1007/s00603-021-02556-2
[49] Xu B, Liu X, Zhou X et al (2022) Investigations on the macro-meso cumulative damage mechanism of the discontinuities in soft-hard interbedded rock mass under pre-peak cyclic shear loading. Int J Rock Mech Min Sci 158:105184. https://doi.org/10.1016/j.ijrmms.2022.105184
[50] Yao Q, Tang C, Xia Z et al (2020a) Mechanisms of failure in coal samples from underground water reservoir. Eng Geol 267:105494. https://doi.org/10.1016/j.enggeo.2020.105494
[51] Yao Q, Zheng C, Tang C et al (2020b) Experimental investigation of the mechanical failure behavior of coal specimens with water intrusion. Front Earth Sci. https://doi.org/10.3389/feart.2019.00348
[52] Yu J, Zhu Y, Yao W et al (2021) Stress relaxation behaviour of marble under cyclic weak disturbance and confining pressures. Measurement 182:109777. https://doi.org/10.1016/j.measurement.2021.109777
[53] Zhang QB, Zhao J (2014) A review of dynamic experimental techniques and mechanical behaviour of rock materials. Rock Mech Rock Eng 47:1411–1478. https://doi.org/10.1007/s00603-013-0463-y
[54] Zhang M, Dou L, Konietzky H et al (2020) Cyclic fatigue characteristics of strong burst-prone coal: experimental insights from energy dissipation, hysteresis and micro-seismicity. Int J Fatigue 133:105429. https://doi.org/10.1016/j.ijfatigue.2019.105429
[55] Zhang C, Wang F, Bai Q (2021a) Underground space utilization of coalmines in China: a review of underground water reservoir construction. Tunn Undergr Sp Technol 107:103657. https://doi.org/10.1016/j.tust.2020.103657
[56] Zhang J, Guo Q, Cai M et al (2021b) Particle flow simulation of the crack propagation characteristics of granite under cyclic load. Chin J Eng 43:636–646
[57] Zhao Y, Konietzky H, Herbst M (2021) Damage evolution of coal with inclusions under triaxial compression. Rock Mech Rock Eng 54:5319–5336. https://doi.org/10.1007/s00603-021-02436-9
[58] Zhuang L, Kim KY, Jung SG et al (2019) Cyclic hydraulic fracturing of pocheon granite cores and its impact on breakdown pressure, acoustic emission amplitudes and injectivity. Int J Rock Mech Min Sci 122:104065. https://doi.org/10.1016/j.ijrmms.2019.104065

About this article

Cite this article

Song, Z., Wu, Y., Zhang, Y. et al. Mechanical responses and acoustic emission behaviors of coal under compressive differential cyclic loading (DCL): a numerical study via 3D heterogeneous particle model.Int J Coal Sci Technol 10, 31 (2023).
  • Received

    21 December 2022

  • Revised

    20 February 2023

  • Accepted

    04 May 2023

  • DOI

    https://doi.org/10.1007/s40789-023-00589-2

  • Share this article

    Copy to clipboard

For Authors

Explore