FLAC3D一个采煤工作面的模拟

set logfile coal_41
title 'sunchun workface model'
;create bound------30m
macro model_size 'size 50 50'
gen zone brick model_size 10 &
p0 0 0 0 p1 142.56 0 0 p2 0 156 0 p3 0 0 30 &
p4 142.56 156 0 p5 0 156 30 p6 142.56 0 102.546 p7 142.56 156 102.546
group bound

inter 1 face range plane dip 27 dd 270 ori 142.56 0 102.546

;create bottom_4 -----2m
gen zone brick model_size 1 &
p0 0 0 30 p1 142.56 0 102.546 p2 0 156 30 p3 0 0 32 &
p4 142.56 156 102.546 p5 0 156 32 p6 142.56 0 104.546 p7 142.56 156 104.546
group bottom_4 range group bound not
inter 2 face range plane dip 27 dd 270 ori 142.56 0 104.546


;create coal_4 ------2m
gen zone brick model_size 1 &
p0 0 0 32 p1 142.56 0 104.546 p2 0 156 32 p3 0 0 34 &
p4 142.56 156 104.546 p5 0 156 34 p6 142.56 0 106.546 p7 142.56 156 106.546
group coal_4 range group bound not group bottom_4 not
inter 3 face range plane dip 27 dd 270 ori 142.56 0 106.546


;create roof_4 -----15.5m
gen zone brick model_size 5 &
p0 0 0 34 p1 142.56 0 106.546 p2 0 156 34 p3 0 0 49.5 &
p4 142.56 156 106.546 p5 0 156 49.5 p6 142.56 0 122.046 p7 142.56 156 122.046
group roof_4 range group bound not group bottom_4 not group coal_4 not
inter 4 face range plane dip 27 dd 270 ori 142.56 0 121.546


;create roof_3----3m

gen zone brick model_size 1 &
p0 0 0 49.5 p1 142.56 0 122.046 p2 0 156 49.5 p3 0 0 52.5 &
p4 142.56 156 122.046 p5 0 156 52.5 p6 142.56 0 125.046 p7 142.56 156 125.046
group roof_3 range group bound not group bottom_4 not group coal_4 not group roof_4 not

inter 5 face range plane dip 27 dd 270 ori 142.56 0 125.046


;create coal_2-----3m
gen zone brick model_size 1 &
p0 0 0 52.5 p1 142.56 0 125.046 p2 0 156 52.5 p3 0 0 55.5 &
p4 142.56 156 125.046 p5 0 156 55.5 p6 142.56 0 128.046 p7 142.56 156 128.046
group coal_2 range group bound not group bottom_4 not group coal_4 not group roof_4 not &
group roof_3 not
inter 6 face range plane dip 27 dd 270 ori 142.56 0 128.046


;create directroof_2 ------8m
gen zone brick model_size 2 &
p0 0 0 55.5 p1 142.56 0 128.046 p2 0 156 55.5 p3 0 0 63.5 &
p4 142.56 156 128.046 p5 0 156 63.5 p6 142.56 0 136.046 p7 142.56 156 136.046
group directroof_2 range group bound not group bottom_4 not group coal_4 not group roof_4 not &
group roof_3 not group coal_2 not
inter 7 face range plane dip 27 dd 270 ori 142.56 0 136.046


;create basicroof_2-------30m
gen zone brick model_size 10 &
p0 0 0 63.5 p1 142.56 0 136.046 p2 0 156 63.5 p3 0 0 93.5 &
p4 142.56 156 136.046 p5 0 156 93.5 p6 142.56 0 166.046 p7 142.56 156 166.046
group basicroof_2 range group bound not group bottom_4 not group coal_4 not group roof_4 not &
group roof_3 not group coal_2 not group directroof_2 not
inter 8 face range plane dip 27 dd 270 ori 142.56 0 166.046

;create ubound---------30m
gen zone brick model_size 10 &
p0 0 0 93.5 p1 142.56 0 166.046 p2 0 156 93.5 p3 0 0 196.046 &
p4 142.56 156 166.046 p5 0 156 196.046 p6 142.56 0 196.046 p7 142.56 156 196.046
group ubound range group bound not group bottom_4 not group coal_4 not group roof_4 not &
group roof_3 not group coal_2 not group directroof_2 not &
group basicroof_2 not

model e
;model m

prop bu 4.86e9 shear 1.02e9


prop den 2000
;prop den 1334 range group coal_4 group coal_2
;prop den 2500 range group bound group ubound
;prop den 2591 range group bottom_4
;prop den 2534 range group roof_4 group roof_3
;prop den 2662 range group directroof_2
;prop den 2558 range group basicroof_2

inter 1 pro ks 0.6e9 kn 0.6e9 fri 15
inter 2 pro ks 0.6e9 kn 0.6e9 fri 15
inter 3 pro ks 0.6e9 kn 0.6e9 fri 15
inter 4 pro ks 0.6e9 kn 0.6e9 fri 15
inter 5 pro ks 0.6e9 kn 0.6e9 fri 15
inter 6 pro ks 0.6e9 kn 0.6e9 fri 15
inter 7 pro ks 0.6e9 kn 0.6e9 fri 15
inter 8 pro ks 0.6e9 kn 0.6e9 fri 15


fix x range x -0.1 0.1
fix x range x 142.5 142.6
fix y range y -0.1 0.1
fix y range y 155.9 156.1
fix z range z -0.1 0.1

set grav 0 0 -10

hist unbal
his gp vel 0 0 196.046
his gp zdis 0 0 196.046

solve
save coal_42.sav