Effects of Freezing and Thawing Cycle on Mechanical Properties and Stability of Soft Rock Slope

1State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining & Technology, Xuzhou 221116, China 2School of Mechanics and Civil Engineering, China University of Mining & Technology, Xuzhou 221116, China 3School of Mines, China University of Mining & Technology, Xuzhou 221116, China 4School of Resources and Geosciences, China University of Mining & Technology, Xuzhou 221116, China 5Department of Mining Engineering & Metallurgical Engineering, Western Australian School of Mines, Curtin University, Kalgoorlie, WA 6430, Australia


Introduction
Slope is a steep hill with a certain angle, and it is composed of soil and rocks.It is the most common rock engineering structure [1].Open-pit mining forms slopes with big geometric size.Safety and stability of open-pit mine slope are of the important factors that restrict safe and high-efficient production of the mine [2][3][4][5].Deformation and failure of big slopes cause huge economic losses and even threaten life safety of workers [6][7][8][9][10][11][12].The shear strength of slope rocks is a basic parameter that determines slope stability, which is influenced greatly by external factors.Temperature fields affect mechanical properties of rocks significantly, especially for rocks with high moisture content and low stiffness [13][14][15][16][17][18][19].In China, coal resources are mainly in northern areas which have more than 10 open-pit coal mines with an annual output of over 20 Mt.These regions belong to seasonal freezing soil areas.Due to temperature rise and drops during season changing, rocks experience physical freezing and thawing cycle, which is easy to cause various slope failures, such as frost heaving, rock mass denudation, and creep deformation.Freezing and thawing cycles cause the development of rock fractures, degrade physical and mechanical properties of rock-soil mass, reduce slope strength, and thus easily result in landslide accidents [20,21].
Many scholars have made long-term researches on freezing and thawing problems of slope.A set of unique mechanical theoretical systems based on damage mechanism of rock and soil have been developed [22].The variation of physical and mechanical properties of rock with defects (e.g., joint fissure) under freezing and thawing cycle and phase changes and moisture and heat transfer characteristics caused by freezing of free water inside the slope were investigated [23].The evolution of macro-and microdamage of rocks caused by different temperature histories and moisture (ice) contents was also studied.Such mechanical theory provides theoretical guidance to the design and construction of rock engineering in cold regions [24].
Existing theoretical and experimental researches concerning the effect of low temperature range (generally less than 0 ∘ C) on physical and mechanical properties of rocks could be divided into three types.The first type, which is based on continuum mechanics and the classical heat transfer theory, studies thermal-hydro-mechanical (THM) coupling of rocks under low temperature and the freezing and thawing cycles [25], as well as distribution laws of frost heaving pressure influenced by frost heaving and temperature on tunnel surrounding rocks in cold regions [26,27].The second type studies basic mechanical properties of rocks under freezing and thawing cycles based on damage mechanics.Nicholson et al. [28] analyzed damage deterioration of 10 kinds of rocks caused by freezing and thawing cycles and studied the macroscopic freezing and thawing damage evolution of rocks through graphic record.Li and others [29] studied freezing and thawing microscopic damage of rocks by using advanced Computed Tomography (CT) scanning technique and attempted to establish a damage constitutive model that uses CT as the freezing and thawing damage variable.Based on related theory of fracture mechanics, the third type studies crack development and failure criterion of cracked rocks under freezing, thawing, and the freezing and thawing cycles [30,31].
China's coal resources are mainly located in the northern areas which belong to the seasonal frozen soil area and experience complete freezing and thawing cycles annually.Due to the mining conditions, open-pit is the preferred mining method in this area.Therefore, a great number of high steep slopes are constructed.Few researches on slope engineering of open-pit mines explore the influence mechanism of freezing and thawing cycles on slope deformation and stability.In seasonal freezing soil areas, strength, moisture content, deformation features, and mechanical parameters of soft rock slope are significantly different than other areas.Therefore, researches on deformation and stability of slope under freezing and thawing conditions shall be based on its mechanical states and corresponding physical and mechanical parameters.Long-term engineering practices have already proved that the freezing and thawing cycles play an important role in influencing the stability failure of slopes.As a result, studying frost heaving characteristics of rocks, physical and mechanical parameters, and stability of soft rock slope during the freezing and thawing cycles has very important engineering value and safety significance to open-pit mines in North China.

