Abstract

The study on the constitutive fatigue model for rocks under cyclic loading has an important significance in rock engineering. In order to study the fatigue properties of rocks under cyclic loading, according to the theory of rheological mechanics and the existing three basic one-dimensional fatigue elements, i.e., elastic, viscous, and plastic fatigue elements, the three-dimensional elastic, viscous, and plastic fatigue elements were constructed. Meanwhile, a fatigue yield criterion for rocks under cyclic loading was proposed, and a three-dimensional nonlinear visco-elastoplastic fatigue constitutive model (NVPFM) for rocks was established by using the flow criterion related to the proposed fatigue yield criterion. Compared with the test results of rocks under cyclic loading, the three-dimensional NVPFM could be used to describe the transient, steady-state, and tertiary creep phases of rocks under cyclic loading.

1. Introduction

When rocks are subjected to cyclic or dynamic loading, different types of rocks have different responses. Some rocks may become stronger under cyclic loading, while others become more fragile [1]. In geotechnical engineering, fatigue problems are often involved such as coal walls under mining stress [25], side slope of the dam under cyclic variation water level [6], tunnels under traffic loading, and the surrounding rocks under earthquake load [712]. The strength of rocks under long-term cyclic loading will be reduced; thus, fatigue damage of rocks is caused, and the engineering accidents may be caused. Therefore, the study of fatigue characteristics of rocks under cyclic loading is of theoretical significance for the long-term stability evaluation and design of geotechnical engineering.

For rocks under long-term cyclic loading, even if the upper limit of cyclic loading is lower than the static strength of the material, the material can also be destroyed; this phenomenon is called fatigue [13]. There is no exception for rocks [1417]. The mechanical properties of rocks under cyclic loading are more complex than those under static load. The mechanical properties of rocks under cyclic loading are affected by many factors. The fatigue life of rocks is related to the stress amplitude of cyclic loading [1821], and the dynamic strength of rocks is related to the stress amplitude of cyclic loading and the stress upper limit of cyclic loading. Similarly, the dynamic strength of rocks is also affected by loading rate and the frequency of cyclic loading [2229]. There is a threshold value for rocks under cyclic loading, also known as critical strength. The deformation of rocks under cyclic loading is similar to the creep of rocks under static load [30]. When the upper limit of cyclic loading is higher than the critical strength, fatigue failure will occur [3133]. However, when the upper limit of cyclic loading is less than the critical strength, fatigue failure does not occur. In order to study the fatigue deformation characteristics of rocks under cyclic loading, some scholars studied the damage process of rocks from the perspective of damage mechanics. The damage of rocks under cyclic loading is divided into three stages [34]. Xiao et al. proposed an inverted-S damage model to describe the damage process of rocks under cyclic loading [35, 36]. Meanwhile, there are many constitutive models based on damage theory [3739]. Additionally, internal variable fatigue constitutive model of rocks was also proposed [40]. However, in recent years, the component model based on rheological mechanics has been applied to build the constitutive model of rocks under cyclic loading as well. Wang et al. [41] established a burgers fatigue model under cyclic confining stress loading with low cyclic frequency. Guo and Huang [42] proposed three one-dimensional fatigue elements, i.e., elastic fatigue element, viscous fatigue element, and plastic fatigue element, and used creep theory to establish a nonlinear fatigue constitutive model of rocks under uniaxial cyclic loading by simplifying the cyclic loading reasonably.

The research on the fatigue constitutive model for rocks is mainly based on the damage theory. Although the established constitutive models can effectively describe the rock fatigue deformation characteristic, most of the proposed models can only reflect the fatigue properties of rocks under uniaxial conditions. Moreover, some proposed models cannot reflect the influence of the properties of cyclic loading and strength of rocks on fatigue deformation of rocks. However, at present, the study of the three-dimensional fatigue constitutive model for rocks under cyclic loading is still lacking. Consequently, in this article, on the basis of rheological mechanics and the previous study [42], three three-dimensional fatigue elements, i.e., elastic, viscous, and plastic elements, are constructed. Moreover, the fatigue yield criterion for rocks under cyclic loading was proposed, and a three-dimensional nonlinear viscoelastic-plastic fatigue model was established using the flow criterion related to the fatigue yield criterion. The rock fatigue constitutive equations under triaxial conditions were deduced. Finally, the model parameters were identified, and model applicability was verified.

