8隧道地层结构计算解读.docx
- 文档编号:15312771
- 上传时间:2023-07-03
- 格式:DOCX
- 页数:86
- 大小:122.14KB
8隧道地层结构计算解读.docx
《8隧道地层结构计算解读.docx》由会员分享,可在线阅读,更多相关《8隧道地层结构计算解读.docx(86页珍藏版)》请在冰点文库上搜索。
8隧道地层结构计算解读
第八章隧道支护地层—结构分析方法
(杨林德、朱合华、丁文其)
第一节一般规定
8.1.1设计模型
1、与地层—结构分析方法相应的隧道设计模型是地层—结构模型。
其设计理念,是认
为围岩具有自支承能力,支护(含衬砌结构,下同)的作用是加固围岩,并与围岩联合组成
共同受力的整体,共同承受荷载的作用。
2、采用地层—结构分析方法进行设计计算时,计算区范围应同时包含支护和围岩地层。
3、采用地层—结构分析方法进行设计计算时,应同时考虑开挖施工步骤的影响。
初步
设计阶段可按常规施工方法选定开挖施工步骤,施工图设计阶段应改按施工组织设计制定的
技术方案确定。
4、采用地层—结构分析方法设计隧道时,应同时检验围岩的稳定性和支护结构的受力
变形状态是否满足强度条件及按使用要求确定的变形量限制条件。
5、采用地层—结构分析方法设计隧道时,内衬结构的工作状态应为弹性受力状态,或
经论证认为仍可保持稳定的弹塑性受力状态;初期支护(含开挖阶段增设的喷射混凝土层)
和围岩的工作状态可为弹塑性受力状态。
8.1.2地层—结构分析方法的适用地质条件
1、地层—结构分析方法的基础理念,是认为围岩具有自支承能力,并可由其与支护结
构共同组成承载体系。
因而这类方法适用于在具有一定自支承能力的围岩中建造的隧道支护
1
结构的计算。
2、V级及V级以上的围岩都具有一定的自支承能力,因而都可采用地层—结构分析方
法进行设计计算。
3、Ⅲ级及Ⅲ级以上的围岩自支承能力强,对在这些级别的围岩中建造的隧道,经验表
明对支护结构根据经验选定设计参数时已可使围岩保持稳定,因而规范规定一般不要求进行
计算。
4、IV级、V级围岩中建造的隧道一般采用复合式支护,对其宜采用地层—结构分析方
法进行设计计算。
但对在V级围岩中建造的浅埋隧道,围岩承载能力较低时仍宜采用荷载-
结构分析方法计算。
5、IV级围岩的自支承能力优于V级围岩,采用地层—结构分析方法进行设计计算时,
宜通过控制荷载释放过程,使IV级围岩中隧道内衬结构经受的荷载相对较小,围岩的自支
承能力可适度充分发挥。
6、VI级围岩的自支承能力差,宜采用荷载—结构分析方法对支护结构进行设计计算。
8.1.3计算方法及其适用场合
1、地层—结构分析方法的计算方法可分为解析解和数值法二类。
其中解析解只适用于
均匀介质中的圆形隧道,且只能计算若干典型工况。
对公路隧道的设计,可供采用的计算方
法通常是数值法。
2、地层—结构分析方法的数值法可分为有限单元法(FEM)、特征单元法(DDA)、边
界单元法(BEM)和有限差分法(FDM)等。
3、有限单元法因有既可模拟各级围岩的性态特征,又能反映断层、节理等地质构造的
影响,并能对开挖施工过程实行动态追踪等显著优点,因而适用于各级围岩(硬岩或软岩)
2
中的公路隧道设计的计算。
同时由于目前已有多种包括前、后处理在内的功能强大的程序软
件可供采用,这类方法是目前最常采用的一类算法。
4、特征单元法可较好模拟块体结构的性态,因而适用于围岩地层为块体状结构的硬岩
地层中的公路隧道的设计。
然因查明块体结构分布的几何特征及合理确定结构面性状的参数
均需开展较多的地质调查工作,公路隧道设计很难满足这一要求,因而这类方法一般仅在规
模较大的大跨度地下结构的设计研究中采用。
5、边界单元法用于均匀介质中的弹性、粘弹性问题的计算时才比有限单元法简捷,一
般仅适用于围岩介质的性态可用弹性、粘弹性模型近似模拟时的公路隧道的设计计算。
6、有限差分法因将控制方程改造为差分方程而具有少占内存等显著优点,目前已为
FLAC等程序吸收,可供各级围岩(硬岩或软岩)中的公路隧道设计计算采用。
8.1.4地层—结构分析方法的荷载
1、采用地层—结构分析方法设计公路隧道时,作用在隧道结构上的荷载可按《公路隧
道设计规范》JTGD70-2004表6.1.1分类,并按规范提出的方法计算。
但其中的围岩压力应
为释放荷载。
2、采用地层-结构分析方法进行隧道设计计算时,对在隧道结构上可能同时出现的荷载,
应按规范规定的原则进行组合,并按最不利组合进行计算和设计。
3、地层-结构分析方法对初期支护和二次衬砌的计算都适用,但在进行具体计算时,对
不同阶段的计算应根据实际情况取用不同的荷载组合。
4、释放荷载与初始地应力、围岩材料的性态、开挖施工步骤及结构施做时机等有关。
工程设计中,释放荷载的计算需按当前地应力(不一定是初始地应力)计算。
各类因素的影
响,则可由根据开挖施工步骤和支护施做时机等设定相应的荷载释放过程体现。
3
5、对初期支护的设计计算,级别相对较高的围岩可取用较大的释放荷载分担比,使初
期支护和围岩承受较大的荷载,结构产生较小的变形;级别相对较低的围岩,则相反。
6、鉴于围岩材料的变形常随时间而增长,由数层喷射混凝土层和内衬结构联合组成复
合式支护时,各层支护结构经受围岩压力作用的程度将有差异。
这类力学现象也可通过控制
荷载释放过程模拟,即将与初期支护(含开挖阶段增设的喷射混凝土层)和围岩的受力状态
相应的释放荷载分担比(即《公路隧道设计规范》JTGD70-2004条文说明中表9-1的第一
列数据),按支护层数(不含内衬结构)合理分配。
4
第二节初始地应力与开挖效应
8.2.1初始地应力
8.2.1.1初始地应力的含义与组成
1、初始地应力是指天然状态下存在于岩体或土体介质内部的应力。
2、在未经扰动的岩层中开挖隧道时,岩体内部存在的应力即为初始地应力,相应的应
力状态可称为初始应力状态。
3、在先期经受过开挖扰动影响的岩层中开挖隧道时,岩体内部存在的应力并不是初始
地应力,而是包含先期开挖扰动影响的合应力,但相应的应力状态仍可称为初始应力状态。
4、初始地应力由初始自重应力及构造应力组成,表达式可写为:
gt
(8-2-1)
式中:
—初始地应力;
g—自重应力分量;
t—构造应力分量。
8.2.1.2初始地应力的确定方法
1、初始地应力的确定方法有水压致裂法、钻孔应力法、位移反分析法和回归分析法等。
前两种方法属于直接测量法,后两种方法属于反分析法。
其中直接测量法通常有需要经费较
多、花费时间较长的显著弱点,回归分析法需在工程所在地区的数个点上测得地应力值后才
能采用,因而通常都仅适用于水电站工程等的设计研究。
对公路隧道,这些方法一般仅在长
度特长,地质条件特复杂的场合才考虑采用。
5
2、水压致裂法通过测量垂直钻孔的孔壁开始出现张裂缝时的破裂水压力和在水泵停开
后使水压裂缝保持张开状态所必须的封井压力,进而得出地应力值。
水压致裂法用于测量深
层岩体的地应力时,测得的地应力即为测点的初始地应力。
适用于初始地应力的一个主应力
为垂直应力的情况。
3、钻孔应力法通过量测套芯应力解除前后钻孔孔径的变化确定地应力。
钻孔深度超过
扰动影响区时测得的地应力即为测点的初始地应力。
适用于测点范围内岩性均匀,且岩芯无
大的裂隙通过的情况。
4、位移反分析法
1)位移反分析法利用在工程现场测得的,由开挖扰动引起的位移量确定初始地应力,
因分析过程(位移→荷载)与常见过程(荷载→位移)相反而得名。
按算法特点可分为正反
分析法和逆反分析法两类。
后者由正分析计算的逆解过程确定初始地应力,因对非线性问题
的分析难于得到解析式,以及需要针对各类具体情况分别编制专用程序而很少采用;前者则
通过正分析计算的优化逆解逼近过程确定初始地应力,因可主要采用常用正算程序计算而显
得简便。
目前采用的方法一般都是正反分析法。
2)采用正反分析法确定初始地应力时,可供采用的用于促使优化过程收敛的算法有
单纯形法、阻尼最小二乘法、遗传算法、遗传模拟退火算法以及混合遗传算法等,用于约束
优化过程的目标函数可统一表示为:
N
[1
uj
]
2
JMin
uj
j
1
(8-2-2)
式中:
J—目标函数;
uj—位移量实测值;
uj—位移量真值;
N—位移量测值个数。
6
3)采用正反分析法确定初始地应力时,需先对初始地应力场的分布规律作假设。
通
常认为在工程活动涉及的岩层内,自重应力自上而下呈线性规律分布(地表为零),构造应
力可假设为均布应力,或沿深度分段均布的应力,或沿深度呈线性规律分布的应力(地表不
一定为零)。
4)采用正反分析法确定初始地应力时,目标未知数宜选为沿计算区域的边界分布的
应力。
求得边界应力后,即可由数值分析的正演方法算得计算区域内各点的初始地应力。
5)按计算区域的几何特征,正反分析法可分为二维平面应变问题和三维空间问题的
反分析方法两类。
前者采用的目标未知数为沿边界线均布或线性分布的线荷载(应力),后
者则为沿边界面均布或线性分布的面荷载(应力)。
6)对计算区域的边界需设定边界条件时,采用的方法宜与数值分析的正演分析法相
同。
7)根据上述原则建立的正反分析法可有许多种,8.2.1.3列出的正算逆解逼近法为其
中之一,适用于可简化为平面应变问题计算的情况,可供参考。
8)鉴于围岩位移量的量测值与初始地应力及工程岩体的弹性参数值(尤其是弹性模
量值或变形模量值)都关系密切,采用位移反分析法确定初始地应力时,应同时确定弹性参
数值。
5、回归分析法
1)回归分析法类属应力反分析法,特点为利用散布在工程所在区域内的数个地点的
初始地应力实测值,借助根据数理统计原理建立的优化过程反演确定区域范围内的初始地应
力场的分布规律,从而得出工程建设地点的初始地应力的估计值。
2)采用回归分析法确定工程所在区域的初始地应力场的分布时,可供采用的用于促
使优化过程收敛的算法可与位移反分析法相同,用于约束优化过程的目标函数则需表示为:
7
N
6
JMin
k
k
2
[wkij
ij]
k
1i,i1
(8-2-3)
k
式中:
ij—应力分量的量测值;
k
ij—应力分量的真值;
N—测点数;
wk—加权系数,一般可令wk=1。
3)采用回归分析法确定初始地应力时,也需先对初始地应力场的分布规律作假设。
其原则,可与位移反分析法相同。
4)采用回归分析法确定初始地应力时,计算方法均属三维空间问题的反分析法,目
标未知数常选为边界面力。
求得边界面力后,再由数值分析的正演分析方法算得工程所在部
位的初始地应力。
5)对计算区域的边界需设定边界条件时,采用的方法宜与数值分析的正演分析法相
同。
6、在丘陵地带建造公路隧道时,围岩的初始地应力场通常即是自重应力场,其分布规
律可借助正演分析的数值方法通过计算确定,也可将垂直应力取为上覆地层重量之和,并按
给定水平侧压力系数法确定侧压力。
后者的计算公式可参见《公路隧道设计规范》JTG
D70-2004中的附录J,或有关文献。
8.2.1.3二维平面应变问题反分析计算的正算逆解逼近法
1、线弹性问题的反分析方法
1)基本假设
(1)隧道围岩的工作状态为线弹性受力状态。
(2)横断面上隧道围岩的受力变形状态符合平面应变假设。
8
(3)某一地点横断面上隧道围岩的初始地应力由线性分布的自重应力和均布构造
应力组成。
2)基本方程
(1)对二维平面应变问题,在满足最小二乘原理条件下,线弹性问题位移反分析
计算的基本方程为:
N
2
[Dk
(dg
ijtdij0)k][dij0]k
0
k1
i,j1
(8-2-4)
(ij)
式中:
N-位移量测值的总数;
Dk-第k个位移量的实测值;
dg-由自重应力引起的,测点在量测方向上的位移量(指由与自重应力相应的释放
荷载引起的位移);
t
ij-均布构造应力分量的量值;
000
dij-由单位均布构造应力ij=1引起的,测点在量测方向上的位移量(指由与ij=
1相应的释放荷载引起的位移)。
(2)式(8-2-4)中Dk为实测值,dg和dij0为可由数值分析的正演方法得出的已知
值,但因dg
、dij0
的量值包含围岩材料弹性性态参数的影响,反分析计算的目标未知数需同
时包括隧道计算断面上的均布构造应力分量及围岩材料的
E、值。
但对具体工程的分析,
可在根据地质条件选定
E、值的基础上,仅将均布构造应力分量作为反分析计算的目标
未知数。
(3)采用这一方法计算时,初始地应力分量的计算式为:
Pij
g
t
0
(8-2-5
ij
ij
ij
)
(i,j
1,2;i
j)
9
式中:
Pij-横断面上隧道围岩中的初始地应力分量;
g
ij-横断面上隧道围岩中与自重相应的地应力分量。
t
ij的含义与式(8-2-4)相同。
3)计算方法
(1)采用这一方法进行计算时,可先算得dg、dij0,将其连同Dk代入式(8-2-4),
tt
得到关于ij的三元一次线性代数方程组,由其解得ij,然后由式(8-2-5)确定初始地应力
Pij。
(2)dg、dij0可采用数值分析的正演方法计算,常用方法为有限元方法。
注意两
0
者的荷载应分别为与自重应力及ij=1相应的释放荷载。
(3)采用这一方法进行计算时,围岩地层的弹性参数需同时确定。
通常可先根据
地质条件选定E、值,通过反分析计算确定Pij,然后根据Pij算得与量测位移相应的位移
量。
如果两者相差较大,则调整E值后重新进行反分析计算,直到两者相差较小。
计算过程
中,用于约束优化过程的目标函数即为式(8-2-2)。
2、弹塑性问题的反分析方法
1)基本假设
除假设隧道围岩的工作状态可为弹塑性受力状态外,其余假设均与线弹性问题的反分析
方法相同。
2)基本方程
(1)对二维平面应变问题,在满足最小二乘原理条件下,弹塑性问题位移反分析
计算的基本方程可选为:
N
2
[Dk
(dg
ijtdijep)k][dijep]k
0
k1
i,j1
(8-2-6)
(ij)
10
ep0
式中dij为弹塑性受力状态下,由单位均布构造应力ij=1引起的,测点在量测方向上的位
0
移量(指由与ij=1相应的释放荷载引起的位移)。
其余符号的含义均与式(8-2-4)相同。
(2)式(8-2-6)中Dk仍为实测值,dg、dijep也仍为可由数值分析的正演方法得出的已知值,但因dg、dijep的量值包含围岩材料弹塑性性态参数的影响,反分析计算的目标
未知数需同时包括隧道计算断面上的均布构造应力分量及围岩材料的E、、c、值。
这
类问题属于高度非线性问题。
对具体工程的分析,仍可在根据地质条件选定E、、c、值
的基础上,仅将均布构造应力作为反分析计算的基本未知数。
(3)采用这一方法计算时,理论上初始地应力的计算式仍为式(8-2-5),然因弹
塑性问题分析的计算方法多为增量迭代法,初始地应力的计算也将包含迭代修正过程。
3)计算方法
(1)与弹性问题的反分析相比较,采用这一方法计算时,仍可先算得dg、dijep,
后由式(8-2-6)解得ijt,及由式(8-2-5)确定初始地应力Pij,区别是这时dg、dijep的计
算包含围岩材料弹塑性性态的影响,使dijep、ijt及Pij的取值均需经过迭代修正,才能最终确
定。
(2)与弹性问题的反分析比较,dg、dijep仍可采用数值分析的正演方法计算,
两者的荷载也仍分别为与自重应力及ij0=1相应的释放荷载,区别是计算过程需考虑围岩地
层材料进入弹塑性状态后对变形的影响,尤应注意计算dijep时围岩地层的应力水平应高于自
重应力。
具体计算方法可有多种。
用于工程问题的分析时,可根据采用的计算程序选用。
(3)采用以上方法计算时,如根据由反分析计算确定的
Pij算得的位移量的计算
值与实测值相比误差持续较大,可重新选定
E、、c、值后再进行反分析计算。
11
8.2.2开挖效应
8.2.2.1开挖效应的含义
隧道开挖后,地层的初始应力平衡状态被破坏,洞周围岩将在沿隧道周边分布的,与初
始应力大小相等、方向相反的不平衡力作用下发生变形,由此产生附加应力场与位移场。
这
类由隧道开挖引起的沿隧道洞周作用的不平衡力习称释放荷载,在释放荷载作用下围岩产生
附加应力场与位移场的现象称开挖效应。
8.2.2.2开挖效应的计算原理与方法
1、开挖效应可通过在洞周边界上设置释放荷载进行计算,常用的计算方法为有限单元
法。
2、释放荷载是与隧道洞周上的初始围岩应力大小相等、方向相反的分布应力。
在未经
扰动的岩体中开挖隧道时,初始围岩应力即为初始地应力;在已扰动过的岩体中开挖隧道时,
应为围岩当前的初始应力。
3、采用有限元方法计算时,释放荷载需转化为等效结点力。
4、释放荷载及其等效结点力的计算方法见,采用有限元方法计算开挖效应的具体
方法见。
8.2.3释放荷载的计算方法
8.2.3.1计算方法的种类及其适用场合
1、释放荷载可采用单元应力法、绕结点平均法或Mana法计算。
2、上述三类方法都可用于计算释放荷载。
但因Mana法在建立具体算法时对边界结点
间围岩应力场变化规律的假设与有限单元法相同,并由此易于编制程序,因而宜优先采用。
12
8.2.3.2单元应力法
1、采用单元应力法计算时,需先根据初始地应力或与前一步开挖相应的应力场,算得
预计开挖边界上各结点的应力,并假定各结点间的应力呈线性分布,然后反转开挖边界上各
结点应力的方向(即改变其符号),得到释放荷载。
2、将开挖边界上呈线性分布的释放荷载转化为等效结点力的计算式为:
Pi
1
[2
i
(bb
)
i1b
2
i1b2
i
(aa
2
)
i1a
2
i1a]
x
6
x
1
2
x
x
1
xz
1
xz
xy1
(8-2-7)
Pi
1[2
i
(a
a
)
i1a
2
i1a
1
2
i
(bb)
i1b
2
i1b]
z
6
z
1
2
z
z
xz
12
xz
xy1
(8-2-8)
a1
xi
1
xi
a2
xi
xi
1
b1
zi
zi
1
b2
zi
1
zi
i
i
i
i
式中:
x、
z、xz—分别为结点
上与初始地应力或与前一步开挖相应的正应力和剪应
力分量(见图
8-2-1);
Pxi、Pzi—分别为结点i在x及z轴方向上的等效结点力;xi、zi—分别为结点i在x及z轴方向上的坐标值。
i+1
i1
x
地层
i
i
x
i
xz
i1
i1
i
x
xz
xz
i-1
i1i1
zz
i
z
图8-2-1等效结点力示意图
13
8.2.3.3绕结点平均法
1、绕结点平均法将围绕该结点的各单元的应力的平均值作为该结点的应力值,并假设
各结点间的应力呈线性分布,然后反转开挖边界上各点应力的方向,将其作为释放荷载。
2、结点应力值的计算式为:
1m
e
i
i
me1
(8-2-9)
式中:
i—结点i的平均应力值;
m—围绕该结点的全部单元的总数。
3、计算结点应力的平均值时,也可引入面积加权系数,即可将式(
8-2-9)改写为:
1
m
e
Ae]
m
i
[
i
Ae
me1
e1
(8-2-10)
式
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 隧道 地层 结构 计算 解读
![提示](https://static.bingdoc.com/images/bang_tan.gif)