Research Article | Open Access

Volume 2019 |Article ID 3065414 | https://doi.org/10.1155/2019/3065414

Ying Chen, Lin Zhang, Hui Xie, Jianfeng Liu, Hong Liu, Baoquan Yang, "Damage Ratio Based on Statistical Damage Constitutive Model for Rock", Mathematical Problems in Engineering, vol. 2019, Article ID 3065414, 12 pages, 2019. https://doi.org/10.1155/2019/3065414

# Damage Ratio Based on Statistical Damage Constitutive Model for Rock

Revised08 Jun 2019
Accepted11 Jul 2019
Published06 Aug 2019

#### Abstract

The study of damage characteristics of rock mass is of great significance to the analysis of rock mass structure. According to the characteristics of the microscopic unit strength of rock with random distribution, the Weibull distribution is widely used as the statistical functions of the strength of the microunit of rock to establish the damage constitutive model. The concepts of damage ratio and damage index are proposed. Damage ratio is mainly used to describe the law of damage evolution in rock. Damage index can be used to evaluate the damage degree of rock. The influence of confining pressure on distribution parameters and damage ratio is analyzed through uniaxial and triaxial compression tests of sandstone. The results show that damage ratio is an index of structural characteristics of rock damage, which can reflect the evolution characteristics of microcracks in rock under spatial stress. Critical damage ratio refers to the damage ratio corresponding to the peak stress of rock and can be used as a parameter to characterize the strength of rock for corresponding to the peak strain one to one. The critical damage ratio is linearly related to the logarithmic function of confining pressure. Its relationship is as follows: . With the increase of , the increasing trend of slows down and gradually tends to a certain value. The larger the damage index is, the more serious the damage of rock is. The smaller the damage index is, the less serious the damage of rock is. Therefore, the damage index can be used to evaluate the damage degree of rock. It will be an important direction of rock damage mechanics research to distinguish the severity of rock damage by using damage index as the limit value.

#### 1. Introduction

Damage research of rock materials has been an important issue in geotechnical engineering. The introduction of damage mechanics theory and damage variables provides a basis for the study of rock damage. At present, the analysis and testing of rock generally assume that it is a homogeneous and isotropic material though there are many subtle and complex defects such as microcracks and cracks inside the rock which have a great influence on the basic mechanical parameters of the rock.

The study of the relationship between rock strength and deformation by using the theory of damage mechanics can accurately reflect the influence of defects, such as microcracks on the mechanical properties of rock, and is considered as one of the most effective methods for studying brittle material, such as rock . In addition, the main difference between rock and rock mass is that rock mass contains discontinuous structural planes such as fissures, joints, bedding, and faults. Rock mass can be regarded as a rock with complex structural planes. And the research on rock damage can promote the development of rock mechanics to a certain extent.

At present, the methods of studying rock damage are mainly as follows. (1) To study the damage of rocks by using the continuum damage mechanics from the macroscopic point of view (e.g., Krajcinovic et al. [2, 3], Marigo et al. , Frantziskonis and Desai , Kawamoto , Cheng and Duseault , and Aubertin , et al.) In these studies, the damage mechanics was introduced into rock mechanics and the rock damage was assumed to be isotropic. The damage mechanism was studied based on the structural characteristics of rocks and the corresponding models were established mainly based on the damage research method of metal materials. Therefore, the model established by this method is quite different from the actual situation of rock materials. (2) To deal with the problem of damage using mesodamage theory from the microscopic point of view (e.g., Xia et al. , Bai et al. , and Horii et al. ). In these studies, they believed that it was a good attempt to combine the physical basis of mesomechanics model with macroscopic continuity theory, for the development of microcracks could be described by the theory of mesostructure and was the only phenomenon on the macroscale. (3) The damage theory of rock is studied based on the microscopic phenomenological theory and the damage variable is assumed to be a kind of distribution (e.g., Krajcinovic et al. , Cao et al. , Wang et al. , Deng and Gu , Zhao et al. , Xu et al. [20, 21], Wang et al. , Chen et al. , and Bian et al. ). Statistical damage theory was introduced in these studies. They quantified the degree of damage in rock by microelement strength and established the statistical damage evolution equation and damage softening constitutive model according to the characteristics of random distribution of damage in rock, which made a breakthrough in the study of rock softening constitutive.

