考虑初始状态影响的膨胀土一维膨胀特性研究

叶云雪; 邹维列; 韩仲; 谢鹏; 张俊峰; 徐永福

doi:10.11779/CJGE202108017

考虑初始状态影响的膨胀土一维膨胀特性研究 English Version

叶云雪^{1, 2,},
邹维列²,
韩仲^2, ,,
谢鹏³,
张俊峰⁴,
徐永福⁵

1.
南昌航空大学土木建筑学院,江西　南昌　330063
2.
武汉大学土木建筑工程学院,湖北　武汉　430072
3.
广东省建筑科学研究院集团股份有限公司,广东　广州　510500
4.
中建三局集团有限公司,湖北　武汉　430074
5.
上海交通大学船舶海洋与建筑工程学院,上海　200240

基金项目:

国家重点研发计划项目 2019YFC1509800

国家自然科学基金项目 51779191

国家自然科学基金项目 51809199

南昌航空大学博士启动基金项目 EA201911396

详细信息

作者简介:
叶云雪(1988— ),男,博士,讲师,主要从事非饱和土方面的科研和教学工作。E-mail：yeyunxue@163.com

通讯作者:
韩仲, E-mail：zhong.han@whu.edu.cn

中图分类号: TU43
计量
- 文章访问数: 267
- HTML全文浏览量: 39
- PDF下载量: 130
出版历程
- 收稿日期: 2020-12-24
- 网络出版日期: 2022-12-02
- 刊出日期: 2021-07-31

One-dimensional swelling characteristics of expansive soils considering influence of initial states

1.
School of Civil Engineering and Architecture, Nanchang Hangkong University, Nanchang 330063, China
2.
School of Civil Engineering, Wuhan University, Wuhan 430072, China
3.
Guangdong Provincial Academy of Building Research Group Co., Ltd., Guangzhou 510500, China
4.
China Construction Third Bureau Co., Ltd., Wuhan 430074, China
5.
School of Naval Architecture, Ocean & Civil Engineering, Shanghai Jiaotong University, Shanghai 200240, China

摘要

摘要: 以南阳中膨胀土为研究对象,通过室内一维膨胀特性试验和土水特性试验探究不同初始状态（包括吸力、含水率和干密度等）影响下非饱和重塑土膨胀率或膨胀力的演化规律。试验结果表明,在不同初始含水率下,无荷膨胀率与初始吸力呈幂函数关系,且与初始干密度无关；膨胀力随初始干密度的增加呈幂函数增大,而膨胀力随初始含水率的增加却是先增大后减小；膨胀力与初始干密度和最终吸水量均呈幂函数关系,给出初始干密度–最终吸水量–膨胀力的三变量模型；构建了全吸力范围内初始含水率–初始干密度–初始吸力–膨胀力的统一模型；所提模型不仅参数少,方程左右量纲统一,且与试验数据具有良好的一致性,便于工程应用。
- 膨胀土 /
- 膨胀特性 /
- 土水特性 /
- 初始状态 /
- 模型
Abstract: Taking the Nanyang medium expansive soil as the research object, the evolution laws of the swelling force and swelling rate of unsaturated remolded soils under different initial conditions (including suction, water content and dry density, etc.) are explored through indoor one-dimensional swelling characteristic tests and soil-water characteristic tests. The test results show that under different initial water contents, the unloaded expansion rate has a power function relationship with the initial suction, and this function relationship has nothing to do with the initial dry density. The swelling force increases as a power function with the increase of the initial dry density, while the swelling force first increases and then decreases with the increase of the initial water content. The swelling force has a power function relationship with the initial dry density and final water absorption. A three-variable model for the initial dry density-final water absorption-swelling force is given. A unified model for the initial water content-initial dry density-initial suction-swelling force within the full suction range is established. The proposed model has few parameters and unified left and right dimensions of the equation, and it has good consistency with the experimental data, which is convenient for engineering application. The research results provide a reliable theoretical method for the analysis of the swelling characteristics and soil-water characteristics of the unsaturated Nanyang expansive soil.
- expansive soil /
- swelling property /
- soil-water characteristic /
- initial state /
- model

HTML全文

0. 引言

天然土体颗粒间通常具有胶结物质,比如碳酸钙、微生物、黏粒、甲烷水合物等,称为结构性土。结构性土可看成是一种特殊的胶结颗粒材料,颗粒间胶结物质的存在使结构性土的宏观力学性质不同于重塑土^[1-2]。为了分析结构性土工程问题,需要建立合理的结构性土本构模型。在临界状态土力学基础上,相关学者引入状态变量,分别建立了边界面模型^[3]、下负荷面模型^[4]、统一硬化模型^[5]等重塑土本构模型。在加荷过程中,结构性土颗粒间的胶结逐渐破坏并退出承担应力,土体向重塑土转化。因此,在重塑土本构模型基础上,通过引入表征结构性的参量及其演化规律建立结构性土本构模型是一条可行的研究途径。不少学者提出了结构性土本构模型,在模型中假定了不同的结构性表征参量。比如：①通过比较结构性土屈服面大小 ${\bar{p}}_{c}$ 和重塑土屈服面大小 $p_{c}$ 的不同,定义 $Δ p_{c} = {\bar{p}}_{c} - p_{c}$ ^[6]或 $R^{*} = p_{c} /$ ${\bar{p}}_{c}$ ^[7]为结构性表征参量；②通过引入胶结组分（相对完整状态）和摩擦组分（完全调整状态）应力分担参量（扰动参量、结构破损参量）D^[8]或B₀^[9-10]表达结构性衰退程度。由于胶结破损演化难以通过试验测得,结构性表征参量的演化通常需要假设,缺乏微观试验结果验证。

