模拟隧道开挖基底反弹很大

hjm460727715

我模拟了一个隧道开挖 在第一循环上台阶的时候 基底反弹竟然有0.7m，拱顶下沉 还算正常，应该不是应力释放的问题。盼各位高手 指点下～

hjm460727715

附上 zdis云图

hjm460727715

拱腰的 净空收敛也 很怪异 土体竟然 是想隧道外移动

hjm460727715

水平位移云图

cloud328

应当有问题

ycding

模型有问题吧？
再仔细检查检查。

hjm460727715

1. n
   
2. plo show
   
3. plo blo group
   
4. conf dyn fluid;设置动态与流体算法
   
5. set dyn off fluid off;关闭动态与流体算法
   
6. 
   
7. 
   
8. ;================上部加固圈====================
   
9. ge zo cs p0 0 0 0 p1 11.45 0 0 p2 0 50 0 p3 0 0 11.45 &
   
10. dim 9.45 9.45 9.45 9.45 size 4 50 15 10 group '加固圈'
    
11. ;================下部加固圈====================
    
12. ge zo cs p0 0 50 0 p1 11.45 50 0 p2 0 0 0 p3 0 50 -11.45 &
    
13. dim 5.95 5.95 5.95 5.95 size 11 50 15 10 group '加固圈'
    
14. ;==================管棚=======================
    
15. ge zo cs p0 0 0 0 p1 9.45 0 0 p2 0 50 0 p3 0 0 9.45 &
    
16. dim 5.95 5.95 5.95 5.95 size 7 50 15 10
    
17. group '管棚加固' rang group '加固圈' not
    
18. ;==========5.5加固圈外一层围岩-上部圆形===========
    
19. ge zo radc p0 0 0 0 p1 30 0 0 p2 0 50 0 p3 0 0 20.95 p4 30 50 0 &
    
20. p5 0 50 20.95 p6 30 0 20.95 p7 30 50 20.95 p8 11.45 0 0 p9 0 0 11.45 &
    
21. p10 11.45 50 0 &
    
22. p11 0 50 11.45 size 5 50 15 10
    
23. ;===========5.5加固圈外一层围岩-下部圆形==========
    
24. ge zo radc p0 0 50 0 p1 30 50 0 p2 0 0 0 p3 0 50 -20.95 &
    
25. p4 30 0 0 p5 0 0 -20.95 p6 30 50 -20.95 &
    
26. p7 30 0 -20.95 p8 11.45 50 0 p9 0 50 -11.45 &
    
27. p10 11.45 0 0 p11 0 0 -11.45 size 5 50 15 10
    
28. ;========================初衬==================
    
29. ge zo cs p0 0 0 0 p1 5.95 0 0 p2 0 50 0 p3 0 0 5.95 &
    
30. dim 5.7 5.7 5.7 5.7 size 1 50 15 6 group '初衬' fill group '隧道' &
    
31. ratio 1 1 1 0.9
    
32. 
    
33. ge zo cs p0 0 50 0 p1 5.95 50 0 p2 0 0 0 p3 0 50 -5.95 &
    
34. dim 5.7 5.7 5.7 5.7 size 1 50 15 6 group '初衬' fill group '隧道' &
    
35. ratio 1 1 1 0.9
    
36. 
    
37. ;====================堤坝======================
    
38. ge zo b p0 0 0 20.95 p1 add 30 0 0 p2 add 0 50 0 &
    
39. p3 add 0 20 6 p4 add 30 50 0 p5 add 0 30 6 &
    
40. p6 add 30 20 6 p7 add 30 30 6 size 8 50 6
    
41. group '围岩' ran group '加固圈' not group '管棚加固' not group '初衬' not group '隧道' not
    
42. gen merge 1e-2
    
43. ;===================赋予材料属性================
    
44. 
    
45. mod mohr
    
46. prop dens 2200 bulk 2.08e8 sh 1.19e7 co 4.5e4 fric 30 tens 2e4 ran z -20.95 5.95
    
47. prop dens 1900 bulk 8e6 sh 5.5e6 co 2e4 fric 18 tens 2e4 ran z 26.45 5.95
    
48. ;===================边界条件====================
    
49. ini szz -9.854e5 grad 0,0,2.2e4 ran z -20.95 5.95
    
50. ini szz -3.895e5 grad 0,0,1.9e4 ran z 5.95 26.45
    
51. 
    
52. ini sxx -6.2e5 grad 0,0,1.4e4 ran z -20.95 5.95
    
53. ini sxx -2.4e5 grad 0,0,1.4e4 ran z 5.95 26.45
    
54. 
    
55. ini syy -6e5 grad 0,0,1.3e3 ran z -20.95 5.95
    
56. ini syy -2.3e5 grad 0,0,1.2e3 ran z 5.95 26.45
    
57. 
    
58. fix z ran z -20.85 -30.05
    
59. fix x ran x -0.1 0.1
    
60. fix x ran x 29.9 30.1
    
61. fix y ran y -.1 .1
    
62. fix y ran y 49.9 50.1
    
63. ;===================形成初始应力场===============
    
64. ini pp 1.3e5 grad 0 0 -1e4 ran z 0 -20.95
    
65. ini pp 1.3e5 grad 0 0 -1e4 ran z 0 13
    
66. model fl_iso
    
67. set grav 0 0 -10
    
68. ini fdens 1000
    
69. 
    
70. 
    
71. 
    
72. 
    
