The precondition of rock stress and deformation analysis is a reasonable rock constitutive model. Most of the previous studies have described the heterogeneous microdamage by Weibull distribution or normal distribution. However, both of them have limitations. Therefore, this paper intends to use the lognormal distribution as the probability distribution model of rock microunit strength. Based on the tensile failure of the single-fractured rock under the hydrodynamic force, the maximum tensile strain failure criterion is used as the distribution parameter of rock microunit strength. And, considering the multiphase properties of the filling fractured rock, the equivalent elastic modulus parameter is adopted in the model. We design a triaxial seepage test for the filled single-fractured rock and analyze the applicability and rationality of the modified lognormal statistical damage model for characterizing the fractured rock by using the test data. According to the comparison of the experimental stress-strain curve and the model stress-strain curve and the analysis of the damage value test curve and the model curve, the rationality of the established statistical damage constitutive model is verified, and the advantages and limitations of the model are proposed.

1. Introduction

The precondition of rock stress and deformation analysis is a reasonable rock constitutive model, which is the theoretical basis for stability prediction and evaluation of the surrounding rock in geotechnical engineering and also the key to the development of rock mechanic science. Dougill [1] studied the damage evolution and applied the damage theory to rock mechanics in 1976. Its main content is to study the damage evolution law of damaged materials and the damage theory. At present, there are two main research methods of the damage constitutive model reflecting the whole process of rock deformation and failure. One is based on the phenomenological statistical damage simulation method [27] (as shown in Table 1) and the other is based on CT, scanning electron microscopy, optical microscopy, laser holographic interferometry, and numerical modeling methods [814], which describe the damage situation of materials by quantitatively describing the damage structure information on the surface or inside of materials. Most of the previous studies have regarded the macrofailure phenomenon as the average effect of microunit destruction and described the heterogeneous microdamage by Weibull distribution or normal distribution, so as to establish a statistical damage model of rock under uniaxial or triaxial compression, describing the damage parameters of the rock and the constitutive relationship curve of the failure process.

However, for brittle materials such as the rock, the Weibull theory cannot reflect the nonlocal characteristics of the damage, and the application of the Weibull theory to brittle materials will inevitably lead to size effect, resulting in a large difference between the theoretical results and the experimental results [15]. So far, researchers still pay more attention to the continued research and improvement of this classic theory. At the same time, considering the variables of normal distribution defined on (,), if a normal distribution is used, it means that the model parameters may have negative values, while the material properties can only take positive values, which is not consistent with the physical meaning of the parameters themselves, so strictly speaking, it is not suitable for describing the material properties [16].

Based on the investigation of the factors affecting the strength of rock material RVE (representative volume element), Zhang et al. [17] propose that the lognormal distribution has a strict probability theorem as the theoretical basis, which may better reflect the properties of rock materials in real structures.

Therefore, on the basis of previous studies, this paper intends to use the lognormal distribution as the probability distribution model of rock microunit strength. According to the failure characteristics of the cracked rock in the seepage-stress field, the strength of the rock microunit is determined by the maximum tensile strain failure criterion. At the same time, considering the multiphase properties of the filling fractured rock, the equivalent elastic modulus parameter is adopted in the model. Finally, for the influence of confining pressure and water pressure, Hooke’s law and effective stress principle are used to modify. The statistical constitutive model of the lognormal distribution of seepage damage in the cracked rock is established. In this paper, the advantages and limitations of the model are obtained by comparing the experimental stress-strain curve with the model stress-strain curve. Based on the analysis of the damage variable test curve and the model curve, the rationality of the established statistical damage constitutive model is verified.

2. Rock Failure

2.1. Hydraulic Failure Mechanism of Single-Cracked Rock

As shown in Figure 1, the water pressure acts on the inner side of the crack and the surrounding rock stress acts on the far side of the fracture. The force on the crack boundary is (rock tensile strength). The failure problem of the single-cracked rock mass can be solved by the Dugdale model.

