Next Article in Journal
Extreme Waves
Previous Article in Journal
Dynamic Appointment Rescheduling of Trucks under Uncertainty of Arrival Time
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Damage Model of Rock Based on Compaction Stage and Post-Peak Shape under Chemical-Freezing-Thawing-Loading

1
School of Resources and Environmental Engineering, Wuhan University of Science and Technology, Wuhan 430081, China
2
Hubei Key Laboratory for Efficient Utilization and Agglomeration of Metallurgic Mineral Resources, Wuhan 430081, China
3
Huangshi Construction Market Management Station, Huangshi 435000, China
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(5), 696; https://0-doi-org.brum.beds.ac.uk/10.3390/jmse10050696
Submission received: 28 March 2022 / Revised: 13 May 2022 / Accepted: 17 May 2022 / Published: 19 May 2022
(This article belongs to the Section Ocean Engineering)

Abstract

:
The deterioration of rock mechanical properties under chemical-freezing-thawing-loading will seriously threaten the stability and safety of engineering rock mass, so the study of its constitutive model has very important theoretical and engineering application significance. In view of the deficiency that the existing statistical damage constitutive model cannot describe the nonlinear characteristics of the compaction stage and the post-peak shape, the compaction index is introduced to measure the stress-strain nonlinear characteristics caused by the iso-compression of the pores in the rock, and the post-peak correction coefficient is introduced to optimize the post-peak shape of rock loading. Assuming that the strength and failure of rock microunits obey the Weibull distribution, on the basis of using the Drucker-Prager strength criterion to measure the strength of rock microunits, a statistical damage model for rock considering the compaction stage and post-peak shape under chemical-freezing-thawing-loading is established. Finally, it is compared with the test curve. The results show that the calculation curve of the established chemical-freezing-thawing-loading rock statistical damage model has the same trend as the rock uniaxial compression test curve. Furthermore, it can better describe the rock stress-strain law with different chemical solutions and different freeze-thaw cycles, which fully reflects the rationality and accuracy of the constructed constitutive model. The research results can provide a theoretical basis for the calculation of deformation and failure of rocks under chemical-freezing-thawing-loading.

1. Introduction

On the one hand, engineering rock masses in cold regions are chemically eroded by water, acid rain and many ions in groundwater, causing changes in the mineral composition and structure inside the rock [1,2]. On the other hand, under the influence of seasons and alternation of day and night, rocks are damaged by freeze-thaw cycles under repeated freezing and thawing of water and chemical solutions [3,4,5], and the rock mass is always under a certain stress field condition. Therefore, under the action environment of chemical erosion, freeze-thaw cycles and loads, the internal microscopic defects of the rock gradually generate and expand, and the damage accumulates, thereby deteriorating the physical and mechanical properties of the rock mass. This affects the safety of rock mass engineering in cold regions at all times and it makes the research on the mechanical damage of rock under chemical-freezing-thawing-loading extremely important in theory and engineering practice. Establishing an objective and reasonable constitutive model is one of the keys to explore and study the mechanical properties of rock under chemical-freezing-thawing-loading.
A large number of scholars first carried out rock mechanics test and damage constitutive theory research under a single chemical erosion or freeze-thaw environment, and made useful progress. In terms of chemical erosion tests, Hutchinson et al. [6], Cai et al. [7], Li et al. [8], and Gong et al. [9] carried out experimental research on the corrosion mechanics of limestone and sandstone with different chemical solutions, and the results showed that chemical erosion will change the mineral composition of the rock and affect the parameters such as uniaxial compressive strength and elastic modulus. Han et al. [10] carried out chemical corrosion sandstone mechanical tests and proposed to use porosity changes to describe the microstructure damage variables. In terms of freeze-thaw tests, Yang et al. [11] studied the damage extension mechanism of rocks, water migration, ice formation and structural damage changes under freeze-thawcycles and discussed the effects of freeze-thaw cycles on the deterioration of physical and mechanical properties of rocks. Zhou et al. [12] first performed freeze-thaw cycles on sandstone, and then conducted NMR tests and impact load tests to analyze the microscopic damage characteristics and dynamic mechanical parameters of sandstone. Yahaghi et al. [13] studied the physical and mechanical properties and failure behavior of Tasmanian sandstone in various freeze-thaw cycles through a P-wave test, uniaxial compression test and a Brazilian split test after freeze-thaw cycles. Li et al. [14] carried out research on rock freeze-thaw cycles and uniaxial compression tests, analyzed the freeze-thaw degradation mechanism, and applied the test results to the stability evaluation of rock mass slopes in open-pit mines. The above laboratory experiments show that chemical erosion and freeze-thaw cycles have damage effects on rocks, which are mainly reflected in the reduction of elastic modulus, mechanical strength, and peak strain increase.
In the study of the constitutive model of chemical erosion, Chen [15] established a microscopic damage evolution constitutive model of rock under uniaxial compression under different pH values, concentrations, and types of chemical solution corrosion by introducing chemical corrosion damage factors. Chen et al. [16] proposed an evolutionary neural network constitutive model to describe the stress-strain characteristics of rocks under chemical corrosion. Lin et al. [17] introduced the damage variable correction coefficient affected by porosity and established a statistical damage constitutive model of rock under chemical corrosion based on the Drucker-Prager strength criterion and the two-parameter Weibull distribution. In terms of freeze-thaw cycles, Fang et al. [3] established a statistical damage constitutive model for rocks under freeze-thaw and load conditions, and clarified the method for determining model parameters, and analyzed the relationship between elastic modulus, peak stress, and freeze-thaw cycles. Huang et al. [18] derived a statistical damage constitutive model for rocks under freezing-thawing-loading conditions based on Weibull distribution and maximum tensile strain yield criterion and applied the model to analyze the stability of tunnels under thermal-water-mechanical coupling conditions in cold regions. Feng et al. [19] carried out the freezing-thawing cycle test of saturated sandstone, established a freezing-thawing damage model based on the energy evolution law, and used this model to perform theoretical calculations on the uniaxial compression test. Zhang et al. [20] proposed the concepts of freeze-thaw damage, load damage and total damage, established a freezing-thawing-loading rock damage model, and described the damage evolution and macroscopic damage behavior of the meso-structure of rock materials through damage variables and constitutive equations.
To sum up, the damage research on single chemical and freeze-thaw factors has achieved relatively fruitful results, but most of the rocks in the actual engineering rock are subject to two or more damage effects, which makes the current theoretical research and engineering applications still have a large gap. In recent years, Chen [21,22,23,24], Wang [25], and Li et al. [26] carried out experimental studies on the mechanical properties of rocks immersed in different chemical solutions and subjected to freeze-thaw cycles. Li et al. [27] studied the pore change law and damage evolution mechanism of sandstone under chemical solution erosion by nuclear magnetic resonance and carried out sandstone freeze-thaw cycle tests with different solutions, indicating that the combined effect of chemical corrosion and freeze-thaw cycles can relatively slow down the damage aggravation. Due to limited test conditions and complex influencing factors, the research on rock mechanics under multi-field coupling is in its infancy, and there is no report on the research on the constitutive model of rock under the combined action of chemical-freezing-thawing-loading. In addition, the traditional constitutive model still has some shortcomings. First, the traditional constitutive model simplifies the compaction stage as a linear feature for processing, and the chemical-freezing-thawing action leads to a large increase in the internal pores of the rock, the compaction effect of rock is much more significant during the mechanical test than in the natural state, so the nonlinear change of compaction stage cannot be ignored for the constitutive model of rock under chemical-freezing-thawing-loading. Second, the degree of fitting of the post-peak shape is relatively low. Cao et al. [28], Zhang et al. [29], Zhang et al. [30], and Bilen et al. [31] carried out research on the evaluation of brittleness and plasticity indicators in the process of rock damage. It is worth noting that the brittleness of the rock is weakened and the plasticity is enhanced after chemical erosion and freeze-thaw cycles. The traditional constitutive model cannot well reflect the post-peak shape changes of the plastic enhancement caused by chemical and freeze-thaw damage.
To sum up, the current rock environment is not considered enough, and the existing constitutive model cannot better describe the whole process of rock deformation and failure, resulting in great differences between the current model and the rock stress-strain curve in engineering practice. Therefore, based on previous research, this paper first considers the coupling effect of chemical erosion and freeze-thaw cycles, assumes that the strength of rock micro units obeys the Weibull distribution, and introduces the compaction index, a post-peak correction coefficient, threshold stress, and threshold strain. A statistical damage constitutive model of rock under chemical-freezing-thawing-loading is established. Then, it is compared with the rock uniaxial compression test curve of the literature [21] to verify the rationality and accuracy of the established model. Finally, the total damage evolution equation under the action of chemical-freezing-thawing-loading is used to describe the rock meso-damage evolution and macro-mechanical performance response, which provides a theoretical reference for rock multi-field coupled damage research.

