Experimental Investigation and Micromechanical Modeling of Hard Rock in Protective Seam Considering Damage–Friction Coupling Effect

: The hard rock in the protective coal seam of the Pingdingshan Mine in China is a typical quasi-brittle material exhibiting complex mechanical characteristics. According to available research on the mechanical property, the inelastic deformation and development of damage are considered related with crack initiation and propagation, which are main causes of the material degradation. In the present study, an original experimental investigation on the rock sample of the Pingdingshan coal mine is ﬁrstly carried out to obtain the basic mechanical responses in a conventional triaxial compression test. Based on the homogenization method and thermodynamic theory, a damage– friction coupled model is proposed to simulate the non-linear mechanical behavior. In the framework of micromechanics, the hard rock in a protective coal seam is viewed as a heterogeneous material composed of a homogeneous solid matrix and a large number of randomly distributed microcracks, leading to a Representative Elementary Volume (REV), i.e., the matrix–cracks system. By the use of the Mori–Tanaka homogenization scheme, the effective elastic properties of cracked material are obtained within the framework of micromechanics. The expression of free energy on the characteristic unitary is derived by homogenization methods and the pairwise thermodynamic forces associated with the inelastic strain and damage variables. The local stress tensor is decomposed to hydrostatic and deviatoric parts, and the effective tangent stiffness tensor is derived by considering both the plastic yield law and a speciﬁc damage criterion. The associated generalized Coulomb friction criterion and damage criterion are introduced to describe the evolution of inelastic strain and damage, respectively. Prepeak and postpeak triaxial response analysis is carried out by coupled damage–friction analysis to obtain analytical expressions for rock strength and to clarify the basic characteristics of the damage resistance function. Finally, by the use of the returning mapping procedure, the proposed damage– friction constitutive model is applied to simulate the deformation of Pingdingshan hard rock in triaxial compression with respect to different conﬁning pressures. It is observed that the numerical results are in good agreement with the experimental data, which can verify the accuracy and show the obvious advantages of the micromechanic-based model.


Introduction
An increasing number of large rock projects such as the Sichuan-tibet Railway tunnel project [1] and Beishan high-level nuclear waste geological disposal underground project [2] are currently being planned and developed worldwide.In order to guarantee the safety and stability of the rock mass engineering, the mechanical properties of the engineering rocks need to be mastered.The high ground stress, high gas and low permeability conditions that exist in protruding mines at a depth of 1000 meters make protrusion accidents a frequent problem [3,4].The rock (as shown in Figure 1) here is drilled and blasted out of the coal mine tunnels and made into standard rock samples.The corresponding mechanical parameters can be obtained by experiment and theoretical simulation.This method can provide reference for the calculation of deformation and stability of the coal mine [5,6].In this paper, the hard rock in the protective coal seam of the Pingdingshan mine is studied experimentally and analytically in order to provided a unified constitutive model for predicting its mechanical behaviors based on the obtained experimental data.The horizontal shafts of the Pingdingshan coal mine are located a thousand meters deep underground.The hard rock is collected from the protective coal seam, which is a typical discontinuous, anisotropic and quasi-brittle material.It has a complex microstructure, with a high number of fractures and microcracks in its internal structure.In recent years, to obtain the basic mechanical properties of various rocks, the conventional triaxial compression tests are accepted as a powerful method by many scholars on this subject [7][8][9].Over the past three decades, significant progress has been made in the modeling of plastic degradation in quasi-brittle materials [10][11][12].Several theoretical frameworks have also been proposed for plastic damage models, including [13][14][15][16][17], only to name a few.For concrete and other related materials, isotropic and anisotropic damage models have been developed with or without the plastic coupling, for example, [18][19][20][21][22][23].For rock materials, certain models [24-27] have been developed.On the other hand, certain discrete plasticity-damage models [28][29][30][31] have been developed with the help of the micro-plane theory and the discrete thermodynamics formulation in order to better understand the consequences of anisotropic distribution of microcracks in brittle materials [32][33][34].Furthermore, certain discrete plasticity-damage models have been developed with the help of the micro-plane theory and the discrete thermodynamics formulation [35][36][37] in order to better understand the consequences of anisotropic distribution of microcracks in brittle materials [38][39][40].Recently, within the framework of micromechanics, elastoplastic damage modeling has been successfully applied on rock materials by considering penny-shaped cracks [41,42].
In this paper, based on the available micromechanic-based modeling method [43][44][45] , a new multiscale constitutive model simulating the mechanical properties of the hard rock in the protective coal seam of the Pingdingshan coal mine is constructed by combining homogenization theory and irreversible thermodynamics.The discussion focuses on strength prediction and parameter determination based on coupled damage-friction analysis, and new damage criteria are proposed based on the characterization of the damage evolution resistance function.The model is illustrated by a returning mapping procedure to simulate the mechanical behavior of the hard rock in conventional triaxial compression tests with respect to different confining pressures.It is emphasized that for the evolution of inelastic strains, the deformation can be still accurately predicted by the constructed model despite the application of the associated flow law.As one of the prominent advantages of the multiscale model, by using the homogenization method to obtain the full expression of the free energy on the characteristics of a unit cell to establish the basic pattern of damage friction coupling, the strengthening of the local stress contains a similar back stress weakening function, and the model parameters are greatly reduced.The physical meaning is clear and provides a convenient process for the parameter calibration and the engineering application, which are also demonstrated in this paper.

