- 积分
- 0
- 注册时间
- 2016-4-12
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-4-24 10:52:36
|
显示全部楼层
来自 河南焦作
new
pl bl gr
;各岩层、模块建立
gen zone brick p0 0 0 0 p1 64 0 0 p2 0 4 0 p3 0 0 8 size 32 4 4 group 粉砂岩
gen zone brick p0 64 0 0 p1 96 0 0 p2 64 4 0 p3 64 0 8 size 32 4 4 group 粉砂岩
gen zone brick p0 96 0 0 p1 160 0 0 p2 96 4 0 p3 96 0 8 size 32 4 4 group 粉砂岩
gen zone brick p0 0 0 8 p1 64 0 8 p2 0 4 8 p3 0 0 14 size 32 4 3 group 中粒砂岩
gen zone brick p0 64 0 8 p1 96 0 8 p2 64 4 8 p3 64 0 14 size 32 4 3 group 中粒砂岩
gen zone brick p0 96 0 8 p1 160 0 8 p2 96 4 8 p3 96 0 14 size 32 4 3 group 中粒砂岩
gen zone brick p0 0 0 14 p1 64 0 14 p2 0 4 14 p3 0 0 20 size 32 4 3 group 细粒砂岩
gen zone brick p0 64 0 14 p1 96 0 14 p2 64 4 14 p3 64 0 20 size 32 4 3 group 细粒砂岩
gen zone brick p0 96 0 14 p1 160 0 14 p2 96 4 14 p3 96 0 20 size 32 4 3 group 细粒砂岩
gen zone brick p0 0 0 20 p1 64 0 20 p2 0 4 20 p3 0 0 24 size 32 4 2 group 泥岩
gen zone brick p0 64 0 20 p1 96 0 20 p2 64 4 20 p3 64 0 24 size 32 4 2 group 泥岩
gen zone brick p0 96 0 20 p1 160 0 20 p2 96 4 20 p3 96 0 24 size 32 4 2 group 泥岩
gen zone brick p0 0 0 24 p1 64 0 24 p2 0 4 24 p3 0 0 31 size 32 2 7 group 煤层
gen zone brick p0 64 0 24 p1 75 0 24 p2 64 4 24 p3 64 0 31 size 11 2 7 group 煤层
gen zone radcylinder p0 80 0 26 p1 85 0 26 p2 80 4 26 p3 80 0 31 dim 2 2 2 2 size 2 2 10 3 ratio 1 1 1 1.2 group 1
gen zone reflect normal (-1,0,0) origin=(80 0 26) range group 1
gen zone cylinder p0 80 0 26 p1 82 0 26 p2 80 4 26 p3 80 0 28 dim 2 2 2 size 2 2 10 group 2
gen zone reflect normal (-1,0,0) origin=(80 ,0,26) range group 2
gen zone brick p0 75 0 24 p1 85 0 24 p2 75 4 24 p3 75 0 26 size 10 2 2 group 煤层
gen zone brick p0 85 0 24 p1 88 0 24 p2 85 4 24 p3 85 0 31 size 3 2 7 group 煤层
gen zone brick p0 89 0 25 p1 97 0 25 p2 89 4 25 p3 89 0 29 size 8 2 4 group 煤层B
gen zone brick p0 97 0 25 p1 161 0 25 p2 97 4 25 p3 97 0 29 size 32 2 4 group 煤层B
interface 1 face range plane norm (0,0,24) origin (88,0,24) dist 0.1 x 92 160
ini x add -1 range group 煤层B
ini z add -1 range group 煤层B
gen zone brick p0 89 0 29 p1 97 0 29 p2 89 4 29 p3 89 0 32 size 8 2 3 group 煤层C
gen zone brick p0 97 0 29 p1 161 0 29 p2 97 4 29 p3 97 0 32 size 32 2 3 group 煤层C
interface 2 face range plane norm (0,0,28) origin (88,0,28) dist 0.1 x 92 160
ini x add -1 range group 煤层C
ini z add -1 range group 煤层C
gen zone brick p0 0 0 31 p1 64 0 31 p2 0 4 31 p3 0 0 35 size 32 4 2 group 砂质泥岩
gen zone brick p0 64 0 31 p1 96 0 31 p2 64 4 31 p3 64 0 35 size 32 4 2 group 砂质泥岩
gen zone brick p0 96 0 31 p1 160 0 31 p2 96 4 31 p3 96 0 35 size 32 4 2 group 砂质泥岩
gen zone brick p0 0 0 35 p1 64 0 35 p2 0 4 35 p3 0 0 45 size 32 4 5 group 细粒砂岩
gen zone brick p0 64 0 35 p1 96 0 35 p2 64 4 35 p3 64 0 45 size 32 4 5 group 细粒砂岩
gen zone brick p0 96 0 35 p1 160 0 35 p2 96 4 35 p3 96 0 45 size 32 4 5 group 细粒砂岩
gen zone brick p0 0 0 45 p1 64 0 45 p2 0 4 45 p3 0 0 65 size 32 4 10 group 基岩
gen zone brick p0 64 0 45 p1 96 0 45 p2 64 4 45 p3 64 0 65 size 32 4 10 group 基岩
gen zone brick p0 96 0 45 p1 160 0 45 p2 96 4 45 p3 96 0 65 size 32 4 10 group 基岩
;链接各点
gen merge 1e-6 range z 0 24 x 0 160 y 0 4
gen merge 1e-6 range z 24 26 x 0 92 y 0 4
gen merge 1e-6 range z 26 31 x 0 92 y 0 4
gen merge 1e-6 range z 28 31 x 92 160 y 0 4
gen merge 1e-6 range z 31 65 x 0 160 y 0 4
;设置本构模型
model mohr
;各岩层赋参数
;基岩
prop dens 2689 bulk 8.86e9 she 4.177e9 coh 3.92e6 fric 32 ten 3.83e6 range group 基岩
;细粒砂岩
prop dens 2689 bulk 4.86e9 she 2e9 coh 2.92e6 fric 35 ten 3.83e6 range group 细粒砂岩
;砂质泥岩
prop dens 2506 bulk 3.71e9 she 1.85e9 coh 1.85e6 fric 28 ten 2.16e6 range group 砂质泥岩
;煤层
prop dens 1420 bulk 1.50e9 she 0.75e9 coh 0.88e6 fric 26 ten 0.98e6 range group 煤层
;煤层B
prop dens 1420 bulk 1.50e9 she 0.75e9 coh 0.88e6 fric 26 ten 0.98e6 range group 煤层B
;煤层C
prop dens 1420 bulk 1.50e9 she 0.75e9 coh 0.88e6 fric 26 ten 0.98e6 range group 煤层C
;泥岩
prop dens 2516 bulk 3.71e9 she 1.85e9 coh 1.85e6 fric 28 ten 1.85e6 range group 泥岩
;细粒砂岩
prop dens 2619 bulk 4.86e9 she 2.17e9 coh 1.92e6 fric 36 ten 3.83e6 range group 细粒砂岩
;中粒砂岩
prop dens 2635 bulk 5.35e9 she 2.69e9 coh 2.7e6 fric 36 ten 3.64e6 range group 中粒砂岩
;粉砂岩
prop dens 2732 bulk 5.41e9 she 2.72e9 coh 1.33e6 fric 33 ten 4.75e6 range group 粉砂岩
;group1
prop dens 1420 bulk 1.50e9 she 0.75e9 coh 0.88e6 fric 26 ten 0.98e6 range group 1
;group 2
prop dens 1420 bulk 1.50e9 she 0.75e9 coh 0.88e6 fric 26 ten 0.98e6 range group 2
interface 1 prop kn 13.68e9 ks 13.68e9 fric 20
interface 2 prop kn 13.68e9 ks 13.68e9 fric 20
;设置边界条件
fix z range z= (-0.1 0.1)
fix x range x= (-0.1 0.1)
fix x range x= (159.9 160.1)
fix y range y= (-0.1 0.1)
fix y range y= (3.9 4.1)
;初始条件
set gravity 0 , 0, -10
apply szz -7.425e6 range z 64.9 65.1
;设置水平和垂直应力
ini szz = -7.425e6
ini sxx = -8.167e6
ini syy = -8.167e6
ini xd 0 zd 0 yd 0
solve
save 10m煤柱加密初始应力.sav
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
mo nu ran x 92 160 y 0 4 z 24 28
set large
solve
save 10m煤柱加密采空区.sav
mo nu ran group 2
mo nu ran x 78 82 y 0 4 z 24 26
set large
solve
save 10m煤柱加密.sav
|
|