Abstract

In cold regions, rock’s load-bearing capacity will be greatly diminished in severe settings. This research investigates the influence of freezing and thawing on the physical and mechanical characteristics of sandstone, as well as the strength under complicated stress conditions. Uniaxial compression of sandstone samples was performed following freeze-thaw cycles, and the changes in elastic modulus and peak stress of dry sandstone and saturated sandstone were investigated under various freeze-thaw cycles. ANOVA was used to analyze whether there were significant differences between the number of freeze-thaw cycles and the peak strength and elastic modulus. A three-parameter strength prediction model based on σC, k0, and m is created based on the critical failure energy function. Experimental data is utilized to assess the model’s correctness, and the model is used to forecast the strength of dry and saturated sandstone during freeze-thaw cycles. The result indicates that peak stress and elastic modulus of dry and saturated sandstone show a steady attenuation pattern as the number of freeze-thaw cycles increases. Brittle failure is the failure mode of dry sandstone, whereas brittle failure to plastic failure is the failure mode of saturated sandstone. After 60 freeze-thaw cycles, peak stress in dry and saturated sandstone was reduced by 55.3% and 56.8%, respectively. The three-parameter model can predict the triaxial compressive strength of specimens under different confining pressures through the uniaxial compressive strength of specimens with an error range of 1% to 13%.

1. Introduction

Rock stability in the process of tunnel excavation in cold regions has long been a hot topic. The mechanical parameters of a rock tunnel are influenced by damage within the rock [13]. In fact, the freeze-thaw (F-T) cycle, as well as moisture content, temperature, and other factors [4, 5] will degrade the rock, resulting in tunnel instability and accidents. In general, under natural conditions, a certain amount of moisture in the rock will improve its plastic performance [6]. However, as a result of the freezing and thawing cycles, additional cracks will result in damage [7]. The growth of interior fissures becomes more perfect as the freezing and thawing times rise, and the bonding force between particles is increasingly weakened by hydration and frost heaving force, diminishing the bearing capacity of rock. F-T action will have a direct impact on tunnel engineering safety [810]. As a result, significant research topics in tunnel engineering are the mechanical properties of rock under the operation of the F-T cycle.

When the water in the pores of the rock turns from liquid to solid, the volume expands by 9%, and the constant migration of water molecules causes the rock damage to worsen, which may cause slope instability [1114]. Fan [15] investigated the features of stress variations over time in a F-T cycle sandstone loading experiment. Based on the combined action of F-T and load. Zhou et al. [16] looked at the sandstone from a microscopic perspective. The porosity rose as the F-T periods increased, and a link between porosity and peak strength was discovered. Gasc Barbier et al. [17] analyzed the crack propagation of rock mass after F-T cycles from the perspective of elastic waves. Zhang et al. [18] found a link between the number of F-T cycles and rock tensile strength (UTS) and anisotropy and discovered that as the number of F-T cycles increased, so did the anisotropy of rock tensile strength. Zhang et al. [19] developed a damage equation based on several types of damage and discovered that increasing the initial load reduces the damage produced by F-T. Jihwan Park et al. [20] investigated the impact of F-T cycles on internal rock structural alterations. Water caused interior particles to loosen and new cracks as a result of volume expansion. They concluded that F-T damage had less of an impact on the deterioration of granite with thick pores. Ke et al. [21] carried out different times of F-T cycles on sandstone to study the characteristics of internal pore development. The pores changed from small pores to medium pores and large pores, and the connectivity of pores led to the decrease of peak stress. Huang et al. [22] monitored the pore structure changes of sandstone under F-T by MIP (mercury intrusion porosimetry). Nanopores and micropores had gradually developed into macropores with F-T. The results showed that, with the increase of F-T cycles, the cohesion of sandstone decreased and the internal friction angle increased. Liu et al. [23] used the improved fatigue damage model to study the F-T degradation mechanism of clay-bearing sandstone under unsaturated state and predicted the strength of clay-bearing red sandstone under different F-T cycles. Kahraman and Ozdemir [24] used four different machine-learning techniques to predict the uniaxial compressive strength (USC), specific gravity (SG), water absorption rate (WA), Los Angeles abrasion (LA), P-wave velocity (Vp), and freeze-thaw resistance (F-T) tests. The results showed that machine-learning could be successfully used to predict F-T resistance. Liu et al. [25] took F-T damage into account and used P-wave velocity to replace other parameters of rock to establish a model to predict UCS of rock with different F-T cycles. However, P-wave velocity cannot precisely forecast the strength of rock because of the anisotropy in the rock. Although the existing research results have analyzed the influence of F-T cycle on rock physical and mechanical parameters in detail, the effect of F-T on each physical and mechanical parameters of sandstone is different. Some scholars only study the influence of F-T cycle on a single physical and mechanical property. They do not further analyze which factors were more sensitive to F-T. In fact, it is very important to find out which mechanical parameters are more sensitive to the F-T cycle so as to take some measures to increase the stability of rock. The strength prediction model based on the mechanical parameters obtained by the uniaxial compression experiment can decrease the strength prediction error caused by rock anisotropy.

