土压平衡盾构渣土浆液配比优化研究

张箭; 金俊杰; 丰土根; 刘江涛

doi:10.11779/CJGE20220024

土压平衡盾构渣土浆液配比优化研究 English Version

1.
河海大学岩土力学与堤坝工程教育部重点实验室, 江苏南京 210098
2.
国网江苏省电力有限公司常州供电分公司, 江苏常州 213000
3.
中铁(上海)投资集团有限公司, 上海 200126

基金项目:

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

中央高校基本科研业务费专项资金项目 B220202016

详细信息

作者简介:
张箭(1989—)，男，博士，副教授，主要从事隧道与地下工程方面的研究工作。E-mail: zhangj0507@163.com

通讯作者:
丰土根, E-mail: fengtugen@hhu.edu.cn

中图分类号: TU441
计量
- 文章访问数: 241
- HTML全文浏览量: 107
- PDF下载量: 81
出版历程
- 收稿日期: 2022-01-05
- 网络出版日期: 2023-04-16
- 刊出日期: 2023-03-31

Optimization of mixture ratio of muck grout by earth pressure balance shield machine

1.
Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering (Hohai University), Nanjing 210098, China
2.
Changzhou Power Supply Branch of State Grid Jiangsu Electric Power Co., Ltd., Changzhou 213000, China
3.
China Railway (Shanghai) Investment Group Co., Ltd., Shanghai 200126, China

摘要

摘要: 粉质黏土是城市浅覆盾构隧道工程中常见渣土，如何实现渣土的就地资源再利用是盾构领域亟待解决的问题。以南京地铁7号线某粉质黏土地层盾构掘进段为工程背景，采用X射线衍射技术及室内试验等手段，结合SPSS多元回归分析及Matlab多目标规划工具，开展符合现场要求的渣土浆液配比优化研究，探讨粉质黏土替代壁后注浆浆液膨润土的可行性，结合现场实际，研究浆液性能及注浆参数等对渣土浆液渗透阶段扩散规律的影响。研究结果表明：通过配比优化，粉质黏土可替代膨润土制备满足现场施工要求的渣土注浆浆液；性能最优配比渣土浆液在7 d和28 d强度、泌水率和28 d结石体收缩率方面有优势，具备较好的胶结强度、结石效果和稳定性，可节省10.0%的成本；成本最优配比浆液比原配比浆液性能更优，且可节省20.2%的成本。相较于原配比浆液，采用性能最优配比浆液时要适当的减小注入率，采用成本最优配比浆液时要适当的增加注入率。
- 盾构隧道 /
- 同步注浆 /
- 粉质黏土 /
- 多目标规划 /
- 扩散机理
Abstract: The silty clay is the common muck in urban shallow tunnels, and how to realize the reuse of on-site resources of muck is an urgent problem to be solved in the field of shield tunneling. A shield tunneling section of Nanjing Metro Line 7 in silty clay formation is taken as the engineering background. Based on the SPSS multiple regression analysis and the Matlab multi-objective planning tools, the X-ray diffraction technology and the laboratory tests are used to carry out the optimization of the mixture ratio of muck grout. The feasibility of the silty clay as a substitute for bentonite of backfill grout is discussed. Based on the actual situation, the influences of grout performance and grouting parameters on the diffusion laws of the grout at the seepage stage are studied. The results show that the silty clay can be used to replace the bentonite to prepare grout that meets the requirements of on-site construction through the optimization of the mixture ratio. The grout with the best performance has advantages in 7 d & 28 d-strength, bleeding rate and 28 d shrinkage rate of bonding body. It has good bonding strength, effect and stability, and 10.0% of the costs can be saved. The grout with the cost-optimized mixture ratio has better performance than the grout with the original mixture ratio, and 20.2% of the costs can be saved. Compared with the original grout, the injection rate should be appropriately reduced when the grout with the best performance is used, and the injection rate should be increased when the grout with cost-optimized mixture ratio is used.
- shield tunnel /
- simultaneous grouting /
- silty clay /
- multi-objective planning /
- diffusion mechanism

HTML全文

0. 引言

顶管技术作为近年来非开挖施工的一项重要工艺，在现代地下空间修建的发展中，起着重要的作用。在地下隧道的建设中，顶管隧道装配式顶进的特性，使其具备广泛的应用前景^[1]。顶管隧道结构纵向十分脆弱，易受到如地震、地层变形等影响，而等效刚度是纵向变形计算的基础参数，且目前鲜有矩形顶管F承插型接头等效刚度的相关研究，因此对于等效刚度的研究是紧迫的。

顶管隧道属于预制装配式混凝土结构，管节之间存在接缝，导致纵向变形的计算不能直接等效为弹性地基梁模型，而应该综合考虑接头的构造与力学行为。朱合华等^[2-3]建立了梁-接头不连续模型，用接头单元来模拟管节间的不连续现象，通过数值模拟与室内接头试验，建立接头刚度模型及其取值方法，研究了顶管的管节接缝张开量问题。丁文其等^[4]提出管节-接头模型，在顶力的角度下，对弯曲性能进行探究。以上几位学者将隧道视作Euler-Bernoulli梁模型，主要对圆形顶管隧道接头张开量进行研究，随着隧道断面的不断增大，Euler-Bernoulli梁模型的局限性也逐渐凸显。

对于大断面隧道，将其视为Timoshenko梁模型能更好的反应隧道变形情况。张冬梅等^[5]提出了接头力学模型，并把接头部位的相互作用离散成由剪切弹簧与法向弹簧组合而成的系列弹簧。闫治国等^[6-8]采用管片接头原型荷载试验的方法，研究了纵缝转角刚度和环缝剪切刚度，同时进行两种类型高刚性盾构管片接头正弯矩试验，对不同类型铸铁件力学性能和高刚性接头破坏模式进行了分析。Salemi等^[9]对混凝土界面进行了大量直剪试验，得到了接触面抗剪能力随接触面轴向力的线性关系。作为装配式地下建筑，接头的刚度是不可忽略的重要参数，但却不利于纵向变形的计算，因此应将接头部位进行折减，从而得到考虑接头变形的纵向等效刚度。

国内外计算等效抗弯刚度主要基于等效连续化模型进行分析^[10-12]，发现其影响因素主要有横向刚度、环缝作用、螺栓预应力^[13-14]。并以此考虑了刚度有效率^[15-16]。以上几位学者提出了等效刚度的影响因素与有效率，但未给出具体的等效刚度的求解方式。

为方便设计施工，等效刚度的解析解尤为重要。部分学者^[17-18]基于截面特点，得到等效抗弯刚度模型与解析解。耿萍等^[19]综合考虑轴力与弯矩对盾构隧道纵向变形的影响，提出5种弯曲模式并建立相应的非线性等效抗弯刚度计算模型。肖时辉等^[20]考虑到管片纵向连接螺栓受剪切作用时也会受弯曲作用，基于Timoshenko理论推导出大直径盾构管片纵向连接抗剪刚度计算公式。

