Abstract
In the current paper, the deformation behaviours of rocks during compression are studied by testing 10 groups of sandstone samples with different porosity characteristics. According to the energy theory, the rock material was divided into two parts: solid skeleton and voids. A statistical damage-based approach was adopted to establish a nonlinear statistical damage constitutive model. The validity of the statistical damage constitutive model is verified by the test data. The statistical damage constitutive model performs well in each stage of rock compression before failure. For different types of rocks, different confining pressures, and different water contents, the statistical damage constitutive model fits well. This model can be applied to most types of rocks and in most engineering environments.
1. Introduction
Previous research [1–7] has shown that compression of rock materials typically involves four stages: (a) initial compression stage, (b) linear elastic deformation stage, (c) yield stage, and (d) failure stage. The stress-strain behaviour of rock is significantly influenced by the size and number of voids [8–10]. During the compression process, the volume changes continuously, leading to variation in the deformation modulus [11, 12]. Therefore, the stress-strain relationship of rock is nonlinear during compression. A thorough understanding of the deformation behaviour of rock during compression is essential for establishing an accurate numerical model for rock engineering.
As an imperfect material, rock includes pores and cracks. The porosity of rock has a significant influence on its mechanical properties and permeability [13–15]. In essence, rock damage is caused by the expansion of cracks and the propagation of voids in the rock leading to a decrease in its strength. At present, the methods for researching rock damage mainly include (a) using meso-damage theory and computed tomography scanning and (b) using continuum damage mechanics from the macroscopic point of view [16–18].
Continuum damage mechanics has been increasingly applied to model the constitutive behaviours of brittle rock [16, 18–23]. Great progress has been made in the incorporation of statistical theory in the establishment of damage constitutive models [16, 18–23]. Many recent statistical damage constitutive models have assumed that the strength of microcells follows a distribution law [16, 17, 24–26]. The most commonly applied distributions describing the microunit intensity distribution of rock are log-normal distributions, normal distributions, and Weibull distributions [17]. Chen et al. [17] explored a statistical damage constitutive model of rock and verified and compared the model with others following different kinds of distributions. The results showed that the Weibull distribution is more reasonable. However, these models assume that the stress-strain relationship of each mesoscopic element obeys Hooke’s law before the load reaches the peak. Thus, these models can simulate only elastoplastic materials, which show a quasilinear elastic stage and a yield stage during the compression of rock and cannot effectively simulate the plastic deformation process in the initial compression stage, which is when the initial voids in the rock are compacted. At the initial compression stage, the initial voids in the rock are diminishing during loading, and plastic deformation occurs. Thus, the stress-strain curve is an upper concave curve that is nonlinear in shape.
Commercial software (such as ANSYS and ABAQUS) for the numerical simulation of geotechnical engineering applications can find nonlinear solutions, as they usually have nonlinear material models. However, as the materials in these models are commonly regarded as a set of units, springs, or beams bonded together, these models can simulate only elastoplastic materials. Fortunately, such software usually offers an interface where user-defined material models can be used, allowing us to construct a constitutive model more suitable for porous rock materials. Such a model would be significant for the numerical simulation of rock engineering applications. Several scholars [27, 28] have tried to import custom functions into numerical simulations.
In this paper, we aimed to propose a statistical damage constitutive model that can effectively simulate the plastic deformation process during the initial compression stage. To study the deformation behaviour of porous rock during compaction, 10 groups of sandstone samples with different porosity characteristics were prepared for uniaxial compression tests. According to the energy principle, the rock material was divided into two parts: solid skeleton and void. Then, it is assumed that the strength of microcells of the solid skeleton follows the Weibull distribution. In this way, we proposed a nonlinear statistical damage constitutive model of rock that considers the initial compression stage.
2. Experimental Test
2.1. Test Materials and Methods
To study the deformation behaviour of rock during compression, we selected 10 groups of sandstone samples with different porosity characteristics for uniaxial compression tests. These sandstones were strongly weathered with high porosity; thus, the initial compression stage was very obvious. The sandstone samples were collected from Longchang County, Sichuan Province, China, and consisted of quartz, chlorite, feldspar, mica, and other auxiliary minerals. The clastic particles were subangular, with particle sizes between 200 and 500 microns.
The porosity of this kind of sandstone varies widely, ranging from 12 to 28%. The sandstone material was processed into standard samples with a size of 50 × 100 mm. According to porosity, the samples were divided into 10 groups, with 12 samples in each group, as shown in Table 1. Through uniaxial compression testing, the stress-strain relationships of the samples were obtained. The test instrument is shown in Figure 1. Axial stress is applied through the hydropress, and the deformation of the rock sample is monitored with strain gauges. The loading rate is 0.5 MPa/s. As the water content is also an important factor affecting the mechanical properties of rocks [29], the samples for the uniaxial compression test were divided into two states: dry and saturated. The test results are shown in Table 1.

2.2. Analyses of Experimental Results
The experimental results show that the porosity of sandstone is closely related to its mechanical properties. Figure 2 shows that there is a linear relationship between porosity and uniaxial compressive strength, whether the samples are dry or saturated.