Sample Preparation and Testing Producure
The tested samples in this study are obtained from the protective coal seam of the Pingdingshan Mine in China.The type of rock is diorite.The overall rock samples are light gray surrounded by black particles (as shown in Figure 1).Rock test specifications require standard specimens to be cylindrical, with a diameter of 50 mm and a range of 48 to 52 mm allowed.The height is 100 mm and the allowable range is 95 to 105 mm.For rocks with a heterogeneous coarse-grained structure, non-standard samples are allowed, but the ratio of height to diameter should be within 2 to 2.5.After mining by blasting method, the rock is processed into standard samples with a diameter of 50 mm and a height of 100 mm.The non-parallelism error of the two end faces of the rock sample is less than 0.05 mm, the end faces are perpendicular to the axis of the rock sample, and the deviation is less than 0.25 • .No obvious cracks are observed on the outer surface, and a good homogeneity is seen.The average density of the samples is 2.6 g/cm 3 .
The conventional triaxial compression tests were carried out at the Multiscale Multi-Field Coupled Rock Mechanics Laboratory of Hohai University, using a triaxial rheometer manufactured by TOP Industrie France (see Figure 2).The equipment mainly consists of a triaxial pressure chamber, an axial pressure servo pump, a perimeter pressure servo pump and a computer control system, which can realize conventional triaxial compression tests on rocks, with a wide range of application and high measurement accuracy.Pressure control adopts a high-precision electronic control servo high pressure pump, and the measurement accuracy can reach 0.01 MPa; two highly sensitive displacement sensors are used for axial displacement measurement, which can directly output the axial displacement value of the tested sample.The utilized loading technique is strain-prescribed load.The strainprescribed load requires the instrument's indenter to press down on the rock at a rate of 0.02 mm/min.The lateral strain deformation measuring device (the right subfigure of Figure 2) here is directly placed on the rubber sleeve outside the rock sample when in use.During the test, the computer can directly read the corresponding deformation of the measuring device.The measurement range is 20 mm, and the measurement accuracy is 10 −3 mm.The system can produce the pressure in many ways, among which the axial pressure can be carried out by axial displacement control, pressure control, flow control, lateral displacement control and other loading methods.
The rock of the protective coal seam is noted to exhibit the characteristic mechanical properties of brittle solids.According to the linear portion of the curves, the mean Young's modulus E = 29 GPa and the Poisson's ratio ν = 0.08 are calculated.

Theoretical Framework of Elastic Damage Model
The initiation, propagation and connection of microcracks within the solid matrix are the main mechanisms leading to the regression of mechanical property and material failure.Quasi-brittle rocks with microcracks can be regarded as a kind of composite material and can be therefore studied as a Representative Elementary Volume (REV) containing a solid matrix and a large number of penny-shaped microcracks.By experimental observation, the fracture of the rock sample in conventional compression tests is mainly caused by the crack initiation and propagation.Thus, only the reduction of effective properties by cracks is considered and that by pores is ignored.Based on this, the REV is viewed as a heterogeneous system with scattered microcracks and a solid matrix that has been weakened by pores (Figure 4).The following derivations are studied within the framework of homogenization methods. of 60 to 120 MPa, the macroscopic dilatant cataclastic flow with no shear bands was observed, resulting from the interactions between microcracking and dislocation activities.The stress concentration at crack tips is relaxed by dislocation activities within the neighboring grains and the ability of cracks is insufficient to propagate to localized into shear bands.In this scenario, shear-oriented interfacial microcracks nucleate and form the grain-scale faults, leading to macroscopic distributed deformation.In addition, with the increase of confining pressure, the dominant mechanism of microcracking changes from tensile cracks to compressive-shear cracks.

