new
; 产生网格体
; 隧道上部混凝土衬砌
gen zon cshell p0 0 0 0 p1 6.3 0 0 p2 0 200 0 p3 0 0 6.3 &
dim 5.6 5.6 5.6 5.6 size 1 200 10
group 'concrete liner'
; 隧道上部
gen zon cyl p0 0 0 0 p1 5.6 0 0 p2 0 200 0 p3 0 0 5.6 &
size 4 200 10
group tunnel range group 'concrete liner' not
; 隧道下部
gen zone brick p0 0 0 -4.5 p1 add 6.3 0 0 p2 add 0 200 0 &
p3 add 0 0 4.5 size 6 200 2
; 周围岩石 (10 primitives)
gen zon radcyl p0 0 0 0 p1 26.3 0 0 p2 0 200 0 p3 0 0 25 &
dim 6.3 6.3 6.3 6.3 size 5 200 10 8 rat 1 1 1 1.3
gen zone brick p0 6.3 0 -4.5 p1 26.3 0 -15 p2 add 0 200 0 p3 6.3 0 0 &
p4 26.3 200 -15 p5 6.3 200 0 p6 26.3 0 0 p7 26.3 200 0 &
size 8 200 3 ratio 1.3 1 1
gen zone brick p0 0 0 -15 p1 add 26.3 0 0 p2 add 0 200 0 &
p3 0 0 -4.5 p4 26.3 200 -15 p5 0 200 -4.5 p6 6.3 0 -4.5 &
p7 6.3 200 -4.5 size 6 200 8 rat 1 1 0.7692307692307692
gen zon brick p0 0 0 25 p1 add 26.3 0 0 p2 add 0 200 0 &
p3 add 0 0 10 size 5 200 2
gen zon bric p0 26.3 0 25 p1 add 16.3 0 0 p2 add 0 200 0 &
p3 add 0 0 10 size 2 200 2 rat 2 1 1
gen zon bric p0 26.3 0 -15 p1 add 16.3 0 0 p2 add 0 200 0 &
p3 add 0 0 40 size 2 200 8 rat 2 1 1
gen zon bric p0 26.3 0 -40 p1 add 16.3 0 0 p2 add 0 200 0 &
p3 add 0 0 25 size 2 200 2 rat 2 1 0.5
gen zon bric p0 0 0 -40 p1 add 26.3 0 0 p2 add 0 200 0 &
p3 add 0 0 25 size 6 200 2 rat 1 1 0.5
group rock range group 'concrete liner' not group tunnel not
;完成视图
plot show
plot add surface yellow
plot add axes black
; 指派 Mohr-Coulomb 材料模型
mo mohr
pro bulk 25e8 she 5.4e8 fric 25 coh 10e4 ten 0 dil 0 &
range z -50 35
sel liner range group 'concrete liner' x 5.6 6.3 y 0 200 z 5.6 6.3
sel liner prop iso=(25e9,0.15) thick=0.7
sel liner prop cs_nk=8e8 cs_sk=8e8 &
cs_ncut=4e6 cs_scoh=4e6 cs_scohres=2e6 cs_sfric=25
fix x y z
sel node fix x y z
; 边界条件
fix x range x -.1 .1
fix x range x 42.5 42.7
fix z range z -40.1 -39.9
fix y range y -.1 .1
fix y range y 199.9 200.1
; 初始应力状态
set grav 0 0 -10
ini density 2200
ini szz -770e3 grad 0 0 22000
ini sxx -770e3 grad 0 0 22000
ini syy -385e3 grad 0 0 11000
; 监测变量
hist id=1 unbal
hist id=2 gp zdisp 0 0 6.3
hist id=3 gp xdisp 6.3 0 0
hist id=4 gp zdisp 0 0 0
hist id=5 gp zdisp 0 0 35
hist id=6 gp zdisp 0 30 6.3
hist id=7 gp xdisp 6.3 30 0
hist id=8 gp zdisp 0 30 0
hist id=9 gp zdisp 0 30 35
hist id=10 gp zdisp 0 12 35
hist id=11 gp zdisp 0 18 35
hist id=12 gp zdisp 0 24 35
hist id=13 gp zdisp 0 36 35
hist id=14 gp zdisp 5.6 30 35
hist id=15 gp zdisp 10 30 35
hist id=16 zone smin 0 0 6.3
hist id=17 zone smin 6.3 0 0
hist id=18 zone smin 0 0 0
hist id=19 zone smin 0 0 35
hist id=20 zone smin 0 30 6.3
hist id=21 zone smin 6.3 30 0
hist id=22 zone smin 0 30 0
hist id=23 zone smin 0 30 35
hist id=24 zone smin 0 12 35
hist id=25 zone smin 0 18 35
hist id=26 zone smin 0 24 35
hist id=27 zone smin 0 36 35
hist id=28 zone smin 5.6 30 35
hist id=29 zone smin 10 30 35
hist id=30 zone ssi 0 0 6.3
hist id=31 zone ssi 6.3 0 0
hist id=32 zone ssi 0 0 0
hist id=33 zone ssi 0 0 35
hist id=34 zone ssi 0 30 6.3
hist id=35 zone ssi 6.3 30 0
hist id=36 zone ssi 0 30 0
hist id=37 zone ssi 0 30 35
hist id=38 zone ssi 0 12 35
hist id=39 zone ssi 0 18 35
hist id=40 zone ssi 0 24 35
hist id=41 zone ssi 0 36 35
hist id=42 zone ssi 5.6 30 35
hist id=43 zone ssi 10 30 35
hist id=44 sel recover sres mx 0 0 6.3 ;开始出现错误
hist id=45 sel recover sres mx 6.3 0 0
hist id=46 sel recover sres mx 0 0 0
hist id=47 sel recover sres mx 0 0 35
hist id=48 sel recover sres mx 0 30 6.3
hist id=49 sel recover sres mx 6.3 30 0
hist id=50 sel recover sres mx 0 30 0
hist id=51 sel recover sres mx 0 30 35
hist id=52 sel recover sres mx 0 12 35
hist id=53 sel recover sres mx 0 18 35
hist id=54 sel recover sres mx 0 24 35
hist id=55 sel recover sres mx 0 36 35
hist id=56 sel recover sres mx 5.6 30 35
hist id=57 sel recover sres mx 10 30 35
hist id=58 sel recover sres qx 0 0 6.3
hist id=59 sel recover sres qx 6.3 0 0
hist id=60 sel recover sres qx 0 0 0
hist id=61 sel recover sres qx 0 0 35
hist id=62 sel recover sres qx 0 30 6.3
hist id=63 sel recover sres qx 6.3 30 0
hist id=64 sel recover sres qx 0 30 0
hist id=65 sel recover sres qx 0 30 35
hist id=66 sel recover sres qx 0 12 35
hist id=67 sel recover sres qx 0 18 35
hist id=68 sel recover sres qx 0 24 35
hist id=69 sel recover sres qx 0 36 35
hist id=70 sel recover sres qx 5.6 30 35
hist id=71 sel recover sres qx 10 30 35
plot create stress
plot set title text '应力云图
plot set back white
plot add con sxx outline on
plot add con szz outline on
plot add con sxz outline on
plot create strain
plot set title text '应变云图
plot set back white
plot add con ssr outline on average
plot add con ssi outline on average
plot create dis-step
plot set title text '监测点位移图
plot set back white
plot add his 2 3 4 5 6 7 8 9 10 11 12 13 14 15 vs step
plot create p-stress-step
plot set title text '监测点最大主应力
plot set back white
plot add his 16 17 18 19 20 21 22 23 24 25 26 27 28 29 vs step
plot create shear-strain-step
plot set title text '监测点剪应变
plot set back white
plot add his 30 31 32 33 34 35 36 37 38 39 40 41 42 43 vs step
plot create moment-step
plot set title text '监测点弯矩
plot set back white
plot add his 44 45 46 47 48 49 50 51 52 53 54 55 56 57 vs step
plot create shear-step
plot set title text '监测点剪力
plot set back white
plot add his 58 59 60 61 62 63 64 65 66 67 68 69 70 71 vs step
