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

38600257
58872
10.1038/s41598-024-58872-4
Article
A damage constitutive model of layered slate under the action of triaxial compression and water environment erosion
Jia Jian 12
Tao Tiejun tjtao@gzu.edu.cn

1
Tian Xingchao 3
Xie Caijin 3
Jian Bingxi 3
Li Guoqing 3
1 https://ror.org/02wmsc916 grid.443382.a 0000 0004 1804 268X College of Mining, Guizhou University, Guiyang, 550025 China
2 Guizhou Power Transmission and Transformation Co., Ltd, Guiyang, 550002 China
3 https://ror.org/02wmsc916 grid.443382.a 0000 0004 1804 268X College of Civil Engineering, Guizhou University, Guiyang, 550025 China
10 4 2024
10 4 2024
2024
14 84454 5 2023
3 4 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Based on the macroscopic structure control theory, The slate with a significant bedding plane is a composite rock mass composed of rock blocks containing microscopic defects, joint surface closure elements, and shear deformation elements. Considering the coupling damage effect of water erosion and triaxial compressive load on bedding structure plane, the transversely isotropic damage constitutive model of slate under triaxial compressive load is derived with the dip angle of bedding and confining pressure as the variable. Firstly, based on the statistical theory of continuous damage mechanics and the maximum tensile strain criterion, the transversely isotropic deformation constitutive model of rock block with micro-defects is given; Secondly, based on the phenomenological theory of closed deformation and shear-slip deformation mechanism of layered structural plane under the coupling action of water erosion and triaxial compression load, the calculation formula of axial deformation of layered structural plane under the coupling action is given; Finally, to verify the accuracy of the established constitutive model, triaxial compression tests are carried out to study the influence of dip angle and confining pressure on the macroscopic mechanical properties and mechanism of slate. The results show that: the established triaxial compression damage constitutive model of bedding slate can accurately describe the stress–strain relationship of bedding slate after water environment erosion. With the increase of bedding dip angle, the strength and deformation capacity of the bedding slate first decreases and then increases, showing a U-shaped distribution as a whole. There are three main types of failure: tension shear composite failure, shear slip failure, and splitting tension failure.

Keywords

Bedding slate
Triaxial compression
Transversely isotropic
Deformation elements
Damage constitutive model
Subject terms

Civil engineering
Geology
http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 52064008 Tao Tiejun issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Various anisotropic rocks, such as slate, limestone, and shale, are often encountered in the excavation of underground projects such as tunnels1–3. Among them, the dip angle of the stratigraphic plane of slab rock is the key factor affecting the macroscopic properties of rock mass mechanics and the stability of surrounding rock. In the water-rich environment, the mineral composition of the rock mass changes after the rock mass is eroded by groundwater, and the microstructure is damaged to varying degrees, resulting in the macroscopic mechanical properties of the rock mass being affected4–6. In the neighborhood of rock engineering, theoretical analysis, laboratory experiments, numerical simulation, and other methods are usually used to explore the mechanical properties and damage mechanisms of rock. The mechanical properties of slate are a hot topic7–15. It is of great theoretical value and engineering guiding significance to consider the influence of the bedding dip angle of slate and the erosion of the water environment on the mechanical properties and damage development law of slate.

The test results of many indoor rock samples show that slate belongs to metamorphic rock, with obvious bedding structure, developed fissures, softening in water, and significant transverse isotropy. Based on uniaxial compression test, triaxial compression test, brazilian test, and set-up for wave velocity measurement, Alejano et al.16,17 studies the failure mechanism and mechanical properties of transversely isotropic rocks, gives an experimental determination method of five elastic parameters in the transversely isotropic elastic model, and verifies its reliability and applicability. Dambly et al.18 proposed a method to directly measurement the shear moduli of transversely isotropic rock mass by uniaxial compression test. Through uniaxial, triaxial, and brazilian tests, Gholami19 studied the influence of the change of bedding dip angle on the strength and deformation capacity of slate under wet and dry conditions. Ding et al.20 studied the influence of bedding inclination on the mechanical properties of slate through a uniaxial compression test, observed the morphology of fracture surface through scanning electron microscope (SEM), and explained its fracture mode and micro failure mechanism. Xu et al.21 carried out uniaxial compressive creep test combined with particle discrete element method simulation to study the creep behavior of phyllite under different water content. Chen et al.22 carried out triaxial compression test to study the influence of bedding dip angle and water content on the anisotropic mechanical behavior of phyllite, using the temporal and spatial distributions of the AE counts acoustic emission experiments to study the process of crack propagation and propagation to failure, the morphology and roughness of the fracture surface were analyzed by scanning electron microscope (SEM).

In theoretical research, based on continuum damage mechanics, fracture mechanics theory, and damage statistical distribution law, many scholars have established a variety of damage constitutive models23–36, and their research content mainly includes four aspects: (1) Select the failure strength criteria of rock microelements, such as maximum principal strain, Drucker–Prager, Hoek–Brown, and Mohr–Coulomb criterion; (2) Determine the distribution law of rock microelement strength, such as power function distribution, Weibull distribution, and lognormal distribution; (3) The determination methods of statistical parameters, such as linear fitting method, inverse analysis method, and peak point method; (4) The accuracy and applicability of the model are verified by indoor experiments and numerical simulation. Based on the Drucker–Prager criterion, Chen28 proposed a statistical damage constitutive model of triaxial compression of acid etched rock under the coupling effect of temperature and confining pressure. Based on the Hoek–Brown criterion, Chen29 established a constitutive model of rock triaxial compression damage, which can better predict the strain softening characteristics of rock and has a simple expression. Based on the modified Hoek–Brown criterion, Zhou30 established a statistical damage constitutive model of rock considering the influence of joint inclination. Fu31 proposed a statistical damage constitutive model that can accurately predict the triaxial compressive strength of anisotropic rocks, which is related to the number of thawing cycles, confining pressure, and bedding direction. Li et al.24,32 considered the strain softening phenomenon of materials and the influence of intermediate principal stress and established the triaxial compression damage constitutive model of soft rock. Liu33 established the triaxial compression damage constitutive model of layered composite rock and determined the model parameters by the curve fitting method.Based on the fracture theory and Weibull distribution, Liu34 established a damage constitutive model of jointed rock mass with micro defects and macro coupling under uniaxial compression. Liu35 proposed a damage constitutive model of non-penetrating jointed rock under biaxial compressive load, considering the influence of joint inclination and confining pressure. Based on the Drucker–Prager criterion, Chen36 established a coupled damage constitutive model of non-penetrating jointed rock mass, which considered the crack propagation length and joint friction effect. Wang37 considered the closure and shear deformation of the bedding structural plane and proposed a mathematical model that can simulate the complete stress–strain curve of a multi-jointed rock mass. Liu38 deduced the micro defects and macro coupling damage variables and established the damage statistical constitutive model of rock through the joint plane, which considered the shear strength of joints and the influence of microcracks on the strength and deformation capacity of a rock mass. Yang et al.39–41 Believed that the rock deformation with significant bedding structure should include the deformation of the joint surface. Taking into account the micro defects of rock blocks and the deformation of the joint surface, they proposed a static constitutive model, which simulates the rock with a significant bedding structure by connecting the rock block damage body with micro defects in series with the bedding deformation element.

