0.   引言

管棚超前预支护技术能在一定范围内较大幅度提高天然土体强度、刚度，有效减小地表沉降、拱顶沉降，确保隧道掘进支护安全，且因具有施工简便、经济、支护距离长和安全效益高的特点而成为经常采用的超前预支护方法。

目前国内外学者关于管棚预支护技术的作用机理和加固效果开展了大量研究。Hisatake等^[1]、Juneja等^[2]通过离心机模型试验，考虑了开挖效果的影响，研究了管棚支护的作用效果。郭璇等^[3]采用设计正交模型试验的方法，研究了管棚预支护的作用效果以及管棚对土压力分布的影响。周顺华等^[4]在原位观测和室内模型离心模型试验的基础上，提出了管棚作用的棚架原理。董新平^[5]则基于棚架体系建立了管棚的空间分析模型，研究了管棚支护作用的主要特征。Chungsik等^[6]、Ali等^[7]通过建立三维数值模型研究了管棚的受力和变形特征以及预支护效果。在理论分析方面，由于基于弹性地基梁理论的管棚分析模型具有简单实用的特点，故而具有很高的进一步研究完善的价值。李忠凯^[8]通过建立基于Winkler地基模型的管棚分析模型，推导出一个开挖循环过程中管棚内力和变形的计算表达式。郑俊杰等^[9]考虑了开挖对基床系数的影响，提出了变基床系数下管棚的地基梁模型。丁祖德等^[10]则对弹性地基梁模型中基床系数的取值展开详细研究。王炳智等^[11]通过改进弹性地基约束形式，采用Pasternak弹性地基理论建立管棚的分析模型并进行了求解。王道远等^[12]通过对比两个地基模型的计算结果与现场试验值的比较，得到基于Pasternak弹性地基理论的管棚模型更具合理性，并通过讨论地层反力系数以及管棚直径对管棚作用效果的影响，得出通过调整注浆改善地层反力系数或直径的方法来控制管棚变形效果甚微的结论。武松等^[13]根据弹性理论的刚度等效法，将管棚的一维弹性地基梁模型拓展至二维弹性地基各向异性板模型，并采用有限元法对模型进行了求解。

综上所述，针对基于弹性地基梁理论的管棚分析模型的研究已较为成熟，但在模型的建立过程中，对荷载的分布范围和管棚约束形式的认识上还存在矛盾，并且通过对模型的求解，并没有得到关于管棚设计和施工参数对管棚支护效果影响的一般性结论。预支护区段开挖支护完成后管棚全长范围内力和变形的最终分布并不明确。本文结合管棚超前预支护技术在实际运用过程中的问题以及国内外学者对管棚预支护技术研究的基础上，考虑掌子面失稳段对管棚约束的弱化以及隧道开挖支护动态施工过程，通过建立基于Pasternak弹性地基梁理论的管棚分析模型，采用理论推导的方法研究管棚预支护技术的支护效果以及受力特点，并系统地研究各设计施工参数对管棚力学响应的影响，从而为管棚的设计施工提供参考和借鉴。

1.   管棚分析模型的建立

根据管棚作用下隧道开挖过程的施工经验以及管棚既有设计方法，可采用Euler-Bernoulli梁来对管棚的受力和变形进行分析计算。根据管棚的受力特点，可将管棚沿纵向分为初期支护段、开挖段和未开挖段3部分来进行分析（见图 1）。位于洞口的导向墙由于具有独立的基础和仅承受管棚作用的特点，可认为导向墙对管棚的约束形式为固定端约束，且无初始位移。在隧道开挖过程中，开挖段在围岩压力的作用下将力传递给导向墙或初期支护以及掌子面前方围岩，管棚发生变形。将该开挖段支护后，管棚则以该变形状态固结在初期支护中，初期支护前端对管棚的约束可认为具有初始位移和初始转角的固定端约束，而此处管棚的变形值即为下个循环开挖分析时管棚起固定端处的初始变形值。在掌子面前方，可认为由于隧道的开挖，掌子面在地应力作用下发生失稳，则该段围岩将无法承受管棚变形引起的附加应力，则可认为掌子面失稳段与开挖段具有相同的性质。失稳段前方围岩承受因管棚变形而产生的附加力，该段管棚受围岩的弹性地基约束。

