;Testing the B value of Poliakov et al. 1994 in Fractals 2. I have equalled lambda and Shear ;modulus, G, by making Bulk=1.667*Shear modulus (K=1.667G), this implies a Poisson's ratio of 0.25 and ;means that only one (lambda or shear) modulus is necessary to define the system. ;this is different from Polidry and Poliakov in which we are now doing a pure shear and not ;simple shear exps new tit FLAC Modelling Input File ; config extra 4 grid 100 100 gen 0.00e0,0.0e0,0.00e0,1.00e3 1.00e3,1.00e3 1.00e3,0.0e3 i=1,101 j=1,101 ; *************************************** model mohr *************************************** ; Properties for each sub-regions: ; Basalt pro dens= 2720 bulk=8.33335e9 shear=5.e9 i=1,101 j=1,101 ; bulk =1.66667 Shear pro coh=3.00e7 tens=1.50e7 i=1,101 j=1,101 pro fric=30 dil= 0 i=1,101 j=1,101 ; this sets stress everywhere to 2kb and gravity is zero ini sxx=-2e8 ini syy=-2e8 ini szz=-2e8 ini sxy=0e6 set gravity=0.0 set large fix x i=1 fix x i= 101 fix y j=1 fix y j= 101 ; Schaubs: make sure to put these commands after the fish function which computes your initial stress hist 1 unbal nstep 1 set small set mech on fl off step 200 ini xd 0 yd 0 ini xvel 0 yvel 0 set large ;initial velocity command should be here ie ini xvel 0.004 etc ; commands to calculate meanstress define meanstress while_stepping loop i (1,izones) loop j (1,jzones) ex_1(i,j)=(sxx(i,j)+syy(i,j)+szz(i,j))/3.0 end_loop end_loop end meanstress ini xvel 0.002 var=-0.004 ,0 i=1, 101 j=1, 101 ini yvel -0.002 Var = 0, 0.004 i=1, 101 j=1, 101 ; step 5000 save Polipure1a.sav step 5000 save Polipure1b.sav ; shear strain= 0.2 step 5000 save Polipure1c.sav ; shear strain= 0.3 step 5000 save Polipure1d.sav ; shear strain= 0.4 step 5000 save Polipure1e.sav ; shear strain= 0.5 step 5000 save Polipure1f.sav ; shear strain= 0.5 return stop