2. Statistical Damage Model of Rock under Chemical-Freezing-Thawing-Loading

2.1. Chemical-Freeze-Thaw Damage Variable and Damage Model

Tang et al. [32] put forward the concept of chemical damage earlier and defined it as the deterioration process of rock (mass) material or structure caused by water-rock chemistry. Water-rock chemistry affects the physical state and micro-structure of rocks, weakens the connection between mineral particles, increases rock mass deformation, decreases strength, and significantly reduces the elastic modulus. Because the elastic modulus of rock is easier to analyze and measure during chemical erosion, the damage variable D C of rock chemical erosion can be expressed as:
D C = 1 E C E 0
where E C is the elastic modulus of rock after chemical erosion damage, E 0 is the elastic modulus of natural rock before chemical erosion.
Due to the freeze-thaw cycles, a large number of meso-cracks are bound to occur in the rock, and the cracks gradually expand with the increase of the number of cycles, and the elastic modulus also decreases significantly, indicating that the freeze-thaw cycles cause damage to the rock and so the elastic modulus is also selected in this paper. To represent the freeze-thaw damage variable, so the freeze-thaw damage variable D T n is defined as
D T n = 1 E T n E 0
where E T n is the elastic modulus after n times of freeze-thaw, and E 0 is the initial elastic modulus without freeze-thaw.
Figure 1 is a schematic diagram of the damage state of the rock, where a is the rock in the initial state, b is the rock after chemical erosion, c is the rock that has undergone chemical-freezing-thawing action, and d is the rock that has been subjected to chemical-freezing-thawing-loading. Based on the generalized Lemaitre strain equivalence principle [33], the natural initial state of the rock is regarded as the first damage state, and the state after chemical damage is regarded as the second damage state, then
σ C = σ 0 ( 1 D C )
The generalized strain equivalence principle is applied again, chemical damage is regarded as the first damage state, and the total damage state caused by chemical-freezing-thawing action is regarded as the second damage state, then
σ T = σ C ( 1 D T n )
The damage model of chemical-freezing-thawing action is
σ T = σ 0 ( 1 D C ) ( 1 D T n ) = E 0 ε ( 1 D C ) ( 1 D T n ) = E ε
where E is the elastic modulus after chemical-freezing-thawing action; σ 0 is the effective stress in the natural initial state; σ C is the effective stress in the chemical damage state; σ T is the effective stress in the chemical-freezing-thawing state.

2.2. Statistical Distribution Function