The rise of groundwater along the fault plane under pressure is called water wedge action. There are two manifestations. One is the splitting failure of rock cracks caused by water pressure. When the surrounding rock is in the state of compressive stress, the water pressure overcomes the surrounding rock stress and the rock tensile strength to induce the crack propagation. The other is the squeezing failure of rock cracks caused by water pressure. When the crack width in the rock mass is large or water-conducting fractures are formed, water pressure squeezes into the rock mass to promote the increase in the crack depth.

Therefore, water pressure enters the cracked rock mass, forming a large pore wall pressure, resulting in the formation of tensile cracks. The tensile strain rate in a certain direction first reaches an allowable value, and finally, the damaged rock in this environment contains more longitudinal tensile cracks.

2.2. Test Design and Results

Scholars carry out some studies on the mechanical properties of seepage in fractured rocks [18, 19], but there are few experiments on seepage in fractured rocks considering the filling material. In order to verify the failure mechanism and mode of cracked rock mass under water pressure, we have conducted the triaxial seepage test of single filled cracked rock.

The natural rock mass structure is complex, and the rock samples taken at the site are random and nonrepeatable. When studying the influence of cracks on the rock mass seepage damage, it is difficult to control the variables to obtain effective conclusions. At present, in the indoor seepage test of the filled fractured rock mass, the joints are prefabricated by cutting and splitting and then filled. In this study, the single crack is made in the intact sandstone by cutting. According to the common types of fillings in natural fractured rock masses recorded in the Engineering Geology Site Manual, cement mortar is used to fill the prefabricated cracks to complete the sample.

First, select the complete sandstone rock sample (). Then, a cutter is used to cut from one end of the cylinder to form a single slit with a width of 6 mm and a length of 25 mm. The single crack is grouted by a slurry prepared by cement, quartz sand (0.5–0.6), and water at a standard of 1 : 1.5 : 0.5. By repeating this process, the single filled cracked rock samples having a crack length of 50 mm and 75 mm are made (as shown in Figure 2).

In an indoor environment, the samples are dried for 18 days and then tested. The experiments are carried out on the MTS815 system. The axial loading is controlled by displacement. The confining pressure design values are 10 MPa and 20 MPa, and the water pressure is 7 MPa. The failure picture of the samples is shown in Figure 3.

It can be seen that under the seepage stress, the macroscopic failure mode of the single-cracked sandstone with filling material shows that the tensile cracks, which develop along the dynamic hydraulic direction, and the mixed tensile shear cracks at the distal end originate from the tensile stress concentration area near the tip of the fracture. The failure is the overall instability caused by the tensile shear mixing cracks that originate from the crack tip and extend to the upper and lower end faces of the rock sample.

3. Constitutive Model of Seepage Damage

3.1. Elastic Modulus

The study of rock structure models has been carried out for a long time, and the proposed models can be roughly divided into three categories [20]. Considering the filling crack characteristics of the test rock samples, the Voigt model is used to calculate the elastic parameters of the equivalent model, as seen in Figure 4.

The model proposed by Voigt assumes that the various minerals that make up the rock are arranged in parallel along the direction of the force. The bulk modulus of each species mineral is , the shear modulus is , and the percentage of the rock volume is (i = 1, 2, …, N). Using the volume averaging method, the multiphase equivalent rock bulk modulus of the Voigt model can be calculated as

According to the relation between the elastic parameters, the calculation equation can be written aswhere is the equivalent elastic modulus and is the Poisson ratio, ignoring its change with damage in this model. The calculation results are shown in Table 2.

3.2. Damage Variables and Constitutive Relations

One definition of the damage variable is a small face is taken inside the object and is evenly damaged by the P force. If the total cross-sectional area of the face is and the damaged area is , the net area of the surface is . Rabotnov [21] introduced a variable that corresponds to a continuous variable called the damage variable :when , it corresponds to an intact state; when , it corresponds to a completely damaged state; when , it corresponds to a different degree of the damage state.

