面板堆石坝混凝土面板配筋计算方法研究

魏匡民; 周恒; 米占宽; 李国英; 邓曌

doi:10.11779/CJGE20220757

面板堆石坝混凝土面板配筋计算方法研究 English Version

魏匡民^{1, 3,},
周恒²,
米占宽^{1, 3},
李国英^{1, 3},
邓曌¹

1.
南京水利科学研究院, 江苏南京 210024
2.
中国电建集团西北勘测设计研究院有限公司, 陕西西安 710065
3.
水利部水库大坝安全重点实验室, 江苏南京 210029

基金项目:

国家重点研发计划 2021YFC3090100

中央级公益性科研院所基本科研业务费项目 Y322003

中央级公益性科研院所基本科研业务费项目 Y320005

国家自然科学基金项目 52109094

江苏省水利科技项目 2018019

详细信息

作者简介:
魏匡民(1985—)，男，工学博士，高级工程师，主要从事土石坝安全保障理论与技术方面的研究工作。E-mail: kmwei@nhri.cn

中图分类号: TU43;TV641
计量
- 文章访问数: 416
- HTML全文浏览量: 73
- PDF下载量: 89
出版历程
- 收稿日期: 2022-06-13
- 网络出版日期: 2023-02-19
- 刊出日期: 2023-05-31

Reinforcement calculation of concrete slab in concrete faced rockfill dams

WEI Kuangmin^{1, 3,},
ZHOU Heng²,
MI Zhankuan^{1, 3},
LI Guoying^{1, 3},
DENG Zhao¹

1.
Nanjing Hydraulic Research Institute, Nanjing 210024, China
2.
Northwest Engineering Corporation Limited, Xi'an 710065, Chian
3.
Key Laboratory of Reservoir and Dam Safety of the Ministry of Water Resources, Nanjing 210029, China

摘要

摘要: 混凝土面板应力变形安全是面板堆石坝工程的关键，然而传统面板坝分析方法大多未考虑面板的配筋效应，实际面板配筋方案制定缺乏理论依据。提出了基于限裂理念的混凝土面板配筋计算方法：采用非协调网格计算技术精细模拟混凝土面板应力，采用钢筋埋置单元模拟钢筋加固作用，采用损伤模型模拟混凝土材料损伤开裂行为，通过引入与几何尺寸、钢筋应力相关的裂缝宽度公式实现混凝土面板限裂配筋计算。提出的方法可定量研究面板限裂配筋，为复杂条件下面板坝配筋设计提供理论依据。文中以正在施工建设的大石峡特高面板坝为例，初步研究了该坝的面板配筋方案。
- 面板堆石坝 /
- 有限元 /
- 面板配筋计算 /
- 精细模拟 /
- 损伤模型
Abstract: The stress and deformation safety of the concrete slab is crucial to the concrete faced rockfill dams (CFRDs), however, most of the traditional methods do not consider the effects of the reinforcement in slab, therefore, the calculated results can not provide enough information to the design of slab reinforcement. In this study, a method for calculating the slab reinforcement based on the concept of crack width limit is proposed: the non-coordinated grid computing technology is used to refinedly simulate the stress of the concrete slab, the embedded steel element is used to simulate the reinforcement effects, and the damage model is used to simulate the cracking of the slab concrete. By introducing a formula for crack width related to geometric size and rebar stress, and setting the crack width limit, the reinforcement calculation is realized. The proposed method can be used to quantitatively study the slab reinforcement and provide a theoretical basis for reinforcement design under complex conditions. The Dashixia CFRD under construction is taken as an example and its slab reinforcement design is preliminarily studied.
- concrete faced rockfill dam /
- finite element /
- reinforcement calculation /
- refined simulation /
- damage model

HTML全文

0. 引言

软黏土地层盾构隧道不均匀沉降问题普遍^[1-2]，部分工程出现列车降速运行、停运加固等问题，严重影响列车运营安全及隧道结构服役寿命。下方地层注浆是盾构隧道过大沉降治理的常用技术手段。然而，注浆扰动地层产生超孔压，注浆结束后超孔压消散引起地基固结和隧道沉降，导致最终抬升量小于注浆过程中的抬升量峰值。例如，为调控上海地铁2号线某区间过大沉降，在隧道下方采取分阶段微扰动注浆，各注浆阶段结束后隧道均发生明显沉降，最终抬升量相比抬升量峰值减少53%^[3]。

目前，工程案例的注浆参数选取主要依赖经验，缺乏定量预测方法。例如，上海地铁2号线盾构隧道下卧地层注浆最大抬升量为10 mm，隧道抬升量与总注浆体积的比值大约为1︰40 mm/L^[4]。盾构隧道注浆抬升过程本质上是注浆使隧道下卧地层隆起，产生竖向抬升力，进而导致隧道抬升。根据地层-隧道相互作用的两阶段分析方法^[5-6]，该过程隧道抬升量与地层隆起变形和基床系数相关。同时，软黏土地层注浆引起的超孔隙水压力分布与土体结构性密切相关。例如，刘维正等^[7]利用考虑软黏土结构性损伤的柱孔扩张模型开展分析，结果表明：土体结构性越强，柱形孔孔壁处的超孔压越大且沿径向衰减越快。因此，考虑到不同地区软黏土在基床系数和结构性相关参数（如灵敏度和屈服应力比）等影响注浆抬升效率方面存在较大差别，各地区经验很难适用其他地区。除此之外，目前鲜有下方注浆引起的隧道抬升-沉降计算方法。付艳斌等^[8]假定隧道侧下方注浆体为圆柱体，将圆形隧道通过面积等效为正方形，建立了盾构隧道注浆抬升量计算模型，但并未考虑注浆后地基固结沉降引起的隧道再沉降。因此，有必要提出软黏土地层注浆引起的上覆盾构隧道抬升量与沉降量计算方法，探明影响盾构隧道抬升效率的主要因素，为软黏土地层盾构隧道的注浆抬升提供参数支撑。

