填充沟对地铁运营引起地基振动隔振效果研究

蔡袁强; 江乾明; 曹志刚; 袁宗浩; 楼佳悦; 张盈哲

doi:10.11779/CJGE202008002

填充沟对地铁运营引起地基振动隔振效果研究 English Version

蔡袁强^{1, 2,},
江乾明¹,
曹志刚^1, ,,
袁宗浩²,
楼佳悦³,
张盈哲³

1.
浙江大学滨海和城市岩土工程研究中心,浙江　杭州　310058
2.
浙江工业大学岩土工程研究所,浙江　杭州310014
3.
浙江华云电力工程设计咨询有限公司,浙江　杭州　310014

基金项目:

国家重点研发计划项目 2016YFC0800200

国家自然科学基金项目 51778571

国家自然科学基金项目 51708503

省属高校基本科研业务费项目（科技类） RF-A2019014

详细信息

作者简介:
蔡袁强(1965—),男,博士,教授,主要从事岩土工程科研与教学工作。E-mail：caiyq@zju.edu.cn

通讯作者:
曹志刚, E-mail：caozhigang2011@zju.edu.cn

中图分类号: TU435
计量
- 文章访问数: 342
- HTML全文浏览量: 43
- PDF下载量: 203
出版历程
- 收稿日期: 2019-09-25
- 网络出版日期: 2022-12-05
- 刊出日期: 2020-07-31

Vibration-isolation effectiveness of in-filled trench on ground vibration induced by subway operation

1.
Research Center of Coastal and Urban Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China
2.
Institute of Geotechnical Engineering, Zhejiang University of Technology, Hangzhou 310014, China
3.
Zhejiang Huayun Electric Power Engineering Design & Consultation Co., Ltd., Hangzhou 310014, China

摘要

摘要: 为分析填充沟对地铁运营引起的地基振动隔振效果,建立了地铁列车–轨道结构–衬砌–地基–填充沟2.5维有限元分析模型。列车荷载模拟为一组符合列车轴载空间分布的常荷载,钢轨和浮置板简化为无限长Euler梁,衬砌与地基模拟为均质弹性体。理论模型通过钢轨与衬砌仰拱处、以及衬砌与地基处应力和位移连续条件进行耦合。混凝土填充沟模拟为地基中异质体,填充沟与地基交界面通过共用节点法处理,并在波数域中进行求解,最后通过快速Fourier逆变换（IFFT）进行波数展开获得三维时域–空间域内的地基动力响应。计算分析了填充沟在不同列车速度、填充沟沟深、沟宽与隧道埋深等情况下的隔振效果。研究表明,填充沟隔振效果随沟深增加而显著改善,但隔振效果改善速率在填充沟超过一定深度后变缓；随着填充沟距地铁距离增加,沟前隔振效果减弱,而沟后隔振效果改善。随着列车速度增加,填充沟隔振效果变好。
- 填充沟 /
- 地铁荷载 /
- 2.5维有限元 /
- 隔振效果 /
- 移动荷载
Abstract: In order to analyze the isolation effectiveness of in-filled trench to the ground vibration caused by subway operation, a coupling model, which is divided into five parts including moving train, track structure, lining, ground and in-filled trench, is established. The train load is simulated by a set of constant loads consistent with the spatial distribution of train axles. The rails and floating plates are simplified as the infinite Euler beams. The lining and ground are simulated as the homogeneous elastic media. The theoretical model is coupled by the continuous conditions for the stresses and displacements between the rail and the lining invert and those between the lining and the ground. The concrete in-filled trench is modelled as the heterogeneous medium in the ground. The interface between the in-filled trench and the ground is treated by the common node method and solved in the wave number domain. Finally, the three-dimensional dynamic response in the time-space domain is obtained by the fast inverse Fourier transform (IFFT). The vibration-isolation effectiveness of the filled trench under different train speeds, in-filled trench depths, trench widths and tunnel depths is calculated and analyzed. The results show that its vibration-isolation effectiveness is improved with the increase of the in-filled trench depth, but when the depth of the in-filled trench exceeds a certain depth, the improvement rate of the vibration-isolation effectiveness slows down. With the increase of the distance between the in-filled trench and the subway, the isolation effectiveness decreases in the area before the in-filled trench and increases in the area after the in-filled trench. With the increase of train speed, the vibration-isolation effectiveness becomes better.
- in-filled trench /
- subway load /
- 2.5D finite element /
- vibration isolation effectiveness /
- moving load

HTML全文

0. 引言

微生物岩土加固技术是利用微生物作用生成的有机、无机物材料填充孔隙和胶结岩土颗粒，以改善岩土体的工程性能的技术方法^[1-2]。相较传统的Portland水泥或高聚物加固方法，该技术具有对土体环境影响小、可持续性强的优点。微生物加固技术中研究最多、应用最广的为微生物诱导碳酸钙沉淀(microbially induced calcite precipitation, MICP)技术。MICP技术中高产脲酶芽孢杆菌使用最为普遍，其基本原理是利用自然界中广泛存在的无毒害的高产脲酶芽孢杆菌将尿素水解成铵根离子和碳酸根离子，然后碳酸根离子与环境中的钙离子结合生成碳酸钙沉淀^[1-2]。目前，MICP技术能够增大土体强度和刚度，封堵裂隙等，因而可广泛应用于提高砂土抗液化能力^[3-4]，降低土体渗透性以改善其抗渗透变形能力^[5-6]，冻土地区残积土加固^[7]，三合土表面耐水性处理^[8]，尾矿黏土处理^[9]，泥炭土和有机质黏土改良^[10-11]，沙丘固化与沙漠治理^[12]以及岩体裂隙加固^[13]等。

