采用FLAC_2D模拟三轴压缩试验的代码如下：
n
grid 25 50
gen 0 0 0 0.1 0.05 0.1 0.05 0
m m
prop dens 2735 shear 27.57e9 bulk 49.57e9 cohes 16.23e6 tens 3.29e6
prop fri 39.51
fix y j=1
fix y j=51
app sxx -6e7 i=1 j=1,51
app sxx -6e7 i=26 j=1,51
solve
ini yvel=-1.67e-7 j=51
def load
sum=0
loop i(1,igp)
sum=sum+yforce(i,jgp)
endloop
load=sum/(x(igp,jgp)-x(1,jgp))
end
def vstrain
vst=-ydisp(1,jgp)/y(1,jgp)
if vst<0 then
vst=0
endif
vstrain=vst
end
his unb
his load
his vstrain
step 4700
plot hold his 2 vs 3
为什么出来的应力应变曲线呈现阶梯状，并不平滑。