2. Basic Theory of Rock Fatigue

2.1. Cyclic Loading

The common types of cyclic loading are sinusoidal wave, cosine wave, rectangular wave, and triangular wave cyclic loading. The cyclic loading can be expressed aswhere is the upper limit of cyclic loading; is the lower limit of cyclic loading; is the frequency of cyclic loading, ; is the period of cyclic loading; is the time of cyclic loading, and the relationship between and is ; is a positive integer; and is the stress amplitude of cyclic loading.

2.2. Fatigue Behavior of Rocks

When the upper limit of cyclic loading is higher than the critical strength, fatigue failure will occur, and the creep behavior presents transient, steady-state, and tertiary creep phases (see II curve in Figures 1(a) and 1(b)) [3033, 44]. The creep phases of under cyclic loading are also called as deceleration, constant velocity, and acceleration creep. However, when the upper limit of cyclic loading is less than the critical strength, fatigue failure does not occur, and the creep behavior presents two creep stages, i.e., transient and steady-state (see I curve in Figure 1(a)) [43, 45].

2.3. Simplified Cyclic Loading

To use the component model based on rheological mechanics, assume that is the direction of the cyclic loading and it is the first principal stress. Under confining pressure, the equivalent deviatoric stress acting on rocks under cyclic loading was simplified into the following stress function in our previous study [46]:where is the average stress of cyclic loading, ; is the upper limit of cyclic loading; is the critical strength of the rock under the complex stress state, which is the maximum stress value from the curve of the rock under cyclic loading when the rock happens to be destroyed; is the peak stress from the curve of the rock under static load; and is the small principal stress on rock specimens.

According to equation (2), when the peak deviatoric stress of cyclic loading , the equivalent deviatoric stress , and the fatigue breakdown of the rock will occur. When the peak deviatoric stress of cyclic loading , the equivalent deviatoric stress , and the rocks do not undergo fatigue breakdown.

When confining pressure , the cyclic loading is simplified into the stress function shown in the following equation [42]:where is the average stress of cyclic loading; is the upper limit of cyclic loading; is the critical strength of the rock under uniaxial cyclic loading; and is the rock uniaxial compressive strength.

Similarly, according to equation (3), when , the simplified stress , and the rock will undergo fatigue breakdown. However, when , , and the rock does not undergo fatigue breakdown.

2.4. Fatigue Yield Criterion for Rocks
2.4.1. Rock Fatigue Yield Criterion under Uniaxial Stress

Under uniaxial cyclic loading, when , namely, , the rock will undergo fatigue damage. When , , and the rock will not undergo fatigue failure.

2.4.2. Rock Fatigue Yield Criterion under Confining Pressure

Regardless of the effect of spherical stress tensor on fatigue deformation, it is assumed that fatigue deformation of the rock is related only to the second invariant of stress deviator J2. According to equation (2), the proposed fatigue yield criterion of the rock under cyclic loading is as follows:where is the second invariant of the equivalent stress deviator, .

When , . When , , and the rock will undergo fatigue breakdown. Conversely, when , , and the rock does not undergo fatigue breakdown. Consequently, the established equation (4) can be used as a fatigue yield function for rocks under cyclic loading, and the yield criterion is a function of of cyclic loading and equivalent stress.

2.5. Basic One-Dimensional Fatigue Elements [42]

The one-dimensional fatigue model consists of three basic fatigue elements [42], i.e., elastic, viscous, and plastic fatigue elements. For the elastic fatigue element (see Figure 2(a)), the stress acting on the material is proportional to the strain. For the viscous fatigue element (see Figure 2(b)), the stress acting on the material is proportional to the strain rate. The plastic fatigue element is represented by friction sheets shown in Figure 2(c). The constitutive equations are as follows:

2.5.1. Elastic Fatigue Element [42]

where is the elastic fatigue coefficient, which can be determined by the numerical fitting method. is the fatigue strain.

2.5.2. Viscous Fatigue Element [42]

where is the viscous fatigue coefficient, which can be determined by the numerical method, is the strain increment of the viscous fatigue element within the time period , and is the fatigue strain increment of the element in each cycle of cyclic loading.

2.5.3. Plastic Fatigue Element [42]

When , that is, , even if the cycle number tends to infinity, the fatigue strain of the element is still 0. When , namely, , rocks will undergo fatigue breakdown, and the fatigue strain of rocks will continue to increase with the increase in the number of cycles of cyclic loading . The equation of the state of the plastic fatigue element iswhere is the plastic strain of the plastic fatigue element.

Based on the above establishment process of one-dimensional fatigue elements and rheological mechanics theory, the three-dimensional fatigue elements under confining pressure can be established in this article.

2.6. Establishment of Three-Dimensional Fatigue Elements
2.6.1. Three-Dimensional (3D) Elastic Fatigue Element

The three-dimensional elastic fatigue element is shown in Figure 3, and its equation of state iswhere is the equivalent deviatoric stress tensor , in which is the Kronecker delta; is the deviatoric strain tensor of the fatigue element; is the first invariant of the stress tensor; is the first invariant of the strain tensor; and is the fatigue volume modulus of the rock, and its physical dimension is stress.

Then, the total strain of the 3D elastic fatigue element under iswhere is the spherical stress tensor, .

2.6.2. Three-Dimensional (3D) Viscous Fatigue Element

For the 3D viscous fatigue element (see Figure 4), without considering the volume strain of the element under the spherical stress tensor, the state equation of the element iswhere is the increment of the deviatoric strain of the fatigue element in period and is the deviatoric strain increment of the viscous fatigue element in each cycle.

The following equation can be obtained by integrating equation (10):where is the deviatoric strain of the viscous fatigue element.

2.6.3. Three-Dimensional (3D) Plastic Fatigue Element

Three-dimensional plastic fatigue element can be represented by the friction sheet shown in Figure 5. When , namely, , and even if the number of cycles tends to infinity, the element fatigue strain is still 0. When , namely, , fatigue failure can occur in rocks, and the fatigue strain of rocks will still increase with the increase in the number of cycles of cyclic loading. The equation of state of the element iswhere is the plastic strain of the plastic fatigue element.

3. Three-Dimensional Nonlinear Visco-Elastoplastic Fatigue Constitutive Model for Rocks

The Burgers fatigue model consists of a series of Maxwell and Kelvin fatigue models, which can describe the transient and steady-state creep phases of the rock but cannot describe the tertiary creep phases of the rock. To fully describe the law of creep of the rock under cyclic loading, a nonlinear visco-elastoplastic fatigue model is connected to the burgers model in series, which is a one-dimensional nonlinear visco-elastoplastic fatigue constitutive model (see Figure 6) proposed by Guo and Huang [42].

Schiessel et al. [47] considered the time accumulation of the viscous elements (ηK) in the Kelvin body and defined the fractional viscous element by fractional derivatives. The creep compliance J(t) of the fractional element is proportional to tβ; β is a constant, and t is the rheological time. According to the aforementioned and Guo and Huang’s study [42], the viscous coefficient of the series-connected viscous fatigue element in the nonlinear visco-elastoplastic fatigue model, shown in Figure 6, iswhere is the parameter that reflects the strain rate of the tertiary creep phases of the rock under cyclic loading, . is the initial value of the viscous coefficient, and ; its physical dimension is stress·time n. Equation (13) can also be further expressed as

3.1. Establishment of the Nonlinear Viscoplastic Fatigue Model