离散单元法^[11]模拟结构性土弥补了室内微观测试在获取胶结破损演化方面的不足,已广泛应用于结构性土体宏微观力学性质研究。Jiang等^[10]通过开展结构性土的二维离散元模拟,对3种结构性土本构模型结构性参量演化假设的合理性进行了验证,基于离散元结果提出了新的胶结破损演化式,并应用到本构模型研究^[12],不足之处在于胶结破损演化式中的参数受加载路径影响。通过PFC软件植入完整胶结接触模型（考虑了颗粒间抗转动和抗扭转能力及胶结尺寸对刚度和强度的影响）建立了结构性土三维离散元试样,开展了侧限压缩、等向压缩、等应力比压缩以及常规三轴和真三轴试验的离散元数值模拟,分析了结构破损参量B₀随等效塑性应变的演化规律,并提出了应力路径相关性低的结构破损参量B_σ。

1. 胶结接触模型

离散元模拟采用软硬复合胶结接触模型框架^[13],考虑了颗粒间抗转动和抗扭转能力^[14]及胶结尺寸对刚度和强度的影响^[15]。对于通过胶结间接传力/力矩的接触（“厚胶结”）直接采用胶结接触公式计算接触力/力矩；对于颗粒直接接触并有胶结传力/力矩的接触（“薄胶结”）,按并联模式通过颗粒接触力/力矩和胶结接触力/力矩之和计算接触力/力矩。

1.1 颗粒接触

假设颗粒间通过一圆形截面接触,半径为

$R_{c} = β R$ ,

(1)

式中, $β$ 为抗转动系数, $R = 2 R_{1} R_{2} / (R_{1} + R_{2})$ , $R_{1}$ , $R_{2}$ 为两颗粒半径。

接触间力/力矩的传递通过下列公式计算：

