第33卷第28期西建筑 Vol. 33No. 28 山2007年10月Oct. 2007SHANXI ARCH ITECTURE
#327#
文章编号:1009-6825(2007) 28-0327-02
地下拱形结构内力分析
张 林
摘 要:根据局部变形地基梁法, 利用改进了的压力拱理论计算地层压力, 对212国道某一隧道衬砌(变截面直墙拱) 进行了内力计算, 同时利用大型有限元进行模拟, 分析了拱形结构的受力特点, 对隧道工程的设计与施工都具有一定的指导作用。
关键词:地下拱形结构, 变截面直墙拱, 局部变形地基梁法, 压力拱理论中图分类号:U 451
当前, 地下空间已开始作为一种重要的自然资源加以开发, 在国民经济各部门都有大量的使用, 如城市地铁、穿越山岭的隧道、地下输油管道以及人防工程等。作为地下空间应用的一种主要结构形式) ) ) 直墙拱, 应用极为广泛。分析地下拱形结构的受力特点, 对地下工程的设计、施工等研究都有理论和工程应用意义。
文献标识码:A
1. 2 计算原理及其计算方法
在进行直墙拱结构内力计算时, 把拱圈和侧墙分开来计算, 拱圈按固定在墙顶上的无铰拱计算, 而侧墙按照竖放的弹性地基
梁计算, 同时考虑拱圈和侧墙的相互制约, 即侧墙墙顶的初始条件应与拱脚的内力和变位一致。
1 局部变形地基梁法1. 1 计算简图
由于隧道的长度远远大于隧道的断面尺寸, 因此可以将其简化成平面问题, 即沿纵向截取单位宽度拱带计算。拱圈和侧墙之间为刚性连接, 侧墙按双向弹性地基梁进行计算, 侧向地层抗力系数为K 1, 底部地层抗力系数为K 2。作用在直墙拱上的荷载主要有地层压力、结构自重及弹性抗力, 前面两者以梯形分布, 后者以抛物线分布, 具体计算简图见图1
。
1. 2. 1 拱圈的计算
由拱顶缺口相对变位为0, 有下列力法方程:
x 1D 11+x 2D 12+$2p +$1R +H 0=0
x 1D 21+x 2D 22+$2p +$2R +H 0f +u 0=0
(1) (2)
其中, D ik 为由x k =1作用使基本结构在拱顶处沿x i 作用方向产生的变位; $ip 为由主动外荷载作用使基本结构在拱顶处沿x i 作用方向产生的变位; $i R 为由被动弹性抗力作用使基本结构在拱顶处沿x i 作用方向产生的变位; H 0, u 0分别为拱脚的转角和水平位移。
由变位叠加原理, 有:
H 0=x 1H 1+x 2(H 2+f H 1) +H p u 0=x 1u 1+x 2(u 2+f u 1) +u p
(3) (4)
其中, H 1, u 1分别为墙顶作用单位力矩时引起的墙顶脚变位和水平位移; H 2, u 2分别为墙顶作用单位水平力时引起的墙顶脚
变位和水平位移; H p , u p 分别为外荷载(包括弹性抗力) 引起的墙顶脚变位和水平位移。
2) 协作体系斜拉桥主梁高度h 取为主跨跨径L 2的1/100~
1/150较为合适。梁高的降低, 大大减小了主梁的弯矩, 主梁变为以承受轴力为主。
3) 边跨长度L 1与主跨长度L 2的比值在1/1. 3~1/1. 5之间较为合理。研究表明, 合理取用协作体系的边跨长度, 其结构受力性能要优于边跨设置辅助墩的普通斜拉桥。当边跨取值接近上述合理范围的下限时, 将端锚索适当外移, 可明显改善结构的受力性能。
4) 合理的塔高, 既影响结构的内力及变形, 也关系到拉索的用钢量。分析表明, 主塔高度H 的合理取值在(1/2. 4~1/2. 8)
L 2左右。参考文献:
[1]林元培. 斜拉桥[M ].北京:人民交通出版社, 1995.
[2]陈德伟, 范立础, 张 权. 独塔斜拉桥的总体布置和参数研究[J].土木工程学报, 1999, 32(3) :34-40.
[3]戴利民, 石志源. 独塔协作体系斜拉桥设计参数分析[J]. 西安公路交通大学学报, 2000, 20(1) :34-38.
[4]戴利民. 独塔协作体系斜拉桥设计参数研究[J].结构工程师,
1999(3) :7-13.
S tu dies of d esign p arameter on cable -stayed brid ge of single -tower cooperative system
FENG Qi WANG Zhuan -feng
Abstract:For cable -stayed bridge of sing le -tow er cooper ative system often used in practice, this paper analyzes the reasonable form of this sys -tem in detail from the length w ithout cable, the ratios of dept h of beam to main span, the ratios o f side span to main span and tow er . s heig ht to main span, and reco mmends reasonable scope of par ameters, and provides r eference for the desig n of such br idge. Key words:cooperative system, cable -stayed bridge, the lengt h w ithout cable, the r atios of depth of beam to main span
收稿日期:2007-04-17
作者简介:张 林(1981-) , 男, 硕士, 助理工程师, 广东省建科建筑设计院, 广东广州 500100
#328#
第33卷第28期
2007年10月
山西建筑
根据变形协调条件, 即拱脚的变位取决于侧墙墙顶的变位, 因此H 1, u 1, H 2, u 2, H p , u p , 可以由局部变形地基梁法计算侧墙得到。
由于R h 为未知数, 由(4) 式所求得的x 1, x 2必须为R h 的函数, 由局部变形理论及叠加原理可得:
u 0=x 1u 1+x 2(u 2+f u 1) +u p =
R h k
(5)
根据该衬砌的受力特点, 参照有关设计规范, 综合考虑后, 选用直墙拱结构作为其支护结构。其净跨l 0=14m, 拱顶厚度d 0=0. 7m, 拱脚厚度d j =1. 0m, 内拱矢高f 0=3. 75m, 侧墙厚度d g =1. 0m, 内墙高h 0=2. 75m, 底板厚度d B =0. 12m, 墙基埋深d m =0. 2m 。
利用改进了的压力拱理论[1], 求得拱顶荷载为:
q =q 地+q 自=113. 5kPa; $q =$q 自=25. 5kPa 。
根据局部变形地基梁法, 可以求得拱身各截面的弯矩、内力见表1。
将(4) 式求得的x 1, x 2的表达式代入, 即可求得R h 。
1. 2. 2 侧墙的计算
侧墙按照双向弹性地基梁法进行计算:墙顶力:M h =X 1+X 2f +M hp +M h R #R j
H h =X 1+H hp +H h R #R j
(6)
(7)
3 数值模拟
文中运用AN SY S 大型有限元程序对该隧道变截面直墙拱衬砌内力进行了数值模拟, 先模拟隧道开挖的地层应力释放, 再分析隧道衬砌内力。
运用A NSYS 进行分析, 可得到各截面的弯矩、剪力, 与理论计算值比较如表2所示。
从表2可以看出, 根据局部变形地基梁法计算出的弯矩和轴力, 与大型AN SY S 有限元分析得出的数据最大偏差为37. 65%, 偏差还是比较显著的, 但是弯矩和轴力分布的变化趋势都相似, 在拱脚处弯矩较大, 侧墙底部的轴力较大。
参照文献[1]第8章弹性地基梁的有关公式求得侧墙任意界面内力。
表1 各截面的M, N 值
截面012345678
M 24. 51496123. 10587018. 86113011. 9567342. 839899-7. 785751-18. 97711-29. 59297-38. 33958
N 69. 59654070. 32441072. 22470375. 37225979. 67225984. 96264591. 01649097. 550816104. 23673
截面[**************]17
M -44. 52912-53. 82247-63. 614430-37. 132598-19. 002794-7. 808139-1. 9028240. 4576131. 046056
N 110. 79039117. 61179104. 969375106. 369723107. 770070109. 170415110. 570763111. 971110113. 371458
4 结语
1) 局部变形地基梁法计算出的数据与大型A NSYS 有限元分析得出的数据有一定的偏差, AN SYS 模拟得出的计算数据, 不能很好地反映隧道的真实受力情况, 但是弯矩和轴力分布的变化趋势与局部变形地基梁法计算结果比较相似, 在拱脚处弯矩较大, 侧墙底部的轴力较大。
2) 从局部变形地基梁法和大型AN SY S 有限元模拟都可以看出, 拱脚处的内力值最大, 即隧道衬砌的最危险点, 因此在设计施工过程中, 拱脚处的开挖、设计施工都要特别注意。
3) 使用局部地基梁法计算隧道衬砌内力时, 没有考虑侧墙脚的水平位移, 实际上侧墙将产生一定的水平位移, 这在计算过程中将导致一定误差, 影响计算结果的准确性。
4) 由于A NSYS 在岩土领域中运用的缺陷, 特别是本构关系中的不足, 导致模拟隧道开挖、隧道衬砌受力时存在偏差, 不很适用于岩土领域。文中采用的模拟方法, 没有考虑隧道开挖过程的影响, 仅适用于全断面开挖时, 隧道衬砌断面整体现浇的应力分布情况。参考文献:
[1]孙 钧, 侯学渊. 地下结构(上册) [M ].北京:科学出版社. 1987. [2]李皓月, 周田朋, 刘相新. ANSY S 工程计算应用教程[M ]. 北京:中国铁道出版社, 2003. [3]周勇超, 郝宪武. 双曲拱桥拓宽加固改造实例及力学性能分析[J].山西建筑, 2006, 32(1) :17-18.
表2 ANSYS 分析值与理论计算值的比较
项次截面
123拱圈
45678910
M(理论) 24. 51496123. 10587018. 86113011. 956734
M (ANSYS) 30. 9795627. 8125422. 9426813. 85307
%26. 3720. 3721. 6415. 86
N (理论) 69. 59654070. 32441072. 22470375. 37225979. 67225984. 96264591. 01649097. 550816104. 23673110. 79039117. 61179
N (ANSYS) 86. 0143685. 7395288. 2080388. 2081597. 96501102. 0826111. 7318122. 7092128. 2946137. 5684141. 8516
%23. 5921. 9222. 1317. 0322. 9620. 1522. 7625. 7923. 0824. 1720. 6122. 0621. 8520. 8319. 4620. 2718. 9517. 62
2. 8398993. 4947823. 06-7. 785751-9. 3140919. 63-18. 97711-23. 093221. 69-29. 59297-37. 22525. 79-38. 33958-48. 369226. 16-44. 52912-57. 135328. 31-53. 82247-64. 226419. 33
11-63. 614430-76. 776320. 69104. 969375128. 125612-37. 132598-50. 778837. 65106. 369723129. 6115侧墙
13-19. 002794-22. 662719. 26107. 770070130. 218614-7. 808139-8. 6170610. 36109. 170415130. 415151617
-1. 902824-2. 1633213. 69110. 570763132. 98350. 4576130. 54245418. 54111. 971110133. 18961. 046056
1. 172524
12. 09113. 371458133. 3475
2 实例计算
该工程为212国道重庆) 合川端中某一隧道工程, 围岩为软质石灰岩, 埋深H =50m, 岩石坚硬系数f i =2, 内摩擦角U =65b , 容重C =2400kg/m 3, 侧向岩层地基系数K 1=2. 5@105kPa/m, 基底岩层地基系数K 2=1. 25K 1=3. 125@10kP a/m 。
5
Analysis of underground arch structural inside force
ZHANG Lin
Abstract:According to partial deformed basement beam method, improved pressure arch theory are used in calculating the pressure of base -ment, inside force of a tunnel along national ro ad 212is also calculated, meantime, simulates by lar ge limited unit, analyzes the forcing fea -tures of arch structur e, has g uiding effect for designing and construction of tunnel pro ject.
Key words:under ground arch structur e, chang ed section dir ect arch, partial deformed basement beam method, pressure arch t heory
第33卷第28期西建筑 Vol. 33No. 28 山2007年10月Oct. 2007SHANXI ARCH ITECTURE
#327#
文章编号:1009-6825(2007) 28-0327-02
地下拱形结构内力分析
张 林
摘 要:根据局部变形地基梁法, 利用改进了的压力拱理论计算地层压力, 对212国道某一隧道衬砌(变截面直墙拱) 进行了内力计算, 同时利用大型有限元进行模拟, 分析了拱形结构的受力特点, 对隧道工程的设计与施工都具有一定的指导作用。
关键词:地下拱形结构, 变截面直墙拱, 局部变形地基梁法, 压力拱理论中图分类号:U 451
当前, 地下空间已开始作为一种重要的自然资源加以开发, 在国民经济各部门都有大量的使用, 如城市地铁、穿越山岭的隧道、地下输油管道以及人防工程等。作为地下空间应用的一种主要结构形式) ) ) 直墙拱, 应用极为广泛。分析地下拱形结构的受力特点, 对地下工程的设计、施工等研究都有理论和工程应用意义。
文献标识码:A
1. 2 计算原理及其计算方法
在进行直墙拱结构内力计算时, 把拱圈和侧墙分开来计算, 拱圈按固定在墙顶上的无铰拱计算, 而侧墙按照竖放的弹性地基
梁计算, 同时考虑拱圈和侧墙的相互制约, 即侧墙墙顶的初始条件应与拱脚的内力和变位一致。
1 局部变形地基梁法1. 1 计算简图
由于隧道的长度远远大于隧道的断面尺寸, 因此可以将其简化成平面问题, 即沿纵向截取单位宽度拱带计算。拱圈和侧墙之间为刚性连接, 侧墙按双向弹性地基梁进行计算, 侧向地层抗力系数为K 1, 底部地层抗力系数为K 2。作用在直墙拱上的荷载主要有地层压力、结构自重及弹性抗力, 前面两者以梯形分布, 后者以抛物线分布, 具体计算简图见图1
。
1. 2. 1 拱圈的计算
由拱顶缺口相对变位为0, 有下列力法方程:
x 1D 11+x 2D 12+$2p +$1R +H 0=0
x 1D 21+x 2D 22+$2p +$2R +H 0f +u 0=0
(1) (2)
其中, D ik 为由x k =1作用使基本结构在拱顶处沿x i 作用方向产生的变位; $ip 为由主动外荷载作用使基本结构在拱顶处沿x i 作用方向产生的变位; $i R 为由被动弹性抗力作用使基本结构在拱顶处沿x i 作用方向产生的变位; H 0, u 0分别为拱脚的转角和水平位移。
由变位叠加原理, 有:
H 0=x 1H 1+x 2(H 2+f H 1) +H p u 0=x 1u 1+x 2(u 2+f u 1) +u p
(3) (4)
其中, H 1, u 1分别为墙顶作用单位力矩时引起的墙顶脚变位和水平位移; H 2, u 2分别为墙顶作用单位水平力时引起的墙顶脚
变位和水平位移; H p , u p 分别为外荷载(包括弹性抗力) 引起的墙顶脚变位和水平位移。
2) 协作体系斜拉桥主梁高度h 取为主跨跨径L 2的1/100~
1/150较为合适。梁高的降低, 大大减小了主梁的弯矩, 主梁变为以承受轴力为主。
3) 边跨长度L 1与主跨长度L 2的比值在1/1. 3~1/1. 5之间较为合理。研究表明, 合理取用协作体系的边跨长度, 其结构受力性能要优于边跨设置辅助墩的普通斜拉桥。当边跨取值接近上述合理范围的下限时, 将端锚索适当外移, 可明显改善结构的受力性能。
4) 合理的塔高, 既影响结构的内力及变形, 也关系到拉索的用钢量。分析表明, 主塔高度H 的合理取值在(1/2. 4~1/2. 8)
L 2左右。参考文献:
[1]林元培. 斜拉桥[M ].北京:人民交通出版社, 1995.
[2]陈德伟, 范立础, 张 权. 独塔斜拉桥的总体布置和参数研究[J].土木工程学报, 1999, 32(3) :34-40.
[3]戴利民, 石志源. 独塔协作体系斜拉桥设计参数分析[J]. 西安公路交通大学学报, 2000, 20(1) :34-38.
[4]戴利民. 独塔协作体系斜拉桥设计参数研究[J].结构工程师,
1999(3) :7-13.
S tu dies of d esign p arameter on cable -stayed brid ge of single -tower cooperative system
FENG Qi WANG Zhuan -feng
Abstract:For cable -stayed bridge of sing le -tow er cooper ative system often used in practice, this paper analyzes the reasonable form of this sys -tem in detail from the length w ithout cable, the ratios of dept h of beam to main span, the ratios o f side span to main span and tow er . s heig ht to main span, and reco mmends reasonable scope of par ameters, and provides r eference for the desig n of such br idge. Key words:cooperative system, cable -stayed bridge, the lengt h w ithout cable, the r atios of depth of beam to main span
收稿日期:2007-04-17
作者简介:张 林(1981-) , 男, 硕士, 助理工程师, 广东省建科建筑设计院, 广东广州 500100
#328#
第33卷第28期
2007年10月
山西建筑
根据变形协调条件, 即拱脚的变位取决于侧墙墙顶的变位, 因此H 1, u 1, H 2, u 2, H p , u p , 可以由局部变形地基梁法计算侧墙得到。
由于R h 为未知数, 由(4) 式所求得的x 1, x 2必须为R h 的函数, 由局部变形理论及叠加原理可得:
u 0=x 1u 1+x 2(u 2+f u 1) +u p =
R h k
(5)
根据该衬砌的受力特点, 参照有关设计规范, 综合考虑后, 选用直墙拱结构作为其支护结构。其净跨l 0=14m, 拱顶厚度d 0=0. 7m, 拱脚厚度d j =1. 0m, 内拱矢高f 0=3. 75m, 侧墙厚度d g =1. 0m, 内墙高h 0=2. 75m, 底板厚度d B =0. 12m, 墙基埋深d m =0. 2m 。
利用改进了的压力拱理论[1], 求得拱顶荷载为:
q =q 地+q 自=113. 5kPa; $q =$q 自=25. 5kPa 。
根据局部变形地基梁法, 可以求得拱身各截面的弯矩、内力见表1。
将(4) 式求得的x 1, x 2的表达式代入, 即可求得R h 。
1. 2. 2 侧墙的计算
侧墙按照双向弹性地基梁法进行计算:墙顶力:M h =X 1+X 2f +M hp +M h R #R j
H h =X 1+H hp +H h R #R j
(6)
(7)
3 数值模拟
文中运用AN SY S 大型有限元程序对该隧道变截面直墙拱衬砌内力进行了数值模拟, 先模拟隧道开挖的地层应力释放, 再分析隧道衬砌内力。
运用A NSYS 进行分析, 可得到各截面的弯矩、剪力, 与理论计算值比较如表2所示。
从表2可以看出, 根据局部变形地基梁法计算出的弯矩和轴力, 与大型AN SY S 有限元分析得出的数据最大偏差为37. 65%, 偏差还是比较显著的, 但是弯矩和轴力分布的变化趋势都相似, 在拱脚处弯矩较大, 侧墙底部的轴力较大。
参照文献[1]第8章弹性地基梁的有关公式求得侧墙任意界面内力。
表1 各截面的M, N 值
截面012345678
M 24. 51496123. 10587018. 86113011. 9567342. 839899-7. 785751-18. 97711-29. 59297-38. 33958
N 69. 59654070. 32441072. 22470375. 37225979. 67225984. 96264591. 01649097. 550816104. 23673
截面[**************]17
M -44. 52912-53. 82247-63. 614430-37. 132598-19. 002794-7. 808139-1. 9028240. 4576131. 046056
N 110. 79039117. 61179104. 969375106. 369723107. 770070109. 170415110. 570763111. 971110113. 371458
4 结语
1) 局部变形地基梁法计算出的数据与大型A NSYS 有限元分析得出的数据有一定的偏差, AN SYS 模拟得出的计算数据, 不能很好地反映隧道的真实受力情况, 但是弯矩和轴力分布的变化趋势与局部变形地基梁法计算结果比较相似, 在拱脚处弯矩较大, 侧墙底部的轴力较大。
2) 从局部变形地基梁法和大型AN SY S 有限元模拟都可以看出, 拱脚处的内力值最大, 即隧道衬砌的最危险点, 因此在设计施工过程中, 拱脚处的开挖、设计施工都要特别注意。
3) 使用局部地基梁法计算隧道衬砌内力时, 没有考虑侧墙脚的水平位移, 实际上侧墙将产生一定的水平位移, 这在计算过程中将导致一定误差, 影响计算结果的准确性。
4) 由于A NSYS 在岩土领域中运用的缺陷, 特别是本构关系中的不足, 导致模拟隧道开挖、隧道衬砌受力时存在偏差, 不很适用于岩土领域。文中采用的模拟方法, 没有考虑隧道开挖过程的影响, 仅适用于全断面开挖时, 隧道衬砌断面整体现浇的应力分布情况。参考文献:
[1]孙 钧, 侯学渊. 地下结构(上册) [M ].北京:科学出版社. 1987. [2]李皓月, 周田朋, 刘相新. ANSY S 工程计算应用教程[M ]. 北京:中国铁道出版社, 2003. [3]周勇超, 郝宪武. 双曲拱桥拓宽加固改造实例及力学性能分析[J].山西建筑, 2006, 32(1) :17-18.
表2 ANSYS 分析值与理论计算值的比较
项次截面
123拱圈
45678910
M(理论) 24. 51496123. 10587018. 86113011. 956734
M (ANSYS) 30. 9795627. 8125422. 9426813. 85307
%26. 3720. 3721. 6415. 86
N (理论) 69. 59654070. 32441072. 22470375. 37225979. 67225984. 96264591. 01649097. 550816104. 23673110. 79039117. 61179
N (ANSYS) 86. 0143685. 7395288. 2080388. 2081597. 96501102. 0826111. 7318122. 7092128. 2946137. 5684141. 8516
%23. 5921. 9222. 1317. 0322. 9620. 1522. 7625. 7923. 0824. 1720. 6122. 0621. 8520. 8319. 4620. 2718. 9517. 62
2. 8398993. 4947823. 06-7. 785751-9. 3140919. 63-18. 97711-23. 093221. 69-29. 59297-37. 22525. 79-38. 33958-48. 369226. 16-44. 52912-57. 135328. 31-53. 82247-64. 226419. 33
11-63. 614430-76. 776320. 69104. 969375128. 125612-37. 132598-50. 778837. 65106. 369723129. 6115侧墙
13-19. 002794-22. 662719. 26107. 770070130. 218614-7. 808139-8. 6170610. 36109. 170415130. 415151617
-1. 902824-2. 1633213. 69110. 570763132. 98350. 4576130. 54245418. 54111. 971110133. 18961. 046056
1. 172524
12. 09113. 371458133. 3475
2 实例计算
该工程为212国道重庆) 合川端中某一隧道工程, 围岩为软质石灰岩, 埋深H =50m, 岩石坚硬系数f i =2, 内摩擦角U =65b , 容重C =2400kg/m 3, 侧向岩层地基系数K 1=2. 5@105kPa/m, 基底岩层地基系数K 2=1. 25K 1=3. 125@10kP a/m 。
5
Analysis of underground arch structural inside force
ZHANG Lin
Abstract:According to partial deformed basement beam method, improved pressure arch theory are used in calculating the pressure of base -ment, inside force of a tunnel along national ro ad 212is also calculated, meantime, simulates by lar ge limited unit, analyzes the forcing fea -tures of arch structur e, has g uiding effect for designing and construction of tunnel pro ject.
Key words:under ground arch structur e, chang ed section dir ect arch, partial deformed basement beam method, pressure arch t heory