
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38862575
62886
10.1038/s41598-024-62886-3
Article
Deterioration mechanisms of tuff with surface fractures under freeze–thaw cycles
Lai Runsen 1
Zhang Zizhao zhangzizhao@xju.edu.cn

12
Zhu Jianhua 3
Xu Zhiguo 4
Wei Xiancheng 3
Liu Xinyu 1
Xiong Binghui 1
1 https://ror.org/059gw8r13 grid.413254.5 0000 0000 9544 7024 College of Geology and Mining Engineering, Xinjiang University, Urumqi, 830017 China
2 https://ror.org/059gw8r13 grid.413254.5 0000 0000 9544 7024 Research Base of Xinjiang University, National Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep-Earth Engineering, Xinjiang University, Urumqi, 830017 Xinjiang China
3 Xinjiang Geological Environment Monitoring Institute, Urumqi, 830099 Xinjiang China
4 Xinjiang Uygur Autonomous Region Geological and Mineral Exploration and Development Bureau, Fourth Geological Brigade, Altay, 836500 Xinjiang China
11 6 2024
11 6 2024
2024
14 1340217 1 2024
22 5 2024
© The Author(s) 2024, corrected publication 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
In practical engineering, the development of surface cracks is one of the most important reasons for the destruction of the rock mass, and the development of complex morphology fractures on the rock mass surface significantly influences rock mass mechanics. This paper addresses the freeze–thaw damage issue in rock mass containing surface fractures in cold regions. Tuffs and a control group are selected as research samples, with the control group being prefabricated surface jointed specimens with an inclination angle of 70° and a fracture depth ratio d (crack depth) /t (sample width) of 0.26. The study analyzes the mass, wave velocity loss, macro-microcosmic fracture damage morphology, and mechanical properties of the two specimen groups through laboratory freeze–thaw cycle tests, uniaxial compression tests, and scanning electron microscopy examinations. The results show an overall decrease in mass, wave velocity, and uniaxial compressive strength as the cycle number increases, with the prefabricated jointed group samples showing more significant changes. However, the two specimen groups exhibit different macroscopic failure fracture states. In addition, scanning electron microscopy images illustrate that after freeze–thaw cycles, the large rock mass particles break into smaller fragments, resulting in looser particle arrangements and a transition from initial surface cementation to point contact., which weakens the compressive strength of the rock mass. The paper also explains the mechanism of the diminishing impact of freeze–thaw cycles on the strength of the rock mass after a certain number of cycles. The research outcomes hold significant reference value for engineering construction in cold regions.

Keywords

Freeze–thaw cycles
Surface fractured rock mass
Prefabricated joint
Mechanical property
Fractures morphology
Subject terms

Natural hazards
Engineering
Formation mechanism and early-warning criteria of loess landslide with ice and snow melt water superimposed with rainfall infiltration in Ili Valley42367021 42367021 42367021 42367021 Lai Runsen Zhang Zizhao Liu Xinyu Xiong Binghui issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The damage of rock mass is often associated with joints, fractures, or other structural surfaces. Natural rock mass inherently contains numerous joints, and these masses experience cyclic loading from frost heaving force under the influence of freeze-thaw cycles. This action causes subsequent expansion and closure of joints and fractures, resulting in the coherent structural surfaces, which will lead to the destruction and destabilization of the rock mass. Studies on rock mass destruction and destabilization primarily show that rock mass damage under freeze-thaw cycles is mainly attributed to cyclic frost heaving force. Wettlaufer et al.1 indicate an interface effect within the rock mass. During the freezing process, a portion of unfrozen water located at joints, pores, or fractures migrates to the bottom due to the frozen water, thereby continually increasing water pressure and leading to downward fracture expansion. Setzer2 discovers that the water-ice phase transition induces a volume expansion of about 9%, with some water in the pores and fractures remaining liquid after freezing. Liu et al.3 propose that frost heaving force results from the common phase transition of in-situ water in rock mass fractures and migrating water, causing fractures to expand and ultimately leading to rock mass damage. The above studies illustrate the failure mechanism of rock mass under freeze-thaw cycles, a critical consideration in rock mass engineering where safety and durability requirements are increasingly stringent. Rock mass and soil experience year-round low temperatures and freeze-thaw cycles in alpine regions; thus, investigating mechanical properties under extreme environmental conditions becomes highly significant.

In practical engineering, rock masses exhibit a higher occurrence of joints, fractures, and various structural surfaces following the construction of tunnels, rocky slopes, etc. The intricate nature of jointed and fractured rock masses significantly influences the safety and stability of rock mass engineering. The focus of research on fractured rock masses primarily includes three categories: natural joints and fractures within the protolith, prefabricated joints and fractures within the protolith, and prefabricated joints and fissures within similar materials. Rodríguez et al.4 conduct loading and unloading splitting tests on marble with microfractures, revealing fracture extension evolution in different conditions through the acoustic emission (AE) method. Bao et al.5 uncover the anisotropic freeze-thaw damage law in rock masses by producing protolithic jointed and fractured rock masses with different bedding angles and conducting freeze-thaw and uniaxial tests. Pan et al.6 investigate the uniaxial compressive mechanical properties of marble with prefabricated fractures containing various fillings. Li et al.7 explore the acoustic emission law and seismic source mechanism of granite specimens with prefabricated fractures during hydraulic fracturing, noting shock source concentration at the tips of prefabricated cracks. Wang et al.8 conduct true triaxial loading and unloading tests on similar materials of rock mass specimens with prefabricated joints and fractures, identifying deformation and damage characteristics. Zhu9,10 determines the change law of the stress intensity factor under compressive loading tests on similar materials of rock mass specimens with prefabricated joints and fractures, which include single and double cracks.

