第10卷第2期水利与建筑工程学报
Vol. 10No. 2
隧道荷载结构法在有限元软件中的实现
王力达
(中铁上海设计院集团有限公司, 上海200071)
摘 要:荷载-结构法是盾构隧道设计中的常用方法, 但其作用力形式在Marc 中无法直接施加, 且地层弹簧的具体布置位置是未知的, 这就限制了该法在Marc 软件中的应用。通过力的转换, 采用Fortran 语言编写子程序施加该力; 采用迭代算法, 通过每次计算后, 令受拉弹簧退出计算, 并与上次计算的地层弹簧位置进行比较, 直到相邻两次计算的地层弹簧位置完全相同时, 结束计算。从而实现了用Marc 进行隧道荷载结构法的计算, 扩大了Marc 软件对地下工程计算分析的应用范围。并通过检验算例说明了开发思路和接口子程序的正确性。
关键词:盾构隧道; 荷载结构法; 有限元法; 二次开发中图分类号:U45
文献标识码:A
文章编号:1672) 1144(2012) 02) 0105) 03
Implementation of Shield Tunnel . s Load S tructure Method in FEM Software
W ANG L-i da
(China Railway Shanghai Design Institute Grou p Co. , Ltd. , Shan ghai 200071, China)
Abstract:Load structure method is commonly used in tunnel design, but the force cannot be applied directly in Marc, and the position of the spring around the tunnel is unkno wn. By means of force c onversion, a subroutine programming with Fortran is adopted to exert the force. An iterative algorithm is also adopted, the tension spring will exit from the ca-l culation, and the comparison for the spring positions in the two steps will be made. If the positions are the sa me, the computation will be finished. The work above would realize the application of the load structure method in Marc software and enlarge the application of Marc software in the calculation of underground engineering. At the same time, the verif-i cation exa mple is given to prove the validity of the interface subroutine.
Keywords:shield tunnel; load structure method; finite element method; secondary development
0 引 言
目前盾构法隧道设计的一个重要环节是采用合适的方法计算衬砌结构的内力
[1-2]
Marc 进行隧道荷载结构法的计算, 扩大了Marc 软件对地下工程计算分析的应用范围, 最后通过检验算例说明了开发思路和接口子程序的正确性。
。地下工程中常
用的计算方法是荷载结构法和地层结构法[3], 荷载结构法仍然是目前隧道结构理论计算用得最多的一种方法, 我国地下铁道及铁路隧道设计规范中均推荐使用此方法[4]。
本文通过力的转换, 采用Fortran 语言编写子程序施加作用于衬砌上的荷载。并采用迭代算法, 通过每次计算后, 令受拉弹簧退出计算, 并与上次计算的地层弹簧位置进行比较, 直到相邻两次计算的地层弹簧位置完全相同时, 结束计算。从而实现了用
1 计算模型及Marc 软件二次开发流程
1. 1 计算模型
图1是衬砌外围的水土压力分布图, 图中:P vs 1为衬砌顶部垂直土压力; P vw 1为衬砌顶部垂直水压力; P v 2为衬砌底部均布土反力; P g 为管片自重; P hs 1为衬砌顶部水平土压力; P hw 1为衬砌顶部水平水压力; P hs 2为衬砌底部水平土压力; P hw 2为衬砌底部水平水压力。
收稿日期:2011-12-20 修稿日期:2012-01-15
作者简介:王力达(1963) ) , 男(汉族) , 甘肃兰州人, 高级工程师, 主要从事岩土工程与地下工程的研究与设计工作。
106 水利与建筑工程学报 第10卷
图中, V 为荷载结构法中所施加的竖向荷载, V 1
为转化为梁单元上的荷载。根据力的等效原理, 则:
V 1@d s =V @d x
根据图2上的角度关系, 上式可以写为:
V 1=V @sin (B )
化为梁单元上的荷载。
h 1=h @cos (B )
化为梁上的荷载。
利用Marc 中提供的子程序施加竖向及水平荷
图1 衬砌外围的水土压力分布
(1) (2)
利用同样的方法可以将荷载结构法中的水平荷载转
(3)
式中:h 为荷载结构法中所施加的水平荷载, h 1为转
载, 其主要思路为:根据边界条件判断是施加水平荷载还是竖向荷载, 如果是施加竖向荷载, 首先根据积分点的坐标计算得到角度B , 然后根据式(2) 进行施加; 如果是施加水平荷载, 首先根据积分点竖向坐标计算得到该点的水平荷载值, 然后根据式(3) 进行施加。主要流程如图4所示。
荷载-结构模型中, 需要重点考虑和解决的问题就是结构与地层的相互作用问题, 即地层抗力如何考虑的问题。当前主要有以下几种方式来考虑地层抗力:(a) 忽略土层的被动抗力; (b) 考虑被动抗力, 被动抗力按假定分布; (c) 用地层弹簧来考虑地基抗力的作用[5]。本文选用地层弹簧来考虑地基抗力的作用。
1. 2 Marc 软件二次开发流程
由图1可以看出, 力的形式无法直接施加到梁单元上, 因此需要做必要的转化。取出一微元段进行分析, 该微元段位于图2中的A 点。A 点的受力示意图见图3
。
图4 子程序施加荷载流程图
利用Marc 程序进行荷载结构法计算的主要思路为:首先在整个衬砌外围施加地层弹簧, 利用上述的方法施加隧道衬砌荷载; 根据计算得到的地层弹簧的受力状态, 确定新的地层弹簧施加范围, 当地层弹簧的受力状态为拉时, 该弹簧退出计算, 为压时保
图2
微元段位置示意图
留; 进行迭代计算, 直到前后两次算的地层弹簧的范围一致时, 结束计算。主要流程如图5所示。
第2期 王力达:隧道荷载结构法在有限元软件中的实现107
2 验证算例
利用上节中的思路对Marc 进行二次开发, 从而实现了用Marc 进行隧道荷载结构法的计算。以下通过与文献[6]中11. 2节中的设计实例比较, 来说明程序开发的有效性。
该设计实例为一圆形隧道, 采用最普遍的常用设计法, 亦即不考虑接头影响的等弯曲刚度均匀圆环法。其中管片外径D =3. 35m, 管片厚为0. 125m 。
管片上力的分布形式如图1所示, 其具体数值见表1。
表1 管片上的荷载值
P vs 153600
P v w 1
P v 2
P g 10210
P hs 120750
P hw 1130630
P hs 239950
本次计算, 管片采用连续框架梁模拟, 周围施加地基弹簧用以模拟地层与结构的相互作用。管片一共划分为720个节点, 720个梁单元。有限元模型见图6
。
单位:Pa
P hw 2162880
[1**********]0
图6 有限元模型
表1中荷载符号的含义在图1中已有说明, 其具体计算可参看文献[6]。地基反力系数为2@107Pa 。
将本文的计算结果与文献[6]中的结果进行比较, 具体见表2。
表2 计算结果比较
角度/b 060120180
弯矩值/Pa
M arc 解6410. 460-4396. 720-274. 5132203. 430
文献[6]中解
6520-4290-2602440
弯矩相对误差
1. 68%2. 49%5. 58%9. 70%
轴力值/N
Marc 解[***********]308526
文献[6]中解[***********]307050
轴力相对误差
0. 83%0. 39%0. 39%0. 48%
通过上表可以看出, 两者的弯矩与轴力的变化趋势是完全一样的, 在数值方面略有差别。弯矩的最大误差为9. 70%, 轴力的最大误差为0. 83%。两者在数值上的差别主要是由于Marc 采用了地层弹簧形式, 而文献[6]中采用的是三角形抗力假定。文献[7]比较了采用各种模型时的弯矩及轴力值, 并得出在假定三角形抗力时, 在竖向荷载和水平土压力相同的情况下, 地基抗力偏小, 从而造成拱顶正弯矩偏大。通过表2可以看出, 该结论与本文得出的结果是一致的。
通过对比可以看出, 两者的弯矩与轴力的变化趋势以及数值基本上是吻合的, 说明了本文对Marc 的二次开发是有效的。
出, 两者的弯矩与轴力的变化趋势是完全一样的, 在数值方面也是基本吻合的, 说明本文对Marc 的二次开发是有效的。
通过二次开发扩大了Marc 软件对地下工程计算分析的应用范围。参考文献:
[1] 刘 琼, 吴雄志, 姚 捷, 等. 盾构隧道管片衬砌内力计
算方法对比分析[J]. 河北工程大学学报(自然科学版) , 2008, 25(3) :26-29.
[2] 关永平, 郭 龙, 李云龙, 等. 城市地铁开挖对相邻地下
管线影响的数值分析[J].水利与建筑工程学报, 2010, 8(2) :11-12.
[3] 孙 钧, 侯学渊. 地下结构[M ]. 北京:科学出版社, 1987. [4] 武振国, 常 翔, 叶 飞. 盾构隧道结构设计模型综述[J].隧道建设, 2008, 28(2) :182-185.
[5] 丁军霞, 冯卫星, 张保俭. 盾构隧道管片衬砌内力计算
[J].石家庄铁道学院学报, 2004, 17(2) :66-69.
[6] 张凤祥, 朱合华, 傅德明. 盾构隧道[M]. 北京:人民交通出版社, 2004. [7] 孙文昊. 荷载-结构模型在盾构管片内力计算中应用
[, -3 结 语
荷载结构法是目前进行隧道结构理论计算用得
最多的一种方法, 通过对Marc 进行二次开发, 采用Fortran 语言编写子程序施加隧道荷载, 并利用迭代算法确定地层弹簧的位置, 从而实现了利用Marc 进行荷载结构法的计算。
第10卷第2期水利与建筑工程学报
Vol. 10No. 2
隧道荷载结构法在有限元软件中的实现
王力达
(中铁上海设计院集团有限公司, 上海200071)
摘 要:荷载-结构法是盾构隧道设计中的常用方法, 但其作用力形式在Marc 中无法直接施加, 且地层弹簧的具体布置位置是未知的, 这就限制了该法在Marc 软件中的应用。通过力的转换, 采用Fortran 语言编写子程序施加该力; 采用迭代算法, 通过每次计算后, 令受拉弹簧退出计算, 并与上次计算的地层弹簧位置进行比较, 直到相邻两次计算的地层弹簧位置完全相同时, 结束计算。从而实现了用Marc 进行隧道荷载结构法的计算, 扩大了Marc 软件对地下工程计算分析的应用范围。并通过检验算例说明了开发思路和接口子程序的正确性。
关键词:盾构隧道; 荷载结构法; 有限元法; 二次开发中图分类号:U45
文献标识码:A
文章编号:1672) 1144(2012) 02) 0105) 03
Implementation of Shield Tunnel . s Load S tructure Method in FEM Software
W ANG L-i da
(China Railway Shanghai Design Institute Grou p Co. , Ltd. , Shan ghai 200071, China)
Abstract:Load structure method is commonly used in tunnel design, but the force cannot be applied directly in Marc, and the position of the spring around the tunnel is unkno wn. By means of force c onversion, a subroutine programming with Fortran is adopted to exert the force. An iterative algorithm is also adopted, the tension spring will exit from the ca-l culation, and the comparison for the spring positions in the two steps will be made. If the positions are the sa me, the computation will be finished. The work above would realize the application of the load structure method in Marc software and enlarge the application of Marc software in the calculation of underground engineering. At the same time, the verif-i cation exa mple is given to prove the validity of the interface subroutine.
Keywords:shield tunnel; load structure method; finite element method; secondary development
0 引 言
目前盾构法隧道设计的一个重要环节是采用合适的方法计算衬砌结构的内力
[1-2]
Marc 进行隧道荷载结构法的计算, 扩大了Marc 软件对地下工程计算分析的应用范围, 最后通过检验算例说明了开发思路和接口子程序的正确性。
。地下工程中常
用的计算方法是荷载结构法和地层结构法[3], 荷载结构法仍然是目前隧道结构理论计算用得最多的一种方法, 我国地下铁道及铁路隧道设计规范中均推荐使用此方法[4]。
本文通过力的转换, 采用Fortran 语言编写子程序施加作用于衬砌上的荷载。并采用迭代算法, 通过每次计算后, 令受拉弹簧退出计算, 并与上次计算的地层弹簧位置进行比较, 直到相邻两次计算的地层弹簧位置完全相同时, 结束计算。从而实现了用
1 计算模型及Marc 软件二次开发流程
1. 1 计算模型
图1是衬砌外围的水土压力分布图, 图中:P vs 1为衬砌顶部垂直土压力; P vw 1为衬砌顶部垂直水压力; P v 2为衬砌底部均布土反力; P g 为管片自重; P hs 1为衬砌顶部水平土压力; P hw 1为衬砌顶部水平水压力; P hs 2为衬砌底部水平土压力; P hw 2为衬砌底部水平水压力。
收稿日期:2011-12-20 修稿日期:2012-01-15
作者简介:王力达(1963) ) , 男(汉族) , 甘肃兰州人, 高级工程师, 主要从事岩土工程与地下工程的研究与设计工作。
106 水利与建筑工程学报 第10卷
图中, V 为荷载结构法中所施加的竖向荷载, V 1
为转化为梁单元上的荷载。根据力的等效原理, 则:
V 1@d s =V @d x
根据图2上的角度关系, 上式可以写为:
V 1=V @sin (B )
化为梁单元上的荷载。
h 1=h @cos (B )
化为梁上的荷载。
利用Marc 中提供的子程序施加竖向及水平荷
图1 衬砌外围的水土压力分布
(1) (2)
利用同样的方法可以将荷载结构法中的水平荷载转
(3)
式中:h 为荷载结构法中所施加的水平荷载, h 1为转
载, 其主要思路为:根据边界条件判断是施加水平荷载还是竖向荷载, 如果是施加竖向荷载, 首先根据积分点的坐标计算得到角度B , 然后根据式(2) 进行施加; 如果是施加水平荷载, 首先根据积分点竖向坐标计算得到该点的水平荷载值, 然后根据式(3) 进行施加。主要流程如图4所示。
荷载-结构模型中, 需要重点考虑和解决的问题就是结构与地层的相互作用问题, 即地层抗力如何考虑的问题。当前主要有以下几种方式来考虑地层抗力:(a) 忽略土层的被动抗力; (b) 考虑被动抗力, 被动抗力按假定分布; (c) 用地层弹簧来考虑地基抗力的作用[5]。本文选用地层弹簧来考虑地基抗力的作用。
1. 2 Marc 软件二次开发流程
由图1可以看出, 力的形式无法直接施加到梁单元上, 因此需要做必要的转化。取出一微元段进行分析, 该微元段位于图2中的A 点。A 点的受力示意图见图3
。
图4 子程序施加荷载流程图
利用Marc 程序进行荷载结构法计算的主要思路为:首先在整个衬砌外围施加地层弹簧, 利用上述的方法施加隧道衬砌荷载; 根据计算得到的地层弹簧的受力状态, 确定新的地层弹簧施加范围, 当地层弹簧的受力状态为拉时, 该弹簧退出计算, 为压时保
图2
微元段位置示意图
留; 进行迭代计算, 直到前后两次算的地层弹簧的范围一致时, 结束计算。主要流程如图5所示。
第2期 王力达:隧道荷载结构法在有限元软件中的实现107
2 验证算例
利用上节中的思路对Marc 进行二次开发, 从而实现了用Marc 进行隧道荷载结构法的计算。以下通过与文献[6]中11. 2节中的设计实例比较, 来说明程序开发的有效性。
该设计实例为一圆形隧道, 采用最普遍的常用设计法, 亦即不考虑接头影响的等弯曲刚度均匀圆环法。其中管片外径D =3. 35m, 管片厚为0. 125m 。
管片上力的分布形式如图1所示, 其具体数值见表1。
表1 管片上的荷载值
P vs 153600
P v w 1
P v 2
P g 10210
P hs 120750
P hw 1130630
P hs 239950
本次计算, 管片采用连续框架梁模拟, 周围施加地基弹簧用以模拟地层与结构的相互作用。管片一共划分为720个节点, 720个梁单元。有限元模型见图6
。
单位:Pa
P hw 2162880
[1**********]0
图6 有限元模型
表1中荷载符号的含义在图1中已有说明, 其具体计算可参看文献[6]。地基反力系数为2@107Pa 。
将本文的计算结果与文献[6]中的结果进行比较, 具体见表2。
表2 计算结果比较
角度/b 060120180
弯矩值/Pa
M arc 解6410. 460-4396. 720-274. 5132203. 430
文献[6]中解
6520-4290-2602440
弯矩相对误差
1. 68%2. 49%5. 58%9. 70%
轴力值/N
Marc 解[***********]308526
文献[6]中解[***********]307050
轴力相对误差
0. 83%0. 39%0. 39%0. 48%
通过上表可以看出, 两者的弯矩与轴力的变化趋势是完全一样的, 在数值方面略有差别。弯矩的最大误差为9. 70%, 轴力的最大误差为0. 83%。两者在数值上的差别主要是由于Marc 采用了地层弹簧形式, 而文献[6]中采用的是三角形抗力假定。文献[7]比较了采用各种模型时的弯矩及轴力值, 并得出在假定三角形抗力时, 在竖向荷载和水平土压力相同的情况下, 地基抗力偏小, 从而造成拱顶正弯矩偏大。通过表2可以看出, 该结论与本文得出的结果是一致的。
通过对比可以看出, 两者的弯矩与轴力的变化趋势以及数值基本上是吻合的, 说明了本文对Marc 的二次开发是有效的。
出, 两者的弯矩与轴力的变化趋势是完全一样的, 在数值方面也是基本吻合的, 说明本文对Marc 的二次开发是有效的。
通过二次开发扩大了Marc 软件对地下工程计算分析的应用范围。参考文献:
[1] 刘 琼, 吴雄志, 姚 捷, 等. 盾构隧道管片衬砌内力计
算方法对比分析[J]. 河北工程大学学报(自然科学版) , 2008, 25(3) :26-29.
[2] 关永平, 郭 龙, 李云龙, 等. 城市地铁开挖对相邻地下
管线影响的数值分析[J].水利与建筑工程学报, 2010, 8(2) :11-12.
[3] 孙 钧, 侯学渊. 地下结构[M ]. 北京:科学出版社, 1987. [4] 武振国, 常 翔, 叶 飞. 盾构隧道结构设计模型综述[J].隧道建设, 2008, 28(2) :182-185.
[5] 丁军霞, 冯卫星, 张保俭. 盾构隧道管片衬砌内力计算
[J].石家庄铁道学院学报, 2004, 17(2) :66-69.
[6] 张凤祥, 朱合华, 傅德明. 盾构隧道[M]. 北京:人民交通出版社, 2004. [7] 孙文昊. 荷载-结构模型在盾构管片内力计算中应用
[, -3 结 语
荷载结构法是目前进行隧道结构理论计算用得
最多的一种方法, 通过对Marc 进行二次开发, 采用Fortran 语言编写子程序施加隧道荷载, 并利用迭代算法确定地层弹簧的位置, 从而实现了利用Marc 进行荷载结构法的计算。