ying0824 发表于 2007-6-3 20:25:29

巷道模拟

;计算范围为50*50M 区域，划分单元10000个
grid 100,100
;设置材料模型
model mohr-coulomb
;生成网格
;直接底
gen   0,00,850,2350,0i=1,101 j=1,25 rat=1,0.9
;煤层
gen   0,80,1650,3150,23i=1,101 j=25,49 rat=1,1
gen   0,16 0,1750,3250,31i=1,101 j=49,52 rat=1,1
gen   0,170,2650,4150,32 i=1,101 j=52,76 rat=1,1
;直接顶/老顶
gen   0,26   0,5050,50 50,41i=1,101 j=76,101 rat=1,1.1

genarc 25,26.4 27.2,26.4 182
genline 22.8,2522.8,26.4
genline 22.8,2527.2,25
genline 27.2,2527.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
;锚杆材料参数设定
struprop=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-5area=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         
计算结束

wgy0207 发表于 2007-11-16 01:35:59

下了，谢谢楼主，试验一下，我也在做巷道的模拟，我的求求号是112316942

bbs_bbs1000 发表于 2007-11-16 14:55:35

强啊佩服

gaosongsong 发表于 2010-1-9 16:33:37

佩服！楼主是高手。

fan123000 发表于 2010-2-1 18:44:51

楼主高手啊

我是新手0554 发表于 2010-2-7 10:36:23

楼主太好了

w123wei 发表于 2010-4-17 18:26:26

楼主是用FLAC直接建几何模型吧，如果用ANSYS来建再导入FLAC会省力气一些。本人FLAC学习ing、、、

phymos 发表于 2010-4-17 20:45:57

楼主很强大请教一下梁式干吗的？

ajin 发表于 2010-4-17 22:38:17

谢谢分享！！！！！

wsongxue0105 发表于 2010-4-18 12:11:05

佩服！楼主是高手。

Visual1988 发表于 2011-3-11 16:47:44

救人于危难之中

syi450643708 发表于 2013-8-9 10:41:05

有没有锚注的模拟啊！！！楼主！！

程广坦 发表于 2013-9-4 10:38:49

谢谢楼主啊

skdwrc 发表于 2013-9-4 20:37:50

xiexie谢谢啦

xy4254 发表于 2013-9-28 20:38:10

先谢谢分享了

qq846765019 发表于 2013-10-9 15:58:00

楼主是好人啊
页: [1]
查看完整版本: 巷道模拟