A large number of rock damage constitutive models have been obtained based on the methods above. These studies describe the mathematical and mechanical models of rock damage from different angles, which promotes the development of rock mechanics experiments and theories. The constitutive model obtained by the third method above is always called a statistical damage constitutive model. According to the existing published results, the research on damage variable and distribution parameters of statistical damage constitutive model of rock is one of the hotspots.

However, the damage degree of different rocks under specific conditions has not yet been recognized as a quantitative index, and the understanding on damage is different. Though damage variable based on statistical damage constitutive model of rock can be used to describe the degree of damage development, it is process parameter, which will change gradually with the process of rock damage and it is difficult to be quantified.

Based on this, the Weibull distribution commonly used by the predecessors is taken as the statistical functions of the strength of the microunit of rock to establish the damage constitutive model. The damage ratio is defined as a new term based on the damage area to describe the damage degree of rock. The damage ratio is defined in rock analogous to the usual formulation of void ratio in soil mechanics. A curve of relationship between damage ratio and confining pressure (comparable to the compression curve or e-p curve in soil mechanics) is fitted and the slope of the curve is used to evaluate the damage degree of rock and is defined to be damage index. The damage index can be used as a quantitative index to define the damage degree of rock under test conditions. The research results are of great significance for evaluating the damage degree of rock. To distinguish the damage severity of rock by using damage index as the limit value may be an important direction of rock damage mechanics.

#### 2. Statistical Damage Constitutive Model for Rock

The strength of the rock is also random due to the fact that there are a large number of microcracks, cracks, and other defects; their distribution law, size, and penetration length of them show obvious randomness, and these randomly distributed defects affect the strength of the rock. Therefore, microunit strength can be regarded as a random variable and its probability distribution can be studied.

##### 2.1. Weibull Distribution and Damage Variable

The Weibull distribution is a distribution function proposed by Swedish scholar Waloddi Weibull to describe the distribution of fracture stress in 1939. It is feasible to use Weibull distribution function to calculate the strength of brittle material. Rock is a kind of brittle material; its strength shows a large dispersion; therefore the Weibull distribution function can be used to statistically process the rock strength. Taking the microunit strength parameter as the mechanical symbol quantity, the density function expression of the Weibull distribution is where and F0 are distribution parameters of the Weibull distribution, which reflect the mechanical properties of brittle material, is the strength of the microunit, and is the probability of microunit destruction when the intensity is , showing the damage rate of microunit in the brittle material like rock during loading.

The definition of damage variable is the core problem to study the damage evolution law of rock and other materials under loading. According to continuum damage mechanics, the damage variable can be defined as the ratio of the damage area to the total area of the material when it is not damaged; that is,where is the Drucker-Prager failure criterion parameter, which is the parameter related to cohesion and internal friction angle , is the first invariant of stress tensor, and is the second variable of stress deviation; they are as follows:where is the effective stress corresponding to the nominal stress , using the definition of damage:

##### 2.2. Definition of Damage Ratio

The damage area of the microunit of the material is , and the area with no damaged of the material microunit is the total area minus the damage area; that is,

Similar to the definition of the void ratio in geotechnical mechanics, the ratio of the damage area to the area with no damaged of the rock microunit is defined as the damage ratio ; that is,

The damage ratio is an indicator of the structural characteristics of rock damage. The damage ratio can reflect the evolution characteristics of microcrack inside the rock under the action of spatial stress and can be used to evaluate the compressive strength and the degree of damage of the rock. Therefore, the study of the damage ratio is also of great significance:

In the process of rock test, internal damage accumulates gradually with the increase of load. When the load stops to increase (that means the stress peak is reached), the internal damage of the rock accumulates to a certain value, and the corresponding damage ratio is defined as the critical damage ratio. Its expression iswhere is the critical damage ratio and is the value corresponding to the peak stress of the rock.

##### 2.3. Establishment of Statistical Damage Constitutive Model

In the triaxial compression experiment, it is assumed that the rock microunit conforms to the generalized Hooke’s law of an isotropic elastomer:where and are the principal strain and stress of the microunit, is the elastic modulus, and is Poisson’s ratio.