针对软黏土地层盾构隧道注浆抬升问题，综合运用源汇法、圆柱孔扩张理论和隧道-地层相互作用理论，提出了盾构隧道底部注浆抬升量-沉降量两阶段计算方法，并应用于实际案例。通过参数分析，探讨了注浆参数、地层参数与隧道参数对盾构隧道注浆抬升效率的影响规律。研究成果可以为软黏土地层盾构隧道注浆抬升预测与控制技术提供支撑。

1. 计算方法

软黏土地层盾构隧道底部注浆竖向位移发展呈现两个阶段（）：①注浆过程隧道抬升。软黏土渗透系数低，注浆挤压地层引起的周围地层变形可以认为满足不排水条件。②注浆后隧道沉降。注浆过程会导致周围土体产生超孔压，注浆后超孔压消散引起隧道下卧地基固结沉降，进而导致隧道沉降。因此，本文将盾构隧道底部注浆最终抬升量 $\Delta \delta$ 分为两个部分进行计算，即先根据注浆引起的体积膨胀计算注浆过程隧道抬升量 $\Delta {\delta _1}$ ，再根据注浆引起的地层超孔压计算注浆后的隧道沉降量 $\Delta {\delta _2}$ 。隧道最终量为 $\Delta \delta$ = $\Delta {\delta _1}$ − $\Delta {\delta _2}$ ，抬升效率η为 $\Delta \delta$ / $\Delta {\delta _1}$ 。

图 1 盾构隧道底部注浆抬升-沉降演化过程示意图

Figure 1. Schematic of evolution process of grouting-induced uplift and settlement of overlying shield tunnel

下载: 全尺寸图片幻灯片

1.1 注浆过程盾构隧道抬升量 $\Delta {\delta _1}$ 计算方法

针对注浆抬升过程较短的工况，注浆过程地层响应满足不排水条件，利用源汇法^[9]计算注浆作用下隧道底部自由场地层的纵向变形，其中注浆体简化为若干个等直径球体膨胀的线性叠加（图 2）。需注意的是，本方法不适用于隧道注浆抬升过程较长的情况。然后，基于文克尔地基上的隧道铁木辛柯梁模型^[1]，建立注浆过程隧道底部附加应力及抬升量的计算公式，通过有限差分法求解。本方法主要计算流程如图 3所示。

图 2 注浆过程隧道抬升量计算模型

Figure 2. Theoretical model for calculating tunnel uplift during underlying grouting

下载: 全尺寸图片幻灯片

图 3 注浆过程盾构隧道抬升量计算流程

Figure 3. Calculation process of uplift of shield tunnel during underlying grouting

下载: 全尺寸图片幻灯片

（1）基本假定

a）土体为半无限空间中各向同性的线弹性体，注浆过程土体响应满足不排水条件。

b）注浆方式为囊袋注浆对应的理想压密注浆，注浆体为圆柱体，浆液不会在土体中发生劈裂，不考虑浆液凝固引起的体积变化，每个注浆体由若干半径相等的子球体叠加而成（图 2）。

c）盾构隧道简化为文克尔地基模型上的一维均质铁木辛柯梁模型，隧道纵向附加应力和变形计算均针对隧道底部。

d）隧道纵向附加应力计算采用的地层变形为注浆引起的隧道拱底位置地层隆起量，附加应力沿纵向分布呈现非均匀性。

（2）自由场土体位移

Segaseta^[9]利用镜像原理，提出了半无限空间中由于地层损失引起任意一点土体位移的计算公式，称为源汇法。该方法被广泛应用于不排水条件下基坑围护墙后地层变形计算^[10]。同时，Segaseta^[9]还提出引入负的地层损失可以计算非渗透材料的膨胀，如本文的低渗透地层注浆问题，可以将隧底注浆简化为若干球体膨胀，进而得到隧底土体隆起变形：

$\begin{gathered} {S_z} = \frac{{{a^3}}}{3}\frac{{z - {z_0}}}{{{r_1}^3}} - \frac{{{a^3}}}{3}\frac{{z{\text{ + }}{z_0}}}{{{r_1}^3}}{\text{ + }}\mathop {\lim }\limits_{b \to \infty } \mathop {\lim }\limits_{c \to \infty } \int_{{y_0} - b}^{{y_0} + b} {\int_{{x_0} - c}^{{x_0} + c} {\frac{{{a^3}}}{\mathtt{π} } \cdot } } \hfill \\ \left\{ {\frac{{{z_0}(u - {x_0})(x - u)}}{{{{[{{(u - {x_0})}^2} + {{(t - {y_0})}^2} + {z_0}^2]}^{\frac{5}{2}}}}}\left[ {\frac{z}{{{R^3}}} + \frac{{1 - 2\mu }}{{R(R + z)}}} \right]{\text{ + }}} \right. \hfill \\ \left. {\frac{{{z_0}(t - {y_0})(y - t)}}{{{{[{{(u - {x_0})}^2} + {{(t - {y_0})}^2} + {z_0}^2]}^{\frac{5}{2}}}}}\left[ {\frac{z}{{{R^3}}} + \frac{{1 - 2\mu }}{{R(R + z)}}} \right]} \right\}dudt \text{，} \end{gathered}$

(1)

$\left.\begin{array}{l} r_1=\left[\left(x-x_0\right)^2+\left(y-y_0\right)^2+\left(z-z_0\right)^2\right]^{1 / 2}, \\ r_2=\left[\left(x-x_0\right)^2+\left(y-y_0\right)^2+\left(z+z_0\right)^2\right]^{1 / 2} ，\\ R=\left[(x-u)^2+(y-t)^2+z^2\right]^{1 / 2} 。 \end{array}\right\}$

