求助！

做一个砂岩的单轴压缩100*100*200的长方体，想出一个应力应变曲线。

但是不知道为什么出来是直线。

用位移控制，不知道命令流具体怎么写。

轴向是z轴，求大神帮助。

附上命令流求助。gr1是填充的材料。gr2 3是砂岩

MO MOHR

prop bul 2.45E+10 shea 1.00E+10 ten 1.17E+06 cohe 1.72E+07 fric 27.8 dens 1860 range gr 1

prop bul 2.09E+10 shea 1.25E+10 ten 1.17E+06 cohe 2.07E+07 fric 31.4 dens 2080 range gr 2

prop bul 2.09E+10 shea 1.25E+10 ten 1.17E+06 cohe 2.07E+07 fric 31.4 dens 2080 range gr 3

fix z range z -0.1 0.1

fix z range z 199.9 200.1

ini zvel 2E-6 range z -0.1 0.1

ini zvel -2E-6 range z 199.8 200.1

hist gp xdisp 50 40 100

hist gp xdisp 50 45 100

hist gp xdisp 50 50 100

hist gp xdisp 50 55 100

hist gp xdisp 50 60 100

hist gp ydisp 50 40 100

hist gp ydisp 50 45 100

hist gp ydisp 50 50 100

hist gp ydisp 50 55 100

hist gp ydisp 50 60 100

hist gp zdisp 50 40 100

hist gp zdisp 50 45 100

hist gp zdisp 50 50 100

hist gp zdisp 50 55 100

hist gp zdisp 50 60 100

hist zone sxx 50 40 100

hist zone sxx 50 45 100

hist zone sxx 50 50 100

hist zone sxx 50 55 100

hist zone sxx 50 60 100

hist zone syy 50 40 100

hist zone syy 50 45 100

hist zone syy 50 50 100

hist zone syy 50 55 100

hist zone syy 50 60 100

hist zone szz 50 40 100

hist zone szz 50 45 100

hist zone szz 50 50 100

hist zone szz 50 55 100

hist zone szz 50 60 100

pl con zd out on

solve

save 145635.sav