73. so elas
    
74. sav hon.sav
    
75. 
    
76. 
    
77. 
    
78. 
    
79. 
    
80. 
    
81. 
    

复制代码

我把命令流发上来吧 请指点

hjm460727715

]我把命令流发上来吧 请指点
n
plo show
plo blo group
conf dyn fluid;设置动态与流体算法
set dyn off fluid off;关闭动态与流体算法

;================上部加固圈====================
ge zo cs p0 0 0 0 p1 11.45 0 0 p2 0 50 0 p3  0 0 11.45 &
dim 9.45 9.45 9.45 9.45 size 4 50 15 10 group '加固圈'
;================下部加固圈====================
ge zo cs p0 0 50 0 p1 11.45 50 0 p2 0 0 0 p3 0 50 -11.45 &
dim 5.95 5.95 5.95 5.95 size 11 50 15 10 group '加固圈'
;==================管棚=======================
ge zo cs p0 0 0 0 p1 9.45 0 0 p2 0 50 0 p3  0 0 9.45 &
dim 5.95 5.95  5.95  5.95  size 7 50 15 10
group '管棚加固'  rang  group  '加固圈'  not
;==========5.5加固圈外一层围岩-上部圆形===========
ge zo radc p0 0 0 0 p1 30 0 0 p2 0 50 0 p3 0 0 20.95 p4 30 50 0 &
p5 0 50 20.95 p6 30 0 20.95 p7 30 50 20.95  p8 11.45 0 0 p9 0 0 11.45 &
p10 11.45 50 0 &
p11 0 50 11.45 size 5 50 15 10
;===========5.5加固圈外一层围岩-下部圆形==========
ge zo radc p0 0 50 0 p1 30 50 0 p2 0 0 0 p3 0 50 -20.95  &
p4 30 0  0  p5 0 0 -20.95 p6 30 50 -20.95 &
p7 30 0  -20.95 p8 11.45 50 0 p9 0 50 -11.45 &
p10 11.45 0 0 p11 0 0 -11.45  size 5 50 15 10  
;========================初衬==================
ge zo cs p0 0 0 0 p1 5.95 0 0 p2 0 50 0 p3 0 0 5.95 &
dim 5.7 5.7 5.7 5.7  size 1 50 15 6 group '初衬' fill group '隧道' &
ratio  1 1 1  0.9
ge zo cs p0 0 50 0 p1 5.95 50 0  p2 0 0 0 p3 0 50 -5.95 &
dim 5.7 5.7 5.7 5.7  size 1 50 15 6 group '初衬' fill group '隧道' &
ratio  1 1 1  0.9
;====================堤坝======================
ge zo b p0 0 0 20.95 p1 add 30 0 0 p2 add  0 50 0 &
p3 add 0 20 6 p4 add 30 50 0 p5 add 0 30 6 &
p6 add 30 20 6 p7 add 30 30 6 size 8 50 6
group '围岩' ran group '加固圈' not  group '管棚加固' not group '初衬' not group '隧道' not
gen merge 1e-2
;===================赋予材料属性================
mod mohr
prop dens 2200 bulk 2.08e8 sh 1.19e7 co 4.5e4 fric 30 tens 2e4 ran z -20.95 5.95
prop dens 1900 bulk  8e6 sh 5.5e6 co 2e4 fric 18 tens  2e4 ran z 26.45 5.95
;===================边界条件====================
ini szz -9.854e5 grad 0,0,2.2e4 ran z -20.95 5.95
ini szz -3.895e5 grad 0,0,1.9e4 ran z 5.95 26.45
ini sxx -6.2e5 grad 0,0,1.4e4 ran z -20.95 5.95
ini sxx -2.4e5 grad 0,0,1.4e4 ran z  5.95 26.45
ini syy -6e5 grad 0,0,1.3e3 ran z -20.95 5.95
ini syy -2.3e5 grad 0,0,1.2e3 ran z 5.95 26.45
fix z ran z -20.85 -30.05
fix x ran x -0.1  0.1
fix x ran x 29.9  30.1
fix y ran y  -.1  .1
fix y ran y 49.9 50.1
;===================形成初始应力场===============
ini pp 1.3e5 grad 0 0 -1e4 ran z 0 -20.95
ini pp 1.3e5 grad 0 0 -1e4 ran z 0 13
model fl_iso
set grav 0 0 -10
ini fdens 1000


so  elas
sav hon.sav

hjm460727715

怎么 会没人搭理我呢？ 急盼 哪位高手给予帮助 谢谢先！

lookcity

大致看了一下，没仔细看。我初步的意见是：要先形成初始应力场之后再进行开挖和支护，而不是一开始就支护好然后再形成初始应力场。

hjm460727715

我自己把它个搞定了 呵呵 我发现 出现较大反弹的只有基底上极小的一层土层 把抗拉强度调大点就可以了哦！ 还有用FALC 建模 开挖后形成的 台阶面总是 参差不齐 我用ANSYS 建模后导入 情况好多了

hjm460727715

还有 要进一步解释一下 为什么开挖中会有土体向侧壁位移的现象
我是分为三台阶开挖 上台阶开挖的时候断面是一个 圆心角小于180的圆弧
在较大的（相对于侧向应力）竖向应力作用下自然会在递补出现拉应力 最明显的
拉出现象发生在底脚处。再开挖了整个断面后这种现象就不存在了 呵呵

zxb911

感谢楼主分享经验，解决了我的部分问题