Rock is a natural material composed of a variety of minerals. Due to the heterogeneity of the microscopic structure inside the rock arising from the geological environment and mineral composition, the mechanical properties of the micro units that make up the rock may have large differences, and the micro cracks and fractures are also randomly distributed. It is a reasonable and effective way to use statistical functions to describe the strength of rock micro units. At present, the commonly used distribution functions in rock statistical damage constitutive models include the Weibull distribution, normal distribution [34], power function distribution, lognormal distribution, etc. According to the literature [35,36], the Weibull probability distribution function has the characteristics of peak effect and simple integral calculation, which can better describe the process of rock deformation and failure. Therefore, this paper assumed that the generation and expansion of cracks under the action of chemical-freezing-thawing obey the Weibull distribution function, and the strength of micro units also obeys the Weibull statistical distribution, and the probability density function p ( x ) is
p ( x ) = m F ( x F ) m 1 exp ( x F ) m
The distribution function P ( x ) is
P ( x ) = 1 exp ( x F ) m
where x is the micro units’ strength value of the rock under chemical-freezing-thawing, m and F are parameters related to the Weibull distribution function affecting the shape and size of rock micro units, which change with the change of chemical and freeze-thaw conditions. Therefore, chemical and freeze-thaw damage will affect parameters m and F . According to the existing research, this paper assumes that m ( C T ) and F ( C T ) are similarly affected by chemical and freeze-thaw damage, and also according to the extended Lemaitre strain equivalence principle mentioned above, we can obtain
m ( C ) = m 0 ( 1 D C ) F ( C ) = F 0 ( 1 D C )
m ( C T ) = m ( C ) ( 1 D T n ) F ( C T ) = F ( C ) ( 1 D T n )
where m 0 and F 0 are the Weibull parameters of the rock without chemical erosion, m ( C ) and F ( C ) are the Weibull parameters of the rock after chemical erosion, m ( C T ) and F ( C T ) are the Weibull parameters of the rock under different chemical and freeze-thaw conditions.
According to Equations (8) and (9), the Weibull parameters of rocks under different chemical erosion and freeze-thaw cycles are
m ( C T ) = m 0 ( 1 D C ) ( 1 D T n ) F ( C T ) = F 0 ( 1 D C ) ( 1 D T n )
Assuming that the total number of micro units in the rock is N , the number of damaged micro units under a certain stress level is N F , and the damage variable D is defined as the ratio of the number of damaged micro units to the total number of micro units, namely
D = N F N
When the stress is loaded to the p ( σ ) state, the micro units start to be destroyed one after another, and the number of micro units that have been destroyed is
N F = 0 p ( σ ) N · p ( σ ) d σ = N 0 p ( σ ) p ( σ ) d σ
From Equations (11) and (12), the rock damage variable D can be obtained as
D = N 0 p ( σ ) p ( σ ) d σ N = 0 p ( σ ) p ( σ ) d σ
Substituting Equations (7) and (10) into Equation (13), we can get
D = 1 exp ( p ( σ ) F 0 ( 1 D C ) ( 1 D T n ) ) m 0 ( 1 D C ) ( 1 D T n )
There is an obvious threshold for rock damage under loading [37]. When the stress state is lower than the threshold stress, the stress-strain relationship of the rock is mainly in the stage of pore compaction and linear elasticity, and the damage caused by the load to the interior of the rock is almost zero and can be ignored. At this time, the total damage is caused by chemical erosion and freezing-thawing, therefore, the damage variable is
D = 1 ( 1 D C ) ( 1 D T n ) = 1 E C E 0 . E T n E 0
Combining Equations (14) and (15), the damage evolution equation in the process of rock loading under chemical-freezing-thawing action can be obtained as
D = 1 ( 1 D C ) ( 1 D T n ) ( 0 ε ε u ) 1 exp ( p ( σ ) F 0 ( 1 D C ) ( 1 D T n ) ) m 0 ( 1 D C ) ( 1 D T n ) ( ε > ε u )
where ε u is the strain at the rock damage threshold point, Equation (16) represents the relationship between the total rock damage D , chemical damage D C , and freeze-thaw damage D T n under loading conditions.
According to Lemaitre’s strain equivalence principle [38] and the concept of effective stress, the strain produced by the nominal stress (stress measured by the test) on the damaged material is equivalent to the strain caused by the effective stress on the non-destructive material. Moreover, the post-peak shapes of the stress-strain curves of different rocks are not exactly the same during the loading process. For example, basalt, granite, and quartzite have high strength, and the stress decreases approximately vertically in the post-peak stage, which is prone to brittle failure. On the other hand, shale and sandstone have low strength, and the stress decreases slowly in the post-peak stage, showing good plasticity [39]. In order to be closer to the actual stress-strain law, this article introduces the post-peak correction coefficient λ to represent the plasticity, and the corrected damage model is
σ i = σ i * ( 1 λ D )     ( i = 1 , 2 , 3 )
where σ i is the nominal stress of the rock (MPa), and σ i * is the effective stress of the rock (MPa).
Due to the obvious elastic characteristics of rock stress-strain under chemical-freezing-thawing action, according to the generalized Hooke’s Law, it can be obtained
ε i = 1 E σ i * μ ( σ j * + σ k * )       ( i , j , k = 1 , 2 , 3 )
where ε i is the strain; μ is the Poisson’s ratio; E is the elastic modulus.
Substitute Equation (17) into Equation (18) to get
ε i = 1 E ( 1 λ D ) σ i μ ( σ j + σ k )         ( i , j , k = 1 , 2 , 3 )
when ε i > ε u , substitute Equation (14) into Equation (19) to get
ε i = 1 E 1 λ + λ exp ( p ( σ ) F 0 ( 1 D C ) ( 1 D T n ) ) m 0 ( 1 D C ) ( 1 D T n ) σ i μ ( σ j + σ k )                   ( i , j , k = 1 , 2 , 3 )
Sorted out
σ i = E ε i 1 λ + λ exp ( p ( σ ) F 0 ( 1 D C ) ( 1 D T n ) ) m 0 ( 1 D C ) ( 1 D T n ) + μ ( σ j + σ k )             ( i , j , k = 1 , 2 , 3 ; ε i > ε u )
Substituting Equation (5) into Equation (21), the rock damage equation considering chemical-freezing-thawing-loading can be obtained as:
σ i = E 0 ( 1 D C ) ( 1 D T n ) ε i 1 λ + λ exp ( p ( σ ) F 0 ( 1 D C ) ( 1 D T n ) ) m 0 ( 1 D C ) ( 1 D T n ) + μ ( σ j + σ k )                   ( i , j , k = 1 , 2 , 3 ; ε i > ε u )

2.3. The Representation of Micro Units’ Strength

At present, many scholars have introduced the commonly used rock strength criteria, such as the maximum tensile strain strength criterion, the Mohr–Coulomb strength criterion, the Drucker–Prager strength criterion and the Hoek-Brown strength criterion, etc. The Drucker–Prager strength criterion considers the effect of intermediate principal stress and the effect of hydrostatic pressure, which makes up for the main defect of the Mohr–Coulomb strength criterion. In view of this, this paper adopts the Drucker–Prager strength failure criterion as the measurement method of the strength of the rock micro units, and the strength of the micro units p ( σ ) is expressed as follows
p ( σ ) = α 0 I 1 + J 2
α 0 = sin φ 9 + 3 sin 2 φ
where φ is the rock internal friction angle; I 1 is the first invariant of the stress tensor; J 2 is the second invariant of the stress deviator.
It can be obtained by combining Equation (19)
σ 1 * = σ 1 1 λ D           = σ 1 ε 1 E σ 1 μ ( σ 2 + σ 3 )
σ 2 * = σ 2 1 λ D           = σ 2 ε 1 E σ 2 μ ( σ 1 + σ 3 )
σ 3 * = σ 3 1 λ D           = σ 3 ε 1 E σ 3 μ ( σ 1 + σ 2 )
Then,
I 1 = σ 1 * + σ 2 * + σ 3 * = E ε 1 ( σ 1 + 2 σ 3 ) σ 1 2 μ σ 3
J 2 = 1 6 ( σ 1 * σ 2 * ) 2 + ( σ 2 * σ 3 * ) 2 + ( σ 3 * σ 1 * ) 2 = E ε 1 ( σ 1 σ 3 ) 3 ( σ 1 2 μ σ 3 )
Substitute Equations (28) and (29) into Equation (23) to get
p ( σ ) = α 0 E ε 1 ( σ 1 + 2 σ 3 ) σ 1 2 μ σ 3 + E ε 1 ( σ 1 σ 3 ) 3 ( σ 1 2 μ σ 3 )

