[0004] 本发明的目的是克服现有技术中不足,提供一种地面堆载引起既有盾构隧道位移计算方法。
[0005] 为了达到上述目的,本发明是通过以下技术方案实现的:
[0006] 本专利利用“剪切错台模型”,将盾构隧道等效为一个由剪切弹簧连接的弹性地基短梁,运用最小势能原理计算地面堆载引起下卧盾构隧道的纵向变形,以及相邻盾构隧道环之间的错台量和环间剪切力的大小,可评估盾构隧道在运营期间的结构安全性和抗渗性。本发明提供一种地面堆载引起既有盾构隧道位移计算方法,其特征在于,先建立地面堆载力学计算模型,如图1和图2所示;令荷载大小为q,单位为kPa;堆载长为L,单位为m;堆载宽为B,单位为m;盾构隧道轴线处埋深为H,单位为m;盾构隧道上部覆土厚h,单位为m;盾构隧道外径为D,单位为m;堆载中心离盾构隧道轴线的水平距离为s,单位为m;以地面堆载中心为原点,取平行盾构隧道方向为y轴,垂直盾构隧道方向为x轴;具体包括如下步骤:
[0007] 步骤1):地面堆载引起的盾构隧道附加应力计算
[0008] 利用Boussinesq公式,在堆载区域取一微单元,微单元的长为ξ,宽为η,作用在微单元上的应力为σ,其单位力表达为σdξdη,通过积分可以求得引起盾构隧道轴线上某一点(x1,y1,z1)的土体附加水平应力σx为:
[0009]
[0010] 盾构隧道轴线上某一点(x1,y1,z1)的土体附加竖向应力σz为:
[0011]
[0012] 式中:
[0013] 步骤2):基于最小势能原理计算盾构隧道纵向变形
[0014] 在分析盾构隧道与土体相互作用时,假定:将盾构隧道衬砌环视为由剪切弹簧连接的弹性地基短梁,地面堆载导致盾构隧道以环间剪切错台的方式进行变形;
[0015] 步骤2.1):运用能量变分法计算盾构隧道纵向位移量
[0016] 步骤2.1.1):盾构隧道的总势能
[0017] 任取盾构隧道一环进行分析,编号为m,其所受到的水平荷载Fx为:
[0018] Fx=Px(y)-kDSx(y)-Kt(ΔWx(m+1)+ΔWx(m)) (3)
[0019] 式中:ΔWx(m+1)为m+1环管片的位移量,ΔWx(m)为m环管片的位移量,[0020] Px(y)=Dσx,Px(y)为附加水平荷载,σx为附加水平应力;kDSx(y)为地基抗力,k为地基基床系数,采用Vesic公式计算, Sx(y)为地基弹簧的位移,根据位移协调条件则Sx(y)=Wx(y),这里Wx(y)为盾构隧道的水平位移,Es为地基土的弹性模量,EtIt为盾构隧道的等效抗弯刚度;kt为盾构隧道的环间剪切刚度;
[0021] kt(ΔWx(m+1)+ΔWx(m))为水平环间剪切力;
[0022] 同理可得到盾构隧道的竖向荷载Fz为:
[0023] FZ=PZ(y)-kDSZ(y)-Kt(ΔWZ(m+1)+ΔWZ(m)) (4)
[0024] 根据盾构隧道每片衬砌环的受荷状况,分析计算得到盾构隧道的总势能,具体分为以下三部分:①地面堆载引起邻近盾构隧道轴线处附加荷载做功Wp;②盾构隧道衬砌环克服地层抗力做功Wk;③衬砌环克服盾构环间剪切力做功Ws;
[0025] 可得到地面堆载引起的临近地铁盾构隧道的总势能为:
[0026] Ep=Wp+Wk+Ws (5)
[0027] 步骤2.1.2):假设盾构隧道衬砌环的位移函数
[0028] 能量变分解法原理是假定合适的位移函数来表示盾构隧道受地面堆载影响的基本变形形状;假设盾构隧道位移函数如下,并按傅里叶级数展开;盾构隧道的水平位移函数为:
[0029]
[0030] 盾构隧道的竖向位移函数为:
[0031]
[0032] 式中:
[0033] δ为盾构隧道环宽,A={a0 a1 ... an}T,A为位移函数中的待定系数矩阵,n为傅里叶级数的展开阶数;
[0034] 步骤2.1.3):变分控制方程
[0035] 基于能量变分法,将总势能Ep对各待定系数取极值,即:
[0036]
[0037] 式中:ξi为矩阵A中的各个元素;
[0038] 对上式求解,可得盾构隧道水平位移控制方程为,式中:ξi为矩阵A中的各项元素:
[0039]
[0040] 同理,可得盾构隧道竖向位移的控制方程为:
[0041]
[0042] 将上式表达为矩阵形式:
[0043] ([Kt]+[Ks]){A}={Pl}T (11)
[0044] 式中:[Kt]为盾构隧道环间刚度矩阵,
[0045]
[0046] [Ks]为土体刚度矩阵,
[0047] 其中:{Pl}T表示自由土体位移和盾构隧道衬砌环的相作用效应(l=x;z),具体表示为:
[0048]
[0049]
[0050] 由式(11)计算可得到待定系数矩阵A,再代入假设的盾构隧道位移函数W(y)即式(6)、(7),可以得到在地面堆载作用下引起的盾构隧道纵向位移值;
[0051] 相邻盾构管片之间位移差值即错台量△W,其中盾构隧道的水平错台量为:
[0052] ΔWx=Wx[(m+1)δ]-Wx(m) (14)
[0053] 盾构隧道的竖向错台量为:
[0054] ΔWz=Wz[(m+1)δ]-Wz(m) (15)
[0055] 相邻盾构管片环之间的水平剪切力为:
[0056] Qx={Wx[(m+1)δ]-Wx(mδ)}·kt (16)
[0057] 盾构隧道的竖向剪切力为:
[0058] Qz={Wz[(m+1)δ]-Wz(mδ)}·kt (17)
[0059] 取10阶的刚度矩阵[Ks]和[Kt]即可满足计算精度,以上算法通过Matlab编程进行数值计算;
[0060] 上述得到由地面堆载引起的盾构隧道纵向位移,以及盾构隧道环之间的错台量、水平剪切力和竖向剪切力,由此判断地面堆载作用下隧道结构的安全性。
[0061] 与现有技术相比,本发明的有益效果如下:
[0062] 本专利考虑了隧道环间的“接头效应”,将隧道等效为一个由剪切弹簧连接的弹性地基短梁,充分考虑了隧道的结构特性和物理特性,运用Matlab可以快速计算得到由于地面堆载引起的盾构隧道纵向位移(包括竖向位移和水平位移),以及盾构隧道环之间的错台量(包括水平错台量和竖向错台量)和剪切力(包括水平剪切力和竖向剪切力),由此可判断地面堆载作用下隧道结构的安全性,与盾构隧道实际的受力和变形较为符合。
[0063] 计算得到的隧道位移值,环间错台量和剪切力值与相关国标或者地方标准进行对比,例如一般隧道位移限制为20mm,如果计算得到的隧道位移值大于20mm,就需要加强监管或者在进行堆载作业时对隧道进行加固;隧道错台量以上海地方标准为例,以4mm为限制,若大于4mm,需要加强监管或者在进行堆载作业时对隧道进行加固;对于剪力值来说,以一般隧道环间有17颗螺栓为例,隧道环间剪切强度极限为665.36kN,若大于665.36kN,需要加强监管或者在进行堆载作业时对隧道进行加固。