相变黏土在三轴压缩下的渗透和应力应变特性研究

杨家琦; 刘东海; 王泽帆

doi:10.11779/CJGE20221143

相变黏土在三轴压缩下的渗透和应力应变特性研究 English Version

天津大学水利工程仿真与安全国家重点实验室, 天津 300350

基金项目:

国家自然科学基金面上项目 52279136

国家自然科学基金面上项目 51979189

详细信息

作者简介:
杨家琦(1993—)，男，博士研究生，主要从事岩土新材料研究、大坝施工与服役性态耦合分析方面的研究工作。E-mail: demyang@163.com

通讯作者:
刘东海, E-mail: liudh@tju.edu.cn

中图分类号: TV641
计量
- 文章访问数: 434
- HTML全文浏览量: 38
- PDF下载量: 134
出版历程
- 收稿日期: 2022-09-18
- 网络出版日期: 2023-04-05
- 刊出日期: 2023-11-30

Permeability and strain-stress characteristics of phase-change clay under triaxial compression

State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300350, China

摘要

摘要: 掺混石蜡基相变材料的相变黏土为寒区黏土心墙冬季施工过程中防冻控温问题的解决提供了新途径，但关系到其筑坝服役性态的渗透和力学特性尚需深入研究。基于三轴渗透-压缩试验，研究了相变黏土在不同应力和变形状态下的力学和渗透特性。研究结果表明：相变黏土的应力-应变曲线为硬化型，剪胀现象不明显，其强度指标与其相同初始干密度条件下的纯黏土较为接近，但刚度指标存在差异，分析认为在相同初始干密度下，非极性石蜡基相变材料掺入会显著增加土体的脆性，但对其他力学指标影响较小；相变黏土的渗透系数总是小于相同条件下的纯黏土，表明疏水性的石蜡基相变材料能够显著提升土料的抗渗性；以孔隙比和偏应力水平为变量，建立的回归模型能够很好预测相变黏土在不同应力和变形状态下的渗透系数，可为相变黏土作为心墙筑坝所用的普通黏土的替代材料提供依据。
- 石蜡基相变材料 /
- 相变黏土 /
- 三轴压缩 /
- 力学特性 /
- 渗透特性
Abstract: The phase-change clay-included paraffin-based phase-change material (PPCM) provides a new way to solve the freeze-thaw problem during winter construction of clay-core walls in cold regions. However, the permeability and mechanical properties related to its service performance as a dam core wall need to be further studied. Accordingly, the permeability and mechanical properties of the phase-change clay under different stress and strain states are studied based on the triaxial permeability and compression tests. The results show that the stress-strain curves of the phase-change clay are strain hardening, and the dilatancy during loading is not obvious. The strength indices of the phase-change clay are close to those of the pure clay with the same initial dry density, but their stiffness indices are different. It is inferred that under the same initial dry density, the addition of non-polar PPCM significantly increases the brittleness of the soil but has small influences on the other mechanical indices. The hydraulic conductivity of the phase-change clay is always lower than that of the pure clay under the same conditions. This indicates that the hydrophobic PPCM can significantly improve the impermeability of soil. Moreover, the regression models established with the void ratio and deviatoric stress level as variables can well predict the hydraulic conductivities of the phase-change clay and the pure clay under different stress and strain states. These results may provide a basis for using the phase-change clay as a substitute of general clay for core-wall construction.
- paraffin-based phase-change material /
- phase-change clay /
- triaxial compression /
- mechanical property /
- permeability

HTML全文

0. 引言

近年来，海上风电迅速成为重点发展的能源形式之一。筒型基础凭借其成本低、工期短且对地基扰动相对较小的优势，已成为中国海上风电的重要基础形式之一。竖向承载力验算是海上风电基础稳定性校核的首要环节，筒型基础通常自重较大且长期承受海洋平台结构自重及一定程度的竖向荷载作用，因此了解其地基竖向承载特性对于实际工程设计而言十分重要。

针对筒型基础地基竖向承载力的研究较多，已有较多国内外的学者运用试验以及数值模拟等方式进行研究。Houlsby等^[1]和Byrne等^[2]在21世纪初期探讨了传统的吸力筒型基础作为海上风电基础的应用前景。Fu等^[3]、Gourvenec等^[4-5]、Chen等^[6]、刘梅梅等^[7]通过开展一系列室内模型试验及离心模型试验，得到了筒型基础的竖向承载模式、承载力p-s曲线、基础各部分载荷分担情况及筒内外土压力分布情况。Hung等^[8]、Vulpe^[9]、Mehravar等^[10]、Mana等^[11]基于大量有限元数值模拟，讨论了多种因素对地基竖向承载力的影响，得到了相关规律及半经验计算公式。然而由于筒型基础筒-土相互作用的特殊性，其地基竖向承载力的理论计算一直未能完全解决。

现有的筒型基础地基竖向承载力计算方法大多是利用极限分析上限法推导，但该方法得到计算结果往往高于地基承载力的真值，在工程应用中因偏于危险而无法推广。极限分析下限法从理论上更符合工程设计中“安全”概念，但传统方式构建容许应力场困难且具有主观性，导致下限法运用范围受限。直到1988年Sloan^[12]提出了适用于平面应变问题的有限元极限分析下限法，为极限分析下限法的应用提供了全新的思路，极大推动了其在岩土工程中的发展与应用。

由于筒型基础相互作用的特殊性，对于筒型基础的有限元极限分析下限法鲜有学者涉足，本文将筒土相互作用转化为筒土间摩擦的不等式约束条件，建立了筒型基础地基竖向承载力的下限解求解程序。利用离心模型试验结果验证了本方法的正确性，在此基础上分析了土体强度不均匀性、筒型基础长径比及筒外壁-土间摩擦对筒型基础地基竖向承载力的影响，提出了综合考虑以上因素的筒型基础地基竖向承载力计算公式。

1. 轴对称有限元极限分析下限法

本节在Sloan^[12]的有限元下限法求解格式基础上，结合轴对称问题相关理论，提出了针对筒型基础这一轴对称结构的单元平衡条件、应力间断面约束条件、土体边界条件、筒-土界面约束条件、土体屈服条件及相应的目标函数求解格式。

1.1 单元平衡约束条件

采用单元-结构面离散方式将结构物离散为有限个单元体，假设单元为刚性单元。相比于平面应变中以（ ${\sigma _x}$ ， ${\sigma _y}$ ， ${\tau _{xy}}$ ）为未知变量，轴对称有限元法引入环向应力 ${\sigma _\theta }$ ，将节点应力（ ${\sigma _r}$ ， ${\sigma _z}$ ， ${\tau _{rz}}$ ， ${\sigma _\theta }$ ）作为基本未知变量，模拟应力场的典型三角形单元如图 1所示。允许沿相邻三角形单元之间的所有界面发生静态允许的应力间断，每个节点都与唯一单元相关联。对于轴对称结构，各子午面上的应力状态均相同，假定当r-z平面（即子午面）上的应力满足极限平衡条件时，地基达到极限承载状态。

图 1 土体有限元离散

Figure 1. Finite element discretization of soil mass

下载: 全尺寸图片幻灯片

对于筒型基础地基竖向承载力分析模型，土体仅受重力作用，柱坐标系下单元分析示意图如图 2所示。

图 2 柱坐标系中单元体

Figure 2. Unit in a bucket coordinate system

