新手求助，不收敛问题。。

struggle222

我是新手，最近在做杆件的拉拔试验模拟，但是做出的模型始终不收敛，是参数设置问题？我已经把一些参数设的挺大了，还是边界约束问题？请前辈帮忙看看，下面是命令流。还想找位前辈学习学习，我的QQ：24812193。。先谢谢大家了。。
new
gen zone radcyl p0 0.15 0.075 0.075 p1 0.15 0.15 0.075 p2 0 0.075 0.075 p3 0.15 0.075 0.15 &
dim 0.0125 0.0125 0.0125 0.0125 size 4 12 8 12 ratio 1 1 1 1.2 group 试块
interface 1 face range cyl end1 (0.15 0.075 0.075) end2 (0 0.075 0.075) rad 0.0126 &
cyl end1 (0.15 0.075 0.075) end2 (0 0.075 0.075) rad 0.0124 not
gen zone cyl p0 0.7 0.075 0.075 p1 0.7 0.0875 0.075 p2 0.2 0.075 0.075 p3 0.7 0.075 0.0875 size 4 40 8
group 杆件 range group 试块 not
gen zone refl norm 0 0 -1 origin 0.15 0.075 0.075
gen zone refl norm 0 -1 0 origin 0.15 0.075 0.075

plot blo gro
pause
model mohr range group 试块
pro bulk 14.8e9 shear 11.1e9 coh 1e6 ten 10e3 range group 试块
model elastic range group 杆件
pro bulk 16.7e9 shear 15.2e9 range group 杆件
ini x add -0.2 range group 杆件
fix y range y -0.001 0.001
fix y range y 0.149 0.151
fix z range z -0.001 0.001
fix z range z 0.149 0.151
set grav 0 0 -9.81
ini den 2000 range group 试块
ini den 2160 range group 杆件
interface 1 prop kn 120e6 ks 120e6 coh 200e3 fri 15
hist unbal
apply nstr=40e3 range x 0.049 0.051
step 5000
save tt.sav
plot set plane ori 0.075 0.075 0.075 norm 0 0 -1
plot con sx plane

struggle222

大侠。。快来。。。

struggle222

额。。。这。。。那我自己研究吧，想问怎么显示轴向应力和剪切应力。。