$F_{n}^{l} = {\begin{cases} k_{n} δ_{n} - F_{a} (δ_{n} \geq 0) \\ 0 (δ_{n} < 0) \end{cases}$ ,

(2)

$F_{s}^{l} = {(F_{s}^{l})}_{o} - k_{s} Δ δ_{s}$ ,

(3)

$M_{r}^{l} = {(M_{r}^{l})}_{o} - k_{r} Δ θ_{r}$ ,

(4)

$M_{t}^{l} = {(M_{t}^{l})}_{o} - k_{t} Δ θ_{t}$ 。

(5)

式中 $F_{a}$ 为粒间引力,比如毛细力、范德华力、电子引力等； $F_{n}^{l}$ , $F_{s}^{l}$ , $M_{r}^{l}$ , $M_{t}^{l}$ 分别为颗粒法向力、切向力、弯矩和扭矩； $k_{n}$ , $k_{s}$ ,k_r,k_t分别为颗粒法向、切向、抗弯、抗扭刚度, $k_{r} = 0.25 k_{n} R_{c}^{2}$ , $k_{t} = 0.5 k_{s} R_{c}^{2}$ ； $δ_{n}$ 为颗粒重叠量, $Δ δ_{s}$ , $Δ θ_{r}$ 和 $Δ θ_{t}$ 分别为剪切位移、转动角和扭转角增量。

颗粒接触的法切向刚度可按PFC软件手册^[16]按下式计算：

$k_{n} = A E^{*} / L$ ,

(6)

$k_{s} = k_{n} / κ^{*}$ ,

(7)

式中, $A = π R^{2}$ , $E^{*}$ 是颗粒接触的等效模量, $κ^{*}$ 是颗粒接触的法切向刚度比, $L$ 为两颗粒半径之和。

接触抗剪强度、抗弯强度和抗扭强度分别为^[14-15]

$‖ F_{s}^{l} ‖ \leq μ (F_{n} + F_{a})$ ,

(8)

$‖ M_{r}^{l} ‖ \leq 0.25 ζ_{c} R_{c} (F_{n} + F_{a})$ ,

(9)

$M_{t}^{l} \leq 0.65 μ R_{c} (F_{n} + F_{a})$ ,

(10)

式中, $μ$ 是摩擦系数, $ζ_{c}$ 取2.1。颗粒接触力和力矩达到以上强度后颗粒将进入滑移、转动和扭转状态。

1.2 胶结计算

假定胶结物质呈圆形截面连接两颗粒,则胶结法向力 $F_{n}^{b}$ ,切向力 $F_{s}^{b}$ ,弯矩 $M_{r}^{b}$ 和扭矩 $M_{t}^{b}$ 分别为^[16]

$F_{n}^{b} = {(F_{n}^{b})}_{o} - k_{n}^{b} Δ δ_{n}^{}$ ,

(11)

$F_{s}^{b} = {(F_{s}^{b})}_{o} - k_{s}^{b} Δ δ_{s}^{}$ ,

(12)

$M_{r}^{b} = {(M_{r}^{b})}_{o} - k_{r}^{b} Δ θ_{r}^{}$ ,

(13)

$M_{t}^{b} = {(M_{t}^{b})}_{o} - k_{t}^{b} Δ θ_{t}$ 。

(14)

式中 $k_{n}^{b}$ , $k_{s}^{b}$ , $k_{r}^{b}$ , $k_{t}^{b}$ 分别为胶结法向、切向、抗弯和抗扭刚度, $k_{t}^{b} = 0.5 k_{s}^{b} R_{b}^{2}$ , $k_{r}^{b} = 0.25 k_{n}^{b} R_{b}^{2}$ ； $Δ δ_{n}^{}$ 是法向位移增量； $R_{b} = λ_{b} R$ 为胶结半径, $λ_{b}$ 为胶结半径系数。

为了考虑胶结厚度和宽度对胶结刚度的影响,用下式计算胶结法向刚度（图1）^[15]：

图 1 胶结刚度计算示意图

Figure 1. Schematic diagram for calculation of bond stiffness

下载: 全尺寸图片幻灯片

$k_{n}^{b} = 1 / (δ_{A 1} + δ_{A 2} + δ_{B 1} + δ_{B 2} + δ_{C 1} + δ_{C 2})$ ,

(15)

式中, $δ_{A 1}$ , $δ_{B 1}$ , $δ_{C 1}$ 分别为颗粒1的A1_,B1和C1三段的线柔度,

$δ_{A 1} = \frac{R_{1} ξ_{E}}{π E_{1} R_{b}^{2}}$ ,

(16)

$δ_{B 1} = {\begin{matrix} \frac{1}{2 π (E_{1} - E_{2}) R_{1} Λ_{E}} (\ln (\frac{Λ_{E} + 1}{Λ_{E} - 1}) - \ln (\frac{Λ_{E} + ξ_{E}}{Λ_{E} - ξ_{E}})) (E_{1} > E_{2}) \\ \frac{R_{1} - R_{1} ξ_{E}}{π E_{1} R_{b}^{2}} (E_{1} = E_{2}) \end{matrix}$ ,

(17)

$δ_{C 1} = δ_{C 2} = \frac{h_{0}^{b} / 2}{π E_{2} R_{b}^{2}}$ ,

(18)

式中, $E_{1}$ , $E_{2}$ 分别为骨架材料和胶结材料的等效变形模量, $R_{1}$ 为颗粒1的半径, $h_{0}^{b}$ 为胶结厚度,同理可计算 $δ_{A 2}$ , $δ_{B 2}$ 。在模型使用中, $E_{1} = {\bar{E}}^{*}$ 为胶结等效模量, $E_{2} = η_{E} {\bar{E}}^{*}_{}$ , $η_{E}$ 为模量折减系数。 $ξ_{E}$ , $Λ_{E}$ 计算如下：

$ξ_{E} = \sqrt{1 - {(R_{b} / R_{1})}^{2}}$ ,

(19)

$Λ_{E} = \sqrt{1 + \frac{E_{2}}{E_{1} - E_{2}} {(R_{b} / R_{1})}^{2}}$ 。

(20)

胶结切向刚度根据法切向刚度比 $κ_{b}$ 计算：

$k_{s}^{b} = k_{n}^{b} / κ_{b}$ 。

(21)

采用考虑胶结厚度和宽度的胶结强度计算方法^[15,17]。胶结抗压和抗拉强度计算如下：

$R_{nc}^{b} = π R_{b}^{2} σ_{c}^{b} \frac{c_{1}}{λ_{b}^{c_{3}}} \exp (\frac{c_{2} h_{0}^{b}}{2 R_{b}})$ ,

(22)

$R_{nt}^{b} = π R_{b}^{2} σ_{t}^{b}$ ,

(23)

式中, $σ_{c}^{b}$ , $σ_{t}^{b}$ 分别为胶结材料抗压和抗拉强度； $c_{1} = 1.45$ , $c_{2} = - 0.7$ 和 $c_{3} = 1 / 6$ 是反映胶结半径和厚度对抗压强度影响的参数（定值）。

在一定的法向力下（未达到胶结抗压和抗拉强度）,胶结可能承受剪切力、弯矩和扭矩,复合荷载作用下胶结强度符合椭球包面：

${(‖ F_{s}^{b} ‖ / R_{s}^{b})}^{2} + {(‖ M_{r}^{b} ‖ / R_{r}^{b})}^{2} + {(M_{t}^{b} / R_{t}^{b})}^{2} = 1$ ,

(24)

式中, $R_{s}^{b}$ , $R_{r}^{b}$ , $R_{t}^{b}$ 分别为胶结抗剪、抗弯和抗扭强度,

$R_{s}^{b} = S_{s}^{b} (R_{nc}^{b} + R_{nt}^{b})$ ,

(25)

$R_{r}^{b} = S_{r}^{b} (R_{nc}^{b} + R_{nt}^{b}) R_{b}$ ,

(26)

$R_{t}^{b} = S_{t}^{b} (R_{nc}^{b} + R_{nt}^{b}) R_{b}$ ,

(27)

式中, $S_{s}^{b}$ , $S_{r}^{b}$ , $S_{t}^{b}$ 为胶结抗剪、抗弯和抗扭包面的形状,

$S_{i}^{b} = m_{i} f_{n}^{b} {[\ln (1 / f_{n}^{b})]}^{3 / 5}$ ,

(28)

式中,下标i是s,r和t, $m_{s} = 0.5$ , $m_{r} = 0.3$ 和 $m_{t} = 0.36$ 为参数（定值）, $f_{n}^{b} = (F_{n}^{b} + R_{nt}^{b}) / (R_{nc}^{b} + R_{nt}^{b})$ 为正则化胶结法向力。

胶结的存在对结构性土体的宏微观力学性质有很大影响。胶结达到强度包面后会发生胶结破坏,胶结破坏后胶结部分不再承担力和力矩,接触力和力矩退化到残余值（颗粒接触部分）。

2. 离散元模拟步骤

2.1 离散元试样制备

以饱和结构性黄土^[15]为基准制备离散元试样和选取接触参数,为提高研究的适用性开展了不同胶结材料强度离散元试样的数值分析,下文仅以一种胶结材料强度为例进行分析。采用分层欠压法^[18]分5层制备离散元试样,颗粒数目42180个,如图2（a）所示。结构性土离散元试样制备采用的颗粒级配曲线如图2（b）所示。试样生成后在12.5 kPa竖向压力下压缩稳定以反映原位K₀应力状态。试样稳定后若两个颗粒之间的间距小于 $g_{c} R_{b}$ 则生成胶结, $g_{c}$ 为胶结临界厚度系数。

图 2 离散元试样及其粒径级配

Figure 2. DEM sample and grain-size distribution curve

下载: 全尺寸图片幻灯片

颗粒间引力作用能使试样以更松散状态稳定赋存^[19],试样制备过程施加弱吸引力以制备松散试样^[15]：

$F_{a} = σ_{a} d_{50}^{2}$ ,

(29)

式中,d₅₀为中值粒径, $σ_{a}$ 为弱吸引力系数,根据典型黄土的最大稳定孔隙比和饱和重塑土黏聚力,通过试算取4 kPa。颗粒与墙体间仅承担法向力作用,忽略摩擦。颗粒密度=2710 kg/m³,颗粒局部阻尼系数=0.7。颗粒之间的接触参数：颗粒等效模量=800 MPa,法切向刚度比=1.5,摩擦系数=0.5,抗转动系数=0.21,胶结等效模量=200 MPa,胶结模量折减系数=0.2,胶结法切向刚度比=2,胶结半径系数=0.35,胶结材料抗压强度=25MPa,胶结材料抗拉强度=2.5MPa,胶结临界厚度系数=0.1。

2.2 不同应力路径离散元模拟步骤

为了研究不同应力路径结构性土的胶结破损演化,对上述结构性土离散元试样进行了分级加载的侧限压缩试验、等向压缩试验和等应力比压缩试验,以及应变加载的常规三轴试验和等p真三轴试验。

对分级加载压缩试验,上下墙体从12.5 kPa开始分级施加荷载,每次施加的荷载约为上级荷载的 $\sqrt{2}$ 倍；侧墙伺服应力目标通过应力比要求计算（侧限压缩试验使4个侧墙固定）。

对三轴试验,首先将离散元试样等向压缩到不同的围压,然后对侧墙进行伺服控制,使上下两墙相向运动进行剪切试验。在常规三轴试验中,侧墙的伺服目标为围压。在真三轴试验中通过平均应力和中主应力系数b不变的条件计算侧墙伺服目标应力。

3. 离散元模拟结果

为验证离散元数值分析结果的有效性,在研究结构性土的胶结破损演化规律之前,需确保离散元试样能够反映结构性土的主要压缩和剪切特性。考虑到本文针对一般结构性土体（不同胶结强度的胶结颗粒材料）进行破损规律数学建模,因此,下面主要进行定性对比,未对某一特定的结构性土进行标定。

图3给出了结构性土离散元试样在不同压缩路径下的压缩曲线。图3中 $η = q / p$ 为应力比, $CSR = σ_{3} / σ_{1}$ 为主应力比。随着平均应力的增加,离散元试样压缩线先缓慢下降而后逐渐转变为快速下降,存在结构屈服应力,离散元试样能够定性再现结构性土体的主要压缩特性^[20]。离散元试样体积变形特征与真实结构性土定量上差别较大,主要是因为颗粒形状、颗粒破碎机制等微观结构的差别。

图 3 结构性土压缩性质

Figure 3. Compression properties of structured soils

下载: 全尺寸图片幻灯片

图4给出了结构性土离散元试样常规三轴试验应力–应变关系和孔隙比–应变关系。离散元分析能够定性再现结构性土体三轴剪切性质^[2,21]。试样在剪切初始段模量很大,近似弹性段。随着围压的增加,试样抗剪强度增加,孔隙比减小,试样从倾向于软化和剪胀向倾向于硬化和剪缩发展。离散元反映结构性土剪切破坏特性的效果优于反映压缩变形特性。

图 4 结构性土常规三轴剪切性质

Figure 4. Conventional triaxial test results of structured soils

下载: 全尺寸图片幻灯片

图5给出了结构性土离散元试样真三轴试验应力-应变关系和孔隙比–应变关系（以p=300 kPa为例）。可见离散元模拟结果能够反映室内试验中中主应力系数b（ $b = (σ_{2} - σ_{3}) / (σ_{1} - σ_{3})$ ）的影响^[22-23]。同样平均应力下随着b的增加,破坏偏应力下降。b=0时试样倾向于剪缩,随着b的增加试样剪缩程度减弱。天然土试样的变形模量相对离散元试样和人工制备结构性土试样较低,这可能与试样扰动有关。

图 5 结构性土真三轴剪切性质

Figure 5. True triaxial test results of structured soils

下载: 全尺寸图片幻灯片

4. 胶结破损演化规律

4.1 胶结颗粒材料微观力学

土体可以看成一种颗粒集合体,重塑土（无黏结）可以看成是一种散粒体材料而结构性土可以看成是一种胶结颗粒材料。对于理想的散粒体材料,土中的应力通过接触力进行传递,而土中的应变由颗粒间的相对变形产生。基于此,学者通过研究土颗粒之间的接触关系和微观变形机理,建立了理想散粒体材料宏观应力和应变与微观接触力和变形之间的基本方程^[24-25]。Jiang等^[10]在此基础上建立了胶结颗粒材料的宏微观理论框架。

对于颗粒材料,荷载通过颗粒之间的接触力进行传递,因此颗粒材料的代表性单元的平均应力张量可根据下式求得

${\bar{σ}}_{i j}^{} = \frac{1}{V^{}} \sum_{k = 1}^{N^{}} F_{i}^{k} l_{j}^{k}$ ,

(30)

式中,V为代表性单元的体积,N为代表性单元内颗粒间接触点的总数, $F_{i}^{k}$ 为第k个接触点处接触力的分量, $l_{j}^{k}$ 为第k个接触点处方向向量的分量。

对于结构性土体,根据颗粒间接触点是否存在胶结,土体代表性单元平均应力张量 ${\bar{σ}}_{i j}^{}$ 可分解为胶结接触分担应力 ${\bar{σ}}_{i j}^{in}$ 与无胶结接触分担应力 ${\bar{σ}}_{i j}^{d}$ 之和,

$\begin{array}{l} {\bar{σ}}_{i j} = \frac{1}{V^{}} (\sum_{k = 1}^{N^{in}} F_{i}^{k} l_{j}^{k} + \sum_{k = 1}^{N^{d}} F_{i}^{k} l_{j}^{k}) \\ = \frac{1}{V^{}} （ V^{in} {\bar{σ}}_{i j}^{in} + V^{d} {\bar{σ}}_{i j}^{d}) = (1 - λ) {\bar{σ}}_{i j}^{in} + λ {\bar{σ}}_{i j}^{d} \end{array}$ ,

(31)

式中, $N^{in}$ , $N^{d}$ 分别为代表性单元内胶结接触点和无胶结接触点的数目,且 $N = N^{in} + N^{d}$ 。 $V^{in}$ , $V^{d}$ 分别为代表性单元内胶结颗粒所占体积与无胶结颗粒所占体积（含相应的孔隙体积）,且 $V = V^{in} + V^{d}$ ； $λ = V^{d} / V$ 为体积破损率,即代表性单元内无胶结组分所占体积与代表性单元总体积的比值； ${\bar{σ}}_{i j}^{in}$ 和 ${\bar{σ}}_{i j}^{d}$ 可由下式计算：

$\begin{array}{l} {\bar{σ}}_{i j}^{in} = \frac{1}{V^{in}} \sum_{k = 1}^{N^{in}} F_{i}^{k} l_{j}^{k}, \\ {\bar{σ}}_{i j}^{d} = \frac{1}{V^{d}} \sum_{k = 1}^{N^{d}} F_{i}^{k} l_{j}^{k} 。 \end{array}}$