The damage ratio is introduced to analyse the effect of confining pressure on rock, for the random distribution of the microcracks inside the rock changes under the influence of confining pressure.

In continuum damage mechanics, the three-dimensional isotropic linear elastic damage constitutive model can be expressed as  where is the nominal stress vector, is the material elastic matrix of the material, is the strain vector, and is the damage variable.

We can get the damage variable according to formula (11):when substituting formula (16) into formula (15), we can get

The statistical damage softening constitutive model of rock is established based on the obedience of microunit strength to Weibull distribution formula (1) in . After introducing damage ratio, the statistical damage constitutive model of rock is as follows:when substituting formula (18) into formula (12), we can get

The formula above is a statistical damage constitutive model of rock based on Weibull distribution. The key of the model is to get the two distribution parameters and F0 in the formula. The common method is to determine the values by curve fitting according to the rock triaxial test curves of rock as the fitting method is complicated and random. According to the method in , the distribution parameters are determined according to the stress strain peaks.

##### 2.4. Determination of Distribution Parameters

According to Drucker-Prager failure criterion,where is the macroscopic strength of the material, that is nominal strength, and are the Drucker-Prager parameters of failure criterion, is the first invariant of the stress tensor, and is the second variable of the stress offset, and the expression of them is as follows:where is the nominal stress of the material, that is, the actually measured stress.

According to the damage definition, formula (8) can be converted as follows:therefore,

In the conventional triaxial test, the confining pressure is as σ2 = σ3, and the strain is as ε2 = ε3; then formula (19) can be simplified aswhen the rock reaches the peak failure under the triaxial stress, there are σ1 =σc, ε1 =εc, where σc is the peak stress, and εc is the strain corresponding to the peak stress, as follows:

The key to solve the distribution parameters is to find the relation of . To assume the σ1 and σ3 as the function of ε1 and ε3, we can get the following:

The following expression can be got after differentiating the two equations of formula (25):

The following expression can be got after differentiating the equation of formula (28):

According to the hypothesis of , and F0 are functions of confining pressure; thenwhen substituting formulas (31) and (32) into (30), we can get the following:the formula above can be simplified aswhen dσ3 is eliminated, the formula above can be simplified aswe can get the following expression with formula (29):

From formula (36), we can get that or . If , then

If there is no dσ3 in formula (37) and it is a contradiction for the above formula is to differentiate σ3. However, it can only be ; then

It can be known from formulas (25) and (30)

Combining formulas (38), (39), (40), and (41), we can get the following:when substituting formula (42) into (27), we can get the following:

The formula above is the expression of the distribution parameter determined by the conventional triaxial test. When it is determined by uniaxial compression test, for σ2=σ3=0, the expression of the parameters determined by the uniaxial compression test is as follows:

In formula, ES is secant modulus of elasticity, i.e., the slope of the straight line connecting the origin and the peak stress point in the stress-strain curve, σc is the peak stress in uniaxial compression test, and εc is the strain corresponding to the peak stress.

#### 3. Uniaxial and Triaxial Tests for Rock

##### 3.1. Test Equipment

The test was carried out on the MTS815 concrete and rock mechanics experiment system of the College of Water Resources and Hydropower, Sichuan University (Figure 1). This equipment is mainly used for testing the mechanical properties of rock and concrete materials. It can collect both high and low speed data, track the whole process of rock failure, and obtain the stress-strain curve of the whole process of rock. It is one of the most complete rock mechanics testing equipment with the highest technical level in the world. The MTS815 experiment system is an all-digital computer automatic control system that collects load, stress, strain, and displacement values in real time. During the loading process, it can control the loading by means of axial force control, displacement control, and circumferential deformation control and can carry out the whole process test of uniaxial compression and conventional triaxial test, etc. The load-displacement, load-axial deformation, and load-circumferential deformation curves can be drawn in real time.

##### 3.2. Test Scheme