In this paper, combined with the effect of F-T on the physical and mechanical properties of sandstone, the peak strength and elastic modulus of sandstone with different F-T cycles were tested, and the variance analysis method (ANOVA) in statistics was used to quantitatively analyze the influence of F-T factors on the peak strength and elastic modulus. Based on the critical failure energy function, the compressive strength of sandstone under complex stress is reasonably predicted.

2. Materials and Methods

2.1. Preparation of Sandstone Samples

The rock samples in this experiment were taken from the same intact rock block without joints, and 20 sandstone samples of Φ50 mm × 100 mm were prepared according to the standards of the International Society of Rock Mechanics. During cutting, the unevenness of the upper and lower ends of the rock sample was not more than 0.05 mm, and the height error of each rock sample was within 0.3 mm. Ultrasonic velocity test was carried out on the samples with complete appearance, and 20 samples with similar quality were selected and divided into 10 groups (2 samples in each group). The first five groups were dried in the oven to remove the excessive moisture in the sample, and the latter five groups were treated with vacuum saturation method to make the sample fully saturated. The preparation of dry and saturated samples is shown in Figure 1.

Firstly, 20 cut samples were put into the oven and dried at a constant temperature of 105° for 48 h. Ten remaining samples were put into the vacuum water saturation device and stored for 48 hours to prepare the saturated samples. The air pressure in the device was controlled at −0.1 MPa during the whole process. The saturated moisture content of the sample is about 3.1%.

The obtained dry and saturated samples were subjected to 5, 10, 20, 40, and 60 F-T cycles at −20°C and 20°C. The sample numbers are shown in Table 1, and the process is shown in Figure 2.

2.2. Test Method and Process

In this experiment, TAW-2000 electrohydraulic servo rock triaxial test instrument was used for uniaxial compression test. The maximum test force of the instrument was 2000 kN. The load and axial deformation of the sample during compression were directly measured by the instrument display system. In the process of uniaxial compression, in order to obtain the full stress-strain curve of the specimen, displacement was controlled at the rate of 0.05 mm/min when loading [26].

3. Results and Discussion

3.1. Physical and Mechanical Properties of Samples after F-T Cycles
3.1.1. Variation of Stress-Strain Curves of Sandstone Samples

Based on the uniaxial compression test after F-T cycles, the stress-strain curves of dry and saturated sandstone after 5, 10, 20, 40, and 60 F-T cycles are recorded as shown in Figure 3.

With the increase of F-T cycles, the peak stress of sandstone decreases gradually. The peak stress of dry sandstone decreases obviously by 11.13 MPa after 40 F-T cycles. The peak stress of saturated sandstone decreased by 13.83 MPa after 10 F-T cycles. After 5 F-T cycles, the peak strength of dry rock samples decreases by 55.3% and that of saturated rock samples decreases by 56.8% after 60 F-T cycles.