综上，现有的研究大多基于盾构接头，且以圆形截面为主，对于矩形顶管隧道的研究尚不充分。在纵向等效刚度的理论研究方面，大部分学者仅对纵向等效抗弯刚度进行研究，未考虑管节之间的剪切作用。为此，本文针对矩形顶管隧道F型承插接头变形展开研究，探索接头的等效刚度解析解。

1. 矩形顶管隧道纵向等效刚度解析解

根据引言所述内容，若要对隧道的纵向变形进行预测，首先要得到隧道的等效刚度，即求得等效刚度的解析解。本节将根据Timoshenko模型，结合矩形顶管的相关特征，推导等效刚度解析解。F承插型接头矩形顶管如图 1所示，顶管与地层的纵向示意图如图 2所示。

图 1 F承插型接头顶管管节图

Figure 1. Pipe joint of pipe jacking with F-type socket joint

下载: 全尺寸图片幻灯片

图 2 顶管-地层纵向示意图

Figure 2. Vertical diagram of pipe jacking-strata

下载: 全尺寸图片幻灯片

1.1 矩形顶管隧道纵向等效抗弯刚度

为求得矩形顶管隧道纵向等效抗弯刚度模型，需要进行以下基本假定：

（1）平截面假定，该模型剪切变形远小于弯曲变形，平截面假定近似可用。

（2）管节以中轴线为平衡线，在管节的拉伸区，钢套环受到均匀的拉应力作用，而在压缩区，只有管节混凝土受到压应力作用，且假定为理想弹性状态。

（3）定义钢套环沿截面均匀分布时的等效弹性模量为 ${E_{\text{s}}}$ ，计算方法如下所示：

。

${E_{\text{s}}} = \frac{{{E_{\text{r}}}{A_{\text{s}}}}}{{{A_{\text{c}}}}}。$

(1)

式中： ${E_{\text{r}}}$ 为钢套环的弹性模量； ${A_{\text{s}}}$ 为钢套环沿隧道环向的横截面面积； ${A_{\text{c}}}$ 为矩形顶管隧道混凝土管节的横截面面积。

基于弹性体转动力学平衡条件和变形协调条件，将相邻管节中线的距离作为计算单位 $l$ ，当管节受到弯矩 $M$ 的作用时，发生弯曲变形，相邻管节中线延伸相交，形成转角 $\theta$ ，管节弯曲变形如图 3所示。

图 3 管节单元受弯变形示意图

Figure 3. Bending deformations of tubular element

下载: 全尺寸图片幻灯片

在弹性条件下，矩形环状截面尺寸如图 4（a）所示，图中m和n分别是隧道管片壁厚截面短边中心线到x轴的距离和长边中心线到y轴的距离； ${t_1}$ 和 ${t_2}$ 分别是管节长短边的厚度（此处统一取为 $t$ ）； $h$ 为隧道矩形环状截面中性轴与矩形截面长轴的距离。取图 4（b）所示计算模型进行分析，其变形协调方程为

$\left.\begin{array}{l}{\delta }_{\text{s(max)}}=(n+h)\theta \text{ }\text{，}\\ {\delta }_{\text{c(max)}}=(n-h)\theta \text{ }。\end{array}\right\}$

(2)

图 4 矩形顶管力学模型简图

Figure 4. Mechanical model for rectangular pipe jacking

下载: 全尺寸图片幻灯片

由上述定义设 ${E_{\text{s}}}$ 为钢套环沿截面分布的等效弹性模量； ${E_{\text{c}}}$ 为管节混凝土弹性模量； ${A_{\text{c}}}$ 为矩形顶管隧道混凝土管节横截面面积；以相邻管节中心线内的长度 $l$ 为一个单元进行计算； ${l_{\text{s}}}$ 为钢套环沿隧道纵向分布的长度。 ${\varepsilon _{{\text{s(max)}}}}$ 与 ${\varepsilon _{{\text{c(max)}}}}$ 为管节截面最大拉应变与压应变。

由管节受力平衡条件得

$\frac{{{E_{\text{s}}}}}{{{l_{\text{s}}}}}\int_{ - n}^h {(h - y)} t{\text{d}}y + \frac{{n + h}}{{{l_{\text{s}}}}}{E_{\text{s}}}mt$

$= \frac{{{E_{\text{c}}}}}{l}\int_h^n {(y - h)} t{\text{d}}y + \frac{{n - h}}{l}{E_{\text{c}}}mt。$

(3)

同理，由弯矩平衡条件可得

$2\frac{{{E_{\text{s}}}}}{{{l_{\text{s}}}}}\theta \int_{ - n}^h {{{(h - y)}^2}} t{\text{d}}y + 2\theta \frac{{n + h}}{{{l_{\text{s}}}}}{E_{\text{s}}}mt(n - h) +$

$2\frac{{{E_{\text{c}}}}}{l}\theta \int_h^n {{{(y - h)}^2}} t{\text{d}}y + 2\theta \frac{{({n^2} - {h^2})}}{l}{E_{\text{c}}}mt{\text{ = }}M 。$

(4)

结合式（2）~（4）可得

$\theta = \frac{M}{{f(h)}} 。$

(5)

式（5）中：

$\begin{gathered} f(h) = 2{E_{\text{s}}}\frac{t}{{{l_{\text{s}}}}}\left[ {\frac{1}{3}({h^3} + {n^3}) - h({h^2} - {n^2}) + {h^2}(h + n)} \right] + \hfill \\ 2\frac{{n + h}}{{{l_{\text{s}}}}}{E_{\text{s}}}mt(n - h) + 2\frac{{{n^2} - {h^2}}}{l}{E_{\text{c}}}mt + 2{E_{\text{c}}}\frac{t}{l} \cdot \hfill \\ \end{gathered}$

$\left[ {\frac{1}{3}({n^3} - {h^3}) - h({n^2} - {h^2}) + {h^2}(n - h)} \right] 。$

由等效连续梁转角公式：

$\theta = \frac{{Ml}}{{{{(EI)}_{{\text{eq}}}}}}。$

(6)

由于进行顶管隧道接头的特殊性，在受弯变形过程中，接头弯矩主要由钢套环承担。结合几何关系，可求得管节间张开角 $\theta$ 表达式为

$\theta = \frac{{M{l_{\text{s}}}}}{{{{(EI)}_{{\text{eq}}}}}}。$

(7)

当隧道发生整体弯曲变形时，其张开量为隧道下底面的张开距离，由于变形较小，将管节之间弯曲张开量视为张开角所对应圆弧的弧长。根据几何关系，弯曲张开量的理论计算公式为

