动静组合加载下颗粒迁移与沉积特性试验研究

动静组合加载下颗粒迁移与沉积特性试验研究 English Version

Experimental study on migration and deposition of particles under alternating dynamic and static loads

0. 引言

1. 泡沫改良粗粒土不排水一维压缩有效应力计算模型

2. 泡沫改良粗粒土残余剪切强度计算

3. 模型验证

4. 讨论

5. 结论

期刊类型引用(3)

1.1 常压下孔隙比及饱和度求解

1.2 不排水一维压缩有效应力计算

2.1 不排水剪切引起的超孔隙压力计算

2.2 率相关残余剪切强度计算

3.1 试验方案

3.2 材料常数确定

3.3 结果验证

作者简介:
高峰(1992—)，男，博士后，主要从事冻土力学与工程、路基冻融翻浆病害机制及防控研究的工作。E-mail：gao-feng@csust.edu.cn

通讯作者:
张升，E-mail: zhang-sheng@csu.edu.cn

高峰; 张军辉; 张升; 郑健龙; 盛岱超

doi:10.11779/CJGE20230168

高峰^{1, 2,},
张军辉¹,
张升^3, ,,
郑健龙¹,
盛岱超⁴

1.
长沙理工大学交通运输工程学院公路养护技术国家工程研究中心, 湖南长沙 410114
2.
青海省交通控股集团有限公司科技信息部, 青海西宁 810000
3.
中南大学土木工程学院高速铁路建造技术国家工程研究中心, 湖南长沙 410075
4.
悉尼科技大学土木与环境工程学院, 悉尼 NSW 2007

基金项目:

国家自然科学基金杰出青年基金项目 52025085

国家重点研发计划项目 2021YFB2600900

国家博士后创新人才支持计划项目 BX20220259

湖南省公路先进建养技术国际科技创新合作基地开放基金重点项目 kfj220801

长沙市自然科学基金青年基金项目 kq2208234

详细信息

中图分类号: TU411
计量
- 文章访问数: 468
- HTML全文浏览量: 44
- PDF下载量: 135
出版历程
- 收稿日期: 2023-02-25
- 网络出版日期: 2023-06-01
- 刊出日期: 2024-04-30

1.
National Engineering Research Center of Highway Maintenance Technology, School of Traffic and Transportation Engineering, Changsha University of Science & Technology, Changsha 410114, China
2.
Department of Science and Information, Qinghai Communications Holding Group Co., Ltd., Xining 810000, China
3.
National Engineering Research Center for High-Speed Railway Construction, School of Civil Engineering, Central South University, Changsha 410075, China
4.
School of Civil & Environmental Engineering, University of Technology Sydney, Sydney NSW 2007, Australia

摘要

摘要: 悬浮颗粒在路基中迁移和沉积运动规律研究对揭示翻浆冒泥病害孕育、发展与致灾机制具有重要的理论意义和应用价值。通过开展动-静组合荷载下分层砾石-砂粉土柱翻浆试验，探究了动荷载加载强度、静荷载间歇时长、以及动静组合加载重复次数等对颗粒迁移和沉积行为的影响规律，并结合试样动水力响应特征分析了诱发颗粒迁移的驱动机制。试验结果表明：颗粒在动水力作用下发生悬浮迁移而在静荷载段受重力作用发生沉降，从而造成砾石层内泥浆浊度呈现波动型增长特征。随动静组合加载重复次数增加，每一个重复周期内的颗粒悬浮增量趋于降低，且悬浮颗粒逐渐堵塞砾石孔隙并在一定程度上抑制翻浆发展。提高动载频率和动应力会进一步降低砂粉土表层的内部稳定性，但相比于增大平均动应力，提高加载频率对增大颗粒迁移量和竖向迁移距离更加显著。
- 翻浆冒泥 /
- 室内试验 /
- 颗粒迁移 /
- 超孔隙水压
Abstract: The study on the migration and deposition of suspended particles in subgrade is important to reveal the generation and disaster-causing mechanism of mud pumping. The effects of magnitude of dynamic loads, intermittent duration of static loads and repetition number of alternating dynamic and static loads on the migration and deposition of particles are studied by carrying out the mud pumping tests on the layered gravel-sandy silt column. According to the hydrodynamic response characteristics of specimens, the driving mechanism of particle migration is analyzed. The test results show that the suspended particles migrate upward under hydrodynamic force and settle under gravity during the interval of static loads, which results in the fluctuating growth of slurry turbidity in the gravel layer. The particle suspension increment tends to decrease with the increasing repetition of alternating dynamic and static loads, and these particles gradually clog gravel pores and inhibit mud pumping to a certain extent. Increasing the dynamic frequency and stress can further reduce the internal stability of the sandy silt surface layer. However, compared with increasing the dynamic stress, increasing the loading frequency is more beneficial to increase the particle migration mass and vertical migration distance.
- mud pumping /
- laboratory test /
- particle migration /
- excess pore pressure

HTML全文

在土压平衡盾构掘进过程中，泡沫作为可压缩性介质能短期稳定地赋存于土舱内粗粒渣土的孔隙中，提高粗粒渣土的压缩性和塑流性，使其便于建立土舱压力，成为更有效支护开挖面的膏状介质^[1]。稳定开挖面所需的支护压力可通过调节改良参数以及渣土流入土舱和沿螺旋输送机排出的流量来实现，从而避免地层过大变形、周围构筑物损坏等风险^[2]。

Bezuijen^[3]提出合适的泡沫注入量能使泡沫改良土的孔隙比大于未改良土的最大孔隙比e_max。进一步地，为了优化改良参数和盾构掘进性能，坍落度试验广泛用于评估盾构渣土的塑流性，结果表明级配、含水率、泡沫注入量是影响泡沫改良渣土塑流性的关键因素^[4]。另有流变试验研究尝试用非牛顿流体的流变参数来表征改良粗粒土的流动性^[5]。然而，实际上盾构渣土是在压力舱内经螺旋输送机排出。由于泡沫具有高压缩性，泡沫改良粗粒土的剪切流动特性与环境压力密切相关。加压十字板剪切试验适用于测试强流动性改良土的带压剪切行为，且能测试较大剪切位移后的残余剪切强度，以反映盾构掘进时连续流动渣土的剪切行为^[6]。图 1展示了假定的圆柱形十字板剪切失效面，假设圆柱侧面剪切应力均匀分布，圆柱顶、底面剪切应力沿径向线性增加。十字板扭矩T的换算公式由积分求得：

。

$T = \pi {D^3}\left( {{\tau _\text{l}} + \frac{1}{6}{\tau _\text{v}}} \right)。$

(1)

图 1 十字板剪切失效面示意

Figure 1. Failure interface of vane shear tests

下载: 全尺寸图片幻灯片

式中： ${\tau _\text{v}}$ 为顶、底面的剪切应力； ${\tau _\text{l}}$ 为侧面的剪切应力。二者大小与所在剪切面的法向有效应力有关。

此外，十字板转速与剪切速率的关系式为^[7]

。

