各位大神，小弟我想模拟下子施加围压状态下，TBM在不同间距侵深下的侵入试验。下面是我模拟的，当然感觉好多问题，比如不知道怎么结束之类的。还请指教，当然也可以有偿指导。
new
;建立单元体
gen z b size 25 20 10 &
p0=(0,0,0) p1=(0.25,0,0) p2=(0,0.2,0) p3=(0,0,0.10)
;定义结构体和参数赋值
model mohr
property young 1.71e10 poisson 0.25 dens 2260
prop cohesion 2.95e7 friction 41.6 tension 9.2e6
;固定z 面
fix x y z range z -0.01 0.01
fix x y z range z 0.09 0.11
;施加起始速度变量
ini zvel -8.3e-6 range x 0.0895 0.0905 z 0.09 0.11
ini zvel -8.3e-6 range x 0.1595 0.1605 z 0.09 0.11
;施加围压
apply syy -2.5e5 range y -0.01 0.01
apply syy -2.5e5 range y 0.19 0.21
apply sxx -3e5 range x -0.01 0.01
apply sxx -3e5 range x 0.249 0.251
hist id 1 gp zdisp 0.08 0.1 0.1
hist id 2 gp zforce 0.08 0.1 0.1
hist id 3 gp zdisp 0.08 0.1 0.092
step 5000
set mech ratio 1e-4
solve
save 5-6.save
new
;建立单元体
gen z b size 25 20 10 &
p0=(0,0,0) p1=(0.25,0,0) p2=(0,0.2,0) p3=(0,0,0.10)
;定义结构体和参数赋值
model mohr
property young 1.71e10 poisson 0.25 dens 2260
prop cohesion 2.95e7 friction 41.6 tension 9.2e6
;固定z 面
fix x y z range z -0.01 0.01
fix x y z range z 0.09 0.11
;施加起始速度变量
ini zvel -8.3e-6 range x 0.0895 0.0905 z 0.09 0.11
ini zvel -8.3e-6 range x 0.1595 0.1605 z 0.09 0.11
;施加围压
apply syy -2.5e5 range y -0.01 0.01
apply syy -2.5e5 range y 0.19 0.21
apply sxx -3e5 range x -0.01 0.01
apply sxx -3e5 range x 0.249 0.251
hist id 1 gp zdisp 0.08 0.1 0.1
hist id 2 gp zforce 0.08 0.1 0.1
hist id 3 gp zdisp 0.08 0.1 0.092
step 5000
set mech ratio 1e-4
solve
save 5-6.save