(32)

为了研究结构性土体的结构性演化,需要在本构模型中引入（或蕴含）胶结破损参量。体积破损率就是一个简单直观的胶结破损参量,考虑到胶结组分和无胶结组分的体积难以划分,可用无胶结接触数量与总接触数量比值近似表示,

$λ = \frac{V^{d}}{V} \approx \frac{N^{d}}{N}$ 。

(33)

无胶结组分分担的应力与总应力的比值也可以作为胶结破损参量,

$B_{0} = \frac{λ {\bar{σ}}_{i j}^{d}}{{\bar{σ}}_{i j}} = \frac{\sum_{k = 1}^{N^{d}} F_{i}^{k} l_{j}^{k}}{\sum_{k = 1}^{N} F_{i}^{k} l_{j}^{k}}$ 。

(34)

胶结破损参量表征胶结破损程度或结构性衰退程度。当胶结破损参量为0时,土体结构性完好,土体应力全部由胶结组分承担；随着胶结破损参量的增加,土体应力由胶结组分和无胶结组分共同分担（胶结组分承担应力逐渐变小）；当胶结破损参量为1时,土体结构性完全破坏,土体应力全部由无胶结组分承担。

胶结破损参量 $λ$ , $B_{0}$ 难以通过室内试验获得,离散单元法可以用来研究 $λ$ , $B_{0}$ 的演化。尽管 $λ$ , $B_{0}$ 参量的大小有一定不同,但二者演化规律是相似的,因此本文以破损参量 $B_{0}$ 为例进行研究。