$\dot \gamma = \frac{{\pi n}}{{5400}} \times \frac{{{D^2}}}{{{D^2} - {d^2}}}。$

(2)

孟庆琳等^[7]发现改良渣土屈服应力和塑性黏度均随围压增加或泡沫注入比减小而增加；获得的压力相关性流变参数已应用于盾构土舱渣土输送数值模型^[8]。Hu等^[9]认为环境压力使泡沫收缩和改良土孔隙比减小，导致流变参数增大。

上述研究尚未考虑有效应力对剪切强度或流变参数的影响。随着环境压力增加，当改良土孔隙比小于最大值时，有效应力产生，而且有实例表明盾构实践中允许土舱内改良渣土存在有效应力^[3]，因此有效应力状态对改良土流动性的影响值得被关注。Mori等^[10]测试发现在有效应力出现前改良土十字板剪切强度很小，随着竖向总应力增加，当改良土孔隙比与最大值之比小于1.2，竖向有效应力开始发展，剪切强度才明显增加，二者的线性关系符合Mohr-Coulomb准则。Zhong等^[11]通过十字板试验探究了不同孔隙比泡沫改良砂的率相关剪切行为，存在有效应力的改良砂在十字板剪切过程中产生少量的超孔隙压力，进入残余状态后有效内摩擦角随剪切率增加而略有增加，残余状态改良砂的剪切行为表明它可被视为具有屈服应力的低黏度流体。

目前，对泡沫改良粗粒土的流动特性分析主要集中在试验研究，尚未有合适的理论模型来计算剪切强度或流变参数。Mori等^[10]基于理想气体方程推导了泡沫受压收缩后改良土的孔隙比公式，但不适用于描述含有效应力改良土的孔隙比；Wang等^[12]结合不排水一维压缩试验结果推导了不同竖向总应力下的孔隙压力理论解。由于既有结果反映了泡沫改良粗粒土剪切强度与有效应力的相关性^[10]，可在不排水压缩泡沫改良粗粒土应力状态计算的基础上进一步开发表征流动性的理论模型。

本次研究目标是填补一维压缩状态下泡沫改良粗粒土大变形流动后的残余剪切强度计算模型的空白，以便于预测不同环境压力和改良参数下大变形泡沫改良粗粒土的宏观流变行为。首先推导不排水一维压缩后的有效应力理论解，然后根据剪切引起的超孔隙压力计算方法和Mohr-Coulomb准则求解残余剪切强度，进而构建残余剪切强度与剪切速率的相关性公式，最后，开展了加压式十字板剪切试验来验证计算模型的准确性，并基于理论模型做了影响因素分析。

土内各组成成分的物理性质和体积分数影响着孔隙状态和有效应力特征，进而决定了宏观的力学行为。土舱内改良渣土在盾壳侧向约束和开挖面压力作用下处于一维压缩状态^[12]，而渣土流经土舱很快排出，可忽略渣土一维压缩下排水效应。此外，改良土的力学行为能在60 min内保持不变^[13]，因此盾构正常掘进时，可忽略渣土中泡沫衰变效应。

加载前泡沫改良粗粒土的孔隙状态可由常压下孔隙比和饱和度表征。图 2展示了粗粒土的各相组成。相比由土相、水相和气相组成的粗粒土（图 2（a）），泡沫改良粗粒土还多了泡沫这一特殊“相”（图 2（b）），它实际由泡沫剂溶液及其所包围的空气组成（图 2（c））。泡沫填充于土的孔隙，其添加量和填充状态会影响孔隙比和饱和度。

图 2 改良前后粗粒土多相示意图

Figure 2. Multiphase schematic diagram of coarse-grained soils with and without foam conditioning

下载: 全尺寸图片幻灯片

泡沫的添加量由泡沫注入比FIR表征，FIR定义为所注入的泡沫与粗粒土的体积比：

$\text{FIR} = \frac{{{V_\text{f}}}}{{{V_\text{ps}}}} \times 100\% 。$

(3)

式中：V_f为所注入泡沫的体积；V_ps为粗粒土的体积。

由可知，注入的泡沫主要起置换和扩大作用。置换作用表现为泡沫占据了粗粒土原孔隙空间，部分原有的气相和水相被排出粗粒土之外；而扩大作用体现在泡沫能阻隔土颗粒接触，一定程度扩大粗粒土的孔隙比。由此提出扩大系数 $\alpha$ 来量化泡沫用于扩大孔隙的体积比：

$\alpha = \frac{{{V_\text{sf}} - {V_\text{ps}}}}{{{V_\text{f}}}}。$

(4)

式中：V_sf为泡沫改良粗粒土总体积。扩大系数位于0~1。当 $\alpha$ 为0，泡沫注入前后粗粒土孔隙比不变；当 $\alpha$ 为1，泡沫全部贡献于扩大粗粒土的孔隙空间。此外，由图 2（c）可知，泡沫内空气与泡沫溶液的体积比可通过发泡倍率FER予以反映：

$\text{FER} = \frac{{{V_\text{f}}}}{{{V_\text{fl}}}} = \frac{{{V_\text{fl}} + {V_\text{fa}}}}{{{V_\text{fl}}}} 。$

(5)

式中：V_fl为泡沫液相体积；V_fa为泡沫气相体积。

因此，常压下泡沫改良粗粒土的孔隙比e_sf-ap可由稳定堆积状态下粗粒土的孔隙比e_ps，FIR和α计算得到，根据图 2的各相体积比例，推导过程为

$\begin{align} {e_\text{sf - ap}} &= \frac{{{V_\text{p - ap}}}}{{{V_\text{s}}}} = \frac{{{V_\text{w}} + {V_\text{a}} + \alpha {V_\text{f}}}}{{{V_\text{s}}}} \hfill \\ &= \frac{{{V_\text{w}} + {V_\text{a}}}}{{{V_\text{s}}}} + \frac{{\alpha \text{FIR}{V_\text{ps}}}}{{{V_\text{s}}}} = {e_\text{ps}} + \alpha \text{FIR}(1 + {e_\text{ps}}) \hfill \\ &= (1 + \alpha \text{FIR})(1 + {e_\text{ps}}) - 1 。 \end{align}$

(6)

式中：V_p-ap为泡沫改良粗粒土的孔隙总体积；V_s，V_w，V_a分别为粗粒土的固相、水相和气相体积。

另外，泡沫改良粗粒土的饱和度表征了被泡沫置换后粗粒土残留孔隙水和泡沫液相体积占孔隙体积的比值。据图 2各相比例可推导常压下的饱和度。