Synthesizing the results of the above research, it is evident that cyclic freeze-thaw force induces gradual deterioration and destabilization of rock masses. The initiation of rock mass body damage often emerges from surface fractures, posing challenges to engineering construction in cold regions. It becomes particularly urgent to delve into the laws governing the degradation and damage of rock masses containing surface fractures under the influence of freeze-thaw cycles. Previous studies show that the test rock mass typically features penetrating prefabricated fractures, assuming fractures that penetrate the entire rock mass. Nonetheless, actual engineering scenarios commonly involve non-penetrating fractures (surface fractures)11, where fracture expansion occurs within the rock mass in three-dimensional space. The suppose of penetrating prefabricated fractures neglects the depth dimension of the fracture, potentially missing crucial three-dimensional information about the rupture12. R.H.C. Wong et al.13,14 establish a close relationship between the expansion mode of surface fractures and the fracture depth ratio (d/t) through the examination of rock-like materials. Nakamura15 measures the surface temperature of tuff fractures, establishing the relationship between frost heaving force and fracture temperature. Huang et al.16 demonstrate that the penetration of three-dimensional surface cracks is more intricate than two-dimensional penetration cracks using PMMA transparent materials. Hence, there is a need to shift the research focus to surface fractures to gain comprehensive insights into the laws governing rock mass fracture extension.

Presently, the researchs about the properties of rock mass with prefabricated joints under freeze-thaw action predominantly center on the mechanical strength of the rock mass16–21,23,24, followed by alterations in the microstructure22–27. Other investigations encompass the angular influence of prefabricated joints17,18, penetration degree17, and the comparison of different lithologies20 etc. The associated research methods primarily revolve around laboratory tests, including uniaxial compression tests16–19, triaxial tests20,21,23, acoustic emission tests19, polarized light microscopy22, scanning electron microscopy23,24, nuclear magnetic resonance (NMR), CT scan, etc.25, and CT scanning26,27, etc. This paper will employ freeze-thaw cycle tests, uniaxial compression tests, and SEM tests to investigate prefabricated surface jointed rock mass samples.

Limited studies have delved into the cracking process and expansion laws of surface fractured rock masses, particularly the correlation between macroscopic mechanical strength characteristics and deformation characteristics of surface fractured rock masses under the extreme conditions of freeze-thaw cycles. The actual projects often involve more joints and fractures on the rock mass surface, significantly influencing the rock's deterioration due to freeze-thaw cycles. This paper fabricates surface jointed specimens and subjects them to uniaxial compression tests after different freeze-thaw cycles. Along with research on macroscopic crack extension and microscopic analysis using scanning electron microscopy, this study aims to unravel the freeze-thaw damage laws of tuff and investigate the mechanical properties and crack extension of fractured rock mass under freeze-thaw cycle conditions.

Laboratory tests

This research studies the rocks extracted from the Camel's Hump on the western side of Altay City's downtown area in Xinjiang. The region experiences an average annual temperature of 4.80 °C, with the highest temperatures occurring in July (averaging around 20 °C) and the lowest temperatures in January (with an average of – 30 ℃ and an extreme minimum of − 41.7 °C). Winter spans a considerable duration of 127–169 days, creating a significant temperature differential between seasons. Given these conditions, the rocks in this area are highly influenced by the freeze-thaw action.

Specimen preparation

The test tuff is characterized mainly by tuff lithology and obtained from the collapsed slope of Camel's Hump Scenic Area in Altay City, Xinjiang. Specimens are prepared into Ф50 × 100 mm rock samples. Utilizing a rock wave velocity detector, samples are screened to ensure an average wave velocity of 4.448 m/s and a rock density of 2786 kg/m3, aligning with the grouping criteria for the tests.Based on the actual survey, the time when the rock is most strongly affected by freezing and thawing is mainly concentrated in the alternation of winter and spring, and the freezing and thawing period is about 1 month.

To prevent damage to the rock mass structure around the joints during the waterjet machining process, an automated machine tool is employed to create non-through joints. These joints, with a thickness of approximately 1.5 mm, length of 20 mm, and a depth (d) of 13 mm (see Fig. 1), are statistically determined based on the most prevalent structural plane at an inclination angle of 37°∠70°, identified through joint statistics of the collapsed hazardous rock mass28. The chosen inclination angle for the prefabricated joints is 70°. Studies indicate that when the crack depth ratio (d/t) exceeds 1/313,14,29, the fracture penetrates the specimen in the thickness direction, creating a stress state close to the two-dimensional state. To simulate the expansion process of surface fractures in a three-dimensional stress state, this study chooses a crack depth ratio of d/t=0.26. This ensures that fracture expansion occurs entirely within the specimen, maintaining a three-dimensional stress state. Following specimen processing, those of superior mass are selected for testing.Figure 1 Schematic diagram of tuff specimen dimensions.

Test scheme

Freeze-thaw cycle test

The test specimens are categorized into two groups: original rock samples (Group A) and prefabricated jointed rock samples (Group B) (Fig. 2a). Each group is further subdivided into seven sets, each comprising three specimens, resulting in a total of 42 standard specimens for testing. Freeze-thaw cycles are conducted 0, 1, 3, 5, 10, 20, and 30 times, with the freeze-thaw temperature set to – 30 °C~20 °C to emulate the climatic conditions of the site. Both Group A and Group B specimens are initially placed in an oven at 105 °C for 24 h (Fig. 2b). Following the recording of dry specimen weights, they are subjected to a freeze-thaw cycler (Fig. 2c), maintaining a temperature of (20 °C ± 1 °C) for 8 hours, followed by stabilization at (−30 °C ± 1 °C) for 8 hours, constituting a 16-hour cycling cycle (Fig. 2d). At the conclusion of each freeze-thaw cycle set, the specimens are retrieved and placed in a drying box. Various parameters such as mass and wave speed are measured, and macroscopic observations of freeze-thaw damage to the specimen's appearance are also recorded.Figure 2 Freeze–thaw cycle test process.

Uniaxial compression test

After each freeze-thaw cycle set and 24-hour drying, three corresponding specimens are subjected to uniaxial compression tests. Using the MTS E45.605 mining rock mechanical property test system (Fig. 3a), the rock samples undergo uniaxial compression tests at a loading rate of 0.1 MPa/s until specimen destruction. This analysis aims to explore the correlation between macroscopic mechanical properties, surface fractures, and the number of freeze-thaw cycles. Abnormal data are excluded, and the remaining results are averaged for each group.Figure 3 Rock specimen test system. (a) Mining rock mechanical property test system. (b) Scanning electron microscopy.

Scanning electron microscopy test

Specimens that completed the freeze-thaw cycle test and uniaxial compression test undergo testing and analysis through the scanning electron microscopy (SEM) method to observe damage in typical rock samples (see Fig. 3b).

Tests results and analysis

