title lrh6lt)
Oedometric test phi=10 psi=10 fwAN9zs
round 0.01 C)z[Blt
block 0 0 0 1 1 1 1 0 $_Q o
gen quad 2 &U xN.vl
prop mat 1 dens 1 VSZ 6;&2^
;change cons 3 2>O2#53ls0
;prop mat 1 bu 200 sh 200 co 1 fric 10dil 10 M Zw%s(lv
zone m m 6EK+] 0
zone bulk 200 sh 200 co 1 fric 10 i&>,aiH @
zone dil 10 @\!wW-:A
;zone dil 0 ".xai.trr
step 0 ^l^_ K)tw*
def d_sigy ^^?DYC
c_k = bulk_mod n/?5[O-D]
c_g = shear_mod oJ8_hk<Va8
e1 = c_k+4.0*c_g/3.0 !nj%n
e2 = c_k-2.0*c_g/3.0 0v0Y(
Mo@
sf = friction * degrad >W'SG3Hmc
nf = sin(sf) (4 {49b
nf = (1.0+nf)/(1.0-nf) U&3*c +B4
sp = dilation*degrad hDlk! #*
np = sin(sp) 9'toj%XQ
np = (1.0+np)/(1.0-np) kFM'?L&
r1 = (e1-e2*nf)/((e1+e2)*nf*np-2.0*e2*(nf+np)+2.0*e1) 31>k3IP&
v_dt = -1.0e-5 \-s) D#Y;r
velv =v_dt/tdel g@Ni!U"_c
dsigy = v_dt*(e1+2.0*r1*(e2*np-e1)) /"CKVQ
step1 = -2.0*cohesion*sqrt(nf)/((e1-e2*nf)*v_dt) Sc6wC H
end & OhKx
def esigy SC2C%.%l`
while_stepping 45MK|4\Y_
if step < step1 then f aO8
&
a_sy = a_sy+e1*v_dt "}SER C7
else &ApJ'uC
a_sy = a_sy+dsigy Xe+&/J5b
end_if <YeF?$S}
n_sy = z_syy(z_pnt) : GK{ JP
end _;B!6cRLps
set bulk_mod 200 shear_mod 200 N@Mea O
set cohesion 1 friction 10 dilation 10 N0v d>b
d_sigy iY|YEi8
def z_pnt qfSoF|
ib = block_head Xb^\{s?b
z_pnt = b_zone(ib) B E"nyTQ
end uaPB M<
z_pnt oLr"8R\d>t
bound xvel 0 yvel 0 dWqFP
bound yvel velv range 0 1 0.9 1.1 5h:SH]tn8]
hist n 50 M@'V4oUz
hist ydisp 0 1 (C-z8R
Z6
hist n_sy Jl]]nO BQ/
hist a_sy xD\Km> |i
label hist 1 doR'=@ W
Vertical displacement }uHc7gTBF7
label hist 2 TLSy+x_gX
Vertical stress_udec 5&X
label hist 3 ZHC sv]l
vertical stress_analysis @i{JqHU"
step 865 3K?0PRg
程序读入后，系统提示“warning：block 202 has no contacts"