下载: 全尺寸图片幻灯片

柱坐标系中单元体平衡方程可表示为

$\left.\begin{array}{l}\frac{\partial \sigma }{{\partial }_{r}}+\frac{\partial {\tau }_{rz}}{{\partial }_{z}}\text{+}\frac{{\sigma }_{r}-{\sigma }_{\theta }}{r}=\text{0 }\text{，}\\ \frac{\partial {\tau }_{rz}}{{\partial }_{r}}\text{+}\frac{\partial {\sigma }_{z}}{{\partial }_{z}}\text{+}\frac{{\tau }_{rz}}{r}=\gamma \text{ }。\end{array}\right\}$

(1)

假设土单元为线性单元，单元中3个节点满足平衡方程即代表单元平衡。在此基础上，将式（1）平衡方程转化为离散形式产生如下等式约束：

${\left[ {A_{{\text{equil}}}^{\text{e}}} \right]_{2 \times 12}}{\left\{ {{\sigma ^{\text{e}}}} \right\}_{12 \times 1}} = {\left\{ {b_{{\text{equil}}}^{\text{e}}} \right\}_{2 \times 1}} 。$

(2)

式中： ${\left[{A}_{\text{equil}}^{\text{e}}\right]}_{\text{2}\times \text{12}}\text{ }$

${\text{ = }}\left[ \begin{gathered} \frac{{{z_{23}}}}{{2A}} + \frac{1}{{3\bar r}}{\text{ 0 }}\frac{{{r_{32}}}}{{2A}}{\text{ }} - \frac{1}{{3\bar r}}{\text{ }} \hfill \\ {\text{ }}0{\text{ }}\frac{{{r_{32}}}}{{2A}}{\text{ }}\frac{{{z_{23}}}}{{2A}} + \frac{1}{{3\bar r}}{\text{ }}0 \hfill \\ \end{gathered} \right.$

$\left.\begin{array}{cc}\frac{{z}_{31}}{2A}+\frac{1}{3\overline{r}}& 0\\ 0& \frac{{r}_{13}}{2A}\end{array}\text{ }\begin{array}{cc}\frac{{r}_{13}}{2A}& -\frac{1}{3\overline{r}}\\ \frac{{z}_{31}}{2A}+\frac{1}{3\overline{r}}& \text{0}\end{array}\text{ }\begin{array}{cc}\frac{{z}_{12}}{2A}+\frac{1}{3\overline{r}}& \text{0}\\ 0& \frac{{r}_{21}}{2A}\end{array}\text{ }\begin{array}{cc}\frac{{r}_{21}}{2A}& -\frac{1}{3\overline{r}}\\ \frac{{z}_{12}}{2A}+\frac{1}{3\overline{r}}& 0\end{array}\right]\text{，} \\ {\left\{ {{\sigma ^{\text{e}}}} \right\}^{\text{T}}} = \{ \sigma _{r{\text{,1}}}^{\text{e}}{\text{ }}\sigma _{z{\text{,1}}}^{\text{e}}{\text{ }}\tau _{rz{\text{,1}}}^{\text{e}}{\text{ }}\sigma _{\theta {\text{,1}}}^{\text{e}}{\text{ }}\sigma _{r{\text{,2}}}^{\text{e}}{\text{ }}\sigma _{z{\text{,2}}}^{\text{e}}{\text{ }}\tau _{rz{\text{,2}}}^{\text{e}}{\text{ }}\sigma _{\theta {\text{,2}}}^{\text{e}}\sigma _{r{\text{,3}}}^{\text{e}}{\text{ }}\sigma _{z{\text{,3}}}^{\text{e}} \\ {\tau }_{rz\text{,3}}^{\text{e}} {\sigma }_{\theta \text{,3}}^{\text{e}}{\}}_{1\times 12} \text{，}{\left\{ {b_{{\text{equil}}}^{\text{e}}} \right\}^{\text{T}}} = {\left\{ {0{\text{ }}\gamma } \right\}_{1 \times 2}}。$

式中： ${r_{ij}} = {r_i} - {r_j},{z_{ij}} = {z_i} - {z_j}$ ； $\bar r = ({r_1} + {r_2} + {r_3})/3$ 。

1.2 应力间断面约束条件

将相邻独立单元的共同面定义为应力间断面，如图 3所示。在应力允许的情况下，控制各节点的正应力和剪应力对应相等。通过一对共轭节点（1, 2）和（3, 4）定义了一个典型的应力间断线。间断线与r轴正向形成的角度ω，定义相对于r轴正向逆时针方向转动为正方向。

图 3 应力间断面

Figure 3. Stress discontinuity

下载: 全尺寸图片幻灯片

为了满足平衡条件，作用于该路径上的正应力（ $\sigma _{\text{n}}^{}$ ）和剪应力（ $\tau _{\text{n}}^{}$ ）的值应始终保持不变，即

${\sigma }_{\text{n,1}}^{\text{a}}={\sigma }_{\text{n,2}}^{\text{b}} \text{，}{\sigma }_{\text{n,3}}^{\text{a}}={\sigma }_{\text{n,4}}^{\text{b}}\text{，}{\tau }_{\text{n,1}}^{\text{a}}={\tau }_{\text{n,2}}^{\text{b}}\text{，}{\tau }_{\text{n,3}}^{\text{a}}={\tau }_{\text{n,4}}^{\text{b}} 。$

(3)

通过r-z平面某点的应力状态( $\sigma _r^{}$ ， $\sigma _z^{}$ 和 $\tau _{rz}^{}$ )，利用应力莫尔圆得到正应力和剪应力的值， $\sigma _{\text{n}}^{}$ 和 $\tau _{\text{n}}^{}$ 的最终表达式如下：

${\sigma _{\text{n}}} = {\sin ^2}\omega {\sigma _r} + {\cos ^2}\omega {\sigma _z} - \sin 2\omega {\tau _{rz}}\text{，}$

(4)

${\tau _{\text{n}}} = - \frac{1}{2}\sin 2\omega {\sigma _r} + \frac{1}{2}\sin 2\omega {\sigma _z} + \cos 2\omega {\tau _{rz}}。$

(5)

将式（3）的约束条件转化为离散形式：

${\left[ {A_{{\text{stat}}}^{{\text{dc}}}} \right]_{4 \times 16}}{\left\{ {{\sigma ^{{\text{dc}}}}} \right\}_{16 \times 1}} = {\left\{ {b_{{\text{stat}}}^{{\text{dc}}}} \right\}_{4 \times 1}} 。$

(6)

式中： ${\left[ {A_{{\text{stat}}}^{{\text{dc}}}} \right]_{{\text{4}} \times {\text{16}}}} = {\left[ {\begin{array}{*{20}{c}} T&{ - T} \\ 0&0 \end{array}{\text{ }}\begin{array}{*{20}{c}} {\text{0}}&{\text{0}} \\ T&{ - T} \end{array}} \right]_{{\text{4}} \times {\text{16}}}}$