The specimens of fine sandstone were used in the test, and they were processed into Φ50 mm×100 mm standard cylindrical specimens. In order to verify the statistical softening constitutive model of rock and explore the relationship between damage ratio and confining pressure, uniaxial and triaxial compression tests on the specimens were carried out, respectively. In triaxial tests, confining pressures were set as 3 MPa, 6 MPa, 9 MPa, 12 MPa, and 24 MPa. To get the general law of damage ratio and confining pressure, the first four groups of confining pressures were set to increase by 3 MPa. And the last group of confining pressure was set large enough to verify the applicability of the law. A total of 18 specimens were tested for 3 specimens were tested for each confining pressure.

##### 3.3. Test Process

The main process of conventional triaxial compression test: (1) Install the prepared specimen on the test bench. (2) Install the axial and circumferential extensometers to measure the deformation of the specimens. (3) Manually apply the initial axial force of 2 kN to the specimen by means of displacement control to make the specimen in full contact with the indenter. (4) Close the triaxial chamber and the confining pressure to be loaded to the target value at a rate of 3 MPa/min by confining pressure loading system. (5) The axial force is controlled by means of axial force control to be loaded at the rate of 20 kN/min until the stress state of the specimen reaches the limit of proportionality; change the axial force to be loaded at a rate of 0.04 mm/min by the means of circumferential deformation control. (6) Stop the test when the axial force does not decrease significantly with the increase of deformation (Figure 3).

The process of uniaxial compression test: The axial force is controlled by means of force to be loaded at the rate of 5 kN/min. In order to ensure the safety of the test equipment, when the displacement reaches 0.28 mm, change the axial force to be loaded by means of displacement control until the specimen was destroyed (Figure 2).

In the whole test process, the data of axial pressure, longitudinal displacement, and transverse displacement in the stress-strain process are collected and recorded automatically by the computer according to the time interval of 3 seconds, the longitudinal displacement of 0.001 mm, and the transverse displacement of 0.005 mm.

##### 3.4. Test Results

The stress strain curve of sandstone under uniaxial stress is shown in Figure 4. It can be seen that the rock quickly destroys after reaching the peak stress and rapidly falls to its residual strength. The brittle failure of rock is obvious under the uniaxial stress. The stress strain curve of sandstone under triaxial stress is shown in Figure 5. It can be seen that the sandstone exhibits different postpeak characteristics under different confining pressures and the brittle failure can be seen under low confining pressure. With the increase of confining pressure, the brittle failure of rock gradually changes into plastic failure. The greater the confining pressures are, the more obvious the plastic deformation is.

In triaxial compression test, the confining pressure is σ2=σ3, so it can be obtained by formula (14):we can get the following expression after simultaneous (45) and (46):

The elastic modulus and Poisson’s ratio of rock can be obtained by using formulas (47) and (48) combined with the data of the straight line before the peak of the stress-strain relationship curve.

The confining pressure can be considered to be 0 MPa in uniaxial compression test. The mechanical parameters of rock under different confining pressures are shown in Table 1. It shows that as the confining pressure increases, the elastic modulus, the peak stress, and the corresponding strain increase as well as the Poisson’s ratio. The confining pressure has a significant influence on the mechanical parameters and deformation characteristics of the rock.

 Confining σ1 ε1/10−3 E Es μ Pressure (MPa) (MPa) (GPa) (GPa) 0 31.34 5.125 7.556 6.115 0.291 3 52.19 12.783 5.525 4.083 0.292 6 75.38 16.573 6.702 4.548 0.301 9 92.40 17.787 8.052 5.195 0.311 12 102.50 19.539 8.351 5.246 0.318 24 151.61 24.780 10.124 6.118 0.324
##### 3.5. Solution of Model Parameters and Damage Ratio

The value can be obtained by substituting the test data in Table 1 into formula (43). and F0 values can be obtained by formulas (24) and (42). In order to obtain value, cohesion C and internal friction angle of rock should be calculated first according to test data.

At present, Mohr-Coulomb criterion is the most widely used strength theory in the study of rock mechanics, and its expression is as follows:

The Mohr-Coulomb criterion can be expressed in another way. The failure criteria of rocks can be expressed in terms of σ1 and σ3. The expression is as follows:

According to the test data, the optimal relation curve is drawn by using least square method with σ3 as abscissa and σ1 as ordinate, as shown in Figure 6. It can be seen that the linear relationship between the axial stress and the confining pressure of sandstone is as follows: σ1=4.8898σ3+40.229.