Under confining pressure stress, the fatigue mechanics parameters of the rock are different from those under one-dimensional stress. The three-dimensional nonlinear viscoelastic-plastic fatigue model consists of a three-dimensional elastic fatigue element, a viscous fatigue element, a viscoelastic fatigue model, and a nonlinear viscoplastic fatigue model (see Figure 7). In the model, G1 is the elastic shear fatigue coefficient, G2 is the viscoelastic shear fatigue coefficient, η1 is the viscous shear fatigue coefficient, η2 is the viscoplastic shear fatigue coefficient, and η30 is the initial value of the viscoplastic shear fatigue coefficient. The physical dimension of model parameters G1 and G2 is stress.

3.2. Viscoelastic Fatigue Model

For the viscoelastic fatigue model shown in Figure 8, the effect of the spherical stress tensor on volumetric strain is not considered. The state equation of the viscoelastic fatigue model, shown in Figure 8, iswhere is the deviator strain increment of the fatigue element in the time period , is the deviator strain increment of the fatigue element in each cycle, and is the deviator strain tensor of the fatigue element.

The following fatigue constitutive equation of the fatigue element can be obtained from equation (15):

3.3. Nonlinear Viscoplastic Fatigue Model

For the viscoplastic fatigue model shown in Figure 9, the viscoplastic shear fatigue coefficient in the model is the same as equation (14), and the plastic strain rate of the fatigue model according to the rheologic basic principle can be written aswhere is the deviator strain increment of the fatigue element in the time period ; is the deviator strain increment of the fatigue element in each cycle; is the fatigue yield function of the rock; and is the plastic potential function of rock fatigue.

The following fatigue constitutive equation can be obtained by substituting equation (14) into equation (17) and integrating :where is the plastic deviator strain tensor of the nonlinear viscoplastic fatigue model. Among them,

In equation (20) [48], Bm, m, and M are constants. In this article, we used and for convenience.

When , the rock does not undergo fatigue breakdown. When , the following fatigue constitutive equation of the fatigue element can be obtained by using the flow rules related to the fatigue yield criterion:

3.4. Constitutive Equation of the Nonlinear Visco-Elastoplastic Fatigue Model

When the rock fatigue yield function , the rock will not be subjected to fatigue failure. At this time, the model degenerates into a burgers model which can describe the two creep stages of the rock under cyclic loading, i.e., transient and steady-state phases. From equations (8)–(20), the constitutive equation of the fatigue model can be obtained aswhere is the total strain tensor of the fatigue model.

When the rock fatigue yield function , the rock will undergo fatigue failure. At this time, the model is a 3D nonlinear visco-elastoplastic fatigue model which can describe the three creep stages of the rock under cyclic loading, i.e., transient, steady-state, and tertiary creep phases. From equations (8)–(21), the following constitutive equation of the fatigue model can be obtained:

4. Fatigue Constitutive Equation for Rocks under Triaxial Compression

The dynamic triaxial test is an indoor method for studying the fatigue properties of rocks under confining pressure stress. The adaptability of the proposed model can be verified by dynamic triaxial test data. Consequently, on the basis of three-dimensional rock fatigue constitutive equation (23), the fatigue constitutive equation under triaxial condition is deduced.

According to equation (2), the equivalent stress state of the rock under dynamic triaxial condition iswhere σ3 is the confining pressure acting on specimens.

Deviatoric stress under the simplified stress state is as follows:where δij is the Kronecker symbol, σm is the average stress, σm = (σ1 + 2σ3)/3, and σmδij is the spherical stress tensor.

The total strain of the rock under dynamic triaxial condition based on the 3D NVPFM is as follows:where , , , and are rock elastic, viscous, viscoelastic, and viscoplastic volume strains, respectively.

Without considering the viscous, viscoelastic, and viscoplastic volume strains of the rock and considering only the elastic volume strain, then  = 0. Consequently, the following equation can be obtained from equation (26):where eij is the total deviatoric strain tensor caused by deviatoric stress.