Let be the nominal stress on the cross-sectional area. Let be the stress of the net section, called the net stress. According to the principle of strain equivalence proposed by Lemaitre [22], it is assumed that the strain caused by the net stress acting on the nondestructive material is equivalent to the strain caused by the stress acting on the damaged material. The constitutive relationship of damaged materials can be obtained by nominal stress in the nondestructive material, as follows:

3.3. Rock Microunit Strength

Assuming that the rock failure criterion iswhere is a constant. With the change of the stress state, if , it means that the rock microunit yields or breaks. It follows that can be used as a distribution variable of rock microunit strength random distribution because it reflects the dangerous degree of rock microunit failure.

Letwhere is a distribution variable of microunit strength random distribution and is the effective stress.

The form of the rock’s failure mechanism and its failure criterion entirely depends on the microunit strength form. In this paper, considering the failure mode of the cracked rock, the theory of maximum tensile strain, which is simple and widely used in the rock media, is used to establish the rock microunit strength , as follows:

The effective stress principle is generally used in the analysis of stress-seepage coupling problems. Terzaghi [23] first proposed an effective stress principle adapted to one-dimensional situation, and Biot [24] revised the effective stress principle and proposed a new one:where is the stress tensor, is the effective stress tensor, is the pore water pressure, is the unit second-order tensor and , otherwise , and is the Biot coefficient and determined by the following:where and are the bulk modulus of the solid particles and the bulk modulus of the skeleton in the rock mass. The results of the research on the Biot coefficient in the literature [2528] indicate that the value of ranges from 0 to 1. For most materials, the compressibility of the skeleton is much smaller than the one of the solid particles themselves, that is, , so it is convenient for research to take .

In cracked rock triaxial flow tests, assuming the rock matrix is impermeable and cracks are the seepage channel, when there is a seepage water pressure () in the crack, the nominal stress () and strain () can be measured, and the corresponding effective stress is .

According to Hooke’s law and equation (7), their relation is

Substituting the above ones into equation (7) to establish the rock’s microunit strength, we get

3.4. Modified Statistical Damage Constitutive Model

This paper uses a lognormal distribution to represent the rock microunit strength:where and are the lognormal distribution parameters.

The number of destroyed microunits under a certain load is and the total number of microunits is . The number of destroyed cells is in any interval . When loading to a level , the number of destroyed cells is

The damage variable can be expressed as

And, it is simplified aswhere and are the distribution parameters, is a standard normal distribution function, and , ().

From equations 10–(12), the following can be obtained:

The damage constitutive equation obtained by substituting (17) into (18) is

According to the triaxial test process, confining pressure and pore water pressure are loaded before deviatoric stress, so that the axial strain is generated, and the strain in this section is not included in the test curve. The axial bias recorded in the test is calculated aswhere is the axial stress and is the confining stress.

And, there is an initial strain before the axial stress is loaded, which can be calculated as follows:

So, the true axial strain iswhere is the measured test strain.

Substituting equations (20)–(22) into equation (19), the constitutive equation of and is

Similarly, by substituting equations (20)–(22), equation (19) can be changed to

The method of calculating the parameters of the rock damage constitutive equation is usually based on the rock triaxial stress-strain test curve using the graph method and the linear regression method. These methods have achieved good results, but the process is complicated, and it is difficult to achieve for some rocks, and the errors caused by insufficient data are large. In addition, the method of multivariate function extremum theory can be used to calculate parameters. According to the characteristic parameters of the test curve, such as peak intensity and peak strain, the parameters can be calculated. This method is not only simple, clear in physical meaning, and high in precision, but also suitable for various types of rocks.

The slope of the peak strength point of the rock stress-strain test curve is 0:

Therefore, the result of deriving on both sides of equation (23) iswhereand let