${\left\{ {{\sigma ^{{\text{dc}}}}} \right\}^{\text{T}}}{\text{ = }}\left\{ {\sigma _{r{\text{,1}}}^{\text{a}}{\text{ }}\sigma _{z{\text{,1}}}^{\text{a}}{\text{ }}\tau _{rz{\text{,1}}}^{\text{a}}{\text{ }}\sigma _{\theta {\text{,1}}}^{\text{a}}{\text{ }}\sigma _{r{\text{,2}}}^{\text{b}}{\text{ }}\sigma _{z{\text{,2}}}^{\text{b}}{\text{ }}\tau _{rz{\text{,2}}}^{\text{b}}{\text{ }}\sigma _{\theta {\text{,2}}}^{\text{b}}} \right.$

${{\sigma }_{r\text{,3}}^{\text{a}}\text{ }{\sigma }_{z\text{,3}}^{\text{a}}\text{ }{\tau }_{rz\text{,3}}^{\text{a}}\text{ }{\sigma }_{\theta \text{,3}}^{\text{a}}\text{ }{\sigma }_{r\text{,4}}^{\text{b}}\text{ }{\sigma }_{z\text{,4}}^{\text{b}}\text{ }{\tau }_{rz\text{,4}}^{\text{b}}\text{ }{\sigma }_{\theta \text{,4}}^{\text{b}}\text{ }\}}_{\text{1}\times \text{16}} \text{，}$

${\left\{ {b_{{\text{stat}}}^{{\text{dc}}}} \right\}^{\text{T}}} = \left\{ {{\text{0 0 0 0}}} \right\}1 \times 4。$

1.3 土体边界条件

在任何土体边界l（如图 4所示）的边界条件可以用以下方式来定义：

${\sigma }_{\text{n,1}}^{\text{l}}={q}_{1}\text{，}{\sigma }_{\text{n,2}}^{\text{l}}={q}_{2}\text{，}{\tau }_{\text{n,1}}^{\text{l}}={t}_{1}\text{，}{\tau }_{\text{n,1}}^{\text{l}}={t}_{2} 。$

(7)

图 4 应力边界条件

Figure 4. Stress boundary conditions

下载: 全尺寸图片幻灯片

将任意土体边界l中两节点所对应公式（7）的4个等式约束转化为离散形式如下：

${\left[ {A_{{\text{bound}}}^{\text{l}}} \right]_{4 \times 8}}{\left\{ {{\sigma ^{\text{l}}}} \right\}_{8 \times 1}} = {\left\{ {b_{{\text{bound}}}^{\text{l}}} \right\}_{4 \times 1}}。$

(8)

式中： ${\left[{A}_{\text{bound}}^{\text{l}}\right]}_{4\times 8}={\left[\begin{array}{c}M 0\\ 0\text{ }M\end{array}\right]}_{4\times 8}$ ，

${\left[M\text{ }\right]}_{\text{2}\times \text{4}}=\left[\begin{array}{cc}{\mathrm{sin}}^{2}\alpha & {\mathrm{cos}}^{2}\alpha \\ -\frac{1}{2}\mathrm{sin}2\alpha & \frac{1}{2}\mathrm{sin}2\alpha \end{array}\text{ }\begin{array}{cc}-\mathrm{sin}2\alpha & 0\\ \mathrm{cos}2\alpha & 0\end{array}\right] \text{，}$

${\left\{{\sigma }^{\text{l}}\right\}}^{\text{T}}={\left\{{\sigma }_{r\text{,1}}^{\text{l}} {\sigma }_{z\text{,1}}^{\text{l}} {\tau }_{rz\text{,1}}^{\text{l}} {\sigma }_{\theta \text{,1}}^{\text{l}} {\sigma }_{r\text{,2}}^{\text{l}} {\sigma }_{z\text{,2}}^{\text{l}} {\tau }_{rz\text{,2}}^{\text{l}} {\sigma }_{\theta \text{,2}}^{\text{l}}\right\}}_{1\times 8} \text{，}$

${\left\{{b}_{\text{bound}}^{\text{l}}\right\}}^{T}={\left\{{q}_{1} {t}_{1} {q}_{2} {t}_{2} \right\}}_{1\times 4} 。$

式中， $\alpha$ 为边界l与r轴之间夹角，定义相对于r轴逆时针方向转动为正方向。

1.4 筒-土界面约束条件

在竖向荷载作用下，与筒壁内外侧接触的土体将产生沿筒壁向上的摩擦力，与筒顶及筒壁端部接触的土体将产生垂直端部向上的土抗力。因此，同时对筒顶-土接触面、筒端-土接触面及筒壁-土接触面的进行边界约束。

定义发生剪切破坏时沿筒顶-土界面、筒端-土界面和筒壁-土界面切向抗剪强度均满足以下条件：

$\left| {{\tau _s}} \right| \leqslant {\alpha _{\text{s}}}{s_{{\text{us}}}}\text{；}\left| {{\tau _{\text{b}}}} \right| \leqslant {\alpha _{\text{b}}}{s_{{\text{ub}}}}\text{；}\left| {{\tau _{{\text{f}}i}}} \right| \leqslant {\alpha _{{\text{f}}i}}{s_{{\text{u}}i}}。$

(9)

式中： ${\alpha _{\text{s}}}$ ， ${\alpha _{\text{b}}}$ 和 ${\alpha _{{\text{f}}i}}$ 分别代表筒顶-土界面、筒端-土界面和筒壁-土界面接触强度的折减系数； ${s_{{\text{us}}}}$ 和 ${s_{{\text{ub}}}}$ 分别为筒顶和筒端处的土体不排水抗剪强度， ${s_{{\text{u}}i}}$ 为筒壁-土界面上任一应力点所在深度处土体的不排水抗剪强度。

式（9）可离散转化为

${\left[ {{A_{{\text{int}}}}} \right]_{{\text{2}} \times {\text{4}}}}{\left\{ \sigma \right\}_{{\text{4}} \times {\text{1}}}} \leqslant {\left\{ {{b_{{\text{int}}}}} \right\}_{{\text{2}} \times {\text{1}}}}。$

(10)

式中： ${\left[ {{A_{{\text{int}}}}} \right]_{{\text{2}} \times {\text{4}}}}{\text{ = }}\left[ {\begin{array}{*{20}{c}} {\text{0}}&{\text{0}} \\ 0&0 \end{array}{\text{ }}\begin{array}{*{20}{c}} { - {\text{1}}}&0 \\ 1&0 \end{array}} \right]$ ，

${\left\{ \sigma \right\}^{\text{T}}}{\text{ = }}\left\{ {{\sigma _{r{\text{,}}i}}{\text{ }}{\sigma _{z{\text{,}}i}}{\text{ }}{\tau _{rz{\text{,}}i}}{\text{ }}{\sigma _{\theta {\text{,}}i}}} \right\}1 \times 4 \text{，}$

${\left\{ {{b_{{\text{int}}}}} \right\}^{\text{T}}}{\text{ = }}{\left\{ {{\alpha _{{\text{su}}}}{\text{ }}{\alpha _{{\text{su}}}}} \right\}_{{\text{1}} \times {\text{2}}}} 。$

图 5展示了筒顶-土界面（EP区域）、筒端-土界面（OQ区域）、筒壁-土界面（PQ区域）的剪应力约束形式，及各筒-土界面的剪应力正方向。

图 5 筒-土边界示意图

Figure 5. Boundary between bucket and soil

下载: 全尺寸图片幻灯片

1.5 土体屈服条件

