[发明专利]一种新建隧道下穿引起既有盾构隧道位移计算方法

摘要

本发明涉及一种新建隧道下穿引起既有盾构隧道位移计算方法，包括如下步骤一、新建隧道下穿引起的既有隧道附加应力计算刀盘附加推力q产生的竖向附加应力σz‑q、盾壳摩擦力f产生的竖向附加应力σz‑f、注浆附加压力p产生的竖向附加应力σz‑p1和σz‑p2、土体损失产生的附加应力；二、基于最小势能原理计算隧道纵向变形盾构隧道的总势能、假设隧道衬砌环的位移函数、变分控制方程。本发明的有益效果是本发明考虑了隧道环间的“接头效应”，充分考虑了隧道的结构特性和物理特性，运用Matlab可以快速计算得到由于新建隧道下穿引起的既有盾构隧道纵向位移，由此可判断新建隧道下穿作用下既有隧道结构的安全性，与盾构隧道实际的受力和变形较为符合。

主权项

一种新建隧道下穿引起既有盾构隧道位移计算方法，其特征在于，包括如下步骤：一、新建隧道下穿引起的既有隧道附加应力计算1.1各个分力产生的附加应力设既有盾构隧道轴线处一点坐标为(x，y，z)，求该点处的附加应力；根据Mindlin公式进行积分求解，可推导得到刀盘附加推力q、盾壳摩擦力f、注浆附加压力p在既有隧道轴线处的附加应力；1.1.1刀盘附加推力q产生的竖向附加应力σz‑q在盾构刀盘处取任一微元体dA＝rdrdθ，运用Mindlin公式积分得到刀盘附加推力q在该点处产生的竖向附加应力σz‑q为：dσz-q=-qrdrdθy8π(1-v){(1-2v)R13-(1-2v)R23-3(z-z0+rsinθ)2R15-3(3-4v)(z+z0-rsinθ)2R25+6(z0-rsinθ)R25[z0-rsinθ+(1-2v)·(1-2v)(z+z0-rsinθ)+5z(z+z0-rsinθ)2R22]}]]>R1=(x-rcosθ)2+y2+(z-z0+rsinθ)2R2=(x-rcosθ)2+y2+(z+z0-rsinθ)2]]>式中：v为土的泊松比；1.1.2盾壳摩擦力f产生的竖向附加应力σz‑f在盾壳壁上取任一微元体dA＝Rsdsdθ，运用Mindlin公式积分得到盾壳摩擦力f在该点处产生的竖向附加应力σz‑f为：dσz-f=-fRsdsdθ(y-s)8π(1-v){(1-2v)R33-(1-2v)R23-3(z-z0+Rssinθ)2R35-3(3-4v)(z+z0-Rssinθ)2R45+6(z0-Rssinθ)R45[z0-Rssinθ+(1-2v)·(1-2v)(z+z0-Rssinθ)+5z(z+z0-Rssinθ)2R42]}]]>R3=(x-Rscosθ)2+(y-s)2+(z-z0+Rssinθ)2R2=(x-rcosθ)2+(y-s)2+(z+z0-Rssinθ)2]]>式中：L为新建盾构机机长；1.1.3注浆附加压力p产生的竖向附加应力σz‑p1和σz‑p2在盾尾注浆处取任一微元体dA＝Rsdsdθ，运用Mindlin公式积分得到：(1)注浆附加压力竖向分力p1dσz-p1=-pRssinθdsdθ8π(1-v)[-(1-2v)(z-z0+Rssinθ)R53+(1-2v)(z-z0+Rssinθ)R63-3(z-z0+Rssinθ)3R55-3(3-4v)z(z-z0+Rssinθ)2R65-3(z0-Rssinθ)(z+z0-Rssinθ)(5z-z0+Rssinθ)R65-30(z0-Rssinθ)z(z+z0-Rssinθ)2R67]]]>R5=(x-Rscosθ)2+(y-L-s)2+(z-z0+Rssinθ)2R6=(x-Rscosθ)2+(y-L-s)2+(z+z0-Rssinθ)2]]>(2)注浆附加压力水平向分力p2dσz-p2=pRscosθdsdθ(x-Rscosθ)8π(1-v){(1-2v)R73-(1-2v)R83-3(z-z0+Rssinθ)2R75-3(3-4v)z(z+z0-Rssinθ)2R85+6(z0-Rssinθ)R85[z0-Rssinθ+(1-2v)(z+z0-Rssinθ)+5z(z+z0-Rssinθ)2R82]}]]>R7=(x-Rscosθ)2+(y-L-s)2+(z-z0+Rssinθ)2R8=(x-Rscosθ)2+(y-L-s)2+(z+z0-Rssinθ)2]]>1.1.4土体损失产生的附加应力得到单线盾构施工由于土体损失引起的管线平面处任一点的土体竖向位移为Uz，进而得到土体损失在既有隧道轴线处产生的附加应力σs为：1.其中，λ=14-gπRsη(y)[arcsin(dRs-Rs(1-η(y)))+1-(dRs-Rs(1-η(y)))2-1]]]>式中：k为地基基床系数；d为土体移动焦点到新建隧道中心点的距离；η为最大土体损失百分率(％)；最终得到隧道下穿引起既有隧道轴线处产生的竖向总附加应力值σz为：σZ＝σz‑q+σz‑f+σz‑p1+σz‑p2+σs  (7)二、基于最小势能原理计算隧道纵向变形在分析盾构隧道与土体相互作用时，假定：将盾构隧道衬砌环视为由剪切弹簧连接的弹性地基短梁，新建隧道下穿导致隧道以环间剪切错台的方式进行变形；2.1运用能量变分法计算隧道纵向位移量2.1.1盾构隧道的总势能任取盾构隧道一环进行分析，编号为m，其所受到的竖向荷载Fz为：Fz＝Pz(x)‑kDSz(x)‑kt(ΔWz(m+1)+ΔWz(m))  (8)式中：Pz(x)＝Dσz，Pz(x)为附加水平荷载；kDSz(x)为地基抗力，D为既有盾构隧道直径，k为地基基床系数，采用Vesic公式计算，Sz(x)为地基弹簧的位移，根据位移协调条件则Sz(x)＝Wz(x)，这里Wz(x)为隧道的水平位移，Es为地基土的弹性模量，EtIt为隧道的等效抗弯刚度；kt为隧道的环间剪切刚度；kt(△Wz(m+1)+△Wz(m))为水平环间剪切力；根据盾构隧道每片衬砌环的受荷状况，分析计算得到盾构隧道的总势能，具体分为以下三部分：①新建隧道下穿引起临近盾构隧道轴线处附加荷载做功Wp；②盾构隧道衬砌环克服地层抗力做功Wk；③衬砌环克服盾构环间剪切力做功Ws；可得到新建隧道下穿引起的临近地铁盾构隧道的总势能为：Ep＝Wp+Wk+Ws  (9)2.1.2假设隧道衬砌环的位移函数能量变分解法原理是假定合适的位移函数来表示盾构隧道受新建隧道影响的基本变形形状；假设隧道位移函数如下，并按傅里叶级数展开,隧道的竖向位移函数为：式中：δ为既有盾构隧道环宽；A＝{a0 a1 ... an}T，A为位移函数中的待定系数矩阵；n为傅里叶级数的展开阶数；2.1.3变分控制方程基于能量变分法，将总势能Ep对各待定系数取极值，即：式中：ξi为矩阵A中的各个元素；对上式求解，可得隧道竖向位移控制方程为：式中：2N为新建隧道施工既有隧道受影响的隧道衬砌环数；将上式表达为矩阵形式：([Kt]+[Ks]){A}＝{Pz}T  (13)式中：[Kt]为隧道环间刚度矩阵，[Ks]为土体刚度矩阵，其中：{Pz}T表示自由土体位移和隧道衬砌环的相作用效应，具体表示为：由式(12)计算可得到待定系数矩阵A，再代入假设的隧道位移函数W(x)即式(10)，可以得到在新建隧道开挖时引起的既有盾构隧道纵向位移值；相邻盾构管片之间位移差值即错台量△W，其中隧道的竖向错台量为：ΔWz＝Wz[(m+1)δ]‑Wz(m)  (15)。