Introduction of Wulin Open-Pit Mine.
The Wulin openpit mine lies in Linkou County, North China (Figure 1).It is surrounded by railways and highways and thus has convenient transport conditions.The whole ore field, in average, is 5.

Effect of Freezing and Thawing Factors.
The area covered within the Wulin mine belongs to the typical seasonal freezing soil area.The surface soil begins to be frozen gradually after November.The freezing soil thickness increases gradually until reaching the maximum.Then, the surface and the maximum underground freezing depth begin to thaw in March.The freezing layer between the surface and the maximum underground freezing depth is called the seasonal freezing layer.November is viewed as the typical freezing period of the mine, while March is viewed as the typical thawing period.Time from November to March is viewed as the frozen period.The seasonal freezing layer experiences both freezing and thawing in one year, thus producing frost heaving and thaw collapse of unfrozen water in soil.Freezing causes heaves at some positions of the slope, increases the overall stiffness, and changes the slope flatness.Thawing causes rough and uneven rock-soil mass on the slope surface and working benches as well as muddy road surface and thus influences production efficiency.After entering into the thawing period, different scales of landslides on the superficial slope layer frequently occur, threatening safety of workers and equipment.Freezing and thawing are the key influencing factors of slope stability in this mine.

Experiments
3.1.The Shearing Experimental System.The whole test system structure is shown in Figure 2. The system is composed of the following three parts.
(1) Multifunction Frozen Soil Simulation System.The frozen soil test system has good heat insulation performance.The  inner wall has a 150 mm thick polyurethane insulation coating and the thermal insulation material on the external wall could reduce temperature influence and ensure high accuracy of test results.
(2) Cooling and Temperature Control System.The cooling system and temperature control system provide the experimental environment.The used cooling system consists of two independent cooling devices, which could provide low temperature environment.The control accuracy and display accuracy are 0.5 ∘ C and 0.1 ∘ C, respectively.
(3) Shear Test System.The test system is equipped with one EDJ-1 strain-controlled direct shear apparatus.The diameter of the cutting ring is 61.8 mm, the height is 20 mm, and the shear velocity is 0.8 mm/min.Since metals are sensitive to temperature, temperature change will directly influence technical parameters of the dynameter ring.Effect of temperature on the dynameter ring was corrected through three times of prepressing with the maximum load under different temperature environment.The correcting algorithm is where  is the corrected value of the dynameter ring (kPa/0.01mm),  0 is the original value of the dynameter ring (kPa/0.01mm), t is the service temperature of the dynameter ring ( ∘ C),  0 is the setting temperature of the dynameter ring ( ∘ C), and  is temperature correction factor (0.0003/ ∘ C).

Shear Test Design and Procedures
(1) Sample Preparation.The specimens were collected from the south end-wall slope of Wulin open-pit mine, and the principal component of the rock specimen is siltstone with the grain size of 0.1 mm.It is known the average unit weight of the rock specimen is about 18.3 kN/m 3 , and the natural moisture content is about 15.3%.In the laboratory, the specimens were made into the standard samples with the diameter of 61.8 mm and 20 mm high, as shown in Figure 3.
It is important to make sure the rock specimens used in the laboratory have the similar properties with the natural state, so that the experimental results found in the presented research could be representative of the real behavior experienced by the materials under consideration.In addition, the required amount of water will be added to the rock specimen during its preparation, so that the moisture content is the same with the original state.
(2) Test Conditions and Steps.Using the control variable method, this experiment prepared freezing rock specimens by changing the initial moisture content, number of freezing and thawing cycle, and freezing and thawing temperature one by one.Then, mechanical strength of prepared specimens was tested through the direct shear test.Based on the statistical