根据Harr-VonKarman假设， ${\sigma _\theta }$ 为在r-z平面上的 ${\sigma _{\text{3}}}$ （最小主应力），即假定环向应力不大于大主应来求解轴对称稳定性问题^[13]。 ${\sigma _\theta }$ 的取值范围在是通过绘制两个莫尔圆确定，破坏面土体受力莫尔圆如所示。为了确保环向应力 ${\sigma _\theta }$ 的值始终接近于 ${\sigma _{\text{3}}}$ 的值，确定以下3个不等式约束：

图 6 破坏面土体受力莫尔圆

Figure 6. Mohr circle of failure surface

下载: 全尺寸图片幻灯片

$\left.\begin{array}{l}{\sigma }_{\theta \text{,}i}\ge {\sigma }_{r\text{,}i}\text{ }\text{，}\\ {\sigma }_{\theta \text{,}i}\ge {\sigma }_{z\text{,}i}\text{ }\text{，}\\ {\sigma }_{\theta \text{,}i}\le {\sigma }_{\text{3f,}i}\text{ }。\end{array}\right\}$

(11)

式中： ${\sigma _{\theta {\text{,}}i}}$ ， ${\sigma _{r{\text{,}}i}}$ 和 ${\sigma _{z{\text{,}}i}}$ 为与节点i相关的应力， ${\sigma _{{\text{3f}}}}$ 为屈服时的小主应力。这些约束可以用矩阵的形式表示为

${\left[ {A_{r\theta }^i} \right]_{{\text{3}} \times {\text{4}}}}{\left\{ {{\sigma ^i}} \right\}_{{\text{4}} \times {\text{1}}}} \leqslant {\left\{ {b_{r\theta }^i} \right\}_{{\text{3}} \times {\text{1}}}}。$

(12)

式中：

${\left[{A}_{r\theta }^{i}\right]}_{\text{3}\times \text{4}}=\left[\begin{array}{cc}1& 0\\ \begin{array}{l}\text{ }0\\ -{\scriptstyle \frac{1}{2}}\text{(}1-\mathrm{sin}\varphi \text{)}\end{array}& \begin{array}{l}\text{ 1}\\ \text{ }-{\scriptstyle \frac{1}{2}}(1-\mathrm{sin}\varphi )\end{array}\end{array}\text{ }\begin{array}{cc}\text{0}& -1\\ \begin{array}{l}0\\ \text{0}\end{array}& \begin{array}{l}-1\\ \text{ 1}\end{array}\end{array}\right]，$

${\left\{{\sigma }^{i}\right\}}^{\text{T}}={\left\{{\sigma }_{r\text{,}i}\text{ }{\sigma }_{z\text{,}i}\text{ }{\tau }_{rz\text{,}i}\text{ }{\sigma }_{\theta \text{,}i} \right\}}_{1\times 4} \text{，}$

${\left\{ {b_{r\theta }^i} \right\}^{\text{T}}} = {\left\{ {0{\text{ }}0{\text{ }}c\cos \varphi } \right\}_{1 \times 3}}。$

在下限解中用内接多边形来逼近Mohr-Coulomb准则，从而方便求解计算的线性规划问题。为保证任一点上的应力状态均满足屈服条件，需在r-z平面上保证满足以下线性化的不等式条件：

${A_k}{\sigma _r} + {B_k}{\sigma _z} + {C_k}{\tau _{rz}} \leqslant \varepsilon {\text{ (}}k = 1,2, \cdots ,p) 。$

(13)

式中： ${A_k} = \cos (2{\rm{\mathsf{π}}}k/p) + \sin \varphi \cos ({\rm{\mathsf{π}}}/p)$ ， ${B_k} = \sin \varphi \cdot$ $\mathrm{cos}({\rm{\mathsf{π}}}/p)-\mathrm{cos}(2{\rm{\mathsf{π}}}k/p)，{C}_{\text{k}}=2\mathrm{sin}(2{\rm{\mathsf{π}}}k/p)，\varepsilon =2c\mathrm{cos}\varphi \cdot$ $\cos ({\rm{\mathsf{π}}}/p)$ 。其中，c为节点i处土体的不排水抗剪强度。

由线性化屈服条件对节点i的4个应力张量施加的不等式约束（数目p）可表示为

${\left[ {A_{{\text{yield}}}^i} \right]_{p \times 4}}{\left\{ {{\sigma ^i}} \right\}_{4 \times 1}} \leqslant {\left\{ {b_{{\text{yield}}}^i} \right\}_{p \times 1}}。$

(14)

式中，

${\left[{A}_{\text{yield}}^{i}\right]}_{p\times 4}=\left[\begin{array}{c}{A}_{1} {B}_{1} {C}_{1} 0\\ {A}_{2} {B}_{2} {C}_{2} 0\\ .\text{ }.\text{ }.\\ {A}_{k} {B}_{k} {C}_{k} 0\\ .\text{ }.\text{ }.\\ {A}_{p} {B}_{p} {C}_{p} 0\end{array} \right] 。$

1.6 最优化函数

通过对筒型基础顶盖部分圆形区域、筒端部分圆环区域的垂直法向应力和筒壁区域切向应力积分。地基竖向承载力（目标函数）可通过以下表达式来定义：

$Q = {Q_{\text{s}}}{\text{ + }}{Q_{\text{f}}}{\text{ + }}{Q_{\text{b}}}{\text{ = }}2{\rm{\mathsf{π}}}\int\limits_{{L_{\text{s}}}} {{\sigma _{\text{s}}}} r{\text{d}}s{\text{ + }}2{\rm{\mathsf{π}}}\int\limits_{{L_{\text{f}}}} {{\tau _{\text{f}}}} z{\text{d}}s{\text{ + }}2{\rm{\mathsf{π}}}\int\limits_{{L_{\text{b}}}} {{\sigma _{\text{b}}}} r{\text{d}}s 。$

(15)

式中：Q为地基竖向承载力，Q_s，Q_f，Q_b分别为筒顶、筒壁及筒端承载部分。 ${\sigma _{\text{s}}}$ 和 ${\sigma _{\text{b}}}$ 分别为作用在筒顶-土和筒端-土边界面上的法向应力， ${\tau _{\text{f}}}$ 为作用在筒壁-土边界面上的切向抗剪强度。由于假设应力在每个单元中为线性变化，筒型基础各部分产生的极限承载力的分量分别为

${Q_{\text{s}}}{\text{ = }}2{\rm{\mathsf{π}}}\frac{{L_{\text{s}}^i}}{2}({\sigma _{{\text{s,1}}}}{\text{ + }}{\sigma _{{\text{s,2}}}})r_{\text{s}}^i \text{，}$

(16)

${Q_{\text{f}}}{\text{ = }}2{\rm{\mathsf{π}}}\frac{{L_{\text{f}}^i}}{2}({\tau _{{\text{f,1}}}}{\text{ + }}{\tau _{{\text{f,2}}}})z_{\text{f}}^i \text{，}$

(17)

${Q_{\text{b}}}{\text{ = }}2{\rm{\mathsf{π}}}\frac{{L_{\text{b}}^i}}{2}({\sigma _{{\text{b,1}}}}{\text{ + }}{\sigma _{{\text{b,2}}}})r_{\text{b}}^i 。$

(18)

式中： $L_{}^i$ 为节点i所在单元边界的长度； $r_{}^i$ 为边界 $L_{}^i$ 面的平均半径； $L_{\text{f}}^i$ 为筒壁节点上i所在单元边界的长度； $z_{\text{f}}^i$ 为边界 $L_{\text{f}}^i$ 面的平均长度。