The research on the damage constitutive model of slate mainly focuses on the influence of failure strength criterion and bedding dip angle on mechanical properties. Many results show that the transversely isotropic damage constitutive model is suitable for bedding slate. For slate in a water-rich environment, few damage constitutive models characterize the influence of bedding dip angle and confining pressure on the macroscopic mechanical properties of bedding slate under triaxial compression load. Based on the structural control theory, the influence of bedding dip angle and confining pressure on the macroscopic mechanical properties and damage evolution of anisotropic slate is considered in this paper. In Cartesian coordinates, the elastic constitutive equation of transversely isotropic is adopted to equivalent bedding slate, and its independent elastic parameters are changed from 2 to 5 when isotropic. Based on phenomenological theory, it is considered that the damage to bedding structure plane under compressive load includes the initial damage caused by water erosion before loading and the damage caused by deformation and instability of bedding during loading. Based on the calculation equations of closure deformation and shear slip deformation of layered structural plane under triaxial compression load, an effective and simple mathematical expression is established to characterize the stress–strain curve of layered slate under triaxial compression load after water erosion. In addition, to verify the accuracy and reliability of the established damage constitutive model and to analyze the influence of changes in confining pressure and bedding on the damage evolution process, mechanical properties, and failure mechanism of slate specimens, triaxial compression tests of laminated slate specimens with different dip angles and confining pressures carried out.

Damage evolution equation

Existing research shows that the significant slate deformation of bedding structures should include rock block deformation, closure deformation of the bedding structural plane, and shear slip deformation39,41. Based on the theory of rock mass structure, this paper uses rock blocks containing micro defects, closed deformation elements of bedding discontinuities, and shear slip deformation elements in series to simulate the deformation model of laminated slate11,41. To establish a relatively complete structural mechanics system, as shown in Fig. 1.Figure 1 Layered slate deformation element combination model diagram. (a) Rock blocks containing micro defects; (b) Closed deformation elements; (c) Shear slip deformation elements.

The axial stress and strain of the bedding slate specimen under triaxial compression load are σz and εz, respectively. The combination relationship shown in Fig. 1, It can be expressed as Eq. (1)1 σz=σr=σjc=σjsεz=εr+εjc+εjs

In Eq. (1), σr is the axial stress of rock block, σjc is the axial stress of joint closure deformation, and σjs is the axial stress of joint shear deformation; εr is the axial strain of rock block, εjc is the axial strain during joint closure deformation, and εjs is the axial strain during joint closure deformation.

The main task of establishing a damage constitutive model is to give mathematical expressions that accurately describe the macroscopic mechanical behavior and damage evolution of three deformation elements under triaxial compression. The subsequent contents of this section give the damage evolution equations of the two deformation elements of the rock damage body and the bedding structural plane under load and solve them according to the combination of Eq. (1).

Damage variable of rock blocks containing micro defects

Constitutive model of transversely isotropic elastic

When the bedding structure plane of the slate is stable, it is feasible to analyze the rock block with a transversely isotropic elastic constitutive model16. Establish the coordinate system as shown in Fig. 2. The global coordinate system XYZ rotates clockwise around the X axis to form a local coordinate system X′Y′Z′. The angle between the normal direction of the isotropic surface (joint dip structure surface) and the Z axis is the anisotropy angle β(the bedding dip angle).Figure 2 Diagram of relationship between joint dip angle and coordinate. (a) Definition of anisotropy angle β; (b) Local X′Y′Z′ and global XYZ coordinate systems.

In the Cartesian coordinate system, the generalized Hooke’s law can be used to describe the transversely isotropic elastic constitutive relation of slate, and its independent elastic parameters are reduced from 21 to 5 in the case of extreme anisotropy. In the local coordinate system, the stress–strain relationship is as follows17,42:2 εr′=S′σr′

In Eq. (2), εr′ is the strain tensor; σr′ is the stress tensor; S′=A00B is the elastic flexibility matrix, where A=11E1E1-v1-v1E1E1-v2-v2E2E2-v1-v1E1E111E1E1-v2-v2E2E2-v2-v2E2E2-v2-v2E2E211E2E2 and B=diag2(1+v1)2(1+v1)E1E1,11G2,11G2G2G2,11G2G2. In the matrix, E1 and v1 are the elastic modulus and Poisson’s ratio in the isotropic plane, respectively; E2, v2 and G2 represent the elastic modulus, Poisson’s ratio, and shear modulus perpendicular to the isotropic plane, respectively.

Similarly, the stress–strain relationship of a transversely isotropic body in the global coordinate system XYZ is as follows:3 εr=Sσr

4 S=S11S12S13S14S15S16S21S22S23S24S25S26S31S32S33S34S35S36S41S42S43S44S45S46S51S52S53S53S55S56S61S62S63S64S65S66

In Eqs. (3): εr represents the strain tensor, σr refers to stress tensor; S is the elastic flexibility matrix, and each component can be obtained by the transformation Eq. (5) between the local coordinate system and the global coordinate system.5 S=RTS′R

In Eq. (5), R is the coordinate transformation matrix, which can be expressed as:6 R=1000000cos2βsin2β0-sin(2β)00sin2βcos2β0sin(2β)0000cosβ0sinβ0sin(2β)/2-sin(2β)/20cos(2β)0000sinβ0-cosβ

Through Eqs. (5) and (6), the elastic flexibility matrix components required in this paper are obtained, such as Eqs. (7)–(9):7 S13=S31=-v1E1sin2β-v2E2cos2β

8 S23=S32=sin2(2β)41E1+1E2-1G2-v2E2(sin4β+cos4β)

9 S33=sin2(2β)41G2-2v2E2+cos4βE2+sin4βE1

The elastic parameters in the elastic flexibility matrix of rock mass are inherent properties18. The five elastic parameters in the above elastic flexibility matrix are inherent properties of slate, and the theoretical values of the five elastic parameters of transversely isotropic rock block can be determined by the triaxial compression test of slate specimens. Slate samples with joint dip angles of 0°, 45° and 90° are selected for static compression test to determine five independent elastic parameters. In Table 1, the calculation Eqs. (10)–(18) are derived from Eqs. (4)–(9).Table 1 Five elastic determination equations of transversely isotropic slate.

The joint dip angle β(°)	Calculating equations	
β=0∘	ΔεxΔσr=-v2E2 (10)	
ΔεyΔσz=-v2E2 (11)	
ΔεrΔσr=1E2 (12)	
β=90∘	ΔεxΔσr=-v1E1 (13)	
ΔεyΔσr=-v2E2 (14)	
ΔεrΔσr=1E1 (15)	
β=45∘	ΔεxΔσr=-12v1E1+v2E2 (16)	
ΔεyΔσr=141E1+1E2-1G2-v22E2 (17)	
ΔεrΔσr=141E2+1E1+1G2-2v2E2 (18)	

Damage evolution equation of rock blocks containing micro defects

Based on the theory of statistical damage mechanics, it is considered that the micro defects in the layered slate sample can be regarded as random damage, and its internal damage development has randomness and complexity under triaxial compression load29. The joint slate is composed of numerous micro-elements, and the macroscopic mechanical properties are the comprehensive reaction of many micro-element combinations. Assuming that the undamaged part of the micro-element bears the load, 38,43 defining the damage variable Dm as the ratio of the number of micro-elements of damage failure Nf to the total number of micro-elements Nt, and the statistical damage variable is expressed as Eq. (19).19 Dm=NfNt

The strength of micro elements obeys the maximum tensile strain failure criterion, and the damage obeys the random distribution and has statistical rules. The loading process of rock block material is a continuous damage process, and the Weibull distribution function30 can be used to quantify the damage degree, as Eq. (20):20 Pεr=mF0εrF0m-1exp-εrF0m,εr>00,εr≤0

