new
config fluid
gen zone brick size 20 1 10
model mohr
prop bulk 8.33e6 shear 3.85e6 fric 15 coh 10e3 tens 1e10;
fix x range x -.1 .1
fix x range x 19.9 20.1
fix x y z range z -.1 .1
fix y
; --- apply load slowly ---
def ramp
ramp = min(1.0,float(step)/200.0)
end
apply nstress = -40e3 hist ramp range x -.1 3.1 z 9.9 10.1
; --- fluid flow model ---
model fl_iso
ini fmod 2e9
; --- pore pressure fixed at zero at the surface ---
fix pp 0 range z 9.9 10.1
; --- settings ---
set fl off
; --- histories ---
hist gp pp 2,.5,9
; --- test ---
step 750
save load.sav
plot set plane ori 0 0.5 0 norm 0 1 0
plot con pp plane ou on
plot add fap red plane