When , the rock will not undergo fatigue failure, and circumferential constitutive equation (28) and axial constitutive equation(29) of rocks under dynamic triaxial condition can be obtained by substituting equations (24)–(27) into equation (22):

When , the rock will undergo fatigue damage, and circumferential constitutive equation (30) and axial constitutive equation (31) under dynamic triaxial condition can be obtained by substituting equations (24)–(27) into equation (23):

5. Validation of NVPFM Adaptability

5.1. Parameter Determination and the Curve Fitting Method for Fatigue Constitutive Equations

For equations (28)–(31), the model parameters can be obtained by [48] regression analysis [49, 50]. 1stOpt is world’s leading software for nonlinear curve fitting and the comprehensive optimization analysis of nonlinear equation parameters. It has been widely recognized in the field of nonlinear regression, nonlinear curve fitting, and parameter estimation of the nonlinear equation. The genetic algorithm (GA) method is a built-in nonlinear optimization control algorithm in 1stOpt, which can solve the parameter value of the multiparameter nonlinear equation and fit the nonlinear curve according to the corresponding data [45]. In this article, the GA method in 1stOpt was used to solve the parameters of the proposed model and fit the curve of test data.

The proposed three-dimensional NVPFM can describe the creep of the rock under cyclic loading. The parameters of the three-dimensional NVPFM equation can be solved. According to the previous study [42], the detailed process of obtaining model parameters and fitting test data is as follows. According to the stress condition of the sample during the test, when rocks undergo fatigue breakdown, equation (30) is selected. However, when rocks do not undergo fatigue breakdown, equation (28) is selected. Then, the test data from the literature and the selected constitutive equation were input into software in the form of a program. Afterwards, the GA method was used to fit the fatigue curve of the rock under cyclic loading. After the program runs, the parameter values and the curve of the model can be obtained. In this study, only the constitutive equation of the rock under fatigue breakdown is studied.

5.2. Verification of the Fatigue Model under Cyclic Loading
5.2.1. Model Parameters and the Curve of Red Sandstone

Zhang et al. [51] used sinusoidal wave cyclic loading with 0.2 Hz frequency to conduct triaxial tests on red sandstone under confining pressure. The fatigue curve of red sandstone under 15 MPa confining pressure is shown in Figure 10. It can be obtained that the static strength of red sandstone under 15 MPa confining pressure was 104.653 MPa, and its critical strength was 98.22 MPa by Origin software. The deviatoric stress acting on the RS-4-3 specimen under cyclic loading was in the range from 51.8 to 99.1 MPa. Constitutive equation (30) based on the proposed fatigue model was used to fit the fatigue curve of the RS-4-3 specimen. The obtained parameters and the curve of the NVPFM are shown in Table 1 and Figure 11, respectively.

5.2.2. Model Parameters and the Curve of the Water-Rich Soft Rock

In order to study the dynamic characteristics and structural safety of the soft rock in the foundation of the high-speed railway tunnel, Ding [52] used the MTS rock dynamic system to carry out fatigue tests on the water-rich soft rock under different dynamic stress ratios ηd (see equation (32)) by using the loading method shown in Figure 12. The effects of dynamic stress amplitude and static deviatoric stress on the fatigue characteristics of the water-rich soft rock were studied by Ding [52]. The critical stress ratio of the water-rich soft rock under 200 kPa confining pressure, 180 kPa static deviatoric stress , 3 Hz loading frequency, and 50 kPa groundwater pressure is 0.3–0.5. In this paper, the critical dynamic stress ratio of the water-rich soft rock is 0.4, and the critical dynamic stress amplitude is 240 kPa. The cyclic loading is expressed by equation (33). Then, the upper limit (34) of cyclic loading can be obtained from equation (33). The critical strength of the water-rich soft rock was . For the dynamic strength rock, Ge [53] found that the critical strength of the rock under confining pressure is 0.85 to 0.9 times of the static strength of the rock. Therefore, the static strength of the water-rich soft rock was taken as in this study. The static strength of the water-rich soft rock under 200 kPa confining pressure was 353 kPa, and its critical strength was 300 kPa. Constitutive equation (28) based on the NVFPM proposed was used to fit the fatigue curve of the water-rich soft rock under cyclic loading with 300 kPa dynamic stress amplitude and dynamic stress ratio of 0.5. The fitting results and model parameters are shown in Table 2 and Figure 13.where σcua is the uniaxial compressive strength of the water-rich soft rock, and its value is 0.6 MPa.