In Eq. (20), Pεr is the probability density function, and εr is the random distribution variable of micro-element strength; m and F0 are the parameters of the Weibull distribution that reflect the degree of concentration and the strength magnitude of the micro-elements, respectively34. With the increase of load, when the strain of the rock sample reaches a certain degree, the statistical expression of the number of rock damage failure micro-elements is Eq. (21).21 Nf(εr)=∫0εNtPεrdε=Nt1-exp-εrF0m

With the increase of load, the internal damage of rock block continues to accumulate, micro-cracks appear and expand, and finally form macro cracks, resulting in the failure of a rock mass. The damage variable Dm is the macro expression of rock micro-defects, with a value between 0 and 1, corresponding to the state that the rock loses its bearing capacity from initial loading to complete failure, substituting Eq. (21) into Eq. (19), the damage variable Dm as Eq. (22):22 Dm=1-exp-εrF0m

Under the action of triaxial compression load34, Hooke’s law is followed when the micro-element of the rock block is loaded, and the transverse isotropic damage constitutive model of the rock block (Fig. 1a) without considering the bedding deformation can be obtained, as Eq. (23):23 εr=Sσr1-Dm=(S13σx+S23σy+S33σr)exp-εrF0m

After the transformation of the damage statistical equation, it is not necessary to determine the number of micro-elements in Eq. (19). Expression of microscopic defects of rock block by macroscopic damage variable D This paper mainly discusses the macroscopic mechanical properties and damage evolution law of slate. For this purpose, the determination of damage parameters F0 and m in Eq. (23) should also be given.

Determination of damage parameters

In Eq. (23), F0 and m in the damage constitutive model are the statistic parameters that characterize the concentration degree and macroscopic mechanical strength of slate micro-elements, which is an inherent property. The common methods to determine the statistical parameters based on the stress–strain curve obtained by the triaxial compression test are a linear fitting method, peak point method, and inversion analysis method34,36,38. The three methods have their advantages and disadvantages. In this paper, the peak point method is selected. The solution process includes complex derivation, but the peak point of the test stress–strain curve is well simulated. The implicit constitutive Eq. (24) is derived to obtain Eq. (25)24 S13σx+S23σy+S33σr=εrexp-εrF0m

25 dσrdεr=S-1expεrF0m1-mεrF0m

According to the geometrical characteristics of the stress–strain curve of slate specimen under triaxial compression test, 43 it can be known that the peak point dσrdεr=0 and the vertical Eqs. (24) and (25) can be used as parameters in the damage constitutive equation of rock block under triaxial compression load, as Eq. (26):26 m=1/lnεrS33σrF0=εrm-1/m

Deformation constitutive equation of bedding structural plane

With the change in ion concentration, pH value, and temperature in groundwater, the physical and chemical reactions of micro-elements in slate are complex and changeable44,45. It is difficult to quantify the influence of water erosion on the strength and deformation ability of the bedding structural plane from the microscopic point of view. Phenomenological theory shows that the change of macroscopic physical properties of the layered structural plane can characterize its deterioration degree46. It is generally considered that the deterioration degree of closed elastic modulus and shear modulus is the same when axial deformation and shear deformation occur in the layered structural plane under compressive load.

The damage of the bedding structure plane of slate under the coupling action of water erosion and triaxial compression load mainly includes two parts: (1) Bedding damage caused by water erosion before loading; (2) damage caused by bedding closure deformation and shear deformation instability during loading. According to the stress–strain curve obtained from the test, the change in the elastic modulus of the bedding structural plane is easy to analyze and measure. Under the action of water environment erosion, the damage variables of the bedding structural plane are defined as Eq. (27):27 Dc=1-EjcEj0

In Eq. (27), Dc is the damage variable of water environment erosion; Ej0 represents the closed elastic modulus of rock joint structural plane in its natural state, and Ejc represents the closed elastic modulus of the bedding structural plane of rocks eroded by the water environment.

When carrying out triaxial compression tests on slate specimens with a joint dip inclination of 0°, The shear slip deformation value of the bedding structural plane is 0. The small stress increment at the initial stage of loading corresponds to the large strain increment, that is, the initial stage of the stress–strain curve is relatively flat. Here, the measurement method of the closed elastic modulus Ej0 and Ejc of the bedding structural plane is described. It is approximately considered that the axial deformation of the specimen at the initial stage of the stress–strain curve is mainly the closed deformation of the bedding structure. Therefore, the slope of the initial gentle stage of the stress–strain curve is approximately regarded as the closed elastic modulus of the bedding structural plane.

The deformation of the stratified slate under triaxial compressive loading includes the deformation of the rock block, the closing deformation of the bedding structural plane, and the shear slip deformation of the upper rock block along the joint plane. The mechanical model of laminated slate under triaxial compression load is shown in Fig. 340. Under the action of three-dimensional stress47, the compressive stress and shear stress on the bedding structural plane of slate can be expressed by the Eq. (28):28 σβ=σjccos2β+σycosαsinβτβ=σjscosβsinβ-σycos2β

Figure 3 Diagram of triaxial compression stress of slate.

The contribution of the closure deformation and shear deformation of the bedding plane to the axial deformation of the specimen is discussed, respectively.

First, consider the axial closure strain εjc generated by the bedding structural plane, as shown in Eq. (29)40,47, and εj0 is the maximum closure strain of the structural plane. Under the triaxial compression load, the elastic modulus of the rock blocks containing micro defects is much larger than the closed elastic modulus of the bedding structural plane for the slate specimen with a joint dip angle of β = 0°. It is approximately considered that the strain value at the initial stage of the stress–strain curve is the maximum closed elastic modulus of the bedding structural plane.29 Ejc=Ej0(1-Dc)εjc=εj01-exp-σjccos2β+σycosβsinβEjcεj0·cosβ

Secondly, the shear slip of the bedding plane is considered to affect the total axial strain of the specimen. According to the constitutive relationship of bedding structural plane under shear load39,46, the deformation calculation expression is derived:30 τβ=Kscδs

In Eq. (30), τβ is the shear stress on the joint surface; Ksc represents the tangential stiffness of the joint surface of rock eroded by the water environment; δs is the shear displacement along the joint surface. The axial strain εjs of the specimen due to shear slip along the joint surface is as follows:31 Ksc=Ks0(1-Dc)εjs=σjscosβsinβ-σycos2βsinβKscL

In Eq. (31), Ks0 represents the tangential stiffness of the rock bedding structural plane in the natural state, and Ksc represents the tangential stiffness of the rock bedding structural plane eroded by the water environment.

When the dip angle of bedding is 45°, the deformation of the layered structural plane under triaxial compression load includes closed deformation and shear slip deformation Similarly, in the initial stage of loading, the axial deformation of the specimen is approximately considered to be composed of the closure deformation and shear slip deformation of the bedding structural plane The shear modulus Ks0 and Ksc can be calculated according to the slope of the initial stage of stress–strain curves of slate specimens with bedding dip angles of 0° and 45°.

In Eqs. (29) and (31), the damage variable of water erosion is given. The next step is to consider the damage evolution law of macroscopic mechanical properties of the bedding structural plane under triaxial compression load and give the calculation expression of damage coupling deformation of the bedding structural plane under compression load after water environment erosion. In this paper, the state of the bedding structural plane after water environment erosion is regarded as the initial state, and the damage state caused by compression load after water environment erosion is the second state. Under the action of triaxial compression load, it is considered that the damage deterioration law of the bedding structural plane is the same as that of rock block mechanical properties, and the coupling damage variable D of the bedding structural plane is expressed by the Eq. (32). Bring Eq. (32) into Eqs. (29) and (31) to obtain an Eq. (33) that can accurately describe the deformation of the bedding structural plane.32 D=(1-Dm)(1-Dc)=(1-Dc)exp-εF0m

