考虑孔隙水微观赋存形态的非饱和粉土有效应力方程及其验证

曾立峰; 邵龙潭; 牛庚; 郭晓霞

doi:10.16030/j.cnki.issn.1000-3665.202110023

考虑孔隙水微观赋存形态的非饱和粉土有效应力方程及其验证

1.
工业装备结构分析国家重点实验室（大连理工大学）, 辽宁大连　116024

2.
大连理工大学工程力学系, 辽宁大连　116024

3.
青岛理工大学理学院, 山东青岛　266033

基金项目: 国家自然科学基金面上项目（52079018；51479023）

详细信息

作者简介: 曾立峰（1986-），男，博士，主要从事非饱和土力学方面的研究工作。E-mail: zenglifeng34@163.com

通讯作者: 邵龙潭（1963-），男，博士研究生，教授，博士生导师，主要从事土和孔隙介质力学基本理论研究、土工结构稳定分析、土工试验测试技术与仪器研发等方面的教学和科研。E-mail: shaolt@dlut.edu.cn

中图分类号: P642.11⁺6

收稿日期: 2021-10-13

修回日期: 2021-11-20

刊出日期: 2022-07-25

An effective stress equation for unsaturated silt considering the microstructure of pore water and its verification

1.
State Key Laboratory of Structural Analysis for Industrial Equipment （Dalian University of Technology）, Dalian, Liaoning　116024, China

2.
Department of Engineering Mechanics, Dalian University of Technology, Dalian, Liaoning　116024, China

3.
School of Science, Qingdao University of Technology, Qingdao, Shandong　266033, China

More Information

Corresponding author: SHAO Longtan, shaolt@dlut.edu.cn

Received Date: 13 October 2021

Revised Date: 20 November 2021

Publish Date: 25 July 2022

摘要

摘要:
有效应力参数的合理确定是非饱和土有效应力研究的重要内容。然而，现有的有效应力参数未能较好地考虑孔隙水的微观赋存形态对有效应力的影响。为此，分析了孔隙水的微观赋存形态，明确了孔隙水可分为收缩膜、吸附水和毛细水，建立了非饱和粉土的扩展三相孔隙介质模型，即孔隙气、毛细水和广义土骨架。基于该模型，采用分相平衡分析法，推导了非饱和粉土的有效应力方程，该方程通过采用毛细水有效饱和度合理地考虑了毛细水对有效应力的影响。依据推导的有效应力方程，提出了非饱和粉土的抗剪强度公式。采用5种不同围压下重塑非饱和粉质黏土试验数据验证了提出的公式。结果表明：与现有的抗剪强度公式相比，本次提出的公式更好地预测了重塑非饱和粉质黏土的抗剪强度随基质吸力的变化规律。因此，运用提出的抗剪强度公式评价非饱和土边坡稳定性、地基承载力和挡土墙稳定性将会取得更好的效果。
- 非饱和粉土 /
- 孔隙水 /
- 微观赋存形态 /
- 有效应力 /
- 抗剪强度 /
- 持水特征曲线
Abstract:
Accurate determination of the effective stress coefficient is an important thing in the study of the effective stress of unsaturated soils. However, the existing effective stress coefficients have not well considered the effect of the microscopic occurrence of pore water in unsaturated soils on the effective stress. For this reason, the microscopic occurrence of pore water is analyzed firstly. It is recognized that the pore water in unsaturated soils can be divided into the contractile skin, adsorbed water and capillary water, and an extended three-phase porous medium model of unsaturated silt is established, considering the pore air, capillary water, and generalized soil skeleton. Based on the established model, an effective stress equation for unsaturated silt is derived by using an independent phase balance approach. The equation reasonably considers the influences of the capillary water on the effective stress by adopting a coefficient (i.e., the effective saturation of capillary water). In light of the derived effective stress equation, a shear strength equation for unsaturated silt is proposed. Finally, the proposed shear strength equation is verified by the experimental data on a reconstituted unsaturated silty clay under five different net confining pressures. The results show that compared with the existing two shear strength equations, the proposed equation better predicts the variation of the shear strength of the reconstituted unsaturated silty clay with matric suction. Therefore, under the unsaturated conditions, using the proposed shear strength formula to evaluate the stability of slopes, the bearing capacity of foundation, and the stability of retaining wall will achieve better results.
- unsaturated silt /
- pore water /
- microscopic occurrence /
- effective stress /
- shear strength /
- SWRC

HTML全文

非饱和土的有效应力方程是非饱和土力学的基本方程^[1]。对于理论问题而言，非饱和土的有效应力方程是研究非饱和土的强度^[2-3]、变形^[4]和水-力耦合^[5]等问题的基础；对于工程问题而言，基于非饱和土有效应力方程的抗剪强度公式因形式简单，物理意义清晰，已被广泛地应用于求解非饱和土的边坡稳定性^[6]、地基承载力^[7]和挡土墙土压力^[8]等工程问题。鉴于非饱和土有效应力方程在理论和工程问题上的重要性，自20世纪50年代以来，大量的学者从不同的角度提出了众多的非饱和土有效应力方程。其中最具代表性的是Bishop方程^[9]：