Constitutive modelling and numerical simulation
In view of the operative deformation mechanism related to microcracking, the damaged sandstone can be modeled as a heterogeneous system composed of an elastic matrix and a number of randomly oriented are concerned, the Mori-Tanaka (MT) homogenization method shows close links with the results established in the Linear Elastic Fracture Mechanics [23,32] .According to the Mori-Tanaka homogenization method [5,29] , ℂ p involved in the stored energy takes the following explicit expression [49] where are two constants only dependent on the Poisson's ratio  m of the solid matrix.
Given the strain free energy (2) , one can derive the relation between the macroscopic stress  and the macroscopic strain and the thermodynamic force associated with the plastic strain  p Consider an isotropic rock matrix where the fourth-order elastic stiffness tensor is expressed as: where k m and µ m are the compressive bulk modulus and shear modulus of the matrix, respectively.By introducing the second-order unit tensor δ, the fourth-order tensor operators J and K, they are denoted as The existence of microcracks leads to the discontinuity of the displacement field in rock material.Therefore, the macroscopic strain ε can be decomposed into two components, namely, the elastic strain ε m and the inelastic strain ε c , corresponding to the solid matrix component and the microcrack component, respectively: Then, the relationship between macroscopic stress σ and strain ε is obtained as follows: The degree of shear expansion of the microcrack and the relative slip of the microcrack face are expressed in terms of the scalar β and the second-order tensor T, respectively.Thereby, one has the inelastic strain tensor ε c as the sum of the hydrostatic and deviatoric parts: The construction of a constitutive model in terms of damage mechanics generally includes three steps: firstly, select a suitable damage variable ϕ to describe the damage state of the material.Budiansky and Connell [46] proposed that the damage variable is related to the fracture density, namely, ϕ = Nd 3 (where N is the number of microcracks per unit volume and d is the radius of the coin crack surface).Secondly, an effective elastic tensor or the expression for free energy of the REV is established, and the thermodynamic forces associated with the internal variables are derived; finally, a suitable damage criterion is proposed and the evolution equation for the damage variables is determined [2].
According to Zhu et al. [40], for damaged solids with microcracks, the free energy W of a single unit can be expressed in the following general form: where W represents the elastic free energy of a solid matrix, and the term of the right side of the equation is the free energy stored in a solid matrix caused by inelastic strain related to crack, and C b is the fourth-order back stress modulus.We are taking into account the internal relationship [37] between the Mori-Tanaka homogenization method (MT) and linear elastic fracture mechanics in dealing with crack problems.In the case of isotropic and open cracks, the following effective elastic tensors can be obtained by the application of the MT method: where α 1 , α 2 are constants only related to the Poisson coefficient ν m of the rock matrix, 1−2ν m , and α 2 = 32 45 . Considering the expression of C hom , the system free energy is expressed as W = ε : C hom ε/2.By combining Equation (5), we can obtain: It should be pointed out that according to the research of Zhu et al. [40], Equation ( 5) is also suitable for crack closure.In this case, the energy dissipation mechanism of crack propagation and sliding friction coupling exist in the REV.The analytical relationship between inelastic strain and macroscopic strain caused by microcracks is no longer valid.
The thermaldynamic force σ c associated with the inelastic strain ε c can be determined from the system free energy, that is, the local stress acting on the crack: Furthermore, σ c is decomposed into two portions: deviatoric stress s c and hydrostatic part: s c = K : σ c and p c = trσ c /3 .In order to describe the inelastic strain due to the sliding friction of closed cracks and to capture the compressive shear damage pattern of quasi-brittle rock materials under compression, a Coulomb-type yielding criterion based on local stress is adopted in the present study: where α is the coefficient of friction of the cracked surface of the rock material.The local stress tensor can be decomposed into deviatoric and hydrostatic parts: Combined with the above equations, Equation ( 9) can be rewritten as: where s = K : σ and p = trσ/3.Considering the expression of free energy and the theory of irreversible thermodynamics, the thermaldynamic force associated with the damage variable ϕ, namely the damage driving force F ϕ , can be derived by following equation: Substituting Equations ( 4), ( 5) and (7) into Equation ( 13), the explicit form of damage driving force can be obtained as follows: The damage evolution criterion based on the strain energy release rate can also be derived where R(ϕ) is the current damage evolution resistance force.