2.4. Statistical Damage Constitutive Model Considering Chemical-Freezing-Thawing Effect in Compaction Stage

The entire deformation process of the rock sample under uniaxial compression can be represented by the full stress-strain curve shown in Figure 2 [40]. This shows that the deformation and failure of the rock is divided into four stages, which are the pore and fracture compaction stage (OA stage), elastic deformation to the stable development stage of micro-elastic cracks (AC stage), among them, the AB stage is the elastic deformation stage, and the BC stage is the stable development stage of micro-elastic cracks, unstable rupture development stage (CD stage), and post-rupture stage (after point D).
Before the peak, stress generally increases with increasing strain, but the rate of stress growth is not exactly the same. In the compaction stage, due to the compaction of the internal pores of the rock, the stress increases slowly with the strain, showing an upward concave shape; then in the elastic deformation stage, the stress increases rapidly with the strain and is basically linear elastic. To establish a constitutive model that conforms to the nonlinear characteristics of rock deformation in Figure 2, the existence of pores in the rock needs to be considered, and under the action of chemical erosion and freeze-thaw, the rock crystal plane is prone to dissolve small cavities, resulting in many secondary cracks and pores. The original pores in the rock sample are further expanded or more micro-cavities are formed, which promotes changes in the deformation and failure mechanism of the rock. To sum up, the compaction effect of rock pores is more obvious under the action of chemical-freezing-thawing, and it is particularly important to simulate the nonlinear characteristics of the compaction stage. In addition, chemical erosion and freeze-thaw cycles will also destroy the connection between the particles in the rock and the crystal itself, so that the rock changes from brittleness to plasticity, that is, the brittleness of the rock decreases and the plasticity increases, and the stress-strain curve of the rock sample fall slowly after the peak.
After multiple fittings, it was found that the exponential function can well describe the stress-strain characteristics considering the pore compaction in the OC stage. This paper defines a parameter α to measure the pore compaction of the rock, which is called the compaction index, and its value is related to the rock type, damage and stress state. It can be expressed as
τ = σ u ( ε ε u ) α
where σ u is the stress value at the ultimate elastic strain ε u .
The influence of the compaction index on the shape of the stress-strain curve is shown in Figure 3. The boundary yield point of the compaction stage, the elastic deformation to the stable development stage of micro-elastic cracks (OC stage) and the unstable rupture development stage (CD stage) is C ( ε u , σ u ) , the peak point is D ( ε s , σ s ) , where the peak strain is ε s and the peak stress is σ s . With the increase of α, the degree of concave on the curve is more obvious, α < 1, indicating that there are fewer pores in the rock, when α = 1, the pores in the rock are moderate, that is, the stress-strain curve increases linearly, on the contrary, when α > 1, there are more pores in the rock.
When the stress state is lower than the threshold stress, the rock is mainly in the stage of pore compaction and linear elasticity, and the total damage is caused by chemical erosion and freeze-thaw cycles, the constitutive relationship of the rock can be well obtained by fitting the above exponential function. When the stress state exceeds the threshold stress, the damage is caused by the combined action of chemistry, freeze-thaw, and load, which can be described by the Weibull distribution function. Combining Equations (22) and (31), the statistical damage constitutive model of rock based on the Weibull distribution function and the Drucker–Prager strength criterion under chemical-freezing-thawing-loading can be obtained as,
σ i = σ u ( ε ε u ) α                                                                                                                                                                                                         ( 0 ε ε u ) E 0 1 D C 1 D T n ε i 1 λ + λ exp ( p σ F 0 1 D C 1 D T n ) m 0 1 D C 1 D T n + μ σ j + σ k               ( ε > ε u )
of which
p ( σ ) = α 0 E ε 1 ( σ 1 + 2 σ 3 ) σ 1 2 μ σ 3 + E ε 1 ( σ 1 σ 3 ) 3 σ 1 2 μ σ 3

3. Determination of Constitutive Model Parameters

3.1. Determine the Parameters m and F

In the triaxial and uniaxial tests of rock under chemical-freezing-thawing, the parameters E and μ in the model proposed in this paper can be measured by the rock mechanics test after chemical-freezing-thawing. As can be seen from the above description, the stress has a distinct peak value. The stress-strain curve satisfies the model (32) at the stress peak, we can get
σ ε = ε s = σ s = E ε s 1 λ + λ exp ( p ( σ s ) F ) m + 2 μ σ 3    
The first derivative of the stress-strain curve is zero at the peak, that is
σ 1 ε 1 σ = σ s ε = ε s = 0 = E 1 λ + λ exp ( p ( σ s ) F ) m + ε s exp ( p ( σ s ) F ) m m ( p ( σ s ) F ) m 1 p ( σ s ) ε 1
of which
p ( σ s ) ε 1 = E 3 α 0 ( σ s + 2 σ 3 ) + ( σ s σ 3 ) 3 ( σ s 2 μ σ 3 )
According to Equations (34) and (35), we can get
m = σ s 2 μ σ 3 σ s 2 μ σ 3 + ( λ 1 ) E ε s ln σ s 2 μ σ 3 + ( λ 1 ) E ε s λ E ε s
F = p ( σ s ) ln λ E ε s σ s 2 μ σ 3 + ( λ 1 ) E ε s 1 / m
of which
p ( σ s ) = α 0 E ε s ( σ s + 2 σ 3 ) σ s 2 μ σ 3 + E ε s ( σ s σ 3 ) 3 ( σ s 2 μ σ 3 )

3.2. Determine the Parameter α

The deformation and failure generated in the rock mechanics test is a process of continuous change. The stress-strain curve has the same change law before the elastic deformation to the stable development stage of micro-elastic cracks and the unstable rupture development stage, which is continuous at point C ( ε u , σ u ) , we can get
λ E ε u exp ( p ( σ u ) F ) m           + 2   μ σ 3 = α σ u ( ε ε u ) α
of which
p ( σ u ) = α 0 E ε 1 ( σ u + 2 σ 3 ) σ u 2 μ σ 3 + E ε 1 ( σ u σ 3 ) 3 ( σ u 2 μ σ 3 )
σ 1 ε 1 σ = σ u ε = ε u = λ E exp ( p ( σ u ) F ) m 1 m F ( p ( σ u ) F ) m   = α σ u ε u ( ε ε u ) α - 1
The parameter α can be determined by combining Equations (40) to (42). At the same time, the parameter α can also be obtained by fitting the exponential function curve and it can be obtained by substituting the above formula for verification.