针对MICP加固砂土的力学性能，目前已有学者做了大量相关试验研究。Feng等^[14]通过三轴排水试验研究了围压和胶结程度对MICP加固石英砂力学特性的影响。结果表明，MICP加固砂的强度、剪胀性和初始弹性模量与胶结程度呈正相关。Liu等^[15]对MICP加固钙质砂进行了相关试验研究，表明碳酸钙含量增加并不会引起被加固土体的峰值内摩擦角改变，但是黏聚力会显著提高，这与MICP加固石英砂^[16-17]的试验结果不同。目前，针对MICP加固砂土的本构理论研究较试验研究少。方祥位等^[18]基于损伤力学理论建立了MICP加固珊瑚砂的损伤本构模型。Gai等^[19]基于次加载面理论建立了微生物加固砂土的本构模型，但其未考虑MICP胶结作用对土体黏聚力的影响。

本文基于临界状态土力学理论框架，采用非关联流动法则，建立状态相关的MICP加固砂土弹塑性本构模型。随后，将建立的模型分别模拟MICP加固石英砂和钙质砂的三轴排水剪切试验。计算结果表明，该模型可以较好地模拟不同加固程度和不同固结围压下试样的应力–应变关系和剪胀特性。

1. MICP加固砂土的力学特性分析

1.1 MICP加固与破坏机理

MICP加固砂土时，生成的碳酸钙晶体一般有3种沉积形式：①碳酸钙晶体附着在砂颗粒表面，使砂颗粒表面变粗糙，示意图和SEM(scanning electron microscopy)图^[20]如图 1(a)和图 1(b)所示；②碳酸钙晶体在砂颗粒间逐渐沉积生长成较大的晶簇，并将砂颗粒连接起来，形成桥接作用，示意图和SEM图^[20]如图 1(c)和图 1(d)所示；③碳酸钙晶体附着在相邻的砂颗粒接触处，将其胶结起来，示意图和SEM图^[20]如图 1(e)和图 1(f)所示。

图 1 MICP加固砂土机理示意图与SEM图^[20]

Figure 1. Schematic diagrams and SEM images^[20] of MICP-treated sands

下载: 全尺寸图片幻灯片

碳酸钙晶体的沉积形式会对MICP加固砂土的力学行为产生显著影响：第一种沉积形式可以填充砂颗粒间孔隙并增加其表面粗糙度，使土体内摩擦角增大，剪胀增加；第二种和第三种碳酸钙沉积形式主要是对砂土颗粒起到胶结作用，可以显著提高土体强度和剪胀性。对MICP加固Ottawa20/30砂进行三轴剪切后重塑装样再剪切试验，发现试样经过多次剪切之后强度与未加固砂基本相同，但是剪胀和刚度较未加固砂仍有提高^[21]。MICP加固砂土的破坏机理可以假设为：剪切过程中碳酸钙的胶结作用逐渐破坏但附着在砂颗粒表面的碳酸钙未被完全磨损掉，同时胶结破坏后的碳酸钙转化为第一种沉积形式。胶结作用退化造成强度降低，出现应变软化现象和剪胀。当胶结作用完全丧失后，附着在砂颗粒表面的碳酸钙仍会使其表面粗糙度增加，刚度和剪胀增加，但是对强度贡献较小。

1.2 破坏包络线

胶结土体的破坏包络线一般通过对同一胶结程度，不同围压下三轴剪切试样的峰值强度进行拟合得到，通过破坏包络线可以确定MICP加固砂土的初始黏聚力和峰值摩擦角。图 2分别给出了MICP加固Ottawa20/30砂^[22]，MICP加固石英砂^[23]和MICP加固钙质砂^[24]在 $p' - q$ 空间中的峰值应力点和破坏包络线。

图 2 试验峰值应力点和破坏包络线

Figure 2. Test peak stress points and failure envelopes

下载: 全尺寸图片幻灯片

中的 ${B_{{\text{ca}}}}$ 和 ${\text{CCC}}$ 分别为MICP加固石英砂和MICP加固钙质砂中生成的碳酸钙的质量分数。对于MICP加固石英砂，采用酸洗法去除生成的碳酸钙并根据质量差值计算得到碳酸钙的质量分数^[23]：

${B_{{\text{ca}}}} = ({m_{\text{i}}} - {m_{\text{a}}})/{m_{\text{i}}} \times 100\% \text{，}$

(1)

式中， ${m_{\text{i}}}$ 为MICP加固后的试样干重， ${m_{\text{a}}}$ 为酸洗后的试样干重。

对于MICP加固钙质砂，由于钙质砂主要成分为碳酸钙，所以不能采用酸洗法计算生成的碳酸钙质量，一般直接用MICP加固试样和未加固试样的干重差值计算得到碳酸钙质量分数^[24]：

${\text{CCC}} = ({m_{{\text{MICP}}}} - {m_{{\text{un}}}})/{m_{{\text{un}}}} \times 100\% \text{，}$

(2)

式中， ${m_{{\text{MICP}}}}$ 为MICP加固后的试样干重， ${m_{{\text{un}}}}$ 为未加固试样干重。

从可以看出，由于碳酸钙的胶结作用，MICP加固之后的Ottawa20/30砂，石英砂和钙质砂的黏聚力均增加。MICP加固钙质砂在 $p' - q$ 空间破坏包络线的斜率减小，峰值摩擦角相较未加固钙质砂明显减小，可能原因为钙质砂颗粒表面存在很多孔隙，且砂颗粒形状多棱角状，MICP加固过程中生成的碳酸钙将钙质砂颗粒表面的内孔隙填充，使其表面更加规整导致摩擦角减小^[24]。

1.3 临界状态线

