res model.sav
fix x range x -0.1 0.1
fix x range x 99.1 100.1
fix z range z -0.1 0.1
fix z range z 49.9 50.1
fix x y z range y -0.1 0.1
model mohr
;
set grav 0 -9.81 0
;
pro den 2500 bu 1.19e10 she 1.1e10
pro coh 2.72e5 fric 44 ten 2e5
ini dens 2650
;
pl set plane dd 0 dip 0 ori 0 0 25
plot set rotation 90 0 0 or 50 50 25
pl add axe
pl add bl grou plane
pl add dis plane
pl show
;
def moviedis_plot
loop ii (1,50)
command
movie snap
step 20
end_command
endloop
end
def moviedis_plot1
loop ii (1,100)
command
movie snap
step 100
end_command
endloop
end
;
set movie avi size 800 600 frameperiod 150 file kwmovie
movie start
;
moviedis_plot
;
;
ini xdisp 0 ydisp 0 zdisp 0
;
delete range group 2
moviedis_plot
;
;solve
;
;
config creep
model cpow
;
; Properties and stresses in Pascal units (not MPa)
model cpow
pro den 2500 bu 1.19e10
prop shear 1.1e10
prop a_1=1e-23 n_1=3
pro coh 2.72e5 fric 44 ten 2e5
;
set creep on
set creep dt 1.0e-4
set creep mindt 1.0e-5
set creep maxdt 1.0e-3
set creep dt auto on
;
moviedis_plot1
;
movie finish
;