title (I7&8$Zl
Oedometric test phi=10 psi=10 DO1 JPeIi
round 0.01 xS~yH [k
block 0 0 0 1 1 1 1 0 !#c[~erNZ
gen quad 2 mI7rx`4H
prop mat 1 dens 1 lbKv
;change cons 3 =nvAOvP{?
;prop mat 1 bu 200 sh 200 co 1 fric 10dil 10 Tw`c6^%^y
zone m m RD6h=n4B
zone bulk 200 sh 200 co 1 fric 10 Q"OV>kl k
zone dil 10 g<2lPH
;zone dil 0 kj{rk^x
step 0 ?Rt 1CDu
def d_sigy T Oco({/_/
c_k = bulk_mod x0u?*5-t
c_g = shear_mod T$n>7X-r
e1 = c_k+4.0*c_g/3.0 of+phMev
e2 = c_k-2.0*c_g/3.0 wWJQ ~i?
sf = friction * degrad I|F~HUzA"
nf = sin(sf) %Rd~|$@>x
nf = (1.0+nf)/(1.0-nf) Jcalf{W6
sp = dilation*degrad kJurUDo
np = sin(sp) J-, H6u
np = (1.0+np)/(1.0-np) {
OxAY_
r1 = (e1-e2*nf)/((e1+e2)*nf*np-2.0*e2*(nf+np)+2.0*e1) ''(fH$pY
v_dt = -1.0e-5 jMf 7J
velv =v_dt/tdel v ?Yd LR
dsigy = v_dt*(e1+2.0*r1*(e2*np-e1)) D%jD 8 p
step1 = -2.0*cohesion*sqrt(nf)/((e1-e2*nf)*v_dt) e7XsyL'|p
end hi {2h04
def esigy r*i$+ Z
while_stepping _H4$$
if step < step1 then kMl @v`
a_sy = a_sy+e1*v_dt PTzp;.
else 6+Wr6'kuH
a_sy = a_sy+dsigy 'YZI>V*
end_if .*EOVo9S
n_sy = z_syy(z_pnt) vZ[ $H
end ~mXZfG/D
set bulk_mod 200 shear_mod 200 xm}q6>jRV
set cohesion 1 friction 10 dilation 10 l:zU_J6
d_sigy vbRrk($`
def z_pnt nE]~E xr
ib = block_head (>rS
_#^
z_pnt = b_zone(ib) x2j /8]'o
end ZL-uwI!`D
z_pnt (o x4K{
bound xvel 0 yvel 0 vh|Tb5W<
bound yvel velv range 0 1 0.9 1.1
NVO9XK
hist n 50 5W[3_P+
hist ydisp 0 1 Jt-X mGULB
hist n_sy z"6ZDC6
hist a_sy [GR]!\!%~
label hist 1 (#j2P0B
Vertical displacement W>` g;[ W
label hist 2 Gut J_2f^9
Vertical stress_udec e 8d5(e
label hist 3 P"R97#C
vertical stress_analysis 9C557$nS^
step 865 _.d}lK3$2
程序读入后，系统提示“warning：block 202 has no contacts"
