非饱和黏性土的广义有效应力三剪弹塑性本构模型

胡小荣; 路祥; 蔡晓锋

doi:10.6052/j.issn.1000-4750.2022.05.0410

非饱和黏性土的广义有效应力三剪弹塑性本构模型

南昌大学建筑工程学院，南昌 330031

基金项目: 国家自然科学基金项目(10862004，11362016)

详细信息

作者简介:
胡小荣(1964−)，男，江西人，教授，博士，主要从事岩土本构理论及其工程应用的科研教学工作研究(E-mail: huxiaorong@ncu.edu.cn)

蔡晓锋(1997−)，男，浙江人，硕士，主要从事岩土本构理论及其工程应用研究(E-mail: caixiaofeng4399@163.com)

通讯作者:
路　祥(1994−)，男，安徽人，硕士，主要从事岩土本构理论及其工程应用研究(E-mail: 2253086841@qq.com)

中图分类号: TU432
计量
- 文章访问数: 195
- HTML全文浏览量: 36
- PDF下载量: 65
出版历程
- 收稿日期: 2022-05-04
- 修回日期: 2022-11-25
- 网络出版日期: 2023-03-05
- 刊出日期: 2024-05-24

THE TRIPLE-SHEAR UNIFIED CONSTITUTIVE MODEL FOR UNSATURATED CLAYS BASED ON THE GENERALIZED EFFECTIVE STRESS CONCEPT

School of Civil Engineering and Architecture, Nanchang University, Nanchang 330031, China

摘要

摘要:
将非饱和土广义有效应力原理与三剪强度准则相结合，提出了非饱和土广义有效应力三剪强度准则。将所提准则作为破坏准则，分别采用等量代换法和坐标平移法推导出新的破坏应力比，并将其与非饱和土修正剑桥模型相结合得到了新的屈服函数。相比于原来修正剑桥模型中的破坏应力比为定值，新的屈服函数可以更好地反应土体全应力状态、中间主应力效应和拉压不等效应。在弹塑性理论的框架下，建立了非饱和土的广义有效应力三剪弹塑性本构模型。以江西正常固结非饱和重塑红黏土作为试验研究对象，进行室内土工试验、土水特征曲线试验、压缩回弹试验、非饱和土常规三轴固结排水试验。将该本构模型计算结果与非饱和土三轴固结排水试验结果进行对比验证。结果表明：数值模拟结果与试验结果吻合较好，验证了该本构模型的正确性。在轴向应变较小时，等量代换法和坐标平移法模拟结果比较接近，随着轴向应变增大直至偏应力达到平稳状态的过程中，等量代换法计算结果要大于坐标平移法计算结果，且更接近于试验值。真三轴计算预测结果表明：在固结排水条件下，初始压实度、净围压、基质吸力相同，中间主应力影响系数越大，则剪应力和体应变越大。b值相同的情况下，轴向应变较小时，等量代换法计算结果和坐标平移法计算结果比较接近，随着轴向应变的增大，二者之间的差值也越来越大。
- 非饱和黏性土 /
- 本构模型 /
- 广义有效应力 /
- 三剪强度准则 /
- 江西红黏土
Abstract:
Combining the triple-shear strength criterion with the generalized effective stress principle of unsaturated soils, the triple-shear strength criterion based on generalized effective stress for unsaturated soils was proposed. The new failure stress ratio was derived by using the equivalent substitution method and coordinate translation method, respectively. The new failure stress ratio was substituted into the modified Cambridge model for unsaturated soils to obtain the new yield function. Compared with the original modified Cambridge model whose failure stress ratio is constant, the new yield function can better reflect the total stress state, intermediate principal stress effect and tension-compression unequal effect of soils. A unified triple-shear elastoplastic constitutive model for unsaturated soils was established based on the generalized effective stress variable method under the framework of elasto-plastic theory. Taking Jiangxi remolded unsaturated clay in normal consolidation as the experimental research object, the paper carried out indoor geotechnical test, soil-water characteristic curve test, compression rebound test and conventional triaxial consolidation drainage test of unsaturated soil. The results of the constitutive model were compared with the results of unsaturated soil triaxial consolidation and drainage test. The results show that the numerical simulation results are in good agreement with the experimental results, which verifies the correctness of the constitutive model. When the axial strain is small, the simulation results of the equivalent substitution method and the coordinate translation method are relatively close. With the increase of the axial strain until the deviant stress reaches the steady state, the calculation results of the equal substitution method are greater than those of the coordinate translation method, and are closer to the experimental values. The prediction results of true triaxial calculation show that the initial compacting degree, net confining pressure and matric suction are the same under consolidation drainage condition, and the shear stress and bulk strain increase with the influence coefficient of intermediate principal stress. When the value of b is the same and the axial strain is small, the calculation results of the equivalent substitution method and the coordinate translation method are relatively close. With the increase of the axial strain, the difference between the two methods becomes larger.
- generalized effective stress /
- unsaturated cohesive soil /
- the triple-shear strength criterion /
- constitutive model /
- Jiangxi clay

HTML全文

ALONSO等^[1]在修正剑桥模型的基础上提出了非饱和土的巴塞罗纳模型(也称为BBM模型)。BBM模型采用净应力和基质吸力作为模型中的应力状态变量，其特点在于提出的LC屈服线可以很好地描述非饱和土的屈服应力随基质吸力变化的性质。SHENG等^[2]基于双应力变量法提出了一个新的本构模型—SFG模型。该模型认为净应力和基质吸力对体应变的影响需要分开考虑，同时该模型采用积分的方法来计算LC屈服面方程。SFG模型可以很好地表现基质吸力对土体压缩性的影响，但该模型无法从非饱和状态退化到饱和状态，且无法考虑饱和度对非饱和土性质的影响。刘艳^[3]基于广义有效应力原理，采用连续介质力学中热力学的方法建立了固液气三相耦合的非饱和土本构模型，其特点在于可综合考虑固液气三相对非饱和土力学性质的影响，还考虑到了饱和度的影响。目前，大多数的非饱和土本构模型都是在临界状态土力学理论框架下建立起来的^[4-6]，而破坏应力比是一个非常重要的参数，由土的强度准则得出，强度准则大多采用基于某个有效应力原理表示的摩尔-库仑强度准则^[7] ，但该准则没有考虑中间主应力的影响。为了克服上述弊端，一些研究者相继提出了非饱和土的SMP准则^[8]和双剪统一强度准则^[9]等。其中，SMP准则也因其模型特征单一不能反映土体的多样性。而双剪统一强度准则虽然能较全面地表征非饱和土的强度特性，但存在双重破坏角的现象^[10]。而三剪统一强度准则^[10-11]则能够较全面地反映土的拉压不等效应、中间主应力效应和应力区间效应等特性，同时还不会出现双重破坏角的问题，而且该强度准则是一个强度准则系列，可以对其他准则作非线性逼近，由此以反映土体的多样性。在非饱和土的有效应力原理方面，采用较多的是Bishop应力变量法^[12]和Fredlund应力变量法^[13]。但是，Bishop应力变量法^[12]存在有效应力系数难以确定并且和饱和度对应关系不唯一的问题^[14-15]，Fredlund应力变量法^[13]最明显的缺点就是无法考虑饱和度的影响^[16]且无法描述土体从非饱和状态到饱和状态间的平稳过渡^[17-18]，原因是双应力变量有效应力原理在非饱和土基质吸力为零时不能与饱和土的有效应力原理相衔接。而赵成刚等^[19-21]提出的非饱和土广义有效应力原理可以克服上述缺点。

综上所述，本文的研究内容主要有：

1)将广义有效应力变量法与三剪强度准则^[22-23]相结合，得到非饱和土广义有效应力三剪强度准则，并进行验证。为了更好地在修正剑桥模型中反映非饱和黏性土黏聚力的影响，采用文献[24-25]提出的等量代换法和坐标平移法推导出相应的破坏应力比。

2)将所得破坏应力比与非饱和土修正剑桥模型相结合得到屈服函数。

3)在弹塑性理论的框架下，建立弹塑性本构模型。

4)用江西正常固结非饱和重塑红黏土对模型进行验证。验证内容为不同压实度、不同净围压、不同基质吸力条件下土在常规三轴固结排水试验中的偏应力和轴向应变、体应变和轴向应变的关系。

5)用所建立的本构模型对该红黏土进行真三轴固结排水数值模拟，并探究中间主应力的影响。

1 非饱和土三剪强度准则与破坏应力比

将非饱和土广义有效应力^[19]代入到三剪强度准则^[10-11]中，得到基于广义有效应力的非饱和土三剪强度准则为：

$\begin{split}& [ {{{\left( {{{\widetilde \sigma }_1} - {{\widetilde \sigma }_3}} \right)}^2} + b{{\left( {{{\widetilde \sigma }_1} - {{\widetilde \sigma }_2}} \right)}^2} + b{{\left( {{{\widetilde \sigma }_2} - {{\widetilde \sigma }_3}} \right)}^2}}] - \\&\qquad\left( {1 + b} \right)\left[ {\left( {{{\widetilde \sigma }_1} + {{\widetilde \sigma }_3}} \right)\left( {{{\widetilde \sigma }_1} - {{\widetilde \sigma }_3}} \right)} \right]\sin\varphi '= \\&\qquad 2c'\left( {1 + b} \right)\left( {{{\widetilde \sigma }_1} - {{\widetilde \sigma }_3}} \right)\cos\varphi ' \end{split}$

(1)