In the early stage of loading process, the sandstone is in the pore compaction stage, and the curve is characterized by concave growth. As the skeleton of the pore compacted sandstone begins to bear the load and the force is within the elastic range, the stress-strain curve is nearly a straight line. When it reaches the yield stress, the sandstone particles produce vertical slip, and the internal microcracks gradually grow and plastic deformation occurs. When the stress reaches the peak value, the microcracks generated in the yield stage are connected to each other to form a fracture surface. At this point, the strain continues to increase, and the instantaneous decrease of stress indicates the failure of sandstone.

3.1.2. Effect of F-T Cycles on Peak Stress and Elastic Modulus of Sandstone with Different Water Content

With the increase of F-T cycles, the stress-strain slopes of all samples gradually slow down. Under the action of freezing and thawing, the thermal expansion and contraction of the cementing material in the dry sandstone lead to the continuous change of pore volume and the gradual decrease of the maximum tension. With the increase of F-T cycles, the pores in the sandstone gradually lose elasticity, resulting in the decrease of the peak strength and elastic modulus of the sandstone. For saturated sandstones containing internal water, freezing and thawing under the action of water-ice phase change result in greater pore volume expansion. The presence of water also impairs the interaction force of internal particles and the degree of compactness, resulting in loss of sandstone peak strength. Since the water content of saturated sandstone is higher than that of dry sandstone, and the volumetric deformation produced in the freezing process is greater than that of dry sandstone, therefore, the bearing capacity of dry sandstone will change significantly only after considerable F-T cycles while the saturated sandstone is sensitive to F-T action and the peak stress decreases quickly even after several cycles. After 60 F-T cycles, the peak stress of saturated sandstone is 1.5% higher than that of dry sandstone.

The elastic modulus of sandstone changes as shown in Figure 4.

As the number of F-T cycles increases, the elastic modulus of dry sandstone decreases linearly while that of saturated sandstone decreases exponentially. The elastic modulus of dry sandstone decreases by 52.8% and that of saturated sandstone decreases by 69.7% after 60 F-T cycles.

With the increase of F-T cycles, the internal structure of dry sandstone changes almost uniformly, and the evolution rate of pore throat is close to uniform, which leads to a certain linear relationship between macroscopic mechanical properties and F-T cycles. In the saturated sandstone, due to the existence of pore water, the pore development is random under the action of F-T. In the initial cycles, the pores in the saturated sandstone expand sharply and thus the elastic modulus decreases rapidly. However, with the internal pores connected, the elastic modulus of the sandstone gradually tends to be stable, and no large cracks but only some small cracks will be produced. Hence, the elastic modulus of the sandstone decreases slowly. Since the water content of saturated sandstone is higher than that of dry sandstone, the internal damage is more severe, and the elastic modulus of saturated sandstone is 6.9% lower than that of dry sandstone after 60 F-T cycles.

3.1.3. Statistical Analysis

The statistical methods for analysis are designed to determine whether moisture content has a significant effect on sandstone peak strength and elastic modulus. This section delves into the details of the issue. Because the average value of the drop in peak strength and elastic modulus after 60 F-T cycles follows the normal distribution in statistics, the number of F-T cycles is used to group the samples and determine whether there are significant variations between groups.

An analysis of variance (ANOVA) was conducted to investigate if there were statistically significant differences in peak strength and elastic modulus between sandstone sample groups with different F-T cycles. When using single factor analysis of variance to study sandstone samples with 100% saturation, there is always a control group for comparison.

The original hypothesis in this experiment is that there is no significant difference in the average μ between groups, indicating that μ1 = μ2 = μ3. Assuming that the explicitness level α = 0.05, the total amount of samples is defined as N, N is divided into r different groups according to different moisture contents, for instance, sandstone samples with five F-T cycles consider three different moisture contents of 0, 2.262%, and 3.77%, each moisture content group contains 10 samples, r is equal to 10, and N is equal to 30.

In the process, the following procedure must be followed [27].

According to F statistics,where N is the total number of sandstone samples; R refers to groups by number of F-T cycles;  is the average value of each group; SSE is the intragroup variation; SSA is the intergroup variation.

Thus, if F ≤ Fα=0.05, the original assumption will not be rejected; then, µ1 = µ2 = µ3, which indicates that the differences between the means are statistically insignificant.