33 εjc=εj01-exp-σjccos2β+σycosβsinβ(1-Dc)exp-εF0mEj0εj0cosβεjs=σjscosβsinβ-σycos2βsinβ(1-Dc)exp-εF0mKs0L

Establish triaxial compression damage constitutive model

Under the triaxial compression load, the axial deformation of laminated slate consists of three parts: the deformation of damaged rock mass containing micro-cracks, the closure deformation of the bedding structural plane, and shear slip deformation. In Sections “Damage variable of rock blocks containing micro defects” and “Deformation constitutive equation of bedding structural plane”, the detailed derivation process of the deformation calculation expression of the three components is given, respectively. Bring Eqs. (23) and (33) into Eq. (1) to obtain the triaxial damage constitutive model of laminated slate after water environment erosion under a triaxial compression load, as shown in Eq. (34).34 εr=(S13σx+S23σy+S33σz)exp-εrF0mεz=εr+εj01-exp-σzcos2β+σycosβsinβ(1-Dc)exp-εrF0mEj0εj0cosβ+σzcosβsinβ-σycos2βsinβ(1-Dc)exp-εrF0mKs0·L

The Eq. (34) contains three parts corresponding to the axial deformation of rock block under the coupling action of water environment erosion and triaxial compression load, the axial deformation caused by the closure deformation of the bedding structure plane, and the shear slip, respectively.

Experimental work and results

To verify the accuracy and applicability of the damage constitutive model established in this paper, and explore the damage evolution law, failure characteristics, compressive strength, and deformation capacity of laminated slate after water environment erosion under triaxial compression load. The indoor triaxial compression experiments of layered slate samples with different confining pressures and bedding angles were carried out to analyze the effects of confining pressures and inclination changes on the damage evolution law, failure characteristics, compressive strength, and deformation capacity of the samples.

Experiment preparation

To reduce the discrete type of rock samples and avoid the differences in composition and structure of rock samples, all slate samples are from the Bayue mountain tunnel of the Tongan expressway. According to the standards of the international society for rock mechanics (ISRM)48, the slate samples are processed into qualified samples. The flatness error of both ends of the specimen shall not exceed 0.02mm, and the side shall be smooth and straight, meeting the requirements of verticality, some rock samples are shown in Fig. 4.Figure 4 Layered slate samples with different dip angles.

To reveal the damage evolution mechanism of water-saturated laminated slate under triaxial compression load. Carry out triaxial compression tests of slate specimens after water environment erosion at confining pressures of 0 MPa, 5 MPa, 10 MPa, and 15 MPa. The statistics of the number of samples are shown in Table 2.Table 2 Statistics of rock samples.

Confining pressure (MPa)	0°	30°	45°	60°	90°	Total	
0	3	3	3	3	3	15	
5	3	3	3	3	3	15	
10	3	3	3	3	3	15	
15	3	3	3	3	3	15	
Total	12	12	12	12	12	60	

The equipment used is a TAJW-2000 microcomputer controlled electro-hydraulic servo rock triaxial test system. Its loading control mode is equal displacement rate, and the loading rate is 0.005 mm/s, as shown in Fig. 5.The system consists of the main frame, hydraulic pump unit, controller, strain gauge, data acquisition module, and computer.Figure 5 Microcomputer controlled electro-hydraulic servo rock triaxial test system.

Experimental results

Determine elastic parameters

Using 9 rock samples to carry out compression experiments with confining pressure of 0 MPa, five elastic parameters of transverse isotropy of rock blocks are determined, and the damage variables of closure modulus and shear modulus of bedding structural plane are obtained.

The method to determine the five independent elastic parameters of slate is as follows. As shown in Fig. 6, take the reciprocal of the slope that is the stress–strain curve at 30–50% of the peak stress, and determine five independent elastic constants in combination with the equation given in Table 1.Figure 6 Stress strain curve of test piece.

In Eq. (33), the closing deformation and shear slip deformation of the bedding plane of the slate specimen change with the change of the angle when the bedding dip angle is different. Current the anisotropy angle β = 0°, cosβ = 1, sinβ = 0, the axial component value of shear slip deformation of bedding structural plane is εjs = 0, E2 and v2 can be solved by Eqs. (11) and (12); Current the anisotropy angle β = 90°, cosβ = 0, sinβ = 1, The axial component value of shear slip deformation of the bedding structural plane is εjs = 0, and the axial component value of closed deformation is εjc = 0. The axial deformation of slate specimen is the axial deformation of Rock blocks containing micro defects, E1 and v1 can be solved by Eqs. (14) and (15); When the anisotropy angle 0° < β < 90°, the deformation of the specimen consists of three parts: the axial deformation of the rock block, the axial closure deformation, and shear slip deformation of the bedding structural plane, G2 can be solved by Eq. (18).The calculation results of the five independent elastic parameters are shown in Table 3:Table 3 Five elastic parameters of the wet slate obtained from uniaxial tests.

Statistics	E1 (GPa)	E2 (GPa)	v1	v2	G2 (GPa)	
Min	30.19	26.12	0.21	0.23	7	
Max	46.94	34.74	0.23	0.28	12	
Mean	37.24	29.33	0.22	0.25	9	

According to the experimental test method given in Section “Deformation constitutive equation of bedding structural plane”, three slate specimens with a bedding dip angle of 0°and 45° after natural state and water environment erosion are selected to carry out the uniaxial compression test, respectively. The maximum closure strain of the structural plane εj0 = 0.00098, the closed elastic modulusof the joint surface Ej0 = 1973 MPa, Ejc = 1712 MPa, the tangential stiffness of the joint surface Ksc = 1727 MPa/m, the tangential stiffness of the joint surface Ks0 = 1989 MPa/m, the value of damage variable Dc is 0.132.

After solving the five independent elastic parameters, bring the average value into Eqs. (7)–(9), and calculate the s component of the elastic flexibility matrix of slate with different bedding dip angles, the results are listed in Table 4.Table 4 Elastic flexibility matrix components.

Anisotropy angle β	S13	S23	S33	
0°	− 8.5237e−06	− 8.5237e−06	3.40948e−05	
30°	− 7.86968e−06	− 1.4733e−05	3.84949e−05	
45°	− 7.21566e−06	− 1.68027e−05	3.87582e−05	
60°	− 6.56164e−06	− 1.4733e−05	3.48848e−05	
90°	− 5.90763e−06	− 8.5237e−06	2.68745e−05	

Stress–strain curve

According to the above experimental scheme, the triaxial compression test of laminated slate is carried out, and the stress–strain relationship curve is obtained and compared with the constitutive model established in this paper. The Stress–strain curve and damage evolution curve can completely reflect the mechanical behavior of laminated slate under a triaxial compression load, as shown in Fig. 7. The rising section of the theoretical calculation curve of layered slate with different bedding dip angles is consistent with the indoor triaxial compression test curve, which can better fit the deformation process of rock before failure, it shows that it is more reasonable to consider the deformation of slab strata. There are still some deviations in the post-peak softening stage due to the influence of rock-hard brittleness. After the specimen reaches the peak stress, the instability failure suddenly occurs. The experimental test equipment has no rigid auxiliary device, and the falling section of the stress–strain curve decreases linearly and rapidly.Figure 7 The stress–strain curves and damage variation curves of specimens under different bedding dip angles and confining pressures.