4.2 胶结破损演化规律

随着塑性应变的发展,胶结逐渐破损。由于弹性应变较小,塑性应变采用如下简单方法计算。离散元试样的弹性能 $E_{k}$ 为颗粒弹性能 $E_{k}^{p}$ 和胶结弹性能 $E_{k}^{b}$ 之和：

$E_{k}^{p} = \frac{1}{2} (\frac{{(F_{n}^{l})}^{2}}{k_{n}} + \frac{{‖ F_{s}^{l} ‖}^{2}}{k_{s}} + \frac{{‖ M_{r}^{l} ‖}^{2}}{k_{r}} + \frac{{(M_{t}^{l})}^{2}}{k_{t}})$ ,

(35)

$E_{k}^{b} = \frac{1}{2} (\frac{{(F_{n}^{b})}^{2}}{k_{n}^{b}} + \frac{{‖ F_{s}^{b} ‖}^{2}}{k_{s}^{b}} + \frac{{‖ M_{r}^{b} ‖}^{2}}{k_{r}^{b}} + \frac{{(M_{t}^{b})}^{2}}{k_{t}^{b}})$ ,

(36)

$E_{k} = E_{k}^{p} + E_{k}^{b} = \frac{1}{2} (p ε_{v}^{e} + q ε_{s}^{e}) V_{DEM}$ ,

(37)

$K = p / ε_{v}^{e}$ ,

(38)

$G = q / ε_{s}^{e}$ ,

(39)

式中, $V_{DEM}$ 为离散元试样的体积, $ε_{v}^{e}$ 为弹性体应变, $ε_{s}^{e}$ 为弹性偏应变,假设体积模量K和剪切模量G相等,联立以上各式可求得 $ε_{v}^{e}$ , $ε_{s}^{e}$ ,从而求出塑性体应变和塑性偏应变。

图6给出了破损参量随塑性体应变和塑性偏应变的变化规律。对于三轴剪切试验,塑性体应变较小且变化复杂,而对于等应力比压缩试验（应力比较小时）,塑性偏应变较小,二者均不适合单独用于描述试样破损参量演化。因此,胶结的破损既与塑性体应变有关,又与塑性偏应变有关。

图 6 不同应力路径下破损参量B₀随塑性体应变和塑性偏应变的演化规律

Figure 6. Evolution of degradation variable B₀ versus plastic volumetric and deviator strain for DEM sample under different stress paths

下载: 全尺寸图片幻灯片

为了考虑两种塑性应变对胶结破损参量的影响,采用下式简单等效塑性应变研究胶结破损演化。

$ε_{d}^{p} = \sqrt{{(ε_{v}^{p})}^{2} + {(ε_{s}^{p})}^{2}}$ ,

(40)

式中, $ε_{v}^{p}$ 为塑性体应变, $ε_{s}^{p}$ 为塑性偏应变。

图7给出结构性土离散元试样胶结破损参量 $B_{0}$ 随等效塑性应变的变化规律。从图7中可见,对不同应力路径,破损参量 $B_{0}$ 随等效塑性应变的增加先较快增加,而后缓慢增加。但不同应力路径的破损参量 $B_{0}$ 演化在数值上有较大不同。对不同应力比压缩试验,随着应力比的增加, $B_{0}$ 表征的胶结破损速率加快；对不同中主应力系数b的真三轴试验,随着b的增加, $B_{0}$ 表征的胶结破损速率减小；对不同围压常规三轴试验,围压也对 $B_{0}$ 表征的胶结破损速率影响显著。

图 7 不同应力路径下破损参量B₀随等效塑性应变的演化规律

Figure 7. Evolution of degradation variable B₀ versus equivalent plastic strain for DEM sample under different stress paths

下载: 全尺寸图片幻灯片