Deformation of equation (28) can give

The derivative result of equation (24) is

Then, substituting equation (30) into equation (26), it can be obtained that

Substituting the peak strength point into equation (23) gives the following:

Calculating equation (32), we get the following:

Substituting (27) into (28), we get

Deforming equation (34), the result is

Substituting (31 and 33) into (35), the following results are obtained:

So, according to the peak strength point of the stress-strain curve under a confining pressure and water pressure, a value can be obtained by equation (33). From the normal distribution function table, the corresponding value can be found.

The parameters and can be calculated using equations (29) and (37), respectively.

4. Discussion

4.1. Stress-Strain Curve

Based on the triaxial seepage tests of filled single-cracked sandstone (Table 3), the theoretical curve can be calculated by using the established constitutive model equation (23). The stress-strain model curve and the test curve of the cracked rock under seepage-stress conditions are shown in Figure 5.

It can be seen that this model has the following advantages: (1) the model can reflect the whole process of rock deformation and rupture. The model curve shows a shape consistent with the test curve, and the elastic stage of the curve agrees well, which fully reflects the evolution process of rock seepage damage. (2) It better reflects the variation of rock strength with confining pressure and crack length. The strength of rock increases with the increase in confining pressure, and it decreases with the increase in crack length. (3) From the postpeak softening stage, the damage model can also better reflect the real softening characteristics of the fractured rock after failure. (4) The model has few parameters, and the parameters are easy to calculate.

However, the model also has limitations: (1) the equivalent elastic modulus parameters used in the model cannot perfectly represent the initial compression stage of the cracked rock. Although the crack part of the rock is filled, the cement mortar still contains more pores, which makes the sample produce a larger strain during the compaction stage. Therefore, during the initial compaction stage of the filled fractured rock, the fit between the model curve and the test curve is poor. (2) The peak strength of the model curve is slightly lower than that of the test curve. In this test, the transient pulse technique is used to simultaneously test the permeability of the sample, causing the rock to creep. The viscoelastic parameters of the rock are not considered in the model, resulting in deviations in the calculated values. However, this error can be eliminated after changing the test method.

4.2. Model Verification

In order to verify the rationality of the rock damage statistical model established in this paper, the change law of the rock damage variable with strain is studied in combination with the test curve to analyze the correctness of the model damage variable (damage factor). Therefore, transforming equation (18), we can get

And, substituting equations (20)–(22) into equation (36), it is obtained as

According to and values of the filled cracked sandstone obtained in the test, the corresponding damage variable can be calculated by using equation (29). Thereby, the curve corresponding to the test curve is obtained, as shown in Figure 6.

According to the rock seepage softening damage model proposed in this paper and its parameter acquisition method, the damage factor in the model is calculated using equation (17).

The corresponding theoretical curve of can be obtained, as shown in Figure 7.

The damage variable D is a physical quantity used to measure the degree of rock damage. The damage variable evolution curve obtained from the test curve is usually s-type and transforms into a monotonous increase trend. The microcracks inside the rock rapidly increase in the nonlinear elastic stage within the strain range of the initial compaction stage. With the hardening of the confining pressure, the crack propagation gradually slows down to closure. After entering the linear elastic stage, the crack continuously extends to form a macroscopic crack, and the rock is cumulatively damaged and failure is at the peak strain.

However, since the elastic parameters used in the calculation model do not well characterize the deformation characteristics of the compaction stage, the damage value of the rock in the compaction stage is not reflected in the theoretical curve, and the damage value in the initial stage is zero.

