zone group 'weiyan' range group 'ZG_001'
zone group 'kuangshi1' range group 'ZG_002'
zone group 'kuangshi2' range group 'ZG_003'
zone group 'kuangshi3' range group 'ZG_004'
zone group 'kuangshi4' range group 'ZG_005'
zone group 'kuangshi5' range group 'ZG_006'
zone group 'kuangshi6' range group 'ZG_007'
zone group 'kuangshi7' range group 'ZG_008'
zone group 'kuangshi8' range group 'ZG_009'
zone group 'kuangshi9' range group 'ZG_010'
zone group 'kuangshi10' range group 'ZG_011'
zone group 'kuangshi11' range group 'ZG_012'
zone group 'kuangshi12' range group 'ZG_013'
zone group 'kuangshi13' range group 'ZG_014'
zone group 'kuangshi14' range group 'ZG_015'
zone group 'ZG_001' remove
zone group 'ZG_002' remove
zone group 'ZG_003' remove
zone group 'ZG_004' remove
zone group 'ZG_005' remove
zone group 'ZG_006' remove
zone group 'ZG_007' remove
zone group 'ZG_008' remove
zone group 'ZG_009' remove
zone group 'ZG_010' remove
zone group 'ZG_011' remove
zone group 'ZG_012' remove
zone group 'ZG_013' remove
zone group 'ZG_014' remove
zone group 'ZG_015' remove
; Name model boundaries
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
zone gridpoint initialize velocity-x 0
zone gridpoint initialize velocity-y 0
zone gridpoint initialize velocity-z 0
zone gridpoint fix velocity-x range position-z -0.1 0.1
zone gridpoint fix velocity-y range position-z -0.1 0.1
zone gridpoint fix velocity-z range position-z -924.9 -925.1
zone gridpoint fix velocity-x range position-x -0.1 0.1
zone gridpoint fix velocity-x range position-x 3207.9 3208.1 ;
zone gridpoint fix velocity-y
model gravity 0 0 -9.81 ;
zone face apply stress-xx -0.1e6 grad 0 0 -0.05e6 range position-z -2000 0
zone cmodel assign mohr-coulomb
zone property young 2.3e9 poisson 0.2 density 3000 cohesion 2.3e6 friction 44 range group 'weiyan'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi1'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi2'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi3'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi4'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi5'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi6'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi7'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi8'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi9'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi10'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi11
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi12'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi13'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi14'
model solve
model save '1'
zone gridpoint initialize displacement (0,0,0)
zone gridpoint initialize velocity (0,0,0)
zone initialize state 0
model gravity 0 0 -9.81
zone face apply stress-xx -0.1e6 grad 0 0 -0.05e6 range position-z -2000 0
zone cmodel assign null range group 'kuangshi8'
model solve
model save '2'
model restore '2'
zone gridpoint initialize displacement (0,0,0)
zone gridpoint initialize velocity (0,0,0)
zone initialize state 0
model gravity 0 0 -9.81 ;
zone face apply stress-xx -0.1e6 grad 0 0 -0.05e6 range position-z -2000 0
; Assign model and properties
zone cmodel assign elastic range group 'weiyan'
zone cmodel assign elastic range group 'kuangshi1'
zone cmodel assign elastic range group 'kuangshi2'
zone cmodel assign elastic range group 'kuangshi3'
zone cmodel assign elastic range group 'kuangshi4'
zone cmodel assign elastic range group 'kuangshi5'
zone cmodel assign elastic range group 'kuangshi6'
zone cmodel assign elastic range group 'kuangshi7'
zone cmodel assign elastic range group 'kuangshi9'
zone cmodel assign elastic range group 'kuangshi10'
zone cmodel assign elastic range group 'kuangshi11'
zone cmodel assign elastic range group 'kuangshi12'
zone cmodel assign elastic range group 'kuangshi13'
zone cmodel assign elastic range group 'kuangshi14'
zone property density 3000 bulk 3e11 shear 1e11 range group 'weiyan'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi1'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi2'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi3'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi4'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi5'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi6'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi7'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi9'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi10'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi11'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi12'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi13'
zone property density 2930 bulk 3e11 shear 1e11 range group 'kuangshi14'
model solve
model save '3'
model restore '3'
zone gridpoint initialize displacement (0,0,0)
zone gridpoint initialize velocity (0,0,0)
zone initialize state 0
model gravity 0 0 -9.81 ;
zone face apply stress-xx -0.1e6 grad 0 0 -0.05e6 range position-z -2000 0
zone cmodel assign mohr-coulomb range group 'weiyan'
zone cmodel assign mohr-coulomb range group 'kuangshi1'
zone cmodel assign mohr-coulomb range group 'kuangshi2'
zone cmodel assign mohr-coulomb range group 'kuangshi3'
zone cmodel assign mohr-coulomb range group 'kuangshi4'
zone cmodel assign mohr-coulomb range group 'kuangshi5'
zone cmodel assign mohr-coulomb range group 'kuangshi6'
zone cmodel assign mohr-coulomb range group 'kuangshi7'
zone cmodel assign mohr-coulomb range group 'kuangshi9'
zone cmodel assign mohr-coulomb range group 'kuangshi10'
zone cmodel assign mohr-coulomb range group 'kuangshi11'
zone cmodel assign mohr-coulomb range group 'kuangshi12'
zone cmodel assign mohr-coulomb range group 'kuangshi13'
zone cmodel assign mohr-coulomb range group 'kuangshi14'
zone property young 2.3e9 poisson 0.2 density 3000 cohesion 2.3e6 friction 44 range group 'weiyan'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi1'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi2'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi3'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi4'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi5'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi6'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi7'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi9'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi10'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi11
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi12'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi13'
zone property young 2e9 poisson 0.25 density 2930 cohesion 2e6 friction 42 range group 'kuangshi14'
model solve
model save '4'