3.3. Determine the Parameter λ

The value of the post-peak correction coefficient λ is the key to ensure that the rock statistical damage constitutive model can accurately describe the post-peak shape of the rock stress-strain curve. For the post-peak stress path of the rock during the loading process, the smaller the mean square error of the stress calculated by the damage constitutive model, the better the practicability of the parameter λ for the rock. The stress mean square error calculated by the established damage model can be expressed as
M S E = 1 m i = 1 m ( σ 1 M σ 1 E ) 2
where σ 1 M is the axial stress obtained by substituting the parameter into Equation (32), and σ 1 E is the axial stress obtained from the test.
Therefore, when the mean square error of the axial stress is the smallest, it can be considered that the post-peak stress-strain curve of the rock obtained by the established damage model is the most accurate, and the corresponding λ at this time is the optimal solution.

4. Model Verification and Analysis

In order to verify the rationality of the statistical damage model of rock under the action of chemical-freezing-thawing-loading, the uniaxial compression test results of granite with different chemical solutions and different freeze-thaw cycles in the literature [21] were selected for comparative analysis. The freeze-thaw equipment in the literature [21] is the CABR-HDK9A rapid freeze-thaw test machine researched and produced by Jianyan Building Materials Co., Ltd. of China Academy of Building Research, Beijing, China and its freeze-thaw temperature is −20 °C~20 °C. The uniaxial compression test adopts the microcomputer-controlled rigid servo triaxial compression testing machine produced by Xi’an Lichuang Material Testing Technology Co., Ltd., Xi’an China, and its maximum load is 2000 kN.
The samples were soaked in H2O, 1 mol L−1 NaOH solution and 1 mol L−1 HNO3 solution with pH ≈ 7.0 for 90 days. Then the sample was put into the rapid freeze-thaw test machine, add H2O with pH ≈ 7.0, 1 mol L−1 NaOH solution and 1 mol L−1 HNO3 solution, respectively, and freeze-thaw cycles at −15 °C~20 °C. The cycle is about 4 h. The uniaxial compression test was carried out on the frozen thawed granite by using a microcomputer-controlled rigid servo triaxial compression tester. The test adopted a stress control method and applied an axial load at a rate of 0.8 MPa/s until the sample failed. The samples were divided into three grades of freeze-thaw cycles: 0, 75, and 100 cycles, with twelve samples in each grade, including four samples in H2O, NaOH solution and HNO3 solution, as shown in Table 1. The obtained uniaxial compression test data are shown in Table 2.
According to the data in Table 2, the variation laws of peak stress, peak strain, and elastic modulus of different freeze-thaw times under the same chemical solutions and different chemical solutions under the same freeze-thaw times can be obtained, as shown in Figure 4. For the parameters of the constitutive model constructed in this paper, the parameters α, m , F and λ of the constitutive model under different chemical solutions and freeze-thaw cycles can be obtained according to the experimental data of the literature [21] and the above parameter calculation process, as shown in Table 3. From Equation (32), the comparison between the theoretical stress-strain curve of granite under different chemical solutions and freeze-thaw cycles and the experimental curve of the literature [21] can be obtained, as shown in Figure 5.
It can be seen from Figure 4 and Figure 5 that the typical total stress-strain curve change law of granite can be roughly divided into four stages as mentioned above: the pore and fracture compaction stage, elastic deformation to the stable development stage of micro-elastic cracks, the unstable rupture development stage, and the post-rupture stage. The compaction stage is mainly caused by the closure of micro-cracks and pores in the rock under uniaxial loading, and the compaction phenomenon is gradually obvious with the increase of the number of freeze-thaw cycles. At the post-rupture stage, with the increase of the number of freeze-thaw cycles, the stress drop phenomenon weakened and the plasticity gradually increased.
After 90 days of immersion in chemical solution without freeze-thaw, the peak stress of granite soaked in H2O was basically unchanged, and the peak stress of granite soaked in HNO3 solution was slightly reduced. However, the mineral components in granite soaked in NaOH solution react with NaOH solution to repair the initial internal micro-cracks and pores but increase the peak stress. The elastic modulus of granite soaked in H2O and NaOH solution did not change significantly, but the elastic modulus of granite soaked in HNO3 solution decreased. To sum up, HNO3 solution aggravated the damage and deterioration of granite, while NaOH solution has a certain inhibitory effect on rock damage. Under the same number of freeze-thaw cycles, the stress-strain curves of granite in H2O and NaOH solutions were close to each other, while the stress-strain curves of rocks soaked in HNO3 solution changed greatly, mainly manifested as inconsistent shapes and uneven curves. Furthermore, its peak strain is quite different from the stress-strain curves in H2O and NaOH solutions. The plasticity of the rock soaked in HNO3 solution was enhanced, the brittleness was weakened, and it was more easily damaged than that soaked in H2O and NaOH solution.
Under the same solution immersion, the elastic modulus, peak stress, and threshold stress of granite in general decrease with the increase of freeze-thaw cycles, while the peak strain and threshold strain increase with the increase of freeze-thaw cycles. It was gradually obvious with the increase of freeze-thaw cycles, indicating the damage and deterioration effect of freeze-thaw cycles on granite. With the increase of the number of freeze-thaw cycles, the stress-strain curve of the rock soaked in HNO3 solution changed significantly, and the plastic deformation was obvious after 75 freeze-thaw cycles; after 100 freeze-thaw cycles, the NaOH solution showed a certain plastic deformation phenomenon. It shows that in HNO3 solution, the effect of freeze-thaw cycle on the mechanical properties of granite is more obvious, while the effect of freeze-thaw cycle in H2O and NaOH solution is relatively weak.
The theoretical results obtained in this paper have the same general trend as the experimental results in the literature [21]. In particular, the granite stress-stress curve has a good degree of fitting in the compaction stage and the post-destruction stage and can reflect the variation law of stress with strain. It shows that the statistical damage model of rock considering the pore compaction and post-peak shape under the action of chemical-freezing-thawing-loading can reflect the uniaxial compression deformation characteristics of rocks under different chemical solution soaking and freeze-thaw cycles. Additionally, it can well characterize the whole process curve of rock stress-strain under different chemical solution immersion and freeze-thaw cycles, however, after the point (at which damage is generated), there is also a certain deviation between the model curve in this paper and the test data, and the deviation is reasonable. The range is due to the different failure criteria selected during the establishment of the equation, which results in the different setting parameters, and the experimental measurement errors.
According to Equation (16), the change law of damage variables under different chemical solutions (H2O solution, HNO3 solution, NaOH solution) and freeze-thaw times (0 times, 75 times, 100 times) can be obtained, as shown in Figure 6. When the rock is not loaded after chemical-freezing-thawing, some initial damage has already occurred. Additionally, under the same conditions of freeze-thaw times, the value of damage variable D is HNO3 solution > H2O solution > NaOH solution. It shows that the HNO3 solution is more likely to cause damage to the granite, followed by the H2O solution, and the NaOH solution has the least damage. Under the same chemical solution erosion conditions, the value of damage variable D is freeze-thaw 100 times > freeze-thaw 75 times > freeze-thaw 0 times, and shows that freeze-thaw cycles exacerbate rock damage. It can be seen from Figure 6 that the damage variables of the rock in the initial stage, namely the compaction stage, elastic deformation to the stable development stage of micro-elastic cracks, are all zero. When the load reaches the yield point, the damage variable starts to accumulate and rapidly increases to one, and the damage variable remains stable at one, that is, the rock is completely damaged.
The deformation and failure laws of siltstone in the literature [25] and sandstone in the literature [41] are similar to the uniaxial compression test characteristics of granite with different chemical solutions and freeze-thaw cycles in the literature [21], and they are all in line with the damage constitutive model established in this paper. The space is limited here and they will not be repeated here. Only the uniaxial compression test of granite with different chemical solutions and different freeze-thaw cycles in the literature [21] is used as a representative for a detailed comparison.