分别给出了在 $e$ – $\ln p'$ 和 $p'$ – $q$ 空间中MICP加固石英砂^[23]和MICP加固钙质砂^[24]的临界状态线。可以看出，在 $e$ – $\ln p'$ 空间中MICP加固石英砂和MICP加固钙质砂的临界状态线随加固程度增加向上移动，且近似平行。造成这种现象的原因可能是由于剪切过程中试样达到临界状态时砂颗粒间仍然存在一定的胶结作用碳酸钙，同时附着在砂颗粒表面的碳酸钙也会产生影响。从可以看出， $p'$ – $q$ 空间中MICP加固石英砂的临界状态线随加固程度增加，斜率逐渐增大，表明被加固砂土的临界状态摩擦角增加，主要是由于加固程度较高时有更多的碳酸钙附着在砂颗粒表面，导致其表面粗糙度增加。同时，加固程度的提高也会使剪切过程中胶结砂颗粒的碳酸钙破碎之后转化为增加砂颗粒表面粗糙度的碳酸钙更多。从图 3(d)可以看出，MICP加固钙质砂的临界状态线与未加固钙质砂的临界状态线的斜率相差不大，表明MICP加固钙质砂的临界状态摩擦角增加程度有限，主要是因为钙质砂的颗粒表面存在内孔隙，生成的碳酸钙将这些内孔隙填充。胶结作用破坏后仍有部分碳酸钙附着在砂颗粒表面，但是其对临界状态摩擦角的贡献并不如MICP加固石英砂明显。

图 3 MICP加固石英砂和MICP加固钙质砂的临界状态线

Figure 3. Critical state lines of MICP-treated silica sand and MICP-treated calcareous sand

下载: 全尺寸图片幻灯片

2. MICP加固砂土本构模型

2.1 屈服面方程

Yao等^[25]基于UH模型中的椭圆屈服面，引入临界状态参数 $\chi$ ，提出了适用于砂土的屈服面形式，如。在 $p'$ – $q$ 平面内屈服面方程表示为

