求助：高手帮我看看蠕变程序问题

飞鱼游鹰

def ini_cons
  c_bu = 1e5    ;弹性体积模量
  c_sh = 1e5    ;弹性剪切模量
  c_vi = 2e5    ;粘性系数
  
  c_ksh=2e5
  c_pr = 1e5    ;作用力
  a1  = 1/(3*c_sh) +1/(9*c_bu)
  c_a=1/3*(c_ksh)+1/(9*c_bu)
  c_b=9*c_ksh*c_bu/(3*c_bu+c_ksh)
  c_c=c_b/(c_vi)
  
end
ini_cons
conf  creep ext=1
g 20 20
mod cvisc
fix x i=1
fix x i=21
fix y j=1
gen  0 0 0 0.2 0.05 0.2 0.05 0   
prop density 1060
prop bu c_bu sh c_sh viscosity c_vi
prop  k_sh c_ksh
prop coh 2.72e5 tension 2e6   ;(没有给出粘性系数？）
def confine                    ;定义围压
  confine=-con
  app_x=con
end
set con=1e5
def vstrain                      ;定义轴向应变
  vst = -ydisp(1,jgp) / y(1,jgp)
  if vst < 0.0 then
    vst = 0.0
  end_if
  vstrain = vst
end
def disp_mag
loop i (1,igp)
  loop j (1,jgp)
    ex_1(i,j) = sqrt(xdisp(i,j)^2 + ydisp(i,j)^2)
  end_loop
end_loop
end                    ; (定义最大位移）
ini sxx confine syy confine  szz confine
apply pressure c_pr   j=21     ;  施加上边界应力(400 kPa)
def ana_eyy
  ana_eyy = -c_pr * (a1 +c_a* (1-exp(-c_c* crtime)))
  
  ana_syy = -c_pr
end
; --- elastic equilibrium ---
step 1000
his sxx  i=1 j=21
his syy  i=1 j=21
his ana_syy
his ydis i=1 j=21
his ana_eyy
his crtime
hist vstrain
set dt=1.e-3
step 1000
set sratio 0.0
set fobl=1.e-8 lmul=2 umul=1 max=0.1 min=1.e-3
set crdt=100 ;(此时进行初始平衡 ，还未进行动力计算\也未进
step 38880


plot hold hist 5 line skip 10 vs 6
plot hold history 4 line skip 10  vs 6
plot hold hist 7 line  skip 10 vs 6

我编的这个蠕变 程序 是在单轴衡载下的蠕变 我希望它能按我导入的应变关系来发展蠕变 就是广义开尔文模型 可是结果不是我想象的那样   应变结果跳来跳去  总的位移也不是我想的那样 希望高手给我一些建议和意见  不胜感激^_^

swordchang

学习了，虽然没有用过蠕变