式中： $c'$ 为饱和有效黏聚力； $\varphi '$ 为饱和有效内摩擦角； $b$ 为中间主应力影响系数； ${\widetilde \sigma '_1}$ 、 ${\widetilde \sigma '_2}$ 、 ${\widetilde \sigma '_3}$ 分别为最大、中间和最小广义有效主应力。

非饱和土广义有效应力具体表达式为^[3]：

${\widetilde \sigma _{ij}} = {\sigma _{ij}} - {p_{\text{a}}}{\delta _{ij}} + {s_{\text{r}}}s{\delta _{ij}}$

(2)

式中： ${\widetilde \sigma _{ij}}$ 为广义有效应力； ${s_{\text{r}}}$ 为饱和度； $s$ 为基质吸力； ${p_{\text{a}}}$ 为孔隙气压力； ${\delta _{ij}}$ 为Kronecker符号。

其广义有效主应力为：

$\left\{ \begin{gathered} {{\widetilde \sigma }_1} = {\sigma _1} - {p_{\text{a}}} + {s_{\text{r}}}s = {\sigma '_1} + {s_{\text{r}}}s \\ {{\widetilde \sigma }_2} = {\sigma _2} - {p_{\text{a}}} + {s_{\text{r}}}s = {\sigma'_2} + {s_{\text{r}}}s \\ {{\widetilde \sigma }_3} = {\sigma _3} - {p_{\text{a}}} + {s_{\text{r}}}s = {\sigma '_3} + {s_{\text{r}}}s \\ \end{gathered} \right.$

(3)

式中： ${\sigma _1}$ 、 ${\sigma _2}$ 和 ${\sigma _3}$ 分别为最大、中间和最小总主应力； ${\sigma '_1}$ 、 ${\sigma '_2}$ 和 ${\sigma '_3}$ 分别为最大、中间和最小净主应力。

将式(3)代入式(1)得到用净应力表示的基于广义有效应力的非饱和土三剪强度准则为：

$\begin{split}& [ {{{\left( {{\sigma '_1} - {\sigma '_3}} \right)}^2} + b{{\left( {{\sigma '_1} - {\sigma '_2}} \right)}^2} + b{{\left( {{\sigma '_2} - {\sigma '_3}} \right)}^2}} ] - \\&\qquad \left( {1 + b} \right)\left[ {\left( {{\sigma '_1} - {\sigma '_3}} \right)\left( {{\sigma '_1} + {\sigma '_3} + 2{s_{\text{r}}}s} \right)} \right]\sin\varphi ' = \\&\qquad 2c'\left( {1 + b} \right)\left( {{\sigma '_1} - {\sigma '_3}} \right)\cos\varphi ' \end{split}$

(4)

净主应力及其相关的应力不变量关系式为^[26-27]：

$\left\{ \begin{gathered} {\sigma '_1} = p' + \sqrt {\frac{2}{3}} \rho \cos\theta = p' + \frac{2}{3}q\cos\theta \\ {\sigma '_2} = p' + \sqrt {\frac{2}{3}} \rho \cos\left( {\frac{2}{3}\pi - \theta } \right) = p' + \frac{2}{3}q\cos\left( {\frac{2}{3}\pi - \theta } \right) \\ {\sigma '_3} = p' + \sqrt {\frac{2}{3}} \rho \cos\left( {\frac{2}{3}\pi + \theta } \right) = p' + \frac{2}{3}q\cos\left( {\frac{2}{3}\pi + \theta } \right) \\ \end{gathered} \right.$

(5)

式中： $p'$ 为平均净主应力； $\rho$ 为 $\pi$ 平面上极限线上的点到静水压力轴的垂直距离； $q$ 为广义剪应力； $\theta$ 为应力角( ${0 \leqslant \theta \leqslant {{60} }}$ °)。

由式(4)和式(5)得到强度准则 $p'$ - $q$ 子午线表达式为：

$\begin{split} q =& A\left[ {\sin\varphi '\left( {p - {p_{\text{a}}} + {s_{\text{r}}}s} \right) + c'\cos\varphi '} \right] = \\& Ap'\sin\varphi ' + A\left( {{s_{\text{r}}}s\sin\varphi ' + c'\cos\varphi '} \right) \end{split}$

(6)

其中：

$A = \frac{{6\left( {1 + b} \right)\cos\left( {\theta - \dfrac{\pi }{6}} \right)}}{\begin{split} & \Bigg\langle 2\sqrt 3 \left[ {{\rm{co}}{{\rm{s}}^2}\left( {\theta - \frac{\pi }{6}} \right) + b{\rm{co}}{{\rm{s}}^2}\left( {\theta + \dfrac{\pi }{6}} \right) + b{\rm{si}}{{\rm{n}}^2}\theta } \right] - \\& \left( {1 + b} \right)\sin\varphi '\cos\left( {2\theta + \frac{\pi }{6}} \right) \Bigg\rangle \end{split} }$

由式(6)得到基于等量代换法的破坏应力比为：

$M\left( {p',\theta ,s,{s_{\rm{r}}}} \right) = \frac{q}{{p'}} = A\sin \varphi ' + \frac{{A\left( {{s_{\rm{r}}}s\sin \varphi ' + c'\cos \varphi '} \right)}}{{p'}}$

(7)

修正剑桥模型中临界状态极限线是过坐标原点的线，黏聚力不等于零的黏性土并不符合其条件。采用坐标平移法使子午线经过原点，即将坐标轴向左平移使其经过原点。基于坐标平移法下新的非饱和土三剪破坏应力比为：

$M\left( \theta \right){\text{ = }}\frac{{\overline q}}{{\overline p'}} = A\sin \varphi '$

(8)

坐标平移前后其它应力状态量间的关系为：

${\overline I'_1} = {I'_1} + 3\left( {c'\cot \varphi ' + {s_{\rm{r}}}s} \right) ,\; {\overline J_2}{\text{ = }}{J_2} ,\; {\overline J_3}{\text{ = }}{J_3} ,\; \overline \theta = \theta$

(9)

式中： ${\overline I'_1}$ 、 ${\overline J_2}$ 、 ${\overline J_3}$ 和 $\overline \theta$ 分别为坐标平移后的非饱和土净应力的第一应力不变量、第二应力偏张量不变量、第三应力偏张量不变量和Lode角。

坐标平移前后的土体压缩指数关系为：

$\overline \lambda {\text{ = }}{\varLambda _1}\lambda$

(10)

式中， $\lambda$ 和 $\overline \lambda$ 分别为坐标平移前和坐标平移后的土体压缩指数。

坐标平移前后的土体回弹指数关系为：

$\overline \kappa {\text{ = }}{\varLambda _1}\kappa$

(11)

式中， $\kappa$ 和 $\overline \kappa$ 分别为坐标平移前和坐标平移后的土体回弹指数。其中：

${\varLambda _1}{\text{ = }}\frac{{\overline \lambda }}{\lambda }{\text{ = }}\dfrac{{\ln \left( {c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}}{{\ln \left( {\dfrac{{2c'\cot \varphi ' + 2{s_{\rm{r}}}s}}{{1 + c'\cot \varphi ' + {s_{\rm{r}}}s}}} \right)}}$

文献[]中的非饱和黄土在 $s = 100\;{\text{kPa}}$ 时所对应的饱和度 ${s_{\text{r}}} = 34.5\text{%}$ ，有效黏聚力 $c' = 5.3\;{\text{kPa}}$ ，有效内摩擦角 $\varphi '{\text{ = }}27.92$ °。图1为本文新建立的非饱和土广义有效应力三剪强度准则计算结果与文献[28]中的真三轴试验数据对比验证图。

由可知，取 $b{\text{ = }}0.25$ 的非饱和土广义有效应力三剪强度准则计算预测结果与试验结果较为吻合，验证了本文所推导的新准则的正确性。新准则采用了广义有效应力原理，具备广义有效应力原理的优点，考虑了饱和度的影响，可以更全面地描述非饱和土的力学特性。

图 1 广义有效应力变量法计算结果和试验数据对比图

Figure 1. Comparison diagram of calculation results of generalized effective stress variable method and test data

下载: 全尺寸图片幻灯片

由可知，将广义有效应力变量法预测的 $p'$ - $q$ 子午线与试验点相比对，广义有效应力变量法预测的 $p'$ - $q$ 子午线基本与试验点吻合，新准则取得很好的预测效果。

图 2 子午线计算结果与试验数据对比图

Figure 2. Comparison diagram of meridian calculation results and experimental data

下载: 全尺寸图片幻灯片

文献[-]中的试验所用的非饱和黏土砂有效黏聚力 $c' = 0\;{\rm{kPa}}$ ，有效内摩擦角 $\varphi ' = 30$ °，基质吸力 $s = 200\;{\rm{kPa}}$ ，饱和度 ${s_{\text{r}}} = 37.5\text{%}$ ，平均净主应力 $p' = 100\;{\rm{kPa}}$ 。图2为本文新建立的非饱和土广义有效应力三剪强度准则计算结果与文献[29-30]中的真三轴试验数据对比验证图。

由图3可知，将本文所推导的破坏准则预测结果与单应力变量法和双应力变量法的三剪强度准则预测结果相比较，可以得到广义有效应力变量法最为吻合试验结果，同时和双应力变量法比较接近。单应力变量法预测误差更大。新准则虽未对双应力变量法三剪强度准则产生明显的优势，但其采用了广义有效应力原理，能更全面地描述非饱和土的力学特性。

图 3 广义有效应力变量法计算结果和试验数据对比图

Figure 3. Comparison diagram of calculation results of generalized effective stress variable method and test data

下载: 全尺寸图片幻灯片

由可知，将广义有效应力变量法预测的 $p'$ - $q$ 子午线与试验点相比对，文献中的试验数据点比较少，试验点基本与 $\theta = {0}$ °的广义有效应力变量法下的 $p'$ - $q$ 子午线吻合，新准则取得很好的预测效果。

图 4 子午线计算结果与试验数据对比图

Figure 4. Comparison diagram of meridian calculation results and experimental data

下载: 全尺寸图片幻灯片

2 非饱和黏性土三剪统一弹塑性本构模型

在弹塑性力学理论框架下，将第1节所推导的破坏应力比引入到非饱和土修正剑桥模型中，建立基于广义有效应力变量法的非饱和土三剪统一弹塑性本构模型。采用相关联的流动法则，由一致性条件和硬化规律推导出本构模型中的塑性模量和弹塑性刚度矩阵各元素的具体表达式。

2.1 屈服函数

饱和土修正剑桥模型表达式为：

$f = {q^2} + {M^2}{p'^2} - {M^2}p'{p'_x}$

(12)

式中： $f$ 为屈服函数； ${p'_x}$ 为饱和土的有效屈服应力； $p'$ 为饱和土的平均有效主应力； $q$ 为偏应力。

用非饱和土的净屈服应力对饱和土修正剑桥模型中的有效屈服应力进行替换得到非饱和土修正剑桥模型为：

$f = {q^2} + {M^2}{p'^2} - {M^2}p'{p'_y}\left( s \right)$

(13)

式中： $p'$ 为非饱和土净平均主应力； $M$ 为破坏应力比； ${p'_y}\left( s \right)$ 为非饱和土的屈服应力。

文献[31]曾提出过非饱和土的LC屈服曲线，但其中土的回弹指数为一定值，没有考虑基质吸力的影响。文献[32]通过试验研究发现回弹指数并非是一定值，而是随基质吸力变化的量，并将LC屈服曲线改写为：

${p'_y}\left( s \right) = {p'_{{n}}}{\left( {\frac{{{p'_y}\left( 0 \right)}}{{{p'_{{n}}}}}} \right)^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}}}}$

(14)

式中： ${p'_y}\left( s \right)$ 为非饱和土的净屈服应力； ${p'_y}\left( 0 \right)$ 为饱和土的有效(净)屈服应力； ${p'_{\rm{n}}}$ 为吸湿不发生湿化变形的应力； $\lambda \left( 0 \right)$ 为饱和土的压缩指数； $\kappa \left( 0 \right)$ 为饱和土的回弹指数； $\lambda \left( s \right)$ 为基质吸力为s的非饱和土的压缩指数； $\kappa \left( s \right)$ 为基质吸力为s的非饱和土的回弹指数。由文献[32]可知，非饱和土压缩指数和回弹指数与基质吸力间的关系为：

$\left\{ \begin{gathered} \lambda \left( s \right) = \lambda \left( 0 \right) - \frac{{{\lambda _{\rm{s}}}s}}{{{p_{\text{1}}} + s}} \\ \kappa \left( s \right) = \kappa \left( 0 \right) + {\kappa _{\rm{s}}}s \\ \end{gathered} \right.$

(15)

式中： ${\lambda _{\rm{s}}}$ 为表达非饱和土压缩指数随基质吸力大小变化的参数； ${\kappa _{\rm{s}}}$ 为表达非饱和土回弹指数随基质吸力大小变化的参数； ${p_{\text{1}}}$ 为大气压。

将广义有效应力变量法得到的非饱和土三剪统一破坏比代入到非饱和土的修正剑桥模型中，得到 $p'$ - $q$ 子午面上的非饱和土屈服线函数为：

$f = {q^2} + M{\left( {p',\theta ,s,{s_{\text{r}}}} \right)^2}p'\left( {p' - {p'_{{n}}}{{\left( {\frac{{{p'_y}\left( 0 \right)}}{{{p'_{{n}}}}}} \right)}^{\frac{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}{{\lambda \left( s \right) - \kappa \left( s \right)}}}}} \right)$

(16)

饱和土的应变硬化函数为：

${\rm{ln}}{p'_y}\left( 0 \right) - {\rm{ln}}{p'_0} = \frac{{1 + {e_0}}}{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}\varepsilon _{\rm{v}}^{\rm{p}}$

(17)

式中： ${p'_0}$ 为初始平均净主应力； $\varepsilon _{\rm{v}}^{\rm{p}}$ 为塑性体应变。

将式(17)代入式(16)，将式(16)写成椭圆方程形式，得到基于广义有效应力的非饱和土屈服函数表达式如下：

1) 采用等量代换法时，

$\begin{split} f = &{\left\{ {p' - \frac{1}{2}{p'_{{n}}}{{\left[ {\frac{{{p'_0}}}{{{p'_{{n}}}}}\exp\left( {\frac{{1 + {e_0}}}{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)} \right]}^{\frac{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}{{\lambda \left( s \right) - \kappa \left( s \right)}}}}} \right\}^2} + \\&\frac{{{q^2}}}{{M{{\left( {p',\theta ,s,{s_{\rm{r}}}} \right)}^2}}} - \\& \frac{1}{4}{\left\{ {{p'_{{n}}}\left[ {\frac{{{p'_0}}}{{{p'_{{n}}}}}\exp{{\left( {\frac{{1 + {e_0}}}{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)}^{\frac{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}{{\lambda \left( s \right) - \kappa \left( s \right)}}}}} \right]} \right\}^2} = 0 \end{split}$

(18)

2) 采用坐标平移法时，

$\begin{split} & f = \left\{ \left( {p' + c'\cot \varphi ' + {s_{\rm{r}}}s} \right) - \frac{1}{2}\left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)\times\right.\\&\left.{{\left[ {\frac{{\left( {{p'_0} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}}{{\left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}}\exp\left( {\frac{{1 + {e_0}}}{{{\varLambda _1}\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)} \right]}^{\frac{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}{{\lambda \left( s \right) - \kappa \left( s \right)}}}} \right\}^2 +\\& \frac{{{q^2}}}{{M{{\left( \theta \right)}^2}}} - \frac{1}{4}\Bigg\{ \left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)\times\Bigg[ \frac{{\left( {{p'_0} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}}{{\left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}} \\&\left.\exp{{\left( {\frac{{1 + {e_0}}}{{{\varLambda _1}\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)}^{\frac{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}{{\lambda \left( s \right) - \kappa \left( s \right)}}}} \Bigg] \right\}^2 = 0 \end{split}$

(19)

2.2 流动法则

采用相关联的流动的流动法则，屈服函数与塑性势函数相同，且塑性流动方向与加载方向一致。

${L_{ij}} = \frac{{\partial f}}{{\partial {\sigma '_{ij}}}}$

(20)

式中， ${L_{ij}}$ 为加载方向。

2.3 塑性模量

塑性应力和应变之间的增量关系为：

$\left\{ \begin{gathered} {\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}} = \left\langle L \right\rangle \frac{{\partial f}}{{\partial p'}} \\ {\rm{d}}\varepsilon _{\rm{q}}^{\rm{p}} = \left\langle L \right\rangle \frac{{\partial f}}{{\partial q}} \\ \end{gathered} \right.$

(21)

式中： ${\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}$ 为塑性体应变增量； ${\rm{d}}\varepsilon _{\rm{q}}^{\rm{p}}$ 为塑性剪应变增量。

加载指数和塑性模量之间的关系为：

$L = \frac{1}{{{K_{\rm{p}}}}}\frac{{\partial f}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}$

(22)

式中， ${K_{\rm{p}}}$ 为塑性模量。

由一致性条件可得：

$\frac{{\partial f}}{{\partial p'}}{\rm{d}}p' + \frac{{\partial f}}{{\partial q}}{\rm{d}}q + \frac{{\partial f}}{{\partial \varepsilon _{\rm{v}}^{\rm{p}}}}{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}} = 0$

(23)

式中： ${\rm{d}}{\rm{}}\varepsilon _{\rm{v}}^{\rm{p}}$ 为塑性体应变增量； ${\rm{d}}p' = \dfrac{{\partial p'}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}$ ； ${\rm{d}}q = \dfrac{{\partial q}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}$ 。

由式(21)和式(22)得：

${\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}} = \frac{1}{{{K_{\rm{p}}}}}\frac{{\partial f}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}\frac{{\partial f}}{{\partial p'}}$

(24)

将式(24)代入式(23)得：

$\frac{{\partial f}}{{\partial p'}}\frac{{\partial p'}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}} + \frac{{\partial f}}{{\partial q}}\frac{{\partial q}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}} + \frac{{\partial f}}{{\partial \varepsilon _{\rm{v}}^{\rm{p}}}}\frac{1}{{{K_{\rm{p}}}}}\frac{{\partial f}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}\frac{{\partial f}}{{\partial p'}} = 0$

(25)

根据复合函数求导法则得：

$\frac{{\partial f}}{{\partial {{\sigma _{ij}'}}}}{\rm{d}}{\sigma{ '_{ij}}}{\text{ = }}\frac{{\partial f}}{{\partial p'}}\frac{{\partial p'}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}} + \frac{{\partial f}}{{\partial q}}\frac{{\partial q}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}$

(26)

将式(26)代入式(25)得：

$\frac{{\partial f}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}} + \frac{{\partial f}}{{\partial \varepsilon _{\rm{v}}^{\rm{p}}}}\frac{1}{{{K_{\rm{p}}}}}\frac{{\partial f}}{{\partial {{\sigma '}_{ij}}}}{\rm{d}}{\sigma '_{ij}}\frac{{\partial f}}{{\partial p'}} = 0$

(27)

由式(27)求得塑性模量的表达式为：

${K_{\rm{p}}} = - \frac{{\partial f}}{{\partial \varepsilon _{\rm{v}}^{\rm{p}}}}\frac{{\partial f}}{{\partial p'}}$

(28)

其中：

1) 使用等量代换法时，

$\begin{split} \frac{{\partial f}}{{\partial p'}} =& 2\left\{ {p' - \frac{1}{2}{p'_{{n}}}{{\left[ {\frac{{{{p'}_0}}}{{{p'_{{n}}}}}\exp\left( {\frac{{1 + {e_0}}}{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)} \right]}^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}}}}} \right\} +\\& \frac{{2{q^2}}}{{M{{\left( {p',\theta ,s,{s_{\rm{r}}}} \right)}^3}}}\frac{{A\left( {c'\cos \varphi ' + {s_{\rm{r}}}s\sin \varphi '} \right)}}{{{{p'}^2}}}, \end{split}$

$\begin{split} \frac{{\partial f}}{{\partial \varepsilon _{\rm{v}}^{\rm{p}}}} =& - p'{p'_{{n}}}{\left( {\frac{{{{p'}_0}}}{{{p'_{{n}}}}}} \right)^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}}}}\frac{{1 + {e_0}}}{{\lambda \left( s \right) - \kappa \left( s \right)}}\\&{\left[ {\exp\left( {\frac{{1 + {e_0}}}{{\lambda \left( 0 \right) - \kappa \left( 0 \right)}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)} \right]^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}}}} 。 \end{split}$

2) 使用坐标平移法时，

$\begin{split} \frac{{\partial f}}{{\partial p'}} =& 2\left\{ \left( {p' + c'\cot \varphi ' + {s_{\rm{r}}}s} \right) - \frac{1}{2}\left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)\times\right.\\&\left.{{\left[ {\frac{{\left( {{p'_0} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}}{{\left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)}} \exp \left( {\frac{{1 + {e_0}}}{{{\varLambda _1}\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)} \right]}^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)}\right]}}}} \right\}, \end{split}$

$\begin{split} \frac{{\partial f}}{{\partial \varepsilon _{\rm{v}}^{\rm{p}}}} = & - \left( {p' + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)\left( {{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)\times\\&{\left( {\frac{{{p'_0} + c'\cot \varphi ' + {s_{\rm{r}}}s}}{{{p'_{{n}}} + c'\cot \varphi ' + {s_{\rm{r}}}s}}} \right)^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}}}}\frac{{1 + {e_0}}}{{{\varLambda _1}\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}} \times\\& {\left[ {\exp\left( {\frac{{1 + {e_0}}}{{{\varLambda _1}\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}\varepsilon _{\rm{v}}^{\rm{p}}} \right)} \right]^{\frac{{\left[ {\lambda \left( 0 \right) - \kappa \left( 0 \right)} \right]}}{{\left[ {\lambda \left( s \right) - \kappa \left( s \right)} \right]}}}} 。 \end{split}$

2.4 弹塑性本构关系

由弹塑性应变增量可知：

${\rm{d}}{\varepsilon _{ij}} = {\rm{d}}\varepsilon _{ij}^{\rm{e}} + {\rm{d}}\varepsilon _{ij}^{\rm{p}}$

(29)

式中： ${\rm{d}}{\varepsilon _{ij}}$ 为应变增量； ${\rm{d}}\varepsilon _{ij}^{\rm{e}}$ 为弹性应变增量； ${\rm{d}}\varepsilon _{ij}^{\rm{p}}$ 为塑性应变增量。

又因为：

${\rm{d}}{\sigma '_{ij}} = {{\boldsymbol{D}}^{\rm{e}}}{\rm{d}}\varepsilon _{ij}^{\rm{e}} = {{\boldsymbol{D}}^{\rm{e}}}( {{\rm{d}}{\varepsilon _{ij}} - {\rm{d}}\varepsilon _{ij}^{\rm{p}}} )$

(30)

式中： ${\rm{d}}{\sigma '_{ij}}$ 为净应力增量； ${{\boldsymbol{D}}^{\rm{e}}}$ 为弹性刚度矩阵。

${\rm{d}}\varepsilon _{ij}^{\rm{p}} = \left\langle L \right\rangle \frac{{\partial f}}{{\partial {\sigma '_{ij}}}}$

(31)

将式(21)进行变换得式(32)：

${K_{\rm{p}}}\left\langle L \right\rangle = \frac{{\partial f}}{{\partial {\sigma '_{ij}}}}{\rm{d}}{\sigma '_{ij}}$

(32)

将式(30)和式(31)代入式(32)得：

${K_{\rm{p}}}\left\langle L \right\rangle = \frac{{\partial f}}{{\partial {\sigma '_{ij}}}}\left[ {{{\boldsymbol{D}}^{\rm{e}}}\left( {{\rm{d}}{\varepsilon _{ij}} - \left\langle L \right\rangle \frac{{\partial f}}{{\partial {\sigma '_{ij}}}}} \right)} \right]$

(33)

由此得到加载指数的表达式为：

$\left\langle L \right\rangle = \dfrac{{\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}{{\boldsymbol{D}}^{\rm{e}}}{\rm{d}}{\varepsilon _{ij}}}}{{{K_{\rm{p}}} + \dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}{{\boldsymbol{D}}^{\rm{e}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}}}$

(34)

将式(31)和式(34)代入式(30)得：

$\begin{split} {\rm{d}}{\sigma '_{ij}} =& {{\boldsymbol{D}}^{\rm{e}}}\left( {{\rm{d}}{\varepsilon _{ij}} - \dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}\left\langle L \right\rangle } \right) = \\& \left( {{{\boldsymbol{D}}^{\rm{e}}} - \frac{{{{\boldsymbol{D}}^{\rm{e}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}{{\boldsymbol{D}}^{\rm{e}}}}}{{{K_{\rm{p}}} + \dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}{{\boldsymbol{D}}^{\rm{e}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}}}} \right){\rm{d}}{\varepsilon _{ij}} = {{\boldsymbol{D}}^{{\rm{ep}}}}{\rm{d}}{\varepsilon _{ij}} \end{split}$

(35)

式中， ${{\boldsymbol{D}}^{{\rm{ep}}}}$ 为弹塑性刚度矩阵。

从而得到弹塑性刚度矩阵为：

$\begin{split} {{\boldsymbol{D}}^{\rm{ep}}} =& {{\boldsymbol{D}}^{\rm{e}}} - \dfrac{{{{\boldsymbol{D}}^{\rm{e}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}{{\boldsymbol{D}}^{\rm{e}}}}}{{{K_{\rm{p}}} + \dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}{{\boldsymbol{D}}^{\rm{e}}}\dfrac{{\partial f}}{{\partial {\sigma '_{ij}}}}}} = \\& \left[ {\begin{array}{*{20}{c}} {D_{11}^{\rm{ep}}}&{D_{12}^{\rm{ep}}}&{D_{13}^{\rm{ep}}}&{D_{14}^{\rm{ep}}}&{D_{15}^{\rm{ep}}}&{D_{16}^{\rm{ep}}} \\ {D_{21}^{\rm{ep}}}&{D_{22}^{\rm{ep}}}&{D_{23}^{\rm{ep}}}&{D_{24}^{\rm{ep}}}&{D_{25}^{\rm{ep}}}&{D_{26}^{\rm{ep}}} \\ {D_{31}^{\rm{ep}}}&{D_{32}^{\rm{ep}}}&{D_{33}^{\rm{ep}}}&{D_{34}^{\rm{ep}}}&{D_{35}^{\rm{ep}}}&{D_{36}^{\rm{ep}}} \\ {D_{41}^{\rm{ep}}}&{D_{42}^{\rm{ep}}}&{D_{43}^{\rm{ep}}}&{D_{44}^{\rm{ep}}}&{D_{45}^{\rm{ep}}}&{D_{46}^{\rm{ep}}} \\ {D_{51}^{\rm{ep}}}&{D_{52}^{\rm{ep}}}&{D_{53}^{\rm{ep}}}&{D_{54}^{\rm{ep}}}&{D_{55}^{\rm{ep}}}&{D_{56}^{\rm{ep}}} \\ {D_{61}^{\rm{ep}}}&{D_{62}^{\rm{ep}}}&{D_{63}^{\rm{ep}}}&{D_{64}^{\rm{ep}}}&{D_{65}^{\rm{ep}}}&{D_{66}^{\rm{ep}}} \end{array}} \right] \end{split}$

(36)

式(36)中的具体表达式见附录。

3 正常固结非饱和江西红黏土试验

采用江西南昌地区非饱和重塑红黏土，通过室内土工试验得到试验土样的相对重度、最大干密度、最优含水率、液限和塑限等土工参数。制备压实度分别为80%、85%和90%三种土样，在不同基质吸力条件下进行压缩回弹试验，得到本文试验土样的压缩指数和回弹指数。通过土水特征曲线试验得到本文试验土样的土水特征曲线方程，由土水特征曲线方程计算饱和度结果为：初始压实度为90%的条件下，基质吸力等于100 kPa和300 kPa时所对应的土体的饱和度分别为83.9%和69.9%；初始压实度为85%的条件下，基质吸力等于50 kPa和100 kPa时所对应的土体的饱和度分别为79.3%和72.3%；初始压实度为80%的条件下，基质吸力等于100 kPa和200 kPa时所对应的土体的饱和度分别为63.7%和51.2%。最后在压实度分别为80%、85%和90%；净围压分别为100 kPa、200 kPa和300 kPa；不同基质吸力条件下进行非饱和土固结排水试验。由试验结果得到非饱和土的有效内摩擦角、有效黏聚力和偏应力和轴向应变、体应变和轴向应变间的关系曲线。将上述试验所得模型计算参数汇总表如表1所示。

表 1 模型计算参数汇总表

Table 1. Summary table of model calculation parameters

压实度/(%)	饱和土压缩指数 $\lambda \left( 0 \right)$	非饱和土材料参数 ${\lambda _{\rm{s}}}$	饱和土回弹指数 $\kappa \left( 0 \right)$	非饱和土材料参数 ${\kappa _{\rm{s}}}$	饱和土有效黏聚力 $c'$ /kPa	饱和土有效内摩擦角 $\varphi '$ /(^o)	泊松比 $\nu$	初始孔隙比 ${e_0}$
90	0.0666	0.01930	0.00639	−2.640×10⁻⁶	26.90	31	0.35	0.56
85	0.1195	0.06602	0.01369	−2.767×10⁻⁵	26.76	29	0.35	0.71
80	0.1285	0.04024	0.01588	−2.965×10⁻⁵	25.10	26	0.35	0.81

下载: 导出CSV

| 显示表格

4 本构模型的常规三轴固结排水试验验证

为验证新建立的弹塑性本构模型的正确性，采用Fortran语言编制本文模型的计算程序，对非饱和土固结排水试验进行数值模拟。将等量代换法和坐标平移法下的本构模型计算结果与非饱和土三轴固结排水试验结果进行对比，进一步验证新建立的本构模型的正确性。对比结果如图5和图6所示。

图 5 剪应力模拟值与试验值对比图

Figure 5. Comparison between simulated shear stress and experimental stress

下载: 全尺寸图片幻灯片

图 6 小变形体应变模拟值与试验值对比图

Figure 6. Comparison between simulated and experimental strain values of small deformed body

下载: 全尺寸图片幻灯片

由图5(a)可知，在广义有效应力变量法下，等量代换法与试验数据差值最大不超过8%，坐标平移法与试验数据差值最大不超过12%；由图5(b)可知，等量代换法与试验数据差值最大不超过10%，坐标平移法与试验数据差值最大不超过15%；由图5(c)可知，等量代换法与试验数据差值最大不超过9%，坐标平移法与试验数据差值最大不超过14%；由图5(d)可知，等量代换法与试验数据差值最大不超过14%，坐标平移法与试验数据差值最大不超过18%。总体来看，两种方法模拟数据与试验结果均符合较好。不同参数下所有的偏应力模拟值与试验值对比图有着类似的变化规律。本构模型模拟曲线和试验点的变化规律是一致的，说明该本构模型能够较好地反映土体在固结排水状态下的强度特性。同时本构模型计算值与试验值比较接近，对试验结果取得了较好的模拟预测，进一步验证了该模型的正确性。等量代换法模拟值和坐标平移法模拟值在轴向应变较小时比较接近，随着轴应变的增大，二者之间的差值逐渐增大，最终等量代换法计算结果大于坐标平移法计算结果。在轴向应变超过4%以后，坐标平移法的模拟效果相对较差。采用等量代换法计算结果比坐标平移法结果更加精确。土样在剪切的过程中，轴向应变较小时，偏应力随着轴向变形的增大而迅速增大，直至后期逐渐达到稳定值后缓慢增长。在初始压实度和基质吸力相同的条件下，净围压越大，土样的抗剪强度越大；在初始压实度和净围压相同的条件下，基质吸力增大，土样抗剪强度增大；在净围压和基质吸力相同的情况下，初始压实度越大，土样的抗剪强度也会越大。导致计算结果产生误差的原因可能与试验结果误差有关，也有可能与模型中的土样参数值存在误差有关。

由图6可知，本构模型模拟曲线和试验点的变化规律是一致的，说明该本构模型能够较好地反映体应变和轴向应变之间的变化规律。同时该模型计算的体应变值接近于试验值，对试验取得了较好的模拟预测结果，进一步验证了该模型的正确性。在轴向应变较小时，等量代换法计算值和坐标平移法计算值比较接近，在体应变逐渐达到平稳状态时，最终等量代换法计算值大于坐标平移法计算值，等量代换法计算值更符合试验值结果。在土体剪切过程中，体应变在轴向变形较小时增长迅速，直至后期逐渐接近平稳状态时增长缓慢。在初始压实度和基质吸力相同的条件下，净围压越大，土样达到稳定时的体应变值越大；在净围压和基质吸力相同的情况下，初始压实度越大，土样达到稳定时的体应变值越小。

5 固结排水条件下真三轴试验数值模拟

为了反映真实的土体受力情况，本节用Fortran语言根据本文所建立的本构模型，采用等量代换法和坐标平移法对长方体试样进行真三轴非饱和固结排水试验条件下的数值模拟。本节所采用的参数仍旧是江西红黏土的相关土性参数，试验应力加载路径为先对试样三个主应力方向上都施加最小主应力固结，等第一阶段固结完成以后逐步增大其中一个水平方向上的主应力增大至中间主应力大小，直至固结完成。固结阶段完成以后保持中间主应力大小和最小主应力大小不变，逐步增大轴向应力，使试样开始剪切，直至产生15%轴向应变发生剪切破坏。其中，中间主应力影响系数分别取0、0.25、0.50、0.75、1.00，最小主应力取200 kPa，最大主应力取300 kPa。真三轴试样尺寸为： $140\;{\rm{mm}} \times 70\;{\rm{mm}} \times 70\;{\rm{mm}}$ 。在真三轴固结排水试验条件下，本构模型模拟所得到的偏应力与轴向应变关系曲线、体应变与轴向应变关系曲线对比图如图7~图10所示。

图 7 压实度90%、基质吸力200 kPa偏应力和轴向应变关系图

Figure 7. Relationship between the 200 kPa deviant stress and axial strain of the matric suction with 90% compactness

下载: 全尺寸图片幻灯片

图 8 压实度85%、基质吸力200 kPa偏应力和轴向应变关系图

Figure 8. Relationship between the 200 kPa deviant stress and axial strain of the matric suction with 85% compactness

下载: 全尺寸图片幻灯片

图 9 压实度85%、基质吸力200 kPa体应变和轴向应变关系图

Figure 9. Relationship between 200 kPa volume strain and axial strain of 85% compacted matric suction

下载: 全尺寸图片幻灯片

图 10 压实度80%、基质吸力200 kPa体应变和轴向应变关系图

Figure 10. Relationship between 200 kPa volume strain and axial strain of 80% compacted matric suction

下载: 全尺寸图片幻灯片

由图7~图10可知，在真三轴固结排水条件下，中间主应力影响系数对偏应力和体应变都产生影响。在其他条件都相同的情况下，中主应力影响系数值越大，所对应的偏应力和体应变就越大，表明中主应力影响系数可以提高土体的抗剪强度。剪切初期，在轴向应变较小时，不同b值模拟的曲线比较接近，此时偏应力和体应变增长较快；随着轴向应变的增加，不同b值计算的偏应力和体应变值间的差值逐渐增大，曲线增长速度减缓逐渐达到平稳状态。

6 结论

本文将非饱和土广义有效应力原理与三剪强度准则相结合，提出了非饱和土广义有效应力三剪强度准则，并推导出了能够反映土体全应力状态的破坏应力比 $M\left( {p',\theta ,s,{s_{\rm{r}}}} \right)$ 和 $M\left( \theta \right)$ ，将新推导的破坏应力比引入到非饱和土修正剑桥模型中，建立了新的非饱和土三剪统一屈服函数。在弹塑性理论的框架下，建立了基于广义有效应力变量法的非饱和土三剪统一弹塑性本构模型。将本构模型计算结果与非饱和土三轴固结排水试验结果进行了对比验证。结论如下：

(1) 新的强度准则能够很好地反映土体的真实情况，有效地拓宽了非饱和土破坏强度理论，同时它也继承了三剪强度准则的优势，是一系列统一破坏准则，通过取不同的b值可对其他准则作非线性逼近。相比于原来破坏应力比为定值的修正剑桥模型中的屈服函数，新的屈服函数可以更好地反映土体全应力状态、中间主应力效应、拉压不等效应。

(2) 新的本构模型该模型适用于正常固结非饱和黏性土，能够体现中间主应力、黏聚力、基质吸力和饱和度对非饱和土强度和变形特性的影响。可以反映单调静荷载条件下非饱和土体的强度特性、变形特性、屈服特性、压硬性、剪缩特性。

(3) 将数值计算结果与固结排水试验结果进行了对比验证，结果表明：模拟曲线和试验点的变化规律是一致的，且计算值与试验值比较接近，说明新的本构模型能够较好地反映土体在固结排水状态下的强度和变形特性，验证了新的本构模型的正确性。在轴向变形较小时，等量代换法和坐标平移法模拟结果比较接近，随着轴向变形逐渐增大直至剪应力达到平稳状态的过程中，等量代换法计算结果要大于坐标平移法模型计算结果，且更接近于试验值，说明等量代换法能更好地反映土体强度和变形特性，而坐标平移法模拟预测效果相对较差。

附录： ${\;\;}$

${{\boldsymbol{D}}^{\rm{e}}} = \left[ {\begin{array}{*{20}{c}} {L + 2G}&L&L&0&0&0 \\ L&{L + 2G}&L&0&0&0 \\ L&L&{L + 2G}&0&0&0 \\ 0&0&0&G&0&0 \\ 0&0&0&0&G&0 \\ 0&0&0&0&0&G \end{array}} \right],$

$\small D_{11}^{{\rm{ep}}} = \left( {L + 2G} \right) - \frac{{{{\left( {L + 2G} \right)}^2}N_1^2 + {L^2}N_2^2 + {L^2}N_3^2 + 2L\left( {L + 2G} \right){N_1}{N_2} + 2{L^2}{N_2}{N_3} + 2L\left( {L + 2G} \right){N_1}{N_3}}}{Q},$

$\small D_{12}^{{\rm{ep}}} = D_{21}^{{\rm{ep}}} = L - \frac{ L( {L + 2G} )N_1^2 + L( {L + 2G} )N_2^2 + {L^2}N_3^2 + [ {{L^2} + {{( {L + 2G})}^2}} ]{N_1}{N_2} + 2L\left( {L + 2G} ){N_2}{N_3} + 2L( {L + G} \right){N_1}{N_3} }{Q} ,$

$\small { D_{13}^{{\rm{ep}}} = D_{31}^{{\rm{ep}}} = L - \frac{{ L( {L + 2G} )N_1^2 + {L^2}N_2^2 + L( {L + 2G} )N_3^2 + 2L( {L + G} ){N_1}{N_2} + 2L( {L + 2G} ){N_2}{N_3} + [ {{L^2} + {{( {L + 2G} )}^2}} ]{N_1}{N_3} }}{Q} },$

$\small D_{14}^{{\rm{ep}}} = D_{41}^{{\rm{ep}}} = - \frac{{\left( {L + 2G} \right)G{N_1}{N_4} + LG{N_2}{N_4} + LG{N_3}{N_4}}}{Q} ,\; D_{15}^{{\rm{ep}}} = D_{51}^{{\rm{ep}}} = - \frac{{\left( {L + 2G} \right)G{N_1}{N_5} + LG{N_2}{N_5} + LG{N_3}{N_5}}}{Q},$

$\small D_{16}^{{\rm{ep}}} = D_{61}^{{\rm{ep}}} = - \frac{{\left( {L + 2G} \right)G{N_1}{N_6} + LG{N_2}{N_6} + LG{N_3}{N_6}}}{Q} ,$

$\small D_{22}^{{\rm{ep}}} = \left( {L + 2G} \right) - \frac{{{L^2}N_1^2 + {{\left( {L + 2G} \right)}^2}N_2^2 + {L^2}N_3^2 + 2L\left( {L + 2G} \right){N_1}{N_2} + 2L\left( {L + 2G} \right){N_2}{N_3} + 2{L^2}{N_1}{N_3}}}{Q},$

$\small D_{23}^{{\rm{ep}}} = D_{32}^{{\rm{ep}}} = L - \frac{{ {L^2}N_1^2 + L( {L + 2G} )N_2^2 + L( {L + 2G} )N_3^2 + 2L( {L + G} ){N_1}{N_2} + [ {{L^2} + {{( {L + 2G} )}^2}} ]{N_2}{N_3} + 2L( {L + G} ){N_1}{N_3} }}{Q} ,$

$\small D_{24}^{{\rm{ep}}} = D_{42}^{{\rm{ep}}} = - \frac{{LG{N_1}{N_4} + ( {L + 2G} )G{N_2}{N_4} + LG{N_3}{N_4}}}{Q} ,\; D_{25}^{{\rm{ep}}} = D_{52}^{{\rm{ep}}} = - \frac{{LG{N_1}{N_5} + ( {L + 2G} )G{N_2}{N_5} + LG{N_3}{N_5}}}{Q},$

$\small D_{26}^{{\rm{ep}}} = D_{62}^{{\rm{ep}}} = - \frac{{LG{N_1}{N_6} + ( {L + 2G} )G{N_2}{N_6} + LG{N_3}{N_6}}}{Q},$

$\small D_{33}^{{\rm{ep}}} = \left( {L + 2G} \right) - \frac{{{L^2}N_1^2 + {L^2}N_2^2 + {{\left( {L + 2G} \right)}^2}N_3^2 + 2{L^2}{N_1}{N_2} + 2L\left( {L + 2G} \right){N_2}{N_3} + 2L\left( {L + 2G} \right){N_1}{N_3}}}{Q},$

$\small D_{34}^{{\rm{ep}}} = D_{43}^{{\rm{ep}}} = - \frac{{LG{N_1}{N_4} + LG{N_2}{N_4} + \left( {L + 2G} \right)G{N_3}{N_4}}}{Q} ,\; D_{35}^{{\rm{ep}}} = D_{53}^{{\rm{ep}}} = - \frac{{LG{N_1}{N_5} + LG{N_2}{N_5} + \left( {L + 2G} \right)G{N_3}{N_5}}}{Q} ,$

$\small D_{36}^{{\rm{ep}}} = D_{63}^{{\rm{ep}}} = - \frac{{LG{N_1}{N_6} + LG{N_2}{N_6} + \left( {L + 2G} \right)G{N_3}{N_6}}}{Q} ,D_{44}^{\rm {ep}} = G - \frac{{{G^2}N_4^2}}{Q} ,\; D_{55}^{\rm {ep}} = G - \frac{{{G^2}N_5^2}}{Q} ,\;$

$D_{45}^{\rm {ep}} = D_{54}^{\rm {ep}} = - \frac{{{G^2}{N_4}{N_5}}}{Q} ,\; D_{56}^{\rm {ep}} = D_{65}^{\rm {ep}} = - \frac{{{G^2}{N_5}{N_6}}}{Q} ,\; D_{46}^{\rm {ep}} = D_{64}^{\rm {ep}} = - \frac{{{G^2}{N_4}{N_6}}}{Q} ,\; D_{66}^{\rm {ep}} = G - \frac{{{G^2}N_6^2}}{Q},$

$\small Q = {K_{\rm{p}}} + [ {( {L + 2G} )( {N_1^2 + N_2^2 + N_3^2} ) + G( {N_4^2 + N_5^2 + N_6^2} ) + 2L( {{N_1}{N_2} + {N_2}{N_3} + {N_3}{N_1}} )} ],$

$\small {N_1} = \frac{{\partial f}}{{\partial {I'_1}}} + 2\frac{{\partial f}}{{\partial {J_2}}}{S_1} + \frac{{\partial f}}{{\partial \theta }}\frac{{\partial \theta }}{{\partial {J_3}}}\left( {\frac{2}{3}{S_2}{S_3} - \frac{1}{3}{S_1}{S_3} - \frac{1}{3}{S_1}{S_2} - \frac{2}{3}\tau _{23}^2 + \frac{1}{3}\tau _{31}^2 + \frac{1}{3}\tau _{12}^2} \right),$

$\small {N_2} = \frac{{\partial f}}{{\partial {I'_1}}} + 2\frac{{\partial f}}{{\partial {J_2}}}{S_2} + \frac{{\partial f}}{{\partial \theta }}\frac{{\partial \theta }}{{\partial {J_3}}}\left( { - \frac{1}{3}{S_2}{S_3} + \frac{2}{3}{S_1}{S_3} - \frac{1}{3}{S_1}{S_2} + \frac{1}{3}\tau _{23}^2 - \frac{2}{3}\tau _{31}^2 + \frac{1}{3}\tau _{12}^2} \right),$

$\small {N_3} = \frac{{\partial f}}{{\partial {I'_1}}} + 2\frac{{\partial f}}{{\partial {J_2}}}{S_3} + \frac{{\partial f}}{{\partial \theta }}\frac{{\partial \theta }}{{\partial {J_3}}}\left( { - \frac{1}{3}{S_2}{S_3} - \frac{1}{3}{S_1}{S_3} + \frac{2}{3}{S_1}{S_2} + \frac{1}{3}\tau _{23}^2 + \frac{1}{3}\tau _{31}^2 - \frac{2}{3}\tau _{12}^2} \right) ,$

$\small {N_4} = 4\frac{{\partial f}}{{\partial {J_2}}}{\tau _{12}} + \frac{{\partial f}}{{\partial \theta }}\frac{{\partial \theta }}{{\partial {J_3}}}\left[ {2\left( {{\tau _{23}}{\tau _{31}} - {S_3}{\tau _{12}}} \right) + \frac{2}{3}{\tau _{12}}} \right] ,\; {N_5} = 4\frac{{\partial f}}{{\partial {J_2}}}{\tau _{23}} + \frac{{\partial f}}{{\partial \theta }}\frac{{\partial \theta }}{{\partial {J_3}}}\left[ {2\left( {{\tau _{31}}{\tau _{12}} - {S_3}{\tau _{23}}} \right) + \frac{2}{3}{\tau _{23}}} \right] ,$

$\small {N_6} = 4\frac{{\partial f}}{{\partial {J_2}}}{\tau _{31}} + \frac{{\partial f}}{{\partial \theta }}\frac{{\partial \theta }}{{\partial {J_3}}}\left[ {2\left( {{\tau _{12}}{\tau _{23}} - {S_3}{\tau _{31}}} \right) + \frac{2}{3}{\tau _{31}}} \right] , \frac{{\partial f}}{{\partial {I'_1}}} = \frac{{\partial f}}{{\partial p'}}\frac{{\partial p'}}{{\partial {I'_1}}} = \frac{1}{3}\frac{{\partial f}}{{\partial p'}} ,\; \frac{{\partial f}}{{\partial {J_2}}} = \frac{{\partial f}}{{\partial q}}\frac{{\partial q}}{{\partial {J_2}}} = \frac{3}{{2q}}\frac{{\partial f}}{{\partial q}} 。$

$\small \dfrac{{\partial f}}{{\partial q}}$ 和 $\small \dfrac{{\partial f}}{{\partial \theta }}$ 的计算式为：

1) 使用等量代换法时，

$\small \frac{{\partial f}}{{\partial q}} = \frac{{2q}}{{M{{\left( {p',\theta ,s,{s_{\rm{r}}}} \right)}^2}}} ,\; \frac{{\partial f}}{{\partial \theta }} = - 2\frac{{{q^2}}}{{M{{\left( {p',\theta ,s,{s_{\rm{r}}}} \right)}^3}}}\frac{{p'\sin \varphi ' + c'\cos \varphi ' + {s_{\rm{r}}}s\sin \varphi '}}{{p'}}\frac{{\partial A}}{{\partial \theta }}。$

2) 使用坐标平移法时，

$\small \frac{{\partial f}}{{\partial q}} = \frac{{2q}}{{M{{\left( \theta \right)}^2}}} ,\; \frac{{\partial f}}{{\partial \theta }} = - 2\frac{{{q^2}}}{{M{{\left( \theta \right)}^3}}}\frac{{\partial A}}{{\partial \theta }}\sin \varphi ',$

$\small \dfrac{{\partial A}}{{\partial \theta }} = \dfrac{ \begin{gathered}\Bigg\{ - 6\left( {1 + b} \right)\sin \left( {\theta - \dfrac{\pi }{6}} \right)\left[ 2\sqrt 3 {\cos ^2}\left( {\theta - \dfrac{\pi }{6}} \right) + 2\sqrt 3 b{\cos ^2}\left( {\theta + \dfrac{\pi }{6}} \right) + 2\sqrt 3 b{\sin ^2}\left( \theta \right) - \left( {1 + b} \right)\cos \left( {2\theta + \dfrac{\pi }{6}} \right)\sin \varphi ' \right] \\ - 6\left( {1 + b} \right)\cos \left( {\theta - \dfrac{\pi }{6}} \right)\left[ - 2\sqrt 3 \sin \left( {2\theta - \dfrac{\pi }{3}} \right) - 2\sqrt 3 b\sin \left( {2\theta + \dfrac{\pi }{3}} \right) + 2\sqrt 3 b\sin \left( {2\theta } \right) + 2\left( {1 + b} \right)\sin \left( {2\theta + \dfrac{\pi }{6}} \right)\sin \varphi ' \right] \Bigg\} \end{gathered} } {{{{\left[ {2\sqrt 3 {{\cos }^2}\left( {\theta - \dfrac{\pi }{6}} \right) + 2\sqrt 3 b{{\cos }^2}\left( {\theta + \dfrac{\pi }{6}} \right) + 2\sqrt 3 b{{\sin }^2}\left( \theta \right) - \left( {1 + b} \right)\cos \left( {2\theta + \dfrac{\pi }{6}} \right)\sin \varphi '} \right]}^2}}} ,$

$\small \frac{{\partial \theta }}{{\partial {J_3}}} = - \frac{{\sqrt 3 }}{{2\sin \theta }}\frac{1}{{J_2^{3/ 2}}} ,\; {I'_1} = {\sigma '_1} + {\sigma '_2} + {\sigma '_3} ,\; {J_2} = \frac{1}{6}[ {{{( {{\sigma '_1} - {\sigma '_2}} )}^2} + {{( {{\sigma '_2} - {\sigma '_3}} )}^2} + {{( {{\sigma '_3} - {\sigma '_1}} )}^2}} ],$

$\small {J_3} = \frac{1}{{27}}( {2{{\sigma '}_1} - {\sigma '_2} - {\sigma '_3}} )( {2{{\sigma '}_2} - {\sigma '_3} - {\sigma '_1}} )( {2{\sigma '_3} - {\sigma '_1} - {\sigma '_2}} ) 。$

$\small {B_1}$ 的计算式为：

1) 使用等量代换法时，

$\small {B_1} = \frac{{1 + {e_0}}}{\kappa }p'$

2) 使用坐标平移法时，

$\small {B_1} = \frac{{1 + {e_0}}}{{{\varLambda _1}\kappa }}\left( {p' + c'\cot \varphi ' + {s_{\rm{r}}}s} \right)\;,$

$\small G = \frac{{3\left( {1 - 2\nu } \right)}}{{2\left( {1 + \nu } \right)}}{B_1} ,\; L = {B_1} - \frac{2}{3}G ,\; p' = \frac{{{{\sigma '_1}} + {\sigma '_2} + {\sigma '_3}}}{3} ,\; q = \frac{1}{{\sqrt 2 }}\sqrt {{{( {{\sigma '_1} - {\sigma '_2}} )}^2} + {{( {{\sigma '_2} - {\sigma '_3}} )}^2} + {{( {{\sigma '_3} - {\sigma '_1}} )}^2}},$

$\small \frac{\partial {{I}^{\prime }}_{1}}{\left\{\partial \sigma \right\}}\text={\left[1,1,1,0,0,0\right]}^{{\rm{T}}} ,\; \frac{{\partial {J_2}}}{{\left\{ {\partial \sigma } \right\}}} = {\left[ {{S_1},{S_2},{S_3},2{\tau _{23}},2{\tau _{31}},2{\tau _{12}}} \right]^{\rm{T}}} ,$

$\small \dfrac{{\partial {J_3}}}{{\left\{ {\partial \sigma } \right\}}} = \left[ {\begin{array}{*{20}{c}} {\dfrac{2}{3}{S_2}{S_3} - \dfrac{1}{3}{S_1}{S_3} - \dfrac{1}{3}{S_1}{S_2} - \dfrac{2}{3}\tau _{23}^2 + \dfrac{1}{3}\tau _{31}^2 + \dfrac{1}{3}\tau _{12}^2} \\ { - \dfrac{1}{3}{S_2}{S_3} + \dfrac{2}{3}{S_1}{S_3} - \dfrac{1}{3}{S_1}{S_2} + \dfrac{1}{3}\tau _{23}^2 - \dfrac{2}{3}\tau _{31}^2 + \dfrac{1}{3}\tau _{12}^2} \\ { - \dfrac{1}{3}{S_2}{S_3} - \dfrac{1}{3}{S_1}{S_3} + \dfrac{2}{3}{S_1}{S_2} + \dfrac{1}{3}\tau _{23}^2 + \dfrac{1}{3}\tau _{31}^2{\text{ - }}\dfrac{2}{3}\tau _{12}^2} \\ {2\left( {{\tau _{23}}{\tau _{12}} - {S_3}{\tau _{12}} + \dfrac{1}{3}{\tau _{12}}} \right)} \\ {2\left( {{\tau _{31}}{\tau _{12}} - {S_3}{\tau _{23}} + \dfrac{1}{3}{\tau _{23}}} \right)} \\ {2\left( {{\tau _{12}}{\tau _{23}} - {S_3}{\tau _{31}} + \dfrac{1}{3}{\tau _{31}}} \right)} \end{array}} \right] ,$

$\small \begin{split}& {S_1} = {\sigma '_1} - \frac{{{\sigma '_1} + {\sigma '_2} + {\sigma '_3}}}{3},\;{S_2} = {\sigma '_2} - \frac{{{\sigma '_1} + {\sigma '_2} + {\sigma '_3}}}{3},\; {S_3} = {\sigma '_3} - \frac{{{\sigma '_1} + {\sigma '_2} + {\sigma '_3}}}{3},\;{\tau _{12}} = \frac{{{\sigma '_1} - {\sigma '_2}}}{2}, \;{\tau _{23}} = \frac{{{{\sigma '}_2} - {\sigma '_3}}}{2},\;{\tau _{31}} = \frac{{{\sigma '_3} - {\sigma '_1}}}{2}。 \end{split}$

图 1 广义有效应力变量法计算结果和试验数据对比图

Figure 1. Comparison diagram of calculation results of generalized effective stress variable method and test data

下载: 全尺寸图片幻灯片

图 2 子午线计算结果与试验数据对比图

Figure 2. Comparison diagram of meridian calculation results and experimental data

下载: 全尺寸图片幻灯片

图 3 广义有效应力变量法计算结果和试验数据对比图

Figure 3. Comparison diagram of calculation results of generalized effective stress variable method and test data

下载: 全尺寸图片幻灯片

图 4 子午线计算结果与试验数据对比图

Figure 4. Comparison diagram of meridian calculation results and experimental data

下载: 全尺寸图片幻灯片

图 5 剪应力模拟值与试验值对比图

Figure 5. Comparison between simulated shear stress and experimental stress

下载: 全尺寸图片幻灯片

图 6 小变形体应变模拟值与试验值对比图

Figure 6. Comparison between simulated and experimental strain values of small deformed body

下载: 全尺寸图片幻灯片

图 7 压实度90%、基质吸力200 kPa偏应力和轴向应变关系图

Figure 7. Relationship between the 200 kPa deviant stress and axial strain of the matric suction with 90% compactness

下载: 全尺寸图片幻灯片

图 8 压实度85%、基质吸力200 kPa偏应力和轴向应变关系图

Figure 8. Relationship between the 200 kPa deviant stress and axial strain of the matric suction with 85% compactness

下载: 全尺寸图片幻灯片

图 9 压实度85%、基质吸力200 kPa体应变和轴向应变关系图

Figure 9. Relationship between 200 kPa volume strain and axial strain of 85% compacted matric suction

下载: 全尺寸图片幻灯片

图 10 压实度80%、基质吸力200 kPa体应变和轴向应变关系图

Figure 10. Relationship between 200 kPa volume strain and axial strain of 80% compacted matric suction

下载: 全尺寸图片幻灯片

表 1 模型计算参数汇总表

Table 1 Summary table of model calculation parameters

压实度/(%)	饱和土压缩指数 $\lambda \left( 0 \right)$	非饱和土材料参数 ${\lambda _{\rm{s}}}$	饱和土回弹指数 $\kappa \left( 0 \right)$	非饱和土材料参数 ${\kappa _{\rm{s}}}$	饱和土有效黏聚力 $c'$ /kPa	饱和土有效内摩擦角 $\varphi '$ /(^o)	泊松比 $\nu$	初始孔隙比 ${e_0}$
90	0.0666	0.01930	0.00639	−2.640×10⁻⁶	26.90	31	0.35	0.56
85	0.1195	0.06602	0.01369	−2.767×10⁻⁵	26.76	29	0.35	0.71
80	0.1285	0.04024	0.01588	−2.965×10⁻⁵	25.10	26	0.35	0.81

下载: 导出CSV

参考文献(32)

[1]	ALONSO E E, GENS A, JOSA A. A constitutive model for partially saturated soils [J]. Géotechnique, 1990, 40(3): 405 − 430.
[2]	SHENG D C, FREDLUND D G, GENS A. A new modelling approach for unsaturated soils using independent stress variables [J]. Canadian Geotechnical Journal, 2008, 45(4): 511 − 534. doi: 10.1139/T07-112
[3]	刘艳. 非饱和土的广义有效应力原理及其本构模型研究[D]. 北京: 北京交通大学, 2010. LIU Yan. Study on the principle of generalized effective stress and constitutive modeling of unsaturated soils [D]. Beijing: Beijing Jiaotong University, 2010. (in Chinese)
[4]	路德春, 罗磊, 王欣, 等. 土与结构接触面土体软/硬化本构模型及数值实现[J]. 工程力学, 2017, 34(7): 41 − 50. doi: 10.6052/j.issn.1000-4750.2016.01.0072 LU Dechun, LUO Lei, WANG Xin, et al. Softening/hardening constitutive model for soil-structure interface and numerical implementation [J]. Engineering Mechanics, 2017, 34(7): 41 − 50. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.01.0072
[5]	邓子千, 陈嘉帅, 王建伟, 等. 基于SFG模型的统一屈服面本构模型与试验研究[J]. 岩土力学, 2020, 41(2): 527 − 534. doi: 10.16285/j.rsm.2019.0165 DENG Ziqian, CHEN Jiashuai, WANG Jianwei, et al. Constitutive model and experimental study of uniform yield surface based on SFG model [J]. Rock and Soil Mechanics, 2020, 41(2): 527 − 534. (in Chinese) doi: 10.16285/j.rsm.2019.0165
[6]	XIONG Y L, YE G L, XIE Y, et al. A unified constitutive model for unsaturated soil under monotonic and cyclic loading [J]. Acta Geotechnica, 2019, 14(2): 313 − 328. doi: 10.1007/s11440-018-0754-2
[7]	刘艳, 赵成刚, 蔡国庆, 等. 非饱和土力学理论的研究进展[J]. 力学与实践, 2015, 37(4): 457 − 465. doi: 10.6052/1000-0879-15-008 LIU Yan, ZHAO Chenggang, CAI Guoqing, et al. Research progress of unsaturated soil mechanics [J]. Mechanics in Engineering, 2015, 37(4): 457 − 465. (in Chinese) doi: 10.6052/1000-0879-15-008
[8]	MATSUOKA H, SUN D A, KOGANE A, et al. Stress–strain behaviour of unsaturated soil in true triaxial tests [J]. Canadian Geotechnical Journal, 2002, 39(3): 608 − 619. doi: 10.1139/t02-031
[9]	张常光, 范文, 赵均海. 非饱和土统一强度理论及真三轴试验结果验证[J]. 岩石力学与工程学报, 2015, 34(8): 1702 − 1711. doi: 10.13722/j.cnki.jrme.2014.1335 ZHANG Changguang, FAN Wen, ZHAO Junhai. Unified strength theory of unsaturated soils and verification with true triaxial test [J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(8): 1702 − 1711. (in Chinese) doi: 10.13722/j.cnki.jrme.2014.1335
[10]	胡小荣, 魏雪英, 俞茂宏. 三轴压缩下岩石强度与破坏面角度的双剪理论分析[J]. 岩石力学与工程学报, 2003, 22(7): 1093 − 1098. doi: 10.3321/j.issn:1000-6915.2003.07.009 HU Xiaorong, WEI Xueying, YU Maohong. Characteristics of rock strength and failure under triaxial compression analyzed by twin shear theory [J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(7): 1093 − 1098. (in Chinese) doi: 10.3321/j.issn:1000-6915.2003.07.009
[11]	胡小荣, 林太清. 三剪屈服准则及其在极限内压计算中的应用[J]. 中国有色金属学报, 2007, 17(2): 207 − 215. doi: 10.3321/j.issn:1004-0609.2007.02.005 HU Xiaorong, LIN Taiqing. Triple-shear yield criterion and its applications to limit pressures for thin and thick wall cylinders and thick wall spherical shells [J]. The Chinese Journal of Nonferrous Metals, 2007, 17(2): 207 − 215. (in Chinese) doi: 10.3321/j.issn:1004-0609.2007.02.005
[12]	BISHOP A W. The principle of effective stress [J]. Tecknisk Ukeblad, 1959, 106(39): 859 − 863.
[13]	FREDLUND D G, MORGENSTERN N R. Stress state variables for unsaturated soils [J]. Journal of the Geotechnical Engineering Division, 1977, 103(5): 447 − 466. doi: 10.1061/AJGEB6.0000423
[14]	BISHOP A W, DONALD I B. The experimental study of partly saturated soils in triaxial apparatus [C]// Proceedings of the 5th International Conference on Soil Mechanics and Foundation Engineering. Paris: Dunod, 1961: 13 − 21.
[15]	KHALILI N, KHABBAZ M H. A unique relationship for χ for the determination of the shear strength of unsaturated soils [J]. Géotechnique, 1998, 48(5): 681 − 687.
[16]	WHEELER S J, SHARMA R S, BUISSON M S R. Coupling of hydraulic hysteresis and stress-strain behaviour in unsaturated soils [J]. Géotechnique, 2003, 53(1): 41 − 54.
[17]	BOCKING K A, FREDLUND D G. Limitations of the axis-translation technique [C]// Proceedings of the 4th International Conference on Expansive Soils. Denver: ASCE, 1980: 117 − 135.
[18]	TARANTINO A. Direct measurement of soil water tension [C]// Proceedings of the 3th International Conference on Unsaturated Soils. Brazil: Balkema Publishers, 2002: 1005 − 1017.
[19]	赵成刚, 蔡国庆. 非饱和土广义有效应力原理[J]. 岩土力学, 2009, 30(11): 3232 − 3236. doi: 10.3969/j.issn.1000-7598.2009.11.002 ZHAO Chenggang, CAI Guoqing. Principle of generalized effective stress for unsaturated soils [J]. Rock and Soil Mechanics, 2009, 30(11): 3232 − 3236. (in Chinese) doi: 10.3969/j.issn.1000-7598.2009.11.002
[20]	赵成刚, 张雪东. 非饱和土中功的表述以及有效应力与相分离原理的讨论[J]. 中国科学: E辑, 2008, 51(9): 1530 − 1541. doi: 10.1007/s11431-008-0107-8 ZHAO Chenggang, ZHANG Xuedong. Derivation of the work expression and discussion on the effective principle and the phase separation theorem in unsaturated soil [J]. Science in China Series E: Technological Sciences, 2008, 51(9): 1530 − 1541. (in Chinese) doi: 10.1007/s11431-008-0107-8
[21]	ZHAO C G, LIU Y, GAO F P. Work and energy equations and the principle of generalized effective stress for unsaturated soils [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2010, 34(9): 920 − 936.
[22]	胡小荣, 俞茂宏. 材料三剪屈服准则研究[J]. 工程力学, 2006, 23(4): 6 − 11. doi: 10.3969/j.issn.1000-4750.2006.04.002 HU Xiaorong, YU Maohong. Research on triple-shear yield criterion for materials [J]. Engineering Mechanics, 2006, 23(4): 6 − 11. (in Chinese) doi: 10.3969/j.issn.1000-4750.2006.04.002
[23]	胡小荣, 俞茂宏. 岩土类介质强度准则新探[J]. 岩石力学与工程学报, 2004, 23(18): 3037 − 3043. doi: 10.3321/j.issn:1000-6915.2004.18.001 HU Xiaorong, YU Maohong. New research on failure criterion for geomaterial [J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(18): 3037 − 3043. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.18.001
[24]	胡小荣, 董肖龙, 陈晓宇. 正常固结原状饱和黏性土的本构模型研究[J]. 应用力学学报, 2019, 36(1): 36 − 46. HU Xiaorong, DONG Xiaolong, CHEN Xiaoyu. Constitutive model for undisturbed saturated clays in normal consolidation [J]. Chinese Journal of Applied Mechanics, 2019, 36(1): 36 − 46. (in Chinese)
[25]	胡小荣, 蔡晓锋, 李春博, 等. 正常固结非饱和黏性土的三剪弹塑性本构模型[J]. 工程力学, 2022, 39(1): 164 − 174. doi: 10.6052/j.issn.1000-4750.2020.12.0910 HU Xiaorong, CAI Xiaofeng, LI Chunbo, et al. Triple-shear elasto-plastic constitutive models of unsaturated clays in normal consolidation [J]. Engineering Mechanics, 2022, 39(1): 164 − 174. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.12.0910
[26]	陈明祥. 弹塑性力学[M]. 北京: 科学出版社, 2007. CHEN Mingxiang. Elasticity and mechanics [M]. Beijing: Science Press, 2007. (in Chinese)
[27]	罗汀, 姚仰平, 侯伟. 土的本构关系[M]. 北京: 人民交通出版社, 2010. LUO Ting, YAO Yangping, HOU Wei. Soil constitutive models [M]. Beijing: China Communications Press, 2010. (in Chinese)
[28]	王娇. 非饱和黄土土水特性与真三轴屈服及强度变形特性试验研究[D]. 西安: 西安理工大学, 2018. WANG Jiao. Research on soil water characteristics and yielding, strength and deformation properties of unsaturated loess by the true triaxial test [D]. Xi’an: Xi’an University of Technology, 2018. (in Chinese)
[29]	HOYOS L R, PÉREZ-RUIZ D D, PUPPALA A J. Modeling unsaturated soil response under suction-controlled true triaxial stress paths [J]. International Journal of Geomechanics, 2012, 12(3): 292 − 308. doi: 10.1061/(ASCE)GM.1943-5622.0000159
[30]	HOYOS L R, PÉREZ-RUIZ D D, PUPPALA A J. Refined true triaxial apparatus for testing unsaturated soils under suction-controlled stress paths [J]. International Journal of Geomechanics, 2012, 12(3): 281 − 291. doi: 10.1061/(ASCE)GM.1943-5622.0000138
[31]	孙德安. 非饱和土的水力和力学特性及其弹塑性描述[J]. 岩土力学, 2009, 30(11): 3217 − 3231. doi: 10.3969/j.issn.1000-7598.2009.11.001 SUN Dean. Hydro-mechanical behaviours of unsaturated soils and their elastoplastic modelling [J]. Rock and Soil Mechanics, 2009, 30(11): 3217 − 3231. (in Chinese) doi: 10.3969/j.issn.1000-7598.2009.11.001
[32]	胡小荣, 蔡晓锋, 陈昊. 正常固结非饱和黏性土的三剪边界面模型研究[J]. 工程力学, 2021, 38(6): 175 − 190. doi: 10.6052/j.issn.1000-4750.2020.07.0462 HU Xiaorong, CAI Xiaofeng, CHEN Hao. Triple-shear bounding surface models for unsaturated clays in normal consolidation [J]. Engineering Mechanics, 2021, 38(6): 175 − 190. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.07.0462

施引文献(3)

期刊类型引用(2)

1.	姚兆明，唐赛，昌语，李鹏辉. 冻结改良土抗压特性分数阶模型可靠性分析. 河南城建学院学报. 2025(01): 44-51 . 百度学术
2.	曹雪叶，张二媛，孙杰龙，邱明明. 降雨入渗条件下非饱和朗肯土压力三剪统一解. 延安大学学报(自然科学版). 2025(01): 20-28 . 百度学术