一个关于隧道开挖的模拟命令流，紧急求助！

liuxiaowenid

本人是初学者，做了一个隧道开挖的模拟，但是在计算的过程中发现很多问题，最主要的问题就是最大主应力云图围岩都是一个颜色，也就是说围岩的压应力都差不多一样大，不知道怎么回事，请那位大侠帮忙看看。

liuxiaowenid

plot show
plot block group
gen zon radcyl p0 0 0 0 p1 8 0 0 p2 0 30 0 p3 0 0 8 &
size 4 30 12 4 dim 4.25 5.25 4.25 5.25 rat 1 1 1 1 group 围岩
gen zone radcyl p0 0 0 0 p1 0 0 -6.3 p2 0 30 0 p3 8 0 0 &
size 4 30 12 4 dim 3.55 4.25 3.55 4.25 rat 1 1 1 1 group 围岩
gen zone brick p0 8 0 0 p1 20 0 0 p2 8 30 0 p3 8 0 8 &
size 11 30 6 group 围岩
gen zone brick p0 8 0 -6.3 p1 20 0 -6.3 p2 8 30 -6.3 p3 8 0 0 &
size 11 30 6 group 围岩
gen zone brick p0 0 0 8 p1 8 0 8 p2 0 30 8 p3 0 0 15 &
size 6 30 7 group 围岩1
gen zone brick p0 8 0 8 p1 20 0 8 p2 8 30 8 p3 8 0 15 &
size 11 30 7 group 围岩1
gen zone brick p0 0 0 15 p1 8 0 15 p2 0 30 15 p3 0 0 17 &
size 6 30 2 group 围岩2
gen zone brick p0 8 0 15 p1 20 0 15 p2 8 30 15 p3 8 0 17 &
size 11 30 2 group 围岩2
gen zone brick p0 0 0 -7 p1 8 0 -7 p2 0 30 -7 p3 0 0 -6.3 &
size 6 30 1 group 围岩3
gen zone brick p0 8 0 -7 p1 20 0 -7 p2 8 30 -7 p3 8 0 -6.3 &
size 11 30 1 group 围岩3
gen zone brick p0 0 0 -16 p1 8 0 -16 p2 0 30 -16 p3 0 0 -7 &
size 6 30 9 group 围岩4
gen zone brick p0 8 0 -16 p1 20 0 -16 p2 8 30 -16 p3 8 0 -7 &
size 11 30 9 group 围岩4
gen zone reflect normal -1 0 0
gen zon cshell p0 0 0 0 p1 4.25 0 0 p2 0 30 0 p3 0 0 5.25 &
size 1 30 12 4 dim 3.95 4.95 3.95 4.95 rat 1 1 1 1 group 初期支护
gen zon cshell p0 0 0 0 p1 0 0 5.25 p2 0 30 0 p3 -4.25 0 0 &
p4 0 30 5.25 p5 -4.25 30 0 &
size 1 30 12 4 dim 4.95 3.95 4.95 3.95 rat 1 1 1 1 group 初期支护
gen zone cshell p0 0 0 0 p1 3.95 0 0 p2 0 30 0 p3 0 0 4.95 &
size 1 30 12 6 dim 3.45 4.45 3.45 4.45 rat 1 1 1 1 group 二衬 fill group 原岩
gen zone cshell p0 0 0 0 p1 0 0 4.95 p2 0 30 0 p3 -3.95 0 0 &
p4 0 30 4.95 p5 -3.95 30 0 &
size 1 30 12 6 dim 4.45 3.45 4.45 3.45 rat 1 1 1 1 group 二衬 fill group 原岩
gen zone cshell p0 0 0 0 p1 0 0 -3.55 p2 0 30 0 p3 4.25 0 0 &
size 1 30 12 4 dim 3.25 3.95 3.25 3.95 rat 1 1 1 1 group 仰拱初期支护
gen zone cshell p0 0 0 0 p1 -4.25 0 0 p2 0 30 0 p3 0 0 -3.55 &
p4 -4.25 30 0 p5 0 30 -3.55 &
size 1 30 12 4 dim 3.95 3.25 3.95 3.25 rat 1 1 1 1 group 仰拱初期支护
gen zone cshell p0 0 0 0 p1 0 0 -3.25 p2 0 30 0 p3 3.95 0 0 &
size 1 30 12 6 dim 2.75 3.45 2.75 3.45 rat 1 1 1 1 group 仰拱二衬 fill group 仰拱原岩
gen zone cshell p0 0 0 0 p1 -3.95 0 0 p2 0 30 0 p3 0 0 -3.25 &
p4 -3.95 30 0 p5 0 30 -3.25 &
size 1 30 12 6 dim 3.45 2.75 3.45 2.75 rat 1 1 1 1 group 仰拱二衬 fill group 仰拱原岩
save model.sav

