Creep Characteristics of Different Saturated States of Red Sandstone after Freeze-Thaw Cycles

To investigate the creep mechanical characteristics of rocks in different saturated states after freeze-thaw cycles, samples with different saturations (30%, 50%, 70%, 90%, and 100%) were selected for nuclear magnetic resonance (NMR), scanning electron microscopy (SEM), and uniaxial compression creep tests. The internal microscopic damage of the rock sample and mechanical characteristics under long-term loading are analyzed after the action of freeze-thaw cycles. The test results show that, as the saturation increases, the T2 spectrum distribution shifts to the right. The spectrum area gradually increases as the porosity increases. The critical saturation of freeze-thaw damage occurs when the saturation increases from 70% to 90%. It can be seen from the SEM image that the number of pores inside the rock samples gradually increases with increases in saturation, leading to the appearance of cracks. Under long-term loading, the saturation has a significant influence on the time-efficiency characteristics of sandstone freeze-thaw. As the saturation increases, the creep deformation gradually increases. After reaching 70%, the axial creep strain increases significantly. The rate of creep is accelerated, the creep failure stress is reduced, and the creep time under the last level of stress is significantly increased. A modified viscous-plastic body is connected in series to the basic Burgers model, the freeze-thaw-damage viscous element is introduced, and the creep parameters are fitted using test data. The results will provide a reference for long-term antifreeze design for rock engineering in cold areas.


Introduction
Rock masses in cold regions are subjected to the thermal action of seasonal changes and day-night cycles for long periods of time, leading to severe damage and deterioration. Geotechnical engineering in mines and tunnels in the cold regions of western China operate in freeze-thaw environments, with significant loads and surrounding rock pressure for extended periods of time. Rock creep under long-term loading significantly affects the long-term stability of rock masses. If the long-term effects of freeze-thaw and load are not considered at the same time, it will bring safety hazards to engineering construction in cold areas. Therefore, studying the creep characteristics of rock-mass engineering under the action of freeze-thaw cycles and investigating the time-efficiency characteristics play an important role in ensuring the long-term stability of rock-mass engineering in cold regions.
In recent years, Li et al. [1] investigated the influence of the number of freeze-thaw cycles and graded loads on creep and established a nonlinear creep constitutive model of freeze-thaw sandstone. Zhou et al. [2] conducted a triaxial unloading creep test on sandstone; the results showed that both the unloading process and the freeze-thaw cycle can improve the creep deformation of the rock sample, where radial deformation is more sensitive. Chen et al. [3] conducted triaxial creep tests on quartz sandstone after different freeze-thaw cycles. The research showed that the microdamage caused by freezing and thawing can lead to strong creep characteristics; the creep mechanical parameters are obviously changed with the increase in freeze-thaw cycles. Yang et al. [4] carried out triaxial creep tests and microscopic tests under the freeze-thaw cycle and systematically analyzed the operating mechanism of the effects of freezing and thawing on the creep characteristics of gneiss. Liu et al. [5] carried out a triaxial creep test on carbon shale under freeze-thaw action and determined the nonlinear characteristics of creep damage. It was concluded that, under the combined action of freeze-thaw and creep, the freeze-thaw action would accumulate rock damage and accelerate creep failure. Zhang et al. [6] proposed the freeze-thaw rock damage viscous element and constructed a constitutive model for granite freeze-thaw shear creep. Zhou et al. [7,8] carried out mechanical and micro-microscopic tests on prefabricated fractured sandstone after freezing and thawing; the interaction mechanism of moisture inside the defective rock mass and the evolution of the fatigue damage zone of the rock mass have been fully studied under the action of freezing and thawing cycles.
It can be seen from the above research that the freezethaw cycle has a significant influence on the long-term mechanical properties of rocks; however, most current creep research is only carried out in the fully saturated state. Owing to the several changes in geology, geography, climate, and hydrology over the years, the water content of rock masses differs significantly.
McGreevy and Whalley [9] believed that the initial water content of a rock mass determines the degree of frost-heave damage, and the water content will fluctuate with the number of freeze-thaw cycles, freezing duration, and seasonal changes. Liu et al. [10] found that the freeze-thaw damage of rock masses is dependent on the initial saturation, which has an important influence on the frost-heave force. Only the rock masses that exceed the critical value of saturation will lead to effective freeze-thaw damage [11,12]. Chen et al. [13] carried out uniaxial compression, longitudinal wave velocity, and porosity-measurement tests on tuffs with different initial saturations after the freeze-thaw process, and the results showed that a saturation of 70% is the critical value for effective freeze-thaw damage. The level of saturation has a significant influence on the long-term mechanical properties of the rock mass. Al-Omari et al. [14] found that the critical saturation is the main factor that dictates the freeze-thaw failure of limestone in the Chambord Castle in France. With the goal of understanding the time-efficiency deterioration of rocks at different saturations, Zhang et al. [15] proposed an elastic-viscous-plastic model based on a stress function of the influence of freeze-thaw cycles on the characteristics of saturated rocks and the rheological theory. Yang et al. [16] conducted a triaxial creep test for soft rock under water-bearing conditions, where the change law of creep deformation was analyzed. Liu et al. [17] tested slabs under different water-bearing conditions using the rock creep test and established the FLAC3D creep constitutive secondary development model, which considers the effect of water-bearing degradation. Zhang et al. [18] took the auxiliary tunnel of Jinping II Hydropower Station as an example to study the time-dependent deformation of rock tunnels and established a novel four-element fractional viscoplastic (FVP) model based on the test results. Wang et al. [19,20] carried out uniaxial long-period cyclic load creep tests on rock salt specimens under different maximum cyclic stresses and cyclic periods and analyzed the law of creep stage, creep rate, and elastic modulus with stress and number of cycles.
The above studies have laid the foundation for understanding rock creep characteristics and the influence of saturation on rock instantaneous mechanical properties under the effect of freeze-thaw cycles. However, there are few reports on the long-term mechanical properties of rocks with different saturation levels under freeze-thaw cycles. This paper conducts NMR, scanning electron microscopy, and uniaxial compression creep tests on red sandstone in different states of saturation after freeze-thaw cycles and analyzes the microscopic and creep damage mechanisms. Based on the Burgers model, a freeze-thaw damage creep model-considering the effect of saturation-is established. Further, the experimental value and fitted theoretical values are compared, finally obtaining the creep parameters. The results of this research will provide a reference and basis for long-term antifreeze design for rock engineering in cold areas.