Coupled Friction-Damage Effect and the Strength Criterion
As mentioned above, in the process of rock deterioration, there are two fundamental pathways for energy dissipation in compressive stress-dominated loading condition.One is damage evolution caused by the development of microcracks, and the other is friction caused by sliding along the fissure surface accompanied by dilatation (volume expansion) caused by the non-smooth crack surface.Damage and inelastic strain constantly rise upon loading in the damage-friction coupling process.
In this study, the evolution of damage variable ω and inelastic strain ε c are determined by the associated flow rule, and the directions of the evolution are determined by the orthogonalization criteria: where V = s c / || s c || (18) in which V represents the flow direction of the partial portion of inelastic strain ε c .λ ϕ and λ c are damage and plasticity multipliers, respectively.
Comparing the expression of Equation ( 18) with that of Equation ( 4), one has Under standard triaxial loading circumstances, the evolution direction V of inelastic shear strain does not change, so the cumulative damage parameter is Λ ϕ = ε λ ϕ and cumulative inelastic variable Λ c = ε λ c .Thus, Rock material in engineering is mainly subjected to compression, and the strength needs to be determined within the damage-friction coupling framework.It is generally believed that inelastic strain causes the strengthening behavior of the material, while damage causes the strain softening after the peak stress.Therefore, two competing nonlinear mechanical mechanisms occur in the coupling process of inelastic strain and damage.It is difficult to determine the analytical form of material strength for plastic damage coupled models.
For the conventional triaxial compression loading path, in the principal stress space, the stress tensor is σ = σ 1 σ 2 σ 3 .At the same time, assuming that σ 1 < σ 2 = σ 3 , the deviatoric stress s is When the applied stress increases monotonously, the flow direction V can be written as According to the relation s c = K : σ c and p c = trσ c /3, we can obtain the following equations Substituting Equation (20) into Equations ( 13) and (15), the damage criteria in the case of crack closure can be obtained as follows: For simplicity, let us define χ = k m α 2 2α 1 + µ m α 2 .Substituting Equations ( 23) and (24) into Equation ( 9), one obtains From f ϕ = 0, it follows that Substituting Equation ( 27) into Equation (26) and using the sign convention in geotechnics, it is derived: Noting that p = (σ 1 + 2σ 3 )/3 and q = (σ 1 − σ 3 ), Equation ( 28) can be reformed into a loading function on the p-q surface as follows: When the damage variable reaches its critical value ϕ = ϕ c , the damage resistance function R(ϕ) reaches its maximum value R(ϕ c ), at which point the axial stress reaches its maximum value and the material reaches its peak strength.Based on the above analysis, the rock strength envelope predicted by the coupled friction-damage model is expressed in the following form: As mentioned above, geotechnical materials have hardening and softening properties.So, for a given loading path, R(ϕ) reaches its maximum value R(ϕ c ) when the damage factor reaches its critical value ϕ c , namely R(ϕ) < R(ϕ c ), and the rock is in the strainhardening phase.Meanwhile, for R(ϕ) > R(ϕ c ), the rock is in the strain-softening phase.The mathematical description of the above properties is R (ϕ) > 0 when 0 < ϕ < ϕ c ; R (ϕ) = 0 when ϕ = ϕ c ; R (ϕ) < 0 when ϕ < ϕ c .For this purpose, the dimensionless parameter ξ = ϕ/ϕ c is defined, and the following expression for the damage resistance R(ϕ) is used:

Description of the Returning Mapping Procedure
After the constitutive model and strength criterion are established, the inelastic strain ε c and damage variable ω are calculated iteratively according to the loading criterion.
The values of the k − 1 loading step variables ε k−1 , T k−1 , β k−1 , ϕ k−1 and σ k−1 are known.The flow for calculating the kth loading step ε k , T k , β k , ϕ k , σ k using strain loading is shown in Figure 5.

