大佬们帮忙看一下，做出来的伪三轴应力应变曲线总是有问题，也不知道问题出在哪里，求指教。
new
gen zone cyl p0 0 0 0 p1 0.025 0 0 p2 0 0.10 0 p3 0 0 0.025 size 20 50 20
gen zone reflect norm 1 0 0
gen zone reflect norm 0 0 1
group weisanzhou
model mohr
prop bulk 2.5e10 shear 1.2e10 den 2748 range group weisanzhou
prop coh 4e6 fric 45 ten 2e6 range group weisanzhoufix x y z range y -0.001 0.001
fix x y z range y 0.099 0.101apply nstress -4e3 range cyl end1 0 0 0 end2 0 0.1 0 radius 0.024 0.026ini yvel -5e-4 range y 0.099 0.101def axi_stress
f_accum=0.0
pnt = gp_head
loop while pnt # null
if gp_ypos(pnt) <0.003 then
f_accum = f_accum + gp_yfunbal(pnt)
endif
pnt = gp_next(pnt)
endloop
axi_stress = f_accum / (3.141592653*0.025*0.025)
enddef lat_stress
g_accum=0.0
pnt = gp_head
loop while pnt # null
if gp_xpos(pnt) <0.00125 then
g_accum = g_accum + gp_xfunbal(pnt)
endif
pnt = gp_next(pnt)
endloop
lat_stress = g_accum / (3.141592653*0.025*2*0.1)
endhist n 1
hist @axi_stress
hist @lat_stress
hist gp ydis 0 0.1 0
hist gp ydis -0.025 0.1 0
hist gp xdisp -0.025 0.05 0
hist gp xdisp 0.025 0.05 0
hist unbalplot hist -1 vs -3
plot hist -1 vs -4
;plot hist -2 vs -5step 1500
solve

这是应力-位移曲线吧，并没有看到你定义应变

解决了吗楼主。。请教一下

楼主 请教一下您这个命令流