将方程（1），（3）和（7）进行整合，得到全局等式约束矩阵形式如下：

$\left[ {{A_1}} \right] = \sum\limits_{{\text{e = }}1}^E {\left[ {A_{_{{\text{equil}}}}^{\text{e}}} \right]} {\text{ + }}\sum\limits_{{\text{dc}}}^{{D_{\text{c}}}} {\left[ {A_{_{{\text{stat}}}}^{{\text{dc}}}} \right]} {\text{ + }}\sum\limits_{{\text{l = }}1}^L {\left[ {A_{{\text{bound}}}^{\text{l}}} \right]} 。$

(19)

将方程（9），（10），（14）进行整合，得到全局不等式约束矩阵成如下形式：

$\left[ {{A_2}} \right] = \sum\limits_{i{\text{ = }}1}^N {\left[ {A_{_{{\text{yield}}}}^i} \right]} {\text{ + }}\sum\limits_{i{\text{ = }}1}^N {\left[ {A_{{\text{s }{\rm{\mathsf{θ}}}}}^i} \right]} {\text{ + }}\sum\limits_{{\text{ns = }}1}^{{N_{\text{s}}}} {\left[ {A_{{\text{ints}}}^{{\text{ns}}}} \right]} {\text{ + }}\sum\limits_{{\text{nf = }}1}^{{N_{\text{f}}}} {\left[ {A_{{\text{intf}}}^{{\text{nf}}}} \right]} +$

$\sum\limits_{{\text{nb = }}1}^{{N_{\text{b}}}} {\left[ {A_{{\text{intb}}}^{{\text{nb}}}} \right]} 。$

(20)

与等式及不等式约束条件对应的结果矩阵 ${b_1}$ 和 ${b_2}$ 形式如下：

$\left[ {{b_1}} \right] = \sum\limits_{{\text{e = }}1}^E {\left\{ {b_{_{{\text{equil}}}}^{\text{e}}} \right\}} {\text{ + }}\sum\limits_{{\text{dc}}}^{{D_{\text{c}}}} {\left\{ {b_{_{{\text{stat}}}}^{{\text{dc}}}} \right\}} {\text{ + }}\sum\limits_{{\text{l = }}1}^L {\left\{ {b_{{\text{bound}}}^{\text{l}}} \right\}} \text{，}$

(21)

$\left[ {{b_2}} \right] = \sum\limits_{i{\text{ = }}1}^N {\left\{ {b_{_{{\text{yield}}}}^i} \right\}} {\text{ + }}\sum\limits_{i{\text{ = }}1}^N {\left\{ {b_{{\text{s }{\rm{\mathsf{θ}}}}}^i} \right\}} {\text{ + }}\sum\limits_{{\text{ns = }}1}^{{N_{\text{s}}}} {\left\{ {b_{{\text{ints}}}^{{\text{ns}}}} \right\}{\text{ + }}\sum\limits_{{\text{nf = }}1}^{{N_{\text{f}}}} {\left\{ {b_{{\text{intf}}}^{{\text{nf}}}} \right\}} } \sum\limits_{{\text{nb = }}1}^{{N_{\text{b}}}} {\left\{ {b_{{\text{intb}}}^{{\text{nb}}}} \right\}} 。$

(22)

调用线性最优化函数对目标函数求解的方法，陈广思^[14]和Vishwas^[15]等已给出介绍，不在此赘述。

2. 离心试验验证

2.1 试验准备

为验证筒型基础地基承载力下限解的准确性，开展了饱和黏土地基离心模型试验。试验在天津大学TLJ-100A型土工离心机完成。模型箱内部尺寸为880 mm×595 mm×400 mm。将高岭土粉按照2倍液限的含水率进行真空搅拌，为模拟天然地基的自然固结过程，将土体在100g离心条件下固结10×12 h，固结完成土体厚度约为280 mm。

模型筒外径100 mm，裙板长度95 mm，壁厚10 mm，顶盖厚5 mm，筒顶处预留连接螺杆，与量程为8 kN的荷载传感器相连，模型示意图如图 7所示。

图 7 筒型基础模型示意图

Figure 7. Model for bucket foundation

下载: 全尺寸图片幻灯片

2.2 试验方案

试验前取样测得软黏土的有效重度 $\gamma '$ =4.83 kN/m³，采用十字板对模型槽内土体的不排水抗剪强度进行了测量，结果如图 8所示，模型土泥面处土体强度s_um=6 kPa且土体不排水抗剪强度以k=0.426 kPa/m随深度线性增加。

图 8 试验土体的不排水强度

Figure 8. Undrained shear strengths of soil

下载: 全尺寸图片幻灯片

试验布置方式如图 9所示，试验采用分级加载，具体加载方式见图 10。

图 9 筒型基础及传感器整体布置

Figure 9. Overall layout of bucket foundation and sensor

下载: 全尺寸图片幻灯片

图 10 荷载时间曲线

Figure 10. Load-time curves

下载: 全尺寸图片幻灯片

2.3 试验结果对比

离心机模型试验得到的p-s曲线如图 11所示，在整个沉降过程中未出现明显的拐点。利用初始直线段与末尾直线段的切线交点对应的q_u值作为地基极限承载力，如图中黑色虚线所示，此时对应的竖向承载力q_u约为93 kPa。

图 11 离心试验p-s曲线

Figure 11. Load-displacement curves

下载: 全尺寸图片幻灯片

根据试验参数建立有限元极限分析下限解分析模型，将下限解结果和试验结果进行对比见表 1。下限解结果与模型试验的结果吻合良好，下限解结果小于API规范结果，且计算结果误差均在可接受范围，验证了极限分析下限解的正确性。

表 1 结果对比

Table 1. Comparison of results

试验结果/kPa	API/kPa	下限解/kPa
93	104.84	92.02

下载: 导出CSV

| 显示表格

3. 筒型基础承载特性分析

3.1 下限分析模型的建立

以海上风电筒型基础作为主要研究对象，基础外径D=10 m，侧壁壁厚t_s=0.01 m。研究不同筒裙长度d筒型基础的承载特性，主要针对长径比d/D=0.25，0.5，0.75，1，1.5和2的基础进行分析。为降低边界效应，对土体模型的尺寸进行灵敏度分析，确定土体直径取D_soil=2.5D，土体高度d_soil=3D+d，网格加密形式及模型示意图见图 12。

图 12 筒型基础下限分析网格示意图

Figure 12. Meshes for lower bound analysis of bucke foundation

下载: 全尺寸图片幻灯片

利用基础横截面积S和筒端处土体不排水土强度s_u0对地基竖向承载力V_ult进行了无量纲化处理得到N_cV，即

$N_{\rm cV}=V_{\rm ult}/S·s_{\rm u0} 。$

(23)

3.2 土体强度变化对竖向承载力的影响

大量现场原位测试结果表明，对于大部分的正常固结土及超固结土而言，其不排水抗剪强度随埋深的增加而呈近似的线性增长趋势，假设黏土的不排水剪切强度随深度呈线性变化：s_uz=s_um+kz。其中k为土体强度随深度的增长系数，s_um为泥面处土体的不排水抗剪强度，s_uz为地基深度为z处土体的不排水抗剪强度。

