new
gen zone cyl p0 0 0 0 p1 0.05 0 0 p2 0 0.2 0 p3 0 0 0.05 size 8 12 8
gen zone reflect norm 1 0 0
gen zone reflect norm 0 0 1
model mohr
prop bulk 6.3e9 shear 4.5e9 coh 1e6 fric 35 ten 0.8e6
fix x y z range y 0.1
fix x y z range y 0
ini yvel -1e-7 range y 0.1
ini yvel 1e-7 rang y 0
def ax_str
str=0
pnt=gp_head
loop while pnt # null
if gp_ypos(pnt)>0.095 then
str=str+gp_yfunbal(pnt)
endif
pnt=gp_next(pnt)
endloop
ax_str=str/pi/0.025/0.025
end
hist n 10
hist gp ydisp 0 0.1 0
hist gp xdisp 0 0.1 0
hist ax_str
plot hist 3 vs -1
step 40000为什么参数的变化，会导致应力应变曲线形式发生很大变化