If F > Fα=0.05, the hypothesis is rejected, then µ1 = µ2 = µ3, and thus the means are significantly different. Table 2 summarizes the results of peak strength and elastic film of saturated sandstone samples with five different F-T cycles. Table 3 presents a summary of the statistical analysis.

3.2. Critical Failure Energy Function

The M-C strength criterion establishes the mathematical relationship between shear stress τ and normal stress σ through the envelope of the limit Mohr circle and accurately describes the ultimate stress state of each point of the fracture surface to a certain extent. Assuming that the friction coefficient between the fracture surfaces is constant, the limit equilibrium state of sliding friction just meets the M-C strength criterion. When rock mass is subjected to force, longitudinal and lateral deformations are generated macroscopically. These two kinds of deformations lead to the destruction of rock mass, namely, shear failure of rock mass and expansion failure of rock mass. According to the functional principle, when the rock mass is unstable, the critical failure energy required for two forms of deformation on the fracture surface is provided by the shear stress and normal stress at this point. The potential energy stored by the shear stress τp in the fracture surface under elastic shear strain γp as well as the potential energy stored by the normal stress σp under volumetric strain can be expressed as follows [28]:where e is the strain energy when the unit breaks; μ and K represent the shear elastic modulus and bulk modulus when the material is damaged; is volume stress:  = σ1 + σ2 + σ3. When loading, the positive stress prevents the formation of rock fracture, so it provides negative energy.

In three-dimensional principal stress space, (σ1, σ2, σ3) is designated as maximum principal stress, intermediate principal stress, and minimum principal stress, respectively. Assuming that the direction of the crack plane is (l, m, n), the shear stress τp and the normal stress σp can be expressed as

According to the criterion of geometric space vector,

Formulas (4)–(6) are substituted into equation (3) to obtain

The parameter k is defined as

When the rock is unstable, the energy e reaches the peak value and cracks are produced. Through , , the cracking direction under the critical energy state is solved:

Substituting equation (9) into (7), the failure energy ed under the critical strain at one point of the fracture surface is

When the rock mass is destroyed, the failure energy ed acts on both sides of the crack in the form of surface energy. According to the law of conservation of energy, it is known thatwhere γs is the unit surface energy of the fracture surface, and the sum of surface energy on both sides of the fracture is 2γs. When pure shear failure occurs on the rock  = 0, equation (11) can be simplified aswhere τ0 is the shear strength of rock. After substituting formulas (12) into (11), we get

Formula (13) describes the critical failure stress state of rock based on the M-C criterion model. When the rock sample is under uniaxial compression, the parameter k = k0, according to formula (13):where σC represents UCS. Using σC to denote the offspring of τ0 into equation (15), we getwhere and represents a reduction factor for Poisson’s ratio at failure. It represents the condition of critical failure energy when σC is reached.

When the rock is subjected to triaxial compression, the volume expansion will occur before failure, and the volume deformation of rock is nonlinear. The deformation relationship between rock volume expansion and failure is expressed by power function:where A represents the empirical constant related to rock properties and σE represents the stress in volume expansion of rock: .

If equation (16) represents the linear relationship between volumetric stress and volumetric strain with respect to volumetric modulus K,

From formulas (8), (16), and (17),

In uniaxial compression σV = σCk = k0, the substitution of equation (18) can be obtained:

Formula (19) can be replaced by formula (15):

Formula (20) can be used to predict the strength of specimens under different confining pressures by determining the three parameters of σC, k0, m. It is a theoretical model of the power function relationship of the rock volumetric strain approximate to formula (16) when the failure occurs.

3.2.1. Model Verification

The uniaxial stress test is used to predict the compressive strength of rock under complex stress state, which is a convenient method for construction in practice. The following analysis is based on the research results published by scholars [29] to validate the accuracy of the theory.

Table 4 shows the physical and mechanical parameters obtained by uniaxial compression test on granite rock. The theoretical predicted strength values of granite under complex stress state can be obtained by the calculation formulas of k0 and m and the experimental data in Table 3. As the granite belongs to type II characteristic rock, Poisson’s ratio reduction factor of 0.2 is taken to calculate k0. The strength prediction results are shown in Figure 5. Through comparison, it can be found that the theoretical value of strength slows down with the increase of confining pressure. The predicted value of theoretical strength is slightly higher than the actual strength value, with the error of 1%–13%, indicating that the established prediction model of this theory can approximately predict the rock strength under complex stress state.