考虑到胶结破损跟胶结强度相关,在宏观上,胶结破损速率跟目前应力状态和强度包线之间的相对关系相关,结构参数 $B_{0}$ 破损的速率随应力比的增加而增加。因此,将应力比和强度包线引入胶结破损参量 $B_{0}$ ,通过应力比与峰值应力比的比值 $η / M_{f}$ 调整胶结破损速率,构建一个新的胶结破损参量如下式,从而蕴含应力路径对胶结破损的影响,

$B_{σ} = \frac{2 B_{0}}{(1 + η^{2} / M_{f}^{2})}$ ,

(41)

式中, $M_{f} = q_{f} / p$ 为峰值应力比,本文通过常规三轴试验结果采用莫尔–库仑准则得到峰值偏应力 $q_{f} =$ $0.9275 p + 123.8 kPa$ 。

图8给出了结构性土离散元试样胶结破损参量 $B_{σ}$ 随等效塑性应变的破损规律。对不同应力路径,胶结破损参量 $B_{σ}$ 随等效塑性应变逐渐增加,但破损速率逐渐降低。除等向压缩试验（ $η = 0$ ）外,试验应力比、围压和中主应力系数对 $B_{σ}$ 表征的胶结破损演化影响很小。由此可见,胶结破损参量 $B_{σ}$ 应力路径相关性明显降低,可作为状态变量引入结构性土体二元介质类本构模型（比如扰动状态概念模型、损伤模型等）表征胶结破损情况,使模型的结构破损规律更为合理。通过数据分析,建议采用指数函数拟合破损参量 $B_{σ}$ 与等效塑性应变的关系,

图 8 不同应力路径下破损参量B_σ随等效塑性应变的演化规律

Figure 8. Evolution of degradation variable B_σversus equivalent plastic strain for DEM sample under different stress paths

下载: 全尺寸图片幻灯片

$B_{σ} = 1 - \exp [- c_{a} {(ε_{d}^{p})}^{c_{b}}]$ ,

(42)

式中, $c_{a}$ , $c_{b}$ 为胶结破损相关的参数,可通过侧限压缩或常规三轴试验确定。

5. 结论

通过将相对完备的三维胶结接触模型植入粒间接触,反映结构性土粒间胶结的关键特征,建立了结构性土离散元模型；开展了侧限压缩、等向压缩、等应力比压缩以及常规三轴、真三轴等不同加载路径下的离散元数值分析；再现了结构性土体的主要宏观力学性质（尤其是结构性土的破坏特性）；分析了胶结破损参量 $B_{0}$ , $B_{σ}$ 随等效塑性应变的演化规律。

（1）胶结破损参量（ $λ$ , $B_{0}$ , $B_{σ}$ ）表征胶结破损程度或结构性衰退程度。当胶结破损参量为0时,土体结构性完好,土体应力全部由胶结组分承担；当胶结破损参量为1时,土体结构性完全破坏,土体应力全部由无胶结组分承担。

（2）胶结破损参量 $B_{0}$ 的演化受应力比、围压和中主应力系数等的影响,具有应力路径相关性。将应力比和强度包线引入胶结破损参量 $B_{0}$ ,构建了一个新的胶结破损参量 $B_{σ}$ ,应力路径相关性明显降低。在本构模型应用中,胶结破损参量 $B_{σ}$ 随等效塑性应变的演化规律可用指数函数表示,并通过侧限压缩或常规三轴试验确定相关参数。

不同应力路径下 $B_{σ}$ 的演化仍然具有一定的离散性,今后将考虑改进等效塑性应变的定义和强度包面公式。因为胶结破损规律 $B_{σ}$ 基于离散元结果提出,适用于理想胶结颗粒材料（理想结构性土体）,对真实结构性土体的适用性需要进一步验证。

图 1 不同初始干密度下的SWCC

Figure 1. SWCCs under different initial dry densities

下载: 全尺寸图片幻灯片

图 2 不同初始干密度下南阳膨胀土在全吸力范围的SWCC

Figure 2. SWCCs of Nanyang expansive soil in full suction range under different initial dry densities

下载: 全尺寸图片幻灯片

图 3 不同因素影响下初始吸力与最大膨胀率的关系

Figure 3. Relationship between initial suction and maximum swelling rate under influence of different factors

下载: 全尺寸图片幻灯片

图 4 不同影响因素下膨胀力的变化规律

Figure 4. Dependences of swelling force on different influencing factors

下载: 全尺寸图片幻灯片

图 5 膨胀力与最终吸水量的关系

Figure 5. Relationship between swelling force and final water absorption

下载: 全尺寸图片幻灯片

图 6 式（8a）的可靠性

Figure 6. Reliability of Formula (8a)

下载: 全尺寸图片幻灯片

表 1 试验方案

Table 1 Testing programs

类别	ρ_d0/(g·cm^-3)	w₀ /%
无荷膨胀率	1.70	12,14,16,18,20,22
	1.60	12,22
	1.50	12,18
	1.40	18,22
无荷膨胀力	1.40,1.45,1.50,1.55,1.60	12,14,16,18,20,22
SWCC	1.40,1.45,1.50,1.55,1.60	20

下载: 导出CSV

表 2 南阳膨胀土SWCC的相关参数

Table 2 Related parameters of SWCCs of Nanyang expansive soil

$ρ_{d0}$ /(g·cm^-3)	VG模型				Webb迭代法
$ρ_{d0}$ /(g·cm^-3)	α /kPa	n	w_r	R²	s^*	w^*
1.40	6.29	1.095	0	0.998	58.16	0.277
1.45	14.19	1.097	0	0.998	128.83	0.259
1.50	37.33	1.103	0	0.997	319.03	0.240
1.55	101.94	1.108	0	0.997	836.04	0.217
1.60	159.78	1.119	0	0.991	1184.21	0.206

下载: 导出CSV

表 3 膨胀力与初始干密度和最终吸水量关系的拟合参数

Table 3 Fitting parameters of relationship among swelling force, initial dry density and final water absorption

