盾构隧道开挖支护
我在模拟隧道开挖时，采用分步开挖，分步支护，开挖、支护第一步没有问题，等到开挖第二步时，网格显示有问题，支护也不完全，请问大家有没有遇到类似问题，大概是什么原因，谢谢！[ 本帖最后由 笑沧海 于 2007-4-30 10:18 编辑 ] lz ,你在吗? 加我QQ 83526695 刚刚试过没计算之前，开挖，支护都没有问题，难道计算一个荷载步以后，网格边界发生变形，不应该啊？ 取消了set large 可以算下去了；
但是发现结果不太理想，总是出现2cm 左右的隆起，而且跟第一次开挖的结果类似，不应该啊？感觉开挖了几个荷载步 跟第一步差不多的结果 ；
虚心学习，天天向上。：） 你结构单元加的可能有问题
你试一试在sel shell(或liner) id n 后面加上 group 'name'
name为你所加结构单元对象的组名
只是和你讨论不知道管不管用 你模拟的是不是上海的隧道工程，双圆盾构，
你要加了盾构面顶推力的话应该有隆起就对了。 原帖由 笑沧海 于 2007-4-18 23:03 发表 http://www.simwe.com/forum/images/common/back.gif
取消了set large 可以算下去了；
但是发现结果不太理想，总是出现2cm 左右的隆起，而且跟第一次开挖的结果类似，不应该啊？感觉开挖了几个荷载步 跟第一步差不多的结果 ；
118732
请教做过盾构隧道模拟的同行指点一下，问题大概出在什么地方，附上开挖命令流；
问题1)开挖分6个荷载步，发现后面得几个荷载步跟第一次开挖结果差不多，变化很小；2）隧道一直表现为隆起 ， 大概在3cm左右,这跟实际不符；3)用ansys算过该过程，虽然刚开始开挖隧道出现隆起，但后期随着盾构推进表现为沉降；不知道问题大概出在哪里，请教大家！论文中......
new
rest gravstress.sav
plot block group
ini xdis=0 ydis=0 zdis=0
;材料参赋值,将计算参数恢复到原始值
model mohr
prop bulk b_mod1 shear s_mod1 fric 19.33 coh 18.56e3 ten 1e20 dil 0 &
range group 1
prop bulk b_mod2 shear s_mod2 fric 35.0 coh 10.0e3 ten 1e20 dil 0 &
range group 2
prop bulk b_mod3 shear s_mod3 fric 16.0 coh 13.0e3 ten 1e20 dil 0 &
range group 3
prop bulk b_mod4 shear s_mod4 fric 11.0 coh 14.0e3 ten 1e20 dil 0 &
range group 4
prop bulk b_mod5 shear s_mod5 fric 14.0 coh 16.0e3 ten 1e20 dil 0 &
range group 5
prop bulk b_mod6 shear s_mod6 fric 24.46 coh 19.49e3 ten 1e20 dil 0 &
range group 6
;注浆体，开挖土体n8,n9按n4土层参数赋值
prop bulk b_mod4 shear s_mod4 fric 11.0 coh 14.0e3 ten 1e20 dil 0 &
range group 7
prop bulk b_mod4 shear s_mod4 fric 11.0 coh 14.0e3 ten 1e20 dil 0 &
range group 8
prop bulk b_mod4 shear s_mod4 fric 11.0 coh 14.0e3 ten 1e20 dil 0 &
range group 9
def derive2 ;定义衬砌,注浆体材料参数
b_mod8=e_mod8/(3.0*(1.0-2.0*p_ratio8))
s_mod8=e_mod8/(2.0*(1.0+p_ratio8))
b_mod9=e_mod9/(3.0*(1.0-2.0*p_ratio9))
s_mod9=e_mod9/(2.0*(1.0+p_ratio9))
end
set e_mod8=1.2e7 p_ratio8=0.2 e_mod9=2.576e9 p_ratio9=0.286
derive2
;隧道分六步开挖
;右侧隧道开挖,设置shell单元模拟衬砌（盾构推进到1～4环）
model null range cylind end1 2.3 0 48 end2 2.3 0 43.1 radi 3.15 x -0.1 5.46
sel shell id=1 range cylind end1 2.3 0 48 end2 2.3 0 43.1 radius 3.15 x -0.1 5.46 &
z 43.2 not
plot add sel geom node=off shrinkfac=0 magfac=0 red red
sel shell prop iso=(3.45e10,0.167) thick 0.3 density 2500
sel node fix z xr yr range z 47.9 48.1
pause
;左侧隧道开挖,生成衬砌
model null range cylind end1 -2.3 0 48 end2 -2.3 0 43.1 radi 3.15 x -5.46 0.1
sel shell id=1 range cylind end1 -2.3 0 48 end2 -2.3 0 43.1 radi 3.15 x -5.46 0.1 &
z 43.2 not
plot add sel geom node=off shrinkfac=0 magfac=0 red red
sel shell prop iso=(3.45e10,0.167) thick 0.3 density 2500
sel node fix z xr yr range z 47.9 48.1
pause
;施加正面推力
apply nstress -0.21e6 range z 43.2 x -0.1 5.46 y -3.16 3.16
apply nstress -0.21e6 range z 43.2 x -5.46 0.1 y -3.16 3.16
model elas range z 43.1 48 group 7;生成初凝注浆体
prop density 1650 bulk b_mod8 shear s_mod8
;set large
hist unbal
plot cont ydisp
plot add hist 1
solve
save excav1.sav ;第一个开挖步
ret
以下过程类似； 原帖由 笑沧海 于 2007-4-23 20:37 发表 http://www.simwe.com/forum/images/common/back.gif
请教做过盾构隧道模拟的同行指点一下，问题大概出在什么地方，附上开挖命令流；
问题1)开挖分6个荷载步，发现后面得几个荷载步跟第一次开挖结果差不多，变化很小；2）隧道一直表现为隆起 ， 大概在3cm左 ...
第一步结果ydisp:
第六步结果ydisp:
你好，能不能将建模的命令流也贴出来啊，我也正在做盾构的模拟，在建模时有点问题，不知道管片是怎么模拟的．谢谢！ 正面推力是如何计算的啊？ 看看短消息吧！ 谢谢回复，前几天调整了一些参数，结果仍不理想；
;右侧隧道开挖,设置shell单元模拟衬砌（盾构推进到1～4环）
model null range cylind end1 2.3 0 48 end2 2.3 0 43.1 radi 3.15 x -0.1 5.46
sel shell id=1 range cylind end1 2.3 0 48 end2 2.3 0 43.1 radius 3.15 x -0.1 5.46 &
z 43.2 not ;生成右侧衬砌及中间立柱
plot add sel geom node=off shrinkfac=0 magfac=0 red red
sel shell prop iso=(3.45e10,0.167) thick 0.3 density 2500
sel node fix z xr yr range z 47.9 48.1
pause
;左侧隧道开挖,生成衬砌
model null range cylind end1 -2.3 0 48 end2 -2.3 0 43.1 radi 3.15 x -5.46 0.1
sel shell id=1 range cylind end1 -2.3 0 48 end2 -2.3 0 43.1 radi 3.15 x -5.46 0.1 &
z 43.2 not ;生成左侧衬砌
plot add sel geom node=off shrinkfac=0 magfac=0 red red
sel shell prop iso=(3.45e10,0.167) thick 0.3 density 2500
sel node fix z xr yr range z 47.9 48.1
pause
通过上面语句，可以生成双圆的衬砌，双圆中间是有立柱的，不太明白楼主的意思，请指点，谢谢！ 原帖由 笑沧海 于 2007-5-10 10:09 发表 http://www.simwe.com/forum/images/common/back.gif
谢谢回复，前几天调整了一些参数，结果仍不理想；
;右侧隧道开挖,设置shell单元模拟衬砌（盾构推进到1～4环）
model null range cylind end1 2.3 0 48 end2 2.3 0 43.1 radi 3.15 x -0.1 5.46
sel shell...
共同学习中,我的理解上肯定有不正确之处，永远谈不上指点。
我说的是：
1、为什么不将左侧隧道的model null range cylind end1 2.3 0 48 end2 2.3 0 43.1 radi 3.15 x-5.46 0.1 中的x= -5.46 0.1
写成x=-5.46 -0.1，听大师说结构单元会自动进行连接的。shell的厚度也给出了，我认为可以模拟中柱的效果。
2、隧道周围的土体的抗拉强度为e20，为什么不调回来。
3、若采用liner单元可通过调整cs_nk，cs_sk等参数得出与shell单元差别很大的结果。 这几天在赶论文，暂未调试；
讨论：1.右侧开挖后，“x= -5.46 0.1 写成x=-5.46 -0.1”效果是相同的，我认为，因为shell只能生成在结构的临空面上；
2.抗拉强度我是人为设置很大，本意不考虑抗拉破坏，不知道会不会有影响；
3.liner单元没有用过，感觉盾构管片可以用shell来模拟，liner用于模拟结构之间的滑移、脱离，shell是共结点的接触；÷
4.看了几个人的模拟，他们好像是一步开挖，一步支护的，不知道是否有影响；
谢谢！ 请问各位老大，谁算过软土浅埋隧道啊，为什么我调了这么久参数结果还是隆起，几乎没信心了?该如何调整才能得到与实际相符合的结果啊，困惑！
回复 #16 笑沧海 的帖子
问题解决了吗？以我的经验来看，FLAC和ansys模拟软土隧道效果都不好，调大岩土参数再试试
回复 #16 笑沧海 的帖子
set e_mod8=1.2e7 p_ratio8=0.2 e_mod9=2.576e9 p_ratio9=0.286个人认为是不是衬砌和注浆体的参数太小了呀。