3.2.2. Strength Prediction

Based on the model established by the critical failure energy function, the strengths of sandstone under complex stress state after F-T cycles are predicted. The experimental results and material parameters of these sandstones are shown in Table 5.

The predicted strength of sand under complex stress state can be obtained by substituting the calculation formula of the parameters k0 and m of each sandstone and the experimental data shown in the table into equation (20). Since the stress-strain curve of sandstone belongs to type I, the reduction coefficient of Poisson ratio is 1/3 when calculating k0 and the strength prediction results are shown in Figure 6. σ1 of dry sandstone (Figure 6(a)) and saturated sandstone (Figure 6(b)) increases with the increase of σ3, and the increase rate of σ1 is rapid at first but slows down with the increase of σ3. When σ3 increases to a certain extent, the growth rate of predicted strength of sandstone tends to be flat. The model can well predict the ultimate strength of rocks in cold regions after F-T cycles.

3.3. Discussion

The UCS and elastic modulus of sandstone samples decrease to varying degrees after F-T cycles. The average peak strength decreases from 32.20 MPa to 14.41 MPa, and the average elastic modulus decreases from 4.53 GPa to 1.42 GPa. The F values of peak strength and elastic modulus calculated by ANOVA are 114.23 and 193.08, respectively, which are much larger than 3.48, indicating that F-T cycles can significantly affect the mechanical properties of sandstone, and the elastic modulus is more sensitive to F-T cycles. H. Yavuz et al. analyzed the effect of F-T cycle on peak strength and P-wave velocity by multiple regression analysis. They believed that the effect of F-T cycle on peak strength was more obvious [30], which was consistent with the conclusions of many scholars [31]. The results of Liu’s research showed that with the increase of confining pressure, the peak stress of coal rock increased, and the ductility of coal rock became more obvious [32], which further verifies the prediction results of the three-parameter model.

4. Conclusions

In this paper, five distinct F-T cycles were performed on two different types of sandstone with variable water content. The stress-strain curve was obtained by uniaxial compression of the specimen with variations in F-T cycles and water content. The relationship between elastic modulus, peak strength, and F-T times was investigated. The apparent differences between five distinct F-T cycles of saturated sandstone were compared using analysis of variance (ANOVA). Finally, the energy failure critical function was used to estimate the specimen’s strength under complex stress conditions, resulting in a strength prediction model.(1)Peak strength and elastic modulus are reduced as a result of the F-T cycle. Cementing force decreases as the physical properties of internal cementing materials change owing to F-T cycles. The internal water-ice phase transition causes the volume of the saturated sandstone to vary, and the frost heaving force between the particles causes the physical and mechanical properties of the entire sample to change.(2)The freezing-thawing cycle can drastically impair the peak strength and elastic modulus of sandstone samples, according to the results of variance analysis (ANOVA). Due to the filling of small pores by water molecules, water molecules migrate during the melting process, and new small cracks are generated between particles during the freezing process. Finally, the pores are joined to produce huge cracks, resulting in macroscale volume expansion as illustrated in Figure 7, which has a major impact on the samples’ peak strength and elastic modulus.(3)A three-parameter strength prediction model based on σC, k0, and m is created based on the critical failure energy theory, assuming that the shear modulus of elasticity is constant and the pressure of bulk modulus changes exponentially. According to experimental results reported by various scholars, the model’s error ranges from 1% to 13%. When using this model to forecast the strength of sandstone samples under complicated stress circumstances, the sample’s σ1 increases with the increase of σ3. The sandstone samples with more F-T cycles and saturated water rise more slowly, although their curves are generally flat.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors appreciate the National Natural Science Fund and the Liaoning Provincial Department of Education’s financial assistance for scientific research endeavors. This work was funded by the National Natural Science Foundation of China (51974145 and 51904144) and the Scientific Research Funding Project of Liaoning Provincial Department of Education (LJKZ0354).