w ₀ /%	式（5a）			式（6）			式（7）
w ₀ /%	k₁ /kPa	λ₁	R²	k₂ /kPa	λ₂	R²	a₂/kPa	b₂	R²
11.11	387.482	9.519	0.880	0.379	-3.040	0.947	0.395	3.022	0.989
13.57	374.323	8.302	0.968	0.535	-2.720	0.979	0.692	2.572	0.985
16.16	564.711	10.515	0.986	0.392	-2.543	0.975	0.345	2.616	0.993
17.23	851.601	11.865	0.986	0.160	-2.854	0.957	0.062	3.317	0.966
18.99	490.285	9.326	0.970	0.759	-1.945	0.985	1.923	1.552	0.979

下载: 导出CSV

表 4 初始干密度–初始吸力–膨胀力关系的拟合参数

Table 4 Fitting parameters of relationship among initial dry density, initial suction and swelling force

w ₀ /%	式（9a）
w ₀ /%	α₂	β₂	γ₂	R²
11.11	157.165	89.231	5.229	0.965
13.57	338.373	103.752	7.821	0.970
16.16	564.711	353.000	10.515	0.986
17.23	0.162	9.117	9.000	0.981
18.99	417.939	87.431	8.562	0.973

下载: 导出CSV

表 5 初始含水率–初始吸力–膨胀力关系的拟合参数

Table 5 Fitting parameters of relationship among initial water content, initial suction and swelling force

$ρ_{d 0}$ /(g·cm^-3)	式（10）
$ρ_{d 0}$ /(g·cm^-3)	α₃	β₃	R²
1.40	54.478	0.873	0.992
1.45	58.683	0.719	0.999
1.50	167.698	1.114	0.985
1.55	319.560	1.226	0.994
1.60	299.223	1.090	0.986

下载: 导出CSV

参考文献(20)