${\sigma {\text{'}}} = (\sigma - {u_{\text{a}}}) + \chi ({u_{\text{a}}} - {u_{\text{w}}})$

(1)

式中： ${\sigma \text{'}}$ ——非饱和土的有效应力/kPa；

$(\sigma - {u_{\text{a}}})$ ——净应力/kPa；

$({u_{\text{a}}} - {u_{\text{w}}})$ ——基质吸力/kPa；

$\chi$ ——有效应力参数，与饱和度有关，其值介于0～1之间。

后来，部分学者^[10-11]认为有效应力参数等于饱和度，在Mohr-Coulomb准则的框架下验证了其对砂土和部分粉土的适用性。另有学者^{[3, 12]}认为有效应力参数等于有效饱和度，采用同样的方法验证了其对砂土和绝大部分粉土的适用性。此外，还有学者^[13-14]认为有效应力参数等于其它形式的饱和度并且从强度的角度验证了各自的有效性。以上研究极大地推动了非饱和土有效应力方程的发展，然而它们未能较好地考虑孔隙水的微观赋存形态对有效应力的影响，这在一定程度上制约了有效应力方程的适用性。

鉴于以上问题，本文将分析非饱和粉土中孔隙水的微观赋存形态；在孔隙水微观分析的基础上，推导非饱和粉土的有效应力方程；从抗剪强度的角度验证推导的有效应力方程。

1. 非饱和粉土中孔隙水的微观赋存形态

一般来讲，非饱和土是由矿物颗粒、水和空气组成的三相系统。然而，由于三相之间存在着极其复杂的物理化学作用，导致非饱和土中矿物颗粒、水和空气已不同于它们各自单独存在时的状态^[15]。对于孔隙水而言，在一个非饱和粉土代表性体积单元（representative volume element）中，根据孔隙水的物理、化学和力学性质可以将其分为收缩膜（contractile skin）、吸附水（adsorbed water）和毛细水（capillary water），如图1所示。这3类形态的孔隙水反映了非饱和粉土中孔隙水的微观赋存特性。

图 1. 非饱和粉土RVE中孔隙水的微观赋存形态

Figure 1. Microstructure of the pore water in an RVE of unsaturated silt

下载: 全尺寸图片幻灯片

1.1 收缩膜

对于非饱和土而言，从力平衡的角度来讲，水-气交界面的具体形态与孔隙水压强、孔隙气压强、水的表面张力及接触角密切相关。当系统平衡后，水-气交界面一般出现在颗粒接触点附近。由于表面张力的存在，水-气交界面会像张紧的弹性薄膜一样将接触点附近的颗粒聚拢在一起，从而增加非饱和土的强度和刚度^[16-17]。为了形象地描述水-气交界面的力学作用，常将其称为收缩膜^[18]，如图1所示。收缩膜对非饱和土力学行为的影响主要取决于收缩膜的多少，而收缩膜的多少可以通过非饱和土中收缩膜的比表面积表征。当土体由饱和状态逐渐变为烘干状态时，收缩膜的比表面积先增大到某一个峰值然后逐渐减小^[19]。

1.2 吸附水

当矿物颗粒与水接触以后，它们之间会形成固-液交界面，其直观的现象是颗粒被一定厚度的吸附水膜包裹。吸附水膜的形成主要源于交界面附近的短程物理化学作用，这种短程作用主要包括电场力作用和范德华力作用^[20]。短程物理化学作用使得吸附水的势能显著降低、正压强明显提高。特别地，短程物理化学作用在黏性土中尤其显著，从而导致黏土矿物表面被包裹着较厚的吸附水膜^[15]。当土体由饱和状态逐渐变为烘干状态时，矿物颗粒表面的吸附水分子会不断地转换为孔隙气中的水分子和毛细水中的水分子，进而通过孔隙气和毛细水排出土体。因此，吸附水膜的厚度会逐渐减小。当土体处于饱和状态时，吸附水膜的厚度达到最大值；当土体接近烘干状态时，吸附水膜的厚度达到最小值^[21]。吸附水对非饱和土力学行为的影响在很大程度上取决于吸附水膜的厚度。

1.3 毛细水

当烘干状态的土体与水接触以后，矿物颗粒表面首先会形成一定厚度的吸附水膜。随着饱和度的增加，颗粒接触点周围会形成互不连通的弯月水，常被称为悬索状毛细水（pendular water）。当饱和度进一步增加时，悬索状毛细水会相互连通，形成纤维状毛细水（funicular water），如图1所示。毛细水的基质势或应力状态主要受水的表面张力、矿物颗粒与水之间的接触角以及矿物颗粒围成的孔隙的几何特性（如孔隙的形状、尺寸、分布及连通性等）控制，在微观上，其可以由Young-Laplace方程表征。当基质吸力大于某一较大值时，毛细水会发生空化，从而使得毛细水逐渐减少直至消失。与吸附水的正压强状态不同的是，毛细水呈现出负压强状态，这种负压强状态有利于毛细水周围的颗粒聚集，从而提高非饱和土的强度和刚度^[22]。尽管从微观上看，毛细水具有不同于常规液态水的属性，但是从宏观上看，毛细水仍然呈现出液态水的流动行为。

2. 考虑孔隙水微观赋存形态的有效应力方程

2.1 非饱和粉土的扩展三相孔隙介质模型

基于孔隙水微观赋存形态的分析结果，可以得知：位于颗粒接触点附近的收缩膜是一层张紧的弹性薄膜，所以其能够承担并传递荷载；位于颗粒表面的吸附水膜具有一系列不同于常规液态水的独特性质，例如高度规律排列的结构^[23]、更高的密度^[24]、更大的黏度^[25]以及更高的剪切模量^[26]。正是由于这些独特的物理力学性质，吸附水膜表现出类似于固体物质的行为，这意味着吸附水可以被视为一种能够承担并传递荷载的结构。与此同时，非饱和土中还存在能够承担并传递荷载的核心结构—矿物颗粒骨架。鉴于收缩膜、吸附水膜（其含量会随饱和度或基质吸力发生变化）和矿物颗粒骨架均能承担并传递荷载，在此将三者组合成一个广义相，称为广义土骨架相（generalized soil skeleton）。于是，广义土骨架、毛细水和孔隙气组成了非饱和土，在此情况下，非饱和粉土可用扩展的三相孔隙介质模型描述，如图2所示。

图 2. 非饱和粉土的扩展三相孔隙介质模型

Figure 2. An extended three-phase porous medium model for unsaturated silt

下载: 全尺寸图片幻灯片

基于扩展的三相孔隙介质模型，在一个非饱和粉土RVE内可以定义一组明确的体积—饱和度关系（图3）。非饱和粉土RVE内各单独相的体积应满足以下关系：

图 3. 非饱和粉土的扩展三相孔隙介质模型的体积—饱和度关系

Figure 3. Schematic diagram of the volume-saturation relations of an extended three-phase porous medium model for unsaturated silt

下载: 全尺寸图片幻灯片

${V_{\text{a}}} + {V_{{\text{cw}}}} + {V_{{\text{gs}}}} = V$

(2)

式中：V_a——孔隙气相的体积/m³；

V_cw——毛细水相的体积/m³；

V_gs——广义土骨架相的体积/m³；

V——非饱和粉土RVE的总体积/m³。

非饱和粉土RVE中孔隙气、毛细水和广义土骨架的体积分数可以表示为：

$\left\{ \begin{gathered} {n_{\text{a}}} = \frac{{{V_{\text{a}}}}}{V} \times 100\% \hfill \\ {n_{{\text{cw}}}} = \frac{{{V_{{\text{cw}}}}}}{V} \times 100\% \hfill \\ {n_{{\text{gs}}}} = \frac{{{V_{{\text{gs}}}}}}{V} \times 100\% \hfill \\ \end{gathered} \right.$

(3)

式中：n_a——孔隙气相的体积分数/%；

n_cw——毛细水相的体积分数/%；

n_gs——广义土骨架相的体积分数/%。

联立式（2）和式（3），可得：

${n_{\text{a}}} + {n_{{\text{cw}}}} + {n_{{\text{gs}}}} = 1$

(4)

在扩展的三相孔隙介质模型下，非饱和粉土的有效孔隙体积是孔隙气的体积与毛细水的体积之和：

$V_{\text{v}}^{{\text{eff}}}{\text{ = }}{V_{\text{a}}} + {V_{{\text{cw}}}}$

(5)

式中： $V_{\text{v}}^{{\text{eff}}}$ ——非饱和粉土的有效孔隙体积/m³。

非饱和粉土的有效孔隙的体积分数（或有效孔隙率）可以表示为：

$n_{\text{v}}^{{\text{eff}}} = \frac{{V_{\text{v}}^{{\text{eff}}}}}{V} \times 100\%$

(6)

式中： $n_{\text{v}}^{{\text{eff}}}$ ——有效孔隙率/%。

在此情况下，饱和度的定义应为毛细水的体积与有效孔隙体积之比：

$S_{{\text{cw}}}^{{\text{eff}}} = \frac{{{V_{{\text{cw}}}}}}{{V_{\text{v}}^{{\text{eff}}}}}$

(7)

式中： $S_{{\text{cw}}}^{{\text{eff}}}$ ——扩展的三相孔隙介质模型的饱和度，本文称其为毛细水的有效饱和度。

2.2 有效应力方程的推导

在扩展的三相孔隙介质模型建立之后，可以在其基础上采用分相平衡分析法^[1]推导非饱和粉土的有效应力方程。

选取一个非饱和土整体微元体并对其进行受力分析，见图4（a）。在此整体微元体中分离出孔隙气相、毛细水相和广义土骨架相并分别分析各单独相的受力情况，见图4（b）—（d）。由于收缩膜已被视为广义土骨架相的组成部分，所以孔隙气相与毛细水相之间并不直接接触，这就意味着孔隙气相与毛细水相之间不存在相互作用力。

图 4. 非饱和土整体及各分相的微元体的受力示意

Figure 4. Illustration for the forces acting on the unsaturated soil, pore air, capillary water and generalized soil skeleton

下载: 全尺寸图片幻灯片

非饱和土整体微元体的力平衡微分方程为：

$\left\{ \begin{gathered} \frac{{\partial {\sigma _x}}}{{\partial x}} + \frac{{\partial {\tau _{\textit{z}x}}}}{{\partial \textit{z}}} = 0 \hfill \\ \frac{{\partial {\tau _{xz}}}}{{\partial x}} + \frac{{\partial {\sigma _\textit{z}}}}{{\partial \textit{z}}} + \rho g = 0 \hfill \\ \end{gathered} \right.$

(8)

式中： ${\sigma _x}$ 、 ${\sigma _{\textit{z}}}$ ——作用在非饱和土整体微元体上的总正应力/kPa；

${\tau _{\textit{z}x}}$ 、 ${\tau _{x\textit{z}}}$ ——作用在非饱和土整体微元体上的总剪应力/kPa；

$\rho$ ——非饱和土的密度/（kg·m⁻³）；

$g$ ——重力加速度/（m·s⁻²）。

孔隙气相微元体的力平衡微分方程为：

$\left\{ \begin{gathered} \frac{{\partial ({n_{\text{a}}}{u_{\text{a}}})}}{{\partial x}} + f_{{\text{a}}-x}^{{\text{gs}}} = 0 \hfill \\ \frac{{\partial ({n_{\text{a}}}{u_{\text{a}}})}}{{\partial \textit{z}}} + f_{{\text{a}}-\textit{z}}^{{\text{gs}}} + {n_{\text{a}}}{\rho _{\text{a}}}g = 0 \hfill \\ \end{gathered} \right.$

(9)

式中： ${u_{\text{a}}}$ ——孔隙气压强/kPa；

$f_{{\text{a}} - x}^{{\text{gs}}}$ 、 $f_{{\text{a}} -\textit{z}}^{{\text{gs}}}$ ——广义土骨架相对孔隙气相的作用力，其是一种体积力/（kN·m⁻³）；

${\rho _{\text{a}}}$ ——孔隙气的密度/（kg·m⁻³）。

毛细水相微元体的力平衡微分方程：

$\left\{ \begin{gathered} \frac{{\partial ({n_{{\text{cw}}}}{u_{\text{w}}})}}{{\partial x}} + f_{{\text{cw}-}x}^{{\text{gs}}} = 0 \hfill \\ \frac{{\partial ({n_{{\text{cw}}}}{u_{\text{w}}})}}{{\partial \textit{z}}} + f_{{\text{cw}-}\textit{z}}^{{\text{gs}}} + {n_{{\text{cw}}}}{\rho _{\text{w}}}g = 0 \hfill \\ \end{gathered} \right.$

(10)

式中： ${u_{\text{w}}}$ ——毛细水压强/kPa；

$f_{{\text{cw}} - x}^{{\text{gs}}}$ 、 $f_{{\text{cw}} - \textit{z}}^{{\text{gs}}}$ ——广义土骨架相对毛细水相的作用力/（kgN·m⁻³）；

${\rho _{\text{w}}}$ ——毛细水的密度/（kg·m⁻³）。

广义土骨架相微元体的力平衡微分方程：

$\left\{ \begin{gathered} \frac{{\partial \sigma _x \text{'}}}{{\partial x}} + \frac{{\partial \tau _{\textit{z}x} \text{'}}}{{\partial \textit{z}}} - f_{{\text{gs} - }x}^{\text{a}} - f_{{\text{gs} - }x}^{{\text{cw}}} + \frac{{\partial (a_{{\text{gs}}}^{\text{a}}{u_{\text{a}}})}}{{\partial x}} + \frac{{\partial (a_{{\text{gs}}}^{{\text{cw}}}{u_{\text{w}}})}}{{\partial x}} = 0 \hfill \\ \frac{{\partial \sigma _\textit{z} \text{'}}}{{\partial \textit{z}}} + \frac{{\partial \tau _{x\textit{z}} \text{'}}}{{\partial x}} - f_{{\text{gs} - }\textit{z}}^{\text{a}} - f_{{\text{gs}} - \textit{z}}^{{\text{cw}}} + {n_{{\text{gs}}}}{\rho _{{\text{gs}}}}g + \frac{{\partial (a_{{\text{gs}}}^{\text{a}}{u_{\text{a}}})}}{{\partial \textit{z}}} + \frac{{\partial (a_{{\text{gs}}}^{{\text{cw}}}{u_{\text{w}}})}}{{\partial \textit{z}}} = 0 \hfill \\ \end{gathered} \right.$

(11)

式中： $\sigma _x \text{'}$ 、 $\sigma _{\textit{z}} \text{'}$ ——作用在广义土骨架相微元体上的有效正应力/kPa；

$\tau _{\textit{z}x} \text{'}$ 、 $\tau _{x\textit{z}} \text{'}$ ——作用在广义土骨架相微元体上的有效剪应力/kPa；

$f_{{\text{gs} - }x}^{\text{a}}$ 、 $f_{{\text{gs} - }x}^{{\text{cw}}}$ 、 $f_{{\text{gs} - }\textit{z}}^{\text{a}}$ 、 $f_{{\text{gs} -{ \textit{z}}}}^{{\text{cw}}}$ ——孔隙气相和毛细水相对广义土骨架相的作用力/（kN·m⁻³）；

${\rho _{{\text{gs}}}}$ ——广义土骨架相的密度/（kg·m⁻³）；

$a_{{\text{gs}}}^{\text{a}}$ 、 $a_{{\text{gs}}}^{{\text{cw}}}$ ——孔隙气相和毛细水相在广义土骨架相微元体单位面积表面上的作用面积。

分开考虑孔隙气压强和毛细水压强作用下广义土骨架相微元体的受力情况，结合式（3）—（7），可以计算得：

$\begin{split}&a_{{\text{gs}}}^{\text{a}} = (1 - S_{{\text{cw}}}^{{\text{eff}}})(1 - n_{\text{v}}^{{\text{eff}}})\\ &a_{{\text{gs}}}^{{\text{cw}}} = S_{{\text{cw}}}^{{\text{eff}}}(1 - n_{\text{v}}^{{\text{eff}}}) \end{split}$

将孔隙气相、毛细水相和广义土骨架相微元体的平衡微分方程式（9）—（11）叠加并考虑到以下条件：

① $f_{{\text{a} - }i}^{{\text{gs}}} = f_{{\text{gs} - }i}^{\text{a}}$ 和 $f_{{\text{cw} - }i}^{{\text{gs}}} = f_{{\text{gs} - }i}^{{\text{cw}}}$ ，其中i=x, z；

② $\rho = {n_{\text{a}}}{\rho _{\text{a}}} + {n_{{\text{cw}}}}{\rho _{\text{w}}} + {n_{{\text{gs}}}}{\rho _{{\text{gs}}}}$ 。

可以得到由各单独相经叠加后的非饱和土整体的平衡微分方程：

$\left\{ \begin{gathered} \frac{{\partial \sigma _x \text{'}}}{{\partial x}} + \frac{{\partial \tau _{\textit{z}x} \text{'}}}{{\partial \textit{z}}} + \frac{{\partial [(1 - S_{{\text{cw}}}^{{\text{eff}}}){u_{\text{a}}}]}}{{\partial x}} + \frac{{\partial (S_{{\text{cw}}}^{{\text{eff}}}{u_{\text{w}}})}}{{\partial x}} = 0 \hfill \\ \frac{{\partial \sigma _\textit{z} \text{'}}}{{\partial \textit{z}}} + \frac{{\partial \tau _{x\textit{z}} \text{'}}}{{\partial x}} + \frac{{\partial [(1 - S_{{\text{cw}}}^{{\text{eff}}}){u_{\text{a}}}]}}{{\partial \textit{z}}} + \frac{{\partial (S_{{\text{cw}}}^{{\text{eff}}}{u_{\text{w}}})}}{{\partial \textit{z}}} + \rho g = 0 \hfill \\ \end{gathered} \right.$

(12)

比较式（12）和式（8），并考虑到通常情况下孔隙气相和毛细水相的抗剪强度可以忽略，可得 ${\tau _{\textit{z}x}} = \tau _{\textit{z}x} \text{'}$ 和 ${\tau _{x\textit{z}}} = \tau _{x\textit{z}} \text{'}$ 。张量方程为：

${{\boldsymbol{\sigma }}\text{'}} = {\boldsymbol{\sigma }} - [S_{{\text{cw}}}^{{\text{eff}}}{u_{\text{w}}} + (1 - S_{{\text{cw}}}^{{\text{eff}}}){u_{\text{a}}}]{{\boldsymbol{\delta}} }$

(13)

式中： ${{\boldsymbol{\sigma }}\text{'}}$ ——有效应力张量；

${\boldsymbol{\sigma }}$ ——总应力张量；

${{\boldsymbol{\delta}}}$ ——Kronecker符号，当i = j时， ${\delta _{ij}} = 1$ ；当i ≠ j时， ${\delta _{ij}} = 0$ 。

式（13）即是非饱和土的有效应力方程。该方程采用毛细水的有效饱和度 $S_{{\text{cw}}}^{{\text{eff}}}$ 作为有效应力参数，较好地考虑了孔隙水的微观赋存形态对有效应力的影响。

3. 有效应力方程的验证

3.1 验证思路

本次将从抗剪强度的角度验证有效应力式（13）。在Mohr-Coulomb准则框架下，基于有效应力式（13），破坏面上的抗剪强度式可表示为：

${\tau _{\text{f}}} = {c\text{'}} + \left\{ {\sigma - [S_{{\text{cw}}}^{{\text{eff}}}{u_{\text{w}}} + (1 - S_{{\text{cw}}}^{{\text{eff}}}){u_{\text{a}}}]} \right\}\tan {\varphi \text{'}}$

(14)

式中： ${\tau _{\text{f}}}$ ——破坏面上的抗剪强度/kPa；

${c\text{'}}$ ——饱和土的有效内黏聚力/kPa

${\varphi \text{'}}$ ——饱和土的有效内摩擦角/（°）。

有效应力式（13）的验证问题转变为抗剪强度式（14）的验证问题。为了验证式（14），首先需要确定计算毛细水有效饱和度 $S_{{\text{cw}}}^{{\text{eff}}}$ 的方法。考察 $S_{{\text{cw}}}^{{\text{eff}}}$ 的表达式（7），改写为如下形式：

$S_{{\text{cw}}}^{{\text{eff}}} = \frac{{{V_{{\text{cw}}}}}}{{V_{\text{v}}^{{\text{eff}}}}} = \frac{{{V_{{\text{cw}}}}}}{{{V_{\text{v}}} - {V_{{\text{aw}}}}}} = \frac{{\dfrac{{{V_{{\text{cw}}}}}}{{{V_{\text{v}}}}}}}{{\dfrac{{{V_{\text{v}}}}}{{{V_{\text{v}}}}} - \dfrac{{{V_{{\text{aw}}}}}}{{{V_{\text{v}}}}}}} = \frac{{{S_{{\text{cw}}}}}}{{1 - {S_{{\text{aw}}}}}}$

(15)

式中：V_aw——吸附水的体积/m³；

V_v——孔隙的体积/m³，V_v = V_a+ V_cw+ V_aw；

S_cw——毛细水的饱和度；

S_aw——吸附水的饱和度。

根据式（15）可知，当毛细水和吸附水的饱和度随基质吸力的变化规律明确以后， $S_{{\text{cw}}}^{{\text{eff}}}$ 随基质吸力的变化规律即可相应地明确。为此，本文采用了Lu^[21]建立的持水特征曲线（SWRC）模型。该模型不仅可以分别量化毛细水和吸附水的饱和度随基质吸力的变化规律，而且也可以考虑空化效应（cavitation）对持水行为的影响，所以其在物理机制上更为合理。这意味着通过Lu^[21]的SWRC模型计算的 $S_{{\text{cw}}}^{{\text{eff}}}$ 具有更高的精度，该模型的基本思想包括3 个方面：

①非饱和土的饱和度可分为毛细水的饱和度和吸附水的饱和度：

$S\left( \psi \right) = {S_{{\text{cw}}}}\left( \psi \right) + {S_{{\text{aw}}}}\left( \psi \right)$

(16)

式中：S——非饱和土的饱和度；

$\psi$ ——基质吸力/kPa。

②毛细水SWRC方程为：

${S_{{\text{cw}}}}\left( \psi \right) = \frac{1}{2}\left[ {1 - {\text{erf}}\left( {\sqrt 2 \frac{{\psi - {\psi _{{\text{cav}}}}}}{{{\psi _{{\text{cav}}}}}}} \right)} \right] \times \left[ {1 - {S_{{\text{aw}}}}\left( \psi \right)} \right] {\left[ {1 + {{\left( {\alpha \psi } \right)}^n}} \right]^{{\textstyle\frac{1}{n}} - 1}}$

(17)

式中： ${\psi _{{\text{cav}}}}$ ——平均空化基质吸力/kPa；

$\alpha$ ——进气值的倒数/kPa⁻¹；

n——与孔隙尺寸分布相关的参数；

erf( )——误差函数。

③吸附水SWRC方程为：

${S_{{\text{aw}}}}\left( \psi \right) = {S_{{\text{aw0}}}}\left\{ {1 - {{\left[ {\exp \left( {\frac{{\psi - {\psi _{\text{d}}}}}{\psi }} \right)} \right]}^l}} \right\}$

(18)

式中：S_aw0——基质吸力为1 kPa时吸附水的饱和度，即吸附水饱和度的最大值；

${\psi _{\text{d}}}$ ——烘干状态对应的基质吸力，即基质吸力的最大值，取为10⁶ kPa；

l——与吸附作用程度有关的参数。

考虑到砂土和一部分粉土试验数据难以有效地评价不同抗剪强度公式的优劣^[27-28]，本文从现有文献[29]中选取了一组重塑粉质黏土（成分为20%的Speswhite 高岭土、10%的伦敦黏土和70%的HPF4 硅微粉）的试验数据。该数据既包含一组广吸力范围内的SWRC数据，也包含5种不同围压（0，50，100，200，400 kPa）下的抗剪强度数据，这为式（14）的验证提供了较为全面和可靠的试验数据。另外，为了全面地评价抗剪强度式（14）的有效性和合理性，在此也给出了其它2 种通用的抗剪强度公式^{[10, 12]}：

${\tau _{\text{f}}} = {c\text{'}} + \left\{ {\sigma - [S{u_{\text{w}}} + (1 - S){u_{\text{a}}}]} \right\}\tan {\varphi \text{'}}$

(19)

$\begin{array}{l} {\tau _{\text{f}}} = {c\text{'}} + \left\{ {\sigma - [{S _{\rm{e}}}{u_{\text{w}}} + (1 - {S _{\rm{e}}}){u_{\text{a}}}]} \right\}\tan {\varphi \text{'}} \\ S_{\rm{e}}= (S-S_{\rm{res}})/(1-S_{\rm{res}})\end{array}$

(20)

式中：S_e——非饱和土的有效饱和度；

S_res——残余饱和度。

饱和度S、S_e和S_res的值均由VG模型确定。

3.2 验证结果

图5（a）显示了VG模型和Lu^[21]模型关于重塑粉质黏土SWRC数据的拟合结果。从图中可以看出，2 种模型对实测数据均具有较好的拟合结果，其中VG模型的拟合相关系数R²=0.9905，Lu^[21]模型的拟合相关系数R²=0.9908。当基质吸力超过28 MPa时，VG模型预测的饱和度随基质吸力的增加基本保持不变，Lu^[21]模型预测的饱和度随基质吸力的增加缓慢减小，直至基质吸力达到最大值10⁶ kPa。现有的实测数据^{[21, 30]}表明：在高吸力阶段，随着基质吸力的增加，饱和度仍然会不断减小。这就意味着Lu^[21]模型更符合实际的物理过程。

图 5. 两种SWRC模型的拟合结果

Figure 5. Fitting results of two SWRC models

下载: 全尺寸图片幻灯片

图5（b）显示了Lu ^[21]模型分解为毛细水SWRC和吸附水SWRC的结果。从图中可以看出，对于毛细水SWRC而言，随着基质吸力的增加毛细水饱和度先大致保持不变；当基质吸力达到800 kPa左右（即进气值附近）时，毛细水开始从土体孔隙中迅速排出，毛细水饱和度相应地迅速减小；当基质吸力达到4 900 kPa左右（即平均空化基质吸力附近）时，毛细水会发生空化现象。随着基质吸力的进一步增加，空化现象会不断地进行。当基质吸力达到10 MPa附近时，毛细水饱和度将趋近于0。然而，由于现有的SWRC模型^{[14, 31]}未能较好地考虑毛细水的空化现象，从而导致高吸力阶段毛细水饱和度普遍被高估。对于吸附水SWRC而言，当基质吸力小于4 900 kPa时，吸附水饱和度随基质吸力的增加基本保持不变；当基质吸力大于4 900 kPa时，吸附水饱和度开始迅速减小。当基质吸力进一步增加至43 MPa时，吸附水饱和度开始缓慢减小，直至基质吸力达到最大值10⁶ kPa。

图6显示了不同净围压下重塑非饱和粉质黏土的实测抗剪强度及3种抗剪强度公式，见式（14）（19）（20），的预测结果。从图中可以看出，对于实测抗剪强度而言，在相同的净围压下，重塑粉质黏土的抗剪强度随基质吸力的变化规律大致可以分为4 个阶段：①抗剪强度快速增加阶段，此阶段的基质吸力范围约为0～500 kPa（近似饱和阶段）；②抗剪强度缓慢增加阶段，此阶段的基质吸力范围约为500～900 kPa；③抗剪强度的峰值阶段，此阶段的基质吸力约为1 000 kPa；④抗剪强度下降阶段，此阶段的基质吸力超过1 000 kPa。与此同时，在相同的基质吸力下，重塑粉质黏土的抗剪强度随净围压的增加而增加。对于重塑粉质黏土的抗剪强度随基质吸力和净围压的变化规律，现有的2 种抗剪强度公式，式（19）（20），未能很好地预测这种变化规律，但是本文提出的抗剪强度式（14）较好地预测了这种变化规律（图6）。

图 6. 不同净围压下重塑非饱和粉质黏土的3种抗剪强度公式的预测结果

Figure 6. Predictions of three shear strength equations for a reconstituted unsaturated silty clay under different net confining pressures

下载: 全尺寸图片幻灯片

与式（14）相比，式（19）（20）在基质吸力超过500 kPa时严重地高估了重塑粉质黏土的抗剪强度，其中式（19）的高估程度更甚。这主要是因为当基质吸力超过500 kPa时，有效应力参数的大小会显著地影响有效应力和抗剪强度的大小。由于式（19）采用了饱和度S作为有效应力参数，认为毛细水和吸附水均对有效应力有贡献，从而高估了抗剪强度的大小。尽管式（20）通过采用有效饱和度S_e分离了部分吸附水，在一定程度上改进了式（19），但仍然高估了抗剪强度的大小。这主要是由于有效饱和度S_e中引入的残余饱和度S_res并不能描述吸附水饱和度随基质吸力的变化。与上述不同的是，式（14）采用了毛细水有效饱和度 $S_{{\text{cw}}}^{{\text{eff}}}$ 作为有效应力参数，从而合理地区分了毛细水和吸附水对有效应力的贡献，然后借助Lu^[21]SWRC模型量化了毛细水和吸附水饱和度随基质吸力的变化规律，最后较为准确地预测了重塑粉质黏土的抗剪强度。

4. 结论

（1）微观分析表明，非饱和粉土中的孔隙水可分为收缩膜、吸附水和毛细水。由于收缩膜和吸附水具有类固体的性质，它们可以与颗粒骨架组合在一起形成广义土骨架。因此，非饱和粉土可以用扩展的三相孔隙介质模型（即孔隙气、毛细水和广义土骨架）表征，其为非饱和粉土有效应力方程的推导提供了物理基础。

（2）通过采用毛细水有效饱和度，本文推导的非饱和土有效应力方程较好地考虑了毛细水对有效应力的影响。借助Lu模型，本文较为精确地确定了毛细水和吸附水饱和度随基质吸力的变化规律，这为准确计算有效应力奠定了基础。

（3）在不同的净围压下，当基质吸力超过500 kPa时，现有的2种抗剪强度公式均严重地高估了重塑非饱和粉质黏土的抗剪强度。但是，本文提出的抗剪强度公式较好地预测了不同净围压下重塑非饱和粉质黏土的抗剪强度随基质吸力的变化规律。囿于篇幅，本文仅从不同净围压的角度进行了验证。从不同土类的角度进行验证将是下一步工作的重点。

图 1 非饱和粉土RVE中孔隙水的微观赋存形态

Figure 1.

下载: 全尺寸图片幻灯片

图 2 非饱和粉土的扩展三相孔隙介质模型

Figure 2.

下载: 全尺寸图片幻灯片

图 3 非饱和粉土的扩展三相孔隙介质模型的体积—饱和度关系

Figure 3.

下载: 全尺寸图片幻灯片

图 4 非饱和土整体及各分相的微元体的受力示意

Figure 4.

下载: 全尺寸图片幻灯片

图 5 两种SWRC模型的拟合结果

Figure 5.

下载: 全尺寸图片幻灯片

图 6 不同净围压下重塑非饱和粉质黏土的3种抗剪强度公式的预测结果

Figure 6.

下载: 全尺寸图片幻灯片

参考文献(31)

[1]	邵龙潭, 郭晓霞. 有效应力新解[M]. 北京: 中国水利水电出版社, 2014 SHAO Longtan, GUO Xiaoxia. New interpretation of effective stress [M]. Beijing: China Water & Power Press, 2014. (in Chinese)
[2]	陈嘉伟, 高游, 付俊杰, 等. 不同类型黏土的强度特性及其预测[J]. 水文地质工程地质,2020,47(3):101 − 106. [CHEN Jiawei, GAO You, FU Junjie, et al. Strength of different clayey soils and its prediction[J]. Hydrogeology & Engineering Geology,2020,47(3):101 − 106. (in Chinese with English abstract) CHEN Jiawei, GAO You, FU Junjie, et al. Strength of different clayey soils and its prediction [J]. Hydrogeology & Engineering Geology, 2020, 47(3): 101-106. (in Chinese with English abstract)
[3]	邵龙潭, 覃亚龙, 郭晓霞, 等. 非饱和土有效应力公式验证[J]. 地下空间与工程学报,2018,14(6):1476 − 1483. [SHAO Longtan, QIN Yalong, GUO Xiaoxia, et al. The validation of the effective stress principle of unsaturated soils[J]. Chinese Journal of Underground Space and Engineering,2018,14(6):1476 − 1483. (in Chinese with English abstract) SHAO Longtan, QIN Yalong, GUO Xiaoxia, et al. The validation of the effective stress principle of unsaturated soils [J]. Chinese Journal of Underground Space and Engineering, 2018, 14(6): 1476-1483. (in Chinese with English abstract)
[4]	盛岱超, 杨超. 关于非饱和土本构研究的几个基本规律的探讨[J]. 岩土工程学报,2012,34(3):438 − 456. [SHENG Daichao, YANG Chao. Discussion of fundamental principles in unsaturated soil mechanics[J]. Chinese Journal of Geotechnical Engineering,2012,34(3):438 − 456. (in Chinese with English abstract) SHENG Daichao, YANG Chao. Discussion of fundamental principles in unsaturated soil mechanics [J]. Chinese Journal of Geotechnical Engineering, 2012, 34(3): 438-456. (in Chinese with English abstract)
[5]	韩博文, 蔡国庆, 李舰, 等. 考虑颗粒黏结效应的非饱和土水–力耦合边界面模型[J]. 岩土工程学报,2020,42(11):2059 − 2068. [HAN Bowen, CAI Guoqing, LI Jian, et al. Hydro-mechanical coupling bounding surface model for unsaturated soils considering bonding effect of particles[J]. Chinese Journal of Geotechnical Engineering,2020,42(11):2059 − 2068. (in Chinese with English abstract) HAN Bowen, CAI Guoqing, LI Jian, et al. Hydro-mechanical coupling bounding surface model for unsaturated soils considering bonding effect of particles [J]. Chinese Journal of Geotechnical Engineering, 2020, 42(11): 2059-2068. (in Chinese with English abstract)
[6]	饶鸿, 王金淑, 赵志明, 等. 基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析[J]. 水文地质工程地质,2021,48(1):154 − 162. [RAO Hong, WANG Jinshu, ZHAO Zhiming, et al. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology,2021,48(1):154 − 162. (in Chinese with English abstract) RAO Hong, WANG Jinshu, ZHAO Zhiming, et al. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model [J]. Hydrogeology & Engineering Geology, 2021, 48(1): 154-162. (in Chinese with English abstract)
[7]	徐永福. 膨胀土地基承载力研究[J]. 岩石力学与工程学报,2000,19(3):387 − 390. [XU Yongfu. Study on bearing capacity of expansive soil basement[J]. Chinese Journal of Rock Mechanics and Engineering,2000,19(3):387 − 390. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-6915.2000.03.030 XU Yongfu. Study on bearing capacity of expansive soil basement [J]. Chinese Journal of Rock Mechanics and Engineering, 2000, 19(3): 387-390. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-6915.2000.03.030
[8]	汪丁建, 童龙云, 邱岳峰. 降雨入渗条件下非饱和土朗肯土压力分析[J]. 岩土力学,2013,34(11):3192 − 3196. [WANG Dingjian, TONG Longyun, QIU Yuefeng. Rankine's earth pressure analysis of unsaturated soil under condition of rainfall infiltration[J]. Rock and Soil Mechanics,2013,34(11):3192 − 3196. (in Chinese with English abstract) WANG Dingjian, TONG Longyun, QIU Yuefeng. Rankine's earth pressure analysis of unsaturated soil under condition of rainfall infiltration [J]. Rock and Soil Mechanics, 2013, 34(11): 3192-3196. (in Chinese with English abstract)
[9]	BISHOP A W. The principle of effective stress[J]. Teknisk Ukeblad,1959,106(39):859 − 863.
[10]	ÖBERG A L, SÄLLFORS G. Determination of shear strength parameters of unsaturated silts and sands based on the water retention curve[J]. Geotechnical Testing Journal,1997,20(1):40 − 48. doi: 10.1520/GTJ11419J
[11]	LU N, LIKOS W J. Suction stress characteristic curve for unsaturated soil[J]. Journal of Geotechnical and Geoenvironmental Engineering,2006,132(2):131 − 142. doi: 10.1061/(ASCE)1090-0241(2006)132:2(131)
[12]	VANAPALLI S K, FREDLUND D G, PUFAHL D E, et al. Model for the prediction of shear strength with respect to soil suction[J]. Canadian Geotechnical Journal,1996,33(3):379 − 392. doi: 10.1139/t96-060
[13]	ALONSO E E, PEREIRA J M, VAUNAT J, et al. A microstructurally based effective stress for unsaturated soils[J]. Géotechnique,2010,60(12):913 − 925.
[14]	ZHOU A, HUANG R, SHENG D. Capillary water retention curve and shear strength of unsaturated soils[J]. Canadian Geotechnical Journal,2016,53(6):974 − 987. doi: 10.1139/cgj-2015-0322
[15]	MITCHELL J K, SOGA K. Fundamentals of soil behavior [M]. 3th ed. Hoboken: John Wiley & Sons Inc, 2005.
[16]	WANG J P, LAMBERT P, DE KOCK T, et al. Investigation of the effect of specific interfacial area on strength of unsaturated granular materials by X-ray tomography[J]. Acta Geotechnica,2019,14(5):1545 − 1559. doi: 10.1007/s11440-019-00765-2
[17]	LOUREN O S D N, GALLIPOLI D, AUGARDE C E, et al. Formation and evolution of water menisci in unsaturated granular media[J]. Géotechnique,2012,62(3):193 − 199.
[18]	FREDLUND D G, RAHARDJO H. Soil mechanics for unsaturated soils [M]. New York: John Wiley & Sons Inc, 1993.
[19]	COSTANZA-ROBINSON M S, BRUSSEAU M L. Air-water interfacial areas in unsaturated soils: Evaluation of interfacial domains[J]. Water Resources Research,2002,38(10):13 − 1.
[20]	LU N, LIKOS W J. Unsaturated soils mechanics [M]. New Jersey: John Wiley & Sons Inc, 2004.
[21]	LU N. Generalized soil water retention equation for adsorption and capillarity[J]. Journal of Geotechnical and Geoenvironmental Engineering,2016,142(10):04016051. doi: 10.1061/(ASCE)GT.1943-5606.0001524
[22]	TOLL D G, ALI RAHMAN Z. Critical state shear strength of an unsaturated artificially cemented sand[J]. Géotechnique,2017,67(3):208 − 215.
[23]	DHOPATKAR N, DEFANTE A P, DHINOJWALA A. Ice-like water supports hydration forces and eases sliding friction[J]. Science Advances,2016,2(8):e1600763. doi: 10.1126/sciadv.1600763
[24]	JACINTO A C, VILLAR M V, LEDESMA A. Influence of water density on the water-retention curve of expansive clays[J]. Géotechnique,2012,62(8):657 − 667.
[25]	GHEZZEHEI T A, OR D. Rheological properties of wet soils and clays under steady and oscillatory stresses[J]. Soil Science Society of America Journal,2001,65(3):624 − 637. doi: 10.2136/sssaj2001.653624x
[26]	BAYUK I O, AMMERMAN M, CHESNOKOV E M. Upscaling of elastic properties of anisotropic sedimentary rocks[J]. Geophysical Journal International,2008,172(2):842 − 860. doi: 10.1111/j.1365-246X.2007.03645.x
[27]	GAO Y, SUN D A, ZHOU A, et al. Predicting shear strength of unsaturated soils over wide suction range[J]. International Journal of Geomechanics,2020,20(2):04019175. doi: 10.1061/(ASCE)GM.1943-5622.0001555
[28]	SHENG D, ZHOU A, FREDLUND D G. Shear strength criteria for unsaturated soils[J]. Geotechnical and Geological Engineering,2011,29(2):145 − 159. doi: 10.1007/s10706-009-9276-x
[29]	CUNNINGHAM M R, RIDLEY A M, DINEEN K, et al. The mechanical behaviour of a reconstituted unsaturated silty clay[J]. Géotechnique,2003,53(2):183 − 194.
[30]	CAMPBELL G S, SHIOZAWA S. Prediction of hydraulic properties of soils using particle-size distribution and bulk density data [C]//VAN GENUCHTEN M T, LEIJ R J, LUND L J. International Workshop on Indirect Methods for Estimating the Hydraulic Properties of Unsaturated Soils. University of California, Riverside. 1992: 317 − 328.
[31]	LEBEAU M, KONRAD J-M. A new capillary and thin film flow model for predicting the hydraulic conductivity of unsaturated porous media[J]. Water Resources Research,2010,46(12):W12554.

施引文献

期刊类型引用(4)

1.	陈林志，罗军尧，王昆. 某地区非饱和粉土抗液化强度的试验研究. 科技通报. 2024(06): 67-71 . 百度学术
2.	金李皓琳，张敏，苑成旺，王占. 孔隙环境特征对高岭土超声波波速影响机理研究. 水文地质工程地质. 2023(05): 89-95 . 本站查看
3.	王康，畅俊斌，李晓科，朱文峰，卢啸，刘慧. 基于水文过程和应力应变耦合的陕北黄土滑坡复活机理分析——以延安二庄科滑坡为例. 中国地质灾害与防治学报. 2023(06): 47-56 . 百度学术
4.	刘佳婷，付昱凯，李同录，赵晨曦，李萍，侯晓坤，胡向阳. 黄土与其矿物颗粒表面水膜类型及其定量表征. 水文地质工程地质. 2022(06): 105-113 . 本站查看

其他类型引用(2)

资源附件(0)

访问统计

图(6)

计量

文章访问数: 1573
PDF下载数: 97
施引文献: 6

1. 非饱和粉土中孔隙水的微观赋存形态
1.1 收缩膜
1.2 吸附水
1.3 毛细水
2. 考虑孔隙水微观赋存形态的有效应力方程
2.1 非饱和粉土的扩展三相孔隙介质模型
2.2 有效应力方程的推导
3. 有效应力方程的验证
3.1 验证思路
3.2 验证结果
4. 结论

考虑孔隙水微观赋存形态的非饱和粉土有效应力方程及其验证

1. 工业装备结构分析国家重点实验室（大连理工大学）, 辽宁 大连 116024 2. 大连理工大学工程力学系, 辽宁 大连 116024 3. 青岛理工大学理学院, 山东 青岛 266033

作者简介: 曾立峰（1986-），男，博士，主要从事非饱和土力学方面的研究工作。E-mail: zenglifeng34@163.com

通讯作者: 邵龙潭（1963-），男，博士研究生，教授，博士生导师，主要从事土和孔隙介质力学基本理论研究、土工结构稳定分析、土工试验测试技术与仪器研发等方面的教学和科研。E-mail: shaolt@dlut.edu.cn