5. Conclusions

In view of the necessity for rock damage research under chemical-freezing-thawing-loading and the defects existing in the traditional rock constitutive model, this paper uses statistical damage theory, assuming that the rock micro units’ strength under chemical-freezing-thawing action obeys the Weibull distribution. It introduces the compaction index and a post-peak correction coefficient that were used to establish a statistical damage model for rock under the action of chemical-freezing-thawing-loading. Analyzing the damage evolution process combined with the rock uniaxial compression test curve, the main conclusions obtained were as follows:
(1)
The compaction index α is introduced to characterize the compaction performance of the internal pores in the rock compaction stage, and fully reflect the nonlinear curve shape of the stress-strain in the undamaged stage. As the compaction index decreases, the compaction effect of the pores in the rock gradually weakens.
(2)
Taking the Drucker–Prager strength criterion as the damage basis of rock micro units, considering the damage threshold and the post-peak correction coefficient, a statistical damage constitutive model of rock under chemical-freezing-thawing-loading based on the Weibull distribution is established, and the determination of each parameter of the model was clarified.
(3)
There is a good correlation between the experimental curve and the theoretical curve calculated by the constitutive model, which can better describe the uniaxial stress-strain relationship under different chemical solutions and freeze-thaw cycles. With the increase in number of freeze-thaw cycles, the rock damage deteriorated seriously, the peak stress, elastic model and threshold stress became smaller, and the increase of plasticity led to the increase of peak strain and threshold strain. Under the same number of freeze-thaw cycles, the HNO3 solution is more likely to cause damage to the granite, followed by the H2O solution, and the NaOH solution has the least damage. From the damage evolution curve, it can be seen that some initial damage has occurred in the rock after chemical-freezing-thawing. During the uniaxial compression process, the damage variable remains unchanged at first, and then begin to change and rapidly increase to one after loading to the yield point. It shows that the model is accurate and reasonable and can provide a theoretical basis for the study of deformation and failure under multi-field coupling conditions.
(4)
The Hoek–Brown strength criterion can be used for further research on the method of measuring the strength of rock micro units, and the model results of the Drucker–Prager strength criterion can be compared and analyzed, hoping to obtain some new discoveries.

Author Contributions