Shear Test Data Processing and Analysis.
To explore the influence of key factors on the strength of frozen rock specimens, several groups of freezing and thawing experiments were conducted with rock specimens collected from the Wulin mine.The orthogonal experimental design was employed to reduce the amount of the experimental work.Experiment results and mechanical parameters of soil sample with 23% water content are shown in Figures 4 and 5 principle of three factors and four levels.The experimental design and the test data are shown in Table 1.
It can be seen from Table 1 that the moisture content and number of the freezing and thawing cycles significantly influence mechanical parameters of the frozen soft rock, while the freezing and thawing temperature is less significant.Therefore, this paper focused on studying the influence of initial moisture content and number of the freezing and thawing cycle.Relationship curves of the initial moisture content (Figure 6) and the number of the freezing and thawing cycles (Figure 7) with mechanical parameters of the frozen soil were gained through additional experiments.
Figure 6(a) reveals that cohesion of the frozen soil sample presents a linear function decrease with the increase of moisture content, which follows the formula as where  is cohesion;  and  are the parameters associated with the temperature difference, where  = 0.008Δ + 1.376  and  = 0.095Δ + 52.77;   is moisture content; Δ is the temperature difference.Such decrease slows down gradually.This is mainly because the increase of moisture content thickens the thin layer of water absorbed by particles in the rock specimens, which increases the distance between internal particles and weakens their connection strength.
In Figure 6(b), a quadratic function decrease of internal friction angle is observed with the increase of moisture content.Such decrease accelerates gradually, which can be described as follows: where  is internal friction angle; , , and  are the parameters associated with the temperature difference, where  = −0.00063Δ+ 0.0355,  = −0.023Δ+ 1.229, and  = 0.189Δ + 10.74;   is moisture content; Δ is the temperature difference.
When pores between particles are filled with water, the thickness of water film increases and the effective contact area of particles decreases, so that the friction coefficient reduces correspondingly.According to test results, frozen soil sample has a critical moisture content.Once exceeding this critical value, the internal friction angle of the soil mass will decrease quickly with the further increase of moisture content.
In Figure 7(a), cohesion of specimens reduces exponentially with the increase of the freezing and thawing cycle.During the freezing period, ice crystals produced in pores destroy internal connection of soil particles and thereby weaken the soil mass structure.Several cycles of freezing and thawing will destroy the internal connection of particles, resulting in continuous decrease of cohesion.In Figure 7(b), the internal friction angle decreases exponentially with the increase of freezing and thawing cycle.The freezing and thawing cycles damage the original porous features in specimens, thus changing the internal skeleton.After the freezing and thawing cycles, both contact points and the friction between internal particles are reduced, thus decreasing the internal friction angle.It can be seen from Figure 7 that the reduction of cohesion and internal friction angle slows down significantly after three freezing and thawing cycles.This demonstrates that the first three freezing and thawing cycles could weaken cohesion and internal friction angle of rock and soil greatly, but they will become stable as the number of cycles further increase.

Computational Simulation Analysis
Temperature field and stress field of the soft rock slope are sensitive to temperature changes.By using the FLAC software, the influences of these factors on the slope stability    of the Wulin mine during freezing and thawing cycles were analyzed.A slope model with 1 : 1 scale of the actual geometric size was established.The slope height and angle are 10 m and 65 ∘ , respectively.The whole model was divided into 1,125 grids (Figure 8).The slope model has displacement constraint and temperature constraint.For displacement boundary condition, the two sides use normal constraint and the bottom use full displacement constraint.For temperature boundary condition, the upper slope surface chooses the first type of boundary condition and temperature is set as the minimum air temperature of the month.In other words, temperature gradient is set as the maximum temperature difference of the month.Sides of the slope are set as thermal insulation conditions, while model bottom considers terrestrial heat effect and is set constant.The involved parameters used in the numerical simulation model are listed in Table 2.
According to the on-site slope structure and tested physical and mechanical parameters of freezing and thawing rocks, a simulation analysis of the slope stability during different stages of the freezing and thawing cycle was conducted.This simulation disclosed the evolution of slope stability coefficient during the freezing and thawing cycle.The original rock gravity (), internal friction angle (), cohesion (C), moisture content, slope angle (), and height (H) are 18.3 kN/m 3 , 21.3 ∘ , 34.5 kPa, 15.3%, 65 ∘ , and 10 m, respectively.Based on the associated flow rule of FLAC that is commonly used to analyze the stability of slope, it is concluded that the stability coefficient of the original slope is 1.51, as shown in Figure 9.