Materials and Methods
2.1. Rock Samples. The rock samples were taken from a tunnel slope in a cold region in western China. Owing to the differences in water content of slopes at different heights, the stability of the slope will gradually weaken under the action of freeze-thaw cycles and long-term loads [21]. The D8Venture X-ray single crystal diffractometer was used to analyze the mineral composition of the rock samples. The mineral composition is shown in Table 1. The internal particles of the rock sample were pore-type cemented, which is a moderately swollen weakly cemented soft rock.
In accordance with the test regulations of the International Society of Rock Mechanics, the rock sample was processed into a standard cylinder with a diameter of 50 mm and a height of 100 mm. A rock sample with good integrity was selected and baked in an oven at 105°C for 24 hours; the primary wave velocity and dry density of rock samples were measured after cooling. The variation range of primary wave velocity and dry density states that (1832 < v < 1865 m/s) and (2:30 < ρ < 2:32 g/cm 3 ), respectively. Three groups were subjected to nuclear magnetic resonance (NMR), conventional uniaxial compression tests, and uniaxial compression creep tests, totaling 15 rock samples, as shown in Figure 1(c). The average values of the initial basic physical parameters for the rock samples are shown in Table 2.

Preparation of Rock Samples with Different Saturated
States. All rock samples were subjected to forced saturation for 24 hours using a vacuum saturating machine to obtain completely saturated rock samples. The rock sample was subsequently naturally air-dried, and its quality was monitored at regular intervals until the corresponding mass of where m t is the corresponding weight of the rock sample in different saturated states, m d is the weight of the dry rock sample, ω represents the saturation of the rock samples, and m s is the weight of the saturated rock sample. After obtaining the required saturation, the rock sample is sealed using plastic wrap and then packed in a plastic bag and sealed in water to keep the rock-sample saturation constant. Figure 1(e) shows the water sealing device for the rock samples.

Freeze-Thaw Cycle
Test. Before the freeze-thaw cycle, the rock samples with different saturated states were placed in the water sealing device for a week to achieve a uniform distribution of water in the rock samples. The freeze-thaw cycle is carried out using a water-sealed device to keep the saturation stable. The rock sample is frozen for 12 hours at -20°C and thawed for 12 hours at 20°C. A freeze-thaw cycle takes place every 24 hours. The number of freeze-thaw cycles is set to 30, the approximate number of days in a month. The equipment for the freeze-thaw test is shown in Figure 1(d).

