蠕变模型应力加载不上

xiaokang_2012

刚接触蠕变，自己建了一个单元小模型，计算过程中监测应力应变，但是应力出现跳动，不知道是什么原因，求大神指教，下面是程序：
new
gen zone b p0  0 0 0 p1 1 0 0 p2 0 1 0 p3 0 0 1 size 1 1 1
plot show
plot add surface yellow
save danyuan.sav
config creep;启动蠕变模块
;材料参数
model cvisc
prop bulk 11.9e9 shear 11e10 coh 2.72e6 fric 34 ten 0.8e6 dilation=10  mvis 1e13 kvis 1e13
;边界条件
fix z range z -0.01 0.01
;加载条件
app szz -10e6 range z 0.99 1.01

set creep mindt=0.1 maxdt=2000
set creep on;开启蠕变计算
set cr dt 1;设置时间步长


;指针
def jiance
  p_gp=gp_near(0,0,1)
  p_z=z_near(0.5,0.5,0.5)
end
jiance


def jisuan
  loop i (1, 3000)
    a=gp_zdisp(p_gp)
    b=z_szz(p_z)
    command
      solve age 1
      hist  a
      hist  b
      hist  crtime
      hist  unbal
      plot show
      plot clear
      plot add hist -2 vs 3
    endcommand
  endloop
end
jisuan

xiaokang_2012

监测的应力应变关系

xiaokang_2012

求高手指点一二

lf_1989

看不出有什么问题额