预测盾构沿曲线路径掘进导致周围土体应力场变化的方法技术

本发明专利技术公开了一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法，包括以下步骤：（1）运用三维源汇法原理，得到半无限体内单位体积空隙引起的空间任一点应力增量解答；（2）根据盾构沿曲线路径掘进特点，确定实际掘进期间的地层损失量；（3）开发盾构沿曲线路径掘进时由实际地层损失量导致的周围土体应力场沿轴线方向及环开挖面方向变化的计算程序，以获得所述盾构沿曲线路径掘进时沿轴线方向变化的地层应力以及环开挖面变化的地层应力。本发明专利技术的优点是：适用于土体损失间隙为三维空间工况，以沿曲线路径掘进为例，施工时需要曲线超挖，运用本发明专利技术所述的预测程序，可准确计算三维超挖间隙导致周围土体应力场的变化。算三维超挖间隙导致周围土体应力场的变化。算三维超挖间隙导致周围土体应力场的变化。

【技术实现步骤摘要】
预测盾构沿曲线路径掘进导致周围土体应力场变化的方法


[0001]本专利技术属于盾构及地下工程
，具体涉及一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法。

技术介绍

[0002]修建地铁时盾构机的掘进路径不一定全为直线段，同时为了满足场地条件限制，盾构沿曲线路径掘进的工况在所难免。盾构隧道开挖过程中，盾尾管片脱空瞬间衬彻外壁与边界土层之间存在有非均匀间隙；特别对于曲线隧道，为使得盾构机顺利掘进，施工过程中又难免产生曲线内侧的超挖间隙。在施工过程中，曲线段超挖导致的间隙难免会对周围地层造成扰动，目前关于盾构掘进对地层扰动的预测程序计算对象为直线隧道，而对于曲线隧道，特别是土体损失导致地层应力的计算程序鲜有报道。更值得注意的是，关于盾构掘进引起地层扰动的预测多局限于其对土体位移的影响，能够预测其对地层应力影响的计算程序或方法较少，而应力状态的变化是引起地层位移的主要原因，因此亟待寻求方便、快捷、可靠的计算程序预测盾构沿曲线路径掘进导致周围土体的应力场变化。

技术实现思路