Conceptualization, B.H. and Z.Z.; methodology, Z.Z.; validation, Z.Z., J.L. and H.X.; formal analysis, Z.Z. and H.X.; data curation, K.C.; writing—original draft preparation, Z.Z.; writing—review and editing, B.H. and J.L.; project administration, B.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under grant number U1802243 and grant number 41672317, in part by the Hubei Province Technical Innovation Special (major projects) Project under grant number 2017ACA184, and in part by the Major Science and Technology Projects of WUST Cultivate Innovation Teams under grant number 2018TDX01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, Z.; He, X.; Fan, J.; Zhou, C. Study on the Softening Mechanism and Control of Red-Bed Soft Rock under Seawater Conditions. J. Mar. Sci. Eng. 2019, 7, 235. [Google Scholar] [CrossRef] [Green Version]
  2. Gao, W.; Ge, M. Fracture of rock affected by chemical erosion environment. Geomech. Eng. 2016, 11, 373–383. [Google Scholar] [CrossRef]
  3. Fang, W.; Jiang, N.; Luo, X. Establishment of damage statistical constitutive model of loaded rock and method for determining its parameters under freeze-thaw condition. Cold Reg. Sci. Technol. 2019, 160, 31–38. [Google Scholar] [CrossRef]
  4. Yuan, W.; Chen, C. Experimental study on mechanical characteristics of freeze-thaw damage of fractured rock mass and its failure mechanism. Fresen. Environ. Bull. 2021, 30, 11388–11394. [Google Scholar]
  5. Abdolghanizadeh, K.; Hosseini, M.; Saghafiyazdi, M. Effect of freezing temperature and number of freeze-thaw cycles on mode I and mode II fracture toughness of sandstone. Theor. Appl. Fract. Mech. 2020, 105, 102428. [Google Scholar] [CrossRef]
  6. Hutchinson, A.; Johnson, J.; Thompson, G.; Wood, G.; Sage, P.; Cooke, M. Stone degradation due to wet deposition of pollutants. Corros. Sci. 1993, 34, 1881–1898. [Google Scholar] [CrossRef]
  7. Cai, Y.; Yu, J.; Fu, G.; Li, H. Experimental investigation on the relevance of mechanical properties and porosity of sandstone after hydrochemical erosion. J. Mt. Sci.-Engl. 2016, 13, 2053–2068. [Google Scholar] [CrossRef]
  8. Li, N.; Zhu, Y.; Zhang, P.; Ge, X. A chemical damage model of sandstone in acid environment. Chin. J. Geotech. Eng. 2003, 25, 395–399. [Google Scholar]
  9. Gong, C.; Wang, W.; Shao, J.; Wang, R.; Feng, X. Quantifying microstructural damage of sandstone after hydrochemical corrosion. Acta Geotech. 2020, 16, 1083–1099. [Google Scholar] [CrossRef]
  10. Han, T.; Shi, J.; Chen, Y.; Cao, X. Quantifying microstructural damage of sandstone after hydrochemical corrosion. Int. J. Geomech. 2018, 18, 04018121. [Google Scholar] [CrossRef]
  11. Yang, G.; Zhang, Q.; Pu, Y. Preliminary study on meso-damage propagation characteristics of rock under condition of freezing temperature. Rock Soil Mech. 2004, 25, 1409–1412. [Google Scholar]
  12. Zhou, K.; Li, B.; Li, J.; Deng, H.; Bin, F. Microscopic damage and dynamic mechanical properties of rock under freeze-thaw environment. Trans. Nonferrous Met. Soc. China 2015, 25, 1254–1261. [Google Scholar] [CrossRef]
  13. Yahaghi, J.; Liu, H.; Li, J.; Chan, A.; Fukuda, D. Experimental and numerical studies on failure behaviours of sandstones subject to freeze-thaw cycles. Transp. Geotech. 2021, 31, 100655. [Google Scholar] [CrossRef]
  14. Li, J.; Zhou, K.; Liu, W.; Zhang, Y. Analysis of the effect of freeze-thaw cycles on the degradation of mechanical parameters and slope stability. Bull. Eng. Geol. Environ. 2018, 77, 573–580. [Google Scholar] [CrossRef]
  15. Chen, S. Meso-Mechanism and Its Constitutive Model of Rock Damage under Chemical Erosion. Ph.D. Dissertation, Northeastern University, Harbin, China, 2003. [Google Scholar]
  16. Chen, B.; Feng, X.; Ding, W.; Yang, C. Evolutionary Neural network constitutive model for complete stress-strain relationship of rock under chemical corrosion. J. Northeast Univ. Nat. Sci. 2004, 25, 695–698. [Google Scholar]
  17. Lin, Y.; Gao, F.; Zhou, K.; Gao, R.; Guo, H. Mechanical properties and statistical damage constitutive model of rock under a coupled chemical-mechanical condition. Geofluids 2019, 2019, 7349584. [Google Scholar] [CrossRef]
  18. Huang, S.; Liu, Q.; Cheng, A.; Liu, Y. A statistical damage constitutive model under freeze-thaw and loading for rock and its engineering application. Cold Reg. Sci. Technol. 2018, 145, 142–145. [Google Scholar] [CrossRef]
  19. Feng, Q.; Jin, J.; Zhang, S.; Li, W.; Yang, X.; Li, W. Study on a damage model and uniaxial compression simulation method of frozen-thawed rock. Rock Mech. Rock Eng. 2021, 55, 187–211. [Google Scholar] [CrossRef]
  20. Zhang, H.; Yang, G. Research on damage model of rock under coupling action of freeze-thaw and load. Chin. J. Rock Mech. Eng. 2010, 29, 471–476. [Google Scholar]
  21. Chen, Y.; Wang, P.; Zhang, X.; Du, X. Experimental research on mechanical properties of granite in chemical dissolution under freeze-thaw cycles. Chin. J. Geotech. Eng. 2014, 36, 2226–2235. [Google Scholar]
  22. Chen, Y.; Xiao, P.; Du, X.; Wang, S.; Wang, Z.; Azzam, R. Study on damage statistical constitutive model of triaxial compression of acid-etched rock under coupling effect of temperature and confining pressure. Materials 2021, 14, 7414. [Google Scholar] [CrossRef] [PubMed]
  23. Liao, H.; Chen, Y.; Li, S.; Yin, Y.; Peng, C.; Li, Z. Mechanical properties of sandstone under freeze-thaw cycling and chemical corrosion. J. Disaster Prev. Mitig. Eng. 2020, 40, 1009–1017. [Google Scholar]
  24. Li, X.; Chen, Y.; Xu, P.; Zhong, M.; Wang, Q. Experiment on mechanical property of triaxial compression of white sandstone after effect of chemical solution and freeze-thaw cycle. J. Water Resour. Water Eng. 2015, 26, 212–218. [Google Scholar]
  25. Wang, L.; Wang, X.; Han, J.; Wang, J. Strength degradation and forecast model of siltstone under combined effect of chemical solutions and freeze-thaw cycles. J. Cent. South Univ. Sci. Technol. 2020, 51, 2361–2372. [Google Scholar]
  26. Li, J.; Kaunda, R.; Zhu, L.; Zhou, K.; Gao, F. Experimental study of the pore structure deterioration of sandstones under freeze-thaw cycles and chemical erosion. Adv. Civ. Eng. 2019, 2019, 9687843. [Google Scholar] [CrossRef]
  27. Li, X.; Qu, D.; Luo, Y.; Ma, R.; Xu, K.; Wang, G. Damage evolution model of sandstone under coupled chemical solution and freeze-thaw process. Cold Reg. Sci. Technol. 2019, 162, 88–95. [Google Scholar] [CrossRef]
  28. Cao, W.; Zhang, S.; Zhao, M. Study on a statistical damage constitutive model with conversion between softening and hardening properties of rock. Eng. Mech. 2006, 23, 110–115. [Google Scholar]
  29. Zhang, M.; Li, Z.; Su, X. Probabilistic volume element modeling in elastic damage analysis of quasi-brittle materials. Chin. J. Rock Mech. Eng. 2005, 24, 4282–4288. [Google Scholar]
  30. Zhang, C.; Yang, Q.; Cao, W. Study of damage constitutive model of brittle rock considering post-peak stress dropping rate. Rock Soil Mech. 2019, 40, 3099–3106. [Google Scholar]
  31. Bilen, C.; Er, S.; Tugrul, A.; Yilmaz, M. A new approach for the prediction of brittleness index based on chemical properties of basaltic rocks. Acta Geodyn. Geomater. 2019, 18, 285–299. [Google Scholar] [CrossRef]
  32. Tang, L.; Wang, S. Analysis on mechanism and quantitative methods of chemical damage in water-rock interaction. Chin. J. Rock Mech. Eng. 2002, 21, 314–319. [Google Scholar]
  33. Zhang, Q.; Yang, G.; Ren, J. New study of damage variable and constitutive equation of rock. Chin. J. Rock Mech. Eng. 2003, 22, 30–34. [Google Scholar]
  34. Chen, X.; He, P.; Qin, Z.; Li, J.; Gong, Y. Statistical damage model of altered granite under dry-wet cycles. Symmetry 2019, 11, 41. [Google Scholar] [CrossRef] [Green Version]
  35. Sanchidrián, J.A.; Ouchterlony, F.; Segarra, P.; Moser, P. Size distribution functions for rock fragments–ScienceDirect. Int. J. Rock Mech. Min. 2014, 71, 381–394. [Google Scholar] [CrossRef]
  36. Chen, S.; Qiao, C.; Ye, Q.; Khan, M. Comparative study on three-dimensional statistical damage constitutive modified model of rock based on power function and Weibull distribution. Environ. Earth Sci. 2018, 77, 108. [Google Scholar] [CrossRef]
  37. Xiang, L.; Cao, W.; Su, Y. A statistical damage constitutive model for softening behavior of rocks. Eng. Geol. 2012, 143–144, 1–17. [Google Scholar]
  38. Lemaitre, J. A continuous damage mechanics model for ductile fracture. J. Eng. Mater.-Technol. ASME 2005, 24, 4282–4288. [Google Scholar] [CrossRef]
  39. Liang, M.; Miao, S.; Cai, M.; Huang, Z.; Yang, J. A damage constitutive model of rock with consideration of dilatation and post-peakshape of the stress-strain curve. Chin. J. Rock Mech. Eng. 2021, 40, 2392–2401. [Google Scholar]
  40. Cai, M.; He, M.; Liu, D. Rock Mechanics and Engineering, 2nd ed.; Science Press: Beijing, China, 2013; pp. 54–57. [Google Scholar]
  41. Han, T.; Shi, J.; Chen, Y. Laboratory investigations on the mechanical properties degradation of sandstone under the combined action between water chemical corrosion and freezing and thawing cycles. J. Hydraul. Eng. 2016, 47, 644–655. [Google Scholar]