$\Delta = \theta r = \frac{{M{l_{\text{s}}}r}}{{{{(EI)}_{{\text{eq}}}}}}。$

(8)

式中：r为隧道变形曲线中性轴到管节下底面的距离。

结合式（5），（6）可求得矩形顶管隧道纵向等效抗弯刚度为

${(EI)_{{\text{eq}}}} = f(h) \cdot l。$

(9)

定义纵向等效刚度有效率：

$\eta = \frac{{{{(EI)}_{{\text{eq}}}}}}{{\frac{{{E_{\text{s}}}{A_{\text{s}}} + {E_{\text{c}}}({A_{\text{c}}} - {A_{\text{s}}})}}{{{A_{\text{c}}}}}{I_{\text{c}}}}}。$

(10)

式中： ${(EI)_{{\text{eq}}}}$ 所示的分母为钢套环与管节混凝土的整体弹性模量； ${I_{\text{c}}}$ 则为截面惯性矩； $\eta$ 为纵向等效刚度有效率，它反映了纵向等效抗弯刚度与管节整体刚度关系。

1.2 矩形顶管隧道纵向等效剪切刚度

为求得矩形顶管隧道纵向等效剪切刚度模型，需要进行以下基本假定：

（1）矩形顶管隧道的横截面是连续的、均匀的。

（2）接头始终处于弹性状态，且假定为理想弹性状态。

（3）接头的剪切力主要由钢套环承受，忽略其它影响因素^[21]。

（4）引入钢套环或管节沿截面均匀分布时的Timoshenko梁截面剪切系数 $\kappa$ ^[22]。

$\kappa =\frac{24{(1+{m}^{2})}^{2}{(1+\nu )}^{2}}{{({m}^{2}+1)}^{2}(32{\nu }^{2}+56\nu +28)+{m}^{2}(76{\nu }^{2}+196\nu +107)}。$

(11)

式中： $\nu$ 为材料的泊松比，m=b/a；a为管节横截面（环状截面）外边界到形心的平均距离，b为管节横截面（环状截面）内边界到形心的平均距离。

基于弹性体剪切力学平衡条件和变形协调条件，取 ${l_{\text{s}}}$ 为单位长度，当管节受到剪力Q时，发生剪切变形，邻近的管节发生了错台变形，形成一个错台量δ，管节受剪变形示意图如图 5所示。

图 5 管节单元受剪变形示意图

Figure 5. Shear deformations of pipe element

下载: 全尺寸图片幻灯片

单元的剪切位移，即相邻管节间的剪切错台量主要由接头处的竖向位移组成，分别是钢套环受剪产生的剪切位移 ${u_{\text{s}}}$ 和相邻管节间的受剪变形位移 ${u_{\text{c}}}$ ^[23]。单元受剪变形位移分析示意图如图 6所示，剪切位移u表示为

$u = {u_{\text{s}}} + {u_{\text{c}}}。$

(12)

图 6 计算单元剪切变形位移示意图

Figure 6. Shear deformation displacement of calculation unit

下载: 全尺寸图片幻灯片

其中管节和接头的剪切位移计算公式如下：

${u_{\text{s}}} = {l_{\text{s}}}\tan {\gamma _{\text{s}}} = {l_{\text{s}}}\tan \frac{Q}{{{\kappa _{\text{s}}}{G_{\text{s}}}{A_{\text{s}}}}} 。$

(13)

式中： ${u_{\text{s}}}$ 为钢套环受剪产生的剪切位移； ${l_{\text{s}}}$ 为钢套环沿隧道纵向分布的长度； ${\gamma _{\text{s}}}$ 为钢套环剪切角； ${\kappa _{\text{s}}}$ 为钢套环的Timoshenko梁截面剪切系数（参考式（11））； ${G_{\text{s}}}$ 为钢套环剪切模量； ${A_{\text{s}}}$ 为钢套环横截面面积。

${u_{\text{c}}} = (l - {l_{\text{s}}})\tan \gamma = (l - {l_{\text{s}}})\tan \frac{Q}{{{\kappa _{\text{c}}}{G_{\text{c}}}{A_{\text{c}}}}} 。$

(14)

式中： ${u_{\text{c}}}$ 为相邻管节间受剪产生的剪切位移； $l$ 为计算单元沿隧道纵向分布的长度； $\gamma$ 为相邻管节间的剪切角； ${\kappa _{\text{c}}}$ 为混凝土管节的Timoshenko梁截面剪切系数； ${G_{\text{c}}}$ 为混凝土管节的剪切模量； ${A_{\text{c}}}$ 为混凝土管节的横截面面积。

由式（12）~（14）可得等效抗剪刚度 ${(\kappa GA)_{{\text{eq}}}}$ ：

${(\kappa GA)_{{\text{eq}}}} = \frac{l}{{\frac{{{l_{\text{s}}}}}{{{\kappa _{\text{s}}}{G_{\text{s}}}{A_{\text{s}}}}} + \frac{{l - {l_{\text{s}}}}}{{{\kappa _{\text{c}}}{G_{\text{c}}}{A_{\text{c}}}}}}} 。$

(15)

上述理论推导是在理想假定条件下进行的，由于矩形顶管隧道管节之间的接头配置相当复杂，因此可根据工程实际按照系数进行折减。

在本篇中，由于混凝土部位的变形远小于接头部位的变形。因此，忽略管节间的剪切变形，以钢套环受剪切变形为主，进行剪切错台计算。剪切错台量理论计算公式如下所示：

$\delta = {u_{\text{s}}} = {l_{\text{s}}}\tan {\gamma _{\text{s}}} = {l_{\text{s}}}\tan \frac{Q}{{{{(\kappa GA)}_{{\text{eq}}}}}} 。$

(16)

2. 试验构件与方案

根据第一节所述内容，为探究等效刚度的影响因素，首先要验证理论结果的合理性，本节将分别开展弯曲与剪切模型试验，采集无地基作用下接头的张开与错台量。

2.1 试件内容

管节横断面尺寸为1.625 m×1.075 m，管节长为1.5 m，混凝土强度等级为C50，混凝土内部设置直径12 mm，HRB400级受力钢筋、构造钢筋，箍筋采用C8@150。管节间通过钢套环与止水圈组成的F型承插接头进行连接，具体尺寸如图 7所示。

图 7 试验管节示意图

Figure 7. Section of test pipe

下载: 全尺寸图片幻灯片

2.2 试验方案

为获得F型承插接头抗弯及抗剪刚度数据，在内蒙古科技大学土木工程安全与耐久重点实验室分别开展了弯曲和剪切模型试验。其中弯曲试验采用两个管节拼装进行，在管节下方左右两侧各设置一个简易支点，同时在侧壁设置限位装置防止侧翻，如图 8所示。