In Table 4, S11 and S23 are the confining pressure coefficients in the x and y directions, respectively. With the increase of bedding dip angle (0° ≤ β ≤ 90°), the values increase first and then decrease, and reach the maximum when β = 60°. S33 represents the coefficient of axial load. With the increase of bedding angle, the value decreases first and then increases, and it is the smallest when β = 60°. The change trend of its value is consistent with the change trend of compressive strength of the sample, as shown in Fig. 7.

The parameter m describes the strength distribution concentration of trace elements in the rock mass and reflects the plasticity and brittleness characteristics of the jointed rock mass. The larger the value is, the better the plasticity is, and 60° is the minimum value. The parameter F0 is the average value of the strength of trace elements in macro statistics, which reflects the strength of jointed rock mass, and the change rule of its value is consistent with the change rule of the compressive strength of the sample, and 60° is the minimum value.

In Fig. 7, the damage evolution curve of bedding slate shows an “S” type distribution law with the whole, which is composed of three segments, the first is a gentle rising segment, the second is a fast rising segment, and the last is a slow rising segment, which tends to 1. When the bedding dip angle is 0°–60°, the initial gentle rising section of slate with the same bedding dip angle becomes longer with the increase of confining pressure; With the increase of bedding dip angle, the initial gentle ascending section becomes shorter; With the increase in bedding dip angle, the steepness of the rising section decreases. When the bedding dip angle is 90°, the rule of damage evolution curve is the same as that of 0°, and the length of the gentle rising section and the steepness of the rapid rising section are between 0° and 30°. The changing trend of the damage curve of layered rock samples is related to many factors, such as bedding dip angle, stress state, failure characteristics, compressive strength, and deformation capacity.

Analysis of failure characteristics of specimens

The layered slate shows obvious brittleness and anisotropy when it fails. Figure 8 shows the failure type and failure section morphology of slate specimen under compression load. The results show that confining pressure and bedding dip angle play a decisive role in the failure type of slate.Figure 8 Failure diagram of specimen under different experimental conditions.

Current bedding dip β = 0°, the specimen has a tensile shear composite failure through the bedding plane, and the failure crack is wide. With the increase of confining pressure, the failure crack becomes narrow, and the failure type is mainly tensile failure. Current bedding dip β = 30°, the fracture surface is a composite failure of tension and shear along the bedding plane and through the bedding plane. With the increase of confining pressure, the crack develops in a stepped manner and expands fully. Current bedding dip β = 45°, the shear failure mainly occurs along the bedding plane, and the shear failure through the bedding occurs with the increase of confining pressure. Current bedding dip β = 60°, the shear slip failure occurred along the bedding plane, and the number of failure planes increased with the increase of confining pressure. Current bedding dip β = 90°, the specimen is fractured and tensioned along the bedding plane and parallel bedding plane. With the increase of confining pressure, the number of wide cracks leading to specimen failure increases, accompanied by more micro cracks.

In addition, it can be observed in Fig. 8 that when the failure type of the sample is mainly tensile, the failure section is dry, and the stratification of the failure section is significant; When the failure type of the specimen is shear slip failure along the bedding plane, the failure section is wet. It shows that when water environment erodes layered slate, the mechanical properties of its bedding structural plane will be greatly affected. Therefore, it is reasonable to consider the coupling effect of water environment erosion and compression load when establishing the damage constitutive model in this paper.

Analysis of compressive strength and deformation capacity

The anisotropy of parameters characterizing mechanical properties is significant with the change of bedding dip angle. In Fig. 9a, under the same confining pressure, with the increase of bedding dip angle, the triaxial compressive strength of laminated slate first decreases and then increases, and the peak strength of the specimen presents a U-shaped parabola with the opening upward, which is consistent with much research results. When the bedding dip angle is 0°, the compressive strength of the slate is the largest, and the compressive strength of the slate with a bedding dip angle of 60° is the smallest. With the increase of confining pressure, the compressive strength of slate with the same bedding angle also increases, and the post-peak residual strength of the specimen is almost zero, with significant brittle failure characteristics.Figure 9 Two-factor change diagram of mechanical properties; (a) The variation of triaxial compressive strength of slate specimens with different bedding angles and confining pressures; (b) The variation of deformation capacity (strain value corresponding to peak stress) of slate specimens with different bedding angles and confining pressures.

In Fig. 9b, shows that the deformation capacity of the specimen is closely related to the compressive strength. It can be seen from Fig. 7 that at the beginning of loading, the stress–strain curves of slate samples with different bedding angles have a nonlinear compaction stage, and then the elastic deformation stage. When the bedding dip angle is 0°–45°, the compaction stage of the stress–strain curve of the rock sample is obvious, and the compaction stage becomes smaller with the increase of the dip angle; When the bedding dip angle is 60°, the compaction stage of the rock sample does not exist, mainly because the slate sample is mainly sheared along the bedding plane under a compression load, and the rock mass is relatively complete at this time; The compaction stage of the 90° rock sample is not obvious, mainly because the influence of the bedding structure of the sample on the compressive strength and deformation capacity is not considered. With the increase of confining pressure, the deformation capacity of slate with the same bedding angle increases, and the brittle failure characteristics are obvious.

Conclusion

The damage constitutive model of slate under triaxial compression load is proposed, which includes the coupling damage variables of water environment erosion and compression load. The triaxial compression experiments of slate with different confining pressures and bedding dip angles were carried out to verify the accuracy and applicability. The influence of confining pressure and bedding dip angle on the damage evolution law, compressive strength, and deformation capacity of laminated slate after water environment erosion is given. The main conclusions are as follows:It is assumed that the statistical law conforms to the Weibull distribution, and considering the influence of the erosion of the water environment, the change of confining pressure, and the bedding dip angle on the stress–strain relationship. The damage constitutive model of transversely isotropic slate under triaxial compression load is given, which is a concise expression that can accurately fit the compaction section and elastic section of the stress–strain curve.

The damage evolution curve of laminated slate is composed of three sections: gentle rising section, fast-rising section, and slow rising section. The overall distribution law is S-shaped. Confining pressure and bedding dip angle are important factors affecting the damage evolution law.

The triaxial compression test results show that when the bedding dip angle changes from 0° to 60°, the failure type of slate specimens develops from tension shear composite failure to shear slip failure; The specimens with a bedding dip angle of 90° undergo splitting tensile failure along the bedding plane and parallel bedding plane. In addition, with the increase of confining pressure, the number of main cracks increases when the slate specimen is damaged, and there are many micro-cracks.

With the change of dip angle, the compressive strength and deformation capacity of laminated slate change in a U-shape. When the bedding dip angle is 0°, the compressive strength and deformation capacity are the largest, and when the bedding dip angle is 60°, the compressive strength and deformation capacity are the smallest. With the increase of confining pressure, the mechanical properties of laminated slate are significantly improved, and the compressive strength and deformation capacity of specimens are improved.

Acknowledgements

This research is supported by the National Natural Science Foundation of China (Grant No.: 52064008).

Author contributions

J.J.: Conceptualization, Writing-original draft; T.T.: Writing-review & editing; C.X.: Data curation, Project administration; X.T.: Investigation, Resources; B.J. and G.L.: Software, Visualization. All authors reviewed the manuscript.

Data availability

All data, models, or codes that support the fifindings of this study are available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

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