[0003]本专利技术的目的是根据上述现有技术的不足之处，提供一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法，该方法以曲线路径掘进为例，运用开发盾构沿曲线路径掘进时由实际地层损失量导致的周围土体应力场沿轴线方向及环开挖面方向变化的计算程序，以获得所述盾构沿曲线路径掘进时沿轴线方向变化的地层应力以及环开挖面变化的地层应力。
[0004]本专利技术目的实现由以下技术方案完成：
[0005]一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法，其特征在于所述方法包括以下步骤：
[0006](1)运用三维源汇法原理，得到半无限体内单位体积空隙引起的空间任一点应力增量解答；
[0007](2)根据盾构沿曲线路径掘进特点，确定实际掘进期间的地层损失量；
[0008](3)开发盾构沿曲线路径掘进时由实际地层损失量导致的周围土体应力场沿轴线方向及环开挖面方向变化的计算程序，以获得所述盾构沿曲线路径掘进时沿轴线方向变化的地层应力以及环开挖面变化的地层应力。
[0009]1.根据权利要求1所述的一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法，其特征在于步骤(1)包括以下步骤：
[0010](1.1)假设土体为无边界的无限体，给出在所述无限体中点F(x0,y0,z0)处的单位体积空隙引起点P(x,y,z)处的位移分量S
k1
：
[0011][0012]式中：
[0013]k为函数自变量，为x或y或z；
[0014](k
–
k0)＝(x
–
x0)、(y
–
y0)、(z
–
z0)中的一种；
[0015]r1＝[(x
‑
x0)2+(y
‑
y0)2+(z
‑
z0)2]1/2
；
[0016](1.2)获得点F(x0,y0,z0)的镜像位置点F
′
(x0,y0,
–
z0)处大小相等的体积膨胀在点P(x,y,z)处的位移分量：
[0017][0018]式中：
[0019]m为x或y；
[0020](m
–
m0)＝(x
–
x0)或(y
–
y0)；
[0021]r2＝[(x
‑
x0)2+(y
‑
y0)2+(z+z0)2]1/2
；
[0022](1.3)已知弹性力学基本方程，得到步骤(1.1)和步骤(1.2)产生的应变、应力解答：
[0023][0024][0025]式中：
[0026]k为函数自变量，为x或y或z；
[0027]S
k2
引用S
m2
的计算公式；
[0028]m＝x，或y；
[0029]G为土体剪切模量；
[0030]μ为土体泊松比；
[0031](1.4)获得以上两步产生沿m轴方向的应力计算程序：
[0032][0033]其中：
[0034]m轴是指x轴或y轴；
[0035](m,n)＝(x,y)，或(y,x)；
[0036]所述两步是指所述无限体内的空隙及其镜像间隙引起任意点应力分量的解答，为符合实际边界条件，即半无限体条件，应将前两步在地表产生的剪应力(Gγ
xz
；Gγ
yz
)反方向作用于地表，可求得其产生的应力分量：
[0037][0038][0039]式中：
[0040]b、c、u、t均为函数自变量；
[0041]r3＝[(x
‑
u)2+(y
‑
t)2+z2]1/2
；
[0042]以上3步解答之和即为半径为1的空隙引起的任一点应力增量解答，从而单位体积空隙产生的应力增量为：
[0043](m＝x或y)。
[0044]步骤(2)包括以下步骤：为满足曲线段转弯目的，所述盾构偏航推进时曲线内侧需要超挖，假设所述盾构的铰接装置的存在将所述盾构的机体分为前盾和后盾，曲线段超挖量满足盾尾覆盖2环管片即可，管片环宽为b1，同时考虑到开挖后盾的及时支护，土体位移受到三维约束，因此应对开挖面超挖量给予折减，根据经验取为1/3，则开挖面超挖间隙ω计算公式为：
[0045][0046]式中：
[0047]Q为曲线隧道轴线的曲率半径；
[0048]R为所述盾构的外径；
[0049]b1为管片环宽。
[0050]步骤(3)包括以下步骤：在曲线超挖影响下，周围土体所引起的应力增量计算程序为：
[0051][0052]式中：
[0053]θ为曲线隧道横断面的圆心及其在z轴上投影点的连线和三维直角坐标系中oxz平面之间的夹角；
[0054]q为函数自变量；
[0055]l为所述盾构的掘进长度；
[0056]h为隧道轴线埋深。
[0057]表征隧道周围应力场的三维空间分布的方法为：假设曲线隧道掘进轴线是水平面内半径为Q的一段圆弧，选取的计算路径如下：一方面，在所述曲线隧道掘进轴线所在的水平面内，取一段与所述曲线隧道掘进轴线为同心圆弧的曲线为计算路径，两所述曲线的径向距离r0为nR，其中，n＝2，3，4，
…
；所述计算路径位于曲线隧道外侧；另一方面，沿所述曲线隧道掘进轴线选取一特定断面，将所述特定断面与所述曲线隧道掘进轴线的交点作为圆心，在所述特定断面画圆，半径r0＝nR，其中，n＝2，3，4，
…
；计算不同r0情况下，圆周各点位置处应力场的环向分布。
[0058]本专利技术的优点是：与现有技术相比，本专利技术适用于土体损失间隙为三维空间工况，以沿曲线路径掘进为例，施工时需要曲线超挖，运用本专利技术所述的预测程序，可准本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法，其特征在于所述方法包括以下步骤：(1)运用三维源汇法原理，得到半无限体内单位体积空隙引起的空间任一点应力增量解答；(2)根据盾构沿曲线路径掘进特点，确定实际掘进期间的地层损失量；(3)开发盾构沿曲线路径掘进时由实际地层损失量导致的周围土体应力场沿轴线方向及环开挖面方向变化的计算程序，以获得所述盾构沿曲线路径掘进时沿轴线方向变化的地层应力以及环开挖面变化的地层应力。2.根据权利要求1所述的一种预测盾构沿曲线路径掘进导致周围土体应力场变化的方法，其特征在于步骤(1)包括以下步骤：(1.1)假设土体为无边界的无限体，给出在所述无限体中点F(x0,y0,z0)处的单位体积空隙引起点P(x,y,z)处的位移分量S
k1
：式中：k为函数自变量，为x或y或z；(k
–
k0)＝(x
–
x0)、(y
–
y0)、(z
–
z0)中的一种；r1＝[(x
‑
x0)2+(y
‑
y0)2+(z
‑
z0)2]
1/2
；(1.2)获得点F(x0,y0,z0)的镜像位置点F
′
(x0,y0,
–
z0)处大小相等的体积膨胀在点P(x,y,z)处的位移分量：式中：m为x或y；(m
–
m0)＝(x
–
x0)或(y
–
y0)；r2＝[(x
‑
x0)2+(y
‑
y0)2+(z+z0)2]
1/2
；(1.3)已知弹性力学基本方程，得到步骤(1.1)和步骤(1.2)产生的应变、应力解答：(1.3)已知弹性力学基本方程，得到步骤(1.1)和步骤(1.2)产生的应变、应力解答：式中：k为函数自变量，为x或y或z；
S
k2
引用S
m2
的计算公式；m＝x，或y；G为土体剪切模量；μ为土体泊松比；(1.4)获得以上两步产生沿m轴方向的应力计算程序：其中：m轴是指x轴或y轴；(m,n)＝(x,y)，或(y,x)；所...

【专利技术属性】
技术研发人员：李少华，黄昌富，姚铁军，赵中华，栾焕强，张帆舸，岳粹洲，田国文，
申请(专利权)人：中铁十五局集团有限公司，
类型：发明
国别省市：