有大佬知道怎么进行模拟直剪实验或者剪切蠕变实验吗？为什么我模拟的时候老是不对劲。
new project
set fish autocreate off
title
sheer test --- cvisc models
config creep
gen zone brick size 8 8 4 &
p0(0,0,0) p1(8,0,0) p2(0,8,0) p3(0,0,4) &
p4(8,8,0) p5(0,8,4) p6(8,0,4) p7(8,8,4) &
group Base
gen zone brick size 8 8 4 &
p0(0,0,8) p1(8,0,8) p2(0,8,8) p3(0,0,12) &
p4(8,8,8) p5(0,8,12) p6(8,0,12) p7(8,8,12) &
group Top
interface 1 face range plane norm 0 0 1 ori 4 4 4 dist 0.1
save int0
plot add zone addlabel "Zone" white
plot add interface line width 2 style solid clear addlabel "1" &
orange alias "1"
ini zpos add -4 range group Top
ini dens=1900 range
def setup
global bulk = 2.5e10
global g_k = 1.6e10
global vis_k = 2.746
global P = -2e6
;global cons = 1./(3.*bulk) + 1./g_k
global alpha = g_k/vis_k
end
@setup
interface 1 prop kn 2.5e10 ks 1.6e10 coh 2e6 fric 45 ten 0
set grav 0 0 -10
model mech cvisc range
prop coh 2e6 tension 1e20 range
fix z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 7.9 8.1
fix x y range x 0 z 4 8
fix x y range x 8 z 4 8
fix x y range y 0 z 4 8
fix x y range y 8 z 4 8
apply syy @P range y 8 z 0 4
def strn_m
global strn_k = P*(1.-exp(-alpha*crtime))/vis_k +P/g_k
end
his add gp ydisp 0 0 4
his add gp ydisp 2 2 4
his add gp ydisp 4 4 4
his add crtime
hist add fish @strn_k
hist add dt
solve
ini xv 0 yv 0 zv 0
set creep dt
set creep dt auto on
set creep mindt 1e-5 maxdt 0.1
solve age 6.
save shear-test
return
算了老半天也没算完
