Effect of Hydrothermal Coupling on Physical and Dynamic Mechanical Properties of Sandstone

Research Center of Mine Underground Engineering, Ministry of Education, Anhui University of Science and Technology, Huainan, Anhui 232001, China School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan, Anhui 232001, China State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mine, Anhui University of Science and Technology, Huainan, Anhui 232001, China Department of Civil, Construction, and Environmental Engineering, (e University of Alabama, Tuscaloosa, AL 35487, USA


Introduction
Heat-cool (H-C) cycle is a typical weathering process, which shows strong effect on the physical and mechanical properties of rock.Rock weathering is a phenomenon of rock breakage, porosity, and secondary change of mineral composition under the combined action of solar radiation, atmosphere, water, or biology [1,2].
e recurrent H-C cycles cause degradation of rock engineering [3,4]; in addition, because of the widely existed volcano eruption, thermal shock, exploitation of geothermal resources, and deep geological repositories for heat-generating radioactive wastes, which have significant effects on degradation degree of rocks, should be considered as an important factor in the H-C cycle process.In addition, the rock is not only in the environment of hydrothermal coupling in the underground and ground engineering but also subjected to the effects of engineering disturbance, earthquake, and other impact loads in the whole construction and operation process [5].Recurrent H-C cycles lead to changes in the microstructure of rocks, which cause generation of new microcracks and micropores inside rocks [3,[6][7][8].Rock tends to experience a deep weathering due to repeated erosion of water and temperature.Resistance ability of the rock to the H-C cycle depends on a complex set of material properties, including mineralogical composition, inner structure, and pore characteristics of rock [9].When the temperature floats down the freezing point of water, a 9% volume expansion occurs inside rock, which causes damage to rock caused by great ice pressure.e frequent freeze-thaw (F-T) of water expands the original fissures inside rock.e freeze-thaw phenomenon also induces new microcracks in rock, which causes great damage to rock engineering [10][11][12].Due to various coefficients for different mineral particles, temperature increment during the H-C cycle leads to secondary thermal stress in rock.Studies show that the main mechanism of H-C cycles on rock is the loss of water and

Specimen Preparation and Experimental Setup
2.1.Sample Preparation.Fine-grained sandstone samples were collected from the Zhujidong coalmine in Huainan city, Anhui province, China.Power X-ray diffraction analysis was adopted to study the mineral composition of rock at 20 °C.According to the XRD spectrum (Figure 1), the main components of sandstone are quartz accounted for 62.8%, besides of which the sandstone was composed of 12% albite, 6.98% kaolinite, and other minerals [22,23].As seen from Figures 2(a) and 2(b) and Table 1, the main elements of sandstone are Si, O, Ca, Fe, and Al.Scanning electron microscope (SEM) was performed on some offcuts of rock specimens, and chemical analyses were also conducted by energy dispersive spectroscopy (EDS).e primary fractures and pores were found by SEM.As seen from the EDS picture and Table 1, the brighter the color block, the more elements it contains, such as the relatively large numbers of elements in block 3. ese primary fissures, pores, and elements have a significant effect on the physical and chemical properties of rocks.e SEM and EDS images of sandstone after different H-C cycles groups are shown in Figures 2(c)-2(h).It can be seen from the diagram that the cracks and pores in rocks increase with the increase of high temperature treatment temperature, showing a phenomenon of high temperature (H-T) group > moderate temperature (M-T) group > low temperature (L-T) group.Additionally, the above figures also illustrate that the damage degree of rock after 12 H-T cycles is obviously higher compared with that after M-T and L-T treatment.Samples with similar P-wave velocity were selected and processed into φ50 mm × 25 mm cylinder with surface parallelism within 0.05 mm and surface flatness within 0.02 mm.Sample processing dimensions and methods comply with the ISRM (International Society for Rock Mechanics) standard [24].Surface planeness is controlled to ±0.05 mm, and the vertical deviation of the upper and lower surfaces is ±0.25 °.   e whole H-C cycles process were shown as follows: sandstone samples were first immersed in water for 12 hours, then put them in a low-temperature test chamber at −20 °C for 6 hours, and finally, heat samples in a high-temperature box at a predetermined temperature for 6 hours.A whole H-C cycle lasts for 24 h.e mass, volume, and P-wave velocity of sandstone samples were tested before and after different H-C cycles.

H-C Cycling and
SHPB equipment has been successfully used to investigate the dynamic behaviors of materials, such as rock and concrete [25].In this test, dynamic impact loading tests were conducted by the SHPB system with a diameter of 50 mm.As shown in Figure 4, it contains launch device, striker, incident bar, transmitted bar, buffer bar, and strain acquisition instrument, and two strain gauges are glued on the incident and transmitted bars to collect the origin strain signals; in addition, a striker with double-tapered shape was adopted to guarantee the stress balance of the rock specimen [26].e SHPB test steps are as follows: (1) check all the test parts on the working state and adjust voltage balance, (2) appropriate lubricant is evenly smeared on the two side surfaces of the rock specimen and then put it between incident and transmitted bars, (3) open the gas pressure switch and collect the test signal, and (4) collect the fragment of the rock specimen.e dynamic stress-strain curves could be calculated according to the "three wave methods" [27][28][29]: where ε I (t), ε R (t), and ε T (t) are the strain of incident wave, reflected wave, and transmitted wave, respectively; E 0 , A 0 , C 0 , A s , and l s are Young's modulus, the cross-sectional area, the elastic wave speed of the bar, and the cross-sectional area, and the height of the rock specimen, respectively; and t is the duration time of the elastic wave.C 0 can be calculated by the following equation: where ρ 0 is the density of the bar material.

Test Results and Analysis
3.1.Variation in P-Wave Velocity.e P-wave velocity results of sandstone measured before and after different H-C cycles are shown in Figure 5.It can be noticed that before H-C cycles, the average wave velocities of the specimens range from 2870 m/s to 3067 m/s.For sandstone specimens after H-C cycles, it shows a decrease tendency of P-wave velocity.After 20 and 4 H-C cycles, the descent rate is relatively slow for L-T and M-T groups and the H-T group, respectively.e descent rate of P-wave velocity from 0 to 20 cycles is 7.93% and 11.52% for the L-Tgroup and M-Tgroup, respectively, while it is 41.09% for the H-T group from 0 to 4 cycles.As a contrast, after 20 cycles, the descent rates are only 1.91% and 4.42% for L-T and M-T groups, respectively, while the descent rate is only 19.42% for the H-T group from 4 to 12 cycles.Hence, 20 cycles is the critical point of P-wave velocity for the L-T group and M-T group, and 4 cycles is the critical point for the H-T group.ere are two reasons for the variation of the P-wave velocity of sandstone.On the one hand, the water in the rock firstly absorbs water to attain saturation and then escapes because of the increase of temperature, which leads to the increase of pore volume due to the repeated absorption and escape of this kind of water.On the other hand, the generation and expansion of internal microcracks and micropores also lead to the attenuation of ultrasonic wave energy after rock damage.

Variation in Density and Porosity.
Rock is naturally a three-phase system of air-water-solid, and rock density refers to the mass of rock in unit volume, which is significance for studying weathering process, stability of rock mass, and prediction of surrounding rock pressure in rock engineering [30].By testing the volume and weight of rock Advances in Civil Engineering 5 samples before and after different H-C cycles, we obtained the variations in density with cycling number, as shown in Figure 6.It can be seen from the diagram that the density of sandstone samples gradually decreases with increasing cycling number for three groups.Figure 6 also shows that the relationship between density and cycling number is similar to that of the P-wave velocity.
e descent rates of density from 0 to 20 cycles are 1.14% and 1.80% for the L-T group and M-T group, respectively, while it is 2.44% for the H-T group from 0 to 4 cycles.After 20 cycles, the descent rates are only 0.27% and 0.33% for the L-T group and M-T group, respectively, while it has a 0.62% descent rate for the H-T group from 4 to 12 cycles.Hence, 20 cycles is the critical point of density for the L-T group and M-T group, and 4 cycles is the critical point for the H-T group.
Rock is a natural polycrystalline material with many defects and pores [31].Porosity has great effect on the hydration and thermal and mechanical properties of rock and rock mass.e existence of pores makes it more vulnerable subject to external load, leading to further deterioration of rock.erefore, the porosity change of rock 6 Advances in Civil Engineering after di erent H-C cycles is measured and analyzed in this test.
Figure 7 illuminates the variation of porosity of sandstone versus H-C cycles.It is clearly notice that (1) under the same cycling number, the porosity of sandstone for the H-T group is much higher than that for L-T and M-T groups.(2) A di erent changing trend is found for three groups, and for L-T and M-T groups, with the increase of cycling number, the porosity exhibits a linear correlation.However, there shows a exponential growth relationship between porosity and cycling number for the H-T group, which attributes to the high temperature damage to the sandstone specimen.(3) e rise rate of porosity is higher from 0 to 20 cycles compared with that from 20 to 40 cycles for L-T and M-T groups, and the values of the rise rate of porosity from 0 to 20 cycles are 39.64% and 92.93% for L-T and M-T groups, respectively, while they are only 10.17% and 18.87% from 20 to 40 cycles.By contrast, the average value of porosity for the H-T group increases from 2.23% to 5.29% with the cycling number increasing from 0 to 4 cycles, with a 137.22% rise rate, while it changes from 5.29% to 7.34% when the cycling number increases from 4 to 12 cycles, with only 38.76% rise rate, which is much lower compared with that from 0 to 4 cycles.

Dynamic Stress-Strain Curve and Relative Elastic
Modulus.Typical dynamic stress-strain curves are presented in Figure 8.By comparison of complete stress-strain curves of di erent groups, we can see that with the increase of axial strain, the stress increase rate gradually decreases with increasing cycling number before peak strain.e peak stress decreases with increasing cycling number.Furthermore, the brittleness of the sandstone samples gradually decreases and the ductility gradually increases as the increase of the cycling number.
e dynamic compressive strength is de ned as the peak stress in the stress-strain curve in this research, and test results of dynamic compressive strength of the sandstone specimen at the strain rate of (190 ± 10) s −1 under di erent H-C cycles are shown in Figure 9(a).
Figure 9(a) results show that the dynamic compressive strength of rock decreases slowly after 40 cycles for L-T and M-T groups, while it decreases sharply after 4 cycles for H-T cycles.Speci cally, the dynamic compressive strength of the control group is 143.17 MPa, and it only decreases to 103.80 MPa and 96.30MPa after 40 cycles for L-T and M-T groups, respectively, and the corresponding descent rates are 27.50% and 32.74%, respectively.However, after 4 cycles for the H-T group, the dynamic compressive strength decreases to 64.38 MPa, with 55.03% descent rate, which is much larger than the other two groups.

Advances in Civil Engineering
Elastic modulus in this research is de ned as the slope corresponding to 40% and 60% of compression strength during the rising phase of the stress-strain curve; the calculated results are shown in Table 2. e results showed that the elastic modulus decreased with the increasing cycling number, in order to describe the drop degree of elastic modulus caused by the H-C cycle, the relative elastic modulus (K r ) is de ned and can be calculated as follows:   Advances in Civil Engineering where K r is the relative elastic modulus and E n and E 0 are dynamic elastic modulus of the rock specimen after 0 and n H-C cycle numbers, respectively.Table 2 and Figure 9(b) illustrate the variation in elastic modulus and K r of sandstone with different cycling numbers.For the L-T group, K r decreases to 0.688 after 20 cycles, and then, it has a slight reduction from 20 to 40 cycles.In addition, K r for the M-T group is larger than that for the L-T group, and after 30 H-C cycles, the values of K r are 0.667 and 0.538 for L-T and M-T groups, respectively.However, the elastic modulus of the specimen for the H-C group has a larger descent rate after 8 cycles, and the corresponding value of K r decreases from 1 to 0.297.e decrease trend of the H-T group is obviously larger than that of L-T and M-T groups.

Variation in Damage.
Using damage mechanics methods to study questions of rock hydrothermal coupling is a new development in rock mechanics.
e damage variable can describe the deformation and failure degree of the sandstone specimen after different H-C cycles.Two defined methods are adopted and compared in this research to describe the damage degree of sandstone after H-C cycles.Various changes occur in the internal structural of rocks after the H-C cycle, which affect the refraction and diffraction phenomena of ultrasonic propagation.Hence, there has a relationship between wave velocity and damage degree of rock after the H-C cycle.Based on this, damage variable is defined by the variation of longitudinal wave velocity [32]: where D p stands for damage variable calculated by variation of wave velocity, V n is the P-wave velocity after n H-T cycle numbers, and V 0 is the P-wave velocity of rock without H-C cycles.e other effective calculation method to describe the damage variable is based on the variation of dynamic strength [18], which can be calculated according to equation (3): where D σ stands for the damage variable calculated by variation of dynamic strength, σ n is the dynamic compressive strength after H-T cycle, and σ 0 is the dynamic compressive strength of the specimen without the H-C cycle.
In addition, the release rate of damage strain energy can describe the damage degree inside rock.e relation between damage and release rate of damage strain energy can be expressed as [33] where θ is the release rate of damage strain energy, σ is the dynamic compressive strength, E is the elastic modulus, and D is the damage variable of rock caused by H-C cycles.e calculated results of damage variable and release rate of damage strain energy of sandstone are listed in Table 3 and Figures 10 and 11.
As can be observed from Table 3 and Figures 10 and 11, both the damage variable and release rate of damage strain energy can describe the damage variation tendency of the sandstone specimen after different H-C cycles.However, compared with the release rate of damage strain energy, damage variable can better reflect the damage degree; the reason can be clarified as follows.With the increase of cycling number, the damage variables increased gradually, which was in accordance with the general rule of damage evolution in the L-T group.However, the release rate of damage strain energy shows a first-increase-and-then-decrease tendency, which is inappropriate according to the damage evolution rule.Additionally, from the test results of P-wave velocity, porosity, density, and dynamic strength, it is clearly noticed that 20 cycles are critical point of density for L-Tand M-Tgroups, and 4 cycles is critical point for the H-T group, which is similar to the relationship between damage variable and cycling number, while the turning points are different for release rate of damage strain energy, as shown in Figure 11.Advances in Civil Engineering e sandstone used in the tests was a typical sedimentary rock with some original pores and cracks in Figure 2.After repeated action of water and temperature, various internal damage degrees happens inside rock, and the main reasons are as follows: (1) Repeated erosion of water is one of the causes of rock damage after H-C cycles.On the saturated process, there are original cracks and pores in the rock, as shown in Figure 1.In the process of water absorption, many minerals are hydrophilic and soluble, which promote the weathering of rocks [13,[34][35][36].e frequent freezing and thawing of pore water inside rock expands the cracks and pores and promotes the development of new microfractures [7,11,37].
(2) Rock specimens are taken out after heat treatment at 50 °C, 100 °C, and 400 °C.e di erent states of water (attached water, bound water, and mineral combined water) existing inside the rock will evaporate after heat treatment.Particularly, the weakly bound water can completely escape at about 150 °C, and the strong bound water can completely escape only at the temperature of about 200 °C to 300 °C.e crystalline water evaporates and escapes at least 400 °C.e loss of component water and crystalline water will lead to the destruction of the crystal structure of sandstone mineral [14,21].(3) Because of various thermal expansion coe cients of mineral composition in rock, induced thermal stress between crystal and ssure leads to aw and exfoliation of rock after heat treatment [32,33].After 400 °C, the ssures between the crystalline grains in the interior of rocks shows a cracking phenomenon, and the crack cannot be recovered [38,39].e mass fractal dimension of fragmentation (D f ) under di erent test conditions is calculated.e fragment distribution parameter can be calculated by the following formulas [40]: where m r is the accumulative mass under sieve of characteristic size r, M is the total mass of the specimen, r m is the particle size, and b is the fragment distribution parameter.e b value is obtained from mass fractal dimension of fragmentation by Formula (10).As shown in Figure 13, with the increase of cycling number, the value of b decreases and the degree of fragmentation increases.
e mass fractal dimension of fragmentation is calculated by the following equation: where D f is the mass fractal dimension of fragmentation.e relationship between D f and cycling number is shown in Figure 14.With the increase of cycling number, the destruction degree of sandstone becomes greater, and the corresponding fragments show more evenly.As seen from Figure 14, the D f of sandstone ranges from 2.024 to 2.304 under di erent test conditions in this research.e values of D f of sandstone increase with increasing cycling number; furthermore, there is a linear relation between mass fractal dimension and cycling number for the L-T group and M-T group, while it exhibits an exponential relation for the H-T Advances in Civil Engineering 11 group.e larger the mass fractal dimension of fragmentation is, the ner the fragmentation is.From the above analysis, it can be noticed that the cycling process promotes the degradation inside the rock specimen, and there exhibits a positive correlation between damage degree and D f , agreeing with the previous studies results [41,42].

Conclusions
(1) With the increase of cycling number, the P-wave velocity and density decrease, while the porosity increases.Additionally, 20 cycles is critical point of P-wave velocity for the L-T group and M-T group, and 4 cycles is critical point for the H-T group.Both the dynamic compressive strength and K r decrease with increasing cycling number for three groups, and the descent rate of them is relatively lower for L-T and M-T groups compared with that for the H-T group under the same cycling number.(2) e damage degree of rock increases with the increase of cycling number.Compared with the release rate of damage strain energy, the damage variable can better describe the damage degree of sandstone after the H-C cycle.(3) e D f of sandstone increase with increasing cycling number; furthermore, there has a linear relation between mass fractal dimension and cycling number for the L-T group and M-T group, while it exhibits an exponential relation for the H-T group.e larger the mass fractal dimension of fragmentation is, the ner the fragmentation is.
Data Availability e datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.Advances in Civil Engineering

Figure 6 :Figure 7 :
Figure 6: Variation of density of sandstone versus H-C cycles.

Figure 8 :
Figure 8: Dynamic stress-strain curve of sandstone specimens after H-C cycles: (a) L-T group, (b) M-T group, and (c) H-T group.

Figure 9 :
Figure 9: E ects of H-C cycles on peak stress and relative elastic modulus under di erent groups: (a) Relation of dynamic compressive strength with cycling number and (b) relation of relative elastic modulus with cycling number.

Figure 10 :Figure 11 :
Figure 10: Relation of damage variable with cycling number.

Figure 14 :
Figure 14: Relation of mass fractal dimension of fragmentation with cycling number.(a) L-T group and M-T group and (b) H-T group.

Table 2 :
Test results on the basic dynamic properties of sandstone after H-C cycles.
55. Fracture of Sandstone.In order to further illuminate the fracture characteristic of sandstone after the H-C cycle, the fracture result is quanti ed by mass fractal dimension of fragmentation.e rock fragments are collected and screened by di erent sizes after the SHPB test.e sizes are 0.15 mm, 0.3 mm, 0.6 mm, 1.18 mm, 2.36 mm, 4.75 mm, 9.5 mm,

Table 3 :
Test results on damage of sandstone after H-C cycles.