The failure processes and stress-strain curves of the sandstone samples with different porosities are very similar. The principle of energy can help us better understand the deformation and failure process of rock. Energy dissipation in the initiation and propagation of voids and cracks in rock during the compression process causes damage in rock [30–35]. According to previous research [33, 35], the total work done by the external force on the rock during compression iswhere is the total work done by the external force; , , and are the principal stresses; and , , and are the principal strains.
During compression, some of the energy is stored in the rock and converted into elastic strain energy, while some energy dissipates due to plastic deformation. The elastic strain energy accumulated in rocks can be calculated by equation (2). The dissipated energy can be calculated bywhere is the elastic strain energy, is the dissipated energy, is Young’s modulus, and is Poisson’s ratio.
The stress-strain curve is an upper concave curve, as shown in Figures 3(a) and 3(b), which shows the energy change curve during the deformation process. At the initial compression stage (OA), a certain number of initial voids are diminished. The rock accumulates little elastic strain energy, and most of the input energy is dissipated. At the elastic stage (AB), the energy absorbed by the rock is stored as the elastic strain energy and dissipates very little. At the yield stage (BC), elastic energy increases at a diminishing rate. Dissipated energy begins to increase as the rock begins to undergo plastic deformation and microcracks begin to form and expand. During this stage, the proportion of elastic strain energy to total energy decreases gradually. Overall, the deformation of the solid skeleton under the load is elastic deformation before reaching the yield stage. Thus, the rock sample can be considered to be composed of voids and the solid skeleton at the macroscopic scale, as shown in Figure 3. The initial length of the sample before loading is defined as , the length of the void is , and the length of the solid skeleton is . The deformation under the axial stress is ∆l. Thus, the deformation of the solid skeleton is , and the deformation of the voids is . Then, the strain of the porous rock can be expressed as (Figure 4).

(a)

(b)

The axial strain of the solid skeleton and the axial strain of the voids can be expressed as (5) and (6), respectively. The relationship between and is defined in (7). The deformation characteristics of the rock sample, solid skeleton, and voids are shown in Figure 5. When the load reaches the critical point A (in Figure 3(a)), compaction of the voids is basically completed. At this time, the axial strain of the voids reaches a maximum of , and then the axial strain will remain unchanged until the rock completely fails, as shown in Figure 5.

The proportional coefficient is defined as , and the coefficient represents the proportion of the void axial strain to the total axial strain. According to our test results, the value of decreases with the increase in axial strain, and the relationship between the axial strain and is roughly linear. Our test materials have different porosities and water contents. From the results, we can see that the variation law of the proportionality coefficient is consistent during compression, as shown in Figure 6,and is not affected by the properties of the rock and the water it contains.

(a)

(b)

(c)

(d)
3. The Establishment of the Nonlinear Statistical Damage Constitutive Model
3.1. The Principle of the Statistical Damage Constitutive Model
Previous studies [16, 26] have shown that using the Weibull distribution to describe the strength of the microscopic unit of rock is reasonable; thus, the statistical damage constitutive model we propose will follow the Weibull distribution. If the strength of mesoscopic elements is assumed to follow the Weibull distribution, then the probability density function can be written as [16, 26]where m is a shape parameter and is a scale parameter.
Let denote the number of elements that have failed and denote the total number of elements. Then, the damage variable can be defined as follows [26]:
The number of failed elements increases by when the stress level increases to . If the external load increases from zero to , then the total number of failed elements can be written as [26]
Substituting (10) into (9), the damage variable can be written as
Based on the continuum damage mechanics, the statistical damage constitutive model of rock can be expressed as [26]where is Young’s modulus of the undamaged rock material; is Poisson’s ratio; D is the damage variable, which takes a value ranging from 0 for the intact or undamaged state (D = 0) and 1 for the fully damaged state.
Assuming that (a) the strength of each mesoscopic element of the rock follows the Weibull distribution during the loading process and (b) the stress-strain relationship of each mesoscopic element obeys Hooke’s law before the load reaches the peak, the statistical damage constitutive model of rock can be expressed as
The complete stress-strain curve of the rock under compression can be simplified as shown in Figure 7. Assume that the confining pressure is 0. According to the geometric conditions shown in Figure 7 and by taking the derivative of (13), we can obtain

Then, the constitutive model of the rock under compression can be written as
Equation (16) is the constitutive model of the rock under compression, which has been widely adopted [17]. Previous studies [16, 17] have proved that this model can describe the stress-strain behaviour of rock during the process of compression. From Figure 8, we can see that this constitutive model fits well for hard rock with few voids. However, for rock with high porosity and low strength, this constitutive model does not fit the test data well. Some scholars [16, 17] have improved this model but have not considered the plastic strain of rock during the initial compression stage (OA). Thus, the coefficient of determination (R2) of the fitting line at the OA stage is not high enough if the rock is not hard. Therefore, to effectively simulate the plastic deformation process during the initial compression stage, the statistical damage constitutive model should be improved.

(a)

(b)

(c)