Nuclear Magnetic Resonance Test.
After the freezethaw cycles, the first set of rock samples were vacuum saturated for 24 hours, and then, the nuclear magnetic resonance test was conducted. The nuclear magnetic resonance test equipment is shown in Figure 1(g). We used an NMR instrument (MacroMR12-150H-I, Suzhou Niumag Analytical Instrument Corporation, Suzhou, China) with a magnetic field strength of 0:3 ± 0:05 T, H proton resonance frequency of 12.77 MHz, and coil radio frequency pulse frequency of 1.499 MHz.

Creep Test.
The conventional uniaxial compression test was performed on the second set of rock samples with different saturated states to obtain the stress level of creep. The loading is controlled by the deformation, with a loading rate of 0.002 mm/s, and the peak strength of rock samples with different saturation was obtained. The third group of rock samples determined the graded loading stress levels according to 30%, 40%, 50%, 60%, 70%, 80%, and 90% of the second group's peak strength. Each level was loaded for 24 h, and the loading rate was 0.02 MPa/s. The surrounding area of the rock samples was smeared with petroleum jelly and wrapped in a thin plastic rubber sleeve before loading to ensure the stability of the internal saturation during the loading process. The test loading device uses the TAW-1000 rock mechanics testing machine, which can perform conventional uniaxial and triaxial compression and creep tests. The testing machine consisted of four main units: conventional triaxial, servo loading, deformation monitoring, and data acquisition units. The maximum axial stress can reach 1000 KN. The maximum confining stress can reach 100 MPa. The loading device is shown in Figure 1(f).

Nuclear Magnetic Resonance Test after Freeze-Thaw Cycle
3.1.1. T 2 Spectrum Distribution Analysis. The rock sample mainly comprises a mineral skeleton and pores. In a fully saturated state, NMR can be used to detect the signal strength of water in each pore and the total porosity, which in turn can determine the nuclear magnetic resonance characteristics. The particle and surface characteristics of the rock samples are reflected by the specific surface area. The lateral relaxation time T 2 of the rock samples is determined by the fluid volume and the specific surface area in the rock pores [12]. Therefore, the T 2 spectrum can clearly and intuitively quantify the internal pore structure of the rock sample. Figure 2 shows the distribution of the T 2 spectrum for rock samples with different saturated states after freezethaw cycles. The figure declares that, after the freeze-thaw cycles, the pore distribution of the rock samples presents a three-peak distribution and two peaks are the main representative forms. The variation trends of the T 2 spectrum distribution curves are very similar. They all gradually move to the upper-right corner as the saturation increases. Based on the principle that the relaxation time T 2 is positively correlated with the pore size, combined with the position where the T 2 spectrum distribution peak appears, different relaxation time intervals can be attributed to different pore types, namely, (r < 0:1 μm, T 2 < 7 ms) micropores, (0:1 < r < 1 μm, 7 < T 2 < 700 ms) mesopores, and (r > 1 μm, T 2 > 700 ms) macropores [22]. From the lateral relaxation time T 2 , we can see that the pore-distribution characteristics change significantly as the saturation reaches 90%. After freezing and thawing, mesopores replace the small pores to a significant extent and become the main pore form, indicating that after the saturation reaches 70%, the rock sample freeze-thaw damage is mainly reflected in the development of micropores, which causes an increase in mesopores.

T 2 Spectrum
Area and Porosity Analysis. The T 2 spectrum area of NMR can reflect the changes in the size and number of pores in the rock, and the peak area is positively correlated with the corresponding pore size and number [23]. The T 2 spectrum area is used to quantify the law governing the changes in the internal structure of rock samples with different saturations under freezing and thawing conditions, facilitating the evaluation of the degree of rock damage. Table 3 shows the total area of the T 2 spectrum and the areas of micropores, mesopores, and macropores under different saturation conditions after the freeze-thaw cycle. It is clearly observed that the total area of the T 2 spectrum gradually increases with increases in saturation.
During the freeze-thaw cycle of the rock samples, there is mutual conversion among the three types of pores. The degree of increase in the volume of each type of pore is not the same. The percentage of differing pore volumes in the total pore volume of the rock samples changes at different saturation levels. The proportion of the total area of the T 2 spectrum is shown in Figure 3. It can be seen from the figure that the inside of the rock sample is mainly composed of micropores and mesopores after freeze-thaw cycles, which account for about 50% of the total pore volume-the ratio of macropores to the total volume of pores is below 2%. The damage of the rock samples caused by the freeze-thaw cycles is mainly in the form of micropores and mesopores. With increases in saturation, the percentage of mesopore area gradually increases, and when the saturation reaches 90%, the peak area of mesopores gradually becomes the main peak. It shows that, when the saturation increases from 70% to 90%, the internal frost-heave force of the rock sample increases continuously, and the critical value of freeze-thaw damage is reached [11].