为得出均质土不排水抗剪强度s_u对筒型基础地基竖向极限承载力的影响，针对s_u=10~40 kPa的均质黏土地基进行计算，分析模型中筒-土界面设置为完全粗糙接触，计算结果如图 13所示。

图 13 土体不排水强度对筒型基础竖向承载力的影响

Figure 13. Effects of s_u on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

由图 13可知，当d/D＜1时，N_cV增长曲线随土体强度增长近似水平，故认为在此范围内筒型基础承载力与s_u无关；而当d/D＞1时，N_cV随强度增长有较为轻微下降趋势。总体来看，对于长径比较大的筒型基础N_cV虽随s_u增长而略有降低，但结合现有筒型基础工程实际情况，此部分降低可忽略。

为得出非均质土体沿深度方向强度不均匀性对筒型基础竖向承载力的影响，针对s_um=4 kPa及s_um=40 kPa且k均为0~2的地基进行了下限解求解，计算结果如所示。引入不同土体强度线性增长梯度下筒型基础地基竖向极限承载力与均质土体（k=0）条件下承载力的比值 ${\eta _k}$ 来评估土体k增长对筒型基础地基极限承载力的增益。 ${\eta _k}$ 定义如下：

${\eta _k} = \frac{{{V_{{\text{ult}}}}(k)}}{{{V_{{\text{ult}}}}(k = 0)}}。$

(24)

由图 14可知，当s_um=4时，随着土体强度增长梯度k由0增长至2， ${\eta _k}$ 增长约120%~410%；当s_um=40时，随着土体强度增长梯度k由0增长至2， ${\eta _k}$ 增长约2%~35%，不同长径比筒型基础V_ult皆随k增长而不同程度的线性增加。单独考虑土体参量s_um或k时，无法通过N_cV来表征竖向承载力特性。为使承载力系数具有普适性，通过无量纲参数κ来定量表征土体强度沿深度方向不均匀性，定义κ=kD/s_um，计算时k和D固定不变，通过改变s_um来改变κ。针对κ=0.5，2，3，5，10，20的非均质饱和黏土地基进行讨论，计算结果如图 15所示。

图 14 k对筒型基础竖向承载力影响

Figure 14. Effects of k on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

图 15 κ对筒型基础竖向承载力的影响

Figure 15. Effects of κ on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

如图 15所示，N_cV随κ的增加而减小。因为κ越高意味沿裙板上平均强度越低，筒壁侧摩阻发挥较低，所以其承载力增长率比基础端部土体强度的增长率低，从而使N_cV降低。

3.3 长径比对筒型基础竖向极限承载力的影响

对均质土中长径比为0.5~2的筒型基础地基承载力进行计算，得到其地基承载力随长径比增长的变化，结果如图 16所示。

图 16 不同长径比的筒型基础在均质土中的承载特性

Figure 16. Vertical bearing capacities of bucket foundation with different length diameter ratios in homogeneous soils

下载: 全尺寸图片幻灯片

由图 16可知，均质土中筒型基础承载力系数N_cV增率随d/D增大呈现非线性增长趋势，且规律随土体强度变化并未发生改变。将均质土中筒型基础地基竖向承载力系数N_cV随土体不排水抗剪强度s_u及长径比d/D的变化规律拟合为

${N_{{\text{cV}}}}{\text{(}}\kappa {\text{ = 0) = 5}}{\text{.91 + 8}}{\text{.41}}\left( {\frac{d}{D}} \right) - {\text{1}}{\text{.12}}{\left( {\frac{d}{D}} \right)^2} 。$

(25)

对于强度沿深度线性增加黏土地基，选取κ=0.5，2，3，5，10，20对不同长径比筒型基础地基竖向承载力进行分析，得到N_cV随长径比增长的变化规律，结果如图 17所示。

图 17 不同长径比筒型基础在非均质土中的承载特性

Figure 17. Vertical bearing capacities of bucket foundation with different length diameter ratios in nonhomogenous soil

下载: 全尺寸图片幻灯片

由图 17可知，对于强度沿深度线性增加的土体，N_cV随d/D的增长而增大。将筒型基础地基竖向承载力系数N_cV随土体强度沿深度不均匀性κ及长径比d/D的变化关系可采用公式（26）进行拟合，由图 17可知，拟合结果良好。

${N_{{\text{cV}}}}{\text{(}}\kappa {\text{) = 3}}{\text{.22 + 2}}{\text{.02}}\left( {\frac{d}{D}} \right){\text{ + 6}}{\text{.05}}{{\text{e}}^{ - 0.45\kappa }} 。$

(26)

为探究筒型基础不同部分对地基承载力的贡献率随长径比的变化，分别计算不同土体强度（s_um=10 kPa，s_um=40 kPa）及不同土体强度线性增长梯度k（k=0，2）条件下不同长径比的筒型基础地基竖向承载力，提取极限状态下筒壁内、外侧摩阻力，筒顶抗力及筒端处抗力。考虑到荷载情况仅为竖向荷载，筒内壁与土体接触良好，认为土体抗力为筒顶抗力和筒内壁摩阻之和，分析各部分在地基竖向承载力中占比，结果绘制于图 18。

图 18 筒基各部分承担荷载比例随长径比的变化

Figure 18. Variation of loading ratio in different parts of bucket foundations with d/D

下载: 全尺寸图片幻灯片

由图 18可知，在不同土体强度条件下，随着基础长径比的增加，筒端抗力贡献率稳定在10%以下，土体抗力的贡献率逐渐降低，筒壁外侧摩阻逐渐增大。分析原因认为，随d/D的增加，筒土接触面积逐步增大，外侧摩阻增长显著，筒端处虽应力较高，但由于筒壁较薄，故对承载力贡献率不高。在相同长径比条件下，随κ增大，土体抗力的贡献率逐步增加。由于κ越高意味着筒端下方参与承载的土体强度越高，筒端以上和筒壁接触的土体强度相对较弱，无法有效地为筒壁-土接触面提供剪切力，故而筒壁外侧摩阻的增长相对较低。

3.4 ${\alpha _{{\text{f}}i}}$ 沿筒壁分布规律

为探究土体达到极限承载状态时筒外壁-土界面接触强度折减程度沿深度方向的实际分布规律，以均质土（s_um=10 kPa，k=0，即κ=0）和强度随深度线性增长黏土（s_um=10 kPa，k=2，即κ=2）地基为例，计算筒外壁-土间完全粗糙接触条件下的地基竖向承载力。提取极限竖向荷载状态下筒外壁-土界面切向抗剪强度 ${\tau _{{\text{f}}i}}$ 及应力点处不排水抗剪强度s_ui，参考式（9）进行归一化处理，即 ${\alpha _{{\text{f}}i}}$ = ${\tau _{{\text{f}}i}}$ /s_ui，计算结果如图 19所示。

图 19

${\alpha _{{\text{f}}i}}$ 沿深度分布规律

Figure 19. Distribution laws of

${\alpha _{{\text{f}}i}}$ along depth

下载: 全尺寸图片幻灯片

图 19中纵坐标w/D为应力点距筒顶的距离w与筒径D之比。分析可知，极限状态时，尽管筒外壁-土间接触关系设定为完全粗糙， ${\alpha _{{\text{f}}i}}$ 沿深度分布也并非定值。

