Advances in Materials Science and Engineering
Volume 2021 (2021), Article ID 6639428, 16 pages
https://doi.org/10.1155/2021/6639428

Research Article

Mechanical Modeling of Frozen Coarse-Grained Materials Incorporating Microscale Investigation

Shan-Zhen Li,1,2 Liang Tang,1,2 Shuang Tian,1,2 Xian-Zhang Ling,1,2 Yang-Sheng Ye,³ and De-Gou Cai³

¹School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China
²Heilongjiang Research Center for Rail Transit Engineering in Cold Regions, Harbin 150090, China
³State Key Laboratory for Track Technology of High-Speed Railway, China Academy of Railway Sciences Corporation Limited, Beijing 100081, China

Correspondence should be addressed to Liang Tang; hit_tl@163.com

Received 9 November 2020; Revised 17 February 2021; Accepted 3 May 2021; Published 13 May 2021

Academic Editor: Carlos R.  Rambo

Copyright © 2021 Shan-Zhen Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In the cold regions of China, coarse-grained materials are frequently encountered or used as backfilling materials in infrastructure construction, such as dams, highways, railways, and mineral engineering structures. Effects of confining pressure (0.2, 0.5, and 1 MPa) and frozen temperature (−2, −5, −10, and −15°C) on the stress-strain response and elastic modulus were investigated using triaxial compression tests. Moreover, the microscale structures of a coarse-grained material were obtained by X-ray computed tomography. The coarse-grained material specimens exhibited strain-softening and significant dilatancy behaviors during shearing. A modified model considering microstructures of the material was proposed to describe these phenomena. The predicted values coincided well with the experimental results obtained in this study and other literatures. The sensitivity analysis of parameters indicated that the model can simulate the initial hardening and post-peak strain-softening behavior of soils. And the transition of volume strain from contraction to dilatancy can also be described using this model. The results obtained in this study can provide a helpful reference for the analysis of frozen coarse-grained materials in geotechnical engineering.

1. Introduction

Civil infrastructures in cold regions, such as dams, highways, railways [1], and mineral engineering structures, have been constructed more and more in recent years. Coarse-grained materials are often encountered or used as backfill. It is important to understand the mechanical features of frozen coarse-grained materials when they are used in cold regions.

Frozen soil is a special kind of geotechnical material, which is composed of ice inclusions, mineral particles, gaseous inclusions, and liquid water [2]. Compared with the unfrozen soils, the mechanical properties of frozen soils are much more complicated. Many research studies [3, 4] have been conducted to investigate the deformation behaviors and strength of frozen soil. Effects of ice content [5], strain rate [6], temperature [7], and confining pressure [8] have been studied on the mechanical features of frozen soils. Generally, those studies mainly focused on clay, sand, and other fine-grained soils. The mechanical properties of frozen soil are significantly influenced by soil types [2]. However, almost none of these works properly explores deformation behaviors and strength characteristics of frozen coarse-grained materials. The possible mechanisms for the mechanical properties of frozen coarse-grained materials are interesting to be clarified.

Many constitutive models (e.g., hyperelastic model, nonlinear elastic model, and viscoelastic-plastic model) have been established to describe the mechanical features of frozen soils [9–14]. In addition, on the basis of continuous damage theory, some constitutive models of frozen soil have also been established. The statistical-damage-based method has been applied to the constitutive response of soils successfully in the past 30 years. It is an important approach to investigate the deformation process and failure mechanisms of geotechnical materials. Until now, in the framework of damage mechanics containing statistical changes, numerous constitutive models for ductile or brittle materials have been established [15–17]. The stress-strain response of geotechnical materials can be predicted by the statistical-damage-based model simply and correctly. However, few models have been established to predict the volume change [18].

The purpose of this paper is to propose a model that can describe the mechanical behavior of frozen coarse-grained material. The effects of confining pressure and temperature on the frozen coarse-grained materials were studied with triaxial tests and X-ray computed tomography (CT). Then, an empirical formula that reflects the relationship between microstructural change and residual strength is proposed. Furthermore, the stress-strain response of frozen coarse-grained material considering microstructures was predicted. At last, the parameter sensitivity of the statistical damage model was discussed.

2. Materials and Methods

2.1. Material Properties and Specimen Preparation

The typical grain size distribution of the tested coarse-grained material is shown in Figure 1. The curve coefficient (Cc) and uniformity coefficient (Cu) are 6.48 and 13.78, respectively. According to the Standard for Soil Test Method [19], the tested coarse-grained material is poorly graded gravel soil. The optimum water content and maximum dry density were determined to be 4.5% and 2.02 g/cm³, respectively. And the detailed characteristics of the tested coarse-grained material are listed in Table 1.

Figure 1: Grain size distribution of the coarse-grained material.

Table 1: Basic physical characteristics of the frozen coarse-grained material.

Parameters
Values

Maximum diameter, Dmax (mm)
10

Diameter at 10% finer, D₁₀ (mm)
0.508

Diameter at 30% finer, D₃₀ (mm)
4.8

Diameter at 60% finer, D₆₀ (mm)
7

Uniformity coefficient, Cu
13.78

Curvature coefficient, Cc
6.48

Maximum dry density, ρd,max (g/cm³)
2.02

Optimum moisture content, opt (%)
4.5

An appropriate amount of distilled water was mixed fully with the dried coarse-grained material to obtain the optimum moisture content (Table 1). A closed container was then used to store the mixture for 12 h to ensure a uniform distribution of moisture. The cylindrical specimens with the maximum dry density of 95% were made by a cylindrical mold. The size of specimens was unified into 125 mm in height and 61.8 mm in diameter on a specially made specimen machine (Figure 2) to decrease the difference in the preparation of specimens. The machine has the advantage of easy to use with a wide operating range. Moreover, it is suitable for various sizes of specimens, and the test precision is high. In order to avoid frost heaving, the specimens fastened by three-piece split molds were placed in an incubator at −30°C for 24 h. Then, the completely frozen specimens were removed from the molds and wrapped with a rubber membrane. After that, they were placed in a refrigerator at the objective temperature for more than 12 h to make the temperature of specimens uniform.

Figure 2: The picture of the specimen machine.

2.2. Testing Apparatus and Scheme

The cryogenic triaxial compression tests were performed on a material test system (MTS-810) exhibited in Figure 3. It is mainly composed of two subsystems: the temperature controlling system and the loading system. In the temperature controlling system, the temperatures at the side and top of the pressure chamber are adjusted by two independent refrigeration circulators (shown in Figure 3). In the loading system, the high confining pressure is withstood by a designed steel airtight pressure cell. The axial load measured by a force sensor is applied to the top of the specimens. The parameters of the triaxial testing apparatus are as follows: the maximum axial load is 100 kN, the maximum axial deformation is 85 mm, the stable temperature ranges from –30°C to 30°C with the precision of ±0.1°C, and the experimental confining pressure is 0.1∼20 MPa. To investigate the effects of temperature and confining pressure on the mechanical features of frozen coarse-grained materials, 12 triaxial tests were conducted. The temperatures ranged from −2°C to −15°C, while the confining pressures varied from 0.2 MPa to 1.0 MPa. And the detailed schemes are listed in Table 2.

(a)

(b)

(c)

Figure 3: Triaxial equipment and tested samples: (a) the triaxial test device; (b) the CT scanner; (c) the specimen before and after testing.

Table 2: Experimental schemes of the triaxial tests.

Samples
ω (%)
ρd (g/cm³)
T (°C)
σ₃ (MPa)
CT scanning

S1
4.5
2.02
−2
0.2
—

S2
0.5
Sections I–V

S3
1
—

S4
−5
0.2
Sections I–V

S5
0.5
Sections I–V

S6
1
Sections I–V

S7
−10
0.2
—

S8
0.5
Sections I–V

S9
1
—

S10
−15
0.2
—

S11
0.5
Sections I–V

S12
1
—

Note. “—” denotes no test.

2.3. Experimental Procedures

The undrained static triaxial test was conducted according to the following procedures: (1) the target temperature was kept for 10 h to ensure the uniform temperature distribution in the pressure chamber; (2) the frozen specimen was transferred from the refrigerator to the triaxial pressure chamber; (3) the confining pressure (σ₃) was set to the specified value and kept for 10 min; and (4) the axial deviatoric stresses (q) were imposed at a strain rate of 1%/min until the axial strain reached 15%.

In the past few years, CT has been widely applied in geotechnical studies as a tool of describing the microstructure due to the advantage of being nondestructive [20]. X-ray CT observations were conducted on some selected specimens to analyze the microstructure variations of frozen coarse-grained materials. The spatial and the density resolution of the CT scanner (PHILIPS Brilliance 16) are 0.208 mm and 0.3%, respectively. The specimens before and after the triaxial test were scanned from the bottom to top in the CT scanner with a spacing of 3 mm. As displayed in Figure 3, five typical sections of a specimen were selected to display the microstructure changes.

3. Results and Discussion

3.1. Stress-Strain Response

The typical relationships under different confining pressures and temperatures between deviator stress (q), axial strain (ε₁), and volumetric strain (εv) are illustrated in Figure 4. Generally, there are roughly three stages in the stress-strain curves, i.e., the initial linear elastic stage, the plastic deformation stage, and the softening stage. In the initial stage, there is remarkable linear relationship between deviator stress and axial strain. When the axial strain increases to larger than 1%, most stress-strain curves present the plastic behavior, which means that their slopes reduce gradually in this plastic deformation stage. Under relatively low confining pressure (i.e., σ₃ = 0.2 MPa), a strain-softening phenomenon is observed at different temperatures. This may attribute to the fact that more water turns into ice with the decrease of temperature. Additionally, ice becomes more brittle as temperature decreases [2, 21], which results in the brittle failure of frozen coarse-grained material. At the same temperature, the deformation behavior of frozen coarse-grained material changes from strain softening to strain hardening with the increase of confining pressures. This indicates that the ductile behavior of material is enhanced at higher pressure. Similar behavior is obtained for frozen sand [22]. It can be concluded that temperature and confining pressure have a remarkable effect on the deformation characteristics of frozen soils.

(a)

(b)

(c)

(d)

Figure 4: Stress-strain behaviors under different testing temperatures and confining pressures: (a) T = −2°C; (b) T = −5°C; (c) T = −10°C; (d) T = −15°C.

It can be seen from Figure 4 that the temperature and confining pressure both have significant effects on the stress-strain curves of frozen coarse-grained materials. Under low confining pressures, the deviatoric stress rises slowly in the initial stage. However, the deviatoric stress increases rapidly with the increase of the confining pressure. With the decrease of temperature, the strength of frozen coarse-grained materials gets higher under the same confining pressure. Generally, the decrease of temperature will improve the strength frozen soil significantly. The frozen water acts as a bonding agent between the soil particles, increasing the stiffness of soil. In addition, lower temperature can significantly influence the strength of both single particle [23, 24] and soils [6, 25], usually enhancing it [26].

The key to determine the volumetric strain of the specimen is to determine its volume change, which is difficult to measure because there is no seepage for the frozen coarse-grained material. In this study, the volume changes of specimens were measured by the volume change caused by the leakage of hydraulic oil in the pressure chamber. The minus sign indicates an increase in volume. With the increase of the axial strain, the volumes of specimens reduce from a starting state to a phase transition state and then expand gradually into a residual state. As observed in Figure 4, the confining pressure has a significant influence on the volumetric strain curves of frozen coarse-grained material. The increase of confining pressure reduces the dilation of specimens.

3.2. Elastic Modulus and Strength Properties

The characteristics of the soil could be reflected by two important parameters: elastic modulus and strength properties. The stress-strain curve slope in the elastic deformation stage can be defined as the elastic modulus, according to a previous report [27]. The initial curve tangent with the axial strain of 1% was used to calculate the elastic modulus in this study [16]. As depicted in Figure 5, the increase of confining pressure and the decrease of temperature lead to the monotonic increase of elastic modulus, indicating that both the confining pressure and temperature have significant influences on stiffness. It can be found that when the tested temperature is −2°C, the elastic modulus of 0.2, 0.5, and 1 MPa is approximately 92.1, 115.3, and 163.9 MPa, respectively. When the temperature drops to −15°C, the values of elastic moduli at 0.2, 0.5, and 1 MPa increase by approximately 234%, 198%, and 126%, respectively. The relationships between elastic modulus and these two parameters can be expressed as follows:

where E, |T|, and pa are the elastic modulus, the absolute value of negative temperature, and the atmospheric pressure (pa = 101 kPa), respectively.

Figure 5: Variation in elastic modulus versus frozen temperature and confining pressure.

The variations of peak strength (qf) and residual strength (qr) at various confining pressures and temperatures are illustrated in Figure 6. It is evident that the variations of peak strength and residual strength are similar to those of elastic modulus. Further, their relationships are as follows:

(a)

(b)

Figure 6: Variation in (a) peak strength and (b) residual strength versus frozen temperature and confining pressure.

The strength of the frozen coarse-grained material enhances as the confining pressure increases, which can be interpreted as the following two aspects. Firstly, the rising confining pressure makes the ice particles and soil particles closer, forming a stronger interlocking friction between these particles. Secondly, the rising confining pressure increases the normal pressure on the shear plane, increasing the resistance of sliding friction.

The angle of the internal friction (φ) and the cohesion (c) were the two shear strength parameters used to analyze the effect of temperature on the shear strength of soil. Figure 7(a) illustrates the effect of temperature on the cohesion of specimens. The cohesion of soil rises with the decrease of temperature in the form of a power function, as shown in the following equation:

(a)

(b)

Figure 7: Variation of (a) cohesion and (b) internal friction angle versus frozen temperature.

The variations of internal friction angle with different frozen temperatures are presented in Figure 7(b). It can be seen that the friction angles fluctuate in a small range near the average value of 34.11°. The results indicate that the frozen temperature has little influence on the tangential interaction between soil and ice particles. On the contrary, the binding force between ice and soil particles is highly dependent on the frozen temperature.

3.3. X-Ray CT Observations

In order to describe the microstructure changes of a specimen before and after the triaxial test, the typical CT images of the specimen under 0.5 MPa and −2°C are presented in Figure 8. Before testing, the particles are evenly distributed, and thus the appearance of specimen slices is almost the same. After testing, from top to bottom, the soil particle distributions are from loose to tight and then to looser.

Figure 8: Typical CT images of the frozen coarse-grained specimen of case S2 before and after testing.

Apart from soil particles, the unfrozen water and internal ice are related to the fabric of frozen soil. The distribution of unfrozen water, internal ice, and soil particles on the cross section of the specimen can be characterized by the average density, which is defined by the CT value (AV) [28]:

where and ui represent the X-ray absorption coefficients of water and material, respectively. Water has a CT value of 0 HU by setting.

The relationships between CT values and section positions are shown in Figure 9. The CT values of specimens after testing are smaller than those before testing. Therefore, the deformation of the specimen belongs to the dilation type, which is coincident with its strain-volumetric behavior presented in Figure 4. It is apparent that the CT values distribute as a bow shape. It decreases to the minimum value at the middle section of the specimen. This means that the particle distributions in the middle part of specimens are looser than those in both ends of specimens, which can be interpreted as the end effect [29]. These results are consistent with the shape of the specimen after the triaxial test, as observed in Figure 3. CT values of tested specimens in five typical sections decrease gradually as the temperature decreases, which is consistent with the strain-softening phenomenon caused by the increase of brittleness from ice (Figure 10(a)). Besides, increasing the confining pressure also leads to an increase in CT values as shown in Figure 10(b).

(a)

(b)

Figure 9: Relationship between CT value and normalized height h/H for specimens before and after testing: (a) confining pressure σ₃ = 0.5 MPa; (b) frozen temperature T = −5°C.

(a)

(b)

Figure 10: CT value under (a) different frozen temperatures and (b) confining pressures.

As reported in a previous study, the gradual failure of cemented structure and interlocking structure gradually reduces the mechanical resistance, resulting in post-peak softening response. The residual strength is closely related to the structure of the shear specimen [30]. The evolution of specimen structure can be reflected by the change of CT values. And the average change ratio of CT values (RAV) can be defined as follows:

where AVbefore and AVafter are the average CT values before and after the test.

Figure 11 shows the relationship between the residual strength and RAV. It can be seen that the residual strength is positively correlated with RAV at −5°C and different confining pressures, which can be expressed as follows:

and the residual strength is negatively correlated with RAV under 0.5 MPa and different temperatures, which is expressed as follows:

Figure 11: Relationship between residual strength and the average change ratio of CT value.

Based on the above analysis, it is speculated that the residual strength of frozen coarse-grained material is linear with RAV, which can be expressed as follows:

where aT is a parameter related to frozen temperature, which is obtained by experimental result.

The relationship between parameter aT and the frozen temperature is presented in Figure 12. Case S8 used for the following validation was not considered in this figure. As depicted in Figure 12, the parameter aT is significantly affected by frozen temperature, which can be expressed as follows:

Figure 12: Relationship between parameter aT and frozen temperature.

Thereafter, a simple empirical formula is proposed by substituting equation (10) into equation (9) to predict the residual strength, which is expressed as follows:

4. Statistical Damage Model Incorporating Microscale Structure

4.1. Prediction Model of Stress-Strain-Volumetric Characteristics

It is assumed that the frozen coarse-grained material consists of microunits. And the material strength of microunits differs. By incorporating the statistical strength theory and the continuous damage theory, the damage variable is described as follows:

where D is the damage variable and Nf and N are the failed and total number of microunits, respectively.

When the Weibull distribution is assumed to be the density function of microdamage, the probability density distribution function f(F) can be expressed as follows:

where Weibull modulus m and F₀ are the material parameters, which correspond to the mechanical features of the frozen coarse-grained material. The variable F corresponds to the random distribution of the microunit strength [31].

The strength criterion [27, 32] or the axial strain [33, 34] of the frozen coarse-grained material can be used to describe the random distribution variable F of the microunit strength. In this study, the axial strain is applied for the expression of random distribution variable of the microunit strength. This form is simple and convenient for application. Therefore, the number of damaged microunits can be determined as follows:

The damage variable based on the Weibull distribution is obtained by substituting equation (14) into equation (12), which is determined as follows:

Then, the stress-strain relationship is described as follows:

where E represents the elastic modulus, ε₁ is the axial strain, and ε₀ is the distribution parameter.

As shown in equation (16), the residual stress cannot be reflected [35]. To describe the failure process of frozen coarse-grained material accurately, a modified model considering residual stress is obtained in this paper. Therefore, the stress-strain relationship can be expressed as follows:

where Ee represents the nominal elastic modulus, which is related to elastic modulus of frozen coarse-grained material.

The prediction model of the stress-strain behavior of the frozen coarse-grained material is obtained by combining equation (15) with equation (17), which is expressed as follows:

All of the parameters are related to the mechanical properties of frozen coarse-grained material. The residual strength qr is determined directly by the stress-strain curves in Figure 4. Other three parameters (Ee, ε₀, and m) can neither be obtained directly from the experimental results, nor be derived by the formulas. Generally, m is determined by the peak point method, which was interpreted by Li et al. in detail [36]. However, this commonly used method cannot reflect the residual stress of frozen coarse-grained material. Actually, m is a parameter of constitutive relation curve, which can be determined by fitting the experimental data according to Huang et al. [37]. Thus, the results of parameters Ee, ε₀, and m in equation (18) are obtained by fitting.

As observed in Figure 13, E is linear with Ee, which can be expressed as follows:

Figure 13: Relationship between E and Ee.

After combining equations (11) and (19) with equation (18), the statistical damage model incorporating the microscale structure is established as follows:

In the binary-medium constitutive model [18], the artificial structural soil is treated as the binary-medium material composed of friction and bonding elements. The volumetric strain is defined by the strains of these two elements. And the damage variable is used to replace the damage rate. Thus, the volume change of frozen coarse-grained material can be expressed as follows:

where and represent the strains of frictional elements and bonded elements, induced by the deviatoric stress and spherical stress, respectively.

and in equation (21) are defined in equations (22) and (23), respectively:

where p₀ and p are the spherical stresses before and during the test, respectively, with the relationship of p = q/3 + p₀; A, B, and C are parameters related to the volume change characteristics of frozen coarse-grained material; and K is the bulk modulus.

After substituting equations (15), (22), and (23) into equation (21), a prediction model of volume change can be obtained as follows:

4.2. Model Verification

Figure 14 shows the comparison of experimental results and the predicted results in this study. According the equations derived in the previous section, the related material parameters are obtained and summarized in Table 3. The predicted values agree well with the experimental results. It can be concluded that the proposed models can predict the stress-strain properties and the volume change characteristics of frozen coarse-grained materials.

(a)

(b)

(c)

(d)

Figure 14: Predicted and experimental deviatoric stress-strain results of frozen coarse-grained materials: (a) T = −2°C; (b) T = −5°C; (c) T = −10°C; (d) T = −15°C.

Table 3: Model parameters of the frozen coarse-grained material in this study.

Samples
Stress-axial strain parameters
Volumetric strain-axial strain parameters

T (°C)
σ₃ (MPa)
qr (MPa)
Ee (MPa)
ε₀
m
R²
K (MPa)
A
B
C
R²

S1
−2
0.2
1.193
76.164
0.019
0.981
0.991
76.667
−0.263
0.693
0.013
0.999

S2
−2
0.5
2.130
96.130
0.025
1.089
0.989
95.833
−0.281
1.059
0.006
0.998

S3
−2
1
3.428
138.049
0.018
1.114
0.988
136.625
−0.408
1.888
0.004
0.943

S4
−5
0.2
1.346
165.000
0.022
0.991
0.944
161.650
−0.518
0.861
0.017
0.999

S5
−5
0.5
2.889
179.791
0.021
0.945
0.994
184.967
−0.419
1.060
0.012
0.997

S6
−5
1
3.969
200.780
0.020
1.126
0.956
209.275
−0.426
1.582
0.009
0.987

S7
−10
0.2
1.440
217.218
0.027
1.135
0.972
192.542
−0.842
1.044
0.016
0.999