图  1  管棚预支护段施工示意图

Figure  1.  Construction diagram of pre-support section of pipe roof

下载: 全尺寸图片 幻灯片

在单个循环进尺的分析中，以初期支护前端作为分析模型的起点，为具有初始竖向位移${\omega _{\text{0}}}$和初始转角${\theta _{\text{0}}}$的固定端。管棚承受开挖段及掌子面失稳段上方的围岩压力。根据莫尔–库仑强度准则，可假设掌子面沿开挖方向失稳的破裂面与铅垂线的夹角为45° $- \varphi {\text{/2}}$，若单个循环开挖进尺长度为$a$，则围岩压力作用范围为$s = a + H\tan ({45^ \circ } - \varphi /2)$，其中$H$为上台阶高度。由于一个循环开挖进尺长度较短，可认为作用于管棚上的围岩压力值在该分析段保持不变。而围岩压力在纵向方向上的连续变化可离散为每个开挖区段上的常围岩压力。破裂面前方围岩以弹性地基约束形式作用于管棚，从而建立单个循环进尺的管棚分析模型（见图 2）。在管棚超前支护作用下隧道的循环开挖，其力学行为等价于单个循环进尺的管棚分析模型的不断移动，通过迭代计算可得到管棚全长范围最终的内力和变形分布。

图  2  单个循环进尺的管棚分析模型

Figure  2.  Analysis model for pipe roof with single cycle footage

下载: 全尺寸图片 幻灯片

通过对管棚受力过程的分析，针对管棚的分析模型提出以下假设：①导向墙或初期支护对管棚的约束假设为固定端约束；②导向墙处管棚无初始挠度和转角；③假设掌子面前方失稳段围岩对管棚不具有约束作用；④掌子面前方岩体对管棚的约束反力采用Pasternak弹性地基理论来确定。

围岩对管棚的弹性地基约束根据管棚长度可分为两种类型。类型Ⅰ，埋入围岩中的管棚长度足够长，该条件下可认为管棚B端的边界条件为：当$x \to \infty$时，$\omega (x) = 0$，$\theta (x) = 0$，分析图示如图 2（a）所示。类型Ⅱ，以类型Ⅰ为基本分析模型，随着在管棚预支护范围内开挖的不断推进，在对预支护段的最后一个循环进尺进行分析时，B端将发生变形。此时B端的边界条件为：当$x = s + l$时，$M(x) = 0$，$Q(x) = 0$，其中$l$为搭接长度，分析如图 2（b）所示。

管棚作为松散破碎地层中开挖隧道所采用的预支护技术，其预支护地层往往是破碎的散体介质，因此在进行管棚结构的分析计算时，采用松散围岩压力计算方法确定作用于管棚结构上的荷载比较符合实际情况。本模型浅埋隧道围岩压力的计算采用全土柱理论，深埋隧道围岩压力的计算采用普氏理论，则围岩压力可表示为统一的表达式：

式中$\gamma$为围岩重度（kN/m³）；$h$为上覆松动围岩的高度（m），浅埋隧道为取隧道埋深，深埋隧道取压力拱高度。

Pasternak弹性地基模型是在Winkler弹性地基模型的基础上假定弹簧之间存在剪切作用，引入地基剪切模量${G_{\text{p}}}$，则弹性地基作用于梁上的约束反力为

式中$k$为地基基床系数（kN/m³）；$\omega (x)$为$x$点处梁的挠度（m）；${G_{\text{p}}}$为地基剪切模量（kN/m）；${b^ * }$为地基剪切层的等效宽度（m），满足${b^ * } = b[1 + \sqrt {{G_{\text{p}}}{\text{/}}k} /b]$，$b$为管棚钢管的管径。Pasternak弹性地基模型相比Winkler地基模型通过引入地基剪切模量，弥补了Winkler模型地基不连续的缺陷，因此更符合实际情况，并且相比弹性半空间模型，Pasternak模型大大降低了计算复杂度。

2.   管棚分析模型的求解