Porosity Changes with Saturation.
Porosity is the most commonly used and most important index for evaluating the characteristics of pores in porous materials. In this study, the change in the porosity of the rock sample after the freeze-thaw cycles was calculated using Equation (2) via the saturation method [24]: where Φ is the porosity of the rock sample, m s is the weight of the saturated sample, m d is the weight of the dry rock sample, ρ w is the density of saturating water, and V b is the bulk volume of the rock sample. Figure 4 shows the change curve for porosity and the T 2 spectrum area of rock samples with respect to saturation after freeze-thaw cycles. The porosity of the rock sample exhibits a nonlinear growth trend with increases in saturation, indicating that the internal volume of the rock sample increases with an increase in the number of freeze-thaw cycles. When the saturation increases from 30% to 70%, the porosity only increases by 13.66%, the total area of the T 2 spectrum increases by 17.92%, and the rock sample damage is low. When the saturation changes from 70% to 100%, the porosity increased by 23.60%, the total area of the T 2 spectrum increased by 28.79%, and critical freeze-thaw damage occurred. Therefore, it is obvious that the pores grow slowly in the first stage and rapidly in the second. The damage inside the rock sample is smallest in the first stage; therefore, the porosity and total area of the T 2 spectrum increase slowly. The increase in porosity is mainly due to the soluble minerals and cement in the rock sample gradually dissolving with the increase in saturation, weakening the connections between particles. This makes the pore volume increase continuously, and the damage caused by frost heave is small. However, the effective freeze-thaw damage occurs in the second stage. As the saturation increases, the frost-heave force increases and the internal fissures of the rock sample expand. The continuous dissolution of soluble minerals and cement by water molecules aggravates the destruction of the rock structure. The pore size inside the rock sample changes significantly, resulting in an increase in the porosity and the total area of the T 2 spectrum [10]. From the above analysis, it can be concluded that the change in the internal microstructure of the rock sample by freeze-thaw cycles is greatly affected by the saturation level.

Scanning Electron Microscope Analysis.
To further analyze the microstructural characteristics of the rock sample affected by saturation, slices were taken from the same position inside the rock sample after the freeze-thaw cycle, and the scanning electron microscopy test under 100x magnification was performed. The scanning result is shown in Figure 5. The PCAS software binarizes the information of the scanning electron microscope [25], and the black portion of the binarized image represents the pores and cracks of the   4 Geofluids scanned slice. From the binarization image, it is clearly observed that as the saturation gradually increases, the number of pores inside the rock sample increases and significant cracks develop. Under low-saturation conditions, owing to the small frost-heave force on the rock sample, even after 30 freeze-thaw cycles, the mineral particles of the rock sample are still tightly bonded, the grain boundary is not obvious, and the pores are fewer. As the saturation increases, particularly when it reaches 90%, the frost-heave force gradually increases, the bonding state of the cement and mineral particles in the rock sample gradually change, and the dissolution of the cement gradually accelerates. As a result, there is a significant loss of cement between particles, weakening their connection. Further, the dissolution holes increase in     5 Geofluids number and penetration cracks through particles appear. The above analysis demonstrates that, when the saturation reaches 70%, freeze-thaw cycle damage gradually appears, a value that can be used as the critical saturation for freezethaw damage.

Uniaxial Strength Test Analysis.
To determine the loading classification, uniaxial compression tests on rock samples with different saturation states were carried out. The stressstrain curve is shown in Figure 6, and the test results are shown in Table 4. In the curve, OA is the compaction stage, AB is the linear-elastic stage of pore and crack development, and BC is the plastic development stage [26]. The table shows that the strength and elastic modulus both decrease with increases in saturation. When the saturation increases from 30% to 70%, the decrease in strength and elastic modulus is smaller, decreasing by 38.86% and 28.64%, respectively; when the saturation is increased from 70% to 100%, it decreases by 51.29% and 62.5%, respectively, which represents a much larger decline than the former case. This shows that, under the action of freeze-thaw cycles, the instantaneous mechanical properties of rock samples are greatly affected by saturation, and the critical value of instantaneous mechanical damage occurs when the saturation increases from 70% to 90%.