(d)
3.2. The Statistical Damage Constitutive Model considering Plastic Deformation at the Initial Compression Stage
In Section 2, we suggest that the rock is composed of voids and a solid skeleton based on energy theory. The proportional coefficient , which represents the ratio of void strain to total strain, was defined. From Figure 6 we can see the value of decreases with the increase in axial strain, and the relationship between the axial strain and is roughly linear. Furthermore, is close to 1 when the axial strain is very low. Therefore, we define that
Coefficient k can be calculated as
In the compression stage (section OA in Figure 3(a)), the relationship between and the axial strain can be written as
Thus, can be calculated by
When reaching critical point A (in Figure 3(a)), the deformation of the voids reaches a maximum of , and the deformation will remain unchanged. can be defined as
According to the evolution law of elastic strain energy, before reaching the yield stage, the deformation of the solid skeleton under the load is elastic. The following is assumed:
Based on the assumption above, the value of can be calculated using (23), and can be calculated using
Combining (19)–(24), the expression of at any time can be written as (25) and (26).
When the axial stress reaches the peak, the corresponding strain on the solid skeleton is
When the load reaches the critical point A (in Figure 3(a)), the deformation of the voids reaches a maximum of , and the deformation remains unchanged. Thus, the stress-strain curve of the sample conforms to the previous statistical damage constitutive model, as defined in (16). It is assumed that (a) the strength of each mesoscopic element of the solid skeleton follows the Weibull distribution during the loading process, and (b) the stress-strain relationship of each mesoscopic element obeys Hooke’s law before failure. By replacing with in (16), the nonlinear statistical damage constitutive model that considers the initial compression stage can be obtained as
The statistical damage constitutive model requires 6 parameters, all of which can be directly obtained by conventional mechanical tests. Therefore, there are no difficulties related to obtaining parameters.
4. Discussion
The statistical damage constitutive model (see (28)) we established considers the plastic deformation during the initial compression stage and can better simulate the process of rock compression than the previous constitutive model (see (16)). To verify the reasonability of this statistical damage constitutive model, we compared the calculated values with the test data. Table 2 shows the error analysis of the statistical damage constitutive model (see (28)) at each stage during compression. It can be seen that, in each stage, the statistical damage constitutive model has a high fitting degree. The statistical damage constitutive model written as (16) is an elastoplastic constitutive model that is widely adopted at present. As it does not consider the plastic deformation during the initial compression stage, the degree of fit is low. If the voids of the rock are not developed or only slightly developed, then the initial compression stage may not be as notable or may not occur. Even if this stage is ignored, the initial compression stage has little impact on the results. However, if the rock develops microcracks or voids, then the stress-strain curve will be an upper concave curve in the initial compression stage. Thus, the statistical damage constitutive model written as (16) cannot effectively simulate the plastic deformation process of porous rocks during the initial compression stage.
To verify the performance of the statistical damage constitutive model on the simulation of other types of rocks, we referenced some published test data [36–38] of different types of rocks, as shown in Figure 9. As shown in Figure 9, for sandstone samples with different porosities, the constitutive model fits well. Thus, the statistical damage constitutive model also performs well for saturated samples. Furthermore, it can be seen that the statistical damage constitutive model also performs well when used to simulate the compression processes of other types of rocks. From Figure 9, we can see that the degree of fit is good even if the confining pressure is considered.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)
In conclusion, the nonlinear statistical damage constitutive model established in this study performs well in each stage of rock compression process before failure. For different types of rocks, different confining pressures, and different water contents, the statistical damage constitutive model fits well. Thus, this model can be applied to most types of rocks and most engineering environments. Since the stress-strain behaviour of rocks in the postpeak stage is too complex, it is not considered in this paper. Therefore, the model proposed in this paper still has room for improvement.
5. Conclusions
In this paper, particular attention is paid to the deformation behaviour of rock in the initial compression stage. A relatively reasonable nonlinear statistical damage constitutive model could be developed with the aid of the statistical damage theory. The following points summarize the present study:(1)Based on energy theory, the rock material was divided into two parts: solid skeleton and void. From the test results, the variation law of strain in each part of the rock during compression was obtained, which has nothing to do with rock properties or test conditions.(2)Based on the theory of damage mechanics, a nonlinear statistical damage constitutive model following the Weibull distribution was developed. This model can simulate the deformation behaviour of rock in the initial compression stage, as the strain of the void and solid skeleton are calculated separately.(3)The validity of the nonlinear statistical damage constitutive model is verified by the existing test data. The statistical damage constitutive model established in this study performs well regardless of which stage of the rock compression process is simulated. For different types of rocks, different confining pressures, and different water contents, the statistical damage constitutive model fits well. Thus, this model can be applied to most types of rocks and in most engineering environments.(4)The statistical damage constitutive model established in this paper is for porous rocks. However, the results of contrastive analysis show that the statistical damage constitutive model also has a high fitting degree for rocks with high stiffness and dense structure, so the statistical damage constitutive model can also be suitable for such rocks.
Data Availability
The test data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was funded by the Natural Science Foundation of China (no. 41602293).