$\begin{align} {S_\text{r}} &= \frac{{{V_\text{wr}} + {V_\text{fl}}}}{{{V_\text{p - ap}}}} = \frac{{\left[ {{V_\text{w}} - \frac{{{V_\text{w}}}}{{{V_\text{a}} + {V_\text{w}}}}(1 - \alpha ){V_\text{f}}} \right] + \frac{{{V_\text{f}}}}{\text{FER}}}}{{{V_\text{a}} + {V_\text{w}} + \alpha {V_\text{f}}}}\, \\ &= \frac{{{V_\text{w}}\left[ {1 - \frac{{(1 - \alpha )\text{FIR}{V_\text{ps}}}}{{{V_\text{ps}} - {V_\text{s}}}}} \right] + \frac{\text{FIR}}{\text{FER}}{V_\text{ps}}}}{{{V_\text{ns}} - {V_\text{s}} + \alpha \text{FIR}{V_\text{ns}}}}\, \\ &= \frac{{\frac{{{V_\text{w}}}}{{{V_\text{s}}}}\left[ {1 - (1 - \alpha )(1 + {e_\text{ps}})\frac{\text{FIR}}{{{e_\text{ps}}}}} \right] + \frac{\text{FIR}}{\text{FER}}\frac{{{V_\text{ps}}}}{{{V_\text{s}}}}}}{{(1 + \alpha \text{FIR})\frac{{{V_\text{ps}}}}{{{V_\text{s}}}} - 1}} \\ &= \frac{{w\frac{{{\rho _\text{s}}}}{{{\rho _\text{w}}}}\left[ {1 - (1 - \alpha )(1 + {e_\text{ps}})\frac{\text{FIR}}{{{e_\text{ps}}}}} \right] + (1 + {e_\text{ps}})\frac{\text{FIR}}{\text{FER}}}}{{(1 + \alpha \text{FIR})(1 + {e_\text{ps}}) - 1}} 。 \end{align}$

(7)

式中：S_r为饱和度；V_wr为加入泡沫后留存在粗粒土孔隙中水的体积；ρ_s，ρ_w分别为土颗粒和孔隙水的密度；w为粗粒土的原始含水率。

因此，常压下泡沫改良粗粒土的饱和度求解需已知土颗粒和水的密度、粗粒土的孔隙比和原始含水率、常压下泡沫注入比、发泡倍率和扩大系数。

由于土颗粒和孔隙水常被假设为不可压缩介质，泡沫改良粗粒土的不排水压缩行为由受压收缩的封闭气相主导。假设不排水条件下粗粒土原有气相与泡沫内气相均为封闭气相，外加荷载后孔隙内封闭气相收缩导致土颗粒间距减小，当土颗粒接触并形成可传递应力的土骨架后，有效应力产生^[10]。因此，存在临界孔隙比e_th，将不排水一维压缩过程划分为无有效应力和存在有效应力两个阶段，e_th需通过试验获得。再者，泡沫改良粗粒土属于含气泡土，而有效应力原理对于含气泡土的适应性已被证实^[14]。

（1）e≥e_th加载阶段

当e≥e_th，土颗粒互不接触，竖向总应力增量全由孔隙压力承担（ ${\sigma _{\text{v}}}$ = u）。一般地，假设封闭气相服从Boyle定律和Henry定律^[15]，即密闭环境下恒温气体的压力与体积成反比关系，而与溶解度成正比关系。由于表面张力很小，封闭气相压力等于孔隙水压力加上大气压^[10]。孔隙压力增量将引起封闭气相收缩，泡沫改良粗粒土内封闭气相体积随孔隙压力的变化为

$\begin{array}{l} {p_\text{ap}}({V_\text{a - ap}} + h{V_\text{l}}) \\ \;\;\;\;\;\;\;= ({p_\text{ap}} + \Delta u)({V_\text{a - ap}} - \Delta {V_\text{a}} + h{V_\text{l}}) 。 \end{array}$

(8)

式中：V_a-ap为常压下泡沫改良粗粒土的气相体积；V_l为泡沫改良粗粒土的液相体积；p_ap为常压绝对值，取101.325 kPa；Δu为孔隙压力增量；ΔV_a为孔隙内气相体积增量；h为溶解度系数（即Henry常数），空气在水中的溶解度系数一般取0.02^[15]。

进而，常压下泡沫改良粗粒土的气相和液相体积可分别为

${V_\text{a - ap}} = {V_\text{p - ap}}(1 - {S_\text{r}}) = {e_\text{sf - ap}}{V_\text{s}}(1 - {S_\text{r}})\text{，}$

(9)

${V_\text{l}} = {e_\text{sf - ap}}{V_\text{s}}{S_\text{r}}。$

(10)

由于竖向总应力等于孔隙压力，将式（9），（10）代入式（8），可推导得到加载后孔隙比的计算式：

$\begin{align} {e_\text{sf}} &= {e_\text{sf - ap}} - \frac{{\Delta {V_\text{a}}}}{{{V_\text{s}}}} = {e_\text{sf - ap}} - \frac{{\Delta u\left[ {{e_\text{sf - ap}}(1 - {S_\text{r}}) + h{e_\text{sf - ap}}{S_\text{r}}} \right]}}{{\Delta u + {p_\text{ap}}}} \\ &= \frac{{(1 - h){\sigma _\text{v}}{e_\text{sf - ap}}{S_\text{r}} + {e_\text{sf - ap}}{p_\text{ap}}}}{{{\sigma _\text{v}} + {p_\text{ap}}}} 。 \end{align}$

(11)

（2）e₀ < e_th加载阶段

当e₀ < e_th，总应力增量由导致土骨架收缩的有效应力增量和引起孔隙封闭气相收缩的孔隙压力增量共同承担。该加载阶段的孔隙压力和封闭气相体积依然遵循Boyle’s定律和Henry’s定律（式（8）），故可推导得到孔隙比与孔隙压力的关系式：

${e}_\text{sf}=\frac{(1-h)u{e}_\text{sf-ap}{S}_\text{r}+{e}_\text{sf-ap}{p}_\text{ap}}{u+{p}_\text{ap}}。$

(12)

在完全侧限条件下，曹文贵等^[16]总结了大量压缩试验数据，认为土体竖向应变ε_v与所受竖向有效应力 ${\sigma '_{\text{v}}}$ 呈双曲线关系：

${\varepsilon _\text{v}} = \frac{{{{\sigma '}_\text{v}}}}{{a + b{{\sigma '}_\text{v}}}}。$

(13)

式中：a，b为具有明确物理含义的压缩常数，可由一维压缩试验结果拟合得到。

由土力学知识可推导一维压缩后孔隙比与竖向应变的关系式：

${e_\text{sf}} = \frac{{(1 - {\varepsilon _\text{v}})({V_\text{p}} + {V_\text{s}})}}{{{V_\text{s}}}} - 1 = (1 - {\varepsilon _\text{v}})(1 + {e_\text{th}}) - 1 。$

(14)

将式（14）代入式（15）可得到孔隙比与竖向有效应力的关系式：

${e_\text{sf}} = {e_\text{th}} - \frac{{{{\sigma '}_\text{v}}(1 + {e_\text{th}})}}{{(a + b{{\sigma '}_\text{v}})}}。$

(15)

不排水压缩过程中，孔隙内封闭气相和土骨架同时收缩，达到一种新的平衡状态，由于改良土与外界无任何气、水交换，土骨架收缩量严格等于孔隙内封闭气相收缩量。因此，由孔隙压力计算得到的孔隙比（式（12））与由竖向有效应力计算得到的孔隙比（式（15））恒定相等。进而根据太沙基有效应力原理，推导得竖向有效应力的二次函数表达式：