Influence of Saturation on Creep Deformation in Rock
Samples. Figure 7 shows the uniaxial creep curves of rock     samples with different saturated states after freeze-thaw cycles. The figure is visible that the rock sample immediately produces instantaneous strain under the action of axial stress, which gradually increases with increases in saturation. Under the action of long-term stress, the creep characteristics of rock samples are relatively obvious. As the stress increases, the creep deformation increases significantly, and the characteristics of decelerated creep and stable creep appear [27]. Under the action of the final load, the characteristics of accelerated creep appear, and the creep rate and creep deformation gradually increase with increases in creep time, until failure. As the saturation increases, the creep deformation gradually increases under the load of the same level, and the creep curve gradually becomes steeper. Looking at Table 5, under the first level of horizontal load, relative to the previous saturation level, the axial creep strain increment of the rock samples with 50%, 70%, 90%, and 100% saturation are 0.0102, 0.0508, 0.2445, and 0.0237, respectively; under the second level of horizontal load, the additional values are 0.105, 0.0369, 0.3332, and 0.029, respectively. From the above analysis, it can be concluded that, under the load of the same level, the axial creep strain of the rock sample gradually increases with the increase of saturation. This is because the internal damage of the rock sample is greatly affected by saturation under the action of a freeze-thaw cycle. According to the analysis in sections, when the saturation reaches 70%, the internal particle cementation of the rock sample weakens, and the number of pores increases, leading to a reduction in creep mechanical properties.

Influence of Saturation on the Creep Rate of Rock
Samples. Figure 8 shows the steady-state creep rates of rock samples with different saturated states at 40%, 60%, and 80% loading stress levels. This is evident from the fitting curve that, under the same stress level, the steady-state creep rate of the rock sample increases nonlinearly with increases in saturation. Taking the 60% stress loading level as an example, when the saturation is low, the steady-state creep rate of the rock sample increases only slightly. When the saturation is 30%, the steady-state creep rate of the rock sample is 4:204 × 10 −7 , and when the saturation is 50%, it is 5:562 × 10 −7 , exhibiting an increase of 32.31%. As the saturation continues to increase, the increase in the steady-state creep rate of the rock sample gradually increases. When the saturation is 70%, the steady-state creep rate is 9:571 × 10 −7 ; when it reaches 90%, the steady-state creep rate is 1:634 × 10 −6 , which is an increase of 70.73%. The same change pattern is also shown at the 40% and 80% loading stress level. It can be concluded that saturation has a great influence on the steady-state creep rate of the rock samples. As the saturation increases, the steady-state creep rate of the rock sample increases gradually. This could be because of the difference in the degree of damage caused by the freeze-thaw cycle to the rock samples in different saturated states. When the saturation reaches 70%, owing to the increase in internal water molecules, the frost-heaving force increases, causing internal pores and microcracks in the rock-sample structure to develop and expand, resulting in increased damage and the gradual deterioration of mechanical properties.
Under the final stage of loading, the rock samples with different saturated states underwent creep failure. The change of the curves of axial strain and creep strain rate under the final stress level with time is shown in Figure 9. To visually analyze the curve characteristics of the accelerated creep stage of a 100% saturation rock sample, only the portion of the curve after the axial-strain change is shown in the figure. Under the action of freezing and thawing, the rock samples with different saturated states experienced the decelerated-creep stage, the stable-creep stage, and the accelerated-creep stage, and rapid destruction occurs in the accelerated creep stage; the creep-failure law is different under the action of the final level of stress. Further, 90% and 100% rock samples have a shorter accelerated creep stage, because when they reach a higher saturation, freezing and thawing cause their plasticity to decrease; as shown in the figure, the creep rate curve is approximately barrelshaped-it first gradually decreases, then remains stable, and finally increases suddenly. The results of a comparison among the creep failure stress level, creep duration, and creep rate of rock samples in different saturated states under the action of freezing and thawing are shown in Table 6. The rock samples with a saturation of 30%, 50%, 70%, 90%, and 100% fail at stresses of 15.47 MPa, 13.05 MPa, 9.46 MPa, 8.2 MPa, and 4.61 MPa, respectively; at low levels of saturation, the creep time gradually decreases with increases in saturation. When the saturation is greater than 70%, the creep time gradually increases; the creep rate of the rock sample increases with increases in saturation.  8 Geofluids the saturation is 30%, only one primary crack appears. As the saturation increases, the failure cracks of the rock sample gradually increase, and microcracks appear. The crack spacing of the rock sample gradually forms, and tension characteristics manifest. This is because, as the saturation increases, the freeze-thaw cycle causes the cementation between the particles of the rock sample to weaken, increasing cracks. The internal cracks in the rock sample structure gradually develop and expand under long-term loading, eventually evolving into multicrack shear failure [28].