Figure 1. Schematic diagram of rock damage state.
Figure 1. Schematic diagram of rock damage state.
Jmse 10 00696 g001
Figure 2. Total stress-strain curve of rock deformation.
Figure 2. Total stress-strain curve of rock deformation.
Jmse 10 00696 g002
Figure 3. Stress-strain relationship of rock with different values of α.
Figure 3. Stress-strain relationship of rock with different values of α.
Jmse 10 00696 g003
Figure 4. Changes of rock parameters under different chemical solutions and freeze-thaw cycles: (a) Elastic modulus; (b) Peak stress; (c) Peak strain; (d) Threshold stress; (e) Threshold strain.
Figure 4. Changes of rock parameters under different chemical solutions and freeze-thaw cycles: (a) Elastic modulus; (b) Peak stress; (c) Peak strain; (d) Threshold stress; (e) Threshold strain.
Jmse 10 00696 g004aJmse 10 00696 g004b
Figure 5. Comparison of the stress-strain curves and the curves of literature [21] under different chemical solutions and freeze-thaw cycles: (a) Freeze-thaw 0 times; (b) Freeze-thaw 75 times; (c) Freeze-thaw 100 times.
Figure 5. Comparison of the stress-strain curves and the curves of literature [21] under different chemical solutions and freeze-thaw cycles: (a) Freeze-thaw 0 times; (b) Freeze-thaw 75 times; (c) Freeze-thaw 100 times.
Jmse 10 00696 g005
Figure 6. Damage evolution curves of uniaxial compression under different chemical solutions and freeze-thaw cycles: (a) Freeze-thaw 0 times; (b) Freeze-thaw 75 times; (c) Freeze-thaw 100 times.
Figure 6. Damage evolution curves of uniaxial compression under different chemical solutions and freeze-thaw cycles: (a) Freeze-thaw 0 times; (b) Freeze-thaw 75 times; (c) Freeze-thaw 100 times.
Jmse 10 00696 g006
Table 1. Uniaxial compression test scheme of granite with different chemical solutions and freeze-thaw cycles.
Table 1. Uniaxial compression test scheme of granite with different chemical solutions and freeze-thaw cycles.
Solution
Category
Soaking Time/dFreeze-Thaw Range/°CFreeze-Thaw Time/hFreeze-Thaw Times
n
Loading Rate/MPa/sNumber of Samples
H2O90−15~20400.84
HNO34
NaOH4
H2O90−15~204750.84
HNO34
NaOH4
H2O90−15~2041000.84
HNO34
NaOH4
Table 2. Uniaxialtest data for different chemical solutions and freeze-thaw cycles.
Table 2. Uniaxialtest data for different chemical solutions and freeze-thaw cycles.
Freeze-Thaw Times
n
Solution
Category
Elastic
Modulus
E/GPa
Peak Stress
σs/MPa
Peak Strain
εs/10−3
Threshold Stress
σu/MPa
Threshold Strain
εu/10−3
0H2O16.25105.8614.6662.6510.59
HNO313.58108.3113.3371.499.82
NaOH16.18115.2914.8476.0111.69
75H2O13.7990.5214.6960.0311.96
HNO311.4583.7113.8757.389.91
NaOH15.8594.9813.6062.4111.00
100H2O12.9884.7415.1456.3012.81
HNO310.2281.4818.0850.9613.42
NaOH14.9093.3315.1062.5212.88
Table 3. Parameter values of constitutive model under different chemical solutions and freeze-thaw cycles.
Table 3. Parameter values of constitutive model under different chemical solutions and freeze-thaw cycles.
Freeze-Thaw CyclesSolution
Category
αmFλ
0H2O2.5115.1541.570.95
HNO31.839.9538.400.90
NaOH2.8214.5032.390.96
75H2O3.851.16126.291.81
HNO32.921.0695.681.49
NaOH3.014.6842.500.72
100H2O2.712.5935.231.02
HNO32.982.5778.992.01
NaOH3.092.5438.630.98
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hu, B.; Zhang, Z.; Li, J.; Xiao, H.; Cui, K. Statistical Damage Model of Rock Based on Compaction Stage and Post-Peak Shape under Chemical-Freezing-Thawing-Loading. J. Mar. Sci. Eng. 2022, 10, 696. https://0-doi-org.brum.beds.ac.uk/10.3390/jmse10050696

AMA Style

Hu B, Zhang Z, Li J, Xiao H, Cui K. Statistical Damage Model of Rock Based on Compaction Stage and Post-Peak Shape under Chemical-Freezing-Thawing-Loading. Journal of Marine Science and Engineering. 2022; 10(5):696. https://0-doi-org.brum.beds.ac.uk/10.3390/jmse10050696

Chicago/Turabian Style

Hu, Bin, Zhen Zhang, Jing Li, Huiping Xiao, and Kai Cui. 2022. "Statistical Damage Model of Rock Based on Compaction Stage and Post-Peak Shape under Chemical-Freezing-Thawing-Loading" Journal of Marine Science and Engineering 10, no. 5: 696. https://0-doi-org.brum.beds.ac.uk/10.3390/jmse10050696

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop