Investigation on the Nonlinear Strength Properties and Damage Statistical Constitutive Model for Frozen Sandy Soils
There are many flaws, such as fissures, cavities, and inclusions, in geomaterials, which make their mechanical properties with great randomness and uncertainty. Upon loading, the soil structure gradually losses the bearing capacity due to the transformation from microdefects to macroscopic breakage bands. Based upon the experimental data of frozen sandy soils, a new nonlinear strength equation between the first and third principal stresses was proposed, and then the nonlinear strength properties for frozen sandy soils in σ-τ plane were analyzed. In addition, by assuming that the microstrength of frozen sandy soil obeys the Weibull distribution function, a statistical damage constitutive model was established based upon the framework of continuum damage mechanics (CDM), with few parameters and a high accuracy. Compared with experimental data, the new model can well grasp the nonlinear strength properties and simulate the stress-strain relationships under different confining pressures for frozen sandy soils.
Frozen soils are geotechnical materials, which are mainly formed in cryogenic environment. Frozen soils are defined as those containing some ice and having a temperature at or lower than 0°C . It is distributed all around the world regularly, such as Russia, America, China, and Canada, which accounts for 24% of the world’s land area . As far as the components are concerned, frozen soils are composite materials consisting of mineral particles, ice inclusions, liquid water, and gaseous inclusions , which are sensitive to temperature, loading rate, external load, soil types, water content, confining pressure, etc. The typical problems in cold regions can be mainly summarized as the frost-heaving effect, thawing-settlement effect, and freeze-thaw cycle effect in these regions. It is well known that the frost-heaving effect causes the structure damage due to the extension of microcracks. Simultaneously, the thawing-settlement effect induces irreversible plastic deformation. Terribly, the interactions of the two effects aggravate the damage and destruction of structures. Hence, up to present, many researchers have studied the physical and mechanical characteristics of permafrost and seasonal frozen soils, such as nonlinear strength properties and damage deformation behaviors.
Currently, a great number of research achievements have been obtained on this field for unfrozen soils and frozen soils. For instance, the mechanical properties of unfrozen soils are obtained with great progress. Matsuoka and Nakai  considered the stress deformation and strength characteristics of soil under three different principal stresses. Liu et al. [5, 6] proposed a general strength criterion for geomaterials, considering isotropic and anisotropic effects. Meanwhile, Liu et al.  employed a nonlinear Drucker–Prager and Matsuoka–Nakai unified failure criterion for geomaterials to describe the nonlinear strength properties with separated stress invariants. Yao et al. [8, 9] established a general nonlinear strength theory, which was extensively employed in concrete, rock, sand, clay, etc. Moreover, Yao et al. [10, 11] proposed a unified strength criterion (USC) for geomaterials to describe the triangle curved shape, which was the combination of SMP criterion and extended Mises criterion in the deviatoric plane, and then he defined a linear interpolation function between the two criteria to describe the shape function. The nonlinear strength envelope mentioned above can be described by hyperbolic, parabolic, exponential, power functions, etc. For structured soils, Liu et al.  proposed a new strength criterion based on binary-medium constitutive model (BMCM), and its expressions on both meridian and deviator planes are given. Due to the peculiar and complicated properties of frozen soils, most of the strength criteria are not completely suitable and should be modified for further applications based upon experimental results. Fish  and Ma et al.  gave a parabolic function to describe the strength properties. Qi and Ma  modified the classical Mohr–Coulomb strength criterion to apply for frozen sandy soils considering the influence of confining pressures. Lai et al.  established a parabolic curve for frozen silt soils in the meridian plane and Lade–Duncan model in the deviatoric plane to describe the shape function. Lai et al.  and Liao et al.  modified the hydrostatic pressure to simulate the strength in the q-p plane and employed a combination of LD and SMP models to describe the deviatoric plane. Lai et al.  proposed a relationship between the first principal σ1 and the third principal stress σ3, and then they established a nonlinear strength formula for frozen sandy soils. The researches mentioned above mainly concentrate on the nonlinear properties of unfrozen soils and frozen soils. The classical and modified strength criterion cannot completely reflect the nonlinear characteristics for frozen soils, so it is necessary to establish an appropriate strength criterion to describe the nonlinear properties of frozen soils.
At the same time, many research results of the stress-strain curves were obtained, especially in rock damage mechanism, mechanical properties, and strength and deformation characteristics [20–29]. On the theoretical framework of continuum damage mechanics (CDM), probability, and statistics theory, the damage constitutive models were established and could well simulate the stress-strain process for geological materials. For instance, Ren et al.  studied microdamage mechanism and damage constitutive model on the basis of uniaxial compression tests and found that the damage mechanism could appropriately describe the variation of stress-strain relationship. Later, a series of experimental study on rock were conducted by Zhang et al.  to analyze the weakened mechanism and damage properties under the freezing-thawing conditions. Ren  carried out computed tomography (CT) tests to study breakage mechanism for frozen cracked rock from a mesoscopic point of view. Lai et al.  and Li et al.  proposed a new and an improved statistical damage constitutive model for warm frozen clay and warm ice-rich frozen clay based upon the experimental results and applied Mohr–Coulomb criterion to judge whether the frozen soil element is damaged or not. Xu et al.  made an investigation on strength and deformation characteristics of ice-saturated frozen sandy soil and then proposed a nonlinear strength criterion and simultaneously employed an improved Duncan–Chang model to simulate the deformation properties. The researches mentioned above are the latest research achievements.
Due to the peculiar and complicated characteristics of frozen soils, some strength criteria, such as Mohr–Coulomb, Drucker–Prager, and Hoek–Brown strength criterion, are not suitable for frozen soils. So, the classical strength criterion should be modified or a new strength criterion may be established to investigate the nonlinear strength properties and breakage mechanism for frozen soils. Therefore, in this paper, a relationship between the first and third principal stresses is proposed, and then the nonlinear strength properties of frozen sandy soils are investigated. Furthermore, a new statistical damage constitutive model is also proposed, and the model parameters are determined by cryogenic triaxial compression tests data. Finally, the applicability of the new constitutive model is validated by comparisons between predicted and experimental data.
2. Test Result Analysis
As is illustrated in Figure 1 from Lai et al. , it can be found that the stress-strain curves go through the linear stage under relatively small axial strain to the elastoplastic stage. In addition, the curves present strain-softening phenomenon under low confining pressures and strain-hardening phenomenon under relatively high confining pressures. From the volumetric strain-axial strain curves, it is indicated that the volume of specimen is compressed at first and then dilated under low confining pressures, while it is merely compressed under high confining pressures. The determination method of peak values of stress-strain curves can be obtained as follows: (a) when the stress-strain curves present strain softening, the maximum value of is taken as the strength of frozen soil in this paper and (b) when the stress-strain curves present strain hardening, the value of at axial strain is taken as the strength of frozen soil . The relationship between deviatoric stress q and mean stress p is depicted in Figure 2, which indicates that the strength increases first and then decreases with increasing mean stress p.
3. Nonlinear Strength Theory of Frozen Sandy Soils
3.1. Existing Classical Strength Criterion
3.1.1. Mohr–Coulomb Strength Criterion
In geomaterial engineering, the Mohr–Coulomb criterion is widely used to predict the strength and deformation properties; hence, one of the classical strength criteria can be described as follows:where and denote the cohesive force and internal frictional angle, respectively. In this paper, for the test data, and .
3.1.2. Hoek–Brown Nonlinear Criterion
In the early years, Hoek  proposed the Hoek–Brown criterion based upon hundreds of experimental results and numerous field test findings for rock materials. It is an empirical strength criterion, and the expression can be presented in the following form:where , , and are material parameters. In this paper, for the test data, , , and .
3.2. The Evolution of the Proposed Strength Criterion
Based upon the experimental data from Figure 1, we can obtain the relationship between and , as illustrated in Figure 3. In order to investigate the nonlinear strength characteristics, we propose an empirical equation as follows:where , , , and are material parameters, respectively, which can be determined based on experimental results. In this paper, for frozen sandy soils, , , , and and the fitting coefficient .
For convenience, in the σ-τ plane, (3) can be rewritten as
According to experimental results, compared with the classical strength criteria, such as Mohr–Coulomb strength criterion and Hoek–Brown nonlinear criterion, the proposed nonlinear strength results are more close to the experimental data, as is depicted in Figure 4. It demonstrates that the new strength criterion can simulate the experimental results of frozen sandy soils well.
3.3. Nonlinear Strength Theory of Frozen Sandy Soils
Figure 5 shows the graphical illustration of the relationship between stress state A(σ, τ) and failure envelop f. That is, when the stress state A(σ, τ) is at the top of Mohr’s circle, the failure envelope is represented by . It should be noted that the failure loci presents a nonlinear changing tendency, which is also called a failure state line. The experimental maximum stress in the failure envelope is tangent to Mohr’s circle, which is the intersection point of A(σ, τ). From Figure 5, it is known that the failure stress state is represented by the point A, which cannot transcend the failure loci . The mathematical relationship of Mohr’s circle can be described as
Due to the independence of and , the relationship of can be calculated as
Differentiating (4) for and , respectively, we can obtain the following:
Furthermore, we can easily get the expression of internal frictional angle based upon the triangle relationship shown in Figure 5 as follows:
Differentiating (4) for , the result is presented as follows:
Finally, substituting (18) into (14), (15), and (16), in the σ-τ plane, the expressions of normal stress σ, shear stress τ, and internal frictional angle φ of the nonlinear strength criterion can be rewritten aswhere the expressions of (19) and (20) can be known as failure envelope in the σ-τ plane for frozen sandy soils.
Based upon experimental data, the nonlinear strength envelope can be obtained from (19) and (20), as is illustrated in Figure 6, in which the same parameters in Section 3.2 are used. The schematic diagram distinctively presents the nonlinear strength properties. Compared with experimental results, the proposed expression of (4) about and can well grasp and simulate the nonlinear strength characteristics of frozen sandy soils. It can be found that the proposed relationship is closer to the experimental results than Mohr–Coulomb strength criterion or Hoek–Brown strength criterion. It is convinced that the proposed equation can appropriately reflect the nonlinear strength properties with the increase of confining pressures.
4. Damage Statistical Constitutive Model for Frozen Sandy Soils
4.1. Formulation of the Damage Statistical Constitutive Model
Similar to unfrozen geomaterials, such as rock, concrete, rockfill material, sand, and clay, the damage process can be mainly accounted for by the same token. That is, the degradation of the material makes the effective areas decrease and the effective stress increase. Accordingly, the intact specimens are easy to damage due to the appearance of fissures, cracks, defects, and shear bands.
Based on the strain equivalent principle , the strain caused by apparent stress applied to a damaged material is equal to the strain caused by equivalent stress acting on the undamaged material. Hence, the tensor expression can be described as follows:where and denote the elastic flexibility matrix tensors of damaged material and undamaged material, respectively; and represent the equivalent stress matrix tensor and apparent stress matrix tensor, respectively; represents the strain matrix tensor, is an identity matrix tensor, and is a damage matrix tensor.
Based on phenomenological method, basic assumptions by Lai et al. [33, 39] are carried out for frozen soils. At first, on a macroscale, the soil specimen, regarded as a representative volume element (abbreviated as RVE), is isotropic and contains the basic information of damage, while it is a microheterogeneous material in the mesolevel. In addition, the linear elastic law is applied when the frozen soil element remains undamaged, and the nonlinearity of the stress-strain relationship is derived from the damage of the material. So, we can take the initial tangent modulus as the elastic modulus of the undamaged material. The damage variable is defined as the ratio of damaged section area to total section area in the mesolevel, and in macroscopic view, it is defined as the ratio of the number of damaged frozen soil elements to the number of all frozen soil elements. So, the damage expression defined is given as follows:where denotes the number of damaged frozen soil elements and stands for the number of all frozen soil elements, respectively.
4.2. Weibull Distribution Function
Compared with the other probability distribution function , it is found that the Weibull distribution can well describe the breakage process of rock [20–22] and warm frozen silt clay [33, 39]. In the paper, it is assumed that the microstrength of frozen sandy soil obeys the Weibull distribution as well. When the stress level reaches the value of F as is illustrated in the later section, the internal flaws and fissures gradually increase and transform to macrocracks or shear bands. So, the expression of probability density for damaged frozen sandy soil elements can be obtained as follows:where represents the probability density function of microstrength for frozen sandy soil, F denotes the yield criterion, and m and F0 are material parameters of the Weibull distribution function.
The yield criterion for frozen sandy soils can be written aswhere denotes the failure strength, which varies with internal state parameters, such as stress level, stress history, stress path, cohesive force, and internal frictional angle.
By integrating the microstrength function of (24), we can get the number of damaged frozen sandy soil elements as follows:
4.3. A New Strength Criterion for Frozen Sandy Soil
Based upon the framework of critical state soil mechanics , for unfrozen soils as noncohesive materials, the strength envelope is a straight line and is passing through the original point of the coordinates in the meridian plane. With respect to frozen soils, the bonding effect of ice crystals and soil particles possesses tensile strength, which will be weakened due to the pressure melting under relatively high confining pressures, similar to the cohesive behaviors of cemented clay. In the previous study [13, 14, 17, 18, 37], it is indicated that the strength of frozen soils increases first and then decreases with the increase of confining pressures. The maximum value qmax of the q-p curve corresponds to failure mean stress pcr. When the stress p is lower than pcr, the q-p curve can be replaced by linear strength criterion in the meridian plane, such as Mohr–Coulomb criterion and Druck–Prager criterion. While the stress p is higher than pcr, the q-p curve starts to bend downward due to pressure melting. In order to simulate evolution laws of the failure state for frozen sandy soils of the phenomenon, Fish , Ma et al. [14, 37], and Lai et al.  suggested a parabolic formula for frozen soils in the meridian plane, and Lai et al.  meanwhile proposed the Lade–Duncan model to simulate the strength characteristics in the deviatoric plane. Later, Nguyen et al.  proposed that the CSL of the cemented clay eventually reduced to an asymptote coinciding with the CSL of natural clay. Lai et al.  and Liao et al.  employed a modified mean effective stress expression to describe the strength criterion for frozen soils.
Hence, in this paper, based on the research results of Nguyen et al. , Lai et al. , and Liao et al. , we propose a modified strength criterion to simulate the variation tendency of the first increasing and then decreasing phenomena for frozen sandy soils. Here, the graphic illustration and the proposed modified strength criterion are given in Figure 8 and (28) to (30).
From Figure 8, the modified failure state line can be written as follows:
Substituting (29) into (28), the expression can be rewritten as follows:where is the deviatoric stress, is the mean stress, and is the intercept of the strength curve in the meridian plane. is the initial failure stress ratio, equal to the slope of the curve. The parameter is related to pressure melting, which reflects the changing rate about with the increase of hydrostatic pressure, and denotes the material parameter.
As is depicted in Figure 8, when the mean stress p is lower than pcr, the modified failure state line can describe the strengthened effect; when the mean stress p is higher than pcr, it denotes the weakened effect for frozen sandy soils.
The mean stress p and the deviatoric stress q can be written in the following form:where , , and are the first, second, and third principal stresses, respectively; J2 is the second deviatoric stress variant; is the stress tensor; is the deviatoric stress tensor; and is the Kronecker function.
Based upon the assumptions, it is easily accepted that the microstrength of frozen soils satisfies the modified failure state line; so, the strength criterion expression of the effective stress can be obtained as follows:where and denote the effective mean stress and deviatoric stress, respectively.
According to (22), the effective stress and apparent stress can be rewritten as follows:
The elastic constants λ and G can be obtained as follows:where λ denotes the Lame constant and K and G are the bulk modulus and shear modulus, respectively.
Based upon the generalized Hooke’s law, the expression of axial strain can be expressed by
Hence, the effective mean stress and deviatoric stress can be obtained as follows:
In order to obtain the material parameters, (43) can be converted to the following form:
4.4. Parameter Determination
4.4.1. Determination of Bulk Modulus K and Shear Modulus G
In order to describe the changing evolution of the bulk modulus K and shear modulus G, the mathematical relationship can be expressed by
where Pa is the atmospheric pressure (Pa = 0.10133 MPa) and , , and , , are material parameters. In this paper, the fitting results are illustrated in Figure 9: , , and , , .
4.4.2. Parameter Determination of the Strength Criterion
According to the cryogenic triaxial compression test data as is illustrated in Figure 1, the values of fitting parameters are shown in Table 2 and the corresponding fitting curve in Figure 10. It can be found that the proposed criterion has a good application for frozen sandy soils, especially the correlation coefficient R2 reaching the value of 0.9971.
4.4.3. Parameter Determination of the Damage Statistical Constitutive Model
The values of fitting parameters of (46) can be obtained from triaxial compression test data. Hence, m and F0 can be found in Table 3. The fitting results of (46) are depicted in Figure 11. It can be found that the relationship between Y and X presents a good linear relationship.
As presented in Figure 12, the relationship between the parameter m and confining pressures can be expressed as follows:where means 1.0 MPa to keep the dimensionless.
4.5.1. Verification of Deviatoric Stress-Axial Strain
After determining the model parameters, the stress under different confining pressures can be calculated by (37). Comparisons of the curves between experimental data and predicted results are depicted in Figure 13. It can be found that this proposed statistical damage constitutive model can simulate the softening phenomenon under low confining pressures and hardening phenomenon under high confining pressures. As a whole, the proposed model has a relatively high accuracy, so it is quite convenient for practical engineering.
4.5.2. Verification of Volumetric Strain-Axial Strain
In order to simulate the volumetric strain of frozen sandy soils, we give a damage relationship from the mesolevel, and the expression is shown as follows:where denotes the initial volume of the specimen, denotes the damage volume of the specimen, and denotes the total volume of the specimen.
Based on the theory framework of continuum damage mechanics (CDM), the damage variable can also be redefined as
So, the strain tensor of (51) can be rewritten as follows:
Under the triaxial symmetric compression condition, the volumetric strain of (53) can be expressed as
We make an assumption of for the initial elastic part and for the damage part; hence, the following expression can be obtained:where and are the material parameters, which can be determined by experimental data.
From Tables 3 and 4, the curves of volumetric strain-axial strain can be obtained as are illustrated in Figure 14. It can be found from Figure 14 that (55) can well simulate the experimental data, which present firstly compression and then dilation phenomenon under relatively low confining pressures and merely compression under high confining pressures.
In this paper, a new nonlinear strength criterion is proposed, and a damage statistical constitutive model is also established. The following conclusions can be reached:(1)The relationship between and is proposed based upon experimental results of frozen sandy soils. It is illustrated that the proposed equation is more close to experimental results than the classical strength criteria, such as Mohr–Coulomb and Hoek–Brown strength criteria. And then, the equation is introduced into the σ-τ plane to investigate the nonlinear strength properties based on failure envelope theorem. Compared with test data, the predicted results have a good accuracy and appropriateness.(2)The microstrength of frozen sandy soils obeys the Weibull distribution function. The strength criterion, which contains the damage properties of microstrength, is chosen as an independent variable in damage variable D. The parameters related to the damage statistical constitutive model are determined by experimental data, and they can well simulate both the curves of deviatoric stress-axial strain and volumetric strain-axial strain.
Highlights. (1) A new strength criterion is proposed based on experimental results, and the nonlinear strength properties are considered. (2) The stress and volumetric strain-axial strain at different confining pressures are validated by damage statistical constitutive model. (3) The damage constitutive model can simulate the strain softening and hardening phenomenon for frozen sandy soils.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
The authors appreciate the funding provided by the CAS Pioneer Hundred Talents Program (Dr. Liu Enlong) and National Science Foundation of China (41771066).
X. Z. Xu, J. C. Wang, and L. X. Zhang, Physics of Frozen Soils, Science press in Chinese, Beijing, China, 2009.
H. M. French, The Periglacial Environment, Essex, London, UK, 2nd edition, 1996.
N. A. Tsytovich, The Mechanics of Frozen Ground, Science Press, Beijing, China, 1985.
A. M. Fish, “Strength of frozen soil under a combined stress state,” Sixth International Symposium on Ground Freezing, vol. 1, pp. 135–145, 1991.View at: Google Scholar
W. Ma, Z. W. Wu, and Y. Sheng, “Effect of confining pressure on strength behaviour of frozen soil,” Chinese Journal Geotechnical Engineering, vol. 17, no. 5, pp. 7–11, 1995.View at: Google Scholar
W. G. Cao, Z. L. Fang, and X. J. Tang, “A study of statistical constitutive model for softening and damage rocks,” Chinese Journal of Rock Mechanics and Engineering, vol. 17, no. 6, pp. 628–633, 1998.View at: Google Scholar
W. G. Cao, M. H. Zhao, and C. X. Liu, “A study on damage statistical strength theory for rock,” Chinese Journal of Geotechnical Engineering, vol. 26, no. 6, pp. 820–823, 2004.View at: Google Scholar
W. G. Cao, P. Li, and M. H. Zhao, “On statistical damage constitutive model and its parameters for rock based on normal distribution,” Chinese Hydrogeology & Engineering Geology, vol. 32, no. 3, pp. 11–14, 2008.View at: Google Scholar
W. Y. Xu and L. D. Wei, “Study on statistical damage constitutive model of rock,” Chinese Journal of Rock Mechanics and Engineering, vol. 21, no. 6, pp. 787–791, 2002.View at: Google Scholar
J. X. Ren and X. R. Ge, “Study of rock meso-damage evolution law and its constitutive model under uniaxial compression loading,” Chinese Journal of Rock Mechanics and Engineering, vol. 20, no. 4, pp. 425–431, 2001.View at: Google Scholar
S. J. Zhang, Y. M. Lai, and X. M. Su, “A laboratory study on the damage propagation of rocks under freeze-thaw cycle condition,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 24, pp. 4105–4111, 2004.View at: Google Scholar
J. X. Ren, “Real-time CT test of damage failure mechanism of frozen cracked rock in loading and unloading condition,” Chinese Journal of Geotechnical Engineering, vol. 26, no. 5, pp. 641–644, 2004.View at: Google Scholar
E. Hoek, “The Hoek-Brown failure criterion,” in Proceedings of 15th Canadian Rock Mechanics Symposium, pp. 31–38, Department of Civil Engineering, University of Toronto, 1988.View at: Google Scholar
W. Weibull, “A statistical distribution function of wide applicability,” Journal of Applied Mechanics, vol. 18, pp. 293–297, 1951.View at: Google Scholar
K. H. Roscoe and J. B. Burland, “On the generalized stress-strain behaviour of wet clay,” in Engineering Plasticity, Cambridge University Press, Cambridge, UK, 1968.View at: Google Scholar