http://s15/mw690/002B1i7lgy6R7UvGbKe3e&690
注意再次打开用fish语言生成的.sav的时候把.sav中多余的0和空格去掉就可以用flac3d rest 成功了，不然可能出现不能打开的情况
**********************************下面的修改（添加int即可改掉上面多余的空格）*************************

;new
;restore zuomianwancheng30-120chongtian.sav
set small
plot con szz gr
plot set plane nor 1 0 0 ori 192.42 0 0
plot con szz plane
plot set rotation ( 0.00, 0.00,90.00)
plot set magnification  1.20e+000
;plot add int red
plot add axes
;采全高
;采右面巷道
   model null range group 8#coal_up x 168.9605 173.9605 y 30.00 120.00;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_mid x 168.9605 173.9605 y 30.00 120.00;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_down x 168.9605 173.9605 y 30.00 120.00;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了

   model null range group 8#coal_up x 210.8940 215.8940   y 30.00 120.00;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_mid x 210.8940 215.8940  y 30.00 120.00;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_down x 210.8940 215.8940 y 30.00 120.00;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
step 500
save 右面大巷.sav
;------------------------------------
;不加interface了，用ten coh 和fri 都很小的mohr模型代替interface传递力
;加intface在全高煤层底面（加在了19-1上表面）
;gen separate  19-1;
;------------------excu------------------------------
def kaiwa_youmian_quangao
loop k(0,8)
  y_beg=10*float(k)+30
  y_end=10*float(k)+40
   y_beg2=y_beg-10;落后10m垮落
   y_end2=y_end-10
name='东欢矿右面全高从30m开挖到'+string(int(y_end))+'.sav'
command

def ini_sxx_syy 
whilestepping               
    _k=0.5                
   p_z=zone_head                
   loop while p_z # null                
   z_sxx(p_z)=_k*z_szz(p_z)                
   z_syy(p_z)=_k*z_szz(p_z)                
    p_z=z_next(p_z)                
  endloop                
 end                
ini_sxx_syy

   model null range group 8#coal_up x 173.9605 210.8940 y y_beg y_end;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_mid x 173.9605 210.8940 y y_beg y_end;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_down x 173.9605 210.8940 y y_beg y_end;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
def congtian
   if k#0 then
     command
    model  mohr range group 8#coal_up x 173.9605 210.8940 y y_beg2  y_end2;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    model  mohr range group 8#coal_mid x 173.9605 210.8940 y y_beg2  y_end2;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
     model  mohr range group 8#coal_down x 173.9605 210.8940 y y_beg2  y_end2;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_up x 173.9605 210.8940 y y_beg2  y_end2; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_mid x 173.9605 210.8940 y y_beg2  y_end2; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_down x 173.9605 210.8940 y y_beg2  y_end2; 20-8#煤   
     end_command  
    end_if
end
congtian
   step 1000

   save name
     end_command
end_loop
end
kaiwa_youmian_quangao
;------------------充填最后一步-------------
     model  mohr range group 8#coal_up  x 173.9605 210.8940  y  110 120;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    model  mohr range group 8#coal_mid  x 173.9605 210.8940  y  110 120;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
     model  mohr range group 8#coal_down  x 173.9605 210.8940  y  110 120;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_up x 173.9605 210.8940  y  110 120; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_mid x 173.9605 210.8940  y  110 120; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_down x 173.9605 210.8940  y  110 120; 20-8#煤   
step 4500
save youmianwancheng30-120chongtian_youhang.sav
call 3-中面开挖全高_30m煤柱有两边巷道.txt

whilestepping               
    _k=0.5                
   p_z=zone_head                
   loop while p_z # null                
   z_sxx(p_z)=_k*z_szz(p_z)                
   z_syy(p_z)=_k*z_szz(p_z)                
    p_z=z_next(p_z)                
  endloop                
 end                
ini_sxx_syy

   model null range group 8#coal_up x 33.16 70.09 y y_beg y_end;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_mid x 33.16 70.09 y y_beg y_end;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
   model null range group 8#coal_down x 33.16 70.09 y y_beg y_end;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
def congtian
   if k#0 then
     command
    model  mohr range group 8#coal_up x 33.16 70.09 y y_beg2  y_end2;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    model  mohr range group 8#coal_mid x 33.16 70.09 y y_beg2  y_end2;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
     model  mohr range group 8#coal_down x 33.16 70.09 y y_beg2  y_end2;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_up x 33.16 70.09 y y_beg2  y_end2; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_mid x 33.16 70.09 y y_beg2  y_end2; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_down x 33.16 70.09 y y_beg2  y_end2; 20-8#煤   
     end_command  
    end_if
end
congtian
   step 1000
   ;solve
   save name
     end_command
end_loop
end
kaiwa_zuomian_quangao
     model  mohr range group 8#coal_up x 33.16 70.09  y  110 120;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    model  mohr range group 8#coal_mid x 33.16 70.09  y  110 120;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
     model  mohr range group 8#coal_down x 33.16 70.09  y  110 120;左工作面50m宽开挖全高，因为倾斜30m边界就不到30了
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_up x 33.16 70.09  y  110 120; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_mid x 33.16 70.09  y  110 120; 20-8#煤
    prop bu 1.00E+09 sh 0.50E+09 coh 2.00E+02 fri 10 ten 1.60E+02 range group 8#coal_down x 33.16 70.09  y  110 120; 20-8#煤   
step 4500
save zuomianwancheng30-120chongtian.sav
call 2-右面开挖充填全高有两边巷道.txt