作为FLAC新手,命令流出了问题不知道怎么改,特来此地请教大神,谢谢~
命令流如下:
new
res 041.sav
pl show
;******************************定义材料属性******************************
mo mohr
group 2 range x -2.4 0 y -2 2 z -0.5 0
group 3 range x -2.4 0 y -2 2 z -1 -0.5
group 4 range x -2.4 0 y -2 2 z -1.5 -1
group 5 range x -2.4 0 y -2 2 z -2 -1.5
group 6 range x -2.4 0 y -2 2 z -2.5 -2
group 7 range x -2.4 0 y -2 2 z -3 -2.5
group 8 range x -2.4 0 y -2 2 z -3.5 -3
group 9 range x -2.4 0 y -2 2 z -4 -3.5
group 10 range x -2.4 0 y -2 2 z -4.5 -4
group 11 range x -2.4 0 y -2 2 z -5 -4.5
group 12 range x -2.4 0 y -2 2 z -5.5 -5
group 13 range x -2.4 0 y -2 2 z -6 -5.5
group 14 range x -2.4 0 y -2 2 z -6.5 -6
group 15 range x -2.4 0 y -2 2 z -7 -6.5
group 16 range x -2.4 0 y -2 2 z -7.5 -7
group 17 range x -2.4 0 y -2 2 z -8 -7.5
group 18 range x -2.4 0 y -2 2 z -8.5 -8
group 19 range x -2.4 0 y -2 2 z -9 -8.5
group 20 range x -2.4 0 y -2 2 z -9.5 -9
group 21 range x -2.4 0 y -2 2 z -10 -9.5
group 22 range x -2.4 0 y -2 2 z -10.5 -10.5
group 23 range x -2.4 0 y -2 2 z -11 -10.5
group 24 range x -2.4 0 y -2 2 z -11.5 -11.5
group 25 range x -2.4 0 y -2 2 z -12 -11.5
group 26 range x -2.4 0 y -2 2 z -12.5 -12.5
group 27 range x -2.4 0 y -2 2 z -13 -12.5
group 28 range x -2.4 0 y -2 2 z -13.5 -13.5
group 29 range x -2.4 0 y -2 2 z -14 -13.5
group 30 range x -2.4 0 y -2 2 z -14.5 -14.5
group 31 range x -2.4 0 y -2 2 z -15 -14.5
initial density 2450 range group 1 ;中风化花岗岩
prop bulk 2800e6 shear 1680.0e6 fric 35 coh 6e5 tension 2.2e5 range group 1
initial density 2450 range x -2.4 0 y -2 2 z -15 0 ;中风化花岗岩
prop bulk 2800e6 shear 1680.0e6 fric 35 coh 6e5 tension 2.2e5 range x -2.4 0 y -2 2 z -15 0
set gravity 0 0 -9.8 ;加载重力
apply nstress -1.75e7 range z 14.9 15.1 ;地面超载(6)
ini sxx -35.48820e4 grad 0 0 2.401e4 range y -25 25
ini syy -11.824473e4 grad 0 0 0.800e4 range y -25 25 ;中风化(44.612)
ini szz -11.824473e4 grad 0 0 2.401e4 range y -25 25
;******************************定义边界条件******************************
fix y range y -24.9 -25.1
fix z range z -15.1 -14.9
fix z range z 14.9 15.1
fix x range x -0.1 0.1
fix x range x 24.9 25.1
;******************************监控变量设置******************************
h nstep=5
h unb
h gp zdis 0,2,0 ;2
h gp zdis 0,2,-0.5 ;3
h gp zdis 0,2,-1
h gp zdis 0,2,-1.5
h gp zdis 0,2,-2
h gp zdis 0,2,-2.5
h gp zdis 0,2,-3
h gp zdis 0,2,-3.5
h gp zdis 0,2,-4
h gp zdis 0,2,-4.5
h gp zdis 0,2,-5
h gp zdis 0,2,-5.5
h gp zdis 0,2,-6
h gp zdis 0,2,-6.5
h gp zdis 0,2,-7
h gp zdis 0,2,-7.5
h gp zdis 0,2,-8
h gp zdis 0,2,-8.5
h gp zdis 0,2,-9
h gp zdis 0,2,-9.5
h gp zdis 0,2,-10
h gp zdis 0,2,-10.5
h gp zdis 0,2,-11
h gp zdis 0,2,-11.5
h gp zdis 0,2,-12
h gp zdis 0,2,-12.5
h gp zdis 0,2,-13
h gp zdis 0,2,-13.5
h gp zdis 0,2,-14
h gp zdis 0,2,-14.5
h gp zdis 0,2,-15 ;32
;******************************求解******************************
set mech force 100
set mech ratio 1e-5
solve ;solve
save 042.sav ;保存初始地应力,只在重力影响
;******************************进行开挖模拟******************************
pl bl gr
pl add axes orange
initial xdisp 0 ydisp 0 zdisp 0
initial xvel 0 yvel 0 zvel 0
solve
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 10.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 10.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 10.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 10.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 10.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 10.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 10.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 10.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 10.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 9.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 9.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 9.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 9.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 9.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 9.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 9.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 9.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 9.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 8.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 8.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 8.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 8.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 8.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 8.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 8.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 8.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 8.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 7.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 7.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 7.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 7.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 7.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 7.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 7.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 7.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 7.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 6.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 6.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 6.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 6.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 6.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 6.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 6.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 6.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 6.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 5.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 5.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 5.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 5.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 5.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 5.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 5.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 5.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 5.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 4.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 4.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 4.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 4.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 4.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 4.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 4.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 4.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 4.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 3.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 3.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 3.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 3.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 3.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 3.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 3.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 3.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 3.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 2.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 2.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 2.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 2.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 2.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 2.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 2.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 2.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 2.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 1.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 1.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 1.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 1.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 1.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 1.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 1.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 1.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 1.000 nseg=4
sel cable id=3 begin -2.400 -1.600 10.000 end -4.900 -1.600 0.000 nseg=4
sel cable id=3 begin -2.400 -0.700 10.000 end -4.900 -0.700 0.000 nseg=4
sel cable id=3 begin -2.400 0.200 10.000 end -4.900 0.200 0.000 nseg=4
sel cable id=3 begin -2.400 1.100 10.000 end -4.900 1.100 0.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -4.900 2.000 0.000 nseg=4
sel cable id=3 begin -2.400 2.000 10.000 end -2.400 4.500 0.000 nseg=4
sel cable id=3 begin -1.600 2.000 10.000 end -1.600 4.500 0.000 nseg=4
sel cable id=3 begin -0.800 2.000 10.000 end -0.800 4.500 0.000 nseg=4
sel cable id=3 begin 0.000 2.000 10.000 end 0.000 4.500 0.000 nseg=4
sel cable id=3 prop xcarea=3.8e-4 emod=200e9 ytension=1e19 gr_k=2e7 gr_per=0.132 gr_coh=2e5 gr_fric=25
model null range group 2
set small
step 5000
save suxing2.sav
model null range group 3
set small
step 5000
save suxing3.sav
model null range group 4
set small
step 5000
save suxing4.sav
model null range group 5
set small
step 5000
save suxing5.sav
model null range group 6
set small
step 5000
save suxing6.sav
model null range group 7
set small
step 5000
save suxing7.sav
model null range group 8
set small
step 5000
save suxing8.sav
model null range group 9
set small
step 5000
save suxing9.sav
model null range group 10
set small
step 5000
save suxing10.sav
model null range group 11
set small
step 5000
save suxing11.sav
model null range group 12
set small
step 5000
save suxing12.sav
model null range group 13
set small
step 5000
save suxing13.sav
model null range group 14
set small
step 5000
save suxing14.sav
model null range group 15
set small
step 5000
save suxing15.sav
model null range group 16
set small
step 5000
save suxing16.sav
model null range group 17
set small
step 5000
save suxing17.sav
model null range group 18
set small
step 5000
save suxing18.sav
model null range group 19
set small
step 5000
save suxing19.sav
model null range group 20
set small
step 5000
save suxing20.sav
model null range group 21
set small
step 5000
save suxing21.sav
model null range group 22
set small
step 5000
save suxing22.sav
model null range group 23
set small
step 5000
save suxing23.sav
model null range group 24
set small
step 5000
save suxing24.sav
model null range group 25
set small
step 5000
save suxing25.sav
model null range group 26
set small
step 5000
save suxing26.sav
model null range group 27
set small
step 5000
save suxing27.sav
model null range group 28
set small
step 5000
save suxing28.sav
model null range group 29
set small
step 5000
save suxing29.sav
model null range group 30
set small
step 5000
save suxing30.sav
model null range group 31
set small
step 5000
save suxing31.sav
model null range group 32
set small
step 5000
save suxing32.sav
save shuzhimoni,sav
set log on
set logfile hist.dat
hist write 1 &
file shuipingweiyi.txt
hist write 2 3 4 5 6 7 8 9 10 &
file shuipingweiyi.txt
hist write 11 12 13 14 15 16 17 18 19 20 &
file shuipingweiyi.txt
hist write 21 22 23 24 25 26 27 28 29 30 31 32 &
file shuipingweiyi.txt
set log off