Based on the results of (17) and (29), the test curves and theoretical curves of are compared:(1)With the increase in the axial strain, when the stress exceeds a certain level, the rock cannot bear the load and begins to produce damage. It can be seen that the peak value of the test curve and the peak value of the theoretical curve correspond to the same strain value, and the final damage variable is the same.(2)As the confining pressure increases, the strain value increases when the rock is damaged; as the crack length increases, the strain value decreases at the same time. This rule is consistent in the test curve and the theoretical curve.(3)In the test curve, when the axial strain is small, the change of the damage variable is not regular, and D is less than 0, which is contrary to its reasonable value interval (0 < D ≤ 1). In the theoretical curve, when the filled fractured sandstone is subjected to a low stress level and is in a nonlinear elastic stage, the damage variable should be constant at zero. This reflects the damage threshold effect of rock material.

Therefore, the model established in this paper can reflect the characteristics of rock damage related to the stress state and can also reflect the damage state affected by the rock structural plane. The damage factor D is nearly 0 in the compaction stage of the theoretical curve, which is consistent with the actual situation and more reasonable, reflecting the influence of the rock damage threshold.

Rock is a nonuniform quasibrittle material with complex mechanical properties. Considering the randomness of rock defect distribution and the nonuniformity of mechanical properties, the common research idea is to analyze the damage process by introducing statistical methods. The established statistical damage model introduces the failure criterion as a physical mechanic parameter into the nonuniformity of the material through a random distribution function. Lognormal distribution is widely used as the structural resistance probability distribution in structural reliability theory, and it can be considered that the random distribution of rock mechanical property parameters is similar to the resistance distribution of structural members. In describing the rock damage deformation curve, the lognormal model curve is more in line with the actual situation, and the curve shape and the final value are more accurate in this paper.

However, in geotechnical engineering, the types of filling in cracked rocks are very rich, such as hydraulic materials that can still play the role of bonding and strength in water like cement mortar, and there are also air-hardening materials that lose their cohesiveness and strength when encountering water. There are self-weight fillings as in this paper’s experiment, and there are also pressure cement fillings. At the same time, the fracture types of rock mass are more diverse, such as single fracture, multifractures, and en echelon joints. This study only reflects the applicability of the established lognormal statistical damage model in describing the seepage mechanical properties of single-fractured sandstone filled with cement mortar, which has little limitations in the experimental design and result data. In the future work, we will consider conducting more experimental studies on the seepage mechanical properties of fractured rocks with different fillings and fracture types under different stress environments to learn more about the applicability of the lognormal distribution statistical model.

5. Conclusion

In this paper, the permeability test of the fractured sandstone filled with cement mortar at a confining pressure of 20 MPa and 10 MPa and the water pressure at 7 MPa is completed. According to the tests data, the lognormal distribution is used as a probability distribution model of rock microunit strength to establish a rock stress-strain relationship constitutive model. According to the failure characteristics of the fractured rock in the seepage-stress field, the strength of the rock microunit is determined based on the maximum tensile strain failure criterion. And, considering the influence of confining pressure and water pressure, Hooke's law and the effective stress principle are used to improve the model parameter. The equivalent elastic modulus of the filled fractured rock is calculated by Voigt’s space average model. Finally, a statistical constitutive model of the lognormal distribution of seepage damage in the fractured rock is established. The following conclusions can be drawn:(1)The established rock damage statistical model agrees well with the experimental curve and can better reflect the stress-strain relationship of the filling fracture rock under different external load conditions and different structural plane types.(2)This model can fully reflect the characteristics of rock damage related to the stress state and also reflect the damage state affected by the characteristics of the rock structural plane.(3)The damage peak of the model curve is consistent with the actual situation, and its initial damage strain reflects the rock damage threshold effect, so the model is reasonable and meaningful.(4)This model has few parameters, simple acquisition, and clear physical meaning, which is extremely convenient in engineering application.

Data Availability

Demand data used to support the findings of this study are available from the first author upon request.

Conflicts of Interest

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


This study was supported by the National Natural Science Funds (41831278 and 51878249) of China, the Guangdong Province Water Resource Science and Technology Innovation Program (2017-30) of China, and the Postgraduate Research and Practice Innovation Program of Jiangsu Province (2018B661X14).