分析Euler-Bernoulli梁在受外荷载作用发生变形时，在梁上任选微段作为分析单元，其中，剪力、弯矩及荷载的正方向均如图 3所示。

图  3  梁单元分析图示

Figure  3.  Diagram of beam element analysis

下载: 全尺寸图片 幻灯片

根据平衡条件可得

再根据梁的挠曲线近似微分方程：

将式（4），（5）代入式（3），即可得到分布荷载作用下，梁的变形控制微分方程为

单个循环进尺的管棚分析模型如图 2所示，根据荷载和约束的分布形式，可分段建立管棚的控制微分方程并求解，得到OA段控制微分方程：

式中，$q = j{q^ * }$，$j$为管棚环向间距（m），即作用于单根梁上围岩压力的宽度。其通解为

根据式（5）~（7）表示的梁的弯矩、剪力和荷载与梁挠度之间的微分关系，将挠度表达式进行求导运算，即可得到管棚转角、弯矩和剪力的分段表达式。其中OA段转角、弯矩和剪力的解析表达式分别为

根据AB段的荷载分布形式，建立控制微分方程：

令${\lambda ^{\text{4}}}{\text{ = }}\frac{{k{b^ * }}}{{4EI}}$，$\alpha = \sqrt {1 + \frac{{{G_{\text{p}}}{\lambda ^2}}}{k}}$，$\beta = \sqrt {1 - \frac{{{G_{\text{p}}}{\lambda ^2}}}{k}}$，求得其通解为

同理求得AB段转角、弯矩和剪力的解析表达式分别为

式中，${B_1}$，${B_2}$，${B_3}$，${B_4}$，${B_5}$，${B_6}$，${B_7}$，${B_8}$均为待求积分常数。

2.1   Ⅰ类管棚分析模型的求解

类型Ⅰ的模型示意图如图 2（a）所示，对于AB段梁，由于B端满足${\omega _2}(x)\left| {_{x \to \infty } = 0} \right.$，${\theta _2}(x)\left| {_{x \to \infty } = 0} \right.$的边界条件，则代入式（13）可求得${B_7}{\text{ = }}{B_{\text{8}}}{\text{ = 0}}$。由边界条件${\omega _1}(x)\left| {_{x = 0} = {\omega _0}} \right.$，${\theta _{\text{1}}}(x)\left| {_{x = 0} = {\theta _0}} \right.$，${\omega _1}(x)\left| {_{x = s} = {\omega _2}(x)} \right.\left| {_{x = s}} \right.$，${\theta _1}(x)\left| {_{x = s} = {\theta _2}(x)} \right.\left| {_{x = s}} \right.$，${M_1}(x)\left| {_{x = s} = {M_2}(x)} \right.\left| {_{x = s}} \right.$，${Q_1}(x)\left| {_{x = s}} \right. =$ ${Q_2}(x)\left| {_{x = s}} \right.$，得方程组：

式中，${\varphi _{35}} = - {{\text{e}}^{ - \alpha \lambda s}}\cos \beta \lambda s$，${\varphi _{3{\text{6}}}} = - {{\text{e}}^{ - \alpha \lambda s}}\sin \beta \lambda s$，${\varphi _{{\text{4}}5}} =$ ${\text{e}}^{-\alpha \lambda s}(\beta \lambda \mathrm{sin}\beta \lambda s+\alpha \lambda \mathrm{cos}\beta \lambda s)，{\phi }_{\text{46}}=-{\text{e}}^{-\alpha \lambda s}(\beta \lambda \mathrm{cos}\beta \lambda s-$$\alpha \lambda \mathrm{sin}\beta \lambda s)，{\phi }_{\text{5}5}=-{\text{e}}^{-\alpha \lambda s}[2\alpha \beta {\lambda }^{2}\mathrm{sin}\beta \lambda s+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\cdot$ $\mathrm{cos}\beta \lambda s]，{\phi }_{\text{56}}=-{\text{e}}^{-\alpha \lambda s}[-2\alpha \beta {\lambda }^{2}\mathrm{cos}\beta \lambda s+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\cdot$ $\mathrm{sin}\beta \lambda s]，{\phi }_{\text{6}5}=-{\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{sin}\beta \lambda s-({\alpha }^{3}-3\alpha {\beta }^{2})\cdot$${\lambda }^{3}\mathrm{cos}\beta \lambda s]，{\phi }_{\text{66}}={\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{cos}\beta \lambda s+({\alpha }^{3}-$$3\alpha {\beta ^2}){\lambda ^3}\sin \beta \lambda s]$，${\psi _3} = - q{s^4}/(24EI)$，${\psi }_{\text{4}}=-q{s}^{\text{3}}/(\text{6}EI)，$${\psi _{\text{5}}} = - q{s^{\text{2}}}/(2EI)$，${\psi _{\text{6}}} = - qs/(EI)$。