From the optimal relation curve of rock and formulas (50) and (51), the cohesion C=9.096 MPa and internal friction angle =41.33° can be obtained.

The distribution parameters F0 and m of sandstone under different confining pressures and the critical damage ratio Der of rock under peak stress can be obtained by incorporating C, , and test data into formulas (4), (5), (13), (24), (42), and (43), as shown in Table 2.

 σ3 (MPa) m F0 (MPa) Der 0 4.726 31.230 0.236 3 2.970 75.256 0.400 6 2.290 129.836 0.548 9 1.997 172.682 0.650 12 1.844 199.821 0.720 24 1.634 320.298 0.844
##### 3.6. Model Verification

Substituting the parameters and F0 with values calculated from tests in formula (19), the statistical damage constitutive model under different confining pressures can be obtained. To verify the model, with the experimental strain as the input value of the model, the stress obtained by using the constitutive model should be compared with the experimental results. Figure 7 shows the comparison of stress strain curves between the proposed model-based curve and the test curve under various confining pressures. It can be seen from the figure that the simulated curve can reflect the process of rock failure before peak. The difference of the postpeak curve between test curve and model curve is caused by the change of loading mode, loading speed, etc. during the test and the model needs to be modified as what has been studied in .

#### 4. Discussion

##### 4.1. Influence of Confining Pressure on Distribution Parameters

According to , is a parameter to characterize the brittleness of rock material, and F0 is a parameter to reflect the macroscopic average strength of rock. From Table 2, it can be found that decreases with the increase of confining pressure, indicating that brittleness of rock is weakening while plasticity is increasing. Moreover, the value subject to the triaxial stress is significantly reduced compared with that in uniaxial stress, indicating that the confining pressure has a significant influence on the failure characteristics of the rock, and even a small confining pressure can significantly improve the postpeak plasticity of the rock.

It has been shown in Table 2 that, as the confining pressure increases, F0 increases, indicating that the macroscopic average strength of rock increases with the increase of confining pressure, for the bearing capacity of the rock is enhanced with the change of the stress of the rock under confining pressure.

##### 4.2. Influence of Confining Pressure on Damage Ratio

Under triaxial stress condition, when the distribution parameters m and F0 are determined, the damage variables can be calculated according to formula (12), and a strain value corresponds to a damage ratio. Therefore, the whole damage ratio curve of sandstone can be obtained by introducing the strain value of the stress-strain curve of sandstone under triaxial compression into (12) (as shown in Figure 7). The curve can reflect the damage development process of the damage of the sandstone under different confining pressures.

From the graph, it can be found that the damage ratio of rock increases with the increase of longitudinal strain, which essentially reflects the process of internal damage expanding and rock failure from loading to gradual failure under the action of axial stress. Furthermore, the increase of longitudinal strain speeds up the increase of the damage ratio. Especially when the rock tends to be damaged, the curve tends to be vertical straight line, the damage ratio develops faster, and the magnitude of increase is more obvious. It can be shown that, when the rock approaches failure, the cumulative area of damage increases and the area of no damage decreases, while the damage ratio increases rapidly as the loading proceeds.

Meanwhile, it can be seen from Figure 8 that, when the damage ratio is 0.4, the tangent value of the curve decreases with the increase of confining pressure, which shows that the increase rate of damage ratio under high confining pressure is less than that under low confining pressure at the same damage ratio. Therefore, the increase of confining pressure can delay the increase of damage ratio, inhibit the development of damage, increase the strength of rock, slow down the failure process of rock, and make the transition from brittle failure to plastic failure.

##### 4.3. Influence of Confining Pressure on Critical Damage Ratio

Critical damage ratio refers to the damage corresponding to the peak stress of rock for it corresponds to peak strain one to one. To be similar to the strength parameters of rock, the critical damage ratio can be used as a parameter to indicate whether rock is damaged or not. The change curve of critical damage ratio and confining pressure is as shown in Figure 9. It can be found that, with the increase of confining pressure, the critical damage ratio increases gradually in curve and the increase range slows down.

In order to further determine the relationship between critical damage ratio and confining pressure, it is found that the relationship between critical damage ratio and logarithmic function of confining pressure is well linear by data fitting, as shown in Figure 10. In the figure, the Y axis is Der and the X axis is lnσ3. Thus, the expression of critical damage ratio and confining pressure can be obtained as follows:

By introducing the data of , the relationship between Der and lnσ3 can be obtained through the methods mentioned above as shown in Figure 11. Therefore, the relationship between critical damage ratio and confining pressure can be obtained as follows:

It can be found that the relationship between critical damage ratio and confining pressure is as follows:

In the formula, and are constants, which can be obtained by linear fitting of experimental data.

The following expression can be getting after differentiating the equation of formula (54):

Formula (55) reflects the tangent of the Der-σ3 curve. It can be seen from formula (55) that when σ3 is infinite, (55) equals zero, which means that, with the increase of σ3, Der tends to a certain value which is related to the category of rock material instead of increasing infinitely.

##### 4.4. Definition and Physical Significance of Damage Index

The Der-lnσ3 curve is obtained from the experimental data in this paper and the data in the literature , as shown in Figure 12. The results of data fitting are shown in Table 3.

 Number Correlation coefficient (R2) 1 0.2170 0.1657 0.9959 2  0.0895 0.2069 0.9455 3  0.2740 0.3913 0.9305 4  0.0761 1.1669 0.9829 5  0.0389 0.0297 0.9965

From the Der-lnσ3 curves shown in Figure 12, we can get the following rule: the curve is close to a straight line. We can define the slope of the curves (the damage index Cc) as

According to slope of test curve and reference curve, it can be seen that the linear slope of rock in different categories is not the same. The larger the value of , the greater the influence of confining pressure on damage ratio.

The existence of confining pressure also changes the failure mechanism and the condition of the microcracks in the rock, which leads to the various closing of microcracks or eases the expanding of fractures.

The microcracks can be regarded as internal damage of the rock. Essentially, the effect of confining pressure is to restrain the development of damage ratio in rock; that is the damage area diffusivity of rock under confining pressure is slower than that without confining pressure, and the damage area diffusivity under high confining pressure is lower than that under low confining pressure. Therefore, greater axial force was required for rock failure when confining pressure exists.

It can also be considered that the more microcracks are distributed in the rock, the more serious the rock damage is, and the more obvious the inhibiting effect of confining pressure on the damage development is, for the damage degree can be greatly changed by only a small confining pressure difference. On the contrary, if the rock interior is complete and the damage is slight, the effect of confining pressure on the damage development is not obvious.

Then, the damage index can reflect the effect of confining pressure on the damage ratio. Therefore, the damage index can be used as an index to evaluate the degree of internal damage of rock. The larger the damage index is, the more serious the damage of rock is. The smaller the damage index is, the slighter the damage degree of rock is. Similar to the compression index that can judge the compressibility of soil, the damage index can be used to judge the degree of rock damage.

Because of the limited test data and limited rock types in this paper, it is impossible to give a clear discriminant index, but we can try to predict the following classification:When Cc > M, rock is seriously damaged.When , rock is moderately damaged.When Cc<N, rock is slightly damaged.

Among them, M and N are the limit values of damage index for judging the damage degree of rock, which need to be determined by a large number of tests. A research idea has been provided in this paper, and the follow-up research can be carried out deeper on the basis of this.

#### 5. Conclusions

The study of rock damage characteristics is of great significance to the analysis of rock mass. In this paper, the main conclusions are as follows.

(1) Damage ratio is an index of damage structure characteristics of rock. It not only can reflect the evolution characteristics of microcracks in rock under spatial stress, but also can be used to evaluate the degree of rock damage and the compressive strength of rock. Based on the study of statistical damage constitutive model and damage variable of rock, the definition and solution of damage ratio are obtained.

(2) Critical damage ratio refers to the damage corresponding to the peak stress of rock. For corresponding to the peak strain one to one, the critical damage ratio can be used as a parameter to characterize whether the rock is damaged or not. The critical damage ratio has a good linear relationship with the logarithmic function of confining pressure. Its relationship is as follows: . With the increase of σ3, the increasing trend of Der slows down and gradually tends to a certain value.

