论坛 › G0201：FLAC3D/FLAC2D › 初始平衡

shasha411 发表于 2012-3-17 00:22:09

初始平衡

做了一个模拟模型，只是加上了一些初始应力，可是显示应力的时候不是均匀分层，附上命令，各位大师给检查一下，问题是出在哪里了？
初始网格

new
plot block state
plot show
set log on

;模型尺寸：200m×120m×68m
;老底，细砂岩，8m
gen zone brick size 50302 &
p0   0 0 0   p1 200 0 0 p2 0 120 0   p3 0 0 8 &
group ld8

;直接底，砂质泥岩，6m
gen zone brick size 100603 &
p0   0 0 8p1 200 0 8p2 0 120 8   p3 0 0 14&
group zjd6

;3上煤层，5m
gen zone brick size 100 605 &
p0   0 0 14p1 200 0 14p2 0 120 14   p3 0 0 19&
group mc5

;直接顶，砂质泥岩，6m
gen zone brick size 100603 &
p0   0 0 19 p1 200 0 19p2 0 120 19   p3 0 0 25&
group zjd6-1

;直接顶，中粒砂岩，6m
gen zone brick size 100603 &
p0   0 0 25p1 200 0 25p2 0 120 25   p3 0 0 31 &
group zjd6-2

;直接顶，砂质泥岩，6m
gen zone brick size 100603 &
p0   0 0 31p1 200 0 31p2 0 120 31   p3 0 0 37&
group zjd6-3

;基本顶，细粒砂岩，8m
gen zone brick size 50 302 &
p0   0 0 37p1 200 0 37p2 0 120 37   p3 0 0 45&
group jbd8

;砂质泥岩，8m
gen zone brick size 50302 &
p0   0 0 45p1 200 0 45p2 0 120 45   p3 0 0 53&
group sfyc8

attach face all

初始平衡

plot block state
plot show

;定义材料模型及参数
model mohr

;细砂岩
prop shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500 range group ld8

;砂质泥岩
prop shear 0.5e9 bulk 0.7e9 cohesion 1.5e6 fric 28 tension 1.2e6 density 2500 range group zjd6

;煤层
prop shear 0.35e9 bulk 0.5e9 cohesion 1.3e6 fric 28 tension 0.9e6 density 1350 range group mc5

;砂质泥岩
prop shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500 range group zjd6-1

;中粒砂岩
prop shear 1.5e9 bulk 1.8e9 cohesion 3.5e6 fric 30 tension 2e6 density 2500 range group zjd6-2

;砂质泥岩
prop shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500 range group zjd6-3

;细粒砂岩
prop shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500 range group jbd8

;砂质泥岩
prop shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500 range group sfyc8

;边界条件
fix z range z-0.1 0.1
fix x range x-0.1 0.1
fix x range x199.9 200.1
fix y range y-0.1 0.1
fix y range y119.9 120.1

set grav 0,0,-9.81


;在上面z=53处施加力
apply szz -5e6 range z 52.9 53.1

;初始条件
ini szz -5e6
ini sxx -3e6
ini syy -3e6

hist unbal
plot hist 1

solve

shasha411 发表于 2012-3-17 00:25:53

初始平衡之后就是这样了，不是应该有规律的分层现象吗？这是出了什么问题啊，会的帮忙解答一下呗，谢了……

whutzll 发表于 2012-3-17 09:27:21

shasha411 发表于 2012-3-17 00:25 static/image/common/back.gif
初始平衡之后就是这样了，不是应该有规律的分层现象吗？这是出了什么问题啊，会的帮忙解答一下呗，谢了… ...

一、用Attach face命令检查模型网格，二、Solve elastic 试算初始应力场
试试看

411shasha 发表于 2012-3-17 10:26:12

whutzll 发表于 2012-3-17 09:27 static/image/common/back.gif
一、用Attach face命令检查模型网格，二、Solve elastic 试算初始应力场
试试看 ...

感觉是建的网格不太对，可是又看不出问题在哪里，很纠结啊……

411shasha 发表于 2012-3-17 11:14:53

这到底是什么问题啊，大师们帮忙检查一下网格建立过程有什么错误啊，这块问题不解决下边的也进行不动，着急啊

还有一个问题就是显示任何一个视图的时候总是比模型小一块，问题到底是在哪里呢？

aaclee 发表于 2012-3-17 18:37:57

估计是网格有问题吧

hclayor 发表于 2012-3-17 23:18:15

;初始条件
ini szz -5e6
ini sxx -3e6
ini syy -3e6
这段话会起到什么作用？

hclayor 发表于 2012-3-17 23:18:49

楼主查查塑性区分布，也许初始状态下就出现塑性破坏了

skzcg 发表于 2012-5-10 11:02:35

请问楼主问题是否解决？

chengcai18 发表于 2013-6-10 19:34:55

搞不懂aplly和ini哪位懂加我394649671 交流一下吧

ancel 发表于 2013-6-10 22:21:43

要是不成水平分布的话，你可以先将tens和coh设置的很大，然后so，再将tens和coh设置成真实的值，在so
你这个可能是网格布均匀的原因

查看完整版本: 初始平衡