;start at disch1d6.sav from disch_fl3, but now I removed again the 2000 steps with mech off, this is similar to the initial run with disch1 but much slower and with the dis similar to that of disch_fl4 (which has a set mech of command). ;Model to test the effects of a discharge pulse at the base of the model. Track fluid propagation tit FLAC Modelling Input File ; config extra 4 gw grid 100 100 gen 0.00e0,0.00e0 0.00e0,1.00e4 1.00e4,1.00e4 1.00e4,0.00e0 i=1,101 j=1,101 mark i=1 mark j=1 mark i = 101 mark j = 101 ; *************************************** model mohr water bulk 2e9 dens 1000 prop poro=0.3 ; to be elaborated below ;; 1 mech step = 1 fluid step; modify if necessary!!! set nmech 1 ngw 1 flow on *************************************** ; ; Properties for each sub-regions: pro dens= 2720 bulk=4.09e9 shear=2.45e9 pro coh=3.00e7 tens=1.5e7 pro fric=30 dil= 5 pro perm=1.00e-13 porosity=0.3 ; ***** ini sxx=-2.66800e8 var=0, 2.66785e8 i=1, 101 j=1,101 ini syy=-2.66800e8 var=0, 2.66785e8 i=1, 101 j=1,101 ini szz=-2.66800e8 var=0, 2.66785e8 i=1, 101 j=1,101 ini sxy=0e6 fix x i=1 fix x i= 101 fix y j=1 ; ; def hydro_pp loop i (1,igp) loop j (1,jgp) if j = jgp Then sss2 = 0 else aaa = abs(y(i, j) - y(i, jgp)) sss2 = 1000 * 9.8 * aaa end_if gpp(i, j) = sss2 end_loop end_loop end hydro_pp * ; ini pp=1e5 j=101 fix pp j=101 ini sat 1 fix sat ; set gravity=9.8 ; hist 1 unbal nstep 1 set small set mech on fl off step 1000 ini xd 0 yd 0 ini xvel 0 yvel 0 plot hold hist 1 set mech on fl on step 200 ; this can be increased if you think the unbalanced forces haven't bottomed out ini xd 0 yd 0 ini xvel 0 yvel 0 set large plot hold hist 1 ;initial velocity command should be here ie ini xvel 0.004 etc; ; ini xvel 0.0002 var=-0.0004 ,0 i=1, 101 j=1, 101 step 50000 save disch1a.sav ; 0.1 % deformation step 50000 save disch1a1.sav ; 0.1 % deformation step 50000 save disch1a2.sav ; 0.1 % deformation step 50000 save disch1a3.sav ; 0.1 % deformation step 50000 save disch1a4.sav ; 0.1 % deformation step 250000 save disch1b.sav ; 0.2 % deformation step 250000 save disch1c.sav ; 0.3 % deformation step 250000 save disch1d.sav ; 0.4 % deformation step 50000 save disch1d1.sav ; 0.1 % deformation step 50000 save disch1d2.sav ; 0.1 % deformation step 50000 save disch1d3.sav ; 0.1 % deformation step 50000 save disch1d4.sav ; 0.1 % deformation step 150000 save disch1d5.sav step 150000 save disch1d6.sav apply dis= 4e-9 i= 1 101 j= 1 set mech on set nmech 1 ngw 1 flow on step 1000 save disch1e1a.sav ; 0.5 % deformation step 1000 save disch1e1b.sav ; 0.5 % deformation step 1000 save disch1e1c.sav ; 0.5 % deformation step 1000 save disch1e1d.sav ; 0.5 % deformation step 1000 save disch1e1f.sav ; 0.5 % deformation step 1000 save disch1e1g.sav ; 0.5 % deformation step 1000 save disch1e1h.sav ; 0.5 % deformation step 1000 save disch1e1i.sav ; 0.5 % deformation step 1000 save disch1e1j.sav ; 0.5 % deformation step 1000 save disch1e1.sav stop step 15000 save disch1e2.sav ; 0.5 % deformation step 15000 save disch1e3.sav ; 0.5 % deformation step 15000 save disch1e4.sav ; 0.5 % deformation step 15000 save disch1e5.sav ; 0.5 % deformation step 15000 save disch1e6.sav ; 0.5 % deformation step 15000 save disch1e7.sav ; 0.5 % deformation step 15000 save disch1e8.sav ; 0.5 % deformation step 15000 save disch1e9.sav ; 0.5 % deformation step 15000 save disch1e10.sav ; 0.5 % deformation step 15000 save disch1e11.sav ; 0.5 % deformation apply dis= 0.0 i= 1 101 j= 1 step 15000 save disch1f.sav ; 0.5 % deformation step 15000 save disch1g.sav ; 0.5 % deformation step 15000 save disch1h.sav ; 0.5 % deformation step 15000 save disch1i.sav ; 0.5 % deformation step 15000 save disch1j.sav ; 0.5 % deformation step 15000 save disch1g.sav ; 0.5 % deformation step 15000 save disch1k.sav ; 0.5 % deformation step 15000 save disch1l.sav ; 0.5 % deformation step 15000 save disch1m.sav ; 0.5 % deformation step 50000 save disch1n.sav ; 0.6 % deformation step 50000 save disch1o.sav ; 0.7 % deformation step 50000 save disch1p.sav ; 0.6 % deformation step 50000 save disch1q.sav ; 0.7 % deformation step 250000 save disch1r.sav ; 0.8 % deformation step 250000 save disch1s.sav ; 0.8 % deformation return stop