5.2.3. Model Parameters and the Curve of the Water-Rich Soft Rock When Fatigue Breakdown Does Not Occur

In order to analyze the effect of dynamic stress amplitude on creep characteristics of the water-rich soft rock, Ding [52] conducted comparative tests on rocks under cyclic loading with dynamic stress amplitudes of 30, 60, 90, 120, and 180 kPa, respectively. The corresponding dynamic stress ratios under the five working conditions were 0.05, 0.1, 0.15, 0.2, 0.3, and 0.5. During testing, the loading frequency of cyclic loading was 3 Hz, the static deviatoric stress was 180 kPa, and the applied confining pressure was 200 kPa. From the data obtained from the test, the model parameters and model curves of the proposed fatigue constitutive equation (28) can be obtained. The obtained parameters and curves of the NVFPM are shown in Table 3 and Figure 14, respectively.

It can be seen from Figure 14 that fatigue constitutive equation (28) established in this paper can reflect the transient, steady-state, creep phases of rocks under cyclic loading. The correlation between model and test data was good. As can be seen from Table 3, when , the rock does not suffer from fatigue failure, and some of the obtained parameters of the model showed apparent regularity. The values of viscoelastic shear fatigue parameters and viscoelastic shear fatigue parameters decay exponentially with the increase in the dynamic strain ratio . The specific attenuation law is shown in Figures 15 and 16. For the specific change of other parameters that cannot be determined for the time being, further study needs to be made in combination with more tests.

From the above regression analysis results, it can be seen that the fatigue constitutive model proposed in this paper can describe the creep characteristics of the rock under cyclic loading. When , rocks do not undergo fatigue breakdown, and the proposed model would degenerate into the burgers model, which can describe the transient, steady-state, and creep phases of the rock. When , rocks undergo fatigue breakdown, and the proposed NVPFM can reflect the transient, steady-state, and tertiary creep phases of the rock under cyclic loading.

6. Conclusions

According to the rheological mechanics, three-dimensional elastic, viscous, and plastic fatigue elements were constructed. Moreover, a fatigue yield criterion for the rock under cyclic loading was proposed based on the simplified cyclic loading in the previous study. A three-dimensional nonlinear visco-elastoplastic fatigue model (NFVPM) was established under confining pressure by using the flow rules related to the proposed fatigue yield criterion. Finally, the rationality of the model is verified. The following conclusions can be drawn from this study.

The established NFVPM can describe the transient, steady-state, and tertiary creep phases of the rock under cyclic loading. Under confining pressure stress state, when the proposed fatigue yield criterion of the rock , the model degenerates into a burgers model that describes the transient and steady-state creep phases of the rock under cyclic loading. However, when and the parameter , the model is a NFVPM that describes the transient, steady-state, and tertiary creep phases of the rock under cyclic loading.

The parameters of the NFVPM can be easily solved using mathematical methods. The parameters and of the NFVPM decreased with the increase in dynamic stress ratio in the way of power function. Moreover, the curves of the NFVPM were in good agreement with the experimental results, and the correlation was above 0.990.

Data Availability

The data used to support the findings of this study are available within this article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the financial support from the National Natural Science Foundation of China (nos. 51968010 and 51608141), the Talents Scientific Research Foundation of Guizhou University (Grant no. (2015)16), and National Science and Technology Support Plan for the 12th Five-Year Plan of China (no. 2011BAJ03B00).