Geofluids
Burgers model can fully reflect the decelerated-creep and stable-creep stage during the creep test, but it cannot describe the accelerated-creep stage. Therefore, in this paper, a modified viscous-plastic body is connected in series on the basis of the Burgers model. By introducing the damage variable to consider the accumulation of damage in the creep process, describing the variation of the viscosity coefficient with the creep process and accounting for the effect of saturation on the viscosity coefficient, a freeze-thaw-damage viscous element was constructed. The viscous components are shown in Figure 10.
The freeze-thaw-damage viscous element is a Newtonian body, and its constitutive relationship is shown in where σ is the stress of the viscous element, _ ε is the strain rate, and η 3 ðN ω , D t Þ is the viscosity coefficient of the freeze-thaw-damage viscous element.      Figure 12: Creep verification curve of 70% saturation rock samples. Geofluids Considering the influence of saturation and long-term loading on the viscosity coefficient, the coefficient of the viscous element is shown in where η 3 ðN ω Þ is the viscosity coefficient after a certain number of freeze-thaw cycles and D t is the damage variable (0 ≤ D t < 1). A significant amount of research has shown that the damage caused by long-term loading in the process of rock creep has a negative exponential function relationship with time [29], as shown in where α is the coefficient related to the material properties and saturation of the rock sample and t is time. Therefore, considering the joint influence of saturation and long-term load, the viscosity coefficient is where αðN ω Þ is the material coefficient of the rock sample at a certain number of freeze-thaw cycles at a degree of saturation. Therefore, the constitutive relationship of the viscous element is The model is composed of four parts connected in series: a water-containing elastic element, a water-containing viscous element, a water-containing viscoelastic element, and a water-containing viscous-plastic body. The creep model is shown in Figure 11, where E 1 , E 2 , η 1 , and η 2 are parameters related to saturation; η 3 is the parameter related to loading time and saturation. The creep equation is obtained as shown in Equations (8) and (9): when σ ≤ σ s , when σ ≥ σ s ,  Figure 12 and Table 8.

Conclusions
The freeze-thaw cycle has a significant influence on the longterm mechanical properties of rocks. However, most current creep research is only carried out in the fully saturated state. This paper conducts NMR, scanning electron microscopy, and uniaxial compression creep tests on red sandstone in different states of saturation after freeze-thaw cycles. The main results of this research are summarized as follows.
Under the action of freeze-thaw cycles, as the saturation of rock samples increases, the T 2 spectrum distribution shifts to the right, the total area of the T 2 spectrum gradually increases, and the porosity increases. Further, when the saturation increased from 70% to 90%, the distribution characteristics of pores changed significantly. After freeze-thaw cycles, mesopores replaced micropores, becoming the main pore form.
By SEM image, it can be seen that the number of pores in the rock samples gradually increases with increasing saturation levels. When the saturation reaches 70%, freeze-thaw cycle damage gradually appears, the number of dissolution holes increases, and penetration fracture occurs in the particle.

Geofluids
Under uniaxial compression, rock samples in different saturated states experience the compaction stage, linearelastic stage, and plastic stage. The instantaneous mechanical properties of rock samples are greatly affected by saturation.
The critical value appears when the saturation increases from 70% to 90%.
Under long-term loading, saturation has a significant influence on the time-efficiency characteristics of the rock sample during the freeze-thaw cycle. With increases in saturation, the creep deformation gradually increases. When it reaches 70%, the axial creep strain increases greatly, the creep rate increases, the creep failure stress decreases, and the creep time under the final level of stress increases significantly. This shows that the saturation level of 70% is the critical saturation for the deterioration of the mechanical properties of the rock sample.
A modified viscous-plastic body is connected in series on the basis of the Burgers model. The damage accumulation and saturation change in the creep process are fully considered. The freeze-thaw-damage viscous element-with the viscosity coefficient η 3 ðN ω , D t Þ-is introduced, and creep parameters were fitted using the experimental data.

Data Availability
The data used to support the findings of this study are included within the article, displayed by figures and tables.

Conflicts of Interest
The authors declare no conflict of interest.