new
config creep
gen zone brick size 4 4 10
title
test programe
def setup
bu = 20.
g_k = 10.
g_m = 15.
vis_k = 10.
vis_m = 30.
P = -10.
cons = 1./(3.*bu) + 1./g_m
alpha = g_k/vis_k
end
setup
model burger
prop coh 1e20 tension 1e20
prop bulk=bu
prop mshear=g_m mvisc=vis_m
prop kshear=g_k kvisc=vis_k
fix z range z 0 0
app szz=p range z 10 10
def strn_m
strn_m = P*(crtime/vis_m + cons)/3.
strn_k = P*((1.-exp(-alpha*crtime))/g_k + cons)/3.
strn_b = P*(crtime/vis_m + (1.-exp(-alpha*crtime))/g_k + cons )/3.
end
strn_m
hist n=10
hist crtime
hist unbal
hist gp zdisp 0 0 10
hist strn_b
hist dt
; --- initial "instantaneous" equilibrum ---
solve
; -----reset velocities----
ini xv 0 yv 0 zv 0
; --- creep test ---
set creep mindt=1e-5 maxdt=0.1
solve age 10
save creep.sav
title
z displacment vs creeptime
plot his -3 vs 1 ymin 0
pause key
print creep
ret
计算结果中，位移蠕变时间曲线不包含加速蠕变阶段，应该如何调整？命令流哪儿还有问题？请高手多多指点！谢谢！
[ 此贴被joidon在2008-05-19 13:41重新编辑 ]