1.
The strain increment dε k was superimposed onto ε k−1 to estimate the strain ε k at step k, and the macroscopic stress σ k was preliminarily calculated according to Equation (3).

2.
The loading and unloading conditions are judged by the yield function f s (σ c ).If f s (σ c ) > 0, then by the consistency condition ḟ (ε, β, T) = 0 and Equation ( 19), we find λ c and hence T k , β k ; f s (σ c ) ≤ 0; then, we update the stress σ k according to ε k only.

3.
Calculate the damage driving force F ϕ(k) according to the updated T k , β k , based on Equation ( 14).
Substitute the updated ε k , T k , β k , ϕ k , σ k into the next loading loop.

Determination of Model Parameters
The coupled damage-friction model proposed in this paper contains only five parameters, E m , ν m , α, R(ϕ c ) and ϕ c , which can all be determined by a set of conventional triaxial mechanical tests.The Young's modulus and Poisson's coefficient are taken as E m = 30 GPa and ν m = 0.1 for the hard rock samples from the Pingdingshan coal mine.The remaining parameters were determined as follows: The parameter ϕ c included in the damage criterion is the damage threshold value, the value of which corresponds to the damage value at the peak stress.According to Lockner [47], the value of ϕ c is approximately linearly related to the perimeter pressure.For simplicity, a constant value of ϕ c = 1.5 is taken here.The effect of the value of ϕ c on the numerical simulation results will be discussed through parametric sensitivity analysis.As mentioned above, when ω = ω c that is, R(ϕ c ) = R(ϕ c ), the material reaches its maximum axial stress (the intensity envelope on the p-q surface is shown in Figure 6).Based on the protected hard rock test data, the parameter values can be determined by applying Equation (30):

Numerical Simulations for Diorite in Pingdingshan Coal Mine
Utilizing the values of the model parameters established in the preceding subsection, the stress-strain data from conventional triaxial compression tests are numerically simulated for different envelope pressure conditions (envelope pressure P c = 0, 10, 20 and 30 MPa), as shown in Figures 7-10.The model more accurately simulates the main macromechanical behavior of the rock, and it better describes the strength and stress-softening characteristics of the rock under different circumferential pressure conditions, especially in the elastic phase.Both in the axial and lateral directions, the stress-strain relationship is more accurately described; after the material enters non-linear deformation, before the peak strength, the model also fits the lateral deformation of the rock relatively well.According to the above simulated results, it can be seen that the constructed damage-friction coupling model in this paper can capture the main properties of the stress-strain curves of Pingdingshan hard rock under different confining pressures (envelope pressure P c = 0, 10, 20 and 30 MPa).Especially for confining pressures of 0 MPa, 10 MPa, and 20 MPa, this model can better simulate the peak pressure and related damage.However, when the confining pressure is 30 MPa, the simulation result of the height strength is smaller than the experimental data.This is because the rock sample in this experiment is obtained by the blasting method, so some microcracks exist inside the body, resulting in a certain dispersion of its peak strength.It must be pointed out that small but noticeable differences on the lateral strain are observed in the above figures.This is due to the fact that the tested tock is collected close to the coal seam, so there is more coal inclusions distributed in the rock solids, resulting in a large discrete property of the rock samples.Since the rock samples in the Pingdingshan Coal mine are obtained by an explosive method, unnatural cracks appear in the rock samples, which may lead to irregular lateral deformation.For most geomatrials, non-associated plastic flow rules are usually adopted to study the lateral and volumetric strains.However, in the present study, for the sake of simplicity, an associated flow law is adopted, which may result in a large difference between the experimental and analytical results of lateral strain.

Sensitivity Analysis of Parameter ϕ c
In order to determine the effect of the parameter ϕ c on the numerical simulation results, a conventional test with an envelope pressure of 10 MPa is used as the subject of the study.Figure 11 compares the results of the simulations for ϕ c values of 1.5, 2.5 and 3.5, respectively.It can be seen that the higher the ϕ c is, the greater the axial and lateral strains at the peak stress point are and the smoother the hardening and softening curve, namely, the slower the rate of hardening and softening, due to the larger the critical value of the damage.In addition, the number of microcracks within the characteristic cell or their radius is greater when the critical value is reached, which is the same for the material deterioration and the axial and lateral strains.The strengthening speed becomes smaller when the inelastic strain between the elastic stage and the peak stress point of REV become larger.As mentioned before, the magnitude of the parameter ϕ c is approximately linearly related to the surrounding pressure.For the sake of simplicity, a constant value was taken in the numerical simulations, which need further improvement of the simulation results for different envelope pressures.At the same time, due to the associated flow rule, the magnitude of lateral strain is measured when the stress peak point is reached at the stress-softening rate.The smaller ϕ c is, the faster the lateral strain reaches the peak, and the faster the stress-softening rate.