(3) The slope of Der-lnσ3 curve is defined as damage index. The larger the damage index, the more serious the damage of rock. The smaller the damage index, the slighter the damage of rock. The damage index can be used as a quantitative index to define the damage degree of rock under test conditions. The research results are of great significance for evaluating the damage degree of rock. To distinguish the damage severity of rock by using damage index as the limit value may be an important direction of rock damage mechanics.

#### Data Availability

(1) The data used to support the findings of this study are included within the Supplementary Information file. (Note that the data above is provided to the editor only to use during the review process. And it is not supported for sharing at this stage.) (2) Previously reported data were used to support this study and are available at CNKI. These prior studies are cited at relevant places within the text as .

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (no. 51609163), National key R & D projects with Grant no. 2016YFC0401908, and Sichuan Science and Technology Program (no. 2018JZ0001).

1. H. Y. Liu, G. W. Wang, and G. Z. Liu, “Research development of rock damage theory characterized by damage variable,” Explosive Materials, vol. 33, no. 6, pp. 25–29, 2004. View at: Google Scholar
2. D. Krajcinovic, M. Basista, and D. Sumarac, “Micromechanically inspired phenomenological damage model,” Journal of Applied Mechanics, vol. 58, no. 2, pp. 305–310, 1991. View at: Publisher Site | Google Scholar
3. D. Krajcinovic and S. Mastilovic, “Some fundamental issues of damage mechanics,” Mechanics of Materials, vol. 21, no. 3, pp. 217–230, 1995. View at: Publisher Site | Google Scholar
4. J. J. Marigo, “Modelling of brittle and fatigue damage for elastic materials by growth of microvoids,” Engineering Fracture Mechanics, vol. 21, no. 4, pp. 861–874, 1985. View at: Publisher Site | Google Scholar
5. G. Frantziskonis and C. S. Desai, “Constitutive model with strain softening,” International Journal of Solids and Structures, vol. 23, no. 6, pp. 733–768, 1987. View at: Publisher Site | Google Scholar
6. T. Kawamoto, Y. Ichikawa, and T. Kyoya, “Deformation and fracturing behaviour of discontinuous rock mass and damage mechanics theory,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 12, no. 1, pp. 1–30, 1988. View at: Publisher Site | Google Scholar
7. H. Cheng and M. B. Dusseault, “Deformation and diffusion behaviour in a solid experiencing damage: a continuous damage model and its numerical implementation,” International Journal of Rock Mechanics & Mining Science & Geomechanics Abstracts, vol. 30, no. 7, pp. 1323–1331, 1993. View at: Google Scholar
8. M. Aubertin, D. E. Gill, and B. Ladanyi, “Constitutive equations with internal state variables for the inelastic behaviour of soft rocks,” Applied Mechanics Reviews, vol. 47, no. 6s, pp. S97–S101, 1994. View at: Publisher Site | Google Scholar
9. M. F. Xia, W. S. Han, F. J. Ke et al., “Statistical meso-scopic damage mechanics and damage evolution induced catastrophe,” Advances in Mechanics, vol. 25, no. 2, pp. 145–173, 1995. View at: Google Scholar
10. Y. L. Bai, M. F. Xia, F. J. Ke et al., “Closed trans-scale statistical microdamage mechanics,” Acta Mechanica Sinica, vol. 18, no. 1, pp. 1–17, 2002. View at: Google Scholar
11. H. Horii and S. Nemat-Nasser, “Compression-induced microcrack growth in brittle solids: Axial splitting and shear failure,” Journal of Geophysical Research: Atmospheres, vol. 90, no. B4, pp. 3105–3125, 1985. View at: Publisher Site | Google Scholar
12. D. Krajcinovic and M. Silva, “Statistical aspects of the continuous damage theory,” International Journal of Solids and Structures, vol. 18, no. 7, pp. 551–562, 1982. View at: Google Scholar
13. W. G. Cao, Z. L. Fang, and X. J. Tang, “A study of statistical constitutive model for soft and damage rocks,” Chinese Journal of Rock Mechanics and Engineering, vol. 17, no. 6, pp. 628–633, 1998. View at: Google Scholar
14. W. G. Cao, M. H. Zhao, and X. J. Tang, “Study on simulation of statistical damage in the full process of rock failure,” Chinese Journal of Geotechnical Engineering, vol. 25, no. 2, pp. 184–187, 2003. View at: Google Scholar
15. W. Cao, X. Li, and H. Zhao, “Damage constitutive model for strain-softening rock based on normal distribution and its parameter determination,” Journal of Central South University of Technology, vol. 14, no. 5, pp. 719–724, 2007. View at: Publisher Site | Google Scholar
16. W. Cao, X. Tan, C. Zhang, and M. He, “Constitutive model to simulate full deformation and failure process for rocks considering initial compression and residual strength behaviors,” Canadian Geotechnical Journal, vol. 56, no. 5, pp. 649–661, 2019. View at: Publisher Site | Google Scholar
17. Z. L. Wang, Y. C. Li, and J. G. Wang, “A damage-softening statistical constitutive model considering rock residual strength,” Computers & Geosciences, vol. 33, pp. 1–9, 2007. View at: Publisher Site | Google Scholar
18. J. Deng and D. S. Gu, “On a statistical damage constitutive model for rock materials,” Computers Geo-Sciences, vol. 37, pp. 122–128, 2011. View at: Publisher Site | Google Scholar
19. H. Zhao, C. Shi, M. Asce, M. Zhao, and X. Li, “Statistical damage constitutive model for rocks considering residual strength,” International Journal of Geomechanics, vol. 16, no. 4, pp. 1–9, 2017. View at: Publisher Site | Google Scholar
20. X. Xu and M. Karakus, “A coupled thermo-mechanical damage model for granite,” International Journal of Rock Mechanics and Mining Sciences, vol. 103, pp. 195–204, 2018. View at: Publisher Site | Google Scholar
21. X. Xu, F. Gao, and Z. Zhang, “Thermo-mechanical coupling damage constitutive model of rock based on the Hoek–Brown strength criterion,” International Journal of Damage Mechanics, vol. 27, no. 8, pp. 1213–1230, 2017. View at: Publisher Site | Google Scholar
22. J. Wang, Z. Song, B. Zhao, X. Liu, J. Liu, and J. Lai, “A study on the mechanical behavior and statistical damage constitutive model of sandstone,” Arabian Journal for Science and Engineering, vol. 43, pp. 5179–5192, 2018. View at: Publisher Site | Google Scholar
23. S. Chen, C. Qiao, Q. Ye, and M. U. Khan, “Comparative study on three-dimensional statistical damage constitutive modified model of rock based on power function and weibull distribution,” Environmental Earth Sciences, vol. 77, no. 3, p. 108, 2018. View at: Publisher Site | Google Scholar
24. K. Bian, J. Liu, W. Zhang, X. Zheng, S. Ni, and Z. Liu, “Mechanical behavior and damage constitutive model of rock subjected to water-weakening effect and uniaxial loading,” Rock Mechanics and Rock Engineering, vol. 52, no. 1, pp. 97–106, 2019. View at: Publisher Site | Google Scholar
25. Z. H. Chen, Y. F. Fu, and C. A. Tang, “Confining pressure effect on acoustic emissions during rock failure,” Chinese Journal of Rock Mechanics and Engineering, vol. 16, no. 1, pp. 65–70, 1997. View at: Google Scholar
26. M. H. Yang, M. H. Zhao, and W. G. Cao, “Method for determining the parameters of statistical damage softening constitutive model for rock,” Journal of Hydraulic Engineering, vol. 36, no. 3, pp. 345–349, 2005 (Chinese). View at: Google Scholar
27. Y. M. Kang, C. W. Liu, Y. Jia et al., “Research on statistical damage constitutive model and critical damage for rock subjected to triaxial stress contidion,” Journal of Sichuan University (Engineering Science Edition, vol. 41, no. 4, pp. 42–47, 2009. View at: Google Scholar
28. X. X. Han, X. B. Yang, and E. L. Liu, “Research on the statistical damage constitutive model of sandstone under triaxial compression,” Journal of China University of Mining & Technology, vol. 3, no. 2, pp. 165–170, 2018. View at: Google Scholar
29. R. L. Cao, S. H. He, J. Wei et al., “Study of modified statistical damage softening constitutive model for rock considering residual strength,” Rock and Soil Mechanics, vol. 34, no. 6, pp. 1652–1660, 2013. View at: Google Scholar

#### More related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.