$A{\sigma '}_{\text{v}}^{\text{2}} + B{{\sigma '}_\text{v}} + C = 0。$

(16)

式中：

$A = \left( {1 - b} \right){e_\text{th}} + \left( {1 - h} \right)b{e_\text{sf - ap}}{S_\text{r}} + 1\text{；}$

(17)

$\begin{align} B =& \left[ {{e_\text{th}}(b - 1) - 1} \right]({\sigma _\text{v}} + {p_\text{ap}}) + \\ &(a - b{\sigma _\text{v}})(1 - h){e_\text{sf - ap}}{S_\text{r}} - b{e_\text{sf - ap}}{p_\text{ap}} - a{e_\text{th}}\text{；} \end{align}$

(18)

$C = ({\sigma _\text{v}} + {p_\text{ap}})a{e_\text{th}} - \left[ {(1 - h){\sigma _\text{v}}{S_\text{r}} + {p_\text{ap}}} \right]a{e_\text{sf - ap}}。$

(19)

通过求根法求解式（16）可得到竖向有效应力，由于竖向总应力恒大于竖向有效应力，竖向有效应力的唯一解为

${\sigma '_\text{v}} = \frac{{ - B - \sqrt {{B^2} - 4AC} }}{{2A}}。$

(20)

综上，当e≥e_th，压缩行为由封闭气相收缩控制，有效应力为0；当e < e_th，压缩行为由有效应力增量引起的土骨架收缩控制，可由式（20）求解某竖向总应力对应的竖向有效应力。竖向有效应力的统一表达式为

$\left. \begin{array}{l}{\sigma }_\text{v, th}=\frac{{p}_\text{ap}({e}_\text{sf-ap}-{e}_\text{th})}{{e}_\text{th}-(1-h){e}_\text{sf-ap}{S}_\text{r}}\text{，}\\ {{\sigma }^{\prime }}_\text{v}=\left\{ \begin{array}{l}0\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;（{e}_\text{sf}\ge {e}_\text{th}, {}^{}{\sigma }_\text{v}\le {\sigma }_\text{v, th}）\\ \frac{-B-\sqrt{{B}^{2}-4AC}}{2A}\;\;\;\;\;\;（{e}_\text{sf} < {e}_\text{th}{, }^{}{\sigma }_\text{v} > {\sigma }_\text{v, th}）\end{array} \right.\text{ }。\end{array}\right\}$

(21)

式中：

$\begin{array}{l} A=(1-b){e}_\text{th}+(1-h)b{e}_\text{sf-ap}{S}_\text{r}+1，B=\left[{e}_\text{th}(b-1)-1\right]\cdot\\ ({\sigma _\text{v}} + {p_\text{ap}}) + (a-b{\sigma }_\text{v})(1-h){e}_\text{sf-ap}{S}_\text{r}-b{e}_\text{sf-ap}{p}_\text{ap}-a{e}_\text{th}\text{，} \\ C = ({\sigma _\text{v}} + {p_\text{ap}})a{e_\text{th}} - \left[ {(1 - h){\sigma _\text{v}}{S_\text{r}} + {p_\text{ap}}} \right]a{e_\text{sf - ap}} 。 \end{array}$

对于存在有效应力的泡沫改良粗粒土，剪切过程中颗粒重排列和泡沫变形会引起孔隙压力的变化，从而影响土体的竖向有效应力和十字板剪切强度。泡沫改良粗粒土进入残余阶段后，有效应力状态趋于稳定，剪切应力不再损失^[11]。盾构改良渣土大变形下的剪切行为可通过泡沫改良粗粒土带压十字板剪切的残余剪切强度来表征^[7]。首先求解不排水剪切引起的超孔隙压力，然后基于Mohr-Coulomb准则建立残余状态下竖向有效应力与剪切强度的关系，最后构建剪切强度与剪切速率的相关性计算公式。

Skempton首先提出了孔隙压力系数来表征加载条件下土中孔隙压力的大小^[17]。不排水不排气三轴应力增量引起的孔隙压力增量Δu的一般表达式为

$\Delta u=\Delta {u}_{\mathrm{Ⅰ}}+\Delta {u}_{\mathrm{Ⅱ}}=B\left[\Delta {\sigma }_{3}+A(\Delta {\sigma }_{1}-\Delta {\sigma }_{3})\right] 。$

(22)

式中：Δu_Ⅰ，Δu_Ⅱ分别为压缩和剪切阶段应力增量引起的孔压增量；B为各向等应力（ $\Delta {\sigma _3}$ ）作用下的孔压系数；A为偏应力 $(\Delta {\sigma _1} - \Delta {\sigma _3})$ 作用下的孔压系数。

Bishop^[17]进一步推导得到孔压系数B：

$B=\frac{\Delta {u}_{\mathrm{Ⅰ}}}{\Delta {\sigma }_{3}}=\frac{1}{1+n\frac{{C}_\text{p}-{C}_\text{s}}{C-{C}_\text{s}}} 。$

(23)

式中：n为孔隙率；C_p为孔隙流体的压缩系数，即单位应力增量引起单位体积孔隙流体的体积变化；C_s为土颗粒的压缩系数；C为土骨架的压缩系数。

另外，孔隙压力系数A反映了土体的剪胀性，与土的性质和物理状态有关，一般由三轴试验确定。

类似地，泡沫改良粗粒土的加压式十字板剪切试验可划分为一维压缩阶段和十字板剪切阶段，尝试建立孔隙压力增量 $\Delta u'$ 的经验公式：

$\Delta u' = \Delta {u_\text{c}} + \Delta {u_\text{s}} = B'(\Delta {\sigma _\text{v}} + A'\Delta \tau )。$

(24)

式中：Δu_c，Δu_s分别为竖向总应力增量和剪切应力增量引起的孔隙压力增量； $A'$ 为十字板剪切作用下的孔压系数； $B'$ 为一维压缩作用下的孔压系数，

$B' = \frac{1}{{1 + \frac{{{e_\text{sf}}({{C'}_\text{p}} - {{C'}_\text{s}})}}{{(1 + {e_\text{sf}})(C' - {{C'}_\text{s}})}}}}。$

(25)

式中： ${C'_\text{s}}$ 为土颗粒的一维压缩系数，因土颗粒不可压缩，一般取0； ${C'_\text{p}}$ ， $C'$ 分别为孔隙流体、土骨架的一维压缩系数，与孔隙内气相体积有关。

当e≥e_th，竖向总应力全转化为孔隙压力， $B'$ = 1；当e < e_th，一维压缩系数随孔隙内气相体积变化而变化，参考式（23）推导得到

${C'_\text{p}} = \frac{{\partial {V_\text{p}}}}{{{V_\text{p}}\partial u}} = - \frac{{\partial {e_\text{sf}}}}{{{e_\text{sf}}\partial u}} = \frac{{{p_\text{ap}}{e_\text{sf - ap}}\left[ {1 - (1 - h){S_\text{r}}} \right]}}{{{e_\text{sf}}{{(u + {p_\text{ap}})}^2}}} .\text{，}$

(26)

$C' = \frac{{\partial V}}{{V\partial {{\sigma '}_\text{v}}}} = \frac{{ - {\text{ }}\partial e}}{{(1 + {e_\text{sf}})\partial {{\sigma '}_\text{v}}}} = \frac{a}{{(a + b{{\sigma '}_\text{v}})\left[ {a + (b - 1){{\sigma '}_\text{v}}} \right]}}。$

(27)

式中：V，V_p分别为改良土的总体积和孔隙体积。

在十字板剪切试验中，试样剪切应力变化特征与初始有效应力相关。因此，在未知剪切应力的前提下，可定义经验函数f（ ${\sigma '_\text{v0}}$ ）来表征剪切过程中初始有效应力对超孔隙压力的贡献，试样进入残余状态产生的总超孔隙压力ΔU_s可表示为

$\Delta {U_\text{s}} = B'\int_0^{{\tau _\text{r}}} {A'\text{d}\tau } = B'f({\sigma '_\text{v0}})。$

(28)

式中：τ_r为残余剪切强度； ${\sigma '_\text{v0}}$ 为剪切前试样的初始竖向有效应力，由式（21）可求得；f $({\sigma '_\text{v0}})$ 需由十字板剪切试验确定。

对于进入残余状态的泡沫改良粗粒土，进一步基于莫尔-库仑强度准则构建垂直于十字板竖向的剪切面对应的剪切强度τ_vr与残余阶段竖向有效应力σ_vr’的计算公式：

${\tau _\text{vr}} = {c_\text{r}} + {\sigma '_\text{vr}}\tan {\varphi '_{\text{r}}} = {c_\text{r}} + \left( {{{\sigma '}_\text{v0}} - \Delta {U_\text{s}}} \right)\tan {\varphi '_{\text{r}}}。$

(29)

式中：c_r， ${\varphi '_{\text{r}}}$ 分别为残余黏聚力和有效内摩擦角。

进一步为了描述泡沫改良粗粒土的连续流动特性，仍需考虑剪切速率对残余剪切强度的影响。对于软黏土或粗粒土材料，大量剪切试验^{[5, 11]}证实了它们的抗剪强度具有率相关性，或者剪切流动时呈现非牛顿流体的流变特性。Galli等^[5]、杨益等^[18]认为泡沫改良粗粒土是具有屈服应力的幂律流体，剪切速率对残余剪切强度的影响可表示为

${\tau _\text{vr}} = {\tau _\text{vr - ref}}\left[ {\delta + \kappa {{\left( {\frac{{\dot \gamma }}{{{{\dot \gamma }_\text{ref}}}}} \right)}^n}} \right]。$

(30)

式中：γ(·)_ref，τ_vr-ref为剪切速率参考值及其对应的残余剪切强度，δ，κ，n分别为表征剪切速率影响的无量纲常数。

其中，残余剪切强度参考值可由式（29）求解得到，所需的f( ${\sigma '_\text{v0}}$ )，c_r和 ${\varphi '_{\text{r}}}$ 由剪切速率参考值下十字板剪切试验结果拟合得到。综上，图 3展示了泡沫改良粗粒土不排水残余剪切强度的计算流程。

图 3 泡沫改良粗粒土残余剪切强度计算流程图

Figure 3. Calculation process of residual shear strength of foam-conditioned coarse-grained soils

下载: 全尺寸图片幻灯片

本次试验粗粒土材料的级配曲线如图 4所示。土样属于级配不良砾砂，其不均匀系数C_u为10.0，曲率系数C_c为0.7，最大孔隙比e_max为0.768，颗粒密度ρ_s为2.66 g/cm³。试验采用符合EFNARC（2005）要求的发泡设备生成泡沫，FER为12.0，半衰期为8.5 min。原始含水率为10%的土样在堆积状态下测得的孔隙比e_ps为0.762，孔隙水密度ρ_w取1 g/cm³。泡沫改良试样的常压下孔隙状态参数和试验工况（包括竖向总应力和十字板转速n）详见表 1。

图 4 试验粗粒土级配曲线

Figure 4. Grain-size distribution curve of testing soils

下载: 全尺寸图片幻灯片

表 1 泡沫改良粗粒土初始状态参数和试验工况

Table 1. Initial parameters and test conditions of specimens

FIR/%	α	e_sf-ap	S_r/%	σ_v/kPa	n/(r·min^-1)
20	0.186	0.827	35.7	50, 100, 200, 300	1/30, 1/5, 1, 3, 5, 25
30	0.366	0.955	32.4
40	0.523	1.130	28.7

下载: 导出CSV

| 显示表格

采用自主研发的加压旋转剪切仪^[6]来开展泡沫改良粗粒土的十字板剪切试验，仪器实物图如图 5所示。仪器可提供量程0~400 kPa，精度0.5%的竖向荷载；十字板扭矩传感器的量程为0~30 N·m，精度为0.1%；孔压表的量程为0~600 kPa，精度为0.5%。

图 5 加压式十字板剪切仪

Figure 5. Pressurized vane shear apparatus

下载: 全尺寸图片幻灯片

考虑到剪切失效面上两类剪切应力（式（1））分别与竖向、侧向有效应力相关，而砂土材料的K₀可由经典Jaky公式^[19]简易求得：

${K_0} = \frac{{{{\sigma '}_\text{l}}}}{{{{\sigma '}_\text{v}}}} = 1 - \sin \varphi '。$

(31)

联立式（1），（29），（31），可得试样进入残余状态后十字板扭矩与竖向有效应力的关系式，结合试验结果可拟合得残余阶段黏聚力和有效内摩擦角。

$\begin{align} {T_\text{r}} &= \left( {\frac{7}{6} - \sin {{\varphi '}_\text{r}}} \right)\pi {D^3}{\tau _\text{vr}} \\ &= \left( {\frac{7}{6} - \sin {{\varphi '}_\text{r}}} \right)\pi {D^3} \cdot \left[ {{c_\text{r}} + ({{\sigma '}_\text{v0}} - \Delta {U_\text{s}})\tan {{\varphi '}_{\text{r}}}} \right]。 \end{align}$

(32)

每次试验均按照以下步骤进行：

（1）将质量相等的泡沫改良粗粒土试样分3层均匀地装入试样腔中，每层厚度约6 cm；安装十字板和顶板后，开启排气阀，施加预压力，以排出试样上方空气，而后关闭排气阀，确保试验的不排水环境。

（2）通过顶板向试样施加目标竖向压力，当监测顶板位置的位移计读数恒定时，十字板开始以表 1设计的速度旋转剪切试样，单次试验的总旋转角度为1440°，以确保试样进入残余状态。试验全程记录试样的竖向变形量、孔隙压力以及十字板扭矩。

以十字板转速1/30 rpm对应的剪切速率为参考值，由式（2）计算得参考剪切速率为7.24×10^-3 s^-1。采用FIR为20%和40%的试样在不同压力下的试验结果来拟合确定计算模型所需的材料常数。图 6展示了一维压缩后初始孔隙比与竖向有效应力的关系，当e₀ ≥ e_max，有效应力未产生；当e₀ < e_max，初始孔隙比与竖向有效应力呈非线性关系。因此，孔隙比临界值e_th等于e_max.。进而通过式（15）拟合e < e_th曲线求解压缩常数a，b，分别为358.66 kPa，8.48。

图 6 一维压缩后初始孔隙比随竖向有效应力的变化

Figure 6. Variation of initial void ratio against vertical effective stress under compression

下载: 全尺寸图片幻灯片

由式（28）可知，需求解剪切引起的超孔隙压力ΔU_s与孔压系数 $B'$ 的比值和初始竖向有效应力的相关性来获得经验函数 $f({\sigma '_\text{v0}})$ 。如所示，两者呈线性关系，因此 $f({\sigma '_\text{v0}})$ 为线性函数，由拟合结果可知f ( ${\sigma '_\text{v0}}$ )=0.662 ${\sigma '_\text{v0}}$ 。图 8展示了试样进入残余状态后十字板扭矩和竖向有效应力的关系。由式（32）拟合得到残余强度参数c_r，φ_r分别为1.03 kPa，12.50°。然后，采用FIR为20%和40%的试样在不同十字板转速下的试验结果来拟合表征剪切速率影响的无量纲常数。展示了 ${\tau _{{\text{vr}}}}/{\tau _{{\text{vr - ref}}}}$ 随 $\gamma/{\gamma _{{\text{ref}}}}$ 的变化。由式（30）拟合得δ，κ，n分别为0.969，0.029，0.269。

图 7 剪切引起超孔隙压力随初始有效应力的变化

Figure 7. Variation of shear-induced excess pore pressure against initial effective stress

下载: 全尺寸图片幻灯片

图 8 残余状态十字板扭矩与竖向有效应力关系

Figure 8. Variation of vane torque against vertical effective stress

下载: 全尺寸图片幻灯片

图 9 无量纲十字板残余剪切强度随剪切速率的变化

Figure 9. Variation of residual strength against shear strain rate

下载: 全尺寸图片幻灯片

综上，图 6~9反映了各拟合公式均具有较高的拟合优度R²（R² > 0.85），且平均绝对误差MAE和均方根误差RMSE远小于目标值，证明了理论模型中新建的式（15），（28）~（30）的合理性。

采用FIR为30%的试样在不同竖向压力和十字板转速下的试验结果来验证计算模型的准确性。图 10对比了残余强度试验值（以空心符号表征）和理论值（以实心符号表征），试验值和理论值的平均误差约为10%，MAE为0.32，RMSE为0.59。表明计算模型能较好预测泡沫改良粗粒土的率相关残余强度，根据特定工况试验结果拟合得到的材料常数可预测任意竖向压力、剪切速率、FIR下的残余强度。

图 10 FIR=30%改良土残余剪切强度试验值与理论值对比

Figure 10. Comparison between analytical and experimental values of residual strength of conditioned soils with FIR=30%

下载: 全尺寸图片幻灯片

进一步地，采用该理论模型来计算竖向总应力和FIR对泡沫改良粗粒土残余剪切强度的影响。对于本研究的试验材料，将粗粒土和泡沫的已知物理参数和拟合得到的材料常数代入理论模型中，另拟定原始含水率取10%、发泡倍率取12、剪切速率取参考值7.24×10^-3 s^-1。应注意到泡沫扩大系数决定了相同FIR泡沫改良粗粒土的孔隙状态，取决于泡沫与粗粒土的混合状态。为聚焦竖向总应力和泡沫注入比的影响，设定扩大系数为0.5。

图 11给出了残余强度随竖向总应力和有效应力的变化情况。由图 11（a）可知，残余强度随着竖向总应力的增加而增加，随着泡沫注入比的增加而减少。而且随着泡沫注入比增加，残余强度随竖向总应力的变化率越小，与试验得到的规律一致。由图 11（b）可知，即便竖向总应力和泡沫注入比变化，残余强度仍与竖向有效应力呈正相关关系。

图 11 残余剪切强度随竖向总应力和有效应力的变化

Figure 11. Variation of analytical residual strength against vertical total stress and effective stress

下载: 全尺寸图片幻灯片

泡沫改良粗粒土的宏观流变行为与剪切变稀流体类似，进一步计算表征残余状态下带压泡沫改良粗粒土宏观流变行为的流变参数，图 12展示了不同竖向总应力和泡沫注入比下的屈服应力和塑性黏度。屈服应力和塑性黏度均随竖向总应力增加而增加，而随泡沫注入比增加而减小，然而相较于屈服应力，塑性黏度变化较小。

图 12 泡沫改良粗粒土残余状态下大变形屈服应力和塑性黏度随竖向总应力的变化

Figure 12. Variation of analytical yield stress and plastic viscosity of conditioned soils in large deformation against vertical total stress

下载: 全尺寸图片幻灯片

构建了预测带压泡沫改良粗粒土率相关残余剪切强度的统一理论模型，可用于定量评估不同土舱压力、改良参数的泡沫改良粗粒土剪切大变形后的流变特性。

（1）在一维加载作用下，泡沫改良粗粒土存在两类孔隙比状态。当e ≥ e_th，有效应力不存在；当e < e_th，基于Boyle定律、Henry定律和土应力-应变关系的双曲线假设，推导竖向有效应力。然后推导不排水剪切引起的超孔隙压力，并依据Mohr-Coulomb准则给出残余状态竖向有效应力与剪切强度的计算公式，最后建立表征残余剪切强度与剪切速率的计算公式。

（2）不排水加压十字板剪切试验结果与理论解对比发现，根据某些特定工况的试验结果拟合得到的材料常数适用于外推预测其它改良参数下泡沫改良粗粒土的残余剪切强度，平均误差约10%。该统一理论模型量化考虑了总应力、有效应力以及泡沫注入量的影响，能克服试验获得的改良土流变参数仅能反映某一特定工况下改良土流动特性的缺点。

（3）采用提出的理论模型进行影响因素分析，不排水残余剪切强度随竖向总应力增加而增加，随泡沫注入比增加而减小，然而它与竖向有效应力的正相关关系不随竖向总应力或泡沫注入比而改变。带压泡沫改良粗粒土的屈服应力和塑性黏度均随竖向总应力增加而增加，随泡沫注入比增加而减小。

相比于既有的由试验拟合的改良土流变模型及参数，本文的计算模型量化考虑了竖向总应力、FIR对泡沫改良粗粒土残余剪切强度和残余状态流变行为的影响，为盾构进排土CFD数值模型在渣土改良策略优化中的应用提供了材料参数计算依据。

图 1 试验材料的颗粒级配曲线

Figure 1. Grain-size distribution curves of test materials

下载: 全尺寸图片幻灯片

图 2 试验装置及分层砾石-砂粉土柱试样

Figure 2. Test devices and layered gravel-sandy silt column specimen

下载: 全尺寸图片幻灯片

图 3 颗粒悬浮和沉积照片

Figure 3. Photos of suspension and deposition of particles

下载: 全尺寸图片幻灯片

图 4 不同加载条件下泥浆浊度时变曲线

Figure 4. Variation of slurry turbidity with time under different loading conditions

下载: 全尺寸图片幻灯片

图 5 颗粒迁移定量化

Figure 5. Quantification of particle migration

下载: 全尺寸图片幻灯片

图 6 试验结束时砾石底层形成的泥化夹层

Figure 6. Generation of interlayer in gravel bottom layer at end of tests

下载: 全尺寸图片幻灯片

图 7 颗粒再悬浮增量随重复次数的变化

Figure 7. Variation of re-suspension increment of particles with repetition number

下载: 全尺寸图片幻灯片

图 8 试样不同深度处超孔压随时间的变化

Figure 8. Variation of excess pore pressure with time at different depths of specimen

下载: 全尺寸图片幻灯片

图 9 不同层深处砂粉土的液化度时变曲线

Figure 9. Variation of liquefaction degree with time at different depths of sandy silt layer

下载: 全尺寸图片幻灯片

图 10 超孔压振荡幅值随重复次数的变化

Figure 10. Variation of oscillation amplitude of excess pore pressure with repetition number

下载: 全尺寸图片幻灯片

图 11 颗粒迁移质量随超孔压振幅的变化

Figure 11. Variation of migration mass of particles with amplitude of excess pore pressure

下载: 全尺寸图片幻灯片

图 12 颗粒迁移质量随平均液化度的变化

Figure 12. Variation of migration mass of particles with average liquefaction degree

下载: 全尺寸图片幻灯片

图 13 不同深度处颗粒悬浮质量随加载强度的变化

Figure 13. Variation of suspension mass of particles with loading strength at different depths

下载: 全尺寸图片幻灯片

表 1 动静组合加载试验方案

Table 1 Test schemes of combined dynamic and static loads

下载: 导出CSV

参考文献(35)

[1]	中华人民共和国交通运输部. 李小鹏出席国新办"权威部门话开局"新闻发布会, 强调: 奋力加快建设交通强国努力当好中国现代化的开路先锋[EB/OL]. 2023-02-24. https://www.mot.gov.cn/jiaotongyaowen/202302/t20230224_3763233.html. Ministry of Transport of the People's Republic of China. Li Xiaopeng attended the press conference of the "authority department speech opening" of the state information office and stressed: strive to speed up the construction of transportation power and strive to be the pioneer of china's modernization. [EB/OL]. 2023-02-24. https://www.mot.gov.cn/jiaotongyaowen/202302/t20230224_3763233.html. ) (in Chinese)
[2]	张军辉, 彭俊辉, 郑健龙. 路基土动态回弹模量预估进展与展望[J]. 中国公路学报, 2020, 33(1): 1-13. doi: 10.3969/j.issn.1001-7372.2020.01.001 ZHANG Junhui, PENG Junhui, ZHENG Jianlong. Progress and prospect of the prediction model of the resilient modulus of subgrade soils[J]. China Journal of Highway and Transport, 2020, 33(1): 1-13. (in Chinese) doi: 10.3969/j.issn.1001-7372.2020.01.001
[3]	彭惠, 马巍, 穆彦虎, 等. 青藏公路普通填土路基长期变形特征与路基病害调查分析[J]. 岩土力学, 2015, 36(7): 2049-2056. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507037.htm PENG Hui, MA Wei, MU Yanhu, et al. Analysis of disease investigation and long-term deformation characteristics of common fill embankment of the Qinghai-Tibet Highway[J]. Rock and Soil Mechanics, 2015, 36(7): 2049-2056. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507037.htm
[4]	NGUYEN T T, INDRARATNA B, KELLY R. Mud pumping under railtracks: mechanisms, assessments and solutions[J]. Australian Geomechanics, 2019, 54(4): 59-80.
[5]	杨新安, 高艳灵. 沪宁铁路翻浆冒泥病害的地质雷达检测[J]. 岩石力学与工程学报, 2004, 23(1): 116-119. doi: 10.3321/j.issn:1000-6915.2004.01.022 YANG Xinan, GAO Yanling. GPR inspection for Shanghai-nanjing railway trackbed[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(1): 116-119. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.01.022
[6]	MNEINA A, SHALABY A. Relating gradation parameters to mechanical and drainage performance of unbound granular materials[J]. Transportation Geotechnics, 2020, 23(6): 100315. http://www.xueshufan.com/publication/2999576606
[7]	边学成, 李书豪, 万章博, 等. 路基注浆对高速铁路轨道-路基体系动力特性的影响[J]. 振动与冲击, 2022, 41(4): 294-302. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ202204038.htm BIAN Xuecheng, LI Shuhao, WAN Zhangbo, et al. Influence of injection remediation on dynamic behaviors of a high-speed railway track-subgrade system[J]. Journal of Vibration and Shock, 2022, 41(4): 294-302. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ202204038.htm
[8]	KERMANI B, XIAO M, STOFFELS S M, et al. Reduction of subgrade fines migration into subbase of flexible pavement using geotextile[J]. Geotextiles and Geomembranes, 2018, 46(4): 377-383. doi: 10.1016/j.geotexmem.2018.03.006
[9]	INDRARATNA B, BABAR SAJJAD M, NGO T, et al. Improved performance of ballasted tracks at transition zones: a review of experimental and modelling approaches[J]. Transportation Geotechnics, 2019, 21: 100260. doi: 10.1016/j.trgeo.2019.100260
[10]	HAYASHI S, SHAHU J T. Mud pumping problem in tunnels on erosive soil deposits[J]. Géotechnique, 2000, 50(4): 393-408. doi: 10.1680/geot.2000.50.4.393
[11]	BEDRIKOVETSKY P, CARUSO N. Analytical model for fines migration during water injection[J]. Transport in Porous Media, 2014, 101(2): 161-189. doi: 10.1007/s11242-013-0238-7
[12]	WANK J R, GEORGE S M, WEIMER A W. Vibro-fluidization of fine boron nitride powder at low pressure[J]. Powder Technology, 2001, 121(2/3): 195-204.
[13]	MAWATARI Y, KOIDE T, TATEMOTO Y, et al. Effect of particle diameter on fluidization under vibration[J]. Powder Technology, 2002, 123(1): 69-74. doi: 10.1016/S0032-5910(01)00432-6
[14]	周健, 杨永香, 刘洋, 等. 循环荷载下砂土液化特性颗粒流数值模拟[J]. 岩土力学, 2009, 30(4): 1083-1088. doi: 10.3969/j.issn.1000-7598.2009.04.039 ZHOU Jian, YANG Yongxiang, LIU Yang, et al. Numerical modeling of sand liquefaction behavior under cyclic loading[J]. Rock and Soil Mechanics, 2009, 30(4): 1083-1088. (in Chinese) doi: 10.3969/j.issn.1000-7598.2009.04.039
[15]	ZHOU L F, WANG J W, GE W, et al. Quantifying growth and breakage of agglomerates in fluid-particle flow using discrete particle method[J]. Chinese Journal of Chemical Engineering, 2018, 26(5): 914-921. doi: 10.1016/j.cjche.2017.05.018
[16]	ZHOU K, HOU J, SUN Q C, et al. A study on particle suspension flow and permeability impairment in porous media using LBM-DEM-IMB simulation method[J]. Transport in Porous Media, 2018, 124(3): 681-698. doi: 10.1007/s11242-018-1089-z
[17]	ALOBAIDI I, HOARE D. Factors affecting the pumping of fines at the subgrade subbase interface of highway pavements: a laboratory study[J]. Geosynthetics International, 1994, 1(2): 221-259. doi: 10.1680/gein.1.0010
[18]	KERMANI B, XIAO M, STOFFELS S M, et al. Measuring the migration of subgrade fine particles into subbase using scaled accelerated flexible pavement testing–a laboratory study[J]. Road Materials and Pavement Design, 2019, 20(1): 36-57. doi: 10.1080/14680629.2017.1374995
[19]	蔡袁强, 严舒豪, 曹志刚, 等. 交通荷载下粉质黏土路基翻浆冒泥机理试验[J]. 吉林大学学报(工学版), 2021, 51(5): 1742-1748. https://www.cnki.com.cn/Article/CJFDTOTAL-JLGY202105022.htm CAI Yuanqiang, YAN Shuhao, CAO Zhigang, et al. Experiments to investigate mechanism of mud pumping of road base on silty clay soil under cyclic loading[J]. Journal of Jilin University (Engineering and Technology Edition), 2021, 51(5): 1742-1748. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JLGY202105022.htm
[20]	韩博文, 蔡国庆, 李舰, 等. 有砟轨道路基翻浆冒泥模型试验系统的研发与应用[J]. 岩土工程学报, 2022, 44(8): 1406-1415. doi: 10.11779/CJGE202208005 HAN Bowen, CAI Guoqing, LI Jian, et al. Development and application of model test system for mud pumping in ballasted track subgrade[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(8): 1406-1415. (in Chinese) doi: 10.11779/CJGE202208005
[21]	白冰, 张鹏远, 宋晓明, 等. 渗透作用下多孔介质中悬浮颗粒的迁移过程研究[J]. 岩土工程学报, 2015, 37(10): 1786-1793. doi: 10.11779/CJGE201510006 BAI Bing, ZHANG Pengyuan, SONG Xiaoming, et al. Transport processes of suspended particles in saturated porous media by column seepage tests[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(10): 1786-1793. (in Chinese) doi: 10.11779/CJGE201510006
[22]	张升, 高峰, 陈琪磊, 等. 砂-粉土混合料在列车荷载作用下细颗粒迁移机制试验[J]. 岩土力学, 2020, 41(5): 1591-1598. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202005016.htm ZHANG Sheng, GAO Feng, CHEN Qilei, et al. Experimental study of fine particles migration mechanism of sand-silt mixtures under train load[J]. Rock and Soil Mechanics, 2020, 41(5): 1591-1598. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202005016.htm
[23]	ZHANG S, GAO F, HE X Z, et al. Experimental study of particle migration under cyclic loading: effects of load frequency and load magnitude[J]. Acta Geotechnica, 2021, 16(2): 367-380. doi: 10.1007/s11440-020-01137-x
[24]	GAO F, HE X Z, ZHANG S. Pumping effect of rainfall-induced excess pore pressure on particle migration[J]. Transportation Geotechnics, 2021, 31: 100669. doi: 10.1016/j.trgeo.2021.100669
[25]	GAO F, ZHANG S, HE X Z, et al. Experimental study on migration behavior of sandy silt under cyclic load[J]. ASCE-Journal of Geotechnical and Geoenvironmental Engineering, 2022, 148(5): 06022003. doi: 10.1061/(ASCE)GT.1943-5606.0002796
[26]	Standard Test Methods for Particle-Size Distribution (Gradation) of Soils Using Sieve Analysis: ASTM D6913/D6913M-17[S]. 2017.
[27]	TASALLOTI A, MARSHALL A M, HERON C M, et al. Geocellular railway drainage systems: physical and numerical modelling[J]. Transportation Geotechnics, 2020, 22: 100299. doi: 10.1016/j.trgeo.2019.100299
[28]	SHENG D, ZHANG S, NIU F, et al. A potential new frost heave mechanism in high-speed railway embankments[J]. Géotechnique, 2014, 64(2): 144-154. doi: 10.1680/geot.13.P.042
[29]	BIAN X C, LI W, HU J, et al. Geodynamics of high-speed railway[J]. Transportation Geotechnics, 2018, 17: 69-76. doi: 10.1016/j.trgeo.2018.09.007
[30]	LAWRENCE M, BULLOCK R, LIU Z M. China's High-Speed Rail Development[M]. Washington D C: World Bank, 2019.
[31]	TERZAGHI K V. The Shearing resistance of saturated soils and the angle between the planes of shear[C]//Proceedings of the 1th International Conference on Soil Mechanics and Foundation Engineering, Cambridge, 1936: 54-56.
[32]	MCDOUGAL W G, TSAI Y T, LIU P L F, et al. Wave-induced pore water pressure accumulation in marine soils[J]. Journal of Offshore Mechanics and Arctic Engineering, 1989, 111(1): 1-11. doi: 10.1115/1.3257133
[33]	BRAY J D, SANCIO R B. Assessment of the liquefaction susceptibility of fine-grained soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(9): 1165-1177. doi: 10.1061/(ASCE)1090-0241(2006)132:9(1165)
[34]	LI J E, ZHANG J H, ZHANG A S, et al. Evaluation on deformation behavior of granular base material during repeated load triaxial testing by discrete-element method[J]. International Journal of Geomechanics, 2022, 22(11): 1-11.
[35]	WINTERWERP J C, DE BOER G J, GREEUW G, et al. Mud-induced wave damping and wave-induced liquefaction[J]. Coastal Engineering, 2012, 64: 102-112. doi: 10.1016/j.coastaleng.2012.01.005

施引文献(4)

平均动应力 ${\sigma _m}$ /kPa

动应力幅值 ${\sigma _d}$ /kPa

动静组合加载重复次数M

动荷载单次循环次数N /10³

静荷载单次间歇时长t/min

1.	韩博文，蔡国庆，苏彦林，单冶鹏，李舰. 间歇荷载–湿化耦合作用下有砟轨道路基翻浆冒泥机制及动力特性试验研究. 岩石力学与工程学报. 2025(01): 69-80 . 百度学术
2.	蒋红光，王新宇，马川义，张宁，刘舜，王川. 高速铁路基床细颗粒动力迁移机制与临界水力梯度研究. 振动与冲击. 2025(06): 263-271 . 百度学术
3.	高峰，曾宪璋，钟闻华，黄生勇，张军辉. 多年冻土区道路工程病害处治技术研究进展与展望. 中外公路. 2024(05): 1-16 . 百度学术