图 8 接头弯曲试验布置图

Figure 8. Layout of bending tests on joints

下载: 全尺寸图片幻灯片

剪切试验为三管节相接而成，在两侧管节下方放置混凝土支撑底座，并在相应位置的上方设置压梁，来限制两侧管节的竖向位移，从而达到仅对中间管节进行加载的效果，如图 9所示。由于本文主要研究对象为接头刚度，且管身变形极小，因此试验忽略管身变形。

图 9 接头剪切试验布置图

Figure 9. Layout of shear tests on joints

下载: 全尺寸图片幻灯片

试验中主要通过位移计来监测管节位移来确定管节的整体变化情况，测量内容为管身竖向位移与侧向水平位移。具体位置情况如图 10，11所示。

图 10 弯曲试验位移计布置图

Figure 10. Layout of displacement meters for bending test

下载: 全尺寸图片幻灯片

图 11 剪切试验位移计布置图

Figure 11. Layout of displacement meters for shear tests

下载: 全尺寸图片幻灯片

2.3 试验结果分析

根据上述试验方案对顶管模型展开了弯曲与剪切试验，获得了各测点的位移。通过对试验数据进行处理，取两端测点平均值，得到接头处 $M{\text{ - }}\theta$ 曲线与 $Q{\text{ - }}S$ 曲线，如图 12，13所示。

图 12 接头

$M{\text{ - }}\theta$ 曲线图

Figure 12.

$M{\text{ - }}\theta$ curves of joints

下载: 全尺寸图片幻灯片

图 13 接头

$Q{\text{ - }}S$ 曲线图

Figure 13.

$Q{\text{ - }}S$ curves of joints

下载: 全尺寸图片幻灯片

为接头 $M{\text{ - }}\theta$ 曲线图。由图 12可知，接头抗弯刚度在初始阶段，约为无穷大，这是因为千斤顶施加的弯矩未达到防水橡胶圈与钢套环之间的静摩擦力极限。从O点开始，接头张开角开始增加，直到达到A点时开始下降，这是由于在弯矩作用下，接头底部钢套环逐渐达到屈服，导致防水橡胶圈与钢套环逐渐脱开。根据试验位移计数据，由弯矩计算公式可得OA段平均抗弯刚度为387.72 kN·m·rad^-1，AB段平均抗弯刚度为-318.93 kN·m·rad^-1。

为 $Q{\text{ - }}S$ 曲线图。由图 13可知，随着管节位移的不断增大，接头的剪切变形大致可分为3个阶段：第一阶段为线弹性变化，钢套环和鹰嘴橡胶在接触后逐步被挤压到一起，最终承插口各个部位压紧，这一阶段曲线斜率最大，抗剪刚度最大；第二阶段为屈服阶段，钢套环出现鼓包等现象，其变形发展速度较第一阶段要慢，抗剪刚度逐渐减小；第三阶段为破坏阶段，钢套环大面积屈服，达到承载能力极限，同时，钢套环发生剪切破坏。由斜率计算公式可得，弹性阶段（OA）50.5 kN/mm、屈服阶段（AB）33.9 kN/mm。

3. 弯曲与剪切理论模型验证

在第一节中，对纵向变形方向上接头刚度进行了理论推导，通过第二节对试验结果的分析，得到了张开量、错台量的数据，本小节对理论公式结合试验进行验证。

3.1 弯曲模型验证

通过改变张开角的值，将其代入式（8）中，可得到 $\Delta {\text{ - }}M$ 曲线。与试验结果进行对比，得到图 14弯曲张开量对比图。

图 14 弯曲张开量对比图

Figure 14. Comparison of opening amounts of bending

下载: 全尺寸图片幻灯片

图 14为接头弯曲张开量对比图。由图 14可知，弯曲张开量的理论值与试验值大致相同，证明前文公式的合理性。由于推导过程中，假设截面拉伸区，钢套环受到均匀的拉应力作用。而试验中，钢套环因装配误差，导致橡胶圈的与钢套环接触不均匀，同时，作为边界条件的侧边限位装置未起到很好的限制作用，使管节受到了非纯弯曲或纯剪切的作用，从而出现对比偏差。

3.2 剪切模型验证

通过改变错台量的值，代入式（16）中，可以得到 $\delta {\text{ - }}Q$ 曲线。与试验结果进行对比，得到图 15剪切错台量对比图。

图 15 剪切错台量对比图

Figure 15. Comparison of shear dislocation quantities

下载: 全尺寸图片幻灯片

图 15为剪切错台量对比图，由图 15可知，当钢套环处于弹性阶段时，剪切错台量的理论值与试验值大致相同，可证明前文公式的合理性。由于试验设置的橡胶圈具有可压缩性，但却不是影响抗剪刚度的主要因素，因此理论推导过程不考虑橡胶圈的性质。

4. 等效刚度影响因素分析

无论是矩形顶管隧道接头的刚度还是隧道整体等效刚度都与接头密切相关，尤其受钢套环的各项参数的影响。同时，管节的截面尺寸也有可能影响隧道整体的等效刚度，下文将结合《矩形顶管工程技术规程：T/CECS716—2020》^[24]中给出的经典矩形顶管细部尺寸，分别对钢套环的纵向分布长度、厚度以及管节的截面尺寸和厚度进行探究，分析其对等效刚度以及有效率的影响。

4.1 管节截面尺寸

为研究管节断面尺寸对接头弯曲和错台变形的影响。基于规程中对管节尺寸和壁厚的要求，保证钢套环尺寸不变，选取截面为7000 mm×5000 mm，7400 mm×4400 mm，7700 mm×4300 mm以及7700 mm×4500 mm（B×H）来进行分析研究，壁厚分别为500，600，700 mm。

为了研究管节厚度对等效抗弯刚度的影响，基于规程，确定截面尺寸为7700 mm×4500 mm，壁厚为500 mm与600 mm，由于工况数量少，不足以总结规律，因此增加一组壁厚700 mm进行分析。

图 16为不同截面尺寸矩形顶管的等效抗弯刚度及有效率变化图，由图 16可知，等效抗弯刚度与有效率的变化曲线近乎一致，说明截面尺寸对等效抗弯刚度以及有效率的影响接近一致，由于规程未按规律给出管节截面尺寸，此处无法按规律进行总结，可以看出的是，增大截面面积可以有效增加等效抗弯刚度，但需要注意的是，不可以增大管节宽度同时减小管节高度，这样反而会降低等效抗弯刚度，最有效的方式是减小管节宽度同时增大管节高度。

图 16 截面尺寸影响下的等效抗弯刚度及纵向等效刚度有效率变化图

Figure 16. Change of effectiveness of equivalent bending stiffness and longitudinal equivalent stiffness under influences of section size

