- 积分
- 0
- 注册时间
- 2007-6-3
- 仿真币
-
- 最后登录
- 1970-1-1
|
;计算范围为50*50M 区域,划分单元10000个
grid 100,100
;设置材料模型
model mohr-coulomb
;生成网格
;直接底
gen 0,0 0,8 50,23 50,0 i=1,101 j=1,25 rat=1,0.9
;煤层
gen 0,8 0,16 50,31 50,23 i=1,101 j=25,49 rat=1,1
gen 0,16 0,17 50,32 50,31 i=1,101 j=49,52 rat=1,1
gen 0,17 0,26 50,41 50,32 i=1,101 j=52,76 rat=1,1
;直接顶/老顶
gen 0,26 0,50 50,50 50,41 i=1,101 j=76,101 rat=1,1.1
gen arc 25,26.4 27.2,26.4 182
gen line 22.8,25 22.8,26.4
gen line 22.8,25 27.2,25
gen line 27.2,25 27.2,26.4
;岩石材料模型描述
;底板(粉砂岩)
prop d=2458 b=9.609e9 s=6.909e9 fric=54 coh=4.3e6 ten=2.58e6 dil=0 &
i=1,101 j=1,25
;煤层
prop d=1366 b=1.293e9 s=0.89e9 fric=40 coh=2.54e6 ten=2.17e6 dil=0 &
i=1,101 j=25,49
;夹矸
prop d=2000 b=5e9 s=3e9 fric=40 coh=3e6 ten=2e6 dil=0 &
i=1,101 j=49,52
;煤层
prop d=1366 b=1.293e9 s=0.89e9 fric=40 coh=2.54e6 ten=2.17e6 dil=0 &
i=1,101 j=52,76
;基本顶(细砂岩)
prop d=2662 b=12.103e9 s=10.642e9 fric=39 coh=9.86e6 ten=8.18e6 dil=0 &
i=1,101 j=76,101
;初始化
;问题描述
set large
set g=9.8
;边界条件设定
fix x y j=1
fix x i=1
fix x i=101
;初始条件设定
initial sxx=-2.0e7
initial syy=-2.0e7
ini xd 0 yd 0
;荷载条件设定
apply syy=-2.0e7 i=1,101 j=101
apply syy=-2.0e7 i=1,101 j=1
apply sxx=-2.0e7 i=1
apply sxx=-2.0e7 i=101
;模拟开挖
model null reg i=51 j=58
;支护结构单元得设定
;锚杆
stru cable begin g 43,55 end g 47,55 prop 1 t=12e4
stru cable begin g 43,59 end g 47,58 prop 1 t=12e4
stru cable begin g 43,62 end g 47,60 prop 1 t=12e4
stru cable begin g 44,66 end g 47,62 prop 1 t=12e4
stru cable begin g 46,68 end g 48,63 prop 1 t=12e4
stru cable begin g 48,69 end g 49,63 prop 1 t=12e4
stru cable begin g 51,63 end g 51,68 prop 1 t=12e4
stru cable begin g 53,62 end g 54,66 prop 1 t=12e4
stru cable begin g 54,60 end g 56,64 prop 1 t=12e4
stru cable begin g 55,59 end g 58,61 prop 1 t=12e4
stru cable begin g 55,57 end g 59,57 prop 1 t=12e4
stru cable begin g 55,55 end g 59,54 prop 1 t=12e4
stru cable begin g 55,52 end g 59,48 prop 1 t=12e4
;锚杆材料参数设定
stru prop=1 e=9.4e10 area=8.04e-4 kbond=2.e8 sbond=8.e4 yield=1.8e5
;锚索
stru cable begin 24,28.4 end 24,36 seg 19 prop 2
stru cable begin 26,28.4 end 26,36 seg 19 prop 2
;锚索材料参数设定
stru prop 2 e 2.0e11 yield 0.3e6 a 9.8125e-4 kbond 66.67e6 sbond 0.40e6
stru prop 2 sfric 20 perimeter 0.4082
;梁的设定与参数确定
stru beam begin g 47,55 end g 47,56 prop 3
stru beam begin g 47,56 end g 47,57 prop 3
stru beam begin g 47,57 end g 47,58 prop 3
stru beam begin g 47,58 end g 47,59 prop 3
stru beam begin g 47,59 end g 47,60 prop 3
stru beam begin g 47,60 end g 47,61 prop 3
stru beam begin g 47,61 end g 47,62 prop 3
stru beam begin g 47,62 end g 48,62 prop 3
stru beam begin g 48,62 end g 48,63 prop 3
stru beam begin g 48,63 end g 49,63 prop 3
stru beam begin g 49,63 end g 50,63 prop 3
stru beam begin g 50,63 end g 51,63 prop 3
stru beam begin g 51,63 end g 52,63 prop 3
stru beam begin g 52,63 end g 52,62 prop 3
stru beam begin g 52,62 end g 53,62 prop 3
stru beam begin g 53,62 end g 53,61 prop 3
stru beam begin g 53,61 end g 54,61 prop 3
stru beam begin g 54,61 end g 54,60 prop 3
stru beam begin g 54,60 end g 54,59 prop 3
stru beam begin g 54,59 end g 55,59 prop 3
stru beam begin g 55,59 end g 55,58 prop 3
stru beam begin g 55,58 end g 55,57 prop 3
stru beam begin g 55,57 end g 55,56 prop 3
stru beam begin g 55,56 end g 55,55 prop 3
stru beam begin g 55,55 end g 55,54 prop 3
stru beam begin g 55,54 end g 55,53 prop 3
stru beam begin g 55,53 end g 55,52 prop 3
stru prop=3 e=2.1e11 i=2.8e-5 area=1.5e-2
;历史值得记录
his unbal
his yd i=51,j=54
his yd i=51,j=63
his xd i=47,j=57
his xd i=55,j=54
;求解过程与计算时步
st 1000
save sjw-pm1c.sav
st 1000
save sjw-pm2c.sav
st 1000
save sjw-pm3c.sav
solve
save sjw-c.sav
return
计算结束 |
|