Based on the outlined test program, both the original rock samples in Group A and the prefabricated jointed rock samples in Group B undergo freeze-thaw cycles and subsequent mechanical tests. The corresponding results and analyses are detailed below.

Changes in specimen mass after freeze–thaw cycles

Post freeze-thaw cycles, alterations in specimen mass, which is a crucial indicator of freeze-thaw damage, are observed. The initial mass of the dried specimen (m0) is weighed after grouping, and the dry mass before and after each cycle is compared to derive the rate of mass change (Rx), characterizing the specimen's mass variation.

The equation for the rate of change in mass for different cycle numbers is given by:1 Rmn=m0-mnm0×100%

where, Rmn represents the rate of change in mass for different cycle numbers, m0 represents the initial rock sample mass in g and mn represents the rock sample mass in g corresponding to n times cycle.

As indicated in Table 1, the mass of both specimen groups gradually decreases with an increase in the number of cycles, albeit with a relatively small overall change. Calculate the freeze-thaw mass rate of change and plot the graph (Fig. 4). The results demonstrate a gradual increase in the mass change rate along with the increasing numbers of freeze-thaw cycles for the two sample groups. The magnitude of mass loss exhibits a trend of initial increase followed by a subsequent decrease.Table 1 Average mass change under freeze–thaw cycles.

Group A of original rock samples	Group B of prefabricated jointed rock samples	
Cycling numbers	Pre-freeze–thaw cycles (g)	Post-freeze–thaw cycles (g)	Mass difference (g)	Mass change rate (Rx)	Cycling numbers	Pre-freeze–thaw cycles (g)	Post-freeze–thaw cycles (g)	Mass difference (g)	Mass change rate (Rx)	
1	580.35	580.23	0.13	0.022%	1	528.32	528.15	0.17	0.033%	
3	563.85	563.69	0.16	0.029%	3	573.94	573.72	0.22	0.039%	
5	555.27	555.09	0.17	0.031%	5	508.30	508.04	0.26	0.052%	
10	513.97	513.73	0.25	0.048%	10	535.96	535.59	0.38	0.070%	
20	573.90	573.61	0.29	0.051%	20	552.36	551.95	0.41	0.075%	
30	539.48	539.16	0.31	0.058%	30	536.73	536.30	0.44	0.081%	

Figure 4 Changes in mass of group A and B specimens after freeze–thaw cycles.

Notably, at the 10th cycle, with a faster increase before and a slower increase afterwards. Post 10 cycles, mass loss exhibits a gradual increase, and mass change rate levels off in the later cycles. This suggests that initial test cycles have a less apparent effect on the dense specimen, with subsequent cycles inducing fatigue damage, particle shedding, and pore fracture expansion. At the 10th cycle, evident mass decline occurs, marking the start of surface particle shedding. Given that the tuff specimen is relatively dense, the freeze-thaw effects on the specimen surface are limited. Therefore, the mass change rate tends to stabilize after the 10th cycle.

Comparative analysis between Groups A and B shows a significantly smaller mass change rate in Group A. Moreover, at the 5th cycle, Group B exhibits an apparent sudden increase in mass change rate, indicating that prefabricated surface joints have caused damage to the rock mass structure. The freeze-thaw impact became more pronounced, causing not only surface particle detachment but also shedding within and around the prefabricated surface joint fractures.

Changes in specimen wave velocity after freeze–thaw cycles

As per the previous discussion, the equation defining the change rate of wave velocity for different cycle numbers is expressed as follows:2 Vvn=P0-PnP0×100\%

where, Vvx represents the change rate of wave velocity for different cycle numbers, P0 represents the initial wave velocity of the rock sample in km/s, and Pn represents the wave velocity of the rock sample in km/s corresponding to n times cycle.(The wave speed here is the longitudinal wave speed.)

Consistent with the law of mass change, the wave velocity of both sample groups gradually decreases (refer to Table 2). As illustrated in Fig. 5, with an increase in the number of cycles, the rate of change in wave velocity for both groups exhibits a gradual increase, rate of change of wave velocity increases levelling off. As the freeze-thaw cycles accumulate, the change rate in wave velocity tends to gently ascend. During freeze-thaw cycles, the rock's internal pores are compressed by frost heaving forces, leading to pore expansion, a decrease in internal densification, and consequently, a gradual reduction in wave velocity. The change rate in wave velocity increases progressively during this process.Table 2 Average wave velocity change under freeze–thaw cycles.

Group A of original rock samples	Group B of prefabricated jointed rock samples	
Cycling numbers	Pre-freeze–thaw cycles (km/s)	Post-freeze–thaw cycles (km/s)	Wave velocity difference (km/s)	Change rate of wave velocity (Vx)	Cycling numbers	Pre-freeze–thaw cycles (km/s)	Post-freeze–thaw cycles (km/s)	Wave velocity difference (km/s)	Change rate of wave velocity (Vx)	
1	4.688	4.626	0.06	1.33%	1	3.704	3.487	0.22	5.85%	
3	4.479	4.363	0.12	2.59%	3	4.027	3.755	0.27	6.77%	
5	4.226	3.981	0.24	5.78%	5	4.288	3.949	0.34	7.91%	
10	4.167	3.875	0.29	7.00%	10	3.776	3.432	0.34	9.10%	
20	4.801	4.405	0.40	8.23%	20	4.511	4.083	0.43	9.49%	
30	4.029	3.641	0.39	9.63%	30	4.085	3.683	0.40	9.85%	

Figure 5 Changes in wave velocity of group A and B specimens after freeze–thaw cycles.

Observing the wave velocity changes in groups A and B, it is evident that the change rate in wave velocity in group A is smaller than in group B. However, the change in wave velocity in group B changes more smoothly. In group A, at the 5th freeze-thaw cycle, which increased rapidly before and then slowly thereafter. This surge is attributed to the initial balance in pore expansion in the original rock samples during the first three cycles, and when the cycles reach five, this balance is disrupted, leading to substantial pore expansion and a consequent decrease in rock mass densification. In contrast, the prefabricated surface joints in group B have disturbed the balanced pore state within the rock mass body. Consequently, the pore spaces in group B are more susceptible to expansion and extension, resulting in a smoother change in wave velocity compared to the more abrupt changes observed in group A.

Apparent damage characteristics after freeze–thaw cycles