下载: 全尺寸图片幻灯片

图 17为不同截面尺寸矩形顶管的等效抗剪刚度变化图，由图 17可知，相较于截面尺寸对等效抗弯刚度的影响，截面尺寸对等效抗剪刚度影响并不大，增大截面面积可以适当增加等效抗剪刚度，但要注意不能减小管节宽度同时增大管节高度，这样反而会降低等效抗剪刚度。

图 17 截面尺寸影响下的等效抗剪刚度变化图

Figure 17. Change of equivalent shear stiffness under influences of section size

下载: 全尺寸图片幻灯片

图 18为壁厚影响下的等效刚度变化图，由图 18可知，在壁厚的影响下，等效抗弯刚度与等效抗剪刚度变化规律恰好相反，但从数值上来看，等效抗弯刚度所受影响较大，等效抗剪刚度受到的影响较小。因此，对于管节壁厚的选取，可着重考虑弯曲变形的影响。

图 18 壁厚影响下的等效刚度变化图

Figure 18. Change of equivalent stiffness under influences of wall thickness

下载: 全尺寸图片幻灯片

4.2 钢套环纵向长度

根据规程，固定管节截面尺寸与钢套环厚度，选取钢套环纵向长度为320，330，340，350，360 mm来进行分析研究。

图 19为等效抗弯刚度以及有效率的影响，由图 19可知。钢套环纵向长度对等效抗弯刚度的影响呈正相关，钢套环纵向长度每增加10 mm，其等效抗弯刚度便增加约3.5%。由于接头钢套环在受弯曲作用时为主要受力构件，当钢套环纵向长度增加时，钢套环与插口混凝土沿纵向的接触长度也将增大，影响了钢套环脱离插口混凝土第一台阶的张开量，从而增强了等效抗弯刚度。

图 19 钢套环纵向长度影响下的等效抗弯刚度及纵向等效刚度有效率变化图

Figure 19. Change of effectiveness of equivalent bending stiffness and longitudinal equivalent stiffness under influences of length of steel ring

下载: 全尺寸图片幻灯片

图 20为钢套环纵向长度影响下的等效抗剪刚度变化图，由图 20可知，钢套环纵向长度对等效抗剪刚度的影响呈正相关，钢套环纵向长度每增加10 mm，其等效抗剪刚度便增加约1.4%。由于接头钢套环在受剪切作用时为主要受剪构件，当钢套环纵向长度增加时，钢套环与插口混凝土的接触面积也将增大，影响了接触面的受力，从而增强了等效抗剪刚度。

图 20 钢套环纵向长度影响下的等效抗剪刚度变化图

Figure 20. Change of equivalent shear stiffness under influences of longitudinal length of steel ring

下载: 全尺寸图片幻灯片

在矩形顶管隧道中，由于承口接头处钢套环与插口混凝土的连接类似于一种“环向包裹接触”的结构，而非使用螺栓固定。因此，在设计管节接头时应重点设计钢套环的各项参数。从钢套环对等效刚度的影响上进行分析，钢套环纵向长度每增加10 mm，等效刚度的增加并不显著，表明仅增加钢套环的纵向长度无法很好的提升等效刚度，需要综合考虑。

4.3 钢套环厚度

根据规范，固定管节截面尺寸与钢套环纵向长度，选取钢套环厚度为12~19 mm的范围进行分析研究，对于等效抗弯刚度以及有效率的影响，如图 21所示。

图 21 钢套环厚度影响下等效抗弯刚度及纵向等效刚度有效率变化图

Figure 21. Change of effectiveness of equivalent bending stiffness and longitudinal equivalent stiffness under influences of steel ring thickness

下载: 全尺寸图片幻灯片

图 21为钢套环厚度对等效抗弯刚度及有效率的影响图，由图 21可知，钢套环厚度对等效抗弯刚度的影响呈正相关，钢套环厚度每增加1 mm，其等效抗剪刚度便增加约8%~10%，影响效果十分显著。

图 22为钢套环厚度影响下等效抗弯刚度增长率变化图，由图 22可知，虽然钢套环厚度对等效抗弯刚度影响显著，但是其增长率呈递减状态，所以建议钢套环厚度取值在13~17 mm。

图 22 钢套环厚度影响下等效抗弯刚度增长率变化图

Figure 22. Change of growth rate of equivalent bending stiffness under influences of steel ring thickness

下载: 全尺寸图片幻灯片

图 23为钢套环厚度影响下对等效抗剪刚度及增长率变化图，由图 23可知，钢套环厚度对等效抗剪刚度的影响呈正相关，钢套环厚度每增加1 mm，其等效抗剪刚度便增加约3%~5%。但其增长率呈递减状态，因此不可以无限制的增加钢套环的厚度，结合规程以及工程经验，建议钢套环厚度取值在13~17 mm。

图 23 钢套环厚度影响下对等效抗剪刚度及增长率变化图

Figure 23. Change of equivalent shear stiffness and growth rate under influences of steel ring thickness

下载: 全尺寸图片幻灯片

5. 结论

本文通过室内接头试验对F型承插接头进行研究，结合接头破坏特征，建立了接头刚度模型，推导了接头等效刚度解析解，验证隧道张开与错台量的理论计算公式，主要得到以下4点结论。

（1）建立矩形顶管隧道的纵向等效抗弯刚度表达式，提出等效抗弯刚度有效率η；结合Timoshenko梁截面剪切系数与顶管隧道特点，形成了矩形顶管隧道的纵向等效抗剪刚度表达式。

（2）在弯曲作用下，通过接头抗弯刚度约为387.72 kN·m·rad^-1。在剪切作用下，接头剪切刚度主要分为3个阶段，即弹性阶段、屈服阶段、破坏阶段，并得出接头的抗剪刚度约为弹性阶段刚度50.55 kN/mm，屈服阶段刚度33.93 kN/mm。

（3）结合矩形顶管隧道结构特性，推导隧道弯曲张开量与剪切错台量的理论计算公式，并与试验所测数据进行对比，验证了所推导公式的合理性。

（4）对管节尺寸、钢套环纵向长度以及厚度等影响因素进行分析。发现减小管节宽度同时增大管节高度可以有效的增加等效抗弯刚度，而截面尺寸以及厚度对等效抗剪刚度影响较小，在设计管节截面尺寸时，可以主要考虑弯曲变形的影响。增大钢套环的厚度可以很好的提升等效刚度，但也不可一直增大，建议钢套环厚度取值为13~17 mm。

图 1 粉质黏土及渣土浆液

Figure 1. Silty clay and muck grout

下载: 全尺寸图片幻灯片

图 2 SEM（×5000）图谱

Figure 2. SEM(×5000)atlas

下载: 全尺寸图片幻灯片