S8
−10
0.5
3.219
227.500
0.024
0.981
0.997
226.225
−0.744
1.300
0.013
0.997

S9
−10
1
5.000
251.082
0.022
1.077
0.948
261.542
−0.487
1.388
0.011
0.993

S10
−15
0.2
1.866
255.460
0.031
1.221
0.899
256.208
−0.931
1.055
0.020
0.999

S11
−15
0.5
3.801
273.700
0.020
0.879
0.996
285.308
−0.955
1.399
0.012
0.997

S12
−15
1
5.313
279.103
0.022
0.996
0.983
309.358
−0.858
1.656
0.010
0.991

Two kinds of frozen soils in previous studies [12, 22] are selected to validate the presented models based on the statistical damage concept. The comparison between experimental and the predicted results by equations (18) and (24) is demonstrated in Figure 15. The parameters used to predict the behaviors of two types of soils are listed in Table 4. It can be seen that the proposed model can correctly simulate the stress-strain features of frozen soils. And the plastic deformation mechanism (i.e., compressibility and expansibility) of soils can be accurately described by the proposed models. This result seems to make it possible for the proposed model to capture the nonlinear mechanical behavior of soil.

(a)

(b)

Figure 15: Predicted and experimental results of frozen soils in previous studies: (a) frozen saline sandy soil [12]; (b) frozen standard sand [22].

Table 4: Model parameters of frozen saline sandy soil and frozen standard sand.

Soil types
Stress-axial strain parameters
Volumetric strain-axial strain parameters

T (°C)
σ₃ (MPa)
qr (MPa)
Ee (MPa)
ε₀
m
R²
K (MPa)
A
B
C
R²

Frozen saline sandy soil [12]
−6
2
10.92
274.1
0.056
1.09
0.984
370.839
−1.017
2.05
0.013
0.982

Frozen standard sand [22]
−15
1
4.304
352.3
0.019
0.895
0.987
489.167
−0.914
1.193
0.012
0.985

4.3. Sensitivity Analysis of the Parameters

The parametric investigation results of the prediction model are analyzed in this section. Figures 16 and 17 present the stress-strain and volume change behavior of frozen coarse-grained materials. These results indicate that the main mechanical characteristics of the frozen coarse-grained materials can be duplicated by the predicted model.

(a)

(b)

(c)

Figure 16: Deviatoric stress-strain curves for different values of parameters (a) ee, (b) ε₀, and (c) m.

(a)

(b)

(c)

Figure 17: Volumetric strain curves for different values of parameters (a) A, (b) B, and (c) C.

In terms of the stress-strain responses, the behavior of specimens transform from strain hardening to strain softening with the increase of Ee. And the peak strength increases sharply. With the increase of ε₀, the peak stress and the corresponding strain increase gradually. It seems that the strength and the plastic performance of frozen coarse-grained material have been enhanced. There is a slight transformation from strain hardening to strain softening with the increase of m. For the volumetric strain, with the increase of A and B, the trend of dilation diminishes. Parameter C determines the rate of contraction. It can be clearly observed that the model can predict the transition from initial hardening to strain softening. Moreover, the dilation response of volume strain can also be described by this model.

5. Conclusions

The effects of confining pressure and temperature on the mechanical property of frozen coarse-grained material were investigated using the triaxial compression test and X-ray CT analysis. Based on the statistical damage theory, a modified model considering the microscale structure was developed. This model is able to describe the strain-softening and dilatancy behaviors of frozen coarse-grained materials. The following conclusions can be drawn:

(1)Temperature and confining pressure have significant influence on the mechanical properties of frozen coarse-grained materials. The specimens exhibit strain-softening and dilatancy behaviors during shearing. With the increase of the confining pressure and temperature, the stress-strain relationship transfers from strain softening to strain hardening.(2)CT values of the tested material specimens decrease gradually with the decrease of the temperature. The increase of the confining pressure leads to the increase of CT values, which is consistent with the dilatancy phenomenon. Moreover, an empirical formula that reflects the relationship between microstructural and residual strength is proposed.(3)A model based on statistical damage theory was proposed to predict the stress-strain response and volume change of frozen coarse-grained material. It was verified with the experimental results obtained by this study and other previous studies. The sensitivity analysis of the model parameters indicated that the model can predict the transition from initial hardening to strain softening. Moreover, the dilation response of volume strain can also be described.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by the National Key Research and Development Project of China (grant nos. 2018YFC1505305 and 2018YFE0207100), the National Major Scientific Instruments Development Project of China (grant no. 41627801), the State Key Program of National Natural Science Foundation of China (grant no. 41731288), and the Technology Research and Development Plan Program of Heilongjiang Province (grant no. GA19A501).

