一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法与流程

1.本发明属于岩土工程位移预测领域，具体涉及一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法。


背景技术：

2.随着隧道工程的快速发展，盾构隧道施工期管片的变形机理及控制技术已逐渐成为工程实践中关注的问题。
3.由各地的盾构施工监测数据可知，无论盾构处于地层条件较差的软土地层，还是处于地层条件较好的砂质、岩质地层，盾构隧道在掘进施工期间都在不同程度上出现了管片上浮问题。
4.盾构隧道管片上浮计算模型的建立、计算参数的选取对于探究管片上浮规律、预测上浮量以及隧道上浮控制等方面至关重要。
5.上浮计算模型方面，浆液的凝固是需要一定时间的，如何考虑浆液的时变性，建立浆液未凝固区的动态变化模型，以及浆液分布与压力的确认还需要进一步的研究。隧道管片变形、内力的计算更多地依靠数据模拟，合理的解析解还有待深入研究。


技术实现要素：

6.本发明的目的是在于提供一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法，本发明易行可靠，优势明显，适用范围广，可用于计算不同土体条件下盾构隧道施工期纵向上浮变形特性。
7.本发明采用以下技术方案实现：
8.一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法，具体如下：
9.盾构开挖掘进期间，注浆浆液或孔隙水产生的上浮力会引起隧道的纵向上浮变形。这是一个土与结构相互作用的问题，为建立考虑浆液凝固特性的盾构隧道管片施工期上浮力学分析模型，需要做以下基本假定：
10.(1)地基模型假定为winkler弹性地基模型，即将周围土体与既有隧道之间的相互作用关系假定为一系列相互独立的弹簧。注浆浆液未凝固区的地基反力系数假定为线性分布，地基反力系数由0按线性递增过渡至浆液凝固区。浆液凝固区的地基系数假定为定值。
11.(2)盾构隧道模型假定为欧拉梁模型，并考虑管片拼装造成的纵向刚度折减。
12.(3)在浆液未凝固区，注浆浆液产生的上浮作用力为线性分布。
13.(4)隧道模型在盾构尾端受到盾尾弹簧约束，以模拟盾构机对管片的约束。由于在远离上浮影响区处，约束形式对结构的变形影响一般较小，以自由端处理。
14.(5)盾构机始终沿隧道设计轴线掘进，不考虑盾构轴线偏移对隧道纵向力学行为的影响。
15.(6)不考虑盾尾同步注浆压力对隧道纵向力学行为的影响。
16.建立考虑浆液凝固特性的盾构隧道管片施工期上浮力学分析模型，具体为：
17.搁置在winkler地基中的欧拉梁，在距盾尾水平距离为x处浆液上浮力f(x)的作用下，可得到关于梁竖向挠曲变形：
[0018][0019]
式中，(ei)
eq
为隧道的等效抗弯刚度，(ei)
eq
＝ξe
l
i,ξ为抗弯刚度折减系数，e
l
为盾构管片的弹性模量，i为管片截面的惯性矩；w(x)为隧道竖向变形；k(x)为地基系数。
[0020]
求解式(1)即可得到在浆液上浮力作用下隧道的纵向变形。由于式(1)为四阶常微分方程，在数学求解上存在一定难度。为简化计算，采用有限差分法求解。隧道离散为n 4个节点单元(其中在隧道两端为2个虚节点单元)，每个节点单元长度为l。根据标准有限差分原理，式(1)中各微分项的有限差分形式如下：
[0021][0022]
式中，wi为节点i竖向位移；ki为节点i的地基系数；fi为节点i的上浮力；h＝l/n，l为隧道长度，i＝0,1,2,3,
…
,n。
[0023]
假定隧道两端为自由端，边界条件则为：
[0024]
m0＝mn＝0
ꢀꢀꢀ
(3)
[0025]
q0＝qn＝0
ꢀꢀꢀ
(4)
[0026]
式(3)和(4)可改写为：
[0027][0028][0029]
求得节点方程如下：
[0030][0031]
式中，{w}
(n 1)
×1为隧道竖向变形列向量，{w}
(n 1)
×1＝{w0,w1,w2,
…
,w
n-1
,wn}
t
；{f}
(n 1)
×1为上浮力列向量，{f}
(n 1)
×1＝{f0,f1,f2,
…
,f
n-1
,fn}
t
；{k
t
}和{ks}为隧道刚度矩阵和地基刚度矩阵。
[0032]
[0033][0034]
上浮力学分析模型参数取值为：
[0035]
(1)盾尾弹簧系数kd[0036]
为模拟盾构盾尾对隧道的约束作用，可在盾尾内的两环管片上假定较大的弹簧系数kd。盾尾弹簧系数kd不小于3000mpa。
[0037]
(2)浆液未凝固区长度s
[0038]
当浆液抗剪强度在时间t后增长至某一临界值，临界抗剪强度的表达式为：
[0039][0040]
式中，ρg为浆液密度，g为重力加速度，s为注浆层厚度，d为隧道外径，g
l
为管片自身重量，可由下式表示：
[0041][0042]
式中，γc为混凝土管片重度，一般取值为26kn
·
m-3
；δc为管片厚度与管片半径的比值，一般取值为0.08～0.10。
[0043]
假定浆液到达临界抗剪强度时即处于凝固状态，根据盾构平均掘进速度vs，可确定浆液未凝固区的长度s为：
[0044]
s＝vst
ꢀꢀꢀ
(16)
[0045]
(3)浆液上浮力f(x)
[0046]
浆液的上浮力主要跟同步注浆浆液配比及浆液的时变特性有关。对于常规抗剪型浆液，浆液上浮力在初期接近阿基米德定律计算得到的浮力，而之后浆液浮力逐渐减小并保持为水的浮力。对于可硬性抗剪型浆液，浆液初期表现与常规抗剪型浆液一致，但在后期浆液硬化之后，浆液本身的重量将作为荷载抵抗来自饱和土中孔隙水压产生的浮力作用。
[0047]
直接利用浆液上浮力时变规律进行模型计算，f(x)为距盾构盾尾水平距离为x处浆液上浮力；
[0048]
则常规抗剪型浆液的上浮力f(x)可由式(17)表示：
[0049][0050]
式中，ρg为浆液密度，ρw为水密度，g为重力加速度，v为管片体积，vs为盾构平均掘进速度，t1为常规抗剪型浆液上浮力时变规律第一阶段线性快速减小段终点对应时间。
[0051]
可硬性抗剪型浆液上浮力f(x)可由式(18)表示：
[0052][0053]
式中，ρg为浆液密度，ρw为水密度，g为重力加速度，v为管片体积，vs为盾构平均掘进速度，t2为可硬性抗剪型浆液上浮力时变规律第一阶段线性快速减小段终点对应时间，t3为可硬性抗剪型浆液上浮力时变规律第二阶段稳定不变段终点对应时间，t4为可硬性抗剪型浆液上浮力时变规律第三阶段二次线性减少段终点对应时间。
[0054]
(4)浆液凝固区地基系数kn[0055]
浆液凝固区的地基系数kn可通过下式得到：
[0056][0057][0058]
其中，η为深度参数；e为土体弹性模量，可取为3倍的土层压缩模量；d为隧道外径；(ei)
eq
为隧道的等效抗弯刚度，z为隧道埋深；ν为土体泊松比。
[0059]
(5)抗弯刚度折减系数ξ
[0060]
盾构隧道纵向抗弯刚度折减系数的取值为1/5-1/7。
[0061]
本发明的有益效果是：
[0062]
可针对分别使用常规抗剪型浆液和可硬性抗剪型浆液条件下，同时考虑盾构隧道的埋深、隧道的直径、盾构的掘进速度以及土体弹性模量等综合因素对盾构隧道施工期上浮性状的影响。通过对该上浮分析计算方法的分析，可得出使用可硬性抗剪型浆液、减小盾构掘进速度、改善上覆土体性质等能够在一定程度上减小盾构隧道施工期上浮量，为工程实践提供指导。
附图说明
[0063]
图1为本发明的盾构隧道管片施工期上浮力学分析模型示意图。
具体实施方式
[0064]
下面将结合实例进一步阐明本发明的内容，但这些实例并不限制本发明的保护范围，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围内。
[0065]
一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法，具体如下：
[0066]
盾构开挖掘进期间，注浆浆液或孔隙水产生的上浮力会引起隧道的纵向上浮变形。这是一个土与结构相互作用的问题，为建立考虑浆液凝固特性的盾构隧道管片施工期
上浮力学分析模型，需要做以下基本假定：
[0067]
(1)地基模型假定为winkler弹性地基模型，即将周围土体与既有隧道之间的相互作用关系假定为一系列相互独立的弹簧。注浆浆液未凝固区的地基反力系数假定为线性分布，地基反力系数由0按线性递增过渡至浆液凝固区。浆液凝固区的地基系数假定为定值。
[0068]
(2)盾构隧道模型假定为欧拉梁模型，并考虑管片拼装造成的纵向刚度折减。
[0069]
(3)在浆液未凝固区，注浆浆液产生的上浮作用力为线性分布。
[0070]
(4)隧道模型在盾构尾端受到盾尾弹簧约束，以模拟盾构机对管片的约束。由于在远离上浮影响区处，约束形式对结构的变形影响一般较小，以自由端处理。
[0071]
(5)盾构机始终沿隧道设计轴线掘进，不考虑盾构轴线偏移对隧道纵向力学行为的影响。
[0072]
(6)不考虑盾尾同步注浆压力对隧道纵向力学行为的影响。
[0073]
建立考虑浆液凝固特性的盾构隧道管片施工期上浮力学分析模型(如图1)，具体为：
[0074]
搁置在winkler地基中的欧拉梁，在距盾尾水平距离为x处浆液上浮力f(x)的作用下，可得到关于梁竖向挠曲变形：
[0075][0076]
式中，(ei)
eq
为隧道的等效抗弯刚度，(ei)
eq
＝ξe
l
i,ξ为抗弯刚度折减系数，e
l
为盾构管片的弹性模量，i为管片截面的惯性矩；w(x)为隧道竖向变形；k(x)为地基系数。
[0077]
求解式(1)即可得到在浆液上浮力作用下隧道的纵向变形。由于式(1)为四阶常微分方程，在数学求解上存在一定难度。为简化计算，采用有限差分法求解。隧道离散为n 4个节点单元(其中在隧道两端为2个虚节点单元)，每个节点单元长度为l。根据标准有限差分原理，式(1)中各微分项的有限差分形式如下：
[0078][0079]
式中，wi为节点i竖向位移；ki为节点i的地基系数；fi为节点i的上浮力；h＝l/n，l为隧道长度，i＝0,1,2,3,
…
,n。
[0080]
假定隧道两端为自由端，边界条件则为：
[0081]
m0＝mn＝0
ꢀꢀꢀ
(3)
[0082]
q0＝qn＝0
ꢀꢀꢀ
(4)
[0083]
式(3)和(4)可改写为：
[0084][0085][0086]
求得节点方程如下：
[0087][0088]
式中，{w}
(n 1)
×1为隧道竖向变形列向量，{w}
(n 1)
×1＝{w0,w1,w2,
…
,w
n-1
,wn}
t
；
{f}
(n 1)
×1为上浮力列向量，{f}
(n 1)
×1＝{f0,f1,f2,
…
,f
n-1
,fn}
t
；{k
t
}和{ks}为隧道刚度矩阵和地基刚度矩阵。
[0089][0090][0091]
上浮力学分析模型参数取值为：
[0092]
(1)盾尾弹簧系数kd[0093]
为模拟盾构盾尾对隧道的约束作用，可在盾尾内的两环管片上假定较大的弹簧系数kd。盾尾弹簧系数kd不小于3000mpa。
[0094]
(2)浆液未凝固区长度s
[0095]
当浆液抗剪强度在时间t后增长至某一临界值，临界抗剪强度的表达式为：
[0096][0097]
式中，ρg为浆液密度，g为重力加速度，s为注浆层厚度，d为隧道外径，g
l
为管片自身重量，可由下式表示：
[0098][0099]
式中，γc为混凝土管片重度，一般取值为26kn
·
m-3
；δc为管片厚度与管片半径的比值，一般取值为0.08～0.10。
[0100]
假定浆液到达临界抗剪强度时即处于凝固状态，根据盾构平均掘进速度vs，可确定浆液未凝固区的长度s为：
[0101]
s＝vst
ꢀꢀꢀ
(16)
[0102]
(3)浆液上浮力f(x)
[0103]
浆液的上浮力主要跟同步注浆浆液配比及浆液的时变特性有关。对于常规抗剪型
浆液，浆液上浮力在初期接近阿基米德定律计算得到的浮力，而之后浆液浮力逐渐减小并保持为水的浮力。对于可硬性抗剪型浆液，浆液初期表现与常规抗剪型浆液一致，但在后期浆液硬化之后，浆液本身的重量将作为荷载抵抗来自饱和土中孔隙水压产生的浮力作用。
[0104]
直接利用浆液上浮力时变规律进行模型计算，f(x)为距盾构盾尾水平距离为x处浆液上浮力；
[0105]
则常规抗剪型浆液的上浮力f(x)可由式(17)表示：
[0106][0107]
式中，ρg为浆液密度，ρw为水密度，g为重力加速度，v为管片体积，vs为盾构平均掘进速度，t1为常规抗剪型浆液上浮力时变规律第一阶段线性快速减小段终点对应时间。
[0108]
可硬性抗剪型浆液上浮力f(x)可由式(18)表示：
[0109][0110]
式中，ρg为浆液密度，ρw为水密度，g为重力加速度，v为管片体积，vs为盾构平均掘进速度，t2为可硬性抗剪型浆液上浮力时变规律第一阶段线性快速减小段终点对应时间，t3为可硬性抗剪型浆液上浮力时变规律第二阶段稳定不变段终点对应时间，t4为可硬性抗剪型浆液上浮力时变规律第三阶段二次线性减少段终点对应时间。
[0111]
(4)浆液凝固区地基系数kn[0112]
浆液凝固区的地基系数kn可通过下式得到：
[0113][0114][0115]
其中，η为深度参数；e为土体弹性模量，可取为3倍的土层压缩模量；d为隧道外径；(ei)
eq
为隧道的等效抗弯刚度，z为隧道埋深；ν为土体泊松比。
[0116]
(5)抗弯刚度折减系数ξ
[0117]
盾构隧道纵向抗弯刚度折减系数的取值为1/5-1/7，本实施例采用的取值为1/7。