According to the freeze-thaw cycle test process, both Groups A and B exhibit noticeable surface changes after five cycles, with no chunk falling out. Therefore, typical samples with freeze-thaw cycles of 5, 10, 20 and 30 times were selected for comparative analysis and research. Figure 6 is a sample without freeze-thaw cycle.Figure 6 Apparent characteristics of Group A and B specimens not subjected to freeze–thaw cycles.

As shown in Fig. 7a, with the increase of freeze-thaw cycle number, the joint cracks on the surface of Group A specimens increase, gradually intertwining and forming damage on the surface. When the cycle number reaches five times, the weak layers in the direction of bedding on the specimen surface begin to emerge, becoming rougher compared to the original specimen. Subsequently, at ten cycles, evident joint fractures emerge, and by the twentieth cycle, these fractures increase, intensifying the freeze-thaw-induced damage. At thirty cycles, these fractures intersect and develop into approximately 0.4 mm wide cracks, indicative of advanced damage. Since Group B consists of prefabricated jointed rock specimens, the freeze-thaw effects are more pronounced than in Group A. In Fig. 7b, minimal debris detachment and small cracks are visible at the edge of the specimens after five cycles. By ten cycles, several groups of joint cracks form, deepening by the twentieth cycle. At thirty cycles, fractures develop at both ends of the specimen along the direction of the prefabricated joints. However, due to the dense nature of tough specimens, the specimen does not show overall damage.Figure 7 Apparent damage characteristics of group A and B after freeze–thaw cycles.

Analysis of mechanical properties of rock mass

Moving to the characterization of uniaxial compressive strength, tests using the MTS E45.605 mining rock mass mechanical property test system are conducted on both Groups A and B for 0, 1, 3, 5, 10, 20, and 30 freeze-thaw cycles. The loading rate was 0.1 MPa/s until specimen destruction, with the test system automatically collecting and storing data. The stress-strain curves obtained (Fig. 8) reveal four stages in the uniaxial compression process: compaction, elastic deformation, plastic deformation, and destruction.Figure 8 Stress–strain curves under different numbers of freeze–thaw cycles. (a) Group A. (b) Group B.

The peak strength of uniaxial compression decreases with increasing freeze-thaw cycles the expansion of microfractures and the increase of soft structures within the rock mass. The compaction phase on the stress-strain curves of both sets of specimens is extended after being subjected to freeze-thaw cycles, indicating that the freeze-thaw cycles create more micropores. The elastic phase is approximately a straight line, and there are more fluctuations during the plastic deformation phase, which is attributed to the more scattered distribution of joints and pores within the rock mass. Curves exhibit vertical trends in the destructive phase, signifying brittleness after freeze-thaw cycles.

Comparing stress-strain curves for Groups A and B, it became evident that prefabricated jointed specimens in Group B are notably affected by freeze-thaw action. Substantial decreases in uniaxial compressive strength are observed in Group A after ten freeze-thaw cycles, while Group B experience this decline after only five cycles. This suggests that prefabricated joint specimens in Group B generate more microporosity in the rock mass after reaching initiation stress. Pores within rocks intersect and penetrate, resulting in damage and peak strength. As the number of freeze-thaw cycles increases, the stress is accompanied by a more minor strain.

Analysis of uniaxial compressive strength characteristics of rock mass

The rock mass rupture is the macroscopic manifestation of the cumulative deformation and destruction of microfabrication during the process of stressing, resulting from the sprouting and expansion of microcracks until their connectivity is achieved30. Due to the difficulty in accurately recording the fracture process of freeze-thawed rock mass cracks during uniaxial compression tests, representative damage morphologies under loading, as shown in Fig. 9, are selected for analysis.Figure 9 Rock mass fractures damage morphology under different numbers of freeze–thaw cycles.

Sketches of rock mass fractures are created to enhance the observation of the broken ring pattern of rock mass fractures (Fig. 10). With an increasing number of freeze-thaw cycles, the fractures on the specimen's surface under uniaxial compression gradually multiply. The fractures on the original specimen's surface are predominantly penetrating cracks, surrounded by numerous secondary fractures mostly oriented vertically in the loading direction. At 20 and 30 freeze-thaw cycles, secondary fractures increase, showing a different extension pattern compared to specimens without freeze-thaw cycles. The influence of freeze-thaw cycles leads to the appearance and aggregation of joints and fractures inside the rock mass, resulting breakage and the formation of multi-group wing fractures.Figure 10 Sketches of rock mass fractures damage morphology.

In contrast, energy release in prefabricated joints reduces the number of secondary fractures around the main fractures. However, the number of cracks increases similarly with the freeze-thaw cycle. Prefabricated jointed specimens show obvious shear damage at the prefabricated joint face under uniaxial compression, resulting in slip damage along the precast joint face and the formation of wing-shaped fractures at the joint's end. However, at 20 and 30 freeze-thaw cycles, fractures display anti-wing-shaped patterns, similar to those observed by Wong, R. H. C. et al.

The two-dimensional digital image processing of the rock mass crack image in Fig. 9 is carried out, and the processed image is transferred into the Matlab software plug-in Fraclab to calculate the fractal dimension. Figure 11 is a two-dimensional digital image processing process.Figure 11 Process diagram of 2D digital image processing of rock mass cracks.

The fractal dimension is used to verify the change rule of cracks in rock specimens with the increase of the number of freeze-thaw cycles, and the increase of fractal dimension represents the increase of the complexity within the image, which is calculated to obtain the change map of the fractal dimension of the rock specimens under different numbers of freeze-thaw cycles (Fig. 12). It can be seen that as the number of freeze-thaw cycles increases, the fractal dimension of the sample gradually increases, and the overall fractal dimension of the prefabricated joint sample group (group B) is larger than that of the original sample (group A), which is consistent with the above crack change law.Figure 12 Variation of fractal dimension with different number of freeze–thaw cycles.

Scanning electron microscopy test results

Scanning electron microscopy (SEM) analyses are employed on the examined specimens after freeze-thaw cycles and uniaxial compression strength tests. Specifically, SEM images at 200x magnification of samples from groups A and B, subjected to 0 and 30 freeze-thaw cycles, respectively, are scrutinized (Fig. 13).Figure 13 Scanning electron microscopy images of rock masses after 0-time and 30-time freeze–thaw cycles.