(2)

式中：a为膨胀半径（m）；μ为土体泊松比；（x₀，y₀，z₀），（x，y，z）分别为膨胀点、计算点坐标；r₁，r₂，R为计算点与源点、像点、地面切应力积分点的距离（m）；S_z(x)为隧底土体隆起量（m）。

（3）隧道附加应力及隧道抬升

为简化计算，假定隧道上覆及下卧地层均被简化为文克尔地基模型。此外，考虑到隧道纵向长度远大于直径且纵向变形由弯曲变形和剪切变形共同组成，本文将盾构隧道简化为铁木辛柯梁模型^[11]。

注浆压力作用下受挤压土体会产生作用于隧道底部、竖直向上的抬升力，该附加应力q(x)为

$q(x) = {k_1}D{S_z}(x) 。$

(3)

式中：D为隧道外径（m）；k₁为下卧地层基床系数（kN/m³）。

已知q(x)后，隧道竖向变形w(x)的微分方程为

$(EI)'\frac{{{\text{d}^4}w(x)}}{{\text{d}{x^4}}} - \frac{{{k_2}D(EI)'}}{{({\text{κGA}})'}}\frac{{{\text{d}^2}w(x)}}{{\text{d}{x^2}}} + {k_2}Dw(x) \\ = q(x) - \frac{{(EI)'}}{{({\text{κGA}})'}}\frac{{{\text{d}^2}q(x)}}{{\text{d}{x^2}}}。$

(4)

式中：k₂为隧道上覆地层基床系数（kN/m³）； $(EI)'$ 为隧道纵向等效抗弯刚度（kN·m²）； $({\text{κGA}})'$ 为隧道等效剪切刚度（kN）。本文采用有限差分法求解。

（4）参数确定

本计算方法所需参数分为隧道参数、土体参数和注浆参数3类。其中，隧道参数包括埋深C、直径D、纵向等效抗弯刚度 $(EI)'$ 以及等效剪切刚度 $({\text{κGA}})'$ 。土体参数包括泊松比μ、基床系数k。注浆参数包括单孔注浆体积V₀、注浆体数量N、注浆体高度h、注浆子层数n_z、注浆体半径a、隧道底-注浆体竖向间距d。主要参数物理意义详见图 2。

a）隧道参数

隧道纵向等效抗弯刚度 $(EI)'$ 及纵向等效剪切刚度 $({\text{κGA}})'$ 是影响隧道变形及内力分布的重要参数，不同隧道的刚度取值可以采用经验法^[12]。

b）土体参数

地基基床系数k计算方法为^[13]

$k = \frac{{3.08}}{{{\eta _\text{p}}D}}\frac{{{E_\text{s}}}}{{1 - {\mu ^2}}}\sqrt[{^8}]{{\frac{{{E_\text{s}}{D^4}}}{{(EI)'}}}} \text{，}$

(5)

$\eta_{\mathrm{p}}=\left\{\begin{array}{ll} 2.18 & (C / D \leqslant 0.5) \\ 1+\frac{1}{1.7 D / h} & (C / D \geqslant 0.5) \end{array},\right.$

(6)

${E_\text{s}} = 350{c_\text{u}}。$

(7)

式中：η_p为地基基床系数深度修正系数；D为隧道外径（m）；C为隧道埋深（m）；μ为土体泊松比；E_s为土体弹性模量（kPa），软黏土可由式（7）^[14]估算；c_u为土体不排水抗剪强度（kPa）。

c）注浆参数

实际施工中，通常将注浆管通过隧道衬砌预留孔插至隧底土体，将每个注浆孔分成若干子层，然后自下而上分层注浆。本文假设每个注浆体为以各子层点为球心的若干球体的叠加（图 2），因而注浆总体积V计算公式为

$V = N{V_0} = N{n_\text{z}}\left( {\frac{4}{3}\mathtt{π} {a^3}} \right) 。$

(8)

式中：N为注浆体数量；V₀为单孔注浆体积；n_z为注浆子层数；a为注浆体半径（m）。

1.2 注浆后的隧道沉降量 $\Delta {\delta _2}$ 计算方法

基于注浆后超孔压消散引起土体固结并导致隧道沉降的基本原理，本节提出了注浆后隧道沉降量的计算方法：利用考虑结构性的柱孔扩张理论^[7]计算塑性区范围的超孔压分布（图 4），通过Mindlin解^[15]以及分层总和法得到隧道下卧土体固结沉降量，再基于文克尔地基上的铁木辛柯梁模型，利用有限差分法求解隧道沉降量。整体计算流程如图 5所示。

图 4 注浆后隧道沉降量理论模型

Figure 4. Theoretical model for tunnel settlement after grouting

下载: 全尺寸图片幻灯片

图 5 注浆后隧道沉降量计算流程

Figure 5. Calculation process of tunnel settlement after grouting

下载: 全尺寸图片幻灯片

（1）基本假定

为简化计算，本计算方法做如下假定：

a）土体固结过程中弹性模量不变，仅考虑塑性区范围内的土体固结沉降。

b）注浆体为规则的圆柱体，内部产生各向同性的注浆压力。注浆体上下边界的超孔压均布且分别等于边界1和边界2在注浆体边缘处（r=0）的塑性区超孔压大小。超孔压消散后将完全转化为隧道下卧地层竖向附加应力。

c）不考虑注浆体凝固引起的体积变形对注浆后上覆盾构隧道沉降的影响。

d）上覆盾构隧道简化为文克尔地基模型上的一维均质铁木辛柯梁模型，隧道纵向附加应力、变形计算均针对隧道底部。

e）隧道纵向附加应力计算采用的地层变形为注浆引起的隧道拱底位置地层沉降量，附加应力沿纵向分布呈现非均匀性。

（2）塑性区土体超孔压分布

刘维正等^[7]针对软黏土的施工扰动及结构性，基于剑桥模型的屈服准则，推导了考虑土结构性损伤的软黏土柱孔扩张弹塑性解，并提出了塑性区的超孔压为

$\begin{array}{l} \Delta u(r) = {\sigma _ \text{r}} - \left[ {\frac{{{q_\text{u}}}}{{\sqrt 3 }} + {{p'}_0}{{\left( {\frac{{ \text{YSR}}}{2}} \right)}^\varLambda }} \right]\left[ {1 - \left( {1 - \frac{1}{{{S_\text{t}}}}} \right)\frac{{\ln \frac{r}{{{r_\text{p}}}}}}{{\ln \frac{{{r_\text{p}}}}{a}}}} \right] \hfill \\ = {\sigma _0} + \frac{{{q_\text{u}}}}{{\sqrt 3 }} - \frac{{2{q_\text{u}}}}{{\sqrt 3 }}\left[ {\ln \frac{r}{{{r_\text{p}}}} + \frac{{1 - \frac{1}{{{S_\text{t}}}}}}{{2\ln \frac{{{r_\text{p}}}}{a}}}{{\left( {\ln \frac{r}{{{r_\text{p}}}}} \right)}^2}} \right] - \hfill \\ \left[ {\frac{{{q_\text{u}}}}{{\sqrt 3 }} + {{p'}_0}{{\left( {\frac{{ \text{YSR}}}{2}} \right)}^\varLambda }} \right]\left[ {1 - \left( {1 - \frac{1}{{{S_\text{t}}}}} \right)\frac{{\ln \frac{r}{{{r_\text{p}}}}}}{{\ln \frac{{{r_\text{p}}}}{a}}}} \right] \text{，} \end{array}$

(9)

${q_\text{u}} = M{p'_0}{\left( {\frac{{ \text{YSR}}}{2}} \right)^\varLambda }\text{，}$

(10)

${r_\text{p}} = a\sqrt {\frac{{\sqrt 3 G}}{{{q_\text{p}}}}} \text{，}$

(11)

${q_\text{p}} = M{p'_0}\sqrt { \text{YSR} - 1} \text{，}$

(12)

$\varLambda = 1 - \frac{\kappa }{\lambda }\text{，}$

(13)

$\text{YSR} = \frac{{{{p'}_{\text{y}}}}}{{{{p'}_0}}}。$

(14)

式中：Δu为超孔压（kPa）；σ₀为初始内应力（kPa）；q_u为极限偏应力（kPa）；q_p为弹塑性交界处应力（kPa）； ${p'_0}$ 为初始有效平均应力（kPa）， ${p'_0}$ =s₀；p_y'为等向固结屈服应力（kPa）；a为注浆体半径（m）；r_p为塑性区半径（m）；S_t为土体灵敏度；M为破坏应力比；G为剪切模量（kPa）；YSR为屈服应力比；Λ为塑性体积比；λ为e-lnp'坐标中正常固结线的斜率；κ为e-lnp'坐标中回弹线的斜率。

（3）隧道底部土体固结沉降

本计算方法将发生固结沉降的土层分为3个部分，包括压缩层1、压缩层2及压缩层3（图 4）。其中，压缩层1为与注浆体同高度的塑性区土体，压缩层2为注浆体顶部与隧道间的土体，压缩层3为注浆体底部以下的土体。注浆后隧道下卧土体的固结沉降量由压缩层1~压缩层3的变形直接叠加得到。

（4）地层超孔压及附加应力计算方法验证

为验证土体超孔压和地层附加应力计算方法的合理性，本文对比了福州淤泥夹砂地层囊袋注浆室内试验^[16]获得的超孔压与计算值。

a）超孔压水平分布规律

本节计算了式（9）在室内试验条件下的理论值，并与室内试验孔压计测得的超孔压结果对比（图 6）。可以发现：下囊袋注浆时，B断面的试验结果与理论计算曲线趋势相同且数值较为接近，证明了式（9）的适用性。在桶壁处两者数据存在差距，主要是模型桶边界效应影响所致。

图 6 B断面超孔压分布曲线

Figure 6. Distribution of excess pore water pressure at Section B

下载: 全尺寸图片幻灯片

b）竖向应力增量云图

针对上述注浆引起的地层竖向应力增量计算假设，本文结合式（9），（16）计算了室内注浆试验条件下地层竖向应力增量并与试验结果进行对比（图 7（a），（b））。提取模型桶对称轴处（d=0）的竖向应力增量值并绘制其沿深度的变化曲线（图 7（c）），可以发现：室内试验与本文方法计算结果趋势相同、数值接近，但两者仍存在一定差异，主要是受模型桶顶板边界效应的影响。

图 7 竖向应力增量计算结果与试验结果对比

Figure 7. Comparison between measured and calculated additional vertical stresses

下载: 全尺寸图片幻灯片

（5）隧道附加应力及沉降

与注浆过程隧道抬升量计算方法一致，注浆后隧道沉降计算中仍将隧道-地层相互作用简化为文克尔地基上的隧道铁木辛柯梁模型。

（6）参数确定

本计算方法涉及隧道参数、土体参数以及注浆参数3类参数。

a）土体参数

土体参数包括剪切模量G₀、有效内摩擦角φ'、临界状态比M、对数硬化模量λ、各向等压膨胀指数k、塑性比体积Λ、等向固结屈服应力 ${p'_{\text{y}}}$ 、等向固结屈服应力比YSR、灵敏度S_t。

其中，剪切模量G₀为^[17]

$\frac{{{G_0}}}{{{p_\text{a}}}} = Sf(e){\left( {\frac{{p'}}{{{p_\text{a}}}}} \right)^n} = S{(1 + {e_0})^{ - 2.4}}{\left( {\frac{{p'}}{{{p_{\text{a}}}}}} \right)^n} 。$

(15)

式中：p_a为参考应力（kPa），取100 kPa；p'为平均有效应力（kPa）， $p'$ =（ ${s'_{\text{v}}}$ + 2 ${s'_{\text{h}}}$ ）/3；n，S为经验参数，对于软黏土，n取0.5，S取24000~30000；e₀为初始孔隙比。

有效内摩擦角φ'和临界状态比M可分别为^[18]

$\varphi ' = 30.4{e_0}^{ - 0.65} \text{，}$

(16)

$M = \frac{{6\sin \varphi '}}{{3 - \sin \varphi '}}。$

(17)

等向固结屈服应力 ${p'_{\text{y}}}$ 可以基于一维固结试验得到的屈服应力p_y并结合剑桥模型的屈服函数计算：

$f \equiv {q^2} + {M^2}p'(p' - {p'_{\text{y}}}) = 0。$

(18)

三轴应力状态下：

$\begin{array}{l}{q}^{\prime }={{\sigma }^{\prime }}_{1}-{{\sigma }^{\prime }}_{3}=(1-{k}_{0}){{\sigma }^{\prime }}_{\text{v}}\text{ }\text{，}\\ {p}^{\prime }=\frac{{{\sigma }^{\prime }}_{1}\text{+2}{{\sigma }^{\prime }}_{3}}{3}=\frac{(1+2{k}_{0}){{\sigma }^{\prime }}_{\text{v}}}{3}。\end{array}\}$

(19)

进而可以得到等向固结屈服应力 ${p'_{\text{y}}}$ ：

${p'_{\text{y}}} = \frac{{(1 + 2{K_0}){p_\text{y}}}}{3} + \frac{{3(1 - {K_0}){p_\text{y}}}}{{{M^2}(1 + 2{K_0})}}。$

(20)

式中：K₀为土体静止侧压力系数；p_y为一维固结试验得到的屈服应力（kPa）；M为临界状态比。

等向固结屈服应力比YSR为等向固结屈服应力 ${p'_{\text{y}}}$ 与初始固结应力p₀之比。

b）隧道参数

隧道参数包括埋深C、直径D、纵向等效抗弯刚度 $({\text{EI}})'$ 以及纵向等效剪切刚度 $({\text{κGA}})'$ 。

c）注浆参数

注浆参数包括初始孔径a₀、注浆体半径a、注浆体高度h、注浆体-隧道竖向间距d。其中，初始孔径a₀为注浆管实际外径，注浆体-隧道竖向间距d为注浆体顶部与隧道底部的距离。

2. 案例应用

本文将上述计算方法应用于宁波地铁2号线某区间盾构隧道注浆抬升工程^[4]。

2.1 工程概况

（1）工程地质条件

宁波市位于中国东部沿海，地貌为第四纪淤积平原，广布深厚淤泥质黏土层。宁波地铁2号线某区间盾构隧道埋深约10~11 m，隧道上部为①₁杂填土、①₂黏土层，隧道主要位于④₁淤泥质黏土中，底部土层为⑤₂粉质黏土，地下水位位于地表以下1 m处，地质剖面如图 8所示。

图 8 宁波地铁2号线地质剖面及注浆方案（166环~174环）

Figure 8. Geological profile and grouting scheme of Ningbo Metro Line 2 (Ring 166 to Ring 174)

下载: 全尺寸图片幻灯片

（2）盾构隧道

盾构隧道采用预制钢筋混凝土管片错缝拼装，外径6.2 m，厚度0.35 m，环宽1.2 m。如图 9所示，隧道环向均布16个预留注浆孔H1~H16，注浆孔直径60 mm。隧道贯通数月后测量发现该区间出现较大差异沉降，最大沉降达214 mm，出现在169环。

图 9 盾构隧道抬升量实测值与计算值对比

Figure 9. Comparison between measured and calculated uplifts of shield tunnel

下载: 全尺寸图片幻灯片

（3）隧道抬升方案

该工程借鉴上海地铁2号线经验，采用微扰动注浆技术抬升隧道，注浆材料为水泥、水玻璃双液浆。

根据隧道沉降情况确定综合治理范围为R166~R175，由管片预留注浆孔H8和H10两个孔竖直注浆（图 8）。根据治理后隧道最小曲率半径目标值及工程经验确定基本注浆参数为：每环设计注浆总量为1900 L，每孔设计注浆量为950 L，将每个注浆孔分成10个子层，每个子层注浆量为95 L，注浆深度范围为隧道底部以下0~2 m。

（4）计算参数

根据隧道结构特征和场地地层条件，本文获得了隧道、土体及注浆计算参数（表 1）。

表 1 宁波地铁2号线盾构隧道注浆抬升工程参数

Table 1. Parameters of grouting uplift project for shield tunnel of Ningbo Metro Line 2

参数名称		参数取值
隧道参数	C/m	10
	D/m	6.2
	(EI)'/(kN·m²)	1.91×10⁸
	(κGA)'/kN	4.14×10⁶
土体参数（淤泥质黏土）	v	0.4
	E_s/MPa	2.83
	S_t	2.22
	k/(kN·m^-3)	4273
	G/MPa	40
	κ	0.01
	λ	0.2
	Μ	1.5
注浆参数	V/L	13300
	V₀/L	950
	N/个	14
	n_z/个	10
	r/m	0.283
	h/m	2
	注浆方式	连续双孔

下载: 导出CSV

| 显示表格

2.2 结果对比

通过对比本文计算值与实测结果（图 9），可以发现本文计算方法得到的最终隧道抬升曲线与实测曲线的趋势基本吻合，数值也较为接近，一定程度上说明了计算方法的有效性。另外，根据本文计算结果，可以得到隧道抬升效率η约51%，与同为淤泥质黏土地层的上海2号线隧道抬升案例^[19]的抬升效率十分吻合。此外，实测结果在治理区段略小于计算值且曲线更平缓，可能是由于实际工程中对盾构隧道内部进行了钢支撑加固，导致该区段纵向刚度变大，隧道更不容易变形且结构差异变形更小。

3. 影响因素分析

本节分别讨论土体屈服应力比、注浆体积以及隧道埋深对隧道最终抬升量及抬升效率的影响规律。

3.1 屈服应力比

图 10（a），（b）分别为地层不同屈服应力比YSR下的隧道最终抬升量δ及抬升效率η纵向分布曲线。可以发现：随着屈服应力比增大，注浆后隧道的最终抬升量及抬升效率减小。根据式（9），这主要是注浆产生的超孔隙水压力大小与土体屈服应力比呈正相关，导致注浆后隧道出现更大沉降，降低隧道抬升效率。另外，随着屈服应力比增大，隧道抬升效率沿纵向分布更加不均匀。

图 10 土体屈服应力比对隧道最终抬升量及抬升效率的影响

Figure 10. Effects of yield stress ratio on final uplift amount and uplifting efficiency of tunnel

下载: 全尺寸图片幻灯片

3.2 注浆体积

图 11（a），（b）分别为不同注浆体积V₀下隧道最终抬升量δ及抬升效率η纵向分布曲线。可以发现：随注浆体积增加，隧道最终抬升量不断增大，隧道抬升效率也不断提高且变化率越来越小。需注意的是，当注浆体积过大时，所需的注浆压力相应增大，可能造成地层劈裂。因此，工程中应综合考虑注浆压力和地层劈裂等影响，确定合理的注浆体积。

图 11 注浆体积对隧道最终抬升量及抬升效率的影响

Figure 11. Effects of grouting volume on final uplift amount and uplifting efficiency of tunnel

下载: 全尺寸图片幻灯片

3.3 隧道埋深比

图 12（a），（b）分别为不同隧道埋深比C/D下的隧道最终抬升量δ及抬升效率η纵向分布曲线。可以发现：随着隧道埋深比增大，注浆范围内的隧道最终抬升量减小，而抬升效率基本不变。因此，结合前述注浆体积对隧道最终抬升量的参数分析结果，若需保证不同埋深比隧道产生相同的抬升量，注浆体积应该随隧道埋深比增大而增大。

图 12 隧道埋深比对隧道最终抬升量及抬升效率的影响

Figure 12. Effects of cover-to-diameter ratio of tunnel on final uplift amount and uplifting efficiency

下载: 全尺寸图片幻灯片

4. 结论

围绕软黏土地层盾构隧道注浆抬升问题，提出了考虑土体结构性的盾构隧道底部注浆抬升量-沉降量两阶段计算方法，得到2点结论。

（1）本文计算方法应用于宁波地铁2号线隧道过大沉降区段注浆抬升工程，计算得到的最终抬升量纵向曲线与实测曲线趋势基本吻合。实测结果在抬升区段略小于计算值且曲线更加平缓，主要是因为抬升前隧道内部钢支撑加固提高了其结构刚度。

（2）隧道最终抬升量及抬升效率与基床系数、注浆体积正相关，与地层结构性负相关。同时，随隧道埋深比增大，注浆治理区域内的隧道最终抬升量减小，抬升效率基本不变。

图 1 水布垭Ⅰ期面板裂缝分布

Figure 1. Cracks in Stage Ⅰ slab of Shuibuya CFRD

下载: 全尺寸图片幻灯片

图 2 HJD大坝左岸水平向裂缝分布

Figure 2. Horizontal cracks in left abutment slab of HJD CFRD

下载: 全尺寸图片幻灯片

图 3 JL面板坝面板裂缝分布

Figure 3. Cracks in slab of JL CFRD

下载: 全尺寸图片幻灯片

图 4 多点约束法三维示例

Figure 4. 3D example of multi-point constraint refinement

下载: 全尺寸图片幻灯片

图 5 钢筋埋置单元

Figure 5. Steel bar embedding element

下载: 全尺寸图片幻灯片

图 6 混凝土单轴应力应变曲线

Figure 6. Uniaxial stress-strain curve of concrete

下载: 全尺寸图片幻灯片

图 7 大石峡面板坝防渗体布置

Figure 7. Arrangement of impervious body of Dashixia CFRD

下载: 全尺寸图片幻灯片

图 8 大石峡大坝材料分区

Figure 8. Material zoning of Dashixia CFRD

下载: 全尺寸图片幻灯片

图 9 基于非协调网格的大石峡面板坝计算模型

Figure 9. Refined FEM model for Dashixia CFRD

下载: 全尺寸图片幻灯片

图 10 面板内钢筋分布

Figure 10. Reinforcements in slab

下载: 全尺寸图片幻灯片

图 11 蓄水运行期面板变形分布

Figure 11. Deformations of slab during operation period

下载: 全尺寸图片幻灯片

图 12 蓄水运行期面板应力分布

Figure 12. Slab stresses during operation period

下载: 全尺寸图片幻灯片

图 13 不同计算模型1640 m高程面板轴向应力比较

Figure 13. Comparison of axial stresses at elevation 1640 m by different models

下载: 全尺寸图片幻灯片

图 14 蓄水运行期面板受拉损伤参数d_t分布

Figure 14. Tensile damage factor d_t of slab during operation period

下载: 全尺寸图片幻灯片

图 15 钢筋拉应力分布

Figure 15. Tensile stresses of reinforcements

下载: 全尺寸图片幻灯片

图 16 面板裂缝分布

Figure 16. Cracks in slabs of Dashixia CFRD

下载: 全尺寸图片幻灯片

图 17 面板配筋加强区示意

Figure 17. Special reinforcement area of Dashixia CFRD

下载: 全尺寸图片幻灯片

表 1 筑坝料“南水”模型参数

Table 1 Parameters of "NHRI" model

坝体分区	${\rho _{\text{d}}}$ /(g·cm^-3)	${\varphi _0}$ /(°)	$\Delta \varphi$ /(°)	$k$	$n$	${R_{\text{f}}}$	c_d/%	n_d	R_d
垫层区（2A）	2.31	49.3	5.9	920	0.29	0.67	0.49	0.45	0.65
主砂砾石区（3BA、3BB）	2.27	50.1	6.3	1294	0.32	0.74	0.30	0.67	0.72
主堆石区（3BC、3C2）	2.22	53.2	9.0	1104	0.22	0.65	0.40	0.72	0.64
下游堆石区（3C1）	2.22	52.6	8.7	1044	0.21	0.66	0.46	0.70	0.66
增模区	2.28	69.1	17.9	0.14	3928.7	0.35	0.08	0.44	0.11

下载: 导出CSV

表 2 筑坝料流变模型参数

Table 2 Parameters of creep model

坝体分区	α	b /%	c/%	d /%	m₁	m₂	m₃
垫层区（2A）	0.006	0.130	0.028	0.277	0.355	0.600	0.700
主砂砾区（3BA、3BB）	0.0062	0.114	0.025	0.241	0.351	0.578	0.696
主堆石料区（3BC、3C2）	0.0055	0.138	0.030	0.363	0.358	0.629	0.783
下游堆石区（3C1）	0.0054	0.140	0.032	0.365	0.358	0.629	0.783

下载: 导出CSV

表 3 混凝土损伤模型参数

Table 3 Parameters for damage model

E_c/10⁴ MPa	f_cr/MPa	ε_cr	α_c	f_tr/MPa	ε_tr	α_t
3.0	20.1	1.470×10^-3	0.74	2.01	9.5×10^-5	1.25

下载: 导出CSV

表 4 不同配筋方案面板裂缝宽度

Table 4 Cracks width under different reinforcements scheme

	配筋加强区Ⅰ		配筋加强区Ⅱ
	水平向配筋率/%	纵向裂缝最大宽度/mm	顺坡向配筋率/%	水平向最大裂缝宽/mm
原方案	0.587	0.32	0.587	0.24
加强方案	0.783	0.20	0.691	0.18

下载: 导出CSV

参考文献(22)

[1]	罗福海, 张保军, 夏界平. 水布垭面板堆石坝施工期裂缝成因及处理措施[J]. 水利水电快报, 2010, 31(12): 5-8. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSK201012002.htm LUO Fuhai, ZHANG Baojun, XIA Jieping. Causes and treatment measures of cracks in Shuibuya concrete face rockfill dam during construction[J]. Express Water Resources & Hydropower Information, 2010, 31(12): 5-8. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSK201012002.htm
[2]	混凝土面板堆石坝设计规范: N/T 10871—2021[S]. 北京: 中国水利水电出版社, 2022. Code for Design of Concrete Face Rockfill Dams: N/T 10871—2021[S]. Beijing: China Water & Power Press, 2022. (in Chinese)
[3]	钟红, 林皋, 胡志强. 有限元计算中疏密网格过渡方法研究[J]. 计算力学学报, 2007, 24(6): 887-891, 898. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJG200706030.htm ZHONG Hong, LIN Gao, HU Zhiqiang. Two methods for transition between coarse and fine finite elements[J]. Chinese Journal of Computational Mechanics, 2007, 24(6): 887-891, 898. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSJG200706030.htm
[4]	周墨臻, 张丙印, 王伟. 高面板堆石坝软缝接触计算模型及其数值实现[J]. 岩石力学与工程学报, 2016, 35(增刊1): 2803-2810. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2016S1024.htm ZHOU Mozhen, ZHANG Bingyin, WANG Wei. Algorithm and simulation methods for the soft longitudinal joint of the concrete faces in high concrete-faced rockfill dams[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S1): 2803-2810. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2016S1024.htm
[5]	周墨臻, 张丙印, 张宗亮, 等. 超高面板堆石坝面板挤压破坏机理及数值模拟方法研究[J]. 岩土工程学报, 2015, 37(8): 1426-1432. doi: 10.11779/CJGE201508010 ZHOU Mozhen, ZHANG Bingyin, ZHANG Zongliang, et al. Mechanisms and simulation methods for extrusion damage of concrete faces of high concrete-faced rockfill dams[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(8): 1426-1432. (in Chinese) doi: 10.11779/CJGE201508010
[6]	孔宪京, 陈楷, 邹德高, 等. 一种高效的FE-PSBFE耦合方法及在岩土工程弹塑性分析中的应用[J]. 工程力学, 2018, 35(6): 6-14. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201806004.htm KONG Xianjing, CHEN Kai, ZOU Degao, et al. An efficient fe-psbfe coupled method and its application to the elasto-plastic analysis of geotechnical engineering structures[J]. Engineering Mechanics, 2018, 35(6): 6-14. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201806004.htm
[7]	邹德高, 刘锁, 陈楷, 等. 基于四叉树网格和多边形比例边界有限元方法的岩土工程非线性静动力分析[J]. 岩土力学, 2017, 38(增刊2): 33-40. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2017S2005.htm ZOU Degao, LIU Suo, CHEN Kai, et al. Nonlinear static and dynamic analysis for geotechnical engineering based on quadtree mesh and polygon scaled boundary finite element method[J]. Rock and Soil Mechanics, 2017, 38(S2): 33-40. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2017S2005.htm
[8]	邹德高, 陈楷, 刘锁, 等. 非线性比例边界有限元在面板坝分析中的应用[J]. 土木与环境工程学报, 2019, 41(3): 11-18. https://www.cnki.com.cn/Article/CJFDTOTAL-JIAN201903002.htm ZOU Degao, CHEN Kai, LIU Suo, et al. Application of nonlinear scaled boundary polygon element method in analysis of concrete face rockfill dam[J]. Journal of Civil and Environmental Engineering, 2019, 41(3): 11-18. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JIAN201903002.htm
[9]	魏匡民, 陈生水, 李国英, 等. 位移多点约束法在面板堆石坝精细模拟中的应用研究[J]. 岩土工程学报, 2020, 42(4): 616-623. doi: 10.11779/CJGE202004003 WEI Kuangmin, CHEN Shengshui, LI Guoying, et al. Application of displacement multi-point constraint refinement method in simulation of concrete-faced rockfill dams[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 616-623. (in Chinese) doi: 10.11779/CJGE202004003
[10]	王勖成. 有限单元法[M]. 北京: 清华大学出版社, 2003. WANG Xucheng. Finite Element Method[M]. Beijing: Tsinghua University Press, 2003. (in Chinese)
[11]	陈国荣. 有限单元法原理及应用[M]. 北京: 科学出版社, 2009. CHEN Guorong. Principle and Application of Finite Element Method[M]. Beijing: Science Press, 2009. (in Chinese)
[12]	孔宪京, 屈永倩, 邹德高, 等. 强震作用下面板堆石坝跨尺度面板开裂演化分析[J]. 岩土工程学报, 2020, 42(6): 989-996. doi: 10.11779/CJGE202006001 KONG Xianjing, QU Yongqian, ZOU Degao, et al. Cross-scale crack evolution analysis for face slab in concrete faced rockfill dams under strong earthquake[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(6): 989-996. (in Chinese) doi: 10.11779/CJGE202006001
[13]	巫昌海, 汪基伟. 混凝土三维钢筋埋置组合式有限单元模型及其网格自动生成[J]. 计算机辅助设计与图形学学报, 2000, 12(10): 761-764. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJF200010011.htm WU Changhai, WANG Jiwei. Three-dimensional embedded finite element model of reinforced concrete and its automatic mesh generation[J]. Jouraal of Computer Aided Design & Computer Graphics, 2000, 12(10): 761-764. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSJF200010011.htm
[14]	XU B, ZOU D G, KONG X J, et al. Dynamic Damage evaluation on the slabs of the concrete faced rockfill dam with the plastic-damage model[J]. Computers and Geotechnics, 2015, 65: 258-265. http://www.onacademic.com/detail/journal_1000037433272510_9a3e.html
[15]	CEN W J, WEN L S, ZHANG Z Q, et al. Numerical simulation of seismic damage and cracking of concrete slabs of high concrete face rockfill dams[J]. Water Science and Engineering, 2016, 9(3): 205-211. http://qikan.cqvip.com/Qikan/Article/Detail?id=7000080633
[16]	混凝土结构设计规范: GB 50010—2010[S]. 北京: 中国建筑工业出版社, 2011. Code for Design of Concrete Structures: GB 50010—2010[S]. Beijing: China Architecture & Building Press, 2011. (in Chinese)
[17]	朱百里, 沈珠江. 计算土力学[M]. 上海: 上海科学技术出版社, 1990. ZHU Baili, SHEN Zhujiang. Computational Soil Mechanics[M]. Shanghai: Shanghai Science and Technology Press, 1990. (in Chinese)
[18]	李国英, 米占宽, 傅华, 等. 混凝土面板堆石坝堆石料流变特性试验研究[J]. 岩土力学, 2004, 25(11): 1712-1716. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200411008.htm LI Guoying, MI Zhankuan, FU Hua, et al. Experimental studies on rheological behaviors for rockfills in concrete faced rockfill dam[J]. Rock and Soil Mechanics, 2004, 25(11): 1712-1716. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200411008.htm
[19]	顾淦臣, 沈长松, 岑威钧. 土石坝地震工程学[M]. 北京: 中国水利水电出版社, 2009. GU Ganchen, SHEN Changsong, CEN Weijun. Earthquake Engineering for Earthrock Dams[M]. Beijing: China Water & Power Press, 2009. (in Chinese)
[20]	高瑞平. 混凝土构件裂缝宽度计算方法研究[D]. 南京: 东南大学, 2008. GAO Ruiping. Study on Calculation Method of Crack Width of Concrete Members[D]. Nanjing: Southeast University, 2008. (in Chinese)
[21]	丁旭柳, 伍鹤皋, 朱忠华. 钢筋混凝土压力隧洞裂缝宽度计算方法的探讨[J]. 武汉大学学报(工学版), 2001, 34(2): 24-27. https://www.cnki.com.cn/Article/CJFDTOTAL-WSDD200102006.htm DING Xuliu, WU Hegao, ZHU Zhonghua. Study on computational method of crack width of reinforced concrete pressure tunnels[J]. Engineering Journal of Wuhan University, 2001, 34(2): 24-27. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-WSDD200102006.htm
[22]	兰伟钦, 肖明. 考虑钢筋滑移效应的高压隧洞衬砌配筋计算方法[J]. 水电与新能源, 2021, 35(4): 66-71. https://www.cnki.com.cn/Article/CJFDTOTAL-HBFD202104016.htm LAN Weiqin, XIAO Ming. On the reinforcement calculation of high-pressure tunnel lining considering the slip effect of reinforcement[J]. Hydropower and New Energy, 2021, 35(4): 66-71. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HBFD202104016.htm