自己用FLAC3D做的一个静三轴试验模拟出了点问题，请高手指点！
这次模拟没有考虑静三轴的固结阶段，用的是摩尔模型。
最后得不到轴向应力-轴向位移曲线，希望有高人现身。
先谢谢了~
new
define change_modulus
sample_young=8.3e7
sample_poisson=0.33
sample_cohesion=21.71
sample_friction=12.76
sample_bulk=sample_young/3*(1-2*sample_poisson)
sample_shear=sample_young/2*(1+sample_poisson)
end
;grid generation
generate zone cylinder p0 0 0 0 p1 19.55e-2 0 0 p2 0 80e-2 0 p3 0 0 19.55e-2 size 5 20 5
generate zone reflect normal 1 0 0
generate zone reflect normal 0 0 1
;choose model and set property
model mohr
change_modulus
property bulk=sample_bulk shear=sample_shear
property cohesion=sample_cohesion friction=sample_friction
;boundary
fix x y z range y=80e-2
fix x y z range y=0
;apply confining pressure
apply nstress -1e5 range cylinder end1 19.55e-2 0 0 end2 19.55e-2 80e-2 0 radius 19.55e-2
apply nstress -1e5 range y=0
apply nstress -1e5 range y=80e-2
;initial condition
initial yvel -2.5e-5 range y=80e-2
define axi_stress
str=0
p_gp=gp_head
loop while p_gp # null
if gp_ypos(p_gp)<0.1 then
str=str+gp_yfunbal(p_gp)
endif
p_gp=gp_next(p_gp)
end_loop
ax_str=4*str/(pi*(39.1e-2)^2)
end
hist n 1
hist id=2 gp ydisp 0 0 0
hist id=3 axi_stress
hist id=4 unbal
plot hist -2 vs 3
step 1500
