一种双圆盾构隧道施工偏转引起土体变形的计算方法技术

本发明专利技术涉及一种双圆盾构隧道施工偏转引起土体变形的计算方法，包括以下步骤：1)绘制双圆盾构隧道开挖的二维断面图；2)确定隧道的收敛半径；3)绘出隧道收敛的形状、大小和位置；4)建立整体直角坐标系；5)建立局部极坐标系；6)建立单位土体损失引起土体变形的函数关系；7)将土体损失分布区域划分为4个子域的代数和；8)计算总的土体变形分量；9)根据土体变形分量确定土体变形量及方向。本方法操作方便，通过简单的设计参数和程序计算就能准确得出双圆盾构隧道施工引起地表及地下土体任一位置的竖向和横向变形，从而可以模拟整个横断面的土体变形场。本方法不仅能计算正常工况下双圆盾构隧道施工引起的土体变形，而且能计算任意偏转工况下的土体变形，并可编写成简单的程序软件，适合推广应用。

【技术实现步骤摘要】
一种双圆盾构隧道施工偏转引起土体变形的计算方法
本专利技术属于隧道工程中的盾构隧道施工
，涉及一种双圆盾构隧道施工偏转引起土体变形的计算方法，适用于软土地层中浅埋双圆盾构隧道施工引起土体变形的预测。
技术介绍
在地铁隧道、地下管网等城市浅埋隧道的施工中，施工引起地表及地下土体变形的预测和控制是一个极为重要的课题。对于传统的单圆隧道施工，其土体变形的预测通常采用经验统计方法，即假定土体沉降曲线为一倒置的正态分布曲线，由最大沉降值和沉降槽宽度系数确定。作为一种新的隧道施工技术，双圆盾构隧道具有施工效率高、工程造价低、节约地下空间资源、对周边环境影响小等诸多优势，而且不存在隧道间横向隧道的施工，避免了一些潜在的危险。但双圆盾构隧道的断面复杂，传统的经验统计方法不再适用于其施工引起的土体变形计算，而且施工中盾构偏转无法避免，由此导致的土体变形计算更加复杂。目前也有采用随机介质理论预测土体变形的计算方法，但这类方法不能计算隧道顶部所在深度以下的土体变形，得到的土体变形场不完整。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种双圆盾构隧道施工偏转引起土体变形的计算方法，该方法基于弹性介质理论，结合镜像法和坐标转换等，能够通过简单的设计参数和程序计算就能准确得出双圆盾构隧道施工引起地表及地下土体任一位置的竖向和横向变形，从而可以模拟整个横断面的土体变形场。为达到上述目的，本专利技术提供如下技术方案：一种双圆盾构隧道施工偏转引起土体变形的计算方法，包括以下步骤：1)绘制双圆盾构隧道开挖形状、大小和位置的二维断面图；2)通过土体损失率确定隧道的收敛半径R2：R2＝R1(1-ε)0.5式中：R1代表隧道的开挖半径，ε为土体损失率；3)结合收敛半径和双圆盾构施工偏转角，绘出隧道收敛的形状、大小和位置；4)以双圆盾构隧道中心正上方的地表位置为坐标原点,垂直于隧道轴线方向为X轴,竖直方向为Z轴,建立整体直角坐标系；5)以任一单位开挖土体所在圆心为坐标原点，X方向为r轴，逆时针方向为θ轴，建立局部极坐标系；6)根据弹性介质理论，采用镜像法和坐标转换，建立单位土体损失引起土体变形的函数关系：竖向变形：横向变形：式中：(x,z)为土体变形的目标点坐标，(r,θ)和(X0,Z0)为单位土体损失所在的位置极坐标及其所在圆心坐标，ν为土体泊松比；7)根据隧道开挖和收敛确定的土体损失分布，将其划分为4个子域的代数和：S＝SI+SII-SIII-SIV；8)通过对4个子域分别求积分计算总的土体变形分量：总的竖向变形：总的横向变形：式中：α代表盾构偏转角，t代表双圆隧道圆心距之半，h代表双圆隧道中心的埋深；9)根据土体变形分量确定土体变形量及方向：变形量：变形方向：本专利技术的有益效果在于：本专利技术所述的方法使用操作方便，通过简单的设计参数和程序计算就能准确得出双圆盾构隧道施工引起地表及地下土体任一位置的竖向和横向变形，从而可以模拟整个横断面的土体变形场。本专利技术不仅能计算正常工况下双圆盾构隧道施工引起的土体变形，而且能计算任意偏转工况下的土体变形，操作简便，计算准确，并可编写成简单的程序软件，适合推广应用。附图说明为了使本专利技术的目的、技术方案和有益效果更加清楚，本专利技术提供如下附图进行说明：图1为双圆盾构隧道施工偏转引起土体变形的计算模型图；图2为土体损失分布的积分子域区划图：(a)沿隧道开挖边界积分子域区划；(b)沿隧道收敛边界积分子域区划；图3为正常工况下(偏转角为0°)土体变形图；图4为偏转工况下(偏转角为0.6°)土体变形图；图中标记为：1-隧道开挖断面；2-隧道收敛断面；3-整体直角坐标系(O-XZ)；4-局部极坐标系(O’-rθ)；5-单位开挖土体；O0-隧道中心；Dmax-最大土体变形；LDmax-局部最大土体变形。具体实施方式下面将结合附图，对本专利技术的优选实施例进行详细的描述。本专利技术所述方法的具体步骤如下：1)根据隧道设计的埋深，以及采用的双圆盾构机断面尺寸，绘制双圆盾构隧道开挖形状、大小和位置的二维断面图。如图1所示。2)根据隧道施工前期的地质调查和试验确定的土体损失率，计算隧道的收敛半径R2：R2＝R1(1-ε)0.5(1)式中：R1代表隧道的开挖半径，ε为土体损失率。3)在偏转工况下，隧道绕其中心O0滚动偏转，结合收敛半径和施工偏转角大小，绘出隧道收敛的形状、大小和位置。如图1所示。4)以双圆盾构隧道中心O0正上方的地表位置为坐标原点O,垂直于隧道轴线方向为X轴,竖直方向为Z轴,建立整体直角坐标系。如图1所示。5)以任一单位开挖土体所在圆心为局部坐标原点O‘，X方向为r轴，逆时针方向为θ轴，建立局部极坐标系。如图1所示。6)根据弹性介质理论，采用镜像法和坐标转换，建立单位土体损失引起土体变形的函数关系：竖向变形：横向变形：式中：(x,z)为土体变形的目标点坐标，(r,θ)和(X0,Z0)为单位土体损失所在的位置极坐标及其所在圆心坐标，ν为土体泊松比。7)为了方便积分，根据隧道开挖和收敛确定的土体损失分布，将其划分为4个子域的代数和：S＝SI+SII-SIII-SIV，每个子域代表一个不完全的圆面积。如图2所示。根据几何关系，各积分子域的数学描述为：SI：{(r,θ)|0≤r≤R1，arccos(t/R1)≤θ≤2π-arccos(t/R1)}∪{(r,θ)|0≤r≤t/cosθ，-arccos(t/R1)≤θ≤arccos(t/R1)}(4)SII：{(r,θ)|0≤r≤R1，-π+arccos(t/R1)≤θ≤π-arccos(t/R1)}∪{(r,θ)|0≤r≤t/cos(π-θ)，π-arccos(t/R1)≤θ≤π+arccos(t/R1)}(5)SIII：{(r,θ)|0≤r≤R2，arccos(t/R2)+α≤θ≤2π-arccos(t/R2)+α}∪{(r,θ)|0≤r≤t/cos(α-θ)，-arccos(t/R2)+α≤θ≤arccos(t/R2)+α}(6)SIV：{(r,θ)|0≤r≤R2，-π+arccos(t/R2)+α≤θ≤π-arccos(t/R2)+α}∪{(r,θ)|0≤r≤t/cos(π-θ+α)，π-arccos(t/R2)+α≤θ≤π+arccos(t/R2)+α}(7)式中：α代表盾构偏转角，t代表双圆隧道圆心距之半；8)通过对4个子域分别求积分计算总的土体变形分量：总的竖向变形：总的横向变形：式中：h代表双圆隧道中心的埋深。9)根据土体变形分量确定土体变形量及方向：变形量：变形方向：实施例：根据隧道设计和地勘资料：采用双圆盾构隧道施工，盾构机断面直径为9.36m，宽度为15.86m，隧道在计算土体变形所在断面埋深为20.68m，土体泊松比为0.38，土体损失率为0.83％。根据盾构机尺寸得到R1＝4.68m，t＝3.25m；代入公式(1)得到R2＝4.66m；代入公式(2)～(3)得到单位土体损失引起土体变形的函数：(1)正常工况下(偏转角α为0°)根据公式(4)～(7)划分4个积分子域，各子域的积分界限为：SI：{(r,θ)|0≤r≤4.68，0.803≤θ≤5.480}∪{(r,θ)|0≤r≤3.25/cosθ，-0.803≤θ≤0.803}SII：{(r,θ)|0≤r≤4本文档来自技高网...

【技术保护点】
一种双圆盾构隧道施工偏转引起土体变形的计算方法，其特征在于：包括以下步骤：1)绘制双圆盾构隧道开挖形状、大小和位置的二维断面图；2)通过土体损失率确定隧道的收敛半径R2：R2＝R1(1‑ε)0.5式中：R1代表隧道的开挖半径，ε为土体损失率；3)结合收敛半径和双圆盾构施工偏转角，绘出隧道收敛的形状、大小和位置；4)以双圆盾构隧道中心正上方的地表位置为坐标原点,垂直于隧道轴线方向为X轴,竖直方向为Z轴,建立整体直角坐标系；5)以任一单位开挖土体所在圆心为坐标原点，X方向为r轴，逆时针方向为θ轴，建立局部极坐标系；6)根据弹性介质理论，采用镜像法和坐标转换，建立单位土体损失引起土体变形的函数关系：竖向变形：uz(x,z;r,θ;X0,Z0)=-1π{z-Z0+rsinθ[(x-X0-rcosθ)2+(z-Z0+rsinθ)2]]]>-(3-4v)(z+Z0-rsinθ)[(x-X0-rcosθ)2+(z+Z0-rsinθ)2]+2z[(x-X0-rcosθ)2-(z+Z0-rsinθ)2][(x-X0-rcosθ)2+(z+Z0-rsinθ)2]2}r]]>横向变形：ux(x,z;r,θ;X0,Z0)=x-X0-rcosθπ{1[(x-X0-rcosθ)2+(z-Z0+rsinθ)2]]]>+3-4v[(x-X0-rcosθ)2+(z+Z0-rsinθ)2]-4z(z+Z0-rsinθ)[(x-X0-rcosθ)2+(z+Z0-rsinθ)2]2}r]]>式中：(x,z)为土体变形的目标点坐标，(r,θ)和(X0,Z0)为单位土体损失所在的位置极坐标及其所在圆心坐标，ν为土体泊松比；7)根据隧道开挖和收敛确定的土体损失分布，将其划分为4个子域的代数和：S＝SI+SII‑SIII‑SIV；8)通过对4个子域分别求积分计算总的土体变形分量：总的竖向变形：Uz(x,z)=∫∫SIuz(x,z;r,θ;-t,h)drdθ+∫∫SIIuz(x,z;r,θ;t,h)drdθ]]>-∫∫SIIIuz(x,z;r,θ;-tcosα,h-tsinα)drdθ-∫∫SIVuz(x,z;r,θ;tcosα,h+tsinα)drdθ]]>总的横向变形：Ux(x,z)=∫∫SIux(x,z;r,θ;-t,h)drdθ+∫∫SIIux(x,z;r,θ;t,h)drdθ]]>-∫∫SIIIux(x,z;r,θ;-tcosα,h-tsinα)drdθ-∫∫SIVux(x,z;r,θ;tcosα,h+tsinα)drdθ]]>式中：α代表盾构偏转角，t代表双圆隧道圆心距之半，h代表双圆隧道中心的埋深；9)根据土体变形分量确定土体变形量及方向：变形量：U(x,z)=Ux2(x,z)+Uz2(x,z)]]>变形方向：β(x,z)=arcsinUz(x,z)U(x,z).]]>...

【技术特征摘要】
1.一种双圆盾构隧道施工偏转引起土体变形的计算方法，其特征在于：包括以下步骤：1)绘制双圆盾构隧道开挖形状、大小和位置的二维断面图；2)通过土体损失率确定隧道的收敛半径R2：R2＝R1(1-ε)0.5式中：R1代表隧道的开挖半径，ε为土体损失率；3)结合收敛半径和双圆盾构施工偏转角，绘出隧道收敛的形状、大小和位置；4)以双圆盾构隧道中心正上方的地表位置为坐标原点,垂直于隧道轴线方向为X轴,竖直方向为Z轴,建立整体直角坐标系；5)以任一单位开挖土体所在圆心为坐标原点，X方向为r轴，逆时针方向为θ轴，建立局部极坐标系；6)根据弹性介质理论，采用镜像法和坐标转换，建立单位土体损失引起土体变形的函数关系：竖向变形：横向变形：式中：(x,z)为土体变形的目...

【专利技术属性】
技术研发人员：黄达，曾彬，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;85