$f = \frac{{({\text{1 + }}\chi ){q^2}}}{{M_{{\text{cs}}}^2{{p'}^2} - \chi {q^2}}} - \frac{{{p_0}}}{{p'}} + 1 = 0 。$

(3)

图 4 砂土与MICP加固砂土的屈服面

Figure 4. Yield surfaces of sands and MICP-treated sands

下载: 全尺寸图片幻灯片

式中 $p'$ 为平均有效应力； $q$ 为剪应力； ${p_0}$ 为屈服面尺寸； ${M_{{\text{cs}}}}$ 为临界状态应力比； $\chi$ 为临界状态参数，可以控制屈服面形状。当 $\chi {\text{ = 0}}$ 时，屈服面形状退化为修正剑桥模型的椭圆，当 ${\text{0 < }}\chi {\text{ < 1}}$ 时，屈服面形状为水滴型，其对屈服面形状的影响可参考Yao等^[25]。

对于MICP加固砂土，其屈服面由于胶结作用相较于未加固砂土的屈服面向 $p'$ 轴的负半轴扩大，表示胶结作用使土体具有一定程度的抗拉强度。MICP加固砂土的屈服面形状如图 4所示。

对式(3)修正之后可得MICP加固砂土的屈服面方程为

$F = \frac{{({\text{1 + }}\chi ){q^2}}}{{M_{{cs} }^2{{(p' + {p_{\text{t}}})}^2} - \chi {q^2}}} - \frac{{{p_0} + {p_{\text{t}}}}}{{p' + {p_{\text{t}}}}} + 1 = 0 \text{，}$

(4)

式中， ${p_{\text{t}}}$ 为胶结作用引起的抗拉强度。

Baudet等^[26]在研究结构性土的本构模型时，认为结构性的丧失同时受塑性体积应变和塑性剪应变的影响，但是两者的影响比例目前还不能通过试验得到，可假设两者的影响相同。Chen等^[27]通过定义塑性损伤应变来表示塑性体积应变和塑性剪应变对胶结作用退化的综合影响：

$\varepsilon _{\text{d}}^{\text{p}} = \sqrt {{{(\varepsilon _{\text{v}}^{\text{p}})}^2} + {{(\varepsilon _{\text{q}}^{\text{p}})}^2}} 。$

(5)

本文采用式(5)的形式并认为抗拉强度 ${p_{\text{t}}}$ 随塑性损伤应变的累计逐渐减小，可以表示为

${p_{\text{t}}} = {p_{{\text{t0}}}}\exp \left( { - \xi \frac{{{\text{1 + }}{e_0}}}{{\lambda - \kappa }}\varepsilon _{\text{d}}^{\text{p}}} \right) \text{，}$

(6)

式中， ${p_{{\text{t0}}}}$ 为初始抗拉强度，可以由 $p'$ – $q$ 空间中的破坏包络线与 $p'$ 轴截距得到， $\xi$ 为胶结退化速率参数， $\lambda$ 和 $\kappa$ 分别为 $e$ – $\ln p'$ 空间中的临界状态线斜率和回弹线斜率。

Chen等^[27]将剪切过程中胶结退化速率与塑性损伤应变建立指数函数关系，表征胶结退化速率随塑性损伤应变的累计逐渐增加。对MICP加固砂土的试验结果表明，胶结作用退化速率与围压也有关。胶结退化速率 $\xi$ 需要综合考虑塑性应变与围压的影响，采用如下表达式：

$\xi = {\xi _0}\exp (\alpha {p_{{\text{ic}}}} + \beta \varepsilon _{\text{d}}^{\text{p}})\text{，}$

(7)

式中， ${\xi _0}$ 为初始胶结退化速率， $\alpha$ 和 $\beta$ 为材料参数， ${p_{{\text{ic}}}}$ 为固结围压。

加载方向的单位向量 $n = {\left[ {{n_{\text{p}}}, {n_{\text{q}}}} \right]^{\text{T}}}$ 的两个分量为

${n_{\text{p}}} = {\text{(}}\partial F/\partial p'{\text{)}}/{L_{\text{F}}} \text{，} {n_{\text{q}}} = {\text{(}}\partial F/\partial q{\text{)}}/{L_{\text{F}}} \text{，}$

(8)

式中， ${L_{\text{F}}}$ 为加载方向向量的模：

${L_{\text{F}}} = \sqrt {{{(\partial F/\partial p')}^2} + {{(\partial F/\partial q)}^2}} 。$

(9)

$\partial F/\partial p'$ 和 $\partial F/\partial q$ 分别为屈服面方程对 $p'$ 和 $q$ 的偏导：

$\frac{{\partial F}}{{\partial p'}} = - \frac{{2M_{{\rm{cs}}}^2(1 + \chi )(p' + {p_{\text{t}}}){q^2}}}{{{{\left[ {M_{{\rm{cs}}}^2{{(p' + {p_t})}^2} - \chi {q^2}} \right]}^2}}} + \frac{{{p_0} + {p_{\text{t}}}}}{{p' + {p_{\text{t}}}}} \text{，}$

(10)

$\frac{{\partial F}}{{\partial q}} = \frac{{2M_{{\rm{cs}}}^2(1 + \chi ){{(p' + {p_{\text{t}}})}^2}q}}{{{{\left[ {M_{{\rm{cs}}}^2{{(p' + {p_{\text{t}}})}^2} - \chi {q^2}} \right]}^2}}} 。$

(11)

2.2 塑性势函数与流动法则

砂土等散粒状土体的力学行为与土体所处的状态有关^[28]，胶结砂土的本构模型研究中同样也可以引入状态参数^[29-31]。Been等^[32]所提出的状态参数应用较为广泛，其可以表示为当前孔隙比 $e$ 与相同平均有效应力下对应的临界状态空隙比 ${e_{{\text{cs}}}}$ 之间的差值：

$\psi = e - {e_{{\text{cs}}}} \text{，}$

(12)

式中，临界状态线CSL上的孔隙比 ${e_{{\text{cs}}}}$ 可表示为

${e_{{\text{cs}}}} = {e_{{\text{cs0}}}} - \lambda \ln p' \text{，}$

(13)

式中， ${e_{{\text{cs0}}}}$ 为 $p' = 1$ kPa时的临界状态孔隙比。

采用非关联流动法则，塑性势函数为

$g = \frac{{({\text{1 + }}\chi ){q^2}}}{{{M_{\text{d}}}^{\text{2}}{{(p' + {p_{\text{t}}})}^2} - \chi {q^2}}} - \frac{{{p_0} + {p_{\text{t}}}}}{{p' + {p_{\text{t}}}}} + 1 = 0 \text{，}$

(14)

式中， ${M_{\text{d}}}$ 为剪胀应力比，其与状态参数有关，采用Li等^[33]提出的指数函数表达式：

${M_{\text{d}}} = {M_{{\text{cs}}}}\exp ({k_{\text{d}}}\psi ) \text{，}$

(15)

式中， ${k_{\text{d}}}$ 为剪胀相关的材料常数。

塑性流动方向单位向量 $m = {[{m_{\text{p}}}, {m_{\text{q}}}]^{\text{T}}}$ 的两个分量可分别表示为

${m_{\text{p}}} = (\partial g/\partial p')/{L_{\text{g}}} \text{，} {m_{\text{q}}} = (\partial g/\partial q)/{L_{\text{g}}} \text{，}$

(16)

式中， ${L_{\text{g}}}$ 为加载方向向量的模：

${L_{\text{g}}} = \sqrt {{{(\partial g/\partial p')}^2} + {{(\partial g/\partial q)}^2}} 。$

(17)

$\partial g/\partial p'$ 和 $\partial g/\partial q$ 分别为塑性势函数对 $p'$ 和 $q$ 的偏导：

$\frac{{\partial g}}{{\partial p'}} = - \frac{{2{M_{\text{d}}}^2(1 + \chi )(p' + {p_t}){q^2}}}{{{{\left[ {{M_{\text{d}}}^{\text{2}}{{(p' + {p_t})}^2} - \chi {q^2}} \right]}^2}}} + \frac{{{p_0} + {p_t}}}{{p' + {p_t}}} \text{，}$

(18)

$\frac{{\partial g}}{{\partial q}} = \frac{{2{M_{\text{d}}}^{\text{2}}(1 + \chi ){{(p' + {p_t})}^2}q}}{{{{\left[ {{M_{\text{d}}}^{\text{2}}{{(p' + {p_{\text{t}}})}^2} - \chi {q^2}} \right]}^2}}} 。$

(19)

2.3 弹塑性增量关系

弹性体积应变分量和剪应变分量为

${\text{d}}\varepsilon _{\text{v}}^{\text{e}} = {\text{d}}p'/K \text{，} {\text{d}}\varepsilon _{\text{q}}^{\text{e}} = {\text{d}}q/(3G) \text{，}$

(20)

式中， $K$ 为体积模量， $G$ 为剪切模量。表达式分别如下：

$K = (1 + {e_0})p/\kappa \text{，} G = 1.5K(1 - 2\nu )/(1 + \nu ) \text{，}$

(21)

式中， ${e_0}$ 为固结完成时的初始孔隙比， $\nu$ 为泊松比。

塑性体积应变分量和剪应变分量为

${\text{d}}\varepsilon _{\text{v}}^{\text{p}} = \left\langle {({n_{\text{p}}}{\text{d}}p' + {n_{\text{q}}}{\text{d}}q)/H} \right\rangle {m_{\text{p}}} \text{，}$

(22)

${\text{d}}\varepsilon _q^p = \left\langle {({n_{\text{p}}}{\text{d}}p' + {n_{\text{q}}}{\text{d}}q)/H} \right\rangle {m_{\text{q}}} \text{，}$

(23)

式中， $H$ 为塑性模量， $\left\langle {} \right\rangle$ 为Macaulay括号：当 $a \leqslant 0$ 时， $\left\langle a \right\rangle = 0$ ；当 $a > 0$ 时， $\left\langle a \right\rangle = a$ 。

2.4 塑性模量

塑性模量一般可以通过屈服面方程的一致性条件求得，但是为了描述土体塑性变形的状态相关性，需要将状态参数引入塑性模量中。本文塑性模量采用如下形式：

$H = {h_0}G{M_{{\text{cs}}}}({M_{\text{b}}} - {\eta ^*})/({M_{\text{b}}}{\eta ^*}) \text{，}$

(24)

式中， ${h_0}$ 为塑性模量参数， ${M_{\text{b}}}$ 为峰值应力比， ${\eta ^*}$ 为修正应力比。

峰值应力比 ${M_{\text{b}}}$ 表达式为^[33]

${M_{\text{b}}} = {M_{{\text{cs}}}}\exp ( - {k_{\text{b}}}\psi ) \text{，}$

(25)

式中， ${k_{\text{b}}}$ 为峰值应力比相关的材料常数。

修正应力比 ${\eta ^*}$ 可以表示为

${\eta ^*} = q/(p' + {p_{\text{t}}}) 。$

(26)

3. 模型预测

3.1 模型参数分析

引入状态参数的MICP胶结砂土边界面本构模型共有13个参数，其中 $\kappa$ 和 $\mu$ 为弹性相关的参数， $\lambda$ ， $\chi$ ， ${M_{{\text{cs}}}}$ 和 ${e_{{\text{cs0}}}}$ 为临界状态参数， ${p_{{\text{t0}}}}$ ， ${\xi _0}$ ， $\alpha$ 和 $\beta$ 为胶结作用及胶结退化有关的参数， ${k_{\text{d}}}$ 为剪胀相关的参数， ${k_{\text{b}}}$ 为峰值状态参数， ${h_0}$ 为塑性模量参数。

参数 $\lambda$ ， $\kappa$ ， $\mu$ ， ${M_{{\text{cs}}}}$ 为修正剑桥模型中的参数，可通过三轴试验和等向压缩试验获得。 ${e_{{\text{cs0}}}}$ 可以通过不同围压三轴试验的临界状态孔隙比在 $e$ – $\ln p'$ 空间中拟合得到^[27]。参数 $\chi$ 通过Yao等^[25]介绍的方法确定，但是由于缺乏MICP加固砂土正常固结线NCL的试验数据，该参数通过拟合试验数据的应力–应变关系得到。 ${p_{{\text{t0}}}}$ 的确定方法前文已经叙述。 ${\xi _0}$ ， $\alpha$ 和 $\beta$ 可以通过拟合三轴排水剪切试验的应力–应变关系得到。 ${k_{\text{d}}}$ 和 ${k_{\text{b}}}$ 的确定方法可以参考Li等^[33]、Xiao等^[34]和孙增春等^[35]。塑性模量参数 ${h_0}$ 可以通过拟合三轴剪切试验的应力–应变曲线得到^[34]。采用中的基准分析参数，分析临界状态参数 $\chi$ ，胶结作用退化速率参数 ${\xi _0}$ ， $\alpha$ 和 $\beta$ 对应力–应变关系，剪胀规律和胶结作用演化的影响。

表 1 模型分析参数

Table 1. Model analysis parameters

初始状态	弹性参数	临界状态参数	胶结作用参数	其他参数
${e_0} = 0.8$ ${p_{{\text{ic}}}} = 100$ kPa	$\mu$ =0.25 $\kappa = 0.005$	$\lambda = 0.08$	${p_{t0}} = 80$ kPa	${k_{\text{d}}} = 1.5$ ${k_{\text{b}}} = 1.0$ ${h_0} = 2.0$
		${M_{{\text{cs}}}} = 1.5$	${\xi _0} = 0.2$
		${e_{{\text{cs0}}}} = 1.2$	$\alpha = 0.01$
		$\chi = 0.1$	$\beta = 1.0$

下载: 导出CSV

| 显示表格

为参数 $\chi$ ， ${\xi _0}$ ， $\alpha$ ， $\beta$ 的变化对应力–应变关系、体变规律和胶结退化的影响。从可知，参数 $\chi$ 对体积应变的变化影响明显，随着 $\chi$ 的增大，体变由剪缩过渡到剪胀。同时，随着 $\chi$ 的增大，胶结退化速率增加，应变软化更加明显，但是峰值偏应力不变，说明 $\chi$ 的变化主要影响峰值强度后的胶结破坏速率，且 $\chi$ 增大会使到达峰值强度的应变减小。可以看出，参数 ${\xi _0}$ 和 $\alpha$ 对模拟结果的影响相似，即随着参数取值增加，胶结退化加快，且峰值强度降低，应变软化和剪胀更加明显，说明上述两个参数影响不仅影响峰值强度后的胶结破坏速率，同时影响峰值强度前的胶结破坏速率，即随着参数取值增大，在达到峰值强度前已有大量胶结作用破坏，使被加固土体整体强度降低。从可以看出，随着参数 $\beta$ 的增大，胶结作用破坏速率增加，同时应变软化和剪胀更加明显，但是峰值强度和剪缩部分的体变几乎不受影响。从式(5)可以看出，参数 $\beta$ 主要控制塑性损伤应变对胶结作用退化速率的影响，随着塑性应变的累计，参数 $\beta$ 对退化速率的影响增大。

图 5 参数对模型预测的影响

Figure 5. Influences of parameters on model performance

下载: 全尺寸图片幻灯片

3.2 模型描述与试验对比

为了验证模型的适用性，选取MICP加固Ottawa 20/30砂、石英砂和钙质砂的三轴排水试验结果与模型计算结果进行对比分析。其中，Lin等^[22]分别对未加固和MICP加固的Ottawa 20/30砂开展了三轴排水剪切试验，试样的初始孔隙比为0.65，初始有效固结围压分别为25，50，100 kPa。Xiao等^[23]分别对未加固和3种MICP加固程度(B_ca=1.8%，B_ca =3.5%和B_ca=5.1%)的石英砂开展了4种有效固结围压(20，50，100，200 kPa)三轴排水剪切试验，试样的初始孔隙比为0.65。Cui等^[24]对不同MICP加固程度的钙质砂进行三轴排水剪切试验，固结围压分别为100，200，300 kPa，试样的初始孔隙比为0.92。模型参数均通过试验数据进行校准，不同材料的模型计算参数取值见表 2。

表 2 模型计算参数

Table 2. Model parameters

参数	Ottawa 20/30砂	MICP 加固Ottawa 20/30砂	MICP 加固石英砂B_ca=1.8%	MICP 加固石英砂B_ca=3.5%	MICP加固钙质砂
$\mu$	0.32	0.30	0.30	0.30	0.25
$\lambda$	0.009	0.001	0.049	0.051	0.047
$\kappa$	0.001	0.015	0.005	0.005	0.002
${M_{{\text{cs}}}}$	1.29	1.43	1.23	1.42	1.43
${e_{{\text{cs0}}}}$	0.714	0.736	1.024	1.054	1.349
$\chi$	0.02	0.1	0.1	0.1	0.1
${k_{\text{d}}}$	2.5	1.5	1.5	1.2	1.5
${k_{\text{b}}}$	1.0	1.5	2.0	0.8	0.1
${h_0}$	0.8	5.0	1.5	2.5	3.5
${p_{{\text{t0}}}}$ /kPa	—	92.5	32.21	36.29	434.5
${\xi _0}$	—	0.28	0.21	0.15	0.27
$\alpha$	—	0.013	-0.016	0.01	-0.005
$\beta$	—	5.0	1.0	2.0	5.0

下载: 导出CSV

| 显示表格

图 6为Ottawa 20/30砂和MICP加固Ottawa 20/30砂^[22]在不同围压下的试验结果与本构模型模拟结果对比。可以看出，模型可以较好地模拟Ottawa 20/30砂在各围压下的应力–应变发展规律和剪胀现象。对于MICP加固Ottawa 20/30砂，本文建立的本构模型计算得到各围压下的应力–应变曲线的峰值应力点对应的轴向应变大于试验结果，模拟的剪胀整体大于试验结果。但是，本文模型可以较好地模拟出MICP加固Ottawa 20/30砂在达到峰值强度后，由于塑性应变的累计，碳酸钙对土体的胶结大量破坏而出现强度陡降的脆性破坏。

图 6 Ottawa 20/30砂和MICP加固Ottawa 20/30砂三轴排水试验结果与模型预测对比

Figure 6. Comparison between test results and model predictions on Ottawa 20/30 sand and MICP-treated Ottawa 20/30 sand under drained triaxial compression

下载: 全尺寸图片幻灯片

图 7为两种加固程度(B_ca=1.8%，B_ca =3.5%)的MICP加固石英砂^[23]在不同初始有效围压下的试验结果与数值模拟结果对比，从图中可以看出模型能够较好地反映MICP加固石英砂的软化特性和体积变化规律。图 8给出了本文模型计算结果与Cui等^[24]对MICP加固钙质砂(CCC=25.5%)的三轴排水剪切试验结果的对比，可以看出本文所建立的本构模型可以较好地模拟出MICP加固钙质砂的应变软化特性，且其胶结退化速率随围压改变而不同，具体为低围压时胶结退化速率较快，软化更加明显；高围压时胶结退化速率较慢。同时，上述结果也验证了本文将胶结退化速率与围压建立关系的正确性。对MICP加固钙质砂体变的变化规律的模拟结果较差，表现为低围压时剪胀较小，且剪胀增长速率较慢，高围压时剪胀过大。

图 7 B_ca=1.8%和B_ca=3.5% MICP加固石英砂三轴排水试验结果与模型预测对比

Figure 7. Comparison between test results and model predictions on MICP-treated silica sand for B_ca=1.8% and B_ca=3.5% under drained triaxial compression

下载: 全尺寸图片幻灯片

图 8 CCC =25.5% MICP加固钙质砂三轴排水试验结果与模型预测对比

Figure 8. Comparison between test results and model predictions on MICP-treated calcareous sand for CCC =25.5% under drained triaxial compression

下载: 全尺寸图片幻灯片

4. 结论

本文针对MICP加固砂土的强度和变形特征，在分析加固和破坏机理的基础上，基于临界状态土力学理论框架，将胶结退化速率与围压和塑性应变建立关系，采用非关联流动法则，引入状态参数和剪胀应力比及峰值应力比，建立了MICP加固砂土状态相关的弹塑性本构模型。将模型计算结果分别与MICP加固Ottawa 20/30砂、MICP加固石英砂和MICP加固钙质砂的三轴试验结果进行对比，得到以下结论：

(1) 通过对MICP加固砂土的破坏包络线进行分析，发现石英砂和钙质砂经MICP加固后土体黏聚力都有所增加，但是MICP加固石英砂的峰值摩擦角增加，MICP加固钙质砂的峰值摩擦角减小。主要是两种砂土颗粒本身的形状、表面粗糙度及钙质砂存在的内孔隙使碳酸钙沉积方式及对力学性能的贡献存在差异导致。

(2) 对MICP加固石英砂和MICP加固钙质砂的临界状态线进行分析，发现临界状态线在 $e$ – $\ln p'$ 空间中随加固程度增加而向上移动，可能是由于被加固土体中的胶结作用未完全破坏或碳酸钙附着在砂颗粒表面的影响。在 $p'$ – $q$ 空间中随加固程度增加MICP加固石英砂的临界状态线的斜率逐渐增大，但是MICP加固钙质砂的斜率增加不明显，主要是由于钙质砂的颗粒表面存在内孔隙。

(3) 通过对不同MICP加固程度的3种砂土的三轴排水剪切试验的模拟，表明本文建立的本构模型能够较好地模拟MICP加固砂土随胶结退化出现的应变软化行为及体积变化规律。

图 1 浮置板轨道结构

Figure 1. Floating-slab track structure

下载: 全尺寸图片幻灯片

图 2 2.5维有限元模型

Figure 2. 2.5D FEM model

下载: 全尺寸图片幻灯片

图 3 本文研究结果与文献[18]对比

Figure 3. Comparison between proposed results and Reference [18]

下载: 全尺寸图片幻灯片

图 4 不同车速下填充沟深度对隔振效果的影响

Figure 4. Influences of h^* on A_r at different train speeds

下载: 全尺寸图片幻灯片

图 5 不同车速下填充沟深度与平均幅值衰减比的关系

Figure 5. Relationship between h^* and $\bar{A_{r}}$ at different train speeds

下载: 全尺寸图片幻灯片

图 6 不同车速下隧道埋深对隔振效果的影响

Figure 6. Influences of h₀ on A_r at different speeds

下载: 全尺寸图片幻灯片

图 7 不同车速下填充沟宽度对隔振效果的影响

Figure 7. Influences of b^* on A_r at different train speeds

下载: 全尺寸图片幻灯片

图 8 不同车速下填充沟与轨道中心水平间距对隔振效果的影响

Figure 8. Influences of d^* on A_r at different train speeds

下载: 全尺寸图片幻灯片

表 1 弹性土和衬砌结构参数

Table 1 Parameters of elastic ground and tunnel structure

介质	密度ρ/(kg·m^-3)	剪切模量μ/Pa	泊松比 $ν$	阻尼系数β
弹性土	1.9×10³	2.0×10⁷	0.30	0.01
衬砌结构	2.5×10³	1.4×10¹⁰	0.25	0.02

下载: 导出CSV

表 2 地铁列车荷载参数

Table 2 Parameters of train load

机车轴重P_E/kN	车厢轴重P_C/kN	一、二组轴载间距w_a/m	二、三组轴载间距w_b/m	机车长度L_E /m	车厢长度L_C/m	列车数量N_T	观察点L_D/m
162	120	2.9	14.8	22.2	24.4	4	0

下载: 导出CSV

表 3 浮置板轨道结构参数

Table 3 Parameters of floating-slab tracks

介质	抗弯刚度EI/(Pa·m⁴)	单位长度质量m/(kg·m^-1)	垫圈弹簧常数k/Pa	垫圈阻尼常数c/(N·S·m^-2)
钢轨	1.00×10⁷	100	4×10⁷	6.30×10³
浮置板	1.43×10⁹	3.5×10³	5×10⁷	4.18×10⁴

下载: 导出CSV

表 4 模型尺寸参数

Table 4 Size parameters of model

隧道中心到地表距离h₀/m	填充沟沟深h/m	填充沟宽度b/m	填充沟到轨道中心水平距离d/m
16	8	2	6

下载: 导出CSV

参考文献(21)

[1]	WOODS R D. Screening of surface waves in soils[J]. Journal of the Soil Mechanics and Foundations Division, 1968, 94(4): 951-979. doi: 10.1061/JSFEAQ.0001180
[2]	AHMAD S, ALHUSSAINI T M. Simplified design for vibration screening by open and in-filled trenches[J]. Journal of Geotechnical Engineering, 1991, 117(1): 67-88. doi: 10.1061/(ASCE)0733-9410(1991)117:1(67)
[3]	高广运, 冯世进, 李伟, 等. 三维层状地基竖向激振波阻板主动隔振分析[J]. 岩土工程学报, 2007, 29(4): 471-476. doi: 10.3321/j.issn:1000-4548.2007.04.001 GAO Guang-yun, FENG Shi-jin, LI Wei, et al. 3D analysis of active vibration isolation with wave impeding block in layered ground under vertical loading[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(4): 471-476. (in Chinese) doi: 10.3321/j.issn:1000-4548.2007.04.001
[4]	ADAM M, VON ESTORFF O. Reduction of train-induced building vibrations by using open and filled trenches[J]. Computers and Structures, 2005, 83(1): 11-24. doi: 10.1016/j.compstruc.2004.08.010
[5]	BESKOS D E, DASGUPTA B, VARDOULAKIS I G. Vibration isolation using open or filled trenches Part l: 2D homogeneous soil[J]. Computational Mechanics, 1986(1): 43-63.
[6]	YANG Y B, HUNG H H. A parametric study of wave barriers for reduction of train-induced vibrations[J]. International Journal for Numerical Methods in Engineering, 1997, 40: 3729-3747. doi: 10.1002/(SICI)1097-0207(19971030)40:20<3729::AID-NME236>3.0.CO;2-8
[7]	邓亚虹, 夏唐代, 陈敬虞. 车辆动荷载作用下隔振沟响应增强区数值分析[J]. 岩土工程学报, 2006, 28(12): 2121-2127. doi: 10.3321/j.issn:1000-4548.2006.12.014 DENG Ya-hong, XIA Tang-dai, CHEN Jing-yu. Numerical analysis of response magnified area of vibration isolation trenches subjected to vehicle dynamic loads[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(12): 2121-2127. (in Chinese) doi: 10.3321/j.issn:1000-4548.2006.12.014
[8]	BESKOS D E, DASGUPTA B, VARDOULAKIS I G. Vibration isolation using open or filled trenches Part 2: 3D homogeneous soil[J]. Computational Mechanics, 1990(7): 137-148.
[9]	YANG Y B, HUNG H H. A 2.5D finite/infinite element approach for modelling visco-elastic bodies subjected to moving loads[J]. International Journal for Numerical Methods in Engineering, 2001, 51(11): 1317-1336. doi: 10.1002/nme.208
[10]	HUNG H H, YANG Y B, CHANG D W. Wave barriers for reduction of train-induced vibrations in soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2004, 130(12): 1283-1291. doi: 10.1061/(ASCE)1090-0241(2004)130:12(1283)
[11]	高广运, 何俊锋, 李宁, 等. 饱和地基上列车运行引起的地面振动隔振分析[J]. 岩土力学, 2011, 32(7): 2191-2198. doi: 10.3969/j.issn.1000-7598.2011.07.045 GAO Guang-yun, HE Jun-feng, LI Ning, et al. Analysis of isolating ground vibration induced by trains running on saturated ground[J]. Rock and Soil Mechanics, 2011, 32(7): 2191-2198. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.07.045
[12]	陈功奇, 高广运. 层状地基中填充沟对不平顺列车动荷载的隔振效果研究[J]. 岩石力学与工程学报, 2014, 33(1): 144-153. doi: 10.13722/j.cnki.jrme.2014.01.021 CHEN Gong-qi, GAO Guang-yun. Vibration screening effect of in-filled trenches on train dynamic loads of geometric irregular track in layered grounds[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(1): 144-153. (in Chinese) doi: 10.13722/j.cnki.jrme.2014.01.021
[13]	VAN HOORICKX C, SCHEVENELS M, LOMBAERT G. Double wall barriers for the reduction of ground vibration transmission[J]. Soil Dynamics and Earthquake Engineering, 2017, 97: 1-13. doi: 10.1016/j.soildyn.2017.02.006
[14]	巴振宁, 梁建文, 王靖雅. 空沟对层状饱和地基中列车移动荷载的隔振性能研究[J]. 岩土工程学报, 2017, 39(5): 848-858. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201705012.htm BA Zhen-ning, LIANG Jian-wen, WANG Jing-ya. Isolation effect of an open trench against train-induced vibrations in a saturated layered ground[J]. Chinese Journal of Geotechnical Engineering,2017, 39(5): 848-858. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201705012.htm
[15]	BALENDRA T, CHUA K H, LO K W, et al. steady-state vibration of subway-soil-building system[J]. Journal of Engineering Mechanics, 1989, 115(1): 145-162. doi: 10.1061/(ASCE)0733-9399(1989)115:1(145)
[16]	ZHANG Sheng-long, WANG Wen-bin, WU Zong-zhen. Effect of periodic pile row in subway vibration isolation[J]. Procedia Engineering, 2017, 199: 302-309. doi: 10.1016/j.proeng.2017.09.043
[17]	DEEKS A J. Axisymmetric time-domain transmitting boundaries[J]. Journal of Engineering Mechanics, 1994, 120(1): 25-42. doi: 10.1061/(ASCE)0733-9399(1994)120:1(25)
[18]	EASON G. The stresses produced in a semi-infinite solid by a moving surface force[J]. International Journal of Engineering Sciences, 1965, 2(6): 581-609. doi: 10.1016/0020-7225(65)90038-8
[19]	HUNG H H, CHEN G H, YANG G H. Effect of railway roughness on soil vibrations due to moving trains by 2.5D finite/infinite element approach[J]. Engineering Structures, 2013, 57: 254-266. doi: 10.1016/j.engstruct.2013.09.031
[20]	HUSSEIN M F M, HUNT H E M. Modelling of floating-slab tracks with continuous slabs under oscillating moving loads[J]. Journal of Sound and Vibration, 2006, 297(1/2): 37-54.
[21]	蔡袁强, 王玉, 曹志刚. 移动荷载作用下饱和地基上板式轨道动力分析[J]. 振动工程学报, 2011, 24(1): 48-54. doi: 10.3969/j.issn.1004-4523.2011.01.009 CAI Yuan-qiang, WANG Yu, CAO Zhi-gang. Dynamic response of slab track and saturated soil subjected to a moving train load[J]. Journal of Vibration Engineering, 2011, 24(1): 48-54. (in Chinese) doi: 10.3969/j.issn.1004-4523.2011.01.009

施引文献(23)

期刊类型引用(8)

1.	刘雪莹，俞缙，周先齐，姚玮，任崇鸿，蔡燕燕. 一种新的岩石各向异性蠕变-渗透率模型. 岩土工程学报. 2025(02): 275-283 . 本站查看
2.	郭建春，张涛，武玺，赵志红，李宗源，曾杰. 煤层气压裂水平井生产动态分析及其渗透率协同演化机制. 煤炭学报. 2025(01): 516-531 . 百度学术
3.	张磊丽，蔡婷婷，石磊，姜玉龙，夏瑾. 不同温度-围压-气体压力下煤体蠕变-渗流演化规律. 煤矿安全. 2024(03): 36-45 . 百度学术
4.	赵一博. 计算多孔介质流线迂曲度及其分形维数的新几何模型. 材料导报. 2024(S2): 282-286 . 百度学术
5.	李钊，倪冠华，杨威，温永瓒，王刚，刘震，王振洋. 不同冲击倾向性煤体孔隙结构及瓦斯吸附特征研究. 煤炭科学技术. 2024(S2): 56-70 . 百度学术
6.	张学博，王豪，杨明，王攀，韩林秀. 抽采钻孔失稳坍塌对瓦斯抽采的影响机制研究及应用. 煤炭学报. 2023(08): 3102-3115 . 百度学术
7.	王科文. 瓦斯抽采钻孔新型封孔工艺试验研究. 山东煤炭科技. 2023(09): 128-130 . 百度学术
8.	许江，饶豪魁，彭守建，陈奕安，陈灿灿，马天宇. 三轴应力条件下凝灰岩广义应力松弛流变模型的构建. 采矿与安全工程学报. 2023(06): 1264-1272 . 百度学术