Freezing Period.
To accurately disclose the influence of a single factor on the freezing and thawing slope stability, the variable-controlled approach was applied.Firstly, the moisture content of the freezing region was kept as 15.3% and the maximum thickness of the freezing soil layer was set to 2 m.The physical mechanical parameters of normal and freezing region are as shown in Table 3.The change of the slope stability coefficient with the freezing depth (0.5 m, 1 m, 1.5 m, and 2 m) is analyzed as shown in Figure 10.

Thawing Period.
Influence of thawing layer thickness (0.5 m, 1 m, 1.5 m, and 2 m) on slope stability was analyzed.Mechanical parameters were chosen according to test results under 15.3%, 18.10%, 20.30%, and 23% moisture contents.The change of the slope stability coefficient with the thawing depth is shown in Figure 11.
According to Figure 11, the stability coefficient is inversely proportional to the thawing depth.Under 15.3% and 18.1% moisture content, the slope stability coefficient shows almost a linear decrease with the increase of thawing depth.Under   20.3% and 23% moisture content, the slope stability coefficient decreases exponentially with the increase of thawing depth.It reveals that higher thickness of the thawing area not only reduces the stability coefficient but also changes the shape and region of the slip mass.The thicker the thawing layer, the larger the area with the risk of sliding surface in the thawing region.

Effect of Number of the Freezing and Thawing Cycles on
Slope Stability.Rock structure as well as physical and mechanical parameters of soft rock slope deteriorate significantly after the freezing and thawing cycle.Each freezing and thawing cycle deteriorates slope rock and stability coefficient to a certain extent.Relationships between the number of the freezing and thawing cycles (1, 2, 3, 4, and 5) and slope stability coefficient are analyzed as shown in Figure 12.
Seen from Figure 12, after experiencing several freezing and thawing cycles, given the same freezing depth and moisture content, the slope stability coefficient shows a negative exponential decrease as the freezing and thawing cycles increase, which discloses an important engineering problem.After the slope is formed, the first freezing and thawing cycle deteriorates most of the overall slope stability.Once the number of freezing and thawing cycle exceeds 3 times, the slope stability coefficient decreases slightly and tends to be stable.

Field Investigation
In Wulin surface mine, the mining direction is from west to east; thus, the high and steep end-wall slopes were formed at both north and south sides of the mine.The slope is located at the eastern part of the end-slope within the scope of 80 m.Under the influence of unloading produced by excavation, the slope appeared as an obvious deformation along the weak layer, as shown in Figure 13.The fragmentation degree and cracks of soft rock in the slope were significantly increased after experienced a freezing and thawing process.This phenomenon embodies the effect of freezing and thawing cycle action on the deformation of newly formed slope.
To estimate the deformation of slope in the mining process, three sets of monitoring equipment were set up on the south end-slope and the west slope.There are three sets of monitoring equipment arranged at the same height, which were named as P2, P4, and P6.In addition, the slopes monitored by P2 and P4 existed for 5 years and 2 years, while the slope under the monitor of P6 is a new slope.The displacement data of the three monitoring points from November 2015 to March 2016 is shown in Figure 14.
Seen from the monitoring data in Figure 14, after one freezing and thawing cycle, the largest displacement appeared at P6, and the smaller displacement appeared at P4. Furthermore the minimum displacement is appeared at P2 after it experienced five times freezing and thawing cycles.The slope displacement distribution is the same as the slope stability changing rule with different times of freezing and thawing cycles.This is why slopes that have existed for years maintain basically the same stability with the increase of freezing and thawing cycle.However, open-pit mines with continuous developing slopes are highly sensitive to freezing and thawing cycle and deserve higher attentions.
The experimental research has concluded that the moisture content as well as the number of the freezing and thawing cycles significantly influences the strength of freezing soft rock specimens.However, there is no explicit conclusion concerning the influence degree of these key factors on rock strength and stability of the soft rock slope.Test results demonstrate that mechanical parameters of soft rock and slope stability coefficient decrease in a negative exponential pattern with the increase of freezing and thawing cycle.But after three freezing and thawing cycles, the slope stability is only influenced to a lesser extent with the further increase of the freezing and thawing cycle.Moisture content significantly affects mechanical parameters of soft rock and slope stability all the time.The higher the moisture content is, the more seriously the soft rock strength and slope stability are influenced.Based on these influence laws, when the number of the freezing and thawing cycles is smaller than three times, the soft rock slope stability is mostly affected.At the same time, when the number of the freezing and thawing cycles is larger than three times, the moisture content influences the soft rock slope stability mostly.To determine the dominant influencing factor of soft rock slope stability, the combination effect of specific rock properties and parameters as well as the number of freezing and thawing cycle the slope have experienced need to be considered.Further researches will focus on the coupling relationship of these two influencing factors, their relative importance, and governing strategies of soft rock slope in seasonal freezing soil areas.