均质土中，筒壁上部w/D=0~0.5范围内，d/D＜1的筒型基础 ${\alpha _{{\text{f}}i}}$ 随深度增加而增大，d/D＞1的筒型基础筒壁侧摩阻力不变且 ${\alpha _{{\text{f}}i}}$ 近似等于1。而当w/D＞0.5，随着应力点靠近筒端 ${\alpha _{{\text{f}}i}}$ 逐渐减小。

对于强度沿深度线性增长的黏土地基， ${\alpha _{{\text{f}}i}}$ 随深度增长而降低，w/D=0~0.5范围内降低速率较小，而距筒顶0.5D后 ${\alpha _{{\text{f}}i}}$ 开始加速降低。对于同一深度的筒壁上各应力点，强度线性增长黏土中 ${\alpha _{{\text{f}}i}}$ 均小于均质土，因此，对于黏土地基，筒壁-土界面处土体不排水抗剪强度越高， ${\alpha _{{\text{f}}i}}$ 越小。

3.5 ${\alpha _{{\text{f}}i}}$ 对筒型基础地基竖向承载力的影响

以均质土（κ=0）和强度沿深度线性增长的黏土地基（κ=20）为例，计算不同 ${\alpha _{{\text{f}}i}}$ 的筒型基础地基竖向承载力，结果绘制于图 20。

图 20 α_fi对筒型基础竖向承载力的影响

Figure 20. Effects of α_fi on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

由可知，对于均质黏土和强度随深度线性增长的黏土地基，N_cV均随 ${\alpha _{{\text{f}}i}}$ 的增大而增大，并且长径比大的筒型基础N_cV增率更为显著。其原因在于长径比大的筒型基础筒外壁-土接触面积更大，筒外壁侧摩阻对承载力的贡献也更大。

综合考虑土体不均匀性κ、筒型基础长径比d/D和筒外壁-土界面接触强度折减系数 ${\alpha _{{\text{f}}i}}$ 的影响，均质黏土及强度沿深度线性增长的黏土地基中筒型基础地基竖向承载力下限解可拟合为

${N}_{\text{cV}}{}_{（\kappa =0）}\text{=}\left(0.72+0.35{\alpha }_{\text{f}}-0.07\frac{d}{D}\right)\cdot$

$\left[ {5.91 + 8.41\frac{d}{D} - 1.12{{\frac{d}{D}}^2}} \right] \text{，}_{ }$

(27)

${N}_{\text{cV}}（\kappa ）=\left(0.79+0.29{\alpha }_{\text{f}}-0.07\frac{d}{D}\right)\cdot$

$\left[ {3.22 + 2.02\frac{d}{D} + 6.05{{\text{e}}^{ - 0.45\kappa }}} \right] \text{，}$

(28)

$V_{{\rm ult}（κ=0）}=N_{{\rm cV}（κ=0）}·As_{\rm u} \text{，}$

(29)

$V_{{\rm ult}（κ）}=N_{{\rm cV}（κ）}·As_{\rm u0} 。$

(30)

4. 结论

本文针对黏性土地基上筒型基础地基竖向承载力问题，采用轴对称有限元极限分析下限法进行分析，主要得到以下4点结论。

（1）将筒型基础简化为轴对称模型并考虑筒-土间摩擦，求得了筒型基础地基竖向承载力下限解，并对计算结果进行了离心模型试验验证。

（2）在饱和黏土地基中，筒型基础地基竖向承载力系数随土体强度不均匀系数的增大而减小，随筒型基础长径比增长而增长，并给出了拟合公式。

（3）土体抗力和侧摩阻力是筒型基础主要承力部位，随长径比的增加或土体强度不均匀性的降低，土体抗力对承载力贡献率逐渐下降，侧壁摩阻贡献逐渐增长。

（4）对于长径比小于1的筒型基础，筒外壁-土间侧摩阻力随深度呈现先增大后减小趋势；对于长径比大于1的筒型基础筒顶下0.5D范围内外壁侧摩阻发挥较高，而后随深度增加而逐渐降低至近筒端处趋于0。

图 1 土料的级配曲线

Figure 1. Grain-size distribution curves of soil

下载: 全尺寸图片幻灯片

图 2 松散状态下的纯黏土与相变黏土

Figure 2. Clay and phase-change clay in loose state

下载: 全尺寸图片幻灯片

图 3 三轴渗透-压缩试验装置

Figure 3. Conventional triaxial permeability and compression test apparatus

下载: 全尺寸图片幻灯片

图 4 纯黏土（ρ₀=1.80 g/cm³）的应力-应变曲线

Figure 4. Stress-strain curves of clay (ρ₀=1.80 g/cm³)

下载: 全尺寸图片幻灯片

图 5 相变黏土（ρ₀=1.71 g/cm³）的应力-应变曲线

Figure 5. Stress-strain curves of phase-change clay (ρ₀=1.71 g/cm³)

下载: 全尺寸图片幻灯片

图 6 土样在不同围压下加载后的典型形态

Figure 6. Typical shapes of samples after loading under different confining pressures

下载: 全尺寸图片幻灯片

图 7 纯黏土（ρ₀=1.71 g/cm³）的应力-应变曲线

Figure 7. Stress-strain curves of clay (ρ₀=1.71 g/cm³)

下载: 全尺寸图片幻灯片

图 8 不同土样的黏聚力和内摩擦角

Figure 8. Fraction angles and cohesions of different samples

下载: 全尺寸图片幻灯片

图 9 不同围压下土样的切线模量和切线模量的变化曲线

Figure 9. Variation curves of E_t and E_s of samples under different confining pressures

下载: 全尺寸图片幻灯片

图 10 土样滑动摩擦角及相关指标的变化规律

Figure 10. Variation of mobilized friction angles of samples and related indices

下载: 全尺寸图片幻灯片

图 11 不同试样渗透系数随轴向变形的变化规律

Figure 11. Variation of hydraulic conductivity of different samples with axial strain

下载: 全尺寸图片幻灯片

图 12 渗透系数与孔隙比的关系

Figure 12. Relationship between void ratio and hydraulic conductivity

下载: 全尺寸图片幻灯片

图 13 渗透系数与偏应力水平的关系

Figure 13. Relationship between deviatoric stress level and hydraulic conductivity

下载: 全尺寸图片幻灯片

图 14 土样渗透系数实测值与预测曲面的对比

Figure 14. Comparison between test results and predicted surfaces of hydraulic conductivity

下载: 全尺寸图片幻灯片

表 1 PPCM的材料特性

Table 1 Material properties of PPCM

类别	指标	数值
物理特性	液体密度/(kg·L^-1)	0.76
	固体密度/(kg·L^-1)	0.85
	折光率	1.43
	相对介电常数	2.03
	表面张力/(mN·m^-1)	26.5
	黏度系数/(Pa·s)	2.18
	溶解性	不溶于水
	沸点/℃	253.70
	闪点/℃	101.00
热学特性	熔融温度范围/℃	4.12~8.92
	凝固温度范围/℃	0.05~3.45
	熔融吸热/(J·g^-1)	193.30
	凝固放热/(J·g^-1)	193.30
	导热系数/(W·m^-1·℃^-1)	0.20
	比热容/(kJ·kg^-1·℃^-1)	2.00

下载: 导出CSV

表 2 黏土与相变黏土的工程特性

Table 2 Engineering properties of clay and phase-change clay