Numerical Simulations for Dagangshan Diabase
For further validation of the proposed model, comparisons between experimental results of the Dagangshan diabase and the numerical simulations are carried out in this subsection.The experimental data are provided by an earlier work [48].Now, we will provide numerical simulations of standard triaxial compression tests conducted with confining pressures of p c = 10, 15, 20, 30, 50 MPa.The specific model parameters used are listed in Table 1 below.The results are illustrated in Figures 12-16 with respect to confining pressures.It can be seen that the predicted numerical solutions for both low and high confining pressure levels match quite well with the experimental data.

Conclusions
Based on homogenization and irreversible thermodynamics theories, a coupled damage-friction model for Pingdingshan hard rock in a protective coal seam and Dagangshan diabase are constructed by the use of an associated flow law.The model successfully simulates the main mechanical behaviors of the tested rock in a conventional triaxial compression test.In order to assess the proposed coupled plastic-damage model, numerical simulations of triaxial compression tests on Pingdingshan rock and Dagangshan diabase have been carried out.By comparing with the experimental data, the proposed model can well describe the peak strengths and the mechanical responses from low to high levels of confining pressure for two types of rocks.However, the lateral deformation around the peak stress is not accurately predicted by the model due to the associated plastic flow rule.Determination of the critical damage parameter is also crucial, which determines the peak strength Determination of the free energy expressions for the RVE of cracked rock based on homogenization theory is the key point of the research results, which are combined with thermodynamic theory to the closed damage-friction coupling and strain strengthening and weakening processes.The proposed multiscale intrinsic model has the advantage of having only five parameters, which have specific physical meanings, and can be easily determined by the experiments.Inelastic deformation and the propagation of cracks are two of the fundamental mechanisms of material damage and destruction.Therefore, the coupled damage-friction model is more suitable for describing the mechanical behavior of such quasi-brittle material.In geotechnics, non-associated flow law is generally considered to be necessary.This model is suitable for quasi-brittle rock and closed fracture.Secondly, the selection of damage parameters in this model is related to the fracture density.After rock is compressed, cracks in its interior will be derived and developed, which is consistent with the material failure mechanism described above.However, in this study, the associated flow law of inelastic strain is adopted to predict the deformation of the tested rock.

Figure 1 .
Figure 1.Hard rock samples obtained from the protective coal seam of the Pingdingshan Mine.

Figure 2 .
Figure 2. Autonomous and auto-compensated multi-field coupling testing system.

Figure 3 .
Figure 3.The stress-strain curves of the hard rock in a protective coal seam with respect to different confining pressures.

Figure 4 .
Figure 4. Representative Elementary Volume of hard rock considering microcracks.

YesNoFigure 5 .
Figure 5. Flow chart of the numerical returning mapping procedure.

Figure 6 .
Figure 6.Strength envelope of the Pingdingshan hard rock under triaxial compressions.

Figure 7 .Figure 8 .Figure 9 .
Figure 7. Simulation of conventional triaxial compression test on hard rock protection layer rock in the Pingdingshan coal mine with the confining pressure of 0 MPa.

Figure 10 .
Figure 10.Simulation of conventional triaxial compression test on hard rock protection layer rock in the Pingdingshan coal mine with the confining pressure of 30 MPa.

Figure 11 .
Figure 11.Sensitivity analysis for parameter ϕ c

Figure 12 .Figure 13 .
Figure 12.Simulation of conventional triaxial compression test on Dagangshan diabase with the confining pressure of 10 MPa.

Figure 14 .
Figure 14.Simulation of conventional triaxial compression test on Dagangshan diabase with the confining pressure of 20 MPa.

Figure 15 .Figure 16 .
Figure 15.Simulation of conventional triaxial compression test on Dagangshan diabase with the confining pressure of 30 MPa.

Table 1 .
Identified parameters in the numerical simulations for Dagangshan diabase.