通过式（17）即可求得积分常数${B_1}$，${B_2}$，${B_3}$，${B_4}$，${B_5}$，${B_6}$，将其回代到挠度、转角、弯矩和剪力表达式中，即可求得满足类型Ⅰ分析模型边界条件下管棚的变形和受力分布。

2.2   Ⅱ类管棚分析模型的求解

类型Ⅱ的模型示意图如图 3（b）所示，由边界条件${\omega _1}(x)\left| {_{x = 0} = {\omega _0}} \right.$，${\theta _{\text{1}}}(x)\left| {_{x = 0} = {\theta _0}} \right.$，${\omega _1}(x)\left| {_{x = s} = {\omega _2}(x)} \right.\left| {_{x = s}} \right.$，${\theta _1}(x)\left| {_{x = s} = {\theta _2}(x)} \right.\left| {_{x = s}} \right.$，${M_1}(x)\left| {_{x = s} = {M_2}(x)} \right.\left| {_{x = s}} \right.$，${Q_1}(x)\left| {_{x = s} = } \right.$ ${Q_2}(x)\left| {_{x = s}} \right.$，${M_{\text{2}}}(x)\left| {_{x = s + l} = } \right.0$，${Q_2}(x)\left| {_{x = s + l}} \right. = 0$，得方程组：