m e
prop den 1680 bulk 4.58e6 shear 2.12e6  range z 17 15 y 0 15
prop den 1640 bulk 3.77e6 shear 1.95e6  range z 15 7 y 0 15
prop den 1750 bulk 4.46e6 shear 2.06e6  range z 7 3 y 0 15
prop den 1650 bulk 4.20e6 shear 2.52e6  range z 3 -4 y 0 15
prop den 1820 bulk 4.27e6 shear 2.44e6  range z -4 -7 y 0 15
prop den 1860 bulk 9.33e6 shear 4.31e6  range z -7 -16 y 0 15
prop den 1680 bulk 4.58e6 shear 2.12e6  range z 17 15 y 18 30
prop den 1640 bulk 3.77e6 shear 1.95e6  range z 15 7 y 18 30
prop den 1750 bulk 4.46e6 shear 2.06e6  range z 7 3 y 18 30
prop den 1650 bulk 4.20e6 shear 2.52e6  range z 3 -4 y 18 30
prop den 1820 bulk 4.27e6 shear 2.44e6  range z -4 -7 y 18 30
prop den 1860 bulk 9.33e6 shear 4.31e6  range z -7 -16 y 18 30
prop den 400  bulk 4.17e6 shear 0.89e6  range  y 15 18
fix z range z -16.1 -15.9
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix y range y 29.9 30.1
fix y range y -0.1 0.1
apply nstress -2e4 range z 14.9 15.1
set grav 0 0 -10
hist unbal
solve
save nature.sav





ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
ini state=0
model mohr
prop den 1680 bulk 4.58e6 shear 2.12e6 cohe 11e3 fric 12.5 range z 17 15 y 0 15
prop den 1640 bulk 3.77e6 shear 1.95e6 cohe 26e3 fric 27   range z 15 7 y 0 15
prop den 1750 bulk 4.46e6 shear 2.06e6 cohe 25e3 fric 22   range z 7 3 y 0 15
prop den 1650 bulk 4.20e6 shear 2.52e6 cohe 30e3 fric 22   range z 3 -4 y 0 15
prop den 1820 bulk 4.27e6 shear 2.44e6 cohe 23e3 fric 22   range z -4 -7 y 0 15
prop den 1860 bulk 9.33e6 shear 4.31e6           fric 32   range z -7 -16 y 0 15
prop den 1680 bulk 4.58e6 shear 2.12e6 cohe 11e3 fric 12.5 range z 17 15 y 18 30
prop den 1640 bulk 3.77e6 shear 1.95e6 cohe 26e3 fric 27   range z 15 7 y 18 30
prop den 1750 bulk 4.46e6 shear 2.06e6 cohe 25e3 fric 22   range z 7 3 y 18 30
prop den 1650 bulk 4.20e6 shear 2.52e6 cohe 30e3 fric 22   range z 3 -4 y 18 30
prop den 1820 bulk 4.27e6 shear 2.44e6 cohe 23e3 fric 22   range z -4 -7 y 18 30
prop den 1860 bulk 9.33e6 shear 4.31e6           fric 32   range z -7 -16 y 18 30
prop den 400  bulk 4.17e6 shear 0.89e6 cohe 5e3  fric 18   range  y 15 18
model null range group 初期支护 y 0 2
model null range group 二衬 y 0 2
model null range group 仰拱初期支护 y 0 2
model null range group 仰拱二衬 y 0 2
model null range group 原岩 y 0 2
model null range group 仰拱原岩 y 0 2
m e range group 初期支护 y 0 2
prop den 2500 bulk 13.7e9 shear 8.2e9 range group 初期支护 y 0 2
m e range group 二衬 y 0 2
prop den 2550 bulk 14.7e9 shear 8.8e9 range group 仰拱二衬 y 0 2
m e range group 仰拱初期支护 y 0 2
prop den 2500 bulk 13.7e9 shear 8.2e9 range group 仰拱初期支护 y 0 2
m e range group 仰拱二衬 y 0 2
prop den 2550 bulk 14.7e9 shear 8.8e9 range group 仰拱二衬 y 0 2
hist unbal
hist gp zdis     0,5,17   
hist gp zdis     4,5,17   
hist gp zdis     8,5,17   
hist gp zdis    12,5,17   
hist gp zdis    16,5,17   
hist gp zdis    20,5,17   
hist gp zdis     0,5,5.25
hist gp zdis     4,5,5.25
hist gp zdis     8,5,5.25
hist gp zdis    12,5,5.25
hist gp zdis    16,5,5.25
hist gp zdis    20,5,5.25
hist gp zdis     0,5,4.45   
hist gp zdis     0,5,-2.75  
hist gp zdis     0,5,-3.55  
hist gp xdis  4.25,5,0      
hist gp xdis  3.45,5,0   
hist gp xdis -4.25,5,0   
hist gp xdis -3.45,5,0   
hist gp zdis     0,17,17  
hist gp zdis     4,17,17  
hist gp zdis     8,17,17  
hist gp zdis    12,17,17  
hist gp zdis    16,17,17  
hist gp zdis    20,17,17  
hist gp zdis     0,17,5.25
hist gp zdis     4,17,5.25
hist gp zdis     8,17,5.25  
hist gp zdis    12,17,5.25  
hist gp zdis    16,17,5.25  
hist gp zdis    20,17,5.25  
hist gp zdis     0,17,4.45  
hist gp zdis     0,17,-2.75  
hist gp zdis     0,17,-3.55
hist gp xdis  4.25,17,0     
hist gp xdis  3.45,17,0     
hist gp xdis -4.25,17,0     
hist gp xdis -3.45,17,0     
hist gp zdis     0,27,17   
hist gp zdis     4,27,17   
hist gp zdis     8,27,17   
hist gp zdis    12,27,17   
hist gp zdis    16,27,17   
hist gp zdis    20,27,17   
hist gp zdis     0,27,5.25   
hist gp zdis     4,27,5.25  
hist gp zdis     8,27,5.25   
hist gp zdis    12,27,5.25  
hist gp zdis    16,27,5.25   
hist gp zdis    20,27,5.25   
hist gp zdis     0,27,4.45   
hist gp zdis     0,27,-2.75   
hist gp zdis     0,27,-3.55   
hist gp xdis  4.25,27,0   
hist gp xdis  3.45,27,0   
hist gp xdis -4.25,27,0   
hist gp xdis -3.45,27,0   
solve
save 2m.sav

liuxiaowenid

真的很急，希望大家能帮帮我

我是白桦树

刚才在群里没人给你解答呀？

stuy_1001823

:P:lol:(

ycgsxdl

你干嘛不把命令流打包发上来？？