1. Wu RJ Li HB Wang DP Full-field deformation measurements from Brazilian disc tests on anisotropic phyllite under impact loads Int. J. Impact Eng. 2020 149 6 103790 10.1016/j.ijimpeng.2020.103790
Wu, R. J., Li, H. B. & Wang, D. P. Full-field deformation measurements from Brazilian disc tests on anisotropic phyllite under impact loads. Int. J. Impact Eng. 149(6), 103790. 10.1016/j.ijimpeng.2020.103790 (2020).10.1016/j.ijimpeng.2020.103790
2. Yiouta-Mitra P Sofianos AI Μulti-jointed stratified hard rock roof analysis and design Int. J. Rock Mech. Min. Sci. 2018 106 96 108 10.1016/j.ijrmms.2018.03.021
Yiouta-Mitra, P. & Sofianos, A. I. Μulti-jointed stratified hard rock roof analysis and design. Int. J. Rock Mech. Min. Sci. 106, 96–108. 10.1016/j.ijrmms.2018.03.021 (2018).10.1016/j.ijrmms.2018.03.021
3. Gegenhuber N Krueger M Linking static and dynamic mechanical properties for metamorphic rocks from Austria including their anisotropic effect Acta Geophys. 2021 69 2 539 546 10.1007/s11600-021-00566-w
Gegenhuber, N. & Krueger, M. Linking static and dynamic mechanical properties for metamorphic rocks from Austria including their anisotropic effect. Acta Geophys. 69(2), 539–546. 10.1007/s11600-021-00566-w (2021).10.1007/s11600-021-00566-w
4. Pan JL Cai MF Li P A damage constitutive model of rock-like materials containing a single crack under the action of chemical corrosion and uniaxial compression J. Cent. South Univ. 2022 29 2 486 498 10.1007/s11771-022-4949-1
Pan, J. L. et al. A damage constitutive model of rock-like materials containing a single crack under the action of chemical corrosion and uniaxial compression. J. Cent. South Univ. 29(2), 486–498. 10.1007/s11771-022-4949-1 (2022).10.1007/s11771-022-4949-1
5. Noiriel C Renard F Doan ML Intense fracturing and fracture sealing induced by mineral growth in porous rocks Chem. Geol. 2010 269 3–4 197 209 10.1016/j.chemgeo.2009.09.018
Noiriel, C. et al. Intense fracturing and fracture sealing induced by mineral growth in porous rocks. Chem. Geol. 269(3–4), 197–209. 10.1016/j.chemgeo.2009.09.018 (2010).10.1016/j.chemgeo.2009.09.018
6. Miao SJ Cai MF Wang GQ Damage effects and mechanisms in granite treated with acidic chemical solutions Int. J. Rock Mech. Min. Sci. 2016 88 77 86 10.1016/j.ijrmms.2016.07.002
Miao, S. J. et al. Damage effects and mechanisms in granite treated with acidic chemical solutions. Int. J. Rock Mech. Min. Sci. 88, 77–86. 10.1016/j.ijrmms.2016.07.002 (2016).10.1016/j.ijrmms.2016.07.002
7. Shen PW Tang HM Zhang BC Investigation on the fracture and mechanical behaviors of simulated transversely isotropic rock made of two interbedded materials Eng. Geol. 2021 286 2 106058 10.1016/j.enggeo.2021.106058
Shen, P. W. et al. Investigation on the fracture and mechanical behaviors of simulated transversely isotropic rock made of two interbedded materials. Eng. Geol. 286(2), 106058. 10.1016/j.enggeo.2021.106058 (2021).10.1016/j.enggeo.2021.106058
8. Chiu CC Weng MC Huang TH Modeling rock joint behavior using a rough-joint model Int. J. Rock Mech. Min. Sci. 2016 89 14 25 10.1016/j.ijrmms.2016.08.001
Chiu, C. C., Weng, M. C. & Huang, T. H. Modeling rock joint behavior using a rough-joint model. Int. J. Rock Mech. Min. Sci. 89, 14–25. 10.1016/j.ijrmms.2016.08.001 (2016).10.1016/j.ijrmms.2016.08.001
9. Cao RH Cao P Lin H Failure characteristics of jointed rock-like material containing multi-joints under a compressive-shear test: Experimental and numerical analyses Arch. Civ. Mech. Eng. 2018 18 3 784 798 10.1016/j.acme.2017.12.003
Cao, R. H. et al. Failure characteristics of jointed rock-like material containing multi-joints under a compressive-shear test: Experimental and numerical analyses. Arch. Civ. Mech. Eng. 18(3), 784–798. 10.1016/j.acme.2017.12.003 (2018).10.1016/j.acme.2017.12.003
10. Zhou JW Xu WY Yang XG A microcrack damage model for brittle rocks under uniaxial compression Mech. Res. Commun. 2010 37 4 399 405 10.1016/j.mechrescom.2010.05.001
Zhou, J. W., Xu, W. Y. & Yang, X. G. A microcrack damage model for brittle rocks under uniaxial compression. Mech. Res. Commun. 37(4), 399–405. 10.1016/j.mechrescom.2010.05.001 (2010).10.1016/j.mechrescom.2010.05.001
11. Liu HY Su TM A dynamic damage constitutive model for a rock mass with non-persistent joints under uniaxial compression Chin. J. Geotech. Eng. 2016 77 12 20 10.1016/j.mechrescom.2016.08.006
Liu, H. Y. & Su, T. M. A dynamic damage constitutive model for a rock mass with non-persistent joints under uniaxial compression. Chin. J. Geotech. Eng. 77, 12–20. 10.1016/j.mechrescom.2016.08.006 (2016).10.1016/j.mechrescom.2016.08.006
12. Le LA Nguyen GD Bui HH Modelling jointed rock mass as a continuum with an embedded cohesive-frictional model-ScienceDirect Eng. Geol. 2017 228 107 120 10.1016/j.enggeo.2017.07.011
Le, L. A. et al. Modelling jointed rock mass as a continuum with an embedded cohesive-frictional model-ScienceDirect. Eng. Geol. 228, 107–120. 10.1016/j.enggeo.2017.07.011 (2017).10.1016/j.enggeo.2017.07.011
13. Yang TY Wang PT Xu T Anisotropic characteristics of jointed rock mass: A case study at Shirengou iron ore mine in China Tunn. Undergr. Sp. Technol. Inc. Trenchless Technol. Res. 2015 48 apr. 129 139 10.1016/j.tust.2015.03.005
Yang, T. Y. et al. Anisotropic characteristics of jointed rock mass: A case study at Shirengou iron ore mine in China. Tunn. Undergr. Sp. Technol. Inc. Trenchless Technol. Res. 48(apr.), 129–139. 10.1016/j.tust.2015.03.005 (2015).10.1016/j.tust.2015.03.005
14. Liu Y Dai F A damage constitutive model for intermittent jointed rocks under cyclic uniaxial compression Int. J. Rock Mech. Min. Sci. 2018 103 289 301 10.1016/j.ijrmms.2018.01.046
Liu, Y. & Dai, F. A damage constitutive model for intermittent jointed rocks under cyclic uniaxial compression. Int. J. Rock Mech. Min. Sci. 103, 289–301. 10.1016/j.ijrmms.2018.01.046 (2018).10.1016/j.ijrmms.2018.01.046
15. Ahmed Z Wang S Hashmi MZ Causes, characterization, damage models, and constitutive modes for rock damage analysis: A review Arab. J. Geosci. 2020 13 16 1 14 10.1007/s12517-020-05755-3
Ahmed, Z. et al. Causes, characterization, damage models, and constitutive modes for rock damage analysis: A review. Arab. J. Geosci. 13(16), 1–14. 10.1007/s12517-020-05755-3 (2020).10.1007/s12517-020-05755-3
16. Cho JW Kim H Jeon S Deformation and strength anisotropy of Asan gneiss, Boryeong shale, and Yeoncheon schist Int. J. Rock Mech. Min. Sci. 2012 50 158 169 10.1016/j.ijrmms.2011.12.004
Cho, J. W. et al. Deformation and strength anisotropy of Asan gneiss, Boryeong shale, and Yeoncheon schist. Int. J. Rock Mech. Min. Sci. 50, 158–169. 10.1016/j.ijrmms.2011.12.004 (2012).10.1016/j.ijrmms.2011.12.004
17. Alejano LR González-Fernández MA Estévez-Ventosa X Anisotropic deformability and strength of slate from NW-Spain Int. J. Rock Mech. Min. Sci. 2021 148 104923 10.1016/j.ijrmms.2021.104923
Alejano, L. R. et al. Anisotropic deformability and strength of slate from NW-Spain. Int. J. Rock Mech. Min. Sci. 148, 104923. 10.1016/j.ijrmms.2021.104923 (2021).10.1016/j.ijrmms.2021.104923
18. Dambly M Nejati M Vogler D On the direct measurement of shear moduli in transversely isotropic rocks using the uniaxial compression test Int. J. Rock Mech. Min. Sci. 2019 113 220 240 10.1016/j.ijrmms.2018.10.025
Dambly, M. et al. On the direct measurement of shear moduli in transversely isotropic rocks using the uniaxial compression test. Int. J. Rock Mech. Min. Sci. 113, 220–240. 10.1016/j.ijrmms.2018.10.025 (2019).10.1016/j.ijrmms.2018.10.025
19. Gholami R Rasouli V Mechanical and elastic properties of transversely isotropic slate Rock Mech. Rock Eng. 2014 47 5 1763 1773 10.1007/s00603-013-0488-2
Gholami, R. & Rasouli, V. Mechanical and elastic properties of transversely isotropic slate. Rock Mech. Rock Eng. 47(5), 1763–1773. 10.1007/s00603-013-0488-2 (2014).10.1007/s00603-013-0488-2
20. Ding CD Hu DW Zhou H Investigations of P-wave velocity, mechanical behavior and thermal properties of anisotropic slate Int. J. Rock Mech. Min. Sci. 2020 127 104176 10.1016/j.ijrmms.2019.104176
Ding, C. D. et al. Investigations of P-wave velocity, mechanical behavior and thermal properties of anisotropic slate. Int. J. Rock Mech. Min. Sci. 127, 104176. 10.1016/j.ijrmms.2019.104176 (2020).10.1016/j.ijrmms.2019.104176
21. Xu GW He C Chen ZQ Transversely isotropic creep behavior of phyllite and its influence on the long-term safety of the secondary lining of tunnels Eng. Geol. 2020 278 105834 10.1016/j.enggeo.2020.105834
Xu, G. W. et al. Transversely isotropic creep behavior of phyllite and its influence on the long-term safety of the secondary lining of tunnels. Eng. Geol. 278, 105834. 10.1016/j.enggeo.2020.105834 (2020).10.1016/j.enggeo.2020.105834
22. Xu GW He C Su A Experimental investigation of the anisotropic mechanical behavior of phyllite under triaxial compression Int. J. Rock Mech. Min. Sci. 2018 104 100 112 10.1016/j.ijrmms.2018.02.017
Xu, G. W. et al. Experimental investigation of the anisotropic mechanical behavior of phyllite under triaxial compression. Int. J. Rock Mech. Min. Sci. 104, 100–112. 10.1016/j.ijrmms.2018.02.017 (2018).10.1016/j.ijrmms.2018.02.017
23. Deng J Gu DS On a statistical damage constitutive model for rock materials Comput. Geosci. 2011 37 2 122 128 10.1016/j.cageo.2010.05.018
Deng, J. & Gu, D. S. On a statistical damage constitutive model for rock materials. Comput. Geosci. 37(2), 122–128. 10.1016/j.cageo.2010.05.018 (2011).10.1016/j.cageo.2010.05.018
24. Li X Cao WG Su YH A statistical damage constitutive model for softening behavior of rocks Eng. Geol. 2012 143–144 1 17 10.1016/j.enggeo.2012.05.005
Li, X., Cao, W. G. & Su, Y. H. A statistical damage constitutive model for softening behavior of rocks. Eng. Geol. 143–144, 1–17. 10.1016/j.enggeo.2012.05.005 (2012).10.1016/j.enggeo.2012.05.005
25. Liu HY Lv SR Zhang LM A dynamic damage constitutive model for a rock mass with persistent joints Int. J. Rock Mech. Min. Sci. 2015 75 132 139 10.1016/j.ijrmms.2015.01.013
Liu, H. Y. et al. A dynamic damage constitutive model for a rock mass with persistent joints. Int. J. Rock Mech. Min. Sci. 75, 132–139. 10.1016/j.ijrmms.2015.01.013 (2015).10.1016/j.ijrmms.2015.01.013
26. Chen HG Zhao C Multiscale damage constitutive model of jointed rock masses based on the influence of crack initiation IOP Conf. Ser. Earth Environ. Sci. 2021 861 2 022006 10.1088/1755-1315/861/2/022006
Chen, H. G. & Zhao, C. Multiscale damage constitutive model of jointed rock masses based on the influence of crack initiation. IOP Conf. Ser. Earth Environ. Sci. 861(2), 022006. 10.1088/1755-1315/861/2/022006 (2021).10.1088/1755-1315/861/2/022006
27. Zhao H Zhang C Cao WG Statistical meso-damage model for quasi-brittle rocks to account for damage tolerance principle Environ. Earth Sci. 2016 10.1007/s12665-016-5681-7
Zhao, H. et al. Statistical meso-damage model for quasi-brittle rocks to account for damage tolerance principle. Environ. Earth Sci.10.1007/s12665-016-5681-7 (2016).10.1007/s12665-016-5681-7
28. Chen YL Xiao P Du X Study on damage statistical constitutive model of triaxial compression of acid-etched rock under coupling effect of temperature and confining pressure Materials 2021 14 23 7414 10.3390/ma14237414 34885571
Chen, Y. L. et al. Study on damage statistical constitutive model of triaxial compression of acid-etched rock under coupling effect of temperature and confining pressure. Materials 14(23), 7414. 10.3390/ma14237414 (2021).34885571 10.3390/ma14237414
29. Chen YF Lin H Wang YX Statistical damage constitutive model based on the Hoek–Brown criterion Arch. Civ. Mech. Eng. 2021 21 3 1 9 10.1007/s43452-021-00270-y
Chen, Y. F. et al. Statistical damage constitutive model based on the Hoek–Brown criterion. Arch. Civ. Mech. Eng. 21(3), 1–9. 10.1007/s43452-021-00270-y (2021).10.1007/s43452-021-00270-y
30. Zhou CT Karakus M Xu CS A new damage model accounting the effect of joint orientation for the jointed rock mass Arab. J. Geosci. 2020 13 7 1 13 10.1007/s12517-020-5274-3
Zhou, C. T. et al. A new damage model accounting the effect of joint orientation for the jointed rock mass. Arab. J. Geosci. 13(7), 1–13. 10.1007/s12517-020-5274-3 (2020).10.1007/s12517-020-5274-3
31. Fu HL Zhang JB Huang Z A statistical model for predicting the triaxial compressive strength of transversely isotropic rocks subjected to freeze-thaw cycling Cold Reg. Sci. Technol. 2018 10.1016/j.coldregions.2017.11.003
Fu, H. L. et al. A statistical model for predicting the triaxial compressive strength of transversely isotropic rocks subjected to freeze-thaw cycling. Cold Reg. Sci. Technol.10.1016/j.coldregions.2017.11.003 (2018).10.1016/j.coldregions.2017.11.003
32. Li HZ Liao HJ Xiong GD A three-dimensional statistical damage constitutive model for geomaterials J. Mech. Sci. Technol. 2015 29 1 71 77 10.1007/s12206-014-1211-2
Li, H. Z. et al. A three-dimensional statistical damage constitutive model for geomaterials. J. Mech. Sci. Technol. 29(1), 71–77. 10.1007/s12206-014-1211-2 (2015).10.1007/s12206-014-1211-2
33. Liu YS Qiu ZZ Zhan XC Study of statistical damage constitutive model of layered composite rock under triaxial compression Appl. Math. Nonlinear Sci. 2021 6 2 299 308 10.2478/amns.2021.2.00048
Liu, Y. S. et al. Study of statistical damage constitutive model of layered composite rock under triaxial compression. Appl. Math. Nonlinear Sci. 6(2), 299–308. 10.2478/amns.2021.2.00048 (2021).10.2478/amns.2021.2.00048
34. Liu H Zhang L A damage constitutive model for rock mass with nonpersistently closed joints under uniaxial compression Arab. J. Sci. Eng. 2015 10.1007/s13369-015-1777-8
Liu, H. & Zhang, L. A damage constitutive model for rock mass with nonpersistently closed joints under uniaxial compression. Arab. J. Sci. Eng.10.1007/s13369-015-1777-8 (2015).10.1007/s13369-015-1777-8
35. Liu HY Yuan XP A compressive damage constitutive model for rock mass with a set of nonpersistently closed joints under biaxial conditions Math. Probl. Eng. 2015 2015 136971 10.1155/2015/136971
Liu, H. Y. & Yuan, X. P. A compressive damage constitutive model for rock mass with a set of nonpersistently closed joints under biaxial conditions. Math. Probl. Eng. 2015, 136971. 10.1155/2015/136971 (2015).10.1155/2015/136971
36. Chen S Qiao CS Composite damage constitutive model of jointed rock mass considering crack propagation length and joint friction effect Arab. J. Geosci. 2018 11 11 283 10.1007/s12517-018-3643-y
Chen, S. & Qiao, C. S. Composite damage constitutive model of jointed rock mass considering crack propagation length and joint friction effect. Arab. J. Geosci. 11(11), 283. 10.1007/s12517-018-3643-y (2018).10.1007/s12517-018-3643-y
37. Wang TT Huang TH A constitutive model for the deformation of a rock mass containing sets of ubiquitous joints Int. J. Rock Mech. Min. Sci. 2009 46 3 521 530 10.1016/j.ijrmms.2008.09.011
Wang, T. T. & Huang, T. H. A constitutive model for the deformation of a rock mass containing sets of ubiquitous joints. Int. J. Rock Mech. Min. Sci. 46(3), 521–530. 10.1016/j.ijrmms.2008.09.011 (2009).10.1016/j.ijrmms.2008.09.011
38. Liu HY Yuan XP A damage constitutive model for rock mass with persistent joints considering joint shear strength Can. Geotech. J. 2015 52 8 3107 3117 10.1139/cgj-2014-0252
Liu, H. Y. & Yuan, X. P. A damage constitutive model for rock mass with persistent joints considering joint shear strength. Can. Geotech. J. 52(8), 3107–3117. 10.1139/cgj-2014-0252 (2015).10.1139/cgj-2014-0252
39. Liu HY Lv SR Zhang LM Dynamic damage constitutive model of penetrating jointed rock mass based on combined model method Geotech. Eng. J. 2014 36 10 1814 1821 10.11779/CJGE201410008
Liu, H. Y., Lv, S. R. & Zhang, L. M. Dynamic damage constitutive model of penetrating jointed rock mass based on combined model method. Geotech. Eng. J. 36(10), 1814–1821. 10.11779/CJGE201410008 (2014) (In Chinese).10.11779/CJGE201410008
40. Yang JM Zhang YH Li QW Dynamic constitutive model of penetrating jointed rock mass based on ZTW model IOP Conf. Ser. Earth Environ. Sci. 2019 10.1088/1755-1315/237/3/032110
Yang, J. M. et al. Dynamic constitutive model of penetrating jointed rock mass based on ZTW model. IOP Conf. Ser. Earth Environ. Sci.10.1088/1755-1315/237/3/032110 (2019).10.1088/1755-1315/237/3/032110
41. Zhao YQ Lv SR Zhang LM Constitutive model of triaxial compression damage of jointed rock mass based on deformation elements J. Cent. South Univ. Nat. Sci. Ed. 2015 46 3 6 10.11817/j.issn.1672-7207.2015.03.029
Zhao, Y. Q. et al. Constitutive model of triaxial compression damage of jointed rock mass based on deformation elements. J. Cent. South Univ. Nat. Sci. Ed. 46(3), 6. 10.11817/j.issn.1672-7207.2015.03.029 (2015) (In Chinese).10.11817/j.issn.1672-7207.2015.03.029
42. Ji DL Zhao H Wang L Analytical solution of stress in a transversely isotropic floor rock mass under distributed loading in an arbitrary direction Appl. Sci. 2021 11 10476 10.3390/app112110476
Ji, D. L. et al. Analytical solution of stress in a transversely isotropic floor rock mass under distributed loading in an arbitrary direction. Appl. Sci. 11, 10476. 10.3390/app112110476 (2021).10.3390/app112110476
43. Li Y Zhang H Chen M Strength criterion of rock mass considering the damage and effect of joint dip angle Sci. Rep. 2022 12 2601 10.1038/s41598-022-06317-1 35173196
Li, Y. et al. Strength criterion of rock mass considering the damage and effect of joint dip angle. Sci. Rep. 12, 2601. 10.1038/s41598-022-06317-1 (2022).35173196 10.1038/s41598-022-06317-1
44. Qiao LP Wang ZC Huang AD Alteration of mesoscopic properties and mechanical behavior of sandstone due to hydro-physical and hydro-chemical effects Rock Mech. Rock Eng. 2017 50 2 255 267 10.1007/s00603-016-1111-0
Qiao, L. P., Wang, Z. C. & Huang, A. D. Alteration of mesoscopic properties and mechanical behavior of sandstone due to hydro-physical and hydro-chemical effects. Rock Mech. Rock Eng. 50(2), 255–267. 10.1007/s00603-016-1111-0 (2017).10.1007/s00603-016-1111-0
45. Pan JL Li P Xi X Research progress on corrosion durability of anchorage structures in underground engineering J. Harbin Inst. Technol. 2019 51 9 1 13
Pan, J. L. et al. Research progress on corrosion durability of anchorage structures in underground engineering. J. Harbin Inst. Technol. 51(9), 1–13 (2019) (in Chinese).
46. Pan JL Ren FH Cai MF Effect of joint density on rockburst proneness of the elastic-brittle-plastic rock mass Shock Vib. 2021 10.1155/2021/5574325
Pan, J. L. et al. Effect of joint density on rockburst proneness of the elastic-brittle-plastic rock mass. Shock Vib.10.1155/2021/5574325 (2021).10.1155/2021/5574325
47. Jaeger JC Shear failure of anisotropic rocks Geol. Mag. 1960 97 1 65 79 10.1017/S0016756800061100
Jaeger, J. C. Shear failure of anisotropic rocks. Geol. Mag. 97(1), 65–79 (1960).10.1017/S0016756800061100
48. GB/T50266-2013 Standard for Test Methods of Engineering Rock Mass 2014 China Planning Press
GB/T50266-2013. Standard for Test Methods of Engineering Rock Mass (China Planning Press, 2014) (in Chinese).
