寒区非圆形隧道冻胀力的解析解

第52卷第1期力学学报V ol.52,No.1 2020年1月Chinese Journal of Theoretical and Applied Mechanics Jan.,2020
生物、工程及交叉力学
寒区非圆形隧道冻胀力的解析解1)
李岩松∗,2)陈寿根†
∗(中国五冶集团有限公司,成都610063)
†(西南交通大学交通隧道工程教育部重点实验室,成都610031)
摘要利用复变函数相应理论引入一种求解考虑衬砌结构的季节性冻土地区非圆形隧道冻胀力和冻胀变形的方
法.方法克服了隧道断面形状为非圆形状并且同时考虑非圆形隧道支护、冻胀圈、未冻围岩时,经典的复变函数
理论不能直接应用求解非圆形隧道应力和位移的问题.方法通过将经典复变函数理论与连续性条件结合,导出
了非圆形隧道衬砌−冻胀圈−未冻围岩系统在正交曲线坐标系ζ平面内的解析式,然后通过保角变换求得直角
坐标系Z平面上考虑衬砌支护的寒区非圆形隧道冻胀力和冻胀变形.将推导的解析式应用于鹧鸪山隧道洞口段
研究中,得到鹧鸪山隧道洞口段冻胀应力和冻胀位移解析解,并将解析解与数值解进行对比,验证解析解的正确
性.由结果可知:围岩冻胀力对衬砌影响明显,在拱顶、拱脚、拱底处因冻胀造成的环向附加应力显著增大;在
拱脚及两侧拱肩处会受到较大的法向附加应力;由于衬砌几何形状的原因,造成了衬砌冻胀变形的不均匀,进而
造成冻胀力分布不均.应用复变函数理论将冻胀圈考虑为马蹄形、直墙拱形等非圆形状较以往将冻胀圈考虑为
单一圆环的冻胀力研究更符合实际,较以往研究更能反映实际工况.研究结果为季节性冻土地区非圆形隧道冻
胀力的弹塑性分析奠定了基础.
关键词季节性冻土隧道,非圆形衬砌,复变函数,冻胀力,解析解
中图分类号:TU91文献标识码:A doi:10.6052/0459-1879-19-226
ANALYTICAL SOLUTION OF FROST HEA VING FORCE IN NON-CIRCULAR
COLD REGION TUNNELS1)
Li Yansong∗,2)Chen Shougen†
∗(China MCC5Group Corp.Ltd,Chengdu610063,China)
†(MOE Key Laboratory of Transportation Tunnel Engineering,Southwest Jiaotong University,Chengdu610031,China)
Abstract A complex variable method is presented of stress and displacement problems for non-circular cold region tunnels.The complex variable method considers the non-circle tunnel support,the frost heaving circle and the unfrozen surrounding rock at the same time,making the problem from the single connected domain to multi connected domain problem,and can not directly apply the classical complex function theory to solve the problem of the stress and dis-placement of the non-circular tunnel.The analytical formula in Zeta plane orthogonal curvilinear coordinate system of non-circular tunnel lining to frozen surrounding rock to unfrozen surrounding rock system is derived by using classical
2019–08–20收稿,2019–12–16录用,2019–12–17网络版发表.
1)中央高校基本科研业务费专项资金项目(SWJTU11ZT33)和教育部创新团队发展计划(IRT0955)资助.
2)李岩松,博士研究生,主要从事隧道及岩石力学方面的研究.E-mail:******************
引用格式:李岩松,陈寿根.寒区非圆形隧道冻胀力的解析解.力学学报,2020,52(1):196-207
Li Yansong,Chen Shougen.Analytical solution of frost heaving force in non-circular cold region tunnels.Chinese Journal of Theoreti-
cal and Applied Mechanics,2020,52(1):196-207
第1期李岩松等:寒区非圆形隧道冻胀力的解析解197
complex variable method and continuity condition with power series and conformal transformation.Then,the frost heav-ing force and frost heaving deformation of the non-circular cold region tunnel are obtained by conformal transformation. This complex variable method is applied to Zhegushan tunnel research,and the analytical solution of Zhegushan tunnel openning frost heave stress and frost heave displacement is paring the analytical solution with the numeri-cal solution to verify the correctness of the analytical solution.It can be seen from the results that frost heaving force has an obvious influence on the lining,and the additional circumferential stresses caused by frost heaving at the vault,arch foot and arch bottom are significantly increased.The additional normal stresses caused by frost heaving at the vault and both sides of arch shoulder are significantly increased.Because of the geometric configurations of lining,the frost heaving deformation is uneven,which leads to uneven distribution of frost heaving pared with previous studies,which consider the frost heaving ring as a single ring,the complex variable method is more practical.And the research results were expected to provide a referenced basis for elasto-plastic analysis of cold-region tunnel.
Key words cold-region tunnel,non-circular lining,complex variables,frost force,analytic solution
引言
在高寒高海拔地区修建隧道,隧址区原有的热力平衡被打破,隧道的通风对流使冷空气在隧道内部的流动,造成了隧道内部温度的下降,从而创造了隧道围岩形成季节冻土或多年冻土的条件.当围岩温度过低围岩中的裂隙、孔隙水冻结时会发生体积膨胀,在这种体积膨胀受到隧道衬砌和未冻岩体的约束时,冻结围岩便会对隧道衬砌产生冻胀力.20世纪70∼80年代,日本对东北寒冷地区隧道的调查分析结果表明,冻胀力是造成寒区隧道变形破坏主要因素.由于我国经济的高速发展,高速铁路、公路建设迎来了高速发展时期,不可避免地需要在高寒高海拔地区修建隧道.隧道常常发生冻害,对隧道的维修养护造成了极大困难,更是对隧道结构和运营安全构成了严重威胁.为保证隧道结构和运营安全,降低维修养护成本,各国学者从不同方向、不同领域对寒区隧道冻胀力进行了研究.
Lai等[1-4]应用无量纲微扰法、Galerkin法、黏弹性理论以及CT技术分别对寒区隧道进行多项研究并取得了显著成果.为寒区隧道施工、运营、维护提供了有力的参考依据;Gao等[5]应用经典弹塑性理论得到了圆形寒区隧道围岩应力和塑性区,发现冻胀力会导致围岩塑性区显著膨胀并在屈服后出现塑性硬化现象;Feng等[6]应用弹塑性理论得到了圆形寒区隧道围岩应力及位移的弹塑性解析解,并且通过参数分析,得出相关参数对冻结圆周内、外边缘塑性半径和围岩应力的影响结果;夏才初等[7]考虑施工过程中机械产生的热量对多年冻土的融化影响,分析了多年冻土融化范围,建立并求解融化作用下多年冻土隧道围岩弹塑性模型,确定了相应控制融化作用的关键因素;张常光等[8]在弹性模量变化、中间主应力效应、围岩应变软化和剪胀条件下,基于统一强度理论推导了深埋圆形隧道塑性位移新解;谭贤君[9]、杨罡[10]、黄继辉等[11-12]、苏林军[13]应用有限元软件对冻土隧道冻胀力进行模拟,得到了一系列令人满意的结果;吴剑[14]通过模型试验结果与数值模拟结果对比,详细研究冻土隧道冻胀力,得出了不同冻胀条件下隧道衬砌变形特征曲线以及冻胀力分布规律;根据Takashi方程Hao等[15]得到了圆柱形冻土试块的冻胀力以及温度分布规律,论证了该方法在三维空间中的适用性;Tan等[16]应用数值模拟建立模型,计算了寒区隧道在多孔介质中冻融循环的温度扩散问题,通过数值模拟得到隧道衬砌表面的建筑保温材料会极大地改变衬砌和围岩温度的结论.
在以上寒区隧道研究中所涉及应用弹性、弹塑性理论计算寒区隧道冻胀力解析解时均采用圆形隧道模型,并未采用与实际隧道断面相似的隧道模型进行计算.由于计算模型的差异,这些研究成果与实际结果存在一定的差异,不能直接应用于非圆形寒区隧道问题中.
当考虑支护为圆形的条件时,寒区隧道冻胀力和位移的弹性解析解通过连续性条件易于给出.而当隧道断面形状为马蹄形、直墙拱形等非圆形状并且考虑隧道支护时,由于考虑了隧道支护,使问题由单连通域问题变为多连通域问题,不能直接应用经典的复变函数理论求解非圆形隧道应力和位移.所以,
198力学学报2020年第52卷
还鲜有学者给出考虑衬砌因素的非圆形寒区隧道冻
胀力和冻胀位移的弹性解析解.综上所述,为得到考
虑衬砌的非圆形寒区隧道冻胀力弹性解析解,应用
经典复变函数理论幂级数解法,结合连续性条件,在
文献[17–18]的基础上,推导了考虑衬砌的非圆形寒
区隧道冻胀力的弹性解析解.为下一步非圆形寒区
隧道冻胀力的弹塑性分析奠定了基础.
1问题描述
根据弹性理论,寒区隧道在深埋条件下,由于隧
道开挖尺寸与埋深相比很小,隧道开挖问题可视为在
无限平面内的孔洞问题,不受边界效应影响,地应力
变化梯度可忽略不计.当考虑冻胀圈和支护结构时,
无限平面内的孔洞问题由单连通域孔洞问题变为多
连通域问题,无限平面被分割为3个区域.图1给出
了无限平面内考虑冻胀圈的隧道孔洞模型.其中在
Z平面中衬砌、冻胀圈和未冻围岩分别表示为S1,S2
和S3区域,所对应的边界L1为衬砌与大气接触面,
它为自由边界,受到的面力为零;L2为衬砌与冻胀圈
接触面,它满足连续边界条件;L3为冻胀圈与未冻围
岩接触面,它满足连续边界条件;围岩的远场应力由
隧道埋深边界条件确定.为了进行计算,必须将Z平
面上的非圆形衬砌、冻胀圈结构通过保角变换映射
到ζ平面上,使得在Z平面上的非圆形结构变换为ζ
平面上的圆环结构.这时Z平面上的区域S3映射为
ζ平面行的区域O3;区域S2映射为ζ平面行的区域
O2;衬砌区域S1映射为圆环区域O1;边界L1,L2及
L3分别映射为圆环R1,R2与R3.其中最一般的映射
z-planeζ-plane
S
3
S
2
S
1L1L
2
L
3
O
3
O
2
O
1
R
1
R
2
R
3
x
y O O'
α
ζ
ρ =1
ρ =r
η
图1Z平面保角变换为ζ平面示意图
Fig.1Conformal mapping of tunnel in z-plane in to two
concentric circles inζ-plane 函数可写成洛朗级数的形式,即
w(ζ)=R
ζ+
∞
∑
k=0
C kζ−k
(1)
其中,R是一个实常数,它的大小反映了隧道孔洞的大小;C k一般为复常数,其中k=1,2,3,···,C k反映了隧道断面形状,在一般情况下C k取前几项时,映射函数的精度已经可以满足工程要求[19].ζ为Z平面上的点映射到ζ平面上的正交曲线坐标,通常可由角度θ与半径ρ表示.
2复变函数解法公式推导及求解
根据Kargar等[20-22]和陈子荫[23]提出的方法,ζ面上未冻围岩O3区域内的应力函数可表示为
φ1(ζ)=Γw(ζ)+φ0(ζ)(2)
ψ1(ζ)=Γ′w(ζ)+ψ0(ζ)(3)
其中
φ0(ζ)=
∞
∑
k=0
h kζ−k(4)
ψ0(ζ)=
∞
∑
k=0
m kζ−k(5)式中,h k和m k为所求复常数,它们均由边界条件确定;φ0(ζ)和ψ0(ζ)满足φ0(∞)=0和ψ0(∞)=0;Lu 等[21-22]提出Γ和Γ′可由无穷远处应力边界条件定出,即
Γ=1+K
4
γH(6)
Γ′=−(1−K)γH(7)
在方程(6)和(7)中K为侧压力系数,γ和H分别表征围岩重度与隧道埋深.
ζ平面上冻胀圈O2区域的应力函数φ2(ζ)和ψ2(ζ)可表示为
φ2(ζ)=
∞
∑
k=0
a kζk+
∞
∑
k=1
b kζ−k(8)
ψ2(ζ)=
∞
∑
k=0
c kζk+
∞
∑
k=1
d kζ−k(9)其中,a k,b k,c k和d k为所求复常数,它们均由边界条件确定.
第1期李岩松等:寒区非圆形隧道冻胀力的解析解199
ζ平面上衬砌O1区域应力函数φ3(ζ)和ψ3(ζ)
可表示为
φ3(ζ)=
∞
∑
k=0
A kζk+
∞
∑
k=1
B kζ−k(10)
ψ3(ζ)=
∞
∑
k=0
C kζk+
∞
∑
k=1
D kζ−k(11)
其中,A k,B k,C k和D k为所求复常数,它们均由边界条件确定.
在复变函数解中[24-26]应力均可以用两个应力函数φ和ψ表达,衬砌、冻胀圈和未冻围岩上应力均可表示为
σθ+σρ=4Re [
φ′(ζ)
w′(ζ)
]
(12)
σθ−σρ+2iτρθ=2ζ2
ρ2
1 w′(ζ)
[
w(ζ)φ′′(ζ)w′(ζ)−φ′(ζ)w′′(ζ)
w′(ζ)w′(ζ)
+ψ′(ζ)
]
(13)
其中,σρ,σθ和τρθ分别为正交曲线坐标下的法向应力、环向应力和剪应力.
假设寒区隧道衬砌S1和未冻围岩S3区域体积不受围岩温度变化影响,相应ζ平面上映射区域O1和O3位移满足复变函数[24-26]应力函数表达,衬砌和未冻围岩上位移可表示为
2G (
uρ+i uθ
)
=
ζ
ρ
w′(ζ)
|w(ζ)|
κφ(ζ)−w(ζ)
w′(ζ)
φ′(ζ)−ψ(ζ)
(14)
其中,uρ和uθ分别表示法向位移与环向位移,G是切变模量,κ由泊松比确定,κ=3−4υ和κ= (3−υ)/(1+υ)分别为平面应变与平面应力状态下κ的取值,υ为泊松比,因隧道纵向长度与断面最大直径相比很大,故问题可视为平面应变问题,κ取3−4υ.
由于冻胀圈S2区域内体积受温度变化影响,故相应在ζ平面上O2区域内位移无法满足复变函数[24-26]应力函数表达,需要重新推导(推导过程详见附录1),经推导冻胀圈上位移可表示为
2G (
uρ+i uθ
)
=
ζ
ρ
w′(ζ)
|w(ζ)|
κφ(ζ)−w(ζ)
w′(ζ)
φ′(ζ)−ψ(ζ)
+
2G
ζ
ρ
w′(ζ)
|w(ζ)|w(ζ)ε(15)
其中ε为冻胀圈线膨胀系数,可由冻胀区单位冻胀体
积变化定出,即
ε=1
3
∆V
其中∆V是冻胀区单位冻胀体积变化.
应力函数φ1(ζ),ψ1(ζ),φ2(ζ)和ψ2(ζ),φ3(ζ)和
ψ3(ζ)必须分别在R2和R3上满足连续边界条件,在
R1上满足应力边界条件.在衬砌与大气接触面上力
的边界条件可表示为
φ3(t1)+w(t1)
w′(t1)
φ′
3
(t1)+ψ3(t1)=0(16)
在计算非圆形隧道围岩应力时,为能更好地反
应实际隧洞变形位移,Lu等[20-22]在连续性位移条
件中引入了开挖卸荷因子.为了能更直观反映冻胀
力对非圆形隧道衬砌的影响,在位移连续边界条件
中并未考虑开挖卸荷因子,如需考虑开挖卸荷效应,
只需在位移边界条件中添加符合实际工况的开挖卸
荷因子即可.根据推导的位移公式,在冻胀圈与衬砌
接触面上位移连续边界条件可表示为
κ2
2G2
φ2(t2)−1
2G2
(
w(t2)
w′(t2)
φ′
2
(t2)+ψ2(t2)
)
+w(ζ)ε=
κ3
2G3
φ3(t2)−1
2G3
(
w(t2)
w′(t2)
φ′
3
(t2)+ψ3(t2)
)
(17)
在冻胀圈与衬砌接触面上力的连续边界条件可
表示为
φ2(t2)+w(t2)
w′(t2)
φ′
2
(t2)+ψ2(t2)=
φ3(t2)+w(t2)
w′(t2)
φ′
2
(t3)+ψ3(t2)(18)
在未冻围岩与冻胀圈接触面上,根据推导的考
虑冻胀作用的位移公式(附录1),位移连续边界条件
可表示为
κ1
2G1
φ0(t3)−1
2G1
(
w(t3)
w′(t3)
φ′
(t3)+ψ0(t3)
)
=
κ2
2G2
φ2(t3)−1
2G2
(
w(t3)
w′(t3)
φ′
2
(t3)+ψ2(t3)
)
+w(ζ)ε(19)
在未冻围岩与冻胀圈接触面上力的连续边界条
件可表示为
φ1(t3)+w(t3)
w′(t3)
φ′
1
(t3)+ψ1(t3)=
φ2(t3)+w(t3)
w′(t3)
φ′
2
(t3)+ψ2(t3)(20)
其中,t1,t2和t3分别表示ζ平面上R1,R2和R3边界
上的点.
200力学学报2020年第52卷
方程(2)和(3)中的表达式Γw(ζ)和Γ′w(ζ)表示围岩的初始地应力和初始位移的分量,所以在方程(19)中计算由于开挖引起的位移时应将其去除[27];方程(17)∼(20)在计算过程中假设接触面之间无滑动.
应用于李岩松等[17-18]幂级数解法,方程(16)可分别展开为正幂次方程组(21)和负幂次方程组(22)
∞
∑k=0A kρ2k−v(k−v+1)L k−v+1σv−
∞
∑
k=1
B kρ−2k−v(k+v−1)L−k−v+1σv+
∞
∑
k=1
A kρ2k+v−2kL k+v−1σv−
∞
∑
k=1
B kρ−2k+v−2kL−k+v−1σv−
∞
∑
k=0
C kρ−v(k+v−1)L−k−v+1σv+
∞
∑
k=1
D kρ−v(k−v+1)L k−v+1σv=0,
v=1,2,3, (21)
∞
∑k=0A kρ2k+v(k+v+1)L k+v+1σ−v−
∞
∑
k=1
B kρ−2k+v(k−v−1)L−k+v+1σ−v+
∞
∑
k=1
A kρ2k−v−2kL k−v−1σ−v−
∞
∑
k=1
B kρ−2k−v−2kL−k−v−1σ−v−
∞
∑
k=0
C kρv(k−v−1)L−k+v+1σ−v+
∞
∑
k=1
D kρv(k+v+1)L k+v+1σ−v=0,
v=1,2,3, (22)
其中,L v可由方程(1)中R和C k计算得出,ζ为Z平面上的点映射到ζ平面上的正交曲线坐标,可表示为ζ=ρσ,这里σ=exp(iθ);ζ的共轭可表示为ζ=ρσ−1.为便于计算,方程(1)可写为
w(ζ)=
+∞
∑
v=−∞
L vζ−v(23)
方程(17)可分别展开为正幂次方程组(24)和负幂次
方程组(25);方程(18)可分别展开为正幂次方程组
(26)和负幂次方程组(27);方程(19)可分别展开为正
幂次方程组(28)和负幂次方程组(29),方程(24)∼(29)
表达式如下
κ2 G2
∞
∑
k=0
a kρ2k−v(k−v+1)L k−v+1σv−
κ2
G2
∞
∑
k=1
b kρ−2k−v(k+v−1)L−k−v+1σv−
1
G2
∞
∑
k=1
a kρ2k+v−2kL k+v−1σv+
1
G2
∞
∑
k=1
b kρ−2k+v−2kL−k+v−1σv+
1
G2
∞
∑
k=0
c kρ−v(k+v−1)L−k−v+1σv−
1
G2
∞
∑
k=1
d kρ−v(k−v+1)L k−v+1σv−
κ3
G3
∞
∑
k=0
A kρ2k−v(k−v+1)L k−v+1σv+
κ3
G3
∞
∑
k=1
B kρ−2k−v(k+v−1)L−k−v+1σv+
1
G3
∞
∑
k=1
A kρ2k+v−2kL k+v−1σv−
1
G3
∞
∑
k=1
B kρ−2k+v−2kL−k+v−1σv−
1
G3
∞
∑
k=0
C kρ−v(k+v−1)L−k−v+1σv+
1
G3
∞
∑
k=1
D kρ−v(k−v+1)L k−v+1σv=
−ε
∞
∑
k=0
L kρ2k−v(k−v+1)L k−v+1σv+ε
∞
∑
k=1
L−kρ−2k−v(k+v−1)L−k−v+1σv,v=1,2,3, (24)
κ2 G2
∞
∑
k=0
a kρ2k+v(k+v+1)L k+v+1σ−v−
k2
G2
∞
∑
k=1
b kρ−2k+v(k−v−1)L−k+v+1σ−v−
1
G2
∞
∑
k=1
a kρ2k−v−2kL k−v−1σ−v+
1
G2
∞
∑
k=1
b kρ−2k−v−2kL−k−v−1σ−v+
1
G2
∞
∑
k=0
c kρv(k−v−1)L−k+v+1σ−v−
1
G2
∞
∑
k=1
d kρv(k+v+1)L k+v+1σ−v−
κ3
G3
∞
∑
k=0
A kρ2k+v(k+v+1)L k+v+1σ−v+
k3
G3
∞
∑
k=1
B kρ−2k+v(k−v−1)L−k+v+1σ−v+
1
G3
∞
∑
k=1
A kρ2k−v−2kL k−v−1σ−v−1
G3
∞
∑
k=1
B kρ−2k−v−2kL−k−v−1σ−v−
1
G3
∞
∑
k=0
C kρv(k−v−1)L−k+v+1σ−v+
1
G3
∞
∑
k=1
D kρv(k+v+1)L k+v+1σ−v=
−ε
∞
∑
k=0
L kρ2k+v(k+v+1)L k+v+1σ−v+ε
∞
∑
k=1
L−kρ−2k+v(k−v−1)L−k+v+1σ−v,v=0,1,2,3, (25)
∞
∑k=0a kρ2k−v(k−v+1)L k−v+1σv−
∞
∑
k=1
b kρ−2k−v(k+v−1)L−k−v+1σv+
∞
∑
k=1
a kρ2k+v−2kL k+v−1σv−
第1期李岩松等:寒区非圆形隧道冻胀力的解析解201∞
∑k=1b kρ−2k+v−2kL−k+v−1σv−
∞
∑
k=0
c kρ−v(k+v−1)L−k−v+1σv+
∞
∑
k=1
d kρ−v(k−v+1)L k−v+1σv−
∞
∑k=0A kρ2k−v(k−v+1)L k−v+1σv+
∞
∑
k=1
B kρ−2k−v(k+v−1)L−k−v+1σv−
∞
∑
k=1
A kρ2k+v−2kL k+v−1σv+
∞
∑k=1B kρ−2k+v−2kL−k+v−1σv+
∞
∑
k=0
C kρ−v(k+v−1)L−k−v+1σv−
∞
∑
k=1
D kρ−v(k−v+1)L k−v+1σv=0,
v=0,1,2,3, (26)
∞
∑k=0a kρ2k+v(k+v+1)L k+v+1σ−v−
∞
∑
k=1
b kρ−2k+v(k−v−1)L−k+v+1σ−v+
∞
∑
k=1
a kρ2k−v−2kL k−v−1σ−v−
∞
∑
k=1
b kρ−2k−v−2kL−k−v−1σ−v−
∞
∑
k=0
c kρv(k−v−1)L−k+v+1σ−v+
∞
∑
k=1
d kρv(k+v+1)L k+v+1σ−v−
∞
∑
k=0
A kρ2k+v(k+v+1)L k+v+1σ−v+
∞
∑
k=1
B kρ−2k+v(k−v−1)L−k+v+1σ−v−
∞
∑
k=1
A kρ2k−v−2kL k−v−1σ−v+
∞
∑
k=1
B kρ−2k−v−2kL−k−v−1σ−v+
∞
∑
k=0
C kρv(k−v−1)L−k+v+1σ−v−
∞
∑
k=1
D kρv(k+v+1)L k+v+1σ−v=0,
v=0,1,2,3, (27)
κ2 G2
∞
∑
k=0
a kρ2k−v(k−v+1)L k−v+1σv−
k2
G2
∞
∑
k=1
b kρ−2k−v(k+v−1)L−k−v+1σv−
1
G2
∞
∑
k=1
a kρ2k+v−2kL k+v−1σv+
1
G2
∞
∑
k=1
b kρ−2k+v−2kL−k+v−1σv+
1
G2
∞
∑
k=0
c kρ−v(k+v−1)L−k−v+1σv−
1
G2
∞
∑
k=1
d kρ−v(k−v+1)L k−v+1σv+
κ1
G1
∞
∑
k=0
h kρ−2k−v(k+v−1)L−k−v+1σv−
1
G1
∞
∑
k=0
h kρ−2k+v−2kL−k+v−1σv+
1
G1
∞
∑
k=0
m kρ−v(k−v+1)L k−v+1σv=
−ε
∞
∑
k=0
L kρ2k−v(k−v+1)L k−v+1σv+ε
∞
∑
k=1
L−kρ−2k−v(k+v−1)L−k−v+1σv,v=0,1,2,3, (28)
κ2 G2
∞
∑
k=0
a kρ2k+v(k+v+1)L k+v+1σ−v−
κ2
G2
∞
∑
k=1
b kρ−2k+v(k−v−1)L−k+v+1σ−v−
1
G2
∞
∑
k=1
a kρ2k−v−2kL k−v−1σ−v+
1
G2
∞
∑
k=1
b kρ−2k−v−2kL−k−v−1σ−v+
1
G2
∞
∑
k=0
c kρv(k−v−1)L−k+v+1σ−v−
1
G2
∞
∑
k=1
d kρv(k+v+1)L k+v+1σ−v+
κ1
G1
∞
∑
k=0
h kρ−2k+v(k−v−1)L−k+v+1σ−v−
1
G1
∞
∑
k=0
h kρ−2k−v−2kL−k−v−1σ−v+
1
G1
∞
∑
k=0
m kρv(k+v+1)L k+v+1σ−v=−ε
∞
∑
k=0
L kρ2k+v(k+v+1)L k+v+1σ−v+ε
∞
∑
k=1
L−kρ−2k+v(k−v−1)L−k+v+1σ−v,v=0,1,2,3, (29)
最后将方程(20)分别展开为正幂次方程组(30)和负幂次方程组(31).应当注意的是:由于边界条件不同,方程(21)和(22)中的极半径ρ与方程(24)∼(27)以及方程(28)∼(31)中的极半径ρ在数值上分别等于r1,1和r2.
结合行列式(21)、式(22)以及行列式(24)∼式(31),可以得到无穷线性方程组[17-18],求解无穷线性方程组得到应力函数(方程(2)和(3),方程(8)和(9),方程(10)和(11))中各项系数.由得到的应力函数计算方程(12)和(13)得到非圆形寒区隧道应力解析解;由方程(14)计算衬砌和未冻围岩位移解析解;最后由方程(15)计算冻胀圈位移解析解.
∞
∑k=0a kρ2k−v(k−v+1)L k−v+1σv−
∞
∑
k=1
b kρ−2k−v(k+v−1)L−k−v+1σv+
∞
∑
k=1
a kρ2k+v−2kL k+v−1σv−∞
∑
k=1
b kρ−2k+v−2kL−k+v−1σv−
∞
∑
k=0
c kρ−v(k+v−1)L−k−v+1σv+
∞
∑
k=1
d kρ−v(k−v+1)L k−v+1σv+
202力学学报2020年第52卷∞
∑k=0h kρ−2k−v(k+v−1)L−k−v+1σv+
∞
∑
k=0
h kρ−2k+v−2kL−k+v−1σv−
∞
∑
k=0
m kρ−v(k−v+1)L k−v+1σv=
2Γ
∞
∑
k=0
L kρ2k−v(k−v+1)L k−v+1σv−2Γ
∞
∑
k=1
L−kρ−2k−v(k+v−1)L−k−v+1σv−
Γ′
∞
∑
k=0
L kρ−v(k+v−1)L−k−v+1σv+Γ′
∞
∑
k=1
L−kρ−v(k−v+1)L k−v+1σv,v=0,1,2,3, (30)
∞
∑k=0a kρ2k+v(k+v+1)L k+v+1σ−v−
∞
∑
k=1
b kρ−2k+v(k−v−1)L−k+v+1σ−v+
∞
∑
k=1
a kρ2k−v−2kL k−v−1σ−v−
∞
∑
k=1
b kρ−2k−v−2kL−k−v−1σ−v−
∞
∑
k=0
c kρv(k−v−1)L−k+v+1σ−v+
∞
∑
k=1
d kρv(k+v+1)L k+v+1σ−v+
∞
∑
k=0
h kρ−2k+v(k−v−1)L−k+v+1σ−v+
∞
∑
k=0
h kρ−2k−v−2kL−k−v−1σ−v−
∞
∑
k=0
m kρv(k+v+1)L k+v+1σ−v=
2Γ
∞
∑
k=0
L kρ2k+v(k+v+1)L k+v+1σ−v−2Γ
∞
∑
k=1
L−kρ−2k+v(k−v−1)L−k+v+1σ−v−
Γ′
∞
∑
k=0
L kρv(k−v−1)L−k+v+1σ−v+Γ′
∞
∑
k=1
L−kρv(k+v+1)L k+v+1σ−v,v=0,1,2,3, (31)
3工程案例应用分析
以新建汶川至马尔康公路鹧鸪山特长隧道洞口段为研究对象,利用上节解析方法,对非圆形寒区隧道应力和变形进行求解,并用Flac有限差分软件进行进行校核,验证解析解的准确性.
根据现场大气、围岩实测温度、围岩及结构热力学试验,估算围岩冻结深度.可知当原岩温度为10◦C 时,冻结深度约为2m[28].据此确定围岩冻胀圈和未冻围岩.
结合鹧鸪山隧道地质报告、室内岩石试验和点荷载试验数据[28].围岩、结构材料物理力学参数见表1.
表1地层及材料参数
Table1Main physical parameter for tunnel calculation
Elasticity Poisson’s Severe/Coefficient
Water
Materials
E/MPa ratio v(kN·m−3)
of lateral
content/% pressure K
slate86000.29200.50.41
lining300000.2250.5 2.85 frozen soil99000.27200.58.26
图2(a)为鹧鸪山隧道K188+480m处断面计算简图.图2(b)为Flac有限差分软件计算模型图,共划分网格3916个,用于模拟深埋非圆形寒区隧道开挖后衬砌、冻胀圈和围岩的应力及变形状态.隧道采用全断面开挖,开挖后立即支护;假定模型为弹性模型;衬砌与冻胀圈,冻胀圈与围岩之间只受法向力;隧道冻胀系数采用等效热胀系数进行模拟;假设冻胀圈温度恒定,无温差;隧道冻胀力和冻胀变形通过控制热力模型和结构模型交替开关进行计算;衬砌与冻胀圈,冻胀圈与围岩之间无滑动.根据式(15)中对线胀系数的说明、水冻结为冰体积增大约9%及冻胀圈含水量计算得冻胀圈线膨胀系数ε≈0.0025.
根据吕爱钟等[19]提出的最优化法求解映射函数的思路,编写相应MATLAB计算程序,通过计算,得到映射函数
w(ζ)=8.0339(ζ+0.3121−0.0697ζ−1+
0.0338ζ−2−0.0087ζ−3−0.0068ζ−4)(32)
在工程计算中方程(1)取前几项就可满足工程精度要求[19],所以为提高计算效率,式(37)取映射函数(方程(1))的前四项;式(37)中ζ=ρσ,其中当ρ=1.225时,表示映射函数所对应未冻围岩与冻胀圈接触面的边界条件;当ρ=1时,表示映射函数所对应冻胀圈与衬砌接触面的边界条件;当ρ=0.8633时,表示映射函数所对应衬砌与大气接触面的边界条件.
第1期李岩松等:寒区非圆形隧道冻胀力的解析解203 x
R9.9
R2.9R
1
6
.9
R7
.43
50
12.05
10.05
9.05
y
z
x
(b)
(a)
图2(a)隧道计算简图(单位:m);(b)隧道计算模型
Fig.2(a)Tunnel calculation diagram(Unit:m)and
(b)Finite mesh calculation mode
根据边界条件(方程(16)∼(20))编写MATLAB 计算程序,求得寒区非圆形隧道冻胀应力及冻胀位移解析解,并将其与Flac有限差分软件所求的数值解进行比较,比较结果如下所示.
由于Flac有限差分软件计算应力为隧道衬砌及冻胀圈的竖向和水平向应力,这与解析解所求法向、环向应力坐标系统不统一,无法进行比较,故将Flac 有限差分软件所求应力结果进行坐标转换后再进行对比.
图3是衬砌内边界环向应力,由图中可看到,数值解与解析解吻合得较好;由于冻胀作用的影响衬砌内边界环向应力显著增大,在拱顶(α=0◦)到拱肩(α=40◦)这一扇形区域衬砌所受的环向应力明显增大;在拱脚(α=120◦)处,环向应力集中明显,拱脚处环向应力较未冻结时增大约25.9%;在拱底(α=180◦)附近30◦的扇形区域范围内衬砌所受环向应力增长显著.
02040
5
4
3
2
1
6080100
analytical solution of frozen
analytical solution of unfrozen
numerical solution of frozen
numerical solution of unfrozen
120140160180200
α/(Ο)
σ
θ
/
M
P
a
图3衬砌内边界环向应力
Fig.3Circumferential stress along the inner lining periphery
图4为衬砌与冻胀圈接触面上法向应力,由图4可以看出:衬砌法向应力较未受冻胀力时有明显波动.由于冻胀作用,衬砌在α=60◦附近法向应力较未受冻胀力时增大约107.5%;在α=110◦附近冻胀后衬砌与未冻胀衬砌相比法向应力减小约82.4%;在α=140◦附近受冻胀力影响法向应力较未冻胀法向应力增大约21.4%;在衬砌与冻胀圈接触面上,数值模拟所求法向应力与解析法所求法向应力吻合得较好,只是在部分路径上有一定的偏差.
02040
5
4
3
2
1
6080100120140160180200
α/(Ο)
σ
ρ
/
M
P
a
analytical solution of frozen
numerical solution of frozen
analytical solution of unfrozen
numerical solution of unfrozen
图4衬砌与冻胀圈接触面上法向应力
Fig.4Normal stress along frost heaving ring-lining interface
由图5可见,冻胀作用对衬砌影响明显,尤其在衬砌拱顶、拱脚、拱底处.由于冻胀力影响,拱顶环
204力学学报2020年第52卷
向应力由2.01MPa 增大至4.55MPa;拱脚处环向应力较未冻胀环向应力增大约38.8%;拱底环向应力由
0.5MPa 增大至2.5MPa;在衬砌外边界上,数值模拟
所求环向应力与解析法所求环向应力吻合得较好,只是在α=100◦应力峰值附近有一定的偏差.
20
40
5
4
3
2
1
60
80
100120140160180200
α/(Ο)
σθ/M P a
analytical solution of frozen numerical solution of frozen
analytical solution of unfrozen numerical solution of unfrozen
图5衬砌外边界环向应力
Fig.5Circumferential stress along outer lining periphery
图6为衬砌与冻胀圈接触面上的法向位移,从图中可以看到:数值法向位移解与解析法向位移解吻合得较好;衬砌受到冻胀力,产生了明显的冻胀变形.衬砌拱顶下沉较未受冻胀力时增大约37.4%.拱底隆起明显,与未受冻胀力时相比增大约3mm.但是,总的来说计算所得未冻胀和冻胀法向位移均较实测值偏小[28],这可能与未考虑开挖卸荷因子和塑性变形有关.
20
3210-1-2-3-4
40
60
80
100120140160180200α/(Ο)
u r /m m
analytical solution of frozen
analytical solution of unfrozen numerical solution of frozen numerical solution of unfrozen
图6衬砌内边界上法向位移
Fig.6Radial displacement along the inner lining periphery
4结论
根据复变函数理论,首次提出了考虑衬砌的非圆形寒区隧道冻胀力和变形方法,并将其应用于鹧鸪山隧道洞口段研究中,得出了以下几点结论:
(1)推导了衬砌−冻胀圈−围岩体系的非圆形隧
道冻胀力弹性解析公式,然后将复变函数解法得到
的解析解与数值分析结果进行对比分析,评价考虑围岩冻胀作用的非圆形隧道应力及位移解析解的合理性和准确性.
(2)解析解与数值解吻合得较好,围岩冻胀力对
衬砌影响明显,在拱顶、拱脚、拱底处因冻胀造成的
环向附加应力显著增大;在拱脚及两侧拱肩处受到较大的法向附加应力;围岩冻胀力造成衬砌受力不均,对隧道的结构安全构成了威胁.
(3)以新建汶川至马尔康高速公路鹧鸪山特长隧
道洞口段为研究对象,结合工程实际,利用MATLAB 计算软件编写计算程序求解,得到鹧鸪山隧道洞口段冻胀应力和冻胀位移解析解,并将所得结果作为其数值计算和安全运行的参考依据.
(4)由于衬砌几何结构的原因造成了隧道冻胀变
形的不均匀,进而造成冻胀力分布不均,这较以往将围岩冻胀圈考虑为单一圆环的围岩冻胀力研究更符合实际.较以往所得解析解更能反映实际工况.
(5)复变函数理论为寒区非圆形隧道冻胀力研究
提出了一种除数值模拟、现场试验外新的计算方法.新方法为下一步寒区非圆形隧道冻胀力的弹塑性分析奠定了基础.
参
考
文
献
1Lai YM,Liu SY ,Wu ZW,et al.Approximate analytical solution for
temperature fields in cold regions circular tunnels.Cold Reg.Sci.Technol ,2002,34:43-49
2Zhang XF,Lai YM,Yu WB,et al.Non-linear analysis for the freezing–thawing situation of the rock surrounding the tunnel in cold regions under the conditions of di fferent construction seasons,initial temperatures and insulations.Tunnelling Underground Space Tech-nol ,2002,17:315-325
3Zhang SJ,Lai YM,Zhang XF,et al.Study on the damage propaga-tion of surrounding rock from a cold-region tunnel under freeze–thaw cycle condition.Tunnelling Underground Space Technol ,2004,19:295-302
4Lai YM,Wu H,Wu ZW,et al.Analytical viscoelastic solution for frost force in cold-region tunnels.Cold Reg.Sci.Technol.,2000,31:227-234
5Gao GY ,Chen QS,Zhang QS,et al.Analytical elasto-plastic solu-tion for stress and plastic zone of surrounding rock in cold region tunnels.Cold Reg.Sci.Technol.,2012,72:50-57
第1期李岩松等:寒区非圆形隧道冻胀力的解析解205
6Feng Q,Zhang BS,Qiang JG,et al.Analytical elasto-plastic solu-tion for stress and deformation of surrounding rock in cold region tunnels.Cold Reg.Sci.Technol.,2014,108:59-68
7夏才初,黄继辉,卞跃威等.融化作用下多年冻土隧道围岩的弹塑性解及其与支护的相互作用分析.岩土力学,2013,34(7):1987-1994(Xia Caichu,Huang Jihui,Bian Yuewei,et al.Elastoplastic analysis of surrounding rock of permafrost tunnel with thawing ef-fect and its interaction with support.Rock and Soil Mechanics,2013, 34(7):1987-1994(in Chinese))
8张常光,徐飞,张庆贺等.岩石隧道塑性位移新解.岩石力学与工程学报,2011,30(增2):3551-3556(Zhang Changguang,Xu Fei, Zhang Qinghe,et al.New plastic solution for displacement of rock tunnles.Chinese Journal of Rock Mechanics and Engineering,2011, 30(supp.2):3551-3556(in Chinese))
9谭贤君.高海拔寒区隧道冻胀机理及其保温技术研究.[博士论文].武汉:中国科学院武汉岩土力学研究所,2010(Tan Xian-jun.Study on the meehanism of frost heave of tunnel in cold region with high altitude and related insulation technology.[PhD Thesis].
Wuhan:Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,2010(in Chinese))
10杨罡.兰新高铁祁连山隧道围岩温度场及冻胀力研究.[硕士论文].西安:西安科技大学,2017(Yang Gang.Study on surrounding rock temperaturefield and frost heave force of Qilian mountain tun-nel of Lanxin railway.[Master Thesis].Xi’An:Xi’an University of Science and Technology,2017(in Chinese))
11黄继辉,夏才初,韩常领等.寒区公路隧道冻胀力荷载的分布形式和简化计算方法.现代隧道技术,2016,53(5):63-70(Huang Jihui,Xia Caichu,Han Changling,et al.The distribution form and simplified calculation method of frost heaving force load in highway tunnels in cold regions.Modern Tunneling Technology,2016,53(5): 63-70(in Chinese))
12黄继辉,夏才初,韩常领等.考虑围岩不均匀冻胀的寒区隧道冻胀力解析解.岩石力学与工程学报,2015,34(增2):3766-3774 (Huang Jihui,Xia Caichu,Han Changling,et al.Analytical solu-tion of frost heave force acting on cold-region tunnel liner consider-ing anisotropy frost heave of surrounding rock.Chinese Journal of Rock Mechanics and Engineering,2015,34(supp.2):3766-3774(in Chinese))
13苏林军.寒区隧道冻害预测与对策研究.[硕士论文].成都:西南交通大学,2007(Su Linjun.Study on forcasting of frost damage and antifreezing strategies of tunnel in cold area.[Master Thesis].
Chengdu:Southwest Jiaotong University,2007(in Chinese))
14吴剑.随道冻害机理及冻胀力计算方法的研究.[硕士论文].成都:西南交通大学,2004(Wu Jian.Study on principle of tunnel freeze injure and calculation method of frost force.[Master Thesis].
Chengdu:Southwest Jiaotong University,2004(in Chinese))
15Zheng H,Shunji KN,Niu FJ,et al.Three-dimensional frost heave evaluation based on practical Takashi’s equation.Cold Regious Sci-ence&Technology,2015,118:30-37
16Tan XJ,Chen WZ,Wu GJ,et al.Numerical simulations of heat transfer with ice–water phase change occurring in porous media and application to a cold-region tunnel.Tunnelling Underground Space Technol.，2013,38:170-179
17李岩松,陈寿根.基于复变函数理论的非圆形隧道解析解.西南交通大学学报,2019,DOI:10.3969/j.issn.0258-2724.20180243(Li Yansong,Chen Shougen.Stress and deformation study of lining and surrounding rock in a non-circular tunnel based on the complex variable method Seepage.Journal of Southwest Jiaotong University, 2019,DOI:10.3969/j.issn.0258-2724.20180243(in Chinese))18李岩松,陈寿根.考虑注浆加固作用下非圆形隧道应力和变形解析解.中国公路学报,2018,31(10):255-265(Li Yansong,Chen Shougen.Analytic solution on stress and deformation of non-circular tunnel considering grouting circle.China Journal of High-way and Transport,2018,31(10):255-265(in Chinese))
19吕爱钟,张路青.地下隧洞力学分析的复变函数方法.北京:科学出版社,2007(L¨u Aizhong,Zhang plex Function Method for Mechanical Analysis of Underground Tunnel.Beijing: Science Press,2007(in Chinese))
20Kargar AR,Rahmannejad R,Hajabasi MA.A semi-analytical elastic solution for stressfield of lined non-circular tunnels at great depth using complex variable method.International Journal of Solids and Structures,2014,51:1475-1482
21Lu AZ,Zhang N,Kuang L.Analytic solutions of stress and displace-ment for a non-circular tunnel at great depth including support delay.
International Journal of Rock Mechanics&Mining Sciences,2014, 70:69-81
22Lu AZ,Zhang N,Qin Y.Analytical solutions for the stress of a lined non-circular tunnel under full-slip contact conditions.International Journal of Rock Mechanics&Mining Sciences,2015,79:183-192 23陈子荫.围岩力学分析中的解析方法.北京:煤炭工业出版社, 1994(Chen Ziyin.Analytical Method in Mechanical Analysis of Surrounding Rock.Beijing:China Coal Industry Publishing House, 1994(in Chinese))
24张治国,杨轩,宫剑飞等.复变函数法分析盾构隧道开挖引起的土体位移和衬砌变形.岩土工程学报,2017,39(9):1626-1635 (Zhang Zhiguo,Yang Xuan,Gong Jianfei,et al.A complex vari-able analysis of soil displacement and liner deformation induced by shield excavation.Chinese Journal of Geotechnical Engineering, 2017,39(9):1626-1635(in Chinese))
25Li YS,Chen SG.A complex variable solution for lining stress and deformation in a non-circular deep tunnel II practical application and verification.Mathematical and Computational Applications,2018, 23,43,doi:10.3390/mca23030043
26覃爱民,骆汉宾.基于Mohr-Coulomb准则寒区隧道围岩应力弹塑性解析.地下空间与工程学报,2018,14(2):395-402(Tan Aimin,Luo Hanbin.Elastoplastic analysis of stress in tunnel sur-rounding rock in cold region based on mohr-coulomb criterion.Chi-nese Journal of Underground Space and Engineering,2018,14(2): 395-402(in Chinese))
27Kargar AR,Rahmannejad R,Hajabasi MA.A semi-analytical elastic solution for stressfield of lined non-circular tunnels at great depth using complex variable method.International Journal of Solids and Structures,2014,51:1475-1482
28赵玉报.高原高寒隧道围岩冻胀行为及施工对策研究.[博士论文].成都:西南交通大学,2011(Zhao Yubao.Study on the frost heaving behavior and construction countermeasure for plateau tun-nel in cold region after earthquake.[PhD Thesis].Chengdu:South-west Jiaotong University,2011(in Chinese))
29黄诗冰,刘泉声,刘艳章等.低温热力耦合下岩体椭圆孔(裂)隙中冻胀力与冻胀开裂特征研究.岩土工程学报,2018,40(3):459-467 (Huang Shibing,Liu Quansheng,Liu Yanzhang,et al.Frost heav-ing pressure and characteristics of frost cracking in elliptical cavity (crack)of rock mass under coupled thermal-mechanical condition at low temperature.Chinese Journal of Geotechnical Engineering, 2018,40(3):459-467(in Chinese))
30夏明.寒区隧道温度场数值分析研究.[硕士论文].淮南:安徽理工大学,2018(Xia Ming.Numerical analysis of tunnel temperature field in cold area.[Master Thesis].Huainan:Anhui University of Science&Technology,2018(in Chinese))。