想模拟一个受竖向应力的柱子的应力分布情况？？
但是结构总是平衡不了，Flac3d一算就没完了，一直都在算。
下面是我的命令行 模型我是在CAD中建好后，到ansys中划分网格，又到flac3d里计算的。
;=================设置边界条件==================
fix x y z ran z -.1 .1
;==============初始地应力的生成============
model elas
prop dens 2000 bulk 3e9 shear 1e9 ;设置物理力学参数用较大值
set gravity 0 0 -9.8
solve
ini xdisp 0 ydisp 0 zdisp 0 ;位移场清零
ini xvel 0 yvel 0 zvel 0 ;速度场清零
;=============安全系数求解===============
model mohr
prop density 2270 bulk 4.4e8 shear 3.0e7 &
coh 80000.0 tens 1.7e5 fric 31.9 ;dila 剪胀角
app szz -1e6 range z 7.5 7.9 ;顶部施加竖向应力
solve
不知道哪里出错了。求指导！！！
