求助：高手帮我看看蠕变程序问题
def ini_consc_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
confcreep ext=1
g 20 20
mod cvisc
fix x i=1
fix x i=21
fix y j=1
gen0 0 0 0.2 0.05 0.2 0.05 0
prop density 1060
prop bu c_bu sh c_sh viscosity c_vi
propk_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 confineszz 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 sxxi=1 j=21
his syyi=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 10vs 6
plot hold hist 7 lineskip 10 vs 6
我编的这个蠕变 程序 是在单轴衡载下的蠕变 我希望它能按我导入的应变关系来发展蠕变 就是广义开尔文模型 可是结果不是我想象的那样 应变结果跳来跳去总的位移也不是我想的那样 希望高手给我一些建议和意见不胜感激^_^ 学习了，虽然没有用过蠕变
页:
[1]