gen zone radcyl p0 0,0,0 p1 24,0,0 p2 0,6,0 p3 0,0,24 &
dim 9.5 9.5 9.5 9.5 size 5 1 12 8 rat 1 1 1 1 group 围岩
gen zone cshell p0 0 0 0 p1 9.5,0,0 p2 0,6,0 p3 0,0,9.5 &
dim 9 9 9 9 size 5 1 12 8 rat 1 1 1 1 group 初期支护
gen zone cshell p0 0 0 0 p1 9,0,0 p2 0,6,0 p3 0,0,9 &
dim 8.6 8.6 8.6 8.6 size 5 1 12 8 rat 1 1 1 1 group 二衬 fill group 原岩
gen zone radcyl p0 0,0,0 p1 0,0,-24 p2 0,6,0 p3 24,0,0 &
dim 1.8 9 1.8 9 size 5 1 12 8 rat 1 1 1 1 group 围岩
gen zone cshell p0 0 0 0 p1 0,0,-1.8 p2 0,6,0 p3 9,0,0 &
dim 1.5 8.8 1.5 8.8 size 5 1 12 8 rat 1 1 1 1 group 仰拱初期支护
gen zone cshell p0 0,0,0 p1 0,0,-1.5 p2 0,6,0 p3 8.6,0,0 &
dim 1.2 8.6 1.2 8.6 size 5 1 12 8 rat 1 1 1 1 group 仰拱二衬 fill group 仰拱原岩
gen zone re dip 90 dd 270 or 0 0 0
plot add sur yellow
plot show
gen zone re dip 90 dd 270 or -24 0 0;隧道区域
gen zone brick p0 -72,0,-52 p1 24,0,-52 p2 -72,6,-52 p3 -72,0,-24 &
size 24 1 7 rat 1 1 1 1;下部扩展区域
attach face range x -24 24 z -23.9 -24.1
gen zone brick p0 24,0,-52 p1 64,0,-52 p2 24,6,-52 p3 24,0,24 &
size 10 1 19 rat 1 1 1 1;右边扩展区域
attach face range x 23.9 24.1 z -52 24
gen zone brick p0 -112,0,-52 p1 -72,0,-52 p2 -112,6,-52 p3 -112,0,24 &
size 10 1 19 rat 1 1 1 1;左边扩展区域
attach face range x -71.9 -72.1 z -52 24
gen zone brick p0 -112,0,24 p1 64,0,24 p2 -112,6,24 p3 -112,0,35 &
size 44 1 4 rat 1 1 1 1;上部扩展区域1
attach face range x -112 64 z 23.9 24.1
gen zone wedge p0 64,0,35 p1 64,0,43 p2 64,6,35 p3 -112,0,35 &
size 3 1 44 rat 1 1 1 1;上部扩展区域2
attach face range x -112 64 z 34.9 35.1
set gravity 0 0 -9.8
fix x y z ran z -52.1 -51.9
fix x ran x -112.1 -111.9
fix x ran x 63.9 64.1
fix y ran y -0.1 0.1
fix y ran y 5.9 6.1
model mohr
ini dens 1700
prop bu 1.11e9 sh 3.7e8 fric 20 coh 5.0e4
prop bu 1.11e9 sh 3.7e8 fric 20 coh 5.0e4
sel pile id=1 beg -30,3,38.727 end -30,3,33.927 nseg=4
sel pile id=2 beg -13,3,39.5 end -13,3,34.7 nseg=4;铁塔基础
sel pile prop dens 2500 Emod 3.0e10 Nu 0.2 XCArea 2.1 Per 5.14
solve
plot cont szz
这是我做的源程序，求指点。。。自己一直找不出原因。。。