式中，${\varphi _{35}} = - {{\text{e}}^{ - \alpha \lambda s}}\cos \beta \lambda s$，${\phi }_{3\text{6}}=-{\text{e}}^{-\alpha \lambda s}\mathrm{sin}\beta \lambda s，{\phi }_{37}=$$-{\mathrm{e}}^{\alpha \lambda s}\mathrm{cos}\beta \lambda s，$${\phi }_{3\text{8}}=-{e}^{\alpha \lambda s}\mathrm{sin}\beta \lambda s，{\phi }_{\text{4}5}={\text{e}}^{-\alpha \lambda s}(\beta \lambda \mathrm{sin}\beta \cdot$ $\lambda s + \alpha \lambda \cos \beta \lambda s)$，${\varphi _{{\text{46}}}} = - {{\text{e}}^{ - \alpha \lambda s}}(\beta \lambda \cos \beta \lambda s - \alpha \lambda \sin \beta \lambda s)$，${\varphi _{{\text{47}}}} = - {{\text{e}}^{\alpha \lambda s}}(\alpha \lambda \cos \beta \lambda s - \beta \lambda \sin \beta \lambda s)$，${\varphi _{{\text{48}}}} = - {{\text{e}}^{\alpha \lambda s}} \cdot (\beta \lambda$ $\mathrm{cos}\beta \lambda s\text{+}\alpha \lambda \mathrm{sin}\beta \lambda s)，{\phi }_{\text{5}5}=-{\text{e}}^{-\alpha \lambda s}[2\alpha \beta {\lambda }^{2}\mathrm{sin}\beta \lambda s+({\alpha }^{2}-$${\beta }^{2}){\lambda }^{2}\mathrm{cos}\beta \lambda s]，{\phi }_{\text{56}}=-{\text{e}}^{-\alpha \lambda s}[-2\alpha \beta {\lambda }^{2}\mathrm{cos}\beta \lambda s+({\alpha }^{2}-{\beta }^{2})$${\lambda ^2}\sin \beta \lambda s]$，${\varphi _{{\text{57}}}} = - {{\text{e}}^{\alpha \lambda s}}[ - 2\alpha \beta {\lambda ^2}\sin \beta \lambda s + ({\alpha ^2} - {\beta ^2}){\lambda ^2}$ $\mathrm{cos}\beta \lambda s]，$${\varphi _{{\text{58}}}} = - {{\text{e}}^{\alpha \lambda s}}[2\alpha \beta {\lambda ^2}\cos \beta \lambda s + ({\alpha ^2} - {\beta ^2}){\lambda ^2} \cdot \sin \beta$$\lambda s]，{\phi }_{\text{6}5}=-{\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{sin}\beta \lambda s-({\alpha }^{3}-3\alpha {\beta }^{2})\cdot {\lambda }^{3}$$\mathrm{cos}\beta \lambda s]，{\phi }_{\text{66}}={\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{cos}\beta \lambda s+({\alpha }^{3}-3\alpha {\beta }^{2})\cdot {\lambda }^{3}$$\mathrm{sin}\beta \lambda s]，{\phi }_{\text{67}}=-{\text{e}}^{\alpha \lambda s}[({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{cos}\beta \lambda s\text{+}({\beta }^{3}-3{\alpha }^{2}\beta )\cdot$${\lambda }^{3}\mathrm{sin}\beta \lambda s]，{\phi }_{\text{68}}=-{\text{e}}^{\alpha \lambda s}[({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{sin}\beta \lambda s-({\beta }^{3}-3{\alpha }^{2}\beta )\cdot$${\lambda ^3}\cos \beta \lambda s]$，${\varphi _{{\text{7}}5}} = {{\text{e}}^{ - \alpha \lambda (s + l)}}[2\alpha \beta {\lambda ^2}\sin \beta \lambda (s + l) + ({\alpha ^2} -$${\beta ^2})$${\lambda ^2}\cos \beta \lambda (s + l)]$，${\varphi _{{\text{76}}}} = {{\text{e}}^{ - \alpha \lambda (s + l)}}[ - 2\alpha \beta {\lambda ^2}\cos \beta \lambda (s +$ $l)+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\mathrm{sin}\beta \lambda (s+l)]，{\phi }_{\text{77}}={\text{e}}^{\alpha \lambda (s+l)}[-2\alpha \beta {\lambda }^{2}\cdot$$\sin \beta \lambda (s + l) + ({\alpha ^2} - {\beta ^2}){\lambda ^2}\cos \beta \lambda (s + l)]$，${\varphi _{{\text{78}}}} = {{\text{e}}^{\alpha \lambda (s + l)}} \cdot$$[2\alpha \beta {\lambda }^{2}\mathrm{cos}\beta \lambda (s+l)+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\mathrm{sin}\beta \lambda (s+l)]，{\phi }_{\text{8}5}={\text{e}}^{-\alpha \lambda (s+l)}\cdot$$[({\beta ^3} - 3{\alpha ^2}\beta ){\lambda ^3}\sin \beta \lambda (s + l) - ({\alpha ^3} - 3\alpha {\beta ^2}) \cdot {\lambda ^3}\cos \beta \lambda (s +$$l)]，{\phi }_{\text{86}}={\text{e}}^{-\alpha \lambda (s+l)}[(3{\alpha }^{2}\beta -{\beta }^{3}){\lambda }^{3}\mathrm{cos}\beta \lambda (s+l)-({\alpha }^{3}-3\alpha {\beta }^{2})\cdot$${\lambda }^{3}\mathrm{sin}\beta \lambda (s+l)]，{\phi }_{\text{87}}={\text{e}}^{\alpha \lambda (s+l)}[({\beta }^{3}-3{\alpha }^{2}\beta )\cdot {\lambda }^{3}\mathrm{sin}\beta \lambda (s+$$l)+({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{cos}\beta \lambda (s+l)]，{\phi }_{\text{88}}={\text{e}}^{\alpha \lambda (s+l)}\cdot [(3{\alpha }^{2}\beta -$${\beta }^{3}){\lambda }^{3}\mathrm{cos}\beta \lambda (s+l)+({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{sin}\beta \lambda (s+l)]，{\psi }_{3}=$$- q{s^4}/(24EI)$，${\psi _{\text{4}}} = - q{s^{\text{3}}}/({\text{6}}EI)$，${\psi _{\text{5}}} = - q{s^{\text{2}}}/(2EI)$，${\psi _{\text{6}}} =$$- qs{\text{/}}(EI)$。

同理，通过式（18）即可求得积分常数${B_1}$，${B_2}$，${B_3}$，${B_4}$，${B_5}$，${B_6}$，${B_7}$，${B_8}$，将其回代到挠度、转角、弯矩和剪力表达式中，即可求得满足类型Ⅱ分析模型边界条件下管棚的变形和受力分布。

3.   案例计算

南宁至崇左铁路引入机场隧道起讫里程DK25+615—DK28+935，全长3320 m。其中DK25+940—DK27+234段长度1294 m为暗挖区间，埋深3.4~15.2 m，施工工法为CRD法（350 m）和三台阶临时仰拱（944 m）；该隧道暗挖区间存在浅埋、开挖下穿既有道路、临近既有构（建）筑物。隧址内岩层上部为杂填土，下部为粉质黏土、泥质砂岩夹砂质泥岩、砾岩，弱风化至全风化，岩溶发育，层面起伏变化大。本隧道进口明暗分界处拱部设置一环35 m长Φ108 mm×6 mm管棚加强支护，管棚环向间距40 cm，具体参数见表 1。

表  1  设计参数取值

Table  1.  Values of design parameters

基床系数/(kN·m-3)	地基剪切模量/(kN·m-1)	隧道埋深/m	内摩擦角/(°)	围岩重度/(kN·m-3)	等效惯性矩/m4	等效弹性模量/kPa	开挖高度/m	管棚间距/m	进尺深度/m
30000	2800	6	45.6	24	6.68×10-6	7.89×107	2.4	0.4	0.6

下载: 导出CSV 

| 显示表格

本文以隧道进口明暗分界处拱顶单根钢管为研究对象，则该钢管的挠曲变形与隧道因开挖引起的拱顶沉降值相同。对管棚分析模型中的参数分别确定取值，并利用MATLAB编程进行应用分析。第一个循环进尺下管棚挠度、转角、弯矩和剪力在开挖方向上的分布曲线如图 4所示。

图  4  管棚的变形及内力分布

Figure  4.  Distribution of deformation and internal force of pipe roof

下载: 全尺寸图片 幻灯片

图 4表明，管棚挠度在OA段逐渐增大，且在A点处达到最大值7.2 mm。AB段挠度逐渐减小，当计算深度为4 m时，该循环开挖对该处管棚的变形无影响。随计算深度的加深，弯矩呈震荡形式趋近于零，其中最大值位于固定端处为-23.5 kN·m。剪力的分布同样随计算深度的增大而收敛于零，最大值同样位于固定端处，大小为62.6 kN。当计算深度达到4 m时，管棚的弯矩和剪力将趋近于零，由此可得一个循环开挖区段对管棚的影响范围长度为4 m。

通过以上单个开挖区段的分析，得到管棚在该分析区段的变形和受力特征。该区段前0.6 m的变形与受力将随初期支护的完成而与初期支护固结在一起，0.6 m处的变形值将成为下一个分析区段的位移初值。通过迭代计算，可得到管棚预支护范围内30 m全部开挖支护完成后隧道拱顶下沉分布如图 5所示，管棚全长范围内内力的分布情况如图 6所示。

图  5  管棚预支护范围内拱顶下沉纵向分布

Figure  5.  Longitudinal distribution of vault settlement within scope of pre-support of pipe roof

下载: 全尺寸图片 幻灯片

图  6  管棚全长范围的内力分布

Figure  6.  Distribution of internal force in whole length of pipe roof

下载: 全尺寸图片 幻灯片

图 5表明，管棚的挠曲变形随掌子面的不断推进而逐渐增大。当拱顶沉降达到22.8 mm时，随掌子面的进一步推进，拱顶沉降将基本保持不变。

王道远等^[12]通过现场试验得出，管棚的变形整体呈凹形分布，在开挖初期，管棚的变形速率较快。而图中管棚变形分布特征符合这一规律，且管棚变形的数量级基本一致，从而证明了该分析模型的合理性和实用性。

图 6表明，管棚的弯矩和剪力的最不利截面均发生在导向墙的约束处，其中弯矩最大值为23.5 kN·m，剪力最大值为62.6 kN。管棚全长范围内弯矩基本保持在12 kN·m，剪力基本保持在52 kN。

综上所述，利用该分析模型可以系统地对管棚的支护效果和受力特征进行预测。通过对算例的数值计算，得到管棚预结构的受力和变形分布特征，从而体现了管棚预支护体系的主要特征。

4.   参数分析

本节在前述理论研究的基础上，通过调整管棚的环向间距、管径、开挖进尺长度以及搭接长度的设计参数，研究各设计参数对管棚受力和变形的影响形式和程度，从而指导管棚的设计和施工。

4.1   环向间距对管棚受力和变形的影响分析

取其他参数不变，当管棚环向间距j为0.40，0.60，0.80，1.00 m时，分别计算并绘出管棚的变形分布曲线（见图 7）和内力分布曲线（见图 8），研究环向间距对管棚变形和受力的影响。

图  7  环向间距对管棚变形的影响

Figure  7.  Influences of circumferential spacing on deformation of pipe roof

下载: 全尺寸图片 幻灯片

图  8  环向间距对管棚受力的影响

Figure  8.  Influence of circumferential spacing on stress of pipe roof

下载: 全尺寸图片 幻灯片

图 7表明，管棚环向间距的变化不改变其变形分布形式，最大挠度${\omega _{\max }}{\text{(mm)}}$和最大转角${\theta _{\max }}$(°)发生的位置不变，其中最大挠度和最大转角均与环向间距之间满足线性增长关系，关系式分别为

图 8表明，管棚环向间距同样不改变管棚的内力分布形式，最不利截面不随环向间距的变化而变化，管棚的最大弯矩${M_{\max }}{\text{(kN}} \cdot {\text{m)}}$和最大剪力${Q_{\max }}{\text{(kN)}}$均与环向间距之间满足线性增长关系，关系式分别为

4.2   管径对管棚受力和变形的影响分析

取其他参数不变，当管径$b$分别为0.089，0.108，0.152，0.180 m时，分别计算并绘制出管棚的挠度分布曲线（见图 9）和弯矩分布曲线（见图 10），研究管径对管棚变形和受力的影响。当管棚的管径发生变化时，管棚的等效截面惯性矩和等效弹性模量均发生变化，根据刚度等效法^[14]求得不同管径下的等效截面惯性矩和等效弹性模量见表 2。

图  9  管径对管棚挠度的影响

Figure  9.  Influence of pipe diameter on deflection of pipe roof

下载: 全尺寸图片 幻灯片

图  10  管径对管棚弯矩的影响

Figure  10.  Influences of pipe diameter on bending moment of pipe roof

下载: 全尺寸图片 幻灯片

表  2  不同管径的截面计算参数

Table  2.  Cross section calculation parameters of different pipe diameters

管径/m	等效惯性矩/(10-6 m4)	等效弹性模量/(107 kPa)
0.089	3.08	9.24
0.108	6.68	7.89
0.152	26.00	5.89
0.180	51.50	5.07

下载: 导出CSV 

| 显示表格

从图 9可以看出，管径越粗，管棚的挠曲变形越平缓，随着管径的增大，最大挠度${\omega _{\max }}$逐渐减小，并且随着管径的不断增大对管棚挠曲变形的影响将越来越小。最大挠度${\omega _{\max }}$与管径b(m)之间满足

从图 10可以看出，随着管径的增大，弯矩的危险截面位置不变，管棚所受的最大弯矩${M_{\max }}$（kN·m）也随着增大。最大弯矩与管径之间满足

4.3   开挖进尺长度对管棚受力和变形的影响分析

取其他参数不变，当开挖进尺长度$a$分别为0.5，0.75，1，1.25，1.5 m时，分别绘制出管棚变形分布曲线（图 11）和内力分布曲线（图 12），研究开挖进尺长度对管棚变形和受力的影响。

图  11  开挖进尺长度对管棚变形的影响

Figure  11.  Influences of excavation footage length on deformation of pipe roof

下载: 全尺寸图片 幻灯片

图  12  开挖进尺长度对管棚受力的影响

Figure  12.  Influences of excavation footage length on stress of pipe roof

下载: 全尺寸图片 幻灯片

图 11表明，随着开挖进尺长度的增大，管棚最大变形位置不断前移，且最大值不断增大。通过开挖进尺长度与最大变形关系曲线图可得到，开挖进尺长度与最大挠度${\omega _{\max }}$和最大转角${\theta _{\max }}$之间呈非线性增长关系，进尺越长，变形增长速度越快，关系式分别为

说明开挖进尺长度是控制管棚变形的关键施工技术指标。

图 12表明，开挖进尺长度的增大同样会使管棚的内力分布前移，但最不利截面总位于固定端处。最大弯矩${M_{\max }}$与开挖进尺长度之间满足非线性关系式：

式中，开挖进尺长度$a$（m），${M_{\max }}$最大弯矩（kN·m）。而最大剪力${Q_{\max }}$（kN）与开挖进尺长度之间满足线性关系式：

4.4   搭接长度对管棚受力和变形的影响分析

由类型Ⅱ分析模型可计算管棚搭接长度对管棚作用效果的影响，搭接长度$D = H\tan ({45^ \circ } - \varphi /$$2) + l$，分别取搭接长度为1.1，1.3，1.5，5 m，分别绘制出管棚挠度分布曲线以及搭接长度与最大挠度关系曲线（图 13）和弯矩分布曲线以及搭接长度与最大弯矩关系曲线（图 14），研究搭接长度对管棚变形和受力的影响。

图  13  搭接长度对管棚挠度的影响

Figure  13.  Influences of lap length on deflection of pipe roof

下载: 全尺寸图片 幻灯片

图  14  搭接长度对管棚弯矩的影响

Figure  14.  Influences of lap length on bending moment of pipe roof

下载: 全尺寸图片 幻灯片

图 13表明，随着搭接长度的增大，管棚发生最大挠度的位置和挠度分布形式基本不变。当搭接长度大于2.5 m时，最大挠度将趋近于7.3 mm而保持不变。

图 14表明，随着搭接长度的增大，管棚的最不利截面位置和弯矩的分布形式基本不变。当搭接长度大于2.5 m时，最大弯矩将趋近于23.5 kN·m而保持不变。综上所述，在该工程条件下，当管棚搭接长度大于2.5 m时，搭接长度的变化将对管棚的变形和受力基本没有影响。

5.   结论

（1）通过建立和求解基于Pasternak弹性地基理论的单个循环进尺下管棚分析模型可得管棚的受力和变形的分布特征。考虑隧道开挖支护动态施工过程，通过迭代计算可有效预测管棚预支护段开挖支护完成后管棚整体的受力和变形分布，且其理论值与已有文献实测结果具有相同的分布特征，从而证明了模型的合理性。

（2）管棚环向间距的变化不改变管棚内力和变形的分布形式，且最不利截面位置保持不变，管棚的受力和变形与环向间距之间呈线性增长关系，随环向间距的增大而增大。

（3）选用的管径越粗，管棚的挠曲变形越平缓，挠曲变形越小，并且随着管径的不断增大，管径对管棚挠曲变形的影响将越来越小。随着管径的增大，弯矩的危险截面位置不变，管棚的最大弯矩与管径之间呈非线性增长关系。

（4）随着开挖进尺长度的增大，管棚最大变形位置不断前移，最大变形值与开挖进尺长度之间呈非线性增长关系，开挖进尺越长，变形增长速度越快。说明开挖进尺长度是控制管棚变形的关键施工技术指标。开挖进尺长度越大，最大弯矩和最大剪力越大。

（5）管棚的最大内力及变形随搭接长度的增长而迅速趋于收敛。在本文案例工况下，当管棚搭接长度大于2.5 m时，搭接长度的变化将对管棚的变形和受力基本没有影响。