After freeze-thaw cycles, large rock mass particles broke into smaller ones, and increased pores are evident (Fig. 13c,d). The particles aggregate into a flocculent structure, and fractures become rougher, indicating looser particle arrangements. The connection between particles changes from original surface cementation contact to point contact under freeze-thaw cycles. Comparing original rock specimens (Group A) and prefabricated surface joints (Group B), it is evident that specimens in Group B mainta in surface cementation contact after freeze-thaw cycles, while Group A specimens mostly exhibit point contacts. The destruction caused by prefabricated surface joints leads specimens to follow existing joint fractures during unidirectional loading, maintaining a surface cementation contact mode.The face cemented contact allows the rock to better resist unidirectional loading, i.e. the face cemented contact is relatively stable and requires more energy to break it when subjected to unidirectional loading, thus the rock is more likely to be damaged along existing joints and fissures when subjected to unidirectional loading.

Discussion

This study explores the uniaxial compressive strength of both original rock masses and prefabricated jointed rock masses in relation to the number of freeze-thaw cycles to delve into the damage characteristics and mechanisms induced by freeze-thaw cycles on rock masses. Anomalies are excluded from each group's test results, and the remaining values are averaged. Figure 14 depicts the relationship between uniaxial compressive strength and the number of freeze-thaw cycles for both groups.Figure 14 Change of uniaxial compressive strength under different numbers of freeze–thaw cycles.

As the number of cycles increases, the rate of strength reduction gradually decreases, as illustrated in Fig. 14. The rock mass strength reduction rate gradually decreases with an increase in the number of cycles, as illustrated in Fig. 14. This trend is attributed to the initial freeze-thaw cycle, during which the pores in the rock mass have not yet been filled and undergo substantial expansion. As a result, the rock mass develops a stronger resistance to loading, leading to greater energy dissipation and faster reduction in strength during the early period of cycling. With an increase in the number of cycles, the addition and expansion of pore space in the rock mass result in a relatively fragile rock mass structure, making it susceptible to destruction by smaller loads during the initial stages, thus slowing down the rate of strength reduction.

The rock mass strength reduction rate gradually decreases with an increase in the number of cycles, as illustrated in Fig. 14. This trend is attributed to the initial freeze-thaw cycle, during which the pores in the rock mass have not yet been filled and undergo substantial expansion. As a result, the rock mass develops a stronger resistance to loading, leading to greater energy dissipation and faster reduction in strength during the early period of cycling. With an increase in the number of cycles, the addition and expansion of pore space in the rock mass result in a relatively fragile rock mass structure, making it susceptible to destruction by smaller loads during the initial stages, thus slowing down the rate of strength reduction.

The equations for the number of freezing and thawing versus uniaxial compressive strength were obtained by fitting the change curves of uniaxial compressive strength of two groups of rock samples under different freezing and thawing cycles (Fig. 15).3 Group A:σ=64.09(1+n)-0.13

4 Group B:σ=58.63(1+n)-0.17

where: n is the number of freeze-thaw cycles.Figure 15 Fitting curve of uniaxial compressive strength under different number of freeze–thaw cycles. (a) Group A. (b) Group B.

The differences and connections between the prefabricated surface jointed rock mass and the original rock mass can be summarized. Prefabricated surface joints significantly deteriorate the rock mass, leading to a peak uniaxial compressive strength that is significantly smaller than that of the original rock mass. The damage cracks in the prefabricated surface jointed rock mass mainly result in sliding damage along the prefabricated joint surface. From a microscopic perspective, particles in the prefabricated surface jointed rock mass are mostly in direct contact with each other, and damage occurs primarily along existing cracks in the rock mass. As the number of freeze-thaw cycles increases, both groups experience a decrease in uniaxial compressive strength; however, after reaching a certain threshold, the effect of freeze-thaw on strength gradually diminishes.

Therefore, the mechanism of the rock mass freeze-thaw cycle is as follows: when water enters the tiny fractures and pores in the rock mass, it will freeze into ice when the temperature drops to the freezing point. When water and ice coexist, the combination of frost heaving force and hydrostatic pressure continues to act on the rock mass (Fig. 16a). As the water in the pores continues to freeze, the ice will produce about 9%2 volume expansion, generating pressure in the pores and further extending the fractures (Fig. 16b), which can be verified by uniaxial compression tests, where the fractures are formed at the ends of the prefabricated joints. When the temperature rises, as the ice melts, the frost heaving force decreases, and water flows into the new fractures (Fig. 16c), after being subjected to a load, the extended fractures will further expand to become through-going. As water freezes again, the process is repeated, and this repetition leads to the cracking or fragmentation of the rock mass.Figure 16 Schematic diagram of the evolution of freeze–thaw damage mechanisms in rocks.

As the number of freeze-thaw cycles increases, the pores in the rock mass gradually increase. When the ice in the pores is completely frozen, volume expansion will cause fractures to appear under the confinement of the pore space. Therefore, pressure will be released, leading to the breakdown of larger pores and thus reducing the effect of frost heaving forces on the rock mass (Fig. 16d). Therefore, macroscopically, it is observed that with the increase in the number of freeze-thaw cycles, the volumetric expansion pressure is gradually released until the ice expansion force becomes less than the pressure for further expansion of pores This will lead to a trend of gradually decreasing influence of the freeze-thaw action on the strength of the rock mass after the number of cycles reaches a certain range.

By observing the samples after 30 freeze-thaw cycles under scanning electron microscopy, samples are taken at the edge of prefabricated surface joints for testing, and typical images are selected for analysis (Fig. 17). Due to the influence of frost heaving force, as depicted in Fig. 17a, the stacking of rock mass particles becomes loose with expanded fractures. Particle spalling occurs, and new micro-fractures appear between particles. With the increase in magnification, a small number of dissolution pores appear on the rock mass surface (Fig. 17b). In the 10,000× image (Fig. 17c), the dissolution and spalling of rock mass particles are obvious, and the phenomenon of fracture expansion is more significant. This is because the volumetric expansion of the water-ice phase change within the crack, and the repeated action of frost heaving force on the micro-cracks cause the fractures to expand, causing deeper frost damage to the more microscopic rock mass particles.Figure 17 Microstructure of rock specimens after 30 cycles of freeze–thaw action.

Conclusions