图 3 养护28 d试块SEM（×10000）图谱

Figure 3. SEM(×1000)atlas of samples cured for 28 days

下载: 全尺寸图片幻灯片

图 4 盾构浆液渗透扩散模型

Figure 4. Diffusion and penetration model for shield grout

下载: 全尺寸图片幻灯片

表 1 土的主要矿物成分

Table 1 Main mineral components of soil

种类	蒙脱石	石英	长石	方解石	伊利石	沸石	其他
膨润土	74%	10%	9%	—	—	7%
粉质黏土	15%	35%	39%	6%	3%	—	2%

下载: 导出CSV

表 2 浆液原材料简要性能指标

Table 2 Brief performance indexes of grout materials

材料名称	性能指标
水	符合国家标准的水；pH值＞4
水泥	P.O42.5普通硅酸盐水泥
粉质黏土黄砂	现场取样烘干破碎细砂，含泥量须≤6%
粉煤灰	Ⅱ级，0.045 mm方孔筛筛余不大于20%~45%，含水率≤5%
膨润土	钠基，过200目筛量超95%，膨胀率范围：13~30 mL/g

下载: 导出CSV

表 3 原配比下浆液性能试验结果

Table 3 Results of grout performance under original ratio mixture

试验号	稠度/cm	抗压强度/MPa		凝结时间/min	泌水率/%	结石体收缩率/%	初始流动度/cm
试验号	稠度/cm	7 d	28 d	凝结时间/min	泌水率/%	结石体收缩率/%	初始流动度/cm
1	12.4	1.27	4.44	1010	7.89	7.83	29.3
2	12.8	1.67	5.16	1050	9.55	8.30	31.1

下载: 导出CSV

表 4 浆液均匀试验设计配比表

Table 4 Design mixture ratios of uniform tests on grout 单位: kg/m³

试验号	列号
试验号	C	F	D	S	W
1	151.40	454.20	72.67	818.37	363.36
2	121.45	485.82	47.37	820.64	394.73
3	172.86	432.15	33.88	817.59	423.51
4	121.97	426.90	82.33	807.16	411.65
5	180.47	360.93	69.30	796.18	433.12
6	137.53	481.34	29.71	910.10	371.32
7	192.07	384.14	89.89	929.37	374.54
8	140.15	420.45	62.79	904.19	392.42
9	110.43	441.73	33.13	890.59	414.12
10	141.43	353.58	95.04	883.94	396.01
11	113.74	454.95	54.59	1015.51	341.21
12	161.00	402.50	43.95	1006.26	366.28
13	113.75	398.14	86.00	1023.78	358.32
14	167.52	335.04	75.38	1005.13	376.92
15	123.20	369.61	47.31	985.63	394.25

下载: 导出CSV

表 5 粉质黏土浆液性能试验结果

Table 5 Test results of grout performance of silty clay

试验号	密度/(g·cm^-3)	稠度/cm	抗压强度MPa		凝结时间/min	泌水率/%	结石体收缩率/%	初始流动度/cm
试验号	密度/(g·cm^-3)	稠度/cm	7 d	28 d	凝结时间/min	泌水率/%	结石体收缩率/%	初始流动度/cm
1	1.86	12.6	2.91	9.69	855	3.06	2.40	28.5
2	1.87	13.7	2.00	7.80	1140*	12.29*	4.30	31.3
3	1.88	13.3	2.34	8.04	1060*	12.41*	7.80	32.1
4	1.85	13.1	1.46	5.53	985*	10.46*	6.10	31.7
5	1.84	13.3	3.72	9.75	900*	17.28*	8.60	33.2
6	1.93	11.4	3.13	8.39	870	8.57	3.40	26.2
7	1.97	11.8	3.48	10.50	700	4.60	4.60	28.5
8	1.92	12.6	2.77	8.14	990*	7.68	5.20	28.9
9	1.89	12.4	1.86	6.74	1183*	9.75	10.60	30.3
10	1.87	12.9	2.79	8.06	1000*	8.30	7.30	34.2
11	1.98	10.7	2.09	7.88	863	5.65	5.20	25.2
12	1.98	11.9	3.00	9.36	780	8.21	8.10	25.3
13	1.98	12.0	1.86	4.92	960*	5.27	6.10	27.3
14	1.96	12.1	3.40	8.02	785	9.05	7.20	28.2
15	1.92	12.2	1.87	5.20	1050*	14.06*	10.30	31.4
注：不满足现场要求的数据右上角标记“*”。

下载: 导出CSV

表 6 拟合效果评价

Table 6 Evaluation of fitting effects

变量	表达式	R²	P
f_稠度	$8.865 - 7.636{X_1} + 81.217{X_3} - 3.615{X_4} + 14.962X_2^2 - 114.868X_3^2 + 5.2{X_1}{X_4} - 68.6{X_2}{X_3}$	0.907	0.04
f_{初始流动度}	$- 11.168 + 58.785{X_1} + 37.382X_2^2 + 23.323X_3^2 - 36.587{X_1}{X_3}$	0.912	0.01
f_凝结时间	$- 121.939 + 3642.742{X_3} + 448.831{X_1}{X_4} + 1243.688X_2^2 - 6660.487{X_2}{X_3} - 30.687X_4^2$	0.902	0.01
f_泌水率	$- 32.735 + 99.226{X_1} - 91.459{X_3} - 20.078{X_1}{X_4} + 1.359X_4^2 + 3.744{X_2}{X_4} + 19.171{X_3}{X_4}$	0.900	0.01
f_{体积收缩率}	$58.{\text{0}} - 256.{\text{2}}{X_{\text{3}}} - 23.{\text{3}}{X_{\text{4}}} - 2{\text{5}}{\text{.0}}X_{\text{2}}^{\text{2}} + 239.{\text{8}}X_{\text{3}}^{\text{2}} + 7.{\text{6}}{X_{\text{1}}}{X_4} + 132.{\text{5}}{X_{\text{2}}}{X_3} + 2.{\text{3}}X_4^2 + 25.{\text{2}}{X_{\text{3}}}{X_4}$	0.966	0.01
f_{7d抗压强度}	$- 3.915 + 9.152{X_1} + 16.127{X_3} + 3.024{X_4} - 5.7{X_3}{X_4} - 4.274{X_1}{X_4}$	0.888	0.02
f_{28d抗压强度}	$2.370 + 42.704{X_3} + 4.670{X_4} + 8.417{X_1}{X_2} - 14.589{X_3}{X_4} - 5.842{X_1}{X_4}$	0.873	0.03

下载: 导出CSV

表 7 优化浆液试验结果验证

Table 7 Validation of test results of optimized grout

