巷道命令流FLAC3D

new
;新建

gen zon radcyl p0 0 0 0 p1 20 0 0 p2 0 10 0 p3 0 0 20 p4 20 10 0 p5 0 10 20 p6 20 0 20 p7 20 10 20 p8 2 0 0 p9 0 0 2 p10 2 10 0 p11 0 10 2 dimension 2 2 2 2 size 2 10 10 18
;建立巷道

plo con szz
;显示垂直应力云图

gen zon reflect ori 0 0 0 norm 1 0 0
;以X轴镜像

gen zon reflect ori 0 0 0 norm 0 0 1
;以Z轴镜像

m m
;模型选摩尔库伦准则


pro bulk 3e9 shear 1.38e9 fric 30 coh 1.8e6 ten 5e6 density 2810 range z -20 20
;赋岩性参数值


fix z range z -20.1 -19.9
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix y range y -0.1 0.1
fix y range y 9.9 10.1
;固定边界，Z上边界加压，不固定

ini szz -12.5e6 grad 0 0 2.5e4
ini sxx -25e6 grad 0,0,5e4
ini syy -25e7 grad 0,0,5e5
;初始应力，所处深度的应力


apply szz -12.5e6 range z 19.9 20.1
apply sxx -25e6 grad 0 0 5e4 range x -20.1 -19.9
apply sxx -25e6 grad 0 0 5e4 range x 19.9 20.1
apply syy -25e7 grad 0 0 5e5 range y -0.1 0.1
apply syy -25e7 grad 0 0 5e5 range y 9.9 10.1
;施加应力于模型上




set gravity 0 0 -10
;设置重力加速度


set large
;设置大变形，可取消



hist id=1 unbal
;记录最大不平衡力


plo sho
;显示


ini xdis 0 ydis 0 zdis 0
;初始速度




solve
;求解


save 001xiugai.sav
;保存，以便下次调用