3.1 Time-Dependent Strain of Medium-grain Sandstone
3.1.1 Behavior of Axial Strain with Time
The creep deformation of saturated and dry medium-grain sandstone with time is shown in Fig. 5. The instantaneous strain of the specimen is obvious at each stage of deviatoric stress, and the instantaneous strain is the largest under the first stage of axial loading. The first-stage instantaneous deformation of saturated rock under different confining pressure differs slightly and remains at about 0.35%, while the first-stage instantaneous deformation of dry rock decreases with the increase of confining pressure, which shows that confining pressure restrict the instantaneous deformation of dry rock more obviously. After the first stage of loading, with the increase of loading level, the instantaneous deformation increases gradually, and the gap between saturated and dry rock axial deformation increases. This also confirms the view that confining pressure has great restrictive effect on dry rock.
When the deviatoric stress is not large enough, such as confining pressure 3.0 MPa and deviatoric stress less than 14 MPa, the rock samples only show instantaneous deformation and decay creep, and the rock deformation reaches a stable state in a short time. When the deviatoric stress increases to a certain extent, the deformation of rock shows obvious time correlation. At this time, it takes a long time for the rock deformation to reach a stable state. With the increase of loading level, the steady creep rate of rock increases gradually until accelerated creep occurs, which eventually leads to rock failure.
The increase of confining pressure improves the maximum deformation and bearing capacity of medium-grain sandstone. When saturated rock is failure, the loading level is low and the axial strain is large, which indicates that the deformation capacity of medium-grain sandstone under the influence of water is increases and the bearing capacity decreases.
3.1.2 Behavior of Lateral Strain with Time
The trend of lateral strain is similar to that of axial strain, as shown in Fig. 5. However, the dynamic principle of lateral deformation and axial deformation is different (Zhang et al. 2016). In the early stage of creep test, the lateral compactness of primary cracks and voids in specimens is generally less than the axial strain. With the increase of deviatoric stress, the original defect is compacted to the minimum, and the restructuring of particles in rock enhances the lateral strain, and the time dependence of lateral strain is greater than that of axial strain. Thereafter, the lateral deformation of rock increases significantly and eventually leads to lateral dilatancy failure in the loading stage near rock failure.
For medium-grain sandstone in coal measures, the instantaneous lateral strain under deviatoric stress is less than the instantaneous axial strain, and the lateral deformation of rock increases faster than the axial deformation after load stabilization. Under the same confining pressure, the difference of lateral strain between saturated and dry rock samples is very small, which indicates that the maximum lateral deformation ability of both samples is similar, and their maximum lateral strain decreases with the increasing of confining pressure.
3.1.3 Behavior of Volumetric Strain with Time
Volume deformation of rock generally includes two stages: volume compaction and creep dilatation (Zhang et al. 2016; Belmokhtar 2017). The volume compaction of medium-grain sandstone is usually reflected in the first or two loading stages, and the subsequent loading stages are characterized by dilatancy behavior. The strain rate in compaction stage is smaller than that in dilatancy stage, and the dilatancy strain increases sharply when rock is near failure.
The maximum volumetric compressive strain of dry rock increases with the confining pressure, and the maximum dilatancy strain decreases with the increasing of confining pressure[]. For saturated rock samples, the difference of volumetric compressive strain under different confining pressure is very small, which remains about 0.2%, while the volume dilatancy under low confining pressure is the most obvious. In addition, the maximum volumetric compression strain of saturated rock is larger than that of dry rock, and the maximum volumetric dilatation strain is smaller than that of dry rock, and these differences gradually decrease with the increase of confining pressure.
3.2 Instantaneous elastic modulus
In multi-axial creep loading, both instantaneous loading and constant axial stress can cause rock damage, and the damage between two adjacent loading levels can be indirectly compared by the elastic modulus between them (Zhao et al. 2017). Based on this, we calculate the elastic modulus of rock samples during each instantaneous loading, as shown in Fig. 6. The elastic modulus of saturated and dry medium-grain sandstone decreases exponentially with the increasing of deviatoric stress. That is to say, with the increase of deviatoric stress, the rock damage increases, which indirectly indicates that the damage of medium-grain sandstone increases exponentially with the increasing of deviatoric stress.
3.3 Strain rate analysis
To study the difference of axial creep rate between saturated and dry medium-grain sandstone under different confining pressures, the creep rate of rock samples under three kinds of deviatoric stress is shown in Fig. 7. At the same deviatoric stress level, the creep rate decreases with the increasing of confining pressure, and the influence becomes more significant with the increasing of deviatoric stress. The reason may be that the confining pressure increases the compressive strength of rock and weakens the creep mechanical properties. Under the same confining pressure, the creep rate of saturated rock sample is higher than that of dry rock sample, and with the increase of confining pressure, the creep rate of saturated rock sample decreases more greatly. In practical engineering, support and reinforcement can effectively improve confining pressure and long-term stability of the project.
Ma et al. (2006) pointed out that the steady state creep rate of tuff lava has a power function relationship with deviatoric stress, and the higher the stress, the more obvious the exponential relationship is. In this paper, the exponential function relationship (y = AeBx) between steady creep rate and deviatoric stress of medium grain sandstone in is shown in Fig. 8. A and B are both experimental parameters, which are greater than 0. The fitted exponential relation is consistent with the conclusion that the elastic modulus decreases exponentially with deviatoric stress. That is, the higher the deviatoric stress, the faster the damage of rock and the faster the creep rate.
3.4 Creep model analysis
The creep experimental data can be analyzed and displayed through constitutive models. The parameters in these models can correspond to the inherent creep characteristics of rock materials. so it is important to quantitatively obtain the creep characteristics of materials for engineering stability prediction and engineering design.
For medium-grained sandstone, the rock presents deceleration creep and steady creep under deviating stress, i.e. linear creep. This creep behavior can be well described by Burgers model. However, according to the above experimental results, the damage of medium-grain sandstone accumulates exponentially under multi-axial creep loading, and the isochronal stress-strain relationship changes nonlinearly. Therefore, we should use the damage Burgers model to describe the creep behavior of medium-grain sandstone, and its expression is as follows (Zhu et al. 2010):
where ε is the strain, σ is the stress, t is time, E0 M, η0 M, E0 K, η0 K are viscoelasticity coefficient of damage Burgers model, εc is viscous strain.
Table 2
Parameters of creep model describing axial deformation of medium-grained sandstone
σ3 (MPa)
|
E0 M (GPa)
|
P-E0 M (GPa)
|
η0 M (GPa h)
|
P-η0 M (GPa h)
|
E0 K (GPa)
|
P-E0 K (GPa)
|
η0 K (GPa h)
|
P-η0 K (GPa h)
|
Saturated
|
1.5
|
6.12
|
6.36
|
5.83
|
6.56
|
52.43
|
61.92
|
158.47
|
146.89
|
3.0
|
6.78
|
7.49
|
8.15
|
7.72
|
68.75
|
72.20
|
159.71
|
169.16
|
4.5
|
9.27
|
8.47
|
8.89
|
8.73
|
91.05
|
81.20
|
188.54
|
188.66
|
6.0
|
10.86
|
9.32
|
10.41
|
9.63
|
113.76
|
92.27
|
203.46
|
209.13
|
Dry
|
1.5
|
9.31
|
9.66
|
9.51
|
9.35
|
96.64
|
105.51
|
138.50
|
139.96
|
3.0
|
10.08
|
10.80
|
9.52
|
10.45
|
110.52
|
116.94
|
153.28
|
153.81
|
4.5
|
12.48
|
11.52
|
11.89
|
11.15
|
121.95
|
124.16
|
164.26
|
162.57
|
6.0
|
14.17
|
12.79
|
13.05
|
12.21
|
134.83
|
135.57
|
177.26
|
174.65
|
It can be seen from Formula (1) that there is a non-linear term (1/e− mεc) in the Burgers model with damage, which is consistent with the conclusion that the damage increases exponentially with the deviatoric stress. m is a coefficient related to rock properties. The greater the m value, the more serious the material damage is under the same viscous strain. Based on the research results of Zhu et al. (2010), the satisfactory fitting result can be obtained by fitting the curve with m = 40.
Based on the creep test results, the Levenberg-Marquardt algorithm is introduced into the least squares method for solving nonlinear problems to identify the parameters in the creep model. The identified parameters of Burgers model with damage are shown in Table 2. The model was found to be suitable to describe the experimental data with good accuracy (R2 = 0.91 − 0.99).
The results in Table 2 show that the saturated rocks with smaller stiffness also have smaller creep parameters than the dry rocks, which results in larger creep deformation. The creep parameters of medium-grain sandstone increase with the confining pressure, and the relationship between rock stiffness and confining pressure is also positive. Therefore, we can infer that the creep model parameters of medium-grained sandstone correspond to its elastic modulus. This inference will make it possible for us to obtain the parameters of creep model based on the instantaneous deformation parameters of rocks.
Maxwell parameter E0 M reflects the instantaneous deformation of rock specimen under axial loading. E0 M in Table 2 shows that the instantaneous deformation of rock before accelerated creep increases with the increase of confining pressure. As for the variation of E0 M with deviatoric stress, taking 4.5 MPa confining pressure as an example, the relationship of E0 M with deviatoric stress as shown in Fig. 9 shows that the instantaneous deformation of medium-grain sandstone increases with the increase of deviatoric stress, which reflects the non-linear characteristics of rock creep, and also confirms the phenomena illustrated by the isochronal stress-strain curve. In addition, the first creep (deformation recoverable) and the second creep (deformation non recoverable) properties of rocks can be reflected by Kelvin parameters η0 K and Maxwell parameter η0 M, respectively. It is found that the Maxwell viscosity η0 M of medium-grain sandstone is usually much smaller than the Kelvin viscosity η0 K, indicating that non recoverable deformation is more easier to produce on medium grain sandstone.
3.5 Relationship between instantaneous and creep parameters
From the above analysis, it can be seen that there is a corresponding relationship between short-term deformation characteristics and creep deformation characteristics (Hamzaa O and Stace R 2018). Therefore, the creep model parameters under different confining pressures are normalized, i.e. the creep model parameters are divided by the elastic modulus (E) under the same confining pressure. The creep parameters after normalization are shown in Fig. 10. The results show that the regularity of creep properties of saturated and dry medium-grain sandstones after normalization is more obvious, which further shows that the instantaneous properties of rocks are proportional to the time-dependent properties.
Therefore, parameters of Burgers model with damage and elastic modulus (E) are plotted in Fig. 11. It is found that the relationship between them can be described by a linear function, and the accuracy of the description is more than 0.89. That is to say, the creep parameters of saturated and dry medium-grain sandstones can be calculated according to the formula parameter = a + bE when the elastic modulus (E) of rocks is known. The values of a and b are given in Table 3.
Table 3
Values of constants a and b used in the equation correlating damage Burgers parameters
|
E0 M (GPa)
|
η0 M (GPa h)
|
E0 K (GPa)
|
η0 K (GPa h)
|
Saturated
|
a
|
-0.1134
|
-0.1346
|
2.4596
|
18.1160
|
b
|
1.7521
|
1.8101
|
16.0702
|
34.8047
|
R2
|
0.9475
|
0.9689
|
0.8979
|
0.9749
|
Dry
|
a
|
-2.3363
|
-2.2295
|
-14.7205
|
-5.9461
|
b
|
2.3252
|
2.2448
|
23.3022
|
28.2757
|
R2
|
0.9652
|
0.9657
|
0.9929
|
0.9991
|