Conclusions
The freezing and thawing cycles affect rock strength of soft rock slope significantly.The experiment shows that moisture content and number of the freezing and thawing cycles are the two important influencing factors for mechanical parameters of soft rock.Their functional relationships with basic mechanical parameters of soft rock are gained.This paper concludes the following: (1) The cohesion of frozen soil specimens presents a quadratic function decrease with the increase of moisture content.Such decrease slows down gradually.The internal friction angle of the frozen soil specimens shows a negative exponential decrease against the increase of moisture content.Such decrease accelerates gradually.As the freezing and thawing cycles increase, both cohesion and internal friction angle of frozen soil specimens decrease exponentially.Such decrease slows down significantly after three freezing and thawing cycles.
(2) Slope stability coefficient is proportional to freezing layer thickness.With the increase of thawing layer thickness, slope stability coefficient decreases linearly when moisture content of the soft rock is smaller than 20.3% but decreases exponentially when moisture content is larger.
(3) After experiencing several freezing and thawing cycles, given the same freezing depth and moisture content, the slope stability coefficient shows a negative exponential decrease as the freezing and thawing cycles increase.After three freezing and thawing cycles, the slope stability coefficient decreases slightly and tends to be stable.Therefore, strip mines with continuous developing slopes are highly sensitive to the number of freezing and thawing cycles.This the key technical problem of strip mines in seasonal freezing soil areas.

Figure 1 :
Figure 1: Location of the study area.
(a) Frozen soil simulation system (b) Cooling system (c) Temperature control system (d) Shear test system

Figure 6 :
Figure 6: Relationship curve between initial moisture content and mechanical parameters.

− 5 ∘− 5 ∘Figure 7 :
Figure 7: The change of cohesion and internal friction angle with the freezing and thawing cycles.

Figure 8 :
Figure 8: Simulation results of temperature field and frost heaving force for slope.

Figure 10 :
Figure 10: Variation of stability coefficient with freezing depth.

Figure 11 :
Figure 11: Variation of slope stability coefficient with thawing depth.

Figure 12 :
Figure 12: Variation of slope stability coefficient with number of the freezing and thawing cycles.

Figure 13 :Figure 14 :
Figure 13: The deformation and fragmentation of new slope after undergoing the freezing and thawing cycle process.
2 km long and 0.82 km wide, covering an area of 4.3 km 2 .The mine lot belongs to hilly slope in terrain and valley plain in landform, showing a relatively simple structure.The lowest temperature of the covering area is in January (−30.3 ∘ C, −17.7 ∘ C in average) and the highest temperature is in July (29.7 ∘ C, 22.6 ∘ C in average).The maximum freezing depth of soil is 2 m.

Table 1 :
Orthogonal experiment design and results.

Table 2 :
The involved parameters used in the numerical simulation model.

Table 3 :
Physicomechanical parameters of normal and freezing region.