This paper conducts freeze-thaw cycle tests, uniaxial compression tests, and scanning electron microscopy examinations on prefabricated surface jointed tuff specimens. The study involves observing macro-micro features before and after freezing and thawing, coupled with mechanical control tests on prefabricated joints and original rock specimens. It focuses on understanding the macro-micro features of fractured rocks post-freezing and thawing and the characteristics of fracture extension damage. The research also delves into the damage mechanisms of rocks due to freeze-thaw cycles, yielding the following key conclusions:Both the mass and wave velocity of the rock samples exhibit a consistent decline with the increasing freeze–thaw cycles. The apparent characteristics of joints and fractures in the rock samples increase. Stress–strain curves for both specimen groups display extended consolidation stages.

The prefabricated jointed specimens are significantly more affected by freeze–thaw action both in fracture numbers and mechanical property, with a substantial decrease in uniaxial compressive strength occurring at the 5th freeze–thaw cycle.

Macroscopic observations reveal a gradual increase in the number of fractures with the cycling. Original rock specimens display multiple sets of wing fractures, while prefabricated surface jointed specimens show anti-wing fractures.

A significant deterioration effect caused by joints on the rock mass is observed. The uniaxial compressive strength of prefabricated jointed specimens is notably smaller than that of the original rock specimens.

After the freeze–thaw cycle, the rock mass's large particles fragment into smaller particles, forming a flocculation structure with aggregated pore particles. The particle connections shift from initial surface cementation contact to point contact under the influence of freeze–thaw cycles.

Acknowledgements

The work was supported by the National Natural Science Foundation of China (42367021), Xinjiang Uygur Autonomous Region Special Program for Key R&D Tasks (2021B03004), Research Base of Xinjiang University, The Innovation Training Program for College Students of Xinjiang Autonomous Region in 2023, National Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep-Earth Engineering.

Author contributions

Conceptualization, R.L. and Z.Z.; Data curation, R.L. and X.L.; Formal analysis, R.L. and J.Z.; Investigation, R.L. and Z.X.; Methodology, R.L. and X.W.; Project administration, Z.Z.; Resources, Z.Z.; Software, J.Z. and X.W.; Supervision, Z.Z.; Validation, Z.X., X.L. and B.X.; Writing—original draft, R.L.; Writing—review & editing, R.L.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

The original online version of this Article was revised: The original version of this Article contained an error in the spelling of the author Runsen Lai which was incorrectly given as Rensen Lai.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Change history

9/18/2024

A Correction to this paper has been published: 10.1038/s41598-024-73024-4
==== Refs
References

