flac使用命令流-源码.zip命令流-源码_.zip

model new zone create brick point 0 0 0 0 point 1 2 0 0 point 2 0 0.5 0 point 3 0 0 3 size 3 1 3 zone create brick point 0 2 0 0 point 1 20 0 0 point 2 2 0.5 0 point 3 2 0 3 ... size 17 1 3 ratio 1.03 1 1 zone create brick point 0 2 0 3 point 1 20 0 3 point 2 2 0.5 3 point 3 12 0 13 point 4 20 0.5 3 ... point 5 12 0.5 13 point 6 20 0 13 point 7 20 0.5 13 ... size 17 1 17 ratio 1.03 1 1 fish define SSR ait1=0.02 k11=1.0 k12=2.0 ks=(k11+k12)/2 loop while (k12-k11)>ait1 coh1=12380/ks fri1=(math.atan((math.tan(20*math.pi/180))/ks))*180/math.pi dila1=20.0 ten1=1e6 grav0=-10 dens1=2000 K1=1e8 G1=3e7 command zone cmodel assign null zone cmodel assign elastic zone property bulk 1e10 shear 3e9 density @dens1 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 -0.1 0.1 zone gridpoint fix velocity-x range position-x 19.9 20.1 zone gridpoint fix velocity-x range position-x -0.1 0.1 zone gridpoint fix velocity-y model gravity 0 0 @grav0 model solve 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 cmodel assign mohr-coulomb zone property bulk @K1 shear @G1 density @dens1 cohesion @coh1 friction @fri1 dilation @dila1 tension @ten1 model solve mechanical ratio 9.8e-6 model solve cycles 10000 endcommand if zone.mech.ratio<1.0e-5 k11=ks k12=k12 else k12=ks k11=k11 endif ks=(k11+k12)/2 endloop fosfile0='_fos'+'.sav' command model save @fosfile0 endcommand end @SSR list @ks