model restore “model_ns.f3sav”
model largestrain on
;add interface
zone cmodel assign null range position (240,100,70) (540,300,72)
zone interface “1” create by-face range position-z 69.99 70.01 position-x 245 540
zone interface “1” node property stiffness-normal 2e10 stiffness-shear 2e10 friction 30 tension 0.5e10 cohesion 0.5e10
zone cmodel assign mohr-coulomb range position (240,100,70) (540,300,72)
zone property bulk 3.858388e9 shear 1.314774e9 cohesion 2.47e6 tension 1.68e6 friction 30 density 1400 range position (240,100,70) (540,300,72)
;assign prop
zone cmodel assign mohr-coulomb
zone property bulk 4.094758e9 shear 2.433307e9 tension 2.56e6 cohesion 3.2e6 friction 33 density 2550 range group “ny”
zone property bulk 7.50879e9 shear 4.227039e9 tension 3.015e6 cohesion 4.8e6 friction 36 density 2690 range group “xsy”
zone property bulk 5.688996e9 shear 4.464116e9 tension 3.05e6 cohesion 4.3e6 friction 38 density 2730 range group “fsy”
zone property bulk 7.50879e9 shear 4.227039e9 tension 3.015e6 cohesion 4.8e6 friction 36 density 2690 range group “zlsy”
zone property bulk 3.858388e9 shear 1.314774e9 tension 1.68e6 cohesion 2.47e6 friction 30 density 1400 range group “mc”
zone property bulk 2.990637e9 shear 1.942822e9 tension 2.4e6 cohesion 2.45e6 friction 34 density 2660 range group “lzny”
zone property bulk 2.5e8 shear 1.69e8 tension 0.48e5 cohesion 0.84e6 friction 26 density 2000 range group “dcd”
zone face apply velocity-normal 0 range group ‘East’ or ‘West’
zone face apply velocity-normal 0 range group ‘North’ or ‘South’
zone face apply velocity-normal 0 range group ‘Bottom’
zone face apply stress-normal -9e6 range group “Top”
model gravity 9.8
zone initialize-stresses ratio 1.2
model solve
model save “ini.sav”