1. Wettlaufer JS Worster MG Premelting dynamics Annu. Rev. Fluid Mech. 2006 38 427 452 10.1146/annurev.fluid.37.061903.175758
Wettlaufer, J. S. & Worster, M. G. Premelting dynamics. Annu. Rev. Fluid Mech. 38, 427–452. 10.1146/annurev.fluid.37.061903.175758 (2006).
2. Setzer, M. J. Development of the micro-ice-lens model. In International RILEM Workshop on Frost Resistance of Concrete (RILEM Publications SARL, 2002), 133–145.
3. Liu Q Huang S Kang Y Fractured rock advance and review on freezing-thawing damage of fractured rock Chin. J. Rock Mech. Eng. 2015 34 03 452 471 10.13722/j.cnki.jrme.2015.03.003
Liu, Q. et al. Fractured rock advance and review on freezing-thawing damage of fractured rock. Chin. J. Rock Mech. Eng. 34(03), 452–471. 10.13722/j.cnki.jrme.2015.03.003 (2015).
4. Rodríguez P Arab PB Celestino TB Characterization of rock cracking patterns in diametral compression tests by acoustic emission and petrographic analysis Int. J. Rock Mech. Min. Sci. 2016 83 73 85 10.1016/j.ijrmms.2015.12.017
Rodríguez, P., Arab, P. B. & Celestino, T. B. Characterization of rock cracking patterns in diametral compression tests by acoustic emission and petrographic analysis. Int. J. Rock Mech. Min. Sci. 83, 73–85. 10.1016/j.ijrmms.2015.12.017 (2016).
5. Bao W Lu H Guo Q Study on the freeze-thaw degradation characteristics of the surrounding rock of the high cold carbonaceous slate tunnel in Xinjiang J. Eng. Geol. 2023 31 04 1213 1224 10.13544/j.cnki.jeg.2023-0102
Bao, W. et al. Study on the freeze-thaw degradation characteristics of the surrounding rock of the high cold carbonaceous slate tunnel in Xinjiang. J. Eng. Geol. 31(04), 1213–1224. 10.13544/j.cnki.jeg.2023-0102 (2023).
6. Pan PZ Miao S Jiang Q The influence of infilling conditions on flaw surface relative displacement induced cracking behavior in hard rock Rock Mech. Rock Eng. 2020 53 4449 4470 10.1007/s00603-019-02033-x
Pan, P. Z. et al. The influence of infilling conditions on flaw surface relative displacement induced cracking behavior in hard rock. Rock Mech. Rock Eng. 53, 4449–4470 (2020).
7. Li BQ da Silva BG Einstein H Laboratory hydraulic fracturing of granite: acoustic emission observations and interpretation Eng. Fract. Mech. 2019 209 200 220 10.1016/j.engfracmech.2019.01.034
Li, B. Q., da Silva, B. G. & Einstein, H. Laboratory hydraulic fracturing of granite: acoustic emission observations and interpretation. Eng. Fract. Mech. 209, 200–220. 10.1016/j.engfracmech.2019.01.034 (2019).
8. Wang M Zhu Z Feng R Experimental study on failure behaviour of cracked rocks around tunnels under true triaxial unloading J. China Coal Soc. 2015 40 02 278 285 10.13225/j.cnki.jccs.2014.0299
Wang, M., Zhu, Z. & Feng, R. Experimental study on failure behaviour of cracked rocks around tunnels under true triaxial unloading. J. China Coal Soc. 40(02), 278–285. 10.13225/j.cnki.jccs.2014.0299 (2015).
9. Zhu ZM An alternative form of propagation criterion for two collinear cracks under compression Math. Mech. Solids 2009 14 8 727 746 10.1177/1081286508090043
Zhu, Z. M. An alternative form of propagation criterion for two collinear cracks under compression. Math. Mech. Solids 14(8), 727–746. 10.1177/1081286508090043 (2009).
10. Zhu ZM Evaluation of the range of horizontal stresses in the earth's upper crust by using a collinear crack model J. Appl. Geophys. 2013 2013 88 114 121 10.1016/j.jappgeo.2012.10.007
Zhu, Z. M. Evaluation of the range of horizontal stresses in the earth’s upper crust by using a collinear crack model. J. Appl. Geophys. 2013(88), 114–121. 10.1016/j.jappgeo.2012.10.007 (2013).
11. Guo Y Wong RHC Zhu W Study on fracture pattern of open surface-flaw in gabbro Chin. J. Rock Mech. Eng. 2007 03 525 531
Guo, Y. et al. Study on fracture pattern of open surface-flaw in gabbro. Chin. J. Rock Mech. Eng. 03, 525–531 (2007).
12. Dyskin AV Germanovich LN Ustinov KB A 3-D model of wing crack growth and interaction Eng. Fract. Mech. 1999 63 1 81 110 10.1016/s0013-7944(96)00115-4
Dyskin, A. V., Germanovich, L. N. & Ustinov, K. B. A 3-D model of wing crack growth and interaction. Eng. Fract. Mech. 63(1), 81–110. 10.1016/s0013-7944(96)00115-4 (1999).
13. Wong RHC Huang ML Jiao MR The mechanisms of crack propagation from surface 3-D fracture under uniaxial compression Key Eng. Mater. 2004 261–263 219 224 10.4028/www.scientific.net/kem.261-263.219
Wong, R. H. C. et al. The mechanisms of crack propagation from surface 3-D fracture under uniaxial compression. Key Eng. Mater. 261–263, 219–224. 10.4028/www.scientific.net/kem.261-263.219 (2004).
14. Wong RHC Law CM Chau KT Crack propagation from 3-D surface fractures in PMMA and marble specimens under uniaxial compression Int. J. Rock Mech. Min. Sci. 2004 41 37 42 10.1016/j.ijrmms.2004.03.016
Wong, R. H. C. et al. Crack propagation from 3-D surface fractures in PMMA and marble specimens under uniaxial compression. Int. J. Rock Mech. Min. Sci. 41, 37–42. 10.1016/j.ijrmms.2004.03.016 (2004).
15. Nakamura D Goto T Suzuki T Ito Y Basic study on the frost heave pressure of rocks: dependence of the location of frost heave on the strength of the rock Cold Regions Eng. 2012 10.1061/978078441247
Nakamura, D. et al. Basic study on the frost heave pressure of rocks: dependence of the location of frost heave on the strength of the rock. Cold Regions Eng.10.1061/978078441247 (2012).
16. Huang M Wong RHC Experimental study on propagation and coalescenceme chanisms of 3D surface cracks Chin. J. Rock Mech. Eng. 2007 2007 09 1794 1799
Huang, M. & Wong, R. H. C. Experimental study on propagation and coalescenceme chanisms of 3D surface cracks. Chin. J. Rock Mech. Eng. 2007(09), 1794–1799 (2007).
17. Liu T Yang R Ding L Study on mechanical properties and meso-initiation mechanismof granite with non-persistent joints Chin. J. Rock Mech. Eng 2022 42 05 1070 1082 10.13722/j.cnki.jrme.2022.0671
Liu, T. et al. Study on mechanical properties and meso-initiation mechanismof granite with non-persistent joints. Chin. J. Rock Mech. Eng 42(05), 1070–1082. 10.13722/j.cnki.jrme.2022.0671 (2022).
18. Wang, C., Wang, L., & Liu, J., et al. Investigations on the mechanical behavior and failure mechanism of singlefractured granite under uniaxial compression. J. Southwest Jiaotong Univ. 1–9 (2024)
19. Meng F Zhou H Wang Z Characteristics of asperity damage and its influence on the shear behavior of granite joints Rock Mech. Rock Eng. 2017 51 2 429 449 10.1007/s00603-017-1315-y
Meng, F. et al. Characteristics of asperity damage and its influence on the shear behavior of granite joints. Rock Mech. Rock Eng. 51(2), 429–449. 10.1007/s00603-017-1315-y (2017).
20. Mu J Pei X Huang R Degradation characteristics of shear strength of joints in three rock types due to cyclic freezing and thawing Cold Reg. Sci. Technol. 2017 138 91 97 10.1016/j.coldregions.2017.03.011
Mu, J. et al. Degradation characteristics of shear strength of joints in three rock types due to cyclic freezing and thawing. Cold Reg. Sci. Technol. 138, 91–97. 10.1016/j.coldregions.2017.03.011 (2017).
21. Lu Y Li X Chan A Damage constitutive model of single flaw sandstone under freeze-thaw and load Cold Reg. Sci. Technol. 2018 10.1016/j.coldregions.2018.11.017
Lu, Y., Li, X. & Chan, A. Damage constitutive model of single flaw sandstone under freeze-thaw and load. Cold Reg. Sci. Technol.10.1016/j.coldregions.2018.11.017 (2018).
22. Wang Z Yan E Influence of material composition and structural characteristics of rock onfreeze-thaw damage and deterioration of schist Chin. J. Geotech. Eng. 2015 37 S2 86 90 10.11779/CJGE2015S2018
Wang, Z. & Yan, E. Influence of material composition and structural characteristics of rock onfreeze-thaw damage and deterioration of schist. Chin. J. Geotech. Eng. 37(S2), 86–90. 10.11779/CJGE2015S2018 (2015).
23. Xiao P Chen Y Du X Mechanical properties of sandstone under freeze-thaw cycles andstudies on meso-damage constitutive model Chin. J. Geotech. Eng. 2023 45 04 805 815 10.11779/CJGE20220219
Xiao, P. et al. Mechanical properties of sandstone under freeze-thaw cycles andstudies on meso-damage constitutive model. Chin. J. Geotech. Eng. 45(04), 805–815. 10.11779/CJGE20220219 (2023).
24. Yang M Zhang J Meng X Experimental study on pore evolution of red sandstone at different scales under freeze-thaw cycles Min. Res. Dev. 2023 43 11 102 106 10.13827/j.cnki.kyyk.2023.11.019
Yang, M., Zhang, J. & Meng, X. Experimental study on pore evolution of red sandstone at different scales under freeze-thaw cycles. Min. Res. Dev. 43(11), 102–106. 10.13827/j.cnki.kyyk.2023.11.019 (2023).
25. Ia H Zhao S Ding S Study on the evolution and influencing factors of frost heaving force of waterbearing cracks during freezing-thawing process Chin. J. Rock Mech. Eng. 2022 41 09 1832 1845 10.13722/j.cnki.jrme.2021.1147
Ia, H. et al. Study on the evolution and influencing factors of frost heaving force of waterbearing cracks during freezing-thawing process. Chin. J. Rock Mech. Eng. 41(09), 1832–1845. 10.13722/j.cnki.jrme.2021.1147 (2022).
26. Liu H Yang G Shen Y CT visual quantitative characterization of meso-damage evolution of sandstoneunder freeze-thaw-loading synergistic effect Chin. J. Rock Mech. Eng. 2023 42 05 1136 1149 10.13722/j.cnki.jrme.2022.0574
Liu, H. et al. CT visual quantitative characterization of meso-damage evolution of sandstoneunder freeze-thaw-loading synergistic effect. Chin. J. Rock Mech. Eng. 42(05), 1136–1149. 10.13722/j.cnki.jrme.2022.0574 (2023).
27. Yang H Liu P Sun B Study on damage mechanisms of the microstructure of sandy conglomerate atMaijishan grottoes under freeze-thaw cycles Chin. J. Rock Mech. Eng. 2021 40 03 545 555 10.13722/j.cnki.jrme.2020.0767
Yang, H. et al. Study on damage mechanisms of the microstructure of sandy conglomerate atMaijishan grottoes under freeze-thaw cycles. Chin. J. Rock Mech. Eng. 40(03), 545–555. 10.13722/j.cnki.jrme.2020.0767 (2021).
28. Lai R Zhang Z Chen D Collapse development characteristics and failure mechanism of Altay Camel Peak in Xinjiang J. Eng. Geol. 2023 31 4 1186 1199 10.13544/j.cnki.jeg.2023-0141
Lai, R. et al. Collapse development characteristics and failure mechanism of Altay Camel Peak in Xinjiang. J. Eng. Geol. 31(4), 1186–1199. 10.13544/j.cnki.jeg.2023-0141 (2023).
29. Teng C Yin X Li S Experimental study of three dimensionalfracture in plate specimens withnonpenetrating crack Chin. J. Geophys. 1987 1987 04 371 378
Teng, C. et al. Experimental study of three dimensionalfracture in plate specimens withnonpenetrating crack. Chin. J. Geophys. 1987(04), 371–378 (1987).
30. Chen S Feng X Li S Effects of chemical erosion on uniaxial compressivestrength and meso-fracturing behaviors of rock Chin. J. Rock Mech. Eng. 2003 2003 04 547 551
Chen, S., Feng, X. & Li, S. Effects of chemical erosion on uniaxial compressivestrength and meso-fracturing behaviors of rock. Chin. J. Rock Mech. Eng. 2003(04), 547–551 (2003).
31. Yu M Xu G Tang G Experimental study on stepped unloading of rock-ike materialswith pre existing fissure Chin. J. Underground Space Eng. 2020 16 06 1672 1681
Yu, M. et al. Experimental study on stepped unloading of rock-ike materialswith pre existing fissure. Chin. J. Underground Space Eng. 16(06), 1672–1681 (2020).
32. Zhao Y Wan W Wang W Compressive-shear rheological fracture of rock-like cracks and subcriticalcrack propagation test and fracture mechanism Chin. J. Geotech. Eng. 2012 34 06 1050 1059
Zhao, Y. et al. Compressive-shear rheological fracture of rock-like cracks and subcriticalcrack propagation test and fracture mechanism. Chin. J. Geotech. Eng. 34(06), 1050–1059 (2012).
33. Hosseini M Khodayari AR Effect of freeze-thaw cycle on strength and rock strength parameters (A Lushan sandstone case study) J. Min. Environ. 2019 10 1 257 270
Hosseini, M. & Khodayari, A. R. Effect of freeze-thaw cycle on strength and rock strength parameters (A Lushan sandstone case study). J. Min. Environ. 10(1), 257–270 (2019).
34. Momeni A Abdilor Y Khanlari GR The effect of freeze–thaw cycles on physical and mechanical properties of granitoid hard rocks Bull. Eng. Geol. Environ. 2016 2016 75 1649 1656 10.1007/s10064-015-0787-9
Momeni, A. et al. The effect of freeze–thaw cycles on physical and mechanical properties of granitoid hard rocks. Bull. Eng. Geol. Environ. 2016(75), 1649–1656 (2016).
35. Lin H Lei D Yong R Analytical and numerical analysis for frost heaving stress distribution within rock joints under freezing and thawing cycles Environ. Earth Sci. 2020 2020 79 1 17
Lin, H. et al. Analytical and numerical analysis for frost heaving stress distribution within rock joints under freezing and thawing cycles. Environ. Earth Sci. 2020(79), 1–17 (2020).
36. Liu H Huang Y Li K Test study of strength and failure mode of pre-existing jointed rock mass Rock Soil Mech. 2013 34 05 1235 1241 10.16285/j.rsm.2013.05.002
Liu, H. et al. Test study of strength and failure mode of pre-existing jointed rock mass. Rock Soil Mech. 34(05), 1235–1241. 10.16285/j.rsm.2013.05.002 (2013).
37. Yanlin Z Wen W Weijun W Fracture experiments on ordered multi-crack body in rock-like materialsunder uniaxial compression and numerical simulation of wing cracks Chin. J. Geotech. Eng. 2013 35 11 2097 2109
Yanlin, Z. et al. Fracture experiments on ordered multi-crack body in rock-like materialsunder uniaxial compression and numerical simulation of wing cracks. Chin. J. Geotech. Eng. 35(11), 2097–2109 (2013).
38. Hao J Qiao L Li Z Study on the fracture behavior of prefabricated fissures granite based on DIC and laser scanning techniques Fatigue Fract. Eng. Mater. Struct. 2021 44 5 1372 1390 10.1111/ffe.13435
Hao, J. et al. Study on the fracture behavior of prefabricated fissures granite based on DIC and laser scanning techniques. Fatigue Fract. Eng. Mater. Struct. 44(5), 1372–1390 (2021).