项目	稠度/cm	初始流动度/cm	凝结时间/h	抗压强度/MPa		2 h泌水率/%	结石体收缩率/%
项目	稠度/cm	初始流动度/cm	凝结时间/h	7 d	28 d	2 h泌水率/%	结石体收缩率/%
性能最优试验值	10.6	24.6	14.4	2.02	6.60	3.97	5.07
性能最优理论值	11.25	26.03	14.7	1.95	6.13	3.35	4.95
误差率%	-6.13	-5.81	-2.08	3.47	7.12	15.62	2.37
成本最优试验值	12.5	32.2	14.8	1.92	6.32	7.43	5.71
成本最优理论值	11.98	30.26	15.07	1.84	6.15	7.06	5.59
误差率%	4.16	6.02	-1.82	4.17	2.69	4.98	2.10
注：误差率指（试验值-理论值）/试验值。

下载: 导出CSV

表 8 浆液水下工作性能

Table 8 Underwater working performance of grout

试验号	抗渗性/MPa	7 d水路强度比	28 d水路强度比
原配比	0.50	0.73	0.66
性能最优	0.65	0.79	0.70
成本最优	0.48	0.70	0.61

下载: 导出CSV

表 9 不同注浆压力下扩散参数

Table 9 Diffusion parameters under different grouting pressures

注浆压力/kPa	扩散半径/cm			对管片产生的压力/kN			对管片产生的压强/kPa
注浆压力/kPa	性能最优	成本最优	原配比	性能最优	成本最优	原配比	性能最优	成本最优	原配比
400	103.67	138.28	109.86	651.40	868.82	690.26	385.83	289.28	364.10
430	107.03	142.80	113.43	722.96	964.51	766.13	401.74	301.13	379.10
460	110.27	147.15	116.86	796.77	1063.25	844.39	417.16	312.61	393.63
490	113.39	151.35	120.17	872.76	1164.93	924.97	432.13	323.75	407.74
520	116.41	155.41	123.38	950.86	1269.44	1007.79	446.70	334.59	421.46
550	119.34	159.35	126.49	1030.99	1376.70	1092.77	460.88	345.15	434.83

下载: 导出CSV

表 10 不同注浆时间下扩散参数

Table 10 Diffusion parameters under different grouting time

注浆时间/min	扩散半径/cm			对管片产生的压力/kN			对管片产生的压强/kPa
注浆时间/min	性能最优	成本最优	原配比	性能最优	成本最优	原配比	性能最优	成本最优	原配比
10	66.88	88.83	70.81	420.21	558.13	444.89	598.10	450.31	564.92
20	88.59	118.00	93.85	556.64	741.39	589.68	451.51	339.00	426.21
30	103.67	138.28	109.86	651.41	868.82	690.26	385.82	289.27	364.10
40	115.29	153.91	122.20	724.41	967.06	767.78	346.94	259.89	327.35
50	124.69	166.56	132.17	783.44	1046.51	830.44	320.80	240.16	302.64
60	132.50	177.07	140.46	832.49	1112.56	882.53	301.90	225.90	284.78

下载: 导出CSV

参考文献(30)

