【复刻】FLAC3D 6.0 真空预压

model new;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;                                                                   ;;      Conference of Hohai University graduate students             ;;          Simple example of vacuum preloading                      ;;                  lakewater@126.com                                ;;                                                                   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;model title 'vacuum_preloading_by_lakewater'model configure fluidzone fluid biot on;zone create brick point 0 0 0 0 point 1 5 0 0 point 2 0 2 0 point 3 0 0 20 size 3 1 10 group 'Face' zone group 'Soft' range position-z 10 18zone group 'Clay' range position-z 0 10;zone cmodel assign mohr-coulombzone property bulk 33e6 shear 7e6 cohesion 0 friction 30 range group 'Face'zone property bulk .83e6 shear .17e6 cohesion 7.0e3 friction 5 range group 'Soft'  zone property bulk 6.66e6 shear 1.42e6 cohesion 8.0e3 friction 20 range group 'Clay';;zone initialize density 1500zone fluid cmodel assign isotropiczone gridpoint initialize fluid-tension -5e5zone initialize fluid-density 1000zone gridpoint initialize biot 4e9zone fluid property permeability 10e-13zone fluid property biot 1zone fluid property permeability 10e-10 range group 'Face';zone gridpoint fix velocity-x range position-z -0.1 0.1zone gridpoint fix velocity-y range position-z -0.1 0.1zone gridpoint fix velocity-z range position-z -0.1 0.1zone gridpoint fix velocity-y range position-y -.01 0.01zone gridpoint fix velocity-y range position-y 1.9 2.1zone gridpoint fix velocity-x range position-x -.01 .01zone gridpoint fix velocity-x range position-x 4.9 5.1;;Initial Stress  K0=1model gravity 0 0 -10zone initialize stress-zz -400e3 gradient 0 0 20e3zone initialize stress-yy -400e3 gradient 0 0 20e3zone initialize stress-xx -400e3 gradient 0 0 20e3zone gridpoint initialize pore-pressure 200e3 gradient 0 0 -10e3;model fluid active offmodel cycle 5000;zone gridpoint initialize displacement-x 0zone gridpoint initialize displacement-y 0zone gridpoint initialize displacement-z 0zone gridpoint initialize velocity-x 0zone gridpoint initialize velocity-y 0zone gridpoint initialize velocity-z 0;for fixed pp on the PVD positionfish define fixPP  loop i(1,10)   pv = -80e3 +(i-1)*10e3   zt = 20.1-(i-1)   zb = 20.1-i   command    zone gridpoint fix pore-pressure @pv range position-x -.1 .1 position-z @zb @zt   end_command  end_loopend@fixpp  
model cycle 5000;model fluid active onzone gridpoint initialize biot 4.3e7 range group 'Soft' notzone gridpoint initialize biot 9.0e6 range group 'Soft';
model mechanical substep 100model mechanical slave onmodel fluid substep 30model fluid slave off; --- histories ---fish define day  day=fluid.time.total/24/3600end@day
fish history @day zone history name '2' displacement-z position 2.5 0 20zone history name '3' displacement-z position 2.5 0 10;; --- solve to 7776000sec, 3 month ---
model solve unbalanced-maximum 15e3 time-total 7.776e6