- 积分
- -2
- 注册时间
- 2011-9-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-9-7 11:41:19
|
显示全部楼层
来自 安徽淮南
我是刚开始学习 请各位看看里面有什么错误 谢谢0 R" L ?4 o3 q* p5 O: u
new8 b' s u4 Q& @0 Z( P: u
plot set title text ""5 h- H% ?0 s5 g
set logfile mn1.log4 ?4 z# l k/ w( _& [ ` X9 ?
set log on. f# E) r' p8 P- {) C& F
plot add axes red7 ?% M2 u) c- Z
9 y4 t9 [% g2 x6 F' s
;巷道及其上下部分! d/ p7 q6 e2 T7 c: R) A, j
gen zone brick size 4 20 2 p0 -1.6 0 0 p1 1.6 0 0 p2 -1.6 100 0 p3 -1.6 0 1.3 ;巷道 0 \7 z5 Z- K% m, j; Z5 J6 a: [7 \
gen zone brick size 4 20 4 p0 -1.6 0 1.3 p1 1.6 0 1.3 p2 -1.6 100 1.3 p3 -1.6 0 6 ;巷道顶部4 s9 t! h3 a8 V5 @7 w. ^
gen zone brick size 4 20 5 p0 -1.6 0 6 p1 1.6 0 6 p2 -1.6 100 6 p3 -1.6 0 11 ;group hl老顶
0 Z/ B8 I/ d, rgen zone brick size 4 20 7 p0 -1.6 0 11 p1 1.6 0 11 p2 -1.6 100 11 p3 -1.6 0 40 ;group hs上覆岩层
3 v+ M* ?& v- k! l% jgen zone brick size 4 20 2 p0 -1.6 0 -1 p1 1.6 0 -1 p2 -1.6 100 -1 p3 -1.6 0 0 ;group hzd直接底
, `% _* C) v. Y3 T* ygen zone brick size 4 20 5 p0 -1.6 0 -5.8 p1 1.6 0 -5.8 p2 -1.6 100 -5.8 p3 -1.6 0 -1 ;group hld老底" z- x w8 V8 t+ k! X1 G
gen zone brick size 4 20 5 p0 -1.6 0 -15 p1 1.6 0 -15 p2 -1.6 100 -15 p3 -1.6 0 -5.8 ;group hx下部岩层. R# Z9 B6 j4 X. y
q! w! `; K4 a% F" U% n. y; ]
;充填体及其上下部分
: p& A0 ?% q& m4 C/ a6 y. Jgen zone brick size 6 20 2 p0 1.6 0 0 p1 4.6 0 0 p2 1.6 100 0 p3 1.6 0 1.3 ;group ctt充填体 5 e( T8 c, b+ T9 ]+ ?" @
gen zone brick size 6 20 4 p0 1.6 0 1.3 p1 4.6 0 1.3 p2 1.6 100 1.3 p3 1.6 0 6 ;group ctz直接顶部+ y& p! y( a- X6 _
gen zone brick size 6 20 5 p0 1.6 0 6 p1 4.6 0 6 p2 1.6 100 6 p3 1.6 0 11 ;group ctl充填老顶
& l3 J1 q$ v6 W7 Fgen zone brick size 6 20 7 p0 1.6 0 11 p1 4.6 0 11 p2 1.6 100 11 p3 1.6 0 40 ;group cts上覆岩层, ?! O8 S# E: n: S9 m/ U
gen zone brick size 6 20 2 p0 1.6 0 -1 p1 4.6 0 -1 p2 1.6 100 -1 p3 1.6 0 0 ;group ctzd直接底
* a j6 c" b7 I5 e; Kgen zone brick size 6 20 5 p0 1.6 0 -5.8 p1 4.6 0 -5.8 p2 1.6 100 -5.8 p3 1.6 0 -1 ;group ctld老底2 D; D' ?1 {1 y
gen zone brick size 6 20 5 p0 1.6 0 -15 p1 4.6 0 -15 p2 1.6 100 -15 p3 1.6 0 -5.8 ;group ctx下部岩层! T5 m8 Q9 i% W) q0 ?4 V+ i L7 A# |* T
;工作面及其上下部分$ I( b5 c1 X# b7 F0 L/ x T
gen zone brick size 20 20 2 p0 4.6 0 0 p1 60 0 0 p2 4.6 100 0 p3 4.6 0 1.3 ;group m 煤体
. c3 ~- n9 P' w- Pgen zone brick size 20 20 4 p0 4.6 0 1.3 p1 60 0 1.3 p2 4.6 100 1.3 p3 4.6 0 6 ;group mz直接顶
1 B9 F9 z% k; agen zone brick size 20 20 5 p0 4.6 0 6 p1 60 0 6 p2 4.6 100 6 p3 4.6 0 11 ;group ml老顶
6 |7 C+ \& A+ Z! R# Rgen zone brick size 20 20 7 p0 4.6 0 11 p1 60 0 11 p2 4.6 100 11 p3 4.6 0 40 ;group ms上覆. P$ V/ o+ u. \: U
gen zone brick size 20 20 2 p0 4.6 0 -1 p1 60 0 -1 p2 4.6 100 -1 p3 4.6 0 0 ;group mzd直接底
; T) `5 j, x" n+ f/ }gen zone brick size 20 20 5 p0 4.6 0 -5.8 p1 60 0 -5.8 p2 4.6 100 -5.8 p3 4.6 0 -1 ;group mld老底7 _, ~; r! Y8 d" s- S8 o y
gen zone brick size 20 20 5 p0 4.6 0 -15 p1 60 0 -15 p2 4.6 100 -15 p3 4.6 0 -5.8 ;group mx下部岩层
/ [! ?& L1 J1 W' L& F;下区段工作面部分
3 ^! A4 N% V6 y& F/ |gen zone brick size 15 20 2 p0 -30 0 0 p1 -1.6 0 0 p2 -30 100 0 p3 -30 0 1.3 ;group m 煤体% _2 N) K* `0 R$ `2 L
gen zone brick size 15 20 4 p0 -30 0 1.3 p1 -1.6 0 1.3 p2 -30 100 1.3 p3 -30 0 6 ;group mz直接顶
+ F, `: \' a4 F$ _( @gen zone brick size 15 20 5 p0 -30 0 6 p1 -1.6 0 6 p2 -30 100 6 p3 -30 0 11 ;group ml老顶; V2 m/ J; @0 ] l- E2 _
gen zone brick size 15 10 7 p0 -30 0 11 p1 -1.6 0 11 p2 -30 100 11 p3 -30 0 40 ;group ms上覆
# ?- P5 K. R h; K. l; X1 @& t* vgen zone brick size 15 20 2 p0 -30 0 -1 p1 -1.6 0 -1 p2 -30 100 -1 p3 -30 0 0 ;group mzd直接底% E6 f' F' B# ]# ^
gen zone brick size 15 20 5 p0 -30 0 -5.8 p1 -1.6 0 -5.8 p2 -30 100 -5.8 p3 -30 0 -1 ;group mld老底
+ |$ F/ K; z7 L; R. }gen zone brick size 15 10 5 p0 -30 0 -15 p1 -1.6 0 -15 p2 -30 100 -15 p3 -30 0 -5.8 ;group mx下部岩层
. C% [! ?) E( h1 h1 u8 ?- ? w5 Battach face( `% P' X+ y' E4 i
plot block state: Y" S* x- K3 v L: u4 {
save moxingjihe1.sav# R5 s( k9 v) ?0 ? J4 _: V& A
mod mohr( W8 O( a" W3 {4 D& C- E/ B9 z% N( [
;煤层9 q- x' ~7 C' B# M
prop bulk=0.4e9 shear=0.133e9 fric=18 coh=2.4e6 tens=1.0e6 range z 0 1.31 ^5 Z8 X1 ?+ O* d8 R1 T$ p
ini density 1300 range z 0 1.3' n* {& C+ {9 l$ |& b% T4 I
;直接顶 泥岩3 J0 W" k# o5 }$ O) x1 M; F
prop bulk=0.8e9 shear=0.4e9 fric=25 coh=6e6 tens=5e6 range z 1.3 6
- P8 e* l5 k, `; M( aini density 2500 range z 1.3 60 E0 ]' S R/ J6 R4 m
;老顶 中砂岩- M. A. ~8 j% @) q5 x
prop bulk=0.55e9 shear=0.32e9 fric=32 coh=7.8e6 tens=9.2e6 range z 6 11
; B+ S# _3 v0 w/ t& Aini density 2550 range z 6 11
% F; M: K- g( ~9 O( ?7 q6 R;上覆岩层 砂岩
' w1 h# j) f- J+ P0 f* v4 Rprop bulk=0.776e9 shear=0.558e9 fric=32 coh=8.0e6 tens=9.0e6 range z 11 40
. N/ }" j) r: g9 p. Pini density 2650 range z 11 401 \. m0 r& P0 O) A- _
;直接底 粉砂岩* E+ O0 K$ J& O& J) P; B
prop bulk=0.55e9 shear=0.32e9 fric=32 coh=7.8e6 tens=9.2e6 range z -1 0, t. y: n, }: q- P& C2 c6 Z
ini density 2550 range z -1 06 o, P% h: s9 E) m
;老底 细砂岩左部
/ m. _, Y# {2 Z( O! r* Gprop bulk=0.536e9 shear=0.369e9 fric=32 coh=8.7e6 tens=13e6 range z -5.8 -1
3 x3 `/ |& y* J7 iini density 2500 range z -5.8 -14 k6 {4 \7 A# C# M0 P: b: `4 H* j
;下部岩层 砂岩* I, F! y" @2 N; J7 m; ?6 ~, z
prop bulk=0.776e9 shear=0.558e9 fric=32 coh=8.0e6 tens=9.0e6 range z -15 -5.8
( Y h2 O2 o& l4 Nini density 2650 range z -15 -5.8
0 |5 o" W9 Y$ X. ] c+ n4 {, Bset grav 0 0 -10
5 N$ r. D' N( efix x range x -60.1 -59.95 d A; [7 {2 `
fix x range x 59.9 60.1
1 {* {- b- ^0 ffix z range z -15.1 -14.9% x; w+ F3 o; l$ l6 Q6 r f
fix y range y -0.1 0.1
, M0 L0 @( [5 B: U1 ?" ^* ^fix y range y 99.9 100.1
5 V! b- ?# V* O* s* k# {;施加应力边界条件$ D5 @) |( B( f
ini szz=-10.6e6 grad 0,0,2.5e4
4 z8 `/ N8 P8 hini sxx=-7.95e6 grad 0,0,1.875e4% J, X: q' j& S
ini syy=-7.95e6 grad 0,0,1.875e4
0 `% h3 M7 f/ _5 w5 i. `5 {% {$ B1 yapply szz=-9.25e6 range z 39.9 40.1
" T6 q: T+ h9 O/ y+ C4 L1 rapply szz=-10.6e6 range z -15.1 -14.96 [% d% }% g! }: i
apply sxx=-7.95e6 grad 0,0,1.875e4 range x -30.1 -29.9
, N1 v4 N# K* [4 x/ m1 Vapply sxx=-7.95e6 grad 0,0,1.875e4 range x 59.9 60.1
* W% T3 ^1 Q6 d8 Japply syy=-7.95e6 grad 0,0,1.875e4 range y -0.1 0.1) X- A, c0 a* P* q- W6 G) o' u9 l- Y
apply syy=-7.95e6 grad 0,0,1.875e4 range y 99.9 100.1$ X8 K: H0 a% @% L
plot block state; b& j7 g/ u/ {
solve
/ s) _- B2 G6 F% P4 psave 初始平衡.sav
5 h; p7 g& G) p$ H' R. _0 |$ ~;位移初始化& N4 `: ^8 [4 W1 o
ini xdis=0
6 R+ c, o; K% m. h3 r; uini ydis=0
K1 x# |- ?- p! m3 S; Zini zdis=0
% Q6 G3 J. j$ x. S4 M0 \- n;挖巷道1 S. D% z3 p8 D9 |) E
model null range x -1.6 1.6 y 0 100 z -1 1.3) @8 ^0 X8 ~6 O- h4 ~
;定义二维数组函数addd，便于确定顶锚杆的施加位置，顶4根锚杆，端部锚杆斜打5 [8 q6 }! A2 U. Z3 r. `
def addd B" N0 u( B9 ^4 x
array xz(2,4) zz(2,4)7 V" R3 r2 a- p# {0 h: s9 i
xz(1,1)=-1.2! d7 v/ D# R/ C+ S
xz(1,2)=-0.4
' t( s. G9 \4 D* a8 Rxz(1,3)=0.4
# h2 B8 v- E9 l- ]: uxz(1,4)=1.2% s& r$ O6 t% s0 W; v( p/ L/ P! e
0 G6 U- C& Z2 |& @zz(1,1)=1.3
C+ d$ g7 Y; V# z9 r1 I( B+ Dzz(1,2)=1.38 v% U+ c% B# j1 {. Z+ D/ T& K0 |
zz(1,3)=1.3
* a% e8 P$ e: ezz(1,4)=1.3
, `( a6 q- U: u8 w4 I7 G& q8 R; H1 ~, G5 | Y
& V/ r' ^( d& \xz(2,1)=-1.9
& c* Y3 ]4 D: V" ]" P0 B O# S! S7 Xxz(2,2)=-0.4. C2 u$ e$ m5 g r1 }, A$ f
xz(2,3)=0.4
0 L' e: C$ M# x: O( jxz(2,4)=1.9
- l; j5 @: A- Y/ B6 ^1 e3 i- M* Y% s
) l8 C: S# k) v% d( tzz(2,1)=2.84
- o# ?( z. o; K. Czz(2,2)=3
2 a ^9 z' z2 T( Q$ }3 Jzz(2,3)=3; V4 V* W, a5 r" c) R: m# m3 L
zz(2,4)=2.848 O) a- ~ }: n9 n+ ~% f+ n; U
' Y) C% L3 A$ S5 Q' }
end
3 ^/ A9 d5 ?6 O3 E0 w
2 u3 E. T5 O; v+ C, p: q0 |def add1% X& u K0 l; p6 w5 S/ t3 }: p1 g
xz1=xz(1,index)% E* ^) n8 J% Y K5 A X3 E
zz1=zz(1,index)
^' g* ~' ^: ~" U
8 D1 r; _/ G0 f9 Y5 D& kxz2=xz(2,index)
4 i+ O8 X) h! P) U# H/ s4 Xzz2=zz(2,index)! w4 y. A5 B$ t. K
end
7 }, F* w2 y9 p# [) W ~0 H% v( ^- `1 W4 w4 w1 T, M2 P6 p
addd* ]. i+ `3 f/ h" ^5 H$ u
;定义函数inscab，采用循环模式模拟计算模型内锚杆的排数9 z, {1 e1 P8 w) Z, @+ x% ^* W
def inscab! Q; O R" Z; j, @: W5 k9 o
6 c y3 i4 P, {4 t0 Q# c) My1=0.2& T" q7 M8 w. j4 F, V
loop m(1,124)
) A- Y: k9 J! ~' }$ n5 b u3 o9 H loop index(1,4)( n$ X4 ^' [2 S, W6 E" r; p
add1
5 o4 V6 K7 ]" O& j/ G6 O command8 I1 u. z% _9 u" \- j- V8 E" m
;设定锚杆的位置
/ R; {6 Q) J4 ?: E5 Z7 C- H! nsel cable id 2 begin xz1 y1 zz1 end xz2 y1 zz2 nseg 4 - I* M1 y. B8 @, L' G* x$ G
;设定锚杆的物理力学参数$ y% \, r, b0 v- v
sel cable prop xca 3.8e-4 e 2e11 yt 1.56e5 gr_k 1.28e10 gr_c 2.25e5 range id 2' ^$ }9 B9 ]0 U( W3 F
end_command% y# ^- c$ r3 J, m4 J6 x& k
end_loop
1 m( ^& `, X0 F;排距0.8m) _2 P0 E+ x3 y
y1=y1+0.8
) g+ C q i# P* v' |! G end_loop
4 V# I8 e4 |* ], s$ k. ?$ P
& T! d5 d" |2 e+ ?$ h. l4 @end
) a. N2 o" ^1 u$ u( u# T;定义addd1数组函数用来确定帮锚杆位置，每帮锚杆3根，帮上上下两根锚杆是倾斜打的
- e4 b6 a: \$ K) wdef addd1
. J+ k6 |( S4 I# q9 ?' ]+ Farray xs(2,3) zs(2,3) xy(2,3) zy(2,3)$ m" c" i F% O) Y9 B0 u# P0 Q9 k
# L3 @& |, @, M4 W' d- R
xs(1,1)=-1.6
4 a) A" N! Y' Cxs(1,2)=-1.6
9 q) k$ E& S3 X; C3 I! ]4 yxs(1,3)=-1.6
2 d# ~( Y9 y9 g% ]' }4 j& M( s: r! C4 T! n
zs(1,1)=-0.5
3 `& o4 |6 u& R5 fzs(1,2)=0.3( G, X9 _+ E: r" N7 G
zs(1,3)=1.1, k1 N3 L7 l5 I# I: m# u3 q
* Q7 k9 F9 x8 }$ N9 s8 cxs(2,1)=-3.32 n! b4 o* J; S j$ Z4 X
xs(2,2)=-3.3
4 e2 ^1 x1 Y! W& \" U3 x0 sxs(2,3)=-3.14
" X2 c& X0 K0 e
" Y* v' n3 o; F2 e6 ~zs(2,1)=-0.5
. Q7 ^) V! l6 F' `+ {* H! vzs(2,2)=0.3" E9 R4 P4 U' R% {; u4 V" ?
zs(2,3)=1.8( W* o- C- t6 Z3 T2 D0 O' X
$ I! ]: `+ f, _" [; k) @& z0 txy(1,1)=1.6
( R. K2 ^% G5 e; ~xy(1,2)=1.6, |- e# W4 z0 O1 x' X4 C' P
xy(1,3)=1.6
% t) a5 N# _; I2 F N8 I- V7 t6 Y$ U8 ]! u. R
zy(1,1)=-0.5
, L; k& X% C4 I xzy(1,2)=0.3/ ~, ]! R& ], t
zy(1,3)=1.1
: n5 [* Q% o) T* U
( D* q9 i+ T$ F) F) w( `xy(2,1)=3.3( }* S, o2 P, ^* k& C, S' }
xy(2,2)=3.3* u4 s d1 [- w" F. a
xy(2,3)=3.14' f4 T' |" [$ V: O& J9 b0 I' c
, q5 q& }+ H' u% k9 _& w/ ]% R a1 b! P
zy(2,1)=-0.5: U n1 T }, I. M( J$ K d! v
zy(2,2)=0.3+ z+ ~, W$ y8 n
zy(2,3)=1.8
" F; `+ X% [6 S# O! K( |/ |
# n* y6 H7 B2 m, o% vend5 _5 Q5 l! F2 Z3 n; g
9 u4 N8 C' \5 qdef add2( L1 z- _, X/ ~3 r/ {
- u: R# d: T0 F( n: j6 S! i2 x# x: C- y
xs1=xs(1,index)7 B6 s+ H W B; V O
zs1=zs(1,index)
8 Z# a m' _0 y( {
* U: m' `+ {* n# Z( w f7 ?xs2=xs(2,index)
' o2 V7 k2 k( ]9 U# _* ~/ {zs2=zs(2,index)
7 M3 C- b5 D& |) l7 m
9 W3 A4 G" l ^! q9 u) dxy1=xy(1,index): t& U; [* }& @1 R) a/ x. i
xy2=xy(2,index)
$ x1 I1 J$ d1 ^/ R4 B
& k F! E* |# u Kzy1=zy(1,index)
+ V! v4 P# n) |9 ]; w1 F! Tzy2=zy(2,index)
; ~9 n6 W! Q# `$ F5 U: H
( g; o9 ^6 U6 s3 lend W' A7 d7 ], o) g. @. s
addd1
9 L! ]# ~7 F" w. B( Gdef inscab1" V8 A3 f" h1 A% w i% |
y0=0.2# p- R4 R* o$ W% E5 U! I$ k q
loop m(1,124)
3 j( Z* t% B9 c; |# C6 @9 D loop index(1,3)
s: z Y6 t) g/ ]5 q& C( W$ f3 j add2' |& U5 y: K n: P- V* m
command
; G) `# L- T; {: e. y; x;施加帮锚杆单元和设定物理力学参数+ D. l0 h, a) _# }, R1 r
sel cable id 1 begin xs1 y0 zs1 end xs2 y0 zs2 nseg 3' L7 t7 p( o b# ]) T; w; z6 r* l5 G7 k
sel cable id 1 begin xy1 y0 zy1 end xy2 y0 zy2 nseg 3
: [3 Y" G& C5 @ sel cable prop xca 2.54e-4 e 2e11 yt 1.04e5 gr_k 1.28e10 gr_c 2.25e5 range id 19 U. m/ U9 x3 b
end_command) z. q# g- T, X4 P5 V% I
end_loop- L, I5 l9 k# k' S( U! G( L: a4 _
# O" R" D0 Y# }$ Q; s, E1 y y0=y0+0.8
" x; n$ S# }8 b. \( o X/ J end_loop
4 W+ _) {( C) x; ?6 `. G( s( Nend- J' j4 k% R4 h1 _& `1 T
7 B* Y+ A. I( i$ Rinscab
4 M) C2 P9 |0 u& \: _% ninscab10 k; t5 B- V& w, O" h( M8 B7 N
;2.5m处巷道变形情况依次为顶板 底板; R; R6 T* w9 b
hist gp zdis 0 2.5 1.3 7 I5 `) T- g v+ d, u5 I
hist gp zdis 0 2.5 -1
, |2 G X9 c) U;50m 处巷道变形情况依次为顶板 底板9 f: V1 B1 @) P2 T3 c# X+ m/ i
hist gp zdis 0 50 1.3
& X, p# u; {5 j" l) x& ~1 Chist gp zdis 0 50 -1
' ~# z7 k5 \ f4 c. ^1 y! l0 M# esolve
" F4 e5 u. A l" M6 ?5 R* usave 巷道成型.sav
' _# a" B( I C0 b;位移初始化
% q* O; n# F- U0 Z, cini xdis=0
* d2 {: M) l3 C% Q2 r; Tini ydis=0+ ^9 ?6 j3 ]4 {
ini zdis=06 | b- Z8 S ]3 E) ^" A
;挖巷道
. Z- E8 i. S! l T% ymodel null range x -1.6 1.6 y 0 100 z -1 1.3
, a% Y9 Q6 z2 L
; O6 Q- A1 t1 K# B( X;定义二维数组函数addd，便于确定顶锚杆的施加位置，顶4根锚杆，端部锚杆斜打
0 q% f$ [2 o, _6 Jdef addd# U9 I1 N. _8 E t7 C0 \' u
array xz(2,4) zz(2,4)4 Y7 M" ?3 q& W4 n1 x
xz(1,1)=-1.2: _% y# f9 \# L, y6 S
xz(1,2)=-0.4
9 r' B' j j# N5 Gxz(1,3)=0.4( x8 n1 H# }, Q( B% D# l2 X, o
xz(1,4)=1.2
, S% ]" i2 v6 E$ p( @
* A/ h- H! o; zzz(1,1)=1.3: l# b! r" c5 R. w l% v
zz(1,2)=1.36 Q0 q+ G' J/ j9 k: \6 I! f' g, A
zz(1,3)=1.37 j( K v4 H. j' `) S
zz(1,4)=1.3* j/ } z8 W* d* p. W
8 I" P+ I- ^+ N/ i* U, ^: F9 J1 j0 r* x6 m& Q& R$ p2 ]" A
xz(2,1)=-1.95 B4 n. v2 [/ ]; x9 l# b6 B) p* x. n4 V
xz(2,2)=-0.4+ O9 J7 C5 \: U* {
xz(2,3)=0.4
/ Z+ n( [3 D3 X% Dxz(2,4)=1.9
- ^* ^& V: Q/ x* i& _9 ~) b3 S
- @1 J! n3 Z. x+ ozz(2,1)=2.842 T$ v5 b/ n9 M0 O- |+ O
zz(2,2)=3, X! g+ y4 _! U) Z% o; S5 F9 I- ^
zz(2,3)=3/ q/ s) @2 m' O* d! j" d; y
zz(2,4)=2.847 `! a T% ]6 |- V, O1 s' Y3 N8 }
( B: r* N+ e' h* Jend
' [5 k& J. m+ Q
% d& P$ \. K" ]- Wdef add1; G' b! {! M* y$ |( o2 [
xz1=xz(1,index)
8 j9 E9 N4 W8 Y& `( G4 _zz1=zz(1,index)8 G: e6 d* j7 b4 D5 m4 v
+ S( o/ `+ J C9 z' |
xz2=xz(2,index)3 ?; f* {9 @/ P9 |! U7 l9 a
zz2=zz(2,index)) ^ o4 _5 z# O, V! S
end
- [3 f- O |8 ` D2 i! j3 q: C
/ m0 H N9 Z& \; B" z& Y' C; iaddd2 q; W; e: B( Q# _
;定义函数inscab，采用循环模式模拟计算模型内锚杆的排数
4 _% @: C! f6 i1 ?6 J) o6 X7 x3 X$ q3 Odef inscab
0 o( z7 l4 ?( f+ I) L& E+ I
4 U! }% v' C" {* I- {: s; dy1=0.2
' @7 G5 V# H" f( Q0 k; v+ @$ x loop m(1,124)% V {& K+ @' U# X
loop index(1,4); x% m% U- G, F! D9 h" T- o' v8 q
add1! ?8 e( @4 k/ Y$ S
command
2 d4 k p! I5 Y/ i! @: f) K;设定锚杆的位置
) U* N7 c# C1 V8 Bsel cable id 2 begin xz1 y1 zz1 end xz2 y1 zz2 nseg 4
$ h& r9 L& L) v8 Y;设定锚杆的物理力学参数- l5 U* b- F7 y3 a8 A* T% ^
sel cable prop xca 3.8e-4 e 2e11 yt 1.56e5 gr_k 1.28e10 gr_c 2.25e5 range id 2
' W; V9 j9 G& M5 b4 H T end_command
1 O# f5 g! a" r! e2 [" t end_loop
) A: ^; R( L7 q; k2 m) ~( p;排距0.8m4 b* X, D: E/ R0 t* r( m6 S
y1=y1+0.8
' j2 D& k7 h0 P. H# n end_loop
1 e! U/ Y! ^4 F6 T C# w% A; E, h! p% M
end
( v$ v: `- ~- Y) `$ ] m- r;定义addd1数组函数用来确定帮锚杆位置，每帮锚杆3根，帮上上下两根锚杆是倾斜打的( l, g# f2 W7 p( p
def addd19 R1 Q7 A% @3 i8 [& A/ B) A+ J L
array xs(2,3) zs(2,3) xy(2,3) zy(2,3)6 d! o* P1 b$ }, S, u
+ k, c2 j3 ^; t
xs(1,1)=-1.6
. c9 D2 n' F% xxs(1,2)=-1.6
. B; [( }7 i; `- k# V6 Gxs(1,3)=-1.69 G0 j; N. b$ y7 t
) `0 w$ l7 o$ S) Z; L5 L$ h% Z7 U% s
zs(1,1)=-0.51 ^4 u* E8 n/ L
zs(1,2)=0.3
6 B- Q) @# n% L e) i; A( ~+ n/ wzs(1,3)=1.12 q7 l# R( G! {; L, A* m5 u' k2 U
" q" E& |5 U' \" {! j2 E/ \
xs(2,1)=-3.3
) l' B, T; p+ s5 A6 y5 Q9 [xs(2,2)=-3.3
$ } @& X3 `. |, r4 [5 gxs(2,3)=-3.14
" b `5 P' Q/ i5 n4 R5 g( _$ D; a0 Q$ `) ]( ?
zs(2,1)=-0.5
' B# H6 d n7 h9 d( czs(2,2)=0.3
/ b& [: X# p# {( \zs(2,3)=1.89 p4 L$ G& d8 J- ~
: e3 n' ^# T* f6 y' V# r: o+ j8 O
xy(1,1)=1.6, _% g( L e( z1 W: n8 S# M
xy(1,2)=1.6
( }! Q h1 y- @! o& X" @xy(1,3)=1.6
; f* _. c J3 I& b3 ~" Z/ O0 u h" M" u u4 F! z' t8 Z
zy(1,1)=-0.5
( F6 b$ @% y1 [% u' `zy(1,2)=0.33 _" D% g1 S; O. N$ t
zy(1,3)=1.16 o6 K9 p3 ?( E, s9 a! Z W5 l- o
' O _1 ]4 [$ ]) v% }! B$ Y
xy(2,1)=3.3
& p! S6 m& T/ K" g# mxy(2,2)=3.3
6 k+ q- _* j! wxy(2,3)=3.148 V' t3 P( v1 l0 o+ U
9 q: Z/ x5 x6 p; A2 n* l# bzy(2,1)=-0.5
4 D) [- D% p D2 bzy(2,2)=0.3) l) T! x+ g1 r* r
zy(2,3)=1.8
7 j: H5 p" Z$ F2 A5 F
! J" p) a5 J2 U' tend6 ?. B: [4 |" c9 y4 Z" L0 c
# U" l" q y0 Y& j5 O! Ndef add2
* p* d$ \; O% u* f6 z% g+ U7 N- g6 I7 Q. k+ g
xs1=xs(1,index)/ h( b: V& b4 g- I" T+ p3 F# U" f
zs1=zs(1,index)
1 N$ \& A: z4 o( f* b4 u' c$ P- x
7 K" R9 F. g W6 @" H; Y# v( S! J5 ~xs2=xs(2,index)5 R L. @8 D; @8 y- [9 C' J
zs2=zs(2,index)
/ N E! i) J9 R
* o; S+ [/ l$ ?' F5 Uxy1=xy(1,index)
+ w0 T0 b6 T+ u. Rxy2=xy(2,index)
. @2 o' G d; l6 w# }7 i
- b$ V8 n2 f% ezy1=zy(1,index)
9 }$ \. x3 s8 ^$ @( h$ V9 N1 izy2=zy(2,index)
5 I2 v( m2 @' x
2 h. b* `4 c0 `, ^end
' P/ ?; I8 u% {+ H; x5 W4 _2 a! Iaddd1, R/ r& z' P7 d" g& b
def inscab1
( r; A. P, O; i3 [ n y0=0.2; Y) m& x+ l7 g8 ~7 y7 Q$ h% F
loop m(1,124)
7 }3 }5 _( o, t" P5 [' O" v loop index(1,3)* x' z$ t8 x4 x$ n, v+ b
add2* f2 Z3 c4 Z3 d7 P
command
; q# D1 {9 S! I) j;施加帮锚杆单元和设定物理力学参数5 s% Z/ O. N" C& i, g$ N
sel cable id 1 begin xs1 y0 zs1 end xs2 y0 zs2 nseg 3
+ M% b) V, F R( p# E% U1 l sel cable id 1 begin xy1 y0 zy1 end xy2 y0 zy2 nseg 3
( J: |, `1 t- {1 \7 I sel cable prop xca 2.54e-4 e 2e11 yt 1.04e5 gr_k 1.28e10 gr_c 2.25e5 range id 1& b# k0 a' |1 E8 r
end_command
: t* `7 x8 a3 e3 b* o end_loop0 B+ |* g" C0 ?' W- L5 [$ `
& j a. X$ W8 j4 {& K, E y0=y0+0.8
; P4 D/ S/ ` Q( N end_loop
l; ]! d N* W3 N$ G1 R8 @end/ n7 z6 q+ P1 c0 |
0 c" v- H' [. d9 X. F& |( v, \8 xinscab
; ~; M! q. @9 linscab1
/ D9 ]0 t5 [# c) Z' @1 j% I;2.5m处巷道变形情况依次为顶板 底板
* z5 ` p! A- G7 `8 k4 A7 xhist gp zdis 0 2.5 1.3 8 [# z, z& a# p% Q
hist gp zdis 0 2.5 -1
3 }( `: Z- z! D- j1 H;50m 处巷道变形情况依次为顶板 底板# V: y z3 Q$ b+ d% D" l8 V3 j: {
hist gp zdis 0 50 1.3
0 y1 f! N2 ~4 J& r- R, i% D thist gp zdis 0 50 -1) k' B* v/ ~+ \3 o: T" B) }
solve
3 j: F+ Z" f" L, Gsave 巷道成型.sav
1 V3 u' D% J# \: Y+ P/ l [, D, L
4 i% }3 G& N4 Z! b# n0 Dhist reset E# \! y/ f. _7 O% L# W
set logfile kaiwatanmo.log
4 m7 X0 I1 M8 n' t" |% `: i L* Yset log on8 g# A; K. Z* R
;位移初始化
, \' A9 O7 G6 L* L. S8 k; a+ gini xdis=0 2 K4 \( C- q5 U2 B
ini ydis=0
B6 Q/ Q/ [( I& C) C) ?ini zdis=0
; O0 p! b; P( L. j: H! E( `9 C' y! m;提前开帮0-5m4 I+ @% J( o5 a7 l! }+ |# X
ini sxx=0 range x 1.6 3.1 y 0 5 z 0 1.3
7 }9 x4 k9 I% J$ B. n; ?ini syy=0 range x 1.6 3.1 y 0 5 z 0 1.3
2 }! n- w) y( lini szz=0 range x 1.6 3.1 y 0 5 z 0 1.35 \/ D# B& ?) r; p. ^
model null range x 1.6 3.1 y 0 5 z 0 1.38 C" \3 |' D0 G, i6 G6 i# X; C
9 a& T) C9 P _' C4 _; B
( \+ K4 u' } @* \2 p; V7 N( ~
! H" u+ L$ U* [$ s7 u;充填
6 P" Q5 j1 J P& t: }model mohr range x 1.6 3.1 y 0 5 z 0 1.39 V( k7 \( M2 I( ]
6 w7 [6 g6 N* w6 U0 R/ q Y `
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 0 5 z 0 1.3 4 |: A4 e d6 I
ini density 3000 range x 1.6 3.1 y 0 5 z 0 1.3 ! _) y* |" |$ I; P/ t; R4 t1 }( }
, ~/ E% p* a9 @+ N" D
step 50
+ S8 l' e s8 a ~ B* x- ^;监测位移% _1 Q' k- P9 M) s3 U
hist n=5
9 A! x# H+ L. ?: |1 ]9 p( |;顶板下沉、底臌变形/ @2 s( A0 t6 k
hist gp zdisp 0 2.5 1.3 ;1 顶沉
% a# D3 @: D8 jhist gp zdisp 0 2.5 -1 ;2 底臌' b7 h0 G/ Y: T
;四角位移 左上角 右上角 左下角 右下角8 q8 T7 s; D- N$ z/ z& S8 b. Q
hist gp zdis -1.6 2.5 1.3 ;3 左上角; B' t. q% p" y
hist gp xdis -1.6 2.5 1.3 ;4 左上角
# T% e5 K7 J* rhist gp zdis 1.6 2.5 1.3 ;5 右上角
5 K. Q$ q) P/ i& S3 C4 X7 zhist gp xdis 1.6 2.5 1.3 ;6 右上角
' v2 Q9 @1 f* [. u+ vhist gp zdis -1.6 2.5 -1 ;7 左下角- I/ u" }; n3 f1 e5 U0 z3 N; i
hist gp xdis -1.6 2.5 -1 ;8 左下角8 E; a, P' V. U7 W& @- C
hist gp zdis 1.6 2.5 -1 ;9 右下角, |& q% v; x% T: x9 v4 x1 y
hist gp xdis 1.6 2.5 -1 ;10 右下角+ ]- [& B. `6 `& k1 u
;两帮移近量 左 右8 k& D# I- N& A4 a2 t: y% C
hist gp xdisp -1.6 2.5 0 ;11 左帮
& E" }, t9 c2 R' ]1 u7 O1 J2 bhist gp xdisp -1.6 2.5 0.6 ;12 左帮
8 b$ j6 ?" S. ?, vhist gp xdisp 1.6 2.5 0 ;13 右帮
P: ]. x/ V" m0 [hist gp xdisp 1.6 2.5 0.6 ;14 右帮
, P& U7 R# E* }% N
# l$ g% n* O+ b;监测应力* \: Q( U: U1 _6 r- T7 E" v4 h
;充填体应力& P m! B4 g* b
hist zone szz 1.7 2.5 1.3 ;15 内侧顶压$ ?7 r5 h7 D. z/ ^
hist zone szz 3.0 2.5 1.3 ;16 外侧顶压. P( h7 M1 e1 \7 S5 U
hist zone szz 1.7 2.5 0 ;17 内侧底压3 q( J1 L1 ?/ ^
hist zone szz 3.0 2.5 0 ;18 外侧底压 j- [# a& J. a
;煤体应力监测
( j0 g+ Q/ R R/ b% T/ s$ dhist zone szz -1.7 2.5 1.3 ;19 近煤壁处顶压力
& N& j. f% ~6 O0 J0 M% w. P) phist zone szz -1.7 2.5 0 ;20 近煤壁处底压
: U! B: C% H O: g4 Ihist zone szz -3.7 2.5 1.3 ;21 2m内部压力
6 h7 V4 L+ h. ^hist zone szz -6.6 2.5 1.3 ;22 5m内部压力
6 u! b& k' C$ G$ ^8 v4 Vhist zone szz -10.6 2.5 1.3 ;23 9m内部压力) @8 x8 p5 r& J
hist zone szz 0 2.5 1.4 ;24 顶板压力
a, X8 [3 X! h0 }3 mhist zone szz 0 2.5 -0.9 ;25 底板压力3 o7 b8 H! `% |! D- b/ _) X% K
hist gp zdis 0 35 1.3 ;26 35m处顶板位移3 i/ G0 q2 |4 S7 t9 k" y- ]
hist gp zdis 0 35 -1 ;27 35m处底板位移8 E& t' j+ r& n* a8 [4 x% M0 _4 ]
hist zone szz 1.7 35 1.3 ;28 35m处工作面侧垂直压力- h4 M y! u# P( ^" @
hist zone szz -1.7 35 1.3 ;29 35m煤壁侧垂直压力
: z5 L2 d, k6 d# s/ h;开挖上工作面0-5m煤层. n o1 }/ q8 C: f; q
model null range x 3.1 60 y 0 5 z 0 1.3% E$ n7 c! v M& l/ f2 |
;采空区处理
" S. B/ g) K X& q( F9 i- hprop bulk=0.4e9 shear=0.2e9 fric=25 coh=8e6 tens=6e6 range x 3.1 60 y 0 5 z 1.3 6
0 T6 Y% w* b9 b6 H$ p6 |ini density 2500 range x 3.1 60 y 0 5 z 1.3 6
+ j$ A, y! r2 b: Rset large% i) P$ B4 M) o, P! T; L, J
solve ratio 4e-48 [; f, z4 s* o& q
save 开挖充填5m.sav
" u8 p+ h- Q% @4 w2 ^3 _/ w2 u3 ]2 J! D1 V' D+ h C
;提前开帮5-10m, M, ^- a" c6 P6 X" ~; p
ini sxx=0 range x 1.6 3.1 y 5 10 z 0 1.3
) V- J3 @0 X* ~: cini syy=0 range x 1.6 3.1 y 5 10 z 0 1.34 |/ V0 _# k8 J
ini szz=0 range x 1.6 3.1 y 5 10 z 0 1.3
& U% M; n4 s: {& ~! [model null range x 1.6 3.1 y 5 10 z 0 1.3% w+ }3 u* n2 D$ V. l5 H, d5 i
% H) S* m n6 O) F% O5 A
model mohr range x 1.6 3.1 y 5 10 z 0 1.3: C% |/ z A) S. A; L5 U
`* f8 e2 t9 E6 F% ]
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 5 10 z 0 1.3
5 W5 u- C2 L2 R8 `6 _- V5 |ini density 3000 range x 1.6 3.1 y 5 10 z 0 1.3 $ h9 S& A' y$ k8 v: d7 I+ L/ x6 w4 {
) ^% r- y ?' N; D; \' G% e) i, h
step 50
4 |' k- f: ]- w: e" o8 k;开挖上工作面5-10煤层2 T/ ~$ D, i- W4 @' Y$ D/ ?
model null range x 3.1 60 y 5 10 z 0 1.3" ^' R: |% {( M
;采空区处理* `& x( [; \( Y! Y! o
prop bulk=0.4e9 shear=0.2e9 fric=25 coh=6e6 tens=2.5e6 range x 3.1 60 y 5 10 z 1.3 6
, }( Q9 a# a! h1 g+ @ini density 2500 range x 3.1 60 y 5 10 z 1.3 6
" H& c. w- X- U# @- F1 e7 zsolve ratio 4e-4
3 q J I3 o' U1 z& R& Usave 开挖充填10m.sav
2 o( y H/ {- O" {% E4 ];提前开帮10-15m' G+ S0 z5 _( F3 W) ?' R' @4 A
ini sxx=0 range x 1.6 3.1 y 10 15 z 0 1.3
& F# ^7 B P) Xini syy=0 range x 1.6 3.1 y 10 15 z 0 1.3' e) _9 ?1 u# {% V+ s( e: ]6 S5 s
ini szz=0 range x 1.6 3.1 y 10 15 z 0 1.38 c( j& P& `7 E8 x! X) z0 W* `
model null range x 1.6 3.1 y 10 15 z 0 1.3
9 K/ H5 ?8 i+ D$ l$ W/ g% Y# j8 O, h. N; j1 ~2 g9 x5 H
;充填
3 _7 r# ]4 S- ^6 D$ V; ]model mohr range x 1.6 3.1 y 10 15 z 0 1.3
2 Y. s( g+ F' ]9 y7 G) V& c# i! mprop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 10 15 z 0 1.3 8 _7 ?4 B- d: J1 K
ini density 3000 range x 1.6 3.1 y 10 15 z 0 1.3
8 S/ Q4 O0 d# E7 L+ f9 Ostep 50
9 a! {( x' I+ F;开挖上工作面10-15煤层2 m" n' I. D# T7 J) m5 A! u* V6 C; B& Q
model null range x 3.1 60 y 10 15 z 0 1.3
# C2 } v d- k; A: E8 q- p;采空区处理9 K% B. o' I' C& K
prop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 10 15 z 1.3 6
+ h8 I; o; c/ `ini density 2500 range x 3.1 60 y 10 15 z 1.3 6
4 b) m6 Q: |. L- K. a5 esolve ratio 4e-4/ `0 t1 H& P4 q6 |8 Z
save 开挖充填15m.sav
$ [+ M. D9 j+ y- ~/ ]7 F' m! O P5 z& [: a
$ P `4 \! _3 E/ [! V
1 ]1 Z' b# a3 y) [0 Z4 B+ p& V! x+ @2 h3 Q& Y
;提前开帮15-20m
2 z, s' g6 s6 a! C& L- ?1 pini sxx=0 range x 1.6 3.1 y 15 20 z 0 1.39 v3 Q9 C0 T# }" S
ini syy=0 range x 1.6 3.1 y 15 20 z 0 1.3
4 Z4 G+ d2 H3 I9 ]- Qini szz=0 range x 1.6 3.1 y 15 20 z 0 1.3
4 o8 I s6 L, b* ?4 J- wmodel null range x 1.6 3.1 y 15 20 z 0 1.3
% f% a" _* c. Q5 k0 Z
6 D9 F0 C8 S: d2 O* U$ U* R;充填( R/ u6 ^7 L/ e
model mohr range x 1.6 3.1 y 15 20 z 0 1.3
# u7 F6 w% c. |2 i! ~( d1 y$ [prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 15 20 z 0 1.3
$ t% J2 H3 H# n7 q0 x3 |( Cini density 3000 range x 1.6 3.1 y 15 20 z 0 1.3
0 N3 X+ {* b, Y: q# t @( S6 @5 a+ G
step 50
: N, }' o9 U& Z B0 p+ U* X;开挖上工作面15-20煤层
% V5 I- y, s, f/ Fmodel null range x 3.1 60 y 15 20 z 0 1.3
$ R p# N4 A. m. E;采空区处理/ x- s9 Q( B) g# Q `8 N0 \1 o$ ^- b2 b
prop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 15 20 z 1.3 6
# d* O4 P7 X3 x2 Mini density 2500 range x 3.1 60 y 15 20 z 1.3 6! ]- Q; e7 `; @* X% w; S+ i+ @
solve ratio 4e-4
4 ^3 Q6 C, `- h! j! z: W8 Lsave 开挖充填20m.sav
6 q& y5 Q5 |5 K" B3 u0 G; ^& I% P. r% d5 p
;提前开帮20-25m7 D/ f, y ]; ~
ini sxx=0 range x 1.6 3.1 y 20 25 z 0 1.3- m/ }+ s0 z4 v9 X1 a ?
ini syy=0 range x 1.6 3.1 y 20 25 z 0 1.36 o* o1 r' u3 ]0 R( N1 @
ini szz=0 range x 1.6 3.1 y 20 25 z 0 1.3
; w+ o) j+ w; Smodel null range x 1.6 3.1 y 20 25 z 0 1.37 R+ k* [* Y6 v* c; P
( B; }* t w4 }) C; ^+ C;充填
+ }+ z. X: z0 ?' G' n+ o% V& i6 {model mohr range x 1.6 3.1 y 20 25 z 0 1.3+ x; F! O, A1 k
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 20 25 z 0 1.3
+ d7 _. [/ w, H2 ?2 Y" {ini density 3000 range x 1.6 3.1 y 20 25 z 0 1.3 , T- D5 r! A1 h6 L9 t
step 50$ U, G0 \# n0 R# m) x3 O Y8 T& q8 u0 O6 h
;开挖上工作面20-25煤层
9 j6 V w+ h/ a+ }) bmodel null range x 3.1 60 y 20 25 z 0 1.3
4 x7 K, g. _- W3 x( R8 Y/ [& H$ V;采空区处理
* n& T+ M a5 u- L- V* l' x, o* hprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 20 25 z 1.3 66 J8 z/ @( n: k0 ~) ]
ini density 2500 range x 3.1 60 y 20 25 z 1.3 6
( S* K* E4 f% bsolve ratio 4e-4
* f, c: |0 w5 ssave 开挖充填25m.sav8 Q# \; X# a: t! e
( {: }* A$ B: \% u
;提前开帮25-30m
" v' q. ^9 F& b4 G# V% Qini sxx=0 range x 1.6 3.1 y 25 30 z 0 1.3
* I. m: c# E o0 w/ B! n# \ini syy=0 range x 1.6 3.1 y 25 30 z 0 1.3& l- I8 V- h4 H4 N+ x& m0 k# c! U
ini szz=0 range x 1.6 3.1 y 25 30 z 0 1.3( ~2 c8 U/ ]: I* z4 ~7 k" o
model null range x 1.6 3.1 y 25 30 z 0 1.38 m8 D* t! o- K
5 c, V; }( g. q" z+ }3 g;充填* N/ p0 ?. i' p4 p: |
model mohr range x 1.6 3.1 y 25 30 z 0 1.3; _( Z. m5 _& U7 R% [( s& @' m1 s
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 25 30 z 0 1.3
9 a7 {* G$ v' a. T, q1 ?9 N: H) L sini density 3000 range x 1.6 3.1 y 25 30 z 0 1.3 ; s n9 y# w' h$ Z) m% G* Y7 ~/ ~
7 _0 e/ | f9 R- K$ k* Q" Z
step 50
* z6 }( p$ l/ p5 A7 r1 N) u4 a;开挖上工作面25-30煤层
I4 i% g# E0 w8 u# q" ]: Tmodel null range x 3.1 60 y 25 30 z 0 1.3
B5 z! f6 f' k/ p0 |/ Q) a; |: \;采空区处理
1 F+ ]! C' F) U! v' U7 k2 Dprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 25 30 z 1.3 6' Y1 v& g- k% m6 e* U& X
ini density 2500 range x 3.1 60 y 25 30 z 1.3 62 }; ~# [* E$ u) v! S: R, F
solve ratio 4e-46 V7 ?5 Z, T9 E
save 开挖充填30m.sav( l9 R" A- j2 M& i9 [( h2 K" F. A& f
. X4 b2 X G3 j& u+ W+ w
' Y7 D# {4 { n: F- f' L6 t- r* W" W& |, p: X9 g
t8 N: D. T6 `1 y
- ~4 F6 l6 T5 j5 s6 N9 O+ u* m4 N8 [3 H
;提前开帮30-35m5 _3 R& b4 R: b) E/ T, e9 E
ini sxx=0 range x 1.6 3.1 y 30 35 z 0 1.3% a r+ m; ` G& u$ V! [+ _
ini syy=0 range x 1.6 3.1 y 30 35 z 0 1.3
. ?6 ]4 v6 i. Q" F5 j! S8 h2 S% ]ini szz=0 range x 1.6 3.1 y 30 35 z 0 1.3, d9 b- _2 x' b' f1 A
model null range x 1.6 3.1 y 30 35 z 0 1.3% l: R" u! m6 r
+ k. x3 R+ i' R
;充填
2 k& `9 T8 Z+ t9 v, Vmodel mohr range x 1.6 3.1 y 30 35 z 0 1.3& r5 V) S0 {: o" r3 @
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 30 35 z 0 1.3
+ {# n4 D9 f/ H! F+ a& y4 aini density 3000 range x 1.6 3.1 y 30 35 z 0 1.3
& D4 u: W* m% t m3 f; b
: _: m7 B2 p4 z* ^- istep 50% l- t- t0 C2 Y5 w) K! H' C
;开挖上工作面30-35煤层; Y1 @& s+ P7 y- s7 n C" {6 h" r
model null range x 3.1 60 y 30 35 z 0 1.38 Y1 h/ c4 d7 r# P; G
;采空区处理4 {: M# z; R! k1 Q, }
prop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 30 35 z 1.3 6
1 t: q& U- r+ e3 e6 |ini density 2500 range x 3.1 60 y 30 35 z 1.3 6) i5 o* J( S2 B8 _4 a) ?
solve ratio 4e-42 L# V! Z- L; v4 o
save 开挖充填35m.sav
9 j9 P- P: p" _) k' ^% V7 a! t- ^2 F: @+ X4 r7 H0 P: _& ^ J
;提前开帮35-40m, G7 {6 z+ _9 G! _2 W- |7 M
ini sxx=0 range x 1.6 3.1 y 35 40 z 0 1.3
6 x* |( X0 i0 q3 ?. B8 {5 {$ I) X6 rini syy=0 range x 1.6 3.1 y 35 40 z 0 1.3
. R8 G" c& G/ i+ o1 k2 d2 d' Gini szz=0 range x 1.6 3.1 y 35 40 z 0 1.3$ X5 \6 ~& F" s2 m- K8 X& c
model null range x 1.6 3.1 y 35 40 z 0 1.3; D2 ?5 S$ \; P- |
1 X! i9 a, M5 k9 O2 g7 y
;充填0 }4 ~" O8 R4 n) Q, b z+ D% ~& w
model mohr range x 1.6 3.1 y 35 40 z 0 1.3
1 b8 f( @* m% q d; b+ i! T E9 Jprop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 35 40 z 0 1.3 $ S1 |4 K& F `- L M- e& R# j
ini density 3000 range x 1.6 3.1 y 35 40 z 0 1.3 + ?, ?- J& E5 L9 A. t, X9 g A& R& [4 {
+ H, W% H9 J- t# \, N* O, t
step 503 x9 n; [9 H3 V6 o* a' f
;开挖上工作面35-40煤层% E3 }5 z, C5 I. L1 A$ r( L
model null range x 3.1 60 y 35 40 z 0 1.3
; c3 }2 P F. N& z2 q;采空区处理
& Q0 ^) n9 c9 m& D3 _; f) Tprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 35 40 z 1.3 64 _+ j. b* ~/ P7 Y' _! l
ini density 2500 range x 3.1 60 y 35 40 z 1.3 6
( f/ ~. O+ x' F1 bsolve ratio 4e-4
% X7 S0 \# B% J( ?8 f) dsave 开挖充填40m.sav
4 _2 i |# v# P
0 ^* H6 ]' d; J1 j5 M2 U;提前开帮40-45m
4 ?, I k* Q3 h* x, {4 e$ y$ Gini sxx=0 range x 1.6 3.1 y 40 45 z 0 1.3
: W9 y0 k4 }; b& Yini syy=0 range x 1.6 3.1 y 40 45 z 0 1.3
, W4 [ q3 n- Wini szz=0 range x 1.6 3.1 y 40 45 z 0 1.3
2 t% q2 }! j! j- R p/ ~model null range x 1.6 3.1 y 40 45 z 0 1.3
4 _, ?6 M$ c( L: R7 z% ?" x3 x4 r6 P3 Y4 m6 S) V2 E" R
;充填+ E3 ] w9 j) I# _4 R- }
model mohr range x 1.6 3.1 y 40 45 z 0 1.3
( K* y- N8 N5 X! l8 Zprop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 40 45 z 0 1.3 ( z; n' N5 h0 N/ M. ^- B
ini density 3000 range x 1.6 3.1 y 40 45 z 0 1.3 , x. n' f5 y. p; g0 ~
+ d4 _' ~0 {7 F% B+ N$ @+ wstep 50
. U$ U1 J) j+ @6 f- A. T;开挖上工作面40-45煤层" q: X% B! k+ V; s6 D/ O. R
model null range x 3.1 60 y 40 45 z 0 1.3 {. U4 E# H4 L6 I! f" F2 d
;采空区处理' t7 F8 c3 W2 ?6 ~* Q8 x: i. X
prop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 40 45 z 1.3 6) B8 m2 C7 u( `, [' V
ini density 2500 range x 3.1 60 y 40 45 z 1.3 6
* k W$ F# R9 X, ?- E+ I) U+ v8 |solve ratio 4e-4
/ x4 [8 I+ @. n' S4 s# wsave 开挖充填45m.sav
, H- R8 O: d* f7 u8 D U
4 P- |" m+ ?/ T! f
, Q! M0 c8 Q6 S7 j;提前开帮45-50m/ _4 z" D* x9 E# E
ini sxx=0 range x 1.6 3.1 y 45 50 z 0 1.38 V7 I. u, s+ d f6 i
ini syy=0 range x 1.6 3.1 y 45 50 z 0 1.3
) P" l$ ^% _/ Q2 X1 `ini szz=0 range x 1.6 3.1 y 45 50 z 0 1.3/ A- k1 z9 J8 R# _6 |0 ]1 C6 N4 U
model null range x 1.6 3.1 y 45 50 z 0 1.3
/ M# Q, E! P, Y5 r5 \
0 R0 [- u2 \! {, l4 ]9 B9 W;充填1 N1 N" S4 C. k: ^8 Q+ y- `
model mohr range x 1.6 3.1 y 45 50 z 0 1.3
6 X: U) N L- a( x, K' aprop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 45 50 z 0 1.3
7 m- |; Z) z; J" z1 p& U( yini density 3000 range x 1.6 3.1 y 45 50 z 0 1.3
8 t1 f2 M- p% Q: O9 t' x# J, m+ r, Y3 T% I& {
step 50# f4 [3 `8 ^1 K& r& l6 d( {% ]
;开挖上工作面45-50煤层( H& `7 v: l8 a2 y: O _8 d: g
model null range x 3.1 60 y 45 50 z 0 1.3
6 j1 ~) f& V% A+ m;采空区处理
, K& c _/ F& @4 @6 Oprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 45 50 z 1.3 61 @' m- W0 \5 c* p* T, m
ini density 2500 range x 3.1 60 y 45 50 z 1.3 6
: w, t- s$ q, n6 v4 xsolve ratio 4e-4; c/ K; K6 d& c4 D
save 开挖充填50m.sav
$ c+ m- G( ?3 _, {4 |. ]5 z/ z x/ q
& E5 [; w/ x Y/ D! b/ L$ P5 |9 a/ M: a7 B9 t6 u3 Y
; t( S! M: ?) J7 s/ B! g" K
;提前开帮50-55m7 c t: Q7 f* U- D% j4 ^% A
ini sxx=0 range x 1.6 3.1 y 50 55 z 0 1.3
1 N& l# C7 @9 {* z! u7 cini syy=0 range x 1.6 3.1 y 50 55 z 0 1.3
2 P" u* W$ Y4 ~9 k1 |- N3 S5 vini szz=0 range x 1.6 3.1 y 50 55 z 0 1.3+ y( a8 [8 t4 R0 G, [
model null range x 1.6 3.1 y 50 55 z 0 1.37 u; H, `# v+ P# l* `
& l: b9 w4 B% ?$ [) s+ p' F8 `8 {;充填+ I4 f# ]* s) s0 j
model mohr range x 1.6 3.1 y 50 55 z 0 1.3
# h0 h& }+ J5 [2 l0 d9 dprop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 50 55 z 0 1.3
( k5 x9 \4 m; N+ {ini density 3000 range x 1.6 3.1 y 50 55 z 0 1.3
/ K9 X8 Y* x3 r2 }) P2 A" ~* a1 L8 h# u5 h$ P0 ^8 k' Y* A+ `
step 50
( G1 ~& z* t# @* s* P;开挖上工作面50-55煤层
/ C& X) h. {- Y ?( S0 @/ j. M/ wmodel null range x 3.1 60 y 50 55 z 0 1.3# V1 m6 g/ t; k: s: r. `
;采空区处理
- Y+ Q P8 N$ Jprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 50 55 z 1.3 6
* _$ p8 ], O7 B4 A2 t7 Qini density 2500 range x 3.1 60 y 50 55 z 1.3 6" S% [ S- d7 e. h9 q' b5 A
solve ratio 4e-42 l5 d" @' M7 G' C% M2 l1 T
save 开挖充填55m.sav
. W! r6 C: s+ }! k. S E9 d9 R D8 H) V
;提前开帮55-60m
! q7 T, D( B o7 l4 j* C7 yini sxx=0 range x 1.6 3.1 y 55 60 z 0 1.3
( F1 J- Z C& o6 K' rini syy=0 range x 1.6 3.1 y 55 60 z 0 1.3
+ @. H# F# f- D+ }ini szz=0 range x 1.6 3.1 y 55 60 z 0 1.3. o) r+ m! N7 x5 h, J1 ~2 G
model null range x 1.6 3.1 y 55 60 z 0 1.3' r( y# c, P1 m/ J. ~6 C5 y
* D @& |4 ?4 e3 J6 d
;充填
% P! N$ Q, B# J/ amodel mohr range x 1.6 3.1 y 55 60 z 0 1.3& s/ U6 m. ` `. ], V1 s; M: y# k5 `# Y
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 55 60 z 0 1.3
[! @4 x: ^$ b6 N2 M" M, X- jini density 3000 range x 1.6 3.1 y 55 60 z 0 1.3 + R1 \+ t( f( [6 a+ ?& l- w" d
;开挖上工作面55-60煤层
6 w/ X K8 C+ p% a8 c7 R: imodel null range x 3.1 60 y 55 60 z 0 1.32 R! J% r# K q$ a
;采空区处理
- @8 e; `& {& [: p3 Y3 t* Bprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 55 60 z 1.3 6
# i" ]! i! C5 g+ U' j% D+ l8 Zini density 2500 range x 3.1 60 y 55 60 z 1.3 6$ }% V% n7 x3 C2 o
6 i: o7 d# I) `6 h6 O
solve ratio 4e-4
6 I1 \7 G2 d& i4 psave 开挖充填60m.sav
9 @' @4 P+ \& B$ ^5 [3 D
( H$ e+ b2 q& k4 h7 C; B;提前开帮60-65m
5 c, B( a% |- U( f+ Pini sxx=0 range x 1.6 3.1 y 60 65 z 0 1.3
1 ^ n$ f3 ?6 R O7 v$ nini syy=0 range x 1.6 3.1 y 60 65 z 0 1.32 A. s: Y! J, T2 i6 b+ L
ini szz=0 range x 1.6 3.1 y 60 65 z 0 1.3
: I8 A5 S# N" ?) B2 _model null range x 1.6 3.1 y 60 65 z 0 1.3
0 b. W0 _# ]) J/ c' Z' E- L0 i6 t% Q% {2 |4 e, O B( R
;充填
. @: D; ~" _& {) {1 Amodel mohr range x 1.6 3.1 y 60 65 z 0 1.3 C& Y- W" z* n
prop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 60 65 z 0 1.3 6 l9 D, L5 r! x$ ]( Z- s( K
ini density 3000 range x 1.6 3.1 y 60 65 z 0 1.3
$ ^# S2 \7 Q0 G* a* l( S2 X1 O" c0 q( Y/ g" [4 F4 o6 ?- f
step 50* d" k/ P! l/ W
;开挖上工作面60-65煤层
2 k( c+ x4 a' q: b/ V# L" Ymodel null range x 3.1 60 y 60 65 z 0 1.3
1 s$ o4 N0 h* B& V" \2 @;采空区处理
% h) U$ ~+ y0 P5 ^2 |prop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 60 65 z 1.3 6
7 s Y) o7 B% A: R! O0 Y. `ini density 2500 range x 3.1 60 y 60 65 z 1.3 6
+ x; v! u" P# z& Gsolve ratio 4e-4
* U6 t6 z$ I$ s; dsave 开挖充填65m.sav% V' S. f( H# M4 |' O. D# a: W
. F1 D9 X) |; g- r+ _' ^5 b
;提前开帮65-70m
8 G Q, d) f- u7 e( Gini sxx=0 range x 1.6 3.1 y 65 70 z 0 1.3$ ?6 J3 z# h9 j! L% m) P
ini syy=0 range x 1.6 3.1 y 65 70 z 0 1.3& u, b7 e4 K8 V0 \# G4 ~9 E# N
ini szz=0 range x 1.6 3.1 y 65 70 z 0 1.3
) x$ T! B( a6 {7 N" K* a- ?model null range x 1.6 3.1 y 65 70 z 0 1.3- l5 v0 O) ]' s- h: B( ~
$ m& g/ s1 O0 [. z
;充填
2 `: I4 a( ] I! Jmodel mohr range x 1.6 3.1 y 65 70 z 0 1.3
I- a5 d* Q% y) W b& \! @* Gprop bulk=1e9 shear=0.8e9 fric=34 coh=10e6 tens=12e6 range x 1.6 3.1 y 65 70 z 0 1.3
0 H! U0 S' u8 z1 t, |0 Jini density 3000 range x 1.6 3.1 y 65 70 z 0 1.3
' T: f3 M9 O, b0 ?' ]- P; W- b- y& ^) b' T! I8 J
step 50
; U3 w8 ?, m! g' p# J* `;开挖上工作面65-70煤层/ b, O$ m& P( s+ x
model null range x 3.1 60 y 65 70 z 0 1.3# U8 Q+ E7 X8 R
;采空区处理
% H9 x' B2 y& q8 Rprop bulk=0.4e9 shear=0.2e9 fric=25 coh=4e6 tens=3e6 range x 3.1 60 y 65 70 z 1.3 6
& v5 H, ~4 s3 m4 L! uini density 2500 range x 3.1 60 y 65 70 z 1.3 6
) u ~( G! J: G- b1 Rsolve ratio 4e-4
+ n/ H9 Y6 t/ U) ]7 i- Esave 开挖充填70m.sav |
|
IP卡
狗仔卡
显身卡
