我做的双向隧道开挖模拟
这是竖直方向的位移图(即沉降),沉降规律根本就不对,请大家指点啊。小弟实在找不出原因了。下面将代码贴出。 怎么没人回答啊?那位大侠帮指点以下啊 new
;建立隧道区间模型,试运行
plot grid
plot add sur green
gen zone radcyl p0 5.58 0 22.43 p1 5.58 0 28.01 p2 5.58 100 22.43 p3 0 0 22.43 size 1 20 4 2 dim 3.5 3.5 3.5 3.5
gen zone cshell p0 5.58 0 22.43 p1 5.58 0 25.93 p2 5.58 100 22.43 p3 2.08 0 22.43 size 1 20 4 1 dim 3 3 3 3 group concrete fill group tunnel
gen zone reflect normal 1 0 0 origin 5.58 0 22.43
gen zone reflect normal 0 0 -1 origin 5.58 0 22.43
gen zone brick p0 11.16 0 16.85 p1 24.28 0 16.85 p2 11.16 100 16.85 p3 11.16 0 28.01 size 6 20 4
gen zone brick p0 24.28 0 16.85 p1 60.28 0 16.85 p2 24.28 100 16.85 p3 24.28 0 28.01 size 12 20 4
gen zone brick p0 0 0 0 p1 11.16 0 0 p2 0 100 0 p3 0 0 16.85 size 4 20 6
gen zone brick p0 11.16 0 0 p1 24.28 0 0 p2 11.16 100 0 p3 11.16 0 16.85 size 6 20 6
gen zone brick p0 24.28 0 0 p1 60.28 0 0 p2 24.28 100 0 p3 24.28 0 16.85 size 12 20 6
gen zone brick p0 0 0 28.01 p1 11.16 0 28.01 p2 0 100 28.01 p3 0 0 40 size 4 20 3
gen zone brick p0 11.16 0 28.01 p1 24.28 0 28.01 p2 11.16 100 28.01 p3 11.16 0 40 size 6 20 3
gen zone brick p0 24.28 0 28.01 p1 60.28 0 28.01 p2 24.28 100 28.01 p3 24.28 0 40 size 12 20 3
gen zone reflect normal -1 0 0 origin 0 0 0
save model.sav
;先用弹性计算
m e
prop b 6.81e6 sh 3.15e6range z 0 16.85
prop b 2.27e7 sh 1.297e7 range z 16.85 28.01
prop b 6.19e6 sh 2.86e6range z 28.01 40
ini dens 1790 range z 0 16.85
ini dens 2050 range z 16.85 28.01
ini dens 1980 range z 28.01 40
set gravity 0 0 -9.8
;施加边界条件
fix x range x -60.38 -60.18
fix x range x 60.18 60.38
fix y range y -0.1 0.1
fix y range y 99.9 100.1
fix x y z range z -0.1 0.1
plot con zdis
solve
save gravity.sav
;施加楼房荷载,由于在隧道附近有一栋楼房,故将视为均布荷载
ini xdisp 0 ydis 0 zdis 0
m m
prop b 6.81e6 sh 3.15e6coh 5.4e4 fric 9ten 1e10 range z 0 16.85
prop b 2.27e7 sh 1.297e7 coh 5e4 fric 35 ten 1e10 range z 16.85 28.01
prop b 6.19e6 sh 2.86e6coh 3.5e4 fric 15 ten 1e10 range z 28.01 40
apply nstress -13.2242e4 range x -24.28 -11.16 y 40 60 z 39.9 40.1 ;不知道这样简化加楼房荷载可以不,请高手指点
set large
solve
save structure.sav
;开挖隧道
ini xdis 0 ydis 0 zdis 0
model null range group tunnel x -11.16 0
model null range group concrete x -11.16 0
step 50
model elas range group concrete x -11.16 0
prop bulk 0.6e9 shear 0.36e9 dens 2500 range group concrete x -11.16 0
solve
save left.sav
model null range group tunnel x 0 11.16
model null range group concrete x 0 11.16
step 50
model elas range group concrete x 0 11.16
prop bulk 0.6e9 shear 0.36e9 dens 2500 range group concrete x 0 11.16
solve
save right.sav 楼主,你的抗拉强度怎么设得那么大,1e10,是故意设得这么大的吗?
没有运行你的命令,apply nstress应该是可以的,不知道可否考虑用apply szz呢? 抗拉强度确实是故意设的很大,但是位移云图明显不对,自己也不知道是什么原因。
回复 #1 xiangruide 的帖子
简单看了一下你的命令流,是不是开挖后应力释放的有点多了?可以将step 50改一下,只是个人意见,不一定对噢 原帖由 xiangruide 于 2007-4-6 14:34 发表 http://www.simwe.com/forum/images/common/back.gifnew
;建立隧道区间模型,试运行
plot grid
plot add sur green
gen zone radcyl p0 5.58 0 22.43 p1 5.58 0 28.01 p2 5.58 100 22.43 p3 0 0 22.43 size 1 20 4 2 dim 3.5 3.5 3.5 3.5
gen zone cshell p ...
大概看了下命令,好像隧道得开挖方向是y 而不是x 啊,请检查一下; ! C:\Documents and Settings\Administrator\桌面\浙江舟山庆丰采石场高边坡稳定性分析及治理研.part3.rar: 必需的压缩卷不存在
! C:\Documents and Settings\Administrator\桌面\浙江舟山庆丰采石场高边坡稳定性分析及治理研.part2.rar: 浙江舟山庆丰采石场高边坡稳定性分析及治理研.nh CRC 失败。文件被破坏 还是有很多热心人,谢谢啊 还是有很多热心人,谢谢啊 好资料 谢谢分享:) henhao ,以后欧普多交流 很好,谢谢分享 楼主,请问您的问题是怎么解决的,求赐教 好好学习,多谢楼主分享 xiangruide 发表于 2007-4-6 14:34
new
;建立隧道区间模型,试运行
plot grid
问下楼主,你这个步数控制是什么原理啊? 开挖隧道时位移又重置为0做什么。 xuexixuexi
页:
[1]