[1]	协会快报: 2020年上半年中国内地城轨交通线路概况[J]. 城市轨道交通, 2020(7): 42. Association express: overview of China urban rail transit lines in the first half of 2020[J]. China Metros, 2020(7): 42. (in Chinese)
[2]	YANG B, LI H, LI H Z, et al. Experimental investigation on the mechanical and hydraulic properties of urease stabilized fine sand for fully permeable pavement[J]. International Journal of Transportation Science and Technology, 2022, 11(1): 60-71. doi: 10.1016/j.ijtst.2020.12.002
[3]	XU Y T, YANG B, LIU X M, et al. Investigation of the medium calcium based non-burnt brick made by red mud and fly ash: durability and hydration characteristics[J]. International Journal of Minerals, Metallurgy, and Materials, 2019, 26(8): 983-991. doi: 10.1007/s12613-019-1814-9
[4]	陈蕊, 杨凯, 肖为, 等. 工程渣土的资源化处理处置分析[J]. 环境工程, 2020, 38(3): 22-26. https://www.cnki.com.cn/Article/CJFDTOTAL-HJGC202003004.htm CHEN Rui, YANG Kai, XIAO Wei, et al. Analysis on recycling treatment and disposal of engineering slag[J]. Environmental Engineering, 2020, 38(3): 22-26. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HJGC202003004.htm
[5]	朱考飞, 张云毅, 薛子斌, 等. 盾构渣土的环境问题与绿色处理[J]. 城市建筑, 2018(29): 108-110. https://www.cnki.com.cn/Article/CJFDTOTAL-JZCS201829032.htm ZHU Kaofei, ZHANG Yunyi, XUE Zibin, et al. Environmental issues and green treatment of shield residues[J]. Urbanism and Architecture, 2018(29): 108-110. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZCS201829032.htm
[6]	姜腾, 姚占虎, 闵凡路. 废弃黏土在泥水盾构泥浆配制中的再利用研究[J]. 隧道建设, 2014, 34(12): 1148-1152. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201412007.htm JIANG Teng, YAO Zhanhu, MIN Fanlu. Study on cyclic use of waste clay in preparing slurry for slurry shield boring[J]. Tunnel Construction, 2014, 34(12): 1148-1152. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201412007.htm
[7]	张亚洲, 夏鹏举, 魏代伟, 等. 南京纬三路过江通道泥水处理及全线路废弃土再利用技术[J]. 隧道建设, 2015, 35(11): 1229-1233. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201511024.htm ZHANG Yazhou, XIA Pengju, WEI Daiwei, et al. Slurry treatment and waste muck recycling use in construction of weisanlu Yangtze River crossing tunnel in Nanjing[J]. Tunnel Construction, 2015, 35(11): 1229-1233. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201511024.htm
[8]	戴勇, 阳军生, 张聪, 等. 泥水盾构弃渣在同步注浆材料中的再利用研究[J]. 华中科技大学学报(自然科学版), 2019, 47(10): 40-45. https://www.cnki.com.cn/Article/CJFDTOTAL-HZLG201910008.htm DAI Yong, YANG Junsheng, ZHANG Cong, et al. Research on recycling of discarded soil produced from slurry shield in synchronous grouting materials[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2019, 47(10): 40-45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HZLG201910008.htm
[9]	郝彤, 李鑫箫, 冷发光, 等. 郑州市地铁粉质黏土层中盾构渣土制备同步注浆材料特性[J]. 长安大学学报(自然科学版), 2020, 40(3): 53-62. https://www.cnki.com.cn/Article/CJFDTOTAL-XAGL202003007.htm HAO Tong, LI Xinxiao, LENG Faguang, et al. Synchronous grouting materials for shield slag in silty clay of Zhengzhou metro[J]. Journal of Chang'an University (Natural Science Edition), 2020, 40(3): 53-62. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XAGL202003007.htm
[10]	ZHANG C, YANG J S, FU J Y, et al. Recycling of discharged soil from EPB shield tunnels as a sustainable raw material for synchronous grouting[J]. Journal of Cleaner Production, 2020, 268: 121947.
[11]	VINAI R, OGGERI C, PEILA D. Soil conditioning of sand for EPB applications: a laboratory research[J]. Tunnelling and Underground Space Technology, 2008, 23(3): 308-317.
[12]	扈萍, 黄茂松, 马少坤, 等. 粉细砂的真三轴试验与强度特性[J]. 岩土力学, 2011, 32(2): 465-470. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201102029.htm HU Ping, HUANG Maosong, MA Shaokun, et al. True triaxial tests and strength characteristics of silty sand[J]. Rock and Soil Mechanics, 2011, 32(2): 465-470. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201102029.htm
[13]	彭芳乐, 白小宇, 雷亮. 多种应力路径条件下上海粉细砂的室内试验[J]. 同济大学学报(自然科学版), 2011, 39(12): 1767-1772. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201112006.htm PENG Fangle, BAI Xiaoyu, LEI Liang. Laboratory experimental study on Shanghai fine sand along different stress paths[J]. Journal of Tongji University (Natural Science), 2011, 39(12): 1767-1772. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201112006.htm
[14]	张保卫. 干湿循环过程中粉细砂改良土路基填料试验研究[J]. 铁道建筑, 2019, 59(2): 54-56. https://www.cnki.com.cn/Article/CJFDTOTAL-TDJZ201902014.htm ZHANG Baowei. Experimental study on fine silt improved soil subgrade filling during dry and wet cycles[J]. Railway Engineering, 2019, 59(2): 54-56. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TDJZ201902014.htm
[15]	钟小春, 左佳, 刘泉维, 等. 地层中粉细砂在盾构壁后注浆中的再利用研究[J]. 岩土力学, 2008, 29(增刊1): 293-296. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2008S1058.htm ZHONG Xiaochun, ZUO Jia, LIU Quanwei, et al. Reuse of excavated fine sand for back grouting of shield tunnelling[J]. Rock and Soil Mechanics, 2008, 29(S1): 293-296. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2008S1058.htm
[16]	PANG C R, WEN X. Experiment of slurry shield tunnel recycling waste sand in sand-pebble layer[J]. Advanced Materials Research, 2014, 919/920/921: 1104-1109.
[17]	ZHOU S H, LI X, JI C, et al. Back-fill grout experimental test for discharged soils reuse of the large-diameter size slurry shield tunnel[J]. KSCE Journal of Civil Engineering, 2017, 21(3): 725-733.
[18]	ABOUSNINA R, MANALO A, LOKUGE W, et al. Properties and structural behavior of concrete containing fine sand contaminated with light crude oil[J]. Construction and Building Materials, 2018, 189: 1214-1231.
[19]	李雪, 黄琦, 王培鑫, 等. 粉细砂地层泥水盾构渣土回收利用及性能优化[J]. 建筑材料学报, 2019, 22(2): 299-307. https://www.cnki.com.cn/Article/CJFDTOTAL-JZCX201902022.htm LI Xue, HUANG Qi, WANG Peixin, et al. Back-fill grouting and proper performance for discharged soils reuse of the slurry shield tunnel on sand stratum[J]. Journal of Building Materials, 2019, 22(2): 299-307. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZCX201902022.htm
[20]	ESTEPHANE P, GARBOCZI E J, BULLARD J W, et al. Using fine sand shape metrics determined from X-ray microcomputed tomography to illustrate the influence of particle shape on the properties of dispersed mortars[J]. Cement and Concrete Composites, 2021, 123: 104176.
[21]	建筑砂浆基本性能试验方法标准: JGJ/T 70—2009[S]. 北京: 中国建筑工业出版社, 2009. Standard for Test Method of Basic Properties of Construction Mortar: JGJ/T 70—2009[S]. Beijing: China Architecture & Building Press, 2009. (in Chinese)
[22]	预应力孔道灌浆剂: GB/T 25182—2010[S]. 北京: 中国标准出版社, 2010. Grouting Admixture for Prestressed Structure: GB/T 25182—2010[S]. Beijing: Standards Press of China, 2010. (in Chinese)
[23]	王星华. 黏土固化浆液固结过程的SEM研究[J]. 岩土工程学报, 1999, 21(1): 34-40. http://cge.nhri.cn/cn/article/id/10251 WANG Xinghua. SEM study of hardening processes of clay hardening grouts[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 34-40. (in Chinese) http://cge.nhri.cn/cn/article/id/10251
[24]	方开泰. 均匀设计与均匀设计表[M]. 北京: 科学出版社, 1994. FANG Kaitai. Uniform Design and Uniform Design Table[M]. Beijing: Science Press, 1994. (in Chinese)
[25]	水下不分散混凝土试验规程: DL/T 5117—2000[S]. 北京: 中国电力出版社, 2001. Test Code on Non-Dispersible Underwater Concrete: DL/T 5117—2000[S]. Beijing: China Electric Power Press, 2001. (in Chinese)
[26]	阮文军. 注浆扩散与浆液若干基本性能研究[J]. 岩土工程学报, 2005, 27(1): 69-73. http://cge.nhri.cn/cn/article/id/11558 RUAN Wenjun. Research on diffusion of grouting and basic properties of grouts[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(1): 69-73. (in Chinese) http://cge.nhri.cn/cn/article/id/11558
[27]	苟长飞. 盾构隧道壁后注浆浆液扩散机理研究[D]. 西安: 长安大学, 2013. GOU Changfei. Study on the Grouts Diffusion Mechanism of Shield Tunnel Back-Filled Grouts[D]. Xi'an: Changan University, 2013. (in Chinese)
[28]	苟长飞, 叶飞, 张金龙, 等. 盾构隧道同步注浆充填压力环向分布模型[J]. 岩土工程学报, 2013, 35(3): 590-598. http://cge.nhri.cn/cn/article/id/15002 GOU Changfei, YE Fei, ZHANG Jinlong, et al. Ring distribution model of filling pressure for shield tunnels under synchronous grouting[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(3): 590-598. (in Chinese) http://cge.nhri.cn/cn/article/id/15002
[29]	KEFAYATI G H R, HUILGOL R R. Lattice Boltzmann method for the simulation of the steady flow of a Bingham fluid in a pipe of square cross-section[J]. European Journal of Mechanics - B/Fluids, 2017, 65: 412-422.
[30]	LIU J, SONG J, ZHANG Z S, et al. Influence of the ground displacement and deformation of soil around a tunnel caused by shield backfilled grouting during construction[J]. Journal of Performance of Constructed Facilities, 2017, 31(3): 04016117.