[1]	陈正汉. 非饱和土与特殊土力学的基本理论研究[J]. 岩土工程学报, 2014, 36(2): 201-272. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201402002.htm CHEN Zheng-han. On basic theories of unsaturated soils and special soils[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(2): 201-272. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201402002.htm
[2]	李献民, 王永和, 杨果林, 等. 击实膨胀土工程变形特征的试验研究[J]. 岩土力学, 2003, 24(5): 826-830. doi: 10.3969/j.issn.1000-7598.2003.05.029 LI Xian-min, WANG Yong-he, YANG Guo-lin, et al. Test study on engineering deformation characteristics of compacted expansive soil[J]. Rock and Soil Mechanics, 2003, 24(5): 826-830. (in Chinese) doi: 10.3969/j.issn.1000-7598.2003.05.029
[3]	徐永福, 孙德安, 董平. 膨润土及其与砂混合物的膨胀试验[J]. 岩石力学与工程学报, 2003, 22(3): 451-455. doi: 10.3321/j.issn:1000-6915.2003.03.022 XU Yong-fu, SUN De-an, DONG Ping. Swelling tests on bentonite and sand-bentonite mixture[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(3): 451-455. (in Chinese) doi: 10.3321/j.issn:1000-6915.2003.03.022
[4]	于海浩, 孙德安, 韦昌富, 等. 氯化钠溶液饱和不同初始含水率膨润土的膨胀特性[J]. 岩土工程学报, 2019, 41(3): 595-600. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201903028.htm YU Hai-bo, SUN De-an, WEI Chang-fu, et al. Swelling characteristics of bentonite with different initial water contents saturated by NaCl solution[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(3): 595-600. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201903028.htm
[5]	石北啸, 陈生水, 韩华强, 等. 考虑吸力变化的膨胀土边坡破坏规律分析[J]. 水利学报, 2014, 45(12): 1499-1506. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201412016.htm SHI Bei-xiao, CHEN Sheng-shui, HAN Hua-qiang, et al. Failure law of the expansive soil slope considering the suction changes[J]. Journal of Hydraulic Engineering, 2014, 45(12): 1499-1506. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201412016.htm
[6]	韩华强, 陈生水, 郑澄锋. 非饱和膨胀土强度及变形特性试验研究[J]. 岩土工程学报, 2008, 30(12): 1872-1876. doi: 10.3321/j.issn:1000-4548.2008.12.017 HAN Hua-qiang, CHEN Sheng-shui, ZHENG Cheng-feng. Experimental study on strength and deformation of unsaturated expansive soils[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(12): 1872-1876. (in Chinese) doi: 10.3321/j.issn:1000-4548.2008.12.017
[7]	贾景超. 膨胀土膨胀机理及细观膨胀模型研究[D]. 大连: 大连理工大学, 2010. JIA Jing-chao. Study on the Swelling Mechanism and Mesomechanical Swelling Model of Expansive Soils[D]. Dalian: Dalian University of Technology, 2010. (in Chinese)
[8]	范秋雁, 梁昕, 韩进仕. 非饱和膨胀岩饱和度及胀缩特性试验研究[J]. 岩石力学与工程学报, 2020, 39(1): 45-56. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202001005.htm FAN Qiu-yan, LIANG Xin, HAN Jin-shi. Experimental study on saturation and swelling-shrinkage characteristics of unsaturated expansive rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(1): 45-56. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202001005.htm
[9]	邹维列, 叶云雪, 韩仲. 基于化学反应动力学方程的一维膨胀时程模型[J]. 岩土工程学报, 2020, 42(4): 737-744. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202004022.htm ZOU Wei-lie, YE Yun-xue, HAN Zhong. Chemical reaction kinetics based models for describing evolution of one-dimensional expansion with time[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 737-744. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202004022.htm
[10]	WEBB S W. A simple extension of two-phase characteristic curves to include the dry region[J]. Water Resources Research, 2000, 36(6): 1425-1430.
[11]	LU S, REN T, GONG Y, et al. Evaluation of three models that describe soil water retention curves from saturation to oven dryness[J]. Soil Science Society of America Journal, 2008, 72(6): 1542-1546.
[12]	SCHNEIDER M, GOSS K U. Prediction of water retention curves for dry soils from an established pedotransfer function: evaluation of the Webb model[J]. Water Resources Research, 2012, 48(6): W06603.
[13]	叶云雪, 邹维列, 袁斐, 等. 基于土壤转换函数(PTF)预测不同初始孔隙比土的土-水特征曲线[J]. 岩土工程学报, 2018, 40(12): 2305-2311. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201812022.htm YE Yun-xue, ZOU Wei-lie, YUAN Fei, et al. Predicating soil-water characteristic curves of soils with different initial void ratios based on a pedotransfer function[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(12): 2305-2311. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201812022.htm
[14]	刘岩, 吴天前, 邹维列. 强膨胀土自然膨胀特性的试验研究[J]. 华中科技大学学报：自然科学版, 2020, 48(2): 29-35. https://www.cnki.com.cn/Article/CJFDTOTAL-HZLG202002006.htm LIU Yan, WU Tian-qian, ZOU Wei-lie. Experimental study on humidification deformation of expansive soil[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2020, 48(2): 29-35. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HZLG202002006.htm
[15]	BURTON G J, SHENG D, CAMPBELL C. Bimodal pore size distribution of a high-plasticity compacted clay[J]. Géotechnique Letters, 2014, 4(2): 88-93.
[16]	谭罗荣, 孔令伟. 膨胀土膨胀特性的变化规律研究[J]. 岩土力学, 2004, 25(10): 1555-1559. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200410009.htm TAN Luo-rong, KONG Ling-wei. Study on variation regularity of swelling behavior of expansive soil[J]. Rock and Soil Mechanics, 2004, 25(10): 1555-1559. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200410009.htm
[17]	MANCA D, FERRARI A, LALOUI L. Fabric evolution and the related swelling behaviour of a sand/bentonite mixture upon hydro-chemo-mechanical loadings[J]. Géotechnique, 2016, 66(1): 41-57.
[18]	ZOU W L, YE J B, HAN Z, et al. Effect of montmorillonite content and sodium chloride solution on the residual swelling pressure of an expansive clay[J]. Environmental Earth Sciences, 2018, 77(19): 1-12.
[19]	卢肇钧, 吴肖茗, 孙玉珍, 等. 膨胀力在非饱和土强度理论中的作用[J]. 岩土工程学报, 1997, 19(5): 20-27. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC705.003.htm LU Zhao-jun, WU Xiao-ming, SUN Yu-zhen, et al. The role of swelling pressure in the shear strength theory of unsaturated soils Chinese Journal of Geotechnical Engineering, 1997, 19(5): 20-27. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC705.003.htm
[20]	TU H, VANAPALLI S K. Prediction of the variation of swelling pressure and 1-D heave of expansive soils with respect to suction using the soil water retention curve as a tool[J]. Canadian Geotechnical Journal, 2016, 53(8): 1213-1234.

施引文献(29)

期刊类型引用(19)

1.	沈书拯，江巍，欧阳晔，谭亿虹，吴剑，陈勇. 基于三维连续-离散耦合模型的门架式双排抗滑桩受力分析. 防灾减灾工程学报. 2025(01): 84-94 . 百度学术
2.	王洪涛，李建华，殷允腾，张红军，宣兆腾，赵明珠. 承压水作用下岩质地层基坑抗突涌安全厚度研究. 建筑结构学报. 2024(06): 213-223 . 百度学术
3.	马勤. 超挖对软土基坑工程稳定性的影响分析. 现代矿业. 2024(04): 121-123 . 百度学术
4.	张瑞程. 新型斜撑-排桩基坑支护体系工程应用分析. 施工技术(中英文). 2024(12): 162-166 . 百度学术
5.	王强，宋歌，廖利，刘国伟，蔡崇德. 全地下空间枢纽深基坑支护结构变形特征控制. 铁路技术创新. 2024(03): 178-185 . 百度学术
6.	赵升峰，汪敏营，邬喜春，李明东. 不同深度基坑中桩加支撑支护特性研究. 江苏建筑. 2024(06): 112-115+129 . 百度学术
7.	赵玉成，代兆宇，夏瑞萌，娄海成. 雄安新区深基坑双排桩支护的实践与探讨. 粉煤灰综合利用. 2024(06): 79-83+95 . 百度学术
8.	蔡红岩. 供水泵站基坑排桩支护局部失效引发连续破坏机理研究. 陕西水利. 2023(01): 100-102 . 百度学术
9.	曾诚. 基坑局部设置双排桩桩长桩径对连续变形控制效果影响研究. 中国新技术新产品. 2023(07): 81-85 . 百度学术
10.	梁文辉. 深基坑双排桩式围护结构支护施工技术研究. 江西建材. 2023(09): 244-245+248 . 百度学术
11.	黄礼明. 某邻近既有建筑的异形软土基坑变形控制. 福建建筑. 2023(11): 84-88 . 百度学术
12.	徐琳，黄程翔，邵根才，杨璞，李宝山. 软土地区基坑超挖对围护结构变形的影响. 建筑安全. 2023(12): 33-36 . 百度学术
13.	魏焕卫，李传斌，种记鑫，郑晓. 局部超挖对内支撑结构的施工效应. 计算机辅助工程. 2023(04): 56-65 . 百度学术
14.	常强. 双排桩治理边坡的数值模拟研究. 黑龙江水利科技. 2022(04): 37-39+179 . 百度学术
15.	郑刚，程雪松，周海祚，张天奇，于晓旋，刁钰，王若展，衣凡，张文彬，郭伟. 岩土与地下工程结构韧性评价与控制. 土木工程学报. 2022(07): 1-38 . 百度学术
16.	唐小军，袁仕贵. 铁路桥台基坑开挖支护方案及稳定性计算. 中国水运. 2022(08): 146-148 . 百度学术
17.	唐小军，袁仕贵. 铁路桥台基坑开挖支护方案及稳定性计算. 中国水运. 2022(15): 146-148 . 百度学术
18.	孙永梅，程海涛，米春荣，孙贝. 土体超挖对深基坑变形影响分析. 山西建筑. 2021(19): 61-63 . 百度学术
19.	曾锦秀. 基于上限法的边坡双排抗滑桩变形计算方法. 福建工程学院学报. 2021(06): 538-544 . 百度学术