References

G. Kweon and T. Hwang, “Deformational characteristics of subgrade soils and subbase materials with freeze-thaw,” KSCE Journal of Civil Engineering, vol. 17, no. 6, pp. 1317–1322, 2013.

M. Kadivar and K. N. Manahiloh, “Revisiting parameters that dictate the mechanical behavior of frozen soils,” Cold Regions Science and Technology, vol. 163, pp. 34–43, 2019.

X. Chang, W. Ma, and D. Wang, “Study on the strength of frozen clay at high confining pressure,” Frontiers of Earth Science in China, vol. 2, no. 2, pp. 240–242, 2008.

Y. Yang, Y. Lai, and X. Chang, “Laboratory and theoretical investigations on the deformation and strength behaviors of artificial frozen soil,” Cold Regions Science and Technology, vol. 64, no. 1, pp. 39–45, 2010.

Y. Yamamoto and S. M. Springman, “Axial compression stress path tests on artificial frozen soil samples in a triaxial device at temperatures just below 0°C,” Canadian Geotechnical Journal, vol. 51, no. 10, pp. 1178–1195, 2014.

J. Lee, Y. S. Kim, D. Chae, and W. Cho, “Loading rate effects on strength and stiffness of frozen sands,” KSCE Journal of Civil Engineering, vol. 20, no. 1, pp. 208–215, 2015.

X. Xu, Y. Wang, Z. Yin, and H. Zhang, “Effect of temperature and strain rate on mechanical characteristics and constitutive model of frozen Helin loess,” Cold Regions Science and Technology, vol. 136, pp. 44–51, 2017.

W. Jinyuan and N. Satoshi, “Interpretation of mechanical behavior of frozen clay through parallel tests of frozen and unfrozen soils,” Japanese Geotechnical Society Special Publication, vol. 5, no. 2, pp. 155–160, 2017.

D. Chen, D. Y. Wang, W. Ma, L. L. Lei, and G. Y. Li, “A strength criterion for frozen clay considering the influence of stress Lode angle,” Canadian Geotechnical Journal, vol. 56, no. 11, pp. 1557–1572, 2018.

S. A. Ghoreishian Amiri, G. Grimstad, M. Kadivar, and S. Nordal, “Constitutive model for rate-independent behavior of saturated frozen soils,” Canadian Geotechnical Journal, vol. 53, no. 10, pp. 1646–1657, 2016.

Y. M. Lai, Y. G. Yang, X. X. Chang, and S. Y. Li, “Strength criterion and elastoplastic constitutive model of frozen silt in generalized plastic mechanics,” International Journal of Plasticity, vol. 26, no. 10, pp. 1461–1484, 2010.

Y. Lai, M. Liao, and K. Hu, “A constitutive model of frozen saline sandy soil based on energy dissipation theory,” International Journal of Plasticity, vol. 78, pp. 84–113, 2016.

S. L. Painter and S. Karra, “Constitutive model for unfrozen water content in subfreezing unsaturated soils,” Vadose Zone Journal, vol. 13, no. 4, pp. 1–8, 2014.

D. Y. Wang, M. Wei, and X. X. Chang, “Analyses of behavior of stress–strain of frozen Lanzhou loess subjected to K0 consolidation,” Cold Regions Science and Technology, vol. 40, no. 1-2, pp. 19–29, 2004.

S. Y. Li, Y. M. Lai, S. J. Zhang, and D. R. Liu, “An improved statistical damage constitutive model for warm frozen clay based on Mohr-Coulomb criterion,” Cold Regions Science and Technology, vol. 57, no. 2-3, pp. 154–159, 2009.

L. Tang, S. Cong, X. Ling, W. Xing, and Z. Nie, “A unified formulation of stress-strain relations considering micro-damage for expansive soils exposed to freeze-thaw cycles,” Cold Regions Science and Technology, vol. 153, pp. 164–171, 2018.

Y. Lai, S. Li, J. Qi, Z. Gao, and X. Chang, “Strength distributions of warm frozen clay and its stochastic damage constitutive model,” Cold Regions Science and Technology, vol. 53, no. 2, pp. 200–215, 2008.

E. L. Liu, H. S. Yu, C. Zhou, Q. Nie, and K. T. Luo, “A binary-medium constitutive model for artificially structured soils based on the disturbed state concept and homogenization theory,” International Journal of Geomechanics, vol. 17, no. 7, Article ID 04016154, 2016.

Ministry of Water Resources of China, Specification of Soil Test (SL237-1999), China Water Power Press, Beijing, China, 1999, in Chinese.

J. L. Qi, Y. M. Lai, and Y. B. Pu, “Experimental study on the micro-fabric of frozen sediment using triaxial deformation and computerized tomography,” Annals of Glaciology, vol. 37, no. 1, pp. 105–107, 2003.

M. A. Rist and S. A. F. Murrell, “Ice triaxial deformation and fracture,” Journal of Glaciology, vol. 40, no. 135, pp. 305–318, 1994.

X. Yao, J. Qi, F. Yu, and L. Ma, “A versatile triaxial apparatus for frozen soils,” Cold Regions Science and Technology, vol. 92, pp. 48–54, 2013.

Y. Xiao, C. S. Desai, A. Daouadji, A. W. Stuedlein, H. Liu, and H. Abuel-Naga, “Grain crushing in geoscience materials-Key issues on crushing response, measurement and modeling: review and preface,” Geoscience Frontiers, vol. 11, no. 2, pp. 363–374, 2020.

M. Meng, Z. Sun, C. Wang, X. He, and Y. Xiao, “Size effect on mudstone strength during freezing-thawing cycle,” Journal of Environmental Geotechnics, vol. 78, no. 15, pp. 1–13, 2020.

S. H. He, H. F. Shan, T. D. Xia, Z. J. Liu, and F. Xia, “The effect of temperature on the drained shear behavior of calcareous sand,” Acta Geotechnica, vol. 16, no. 3, pp. 613–633, 2020.

C. Cai, W. Ma, Z. Zhou et al., “Laboratory investigation on strengthening behavior of frozen China standard sand,” Acta Geotechnica, vol. 39, no. 1-2, pp. 179–192, 2018.

H. Li, H. Liao, G. Xiong, B. Han, and G. Zhao, “A three-dimensional statistical damage constitutive model for geomaterials,” Journal of Mechanical Science and Technology, vol. 29, no. 1, pp. 71–77, 2015.

C. Kawaragi, T. Yoneda, T. Sato, and K. Kaneko, “Microstructure of saturated bentonites characterized by X-ray CT observations,” Engineering Geology, vol. 106, no. 1-2, pp. 51–57, 2009.

R. Pezo, M. Picornell, and S. Nazarian, “Evaluation of strain variation within a triaxial specimen due to end effects,” Transportation Research Record: Journal of the Transportation Research Board, vol. 1614, no. 1, pp. 35–42, 1998.

R. C. K. Wong, “Strength of two structured soils in triaxial compression,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 25, no. 2, pp. 131–153, 2001.

X. Yang, S. Zengchun, S. Chandrakant, and M. Desai, “Strength and surviving probability in grain crushing under acidic erosion and compression,” International Journal of Geomechanics, vol. 19, no. 11, pp. 1–15, 2019.

W.-G. Cao, H. Zhao, X. Li, and Y.-J. Zhang, “Statistical damage model with strain softening and hardening for rocks under the influence of voids and volume changes,” Canadian Geotechnical Journal, vol. 47, no. 8, pp. 857–871, 2010.

C. A. Tang, H. Liu, P. K. K. Lee, Y. Tsui, and L. G. Tham, “Numerical studies of the influence of microstructure on rock failure in uniaxial compression - Part I: effect of heterogeneity,” International Journal of Rock Mechanics and Mining Sciences, vol. 37, no. 4, pp. 555–569, 2000.

H. Liu and L. Zhang, “A damage constitutive model for rock mass with nonpersistently closed joints under uniaxial compression,” Arabian Journal for Science and Engineering, vol. 40, no. 11, pp. 3107–3117, 2015.

J. Wang, Z. Song, B. Zhao, X. Liu, J. Liu, and J. Lai, “A study on the mechanical behavior and statistical damage constitutive model of sandstone,” Arabian Journal for Science and Engineering, vol. 43, no. 10, pp. 5179–5192, 2018.

X. Li, W. G. Cao, and Y. H. Su, “A statistical damage constitutive model for softening behavior of rocks,” Engineering Geology, vol. 143-144, pp. 1–17, 2012.

S. Huang, Q. Liu, A. Cheng, and Y. Liu, “A statistical damage constitutive model under freeze-thaw and loading for rock and its engineering application,” Cold Regions Science and Technology, vol. 145, pp. 142–150, 2017.