指标	黏土	相变黏土
PPCM掺量/%	0	8
最优含水率/%	15.15	9.67
最大干密度/(g·cm^-3)	1.84	1.74
压缩模量^①/kPa	13.79	5.41
渗透系数^②/(10⁷ cm·s^-1)	27.4	3.16
无侧限抗压强度^③/kPa	222.08	225.49
注：①单向固结100~200 kPa；②饱和试样；③非饱和试样。

下载: 导出CSV

表 3 试验方案

Table 3 Test schemes

组号	土料	PPCM掺量/%	初始干密度/(g·cm^-3)	围压/kPa	渗透系数测试点-轴向变形/%
1	黏土	0	1.80	25，50，100，200，400	0，3，6，9，12，15，18，21
2	相变黏土	8	1.71	25，50，100，200，400	0，3，6，9，12，15，18，21
3	黏土	0	1.71	100，200，400	0，3，6，9，12，15，18，21

下载: 导出CSV

表 4 渗透系数响应面的求解结果

Table 4 Solved results of response surface of hydraulic.conductivity

参数	土样类型
参数	0%	8%	0%-1.71
A	205.396	87.742	59.756
B	-2.131	-1.876	-0.326
C	-11.205	3.246	-17.182
D	-147.765	-75.164	-27.662
E	5.453	-0.267	5.932
F	8.794	-2.258	-15.431
R²	0.955	0.981	0.990

下载: 导出CSV

参考文献(17)

[1]	穆彦虎, 朱忻怡, 岳攀, 等. 寒区大坝心墙土料冬季冻融与防控监测[J]. 冰川冻土, 2018, 40(4): 756-763. https://www.cnki.com.cn/Article/CJFDTOTAL-BCDT201804012.htm MU Yanhu, ZHU Xinyi, YUE Pan, et al. Monitoring investigation on winter freezing-thawing of dam core wall soils in cold regions[J]. Journal of Glaciology and Geocryology, 2018, 40(4): 756-763. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BCDT201804012.htm
[2]	王效宾, 杨平, 王海波, 等. 冻融作用对黏土力学性能影响的试验研究[J]. 岩土工程学报, 2009, 31(11): 1768-1772. http://www.cgejournal.com/cn/article/id/8440 WANG Xiaobin, YANG Ping, WANG Haibo, et al. Experimental study on effects of freezing and thawing on mechanical properties of clay[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(11): 1768-1772. (in Chinese) http://www.cgejournal.com/cn/article/id/8440
[3]	LIU D H, WANG Y L, LIANG J Y. Potential applications of phase change materials to extend the winter construction time of earth-rock dam in cold regions[J]. Journal of Materials in Civil Engineering, 2021, 33(8): 4021194. doi: 10.1061/(ASCE)MT.1943-5533.0003818
[4]	刘东海, 戴怀建, 郑涵. 心墙相变砾质土工程特性及寒区冬季施工防冻控温研究[J]. 水利学报, 2022, 53(8): 914-925. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202208003.htm LIU Donghai, DAI Huaijian, ZHENG Han. Engineering characteristics and temperature control of phase change material and gravel mixed soil for core wall anti-freezing[J]. Journal of Hydraulic Engineering, 2022, 53(8): 914-925. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202208003.htm
[5]	LIU D H, LIANG J Y, WANG Y L. Numerical simulation on anti-freezing performance of PCM-Clay in core wall during winter construction[J]. Applied Thermal Engineering, 2022, 215: 118951. doi: 10.1016/j.applthermaleng.2022.118951
[6]	刘东海, 郑涵, 杨家琦. 松铺覆盖下心墙相变土防冻控温性能试验研究[J]. 水力发电学报, 2022, 41(7): 38-46. https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB202207005.htm LIU Donghai, ZHENG Han, YANG Jiaqi. Investigation on anti-freezing performance of core wall phase change material-incorporated clay with an upper loose-covering layer[J]. Journal of Hydroelectric Engineering, 2022, 41(7): 38-46. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB202207005.htm
[7]	HUSSAINI S M S, TOUFIGH V. Strength and fracture behavior of rammed-earth materials[J]. Journal of Materials in Civil Engineering, 2019, 31(10): 4019228. doi: 10.1061/(ASCE)MT.1943-5533.0002876
[8]	SERRANO S, BARRENECHE C, RINCÓN L, et al. Optimization of three new compositions of stabilized rammed earth incorporating PCM: thermal properties characterization and LCA[J]. Construction and Building Materials, 2013, 47: 872-878. doi: 10.1016/j.conbuildmat.2013.05.018
[9]	RAO Y Z, LIU J K, CHANG D, et al. Effects of microencapsulated phase change material characteristics on the thermal performance and mechanical behaviour of silty clay[J]. Transportation Geotechnics, 2021, 29: 100584. doi: 10.1016/j.trgeo.2021.100584
[10]	黄英豪, 陈永, 朱洵, 等. 相变材料改良膨胀土冻融性能试验研究及微观机理分析[J]. 岩土工程学报, 2021, 43(11): 1994-2002. doi: 10.11779/CJGE202111005 HUANG Yinghao, CHEN Yong, ZHU Xun, et al. Experimental study and micro-mechanism analysis of freeze-thaw performance of expansive soils improved by phase-change materials[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(11): 1994-2002. (in Chinese) doi: 10.11779/CJGE202111005
[11]	YONG R N, RAO S M. Mechanistic evaluation of mitigation of petroleum hydrocarbon contamination by soil medium[J]. Canadian Geotechnical Journal, 1991, 28(1): 84-91.
[12]	FERNANDEZ F, QUIGLEY R M. Hydraulic conductivity of natural clays permeated with simple liquid hydrocarbons[J]. Canadian Geotechnical Journal, 1985, 22(2): 205-214.
[13]	土工试验方法标准GB/T 50123—2019[S]. 北京: 中国计划出版社, 2019. Standard for Geotechnical Testing Method GB/T 50123—2019[S]. Beijing: China Planning Press, 2019. (in Chinese)
[14]	碾压式土石坝设计规范: SL 274—2020[S]. 北京: 中国水利水电出版社, 2020. Design Specification for Rolled Earth-Rockfill Dams DL/T 5395—2020[S]. Beijing: China Water Power Press, 2020. (in Chinese)
[15]	朱建华. 土坝心墙原状土的三轴渗透试验[J]. 岩土工程学报, 1989, 11(4): 57-63. http://www.cgejournal.com/cn/article/id/9292 ZHU Jianhua. Triaxial permeability test on undisturbed soils of earth dam core[J]. Chinese Journal of Geotechnical Engineering, 1989, 11(4): 57-63. (in Chinese) http://www.cgejournal.com/cn/article/id/9292
[16]	雷红军. 高土石坝黏性土大剪切变形条件下渗透特性研究[D]. 北京: 清华大学, 2010. LEI Hongjun. A Study on Seepage Characteristics of Clayey Soil of High Earth-Rockfill Dam with Large Shear Deformation[D]. Beijing: Tsinghua University, 2010. (in Chinese)
[17]	王刚, 韦林邑, 魏星, 等. 压实黏土三轴压缩变形过程中的渗透性变化规律[J]. 岩土力学, 2020, 41(1): 32-38. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202001005.htm WANG Gang, WEI Linyi, WEI Xing, et al. Permeability evolution of compacted clay during triaxial compression[J]. Rock and Soil Mechanics, 2020, 41(1): 32-38. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202001005.htm