命令流
new
gen zone brick size 8 1 22 p0 0 0 -6 p1 4.2 0 -6 p2 0 1 -6 p3 0 0 5
group soil
fix x range x 4.2
fix x range x 0
fix y
fix z range z -6
config creep
def burcal
bkk=bme/(3*(1-2*b1)) ;‘bkk——k(体积模量);bme——E（弹性模量）；b1——泊松比；
bk1=bke/(2*(1+b1)) ;‘bke——
bk2=bme/(2*(1+b1)) ; ‘bk2——G（剪切模量）；bme——E（弹性模量）；b1——泊松比；kshear 开尔文抗剪模量, GK
bvis1=vis1 ;‘vis——kviscosity 开尔文粘性系数, ηK
bvis2=vis2
end
set bme=1.35714e9
set bke=3.25714e9
set b1=0.35
set vis1=1.35714e9;d
set vis2=19.2857e9;d
burcal
model burger range group soil any
prop b b1 ksh bke msh bme kvis vis1 mvis vis2 range group soil any
ini dens=2600 range group soil any
set grav 0 0 -10;；重力加速度
;ini syy -4.5e4 grad 0 0 0.9e4;；初始应力平衡
;ini sxx -4.5e4 grad 0 0 0.9e4
;ini szz -9e4 grad 0 0 1.8e4
solve
pause;初试平衡求解结果
ini xdis 0 ydis 0 zdis 0 ;；初始位移为零
set creep mindt=1e-2 maxdt=1
set creep on;开启蠕变计算
solve age 6;；蠕变计算
pause
solve age 120 ;开挖第2步求解结果
