研一学了flac3d，断断续续的学了好久，但是一直好像都没入门。一年过去了，现在又重新捡起flac，可是为什么爱你就这么难呢？
正题：
new
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;生成网格
generate zone brick size  5 5 5  p0 0 0 37.5 p1 add 100 0 0 p2 add 0 50 0  p3 add 0 0 27.6 p4 add 100 50 0  p6 add 100 0 27.6 p5 add 0 50 27.6 p7 add 100 50 27.6 group 粗粒砂岩
generate zone brick size 25 25 5 p0 0 0 9 p1 100 0 9 p2 0 50 9 p3 0 0 37.5 p4 100 50 9 p6 100 0 37.5 p5 0 50 37.5 p7 100 50 37.5  group 细粒砂岩
generate zone brick size 50 25 9 p0 0 0 0 p1 100 0 0 p2 0 50 0 p3 0 0 9 p4 100 50 0 p6 100 0 9 p5 0 50 9  p7 100 50 9 group 煤
generate zone brick size 25 25 2 p0 0 0 -2.5 p1 100 0 -2.5 p2 0 50 -2.5 p3 0 0 0 p4 100 50 -2.5  p6 100 0 0  p5 0 50 0  p7 100 50 0  group 泥岩
generate zone brick size 25 25 3 p0 0 0 -6.7 p1 100 0 -6.7 p2 0 50 -6.7 p3 0 0 -2.5 p4 100 50 -6.7 p6 100 0 -2.5 p5 0 50 -2.5 p7 100 50 -2.5  group 细粒砂岩

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;参数值
model moh
prop shear 1.92e10 bulk 3.2e10 coh 7.5e4 tension 2.5e6 friction 33 range group 粗粒砂岩
prop shear 2e10 bulk 3.33e10 coh 8e4 tension 2.6e6 friction 33 range group 细粒砂岩
prop shear 0.49 bulk 1.27e10 coh 2e4 tension 1.2e6 friction 26 range group 煤
prop shear 1.35e10 bulk 1.92e10 coh 2.5e4 tension 1.5e6 friction 28 range group 泥岩

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;初始化
set grav 0 0 -9.81

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;原岩应力
ini sxx -16.65e6 grad 0 0 2.5e4 ran x 0 100
ini syy -16.65e6 grad 0 0 2.5e4 ran y 0 50
ini szz -16.65e6 grad 0 0 2.5e4 ran z -6.7 65.1
app szz -16.65e6 ran z 56.1

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;边界条件
fix x range x -.1 .1
fix x range y 49.9 50.1
fix y range y -.1 .1
fix y range x 99.9 100.1
fix z range z -6.6  -6.8
set mech ratio 1e-4
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;连接面
attach face range z 37.4 37.6
attach face range z 8.9 9.1
attach face range z -.1 .1
attach face range z -2.6 -2.4
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;平衡
solve
save PINGHENG.sav

平衡结果为什么是这样的呢？好纠结啊？？是哪里错了？？

诶 兄弟看你和我专业这么相近我帮你看看吧
1.生成的网格乱七八糟的。你每个层位接触面对应的网格数目都不一样，网格都不连续，怎么attachface。
generate zone brick size  25 25 5  p0 0 0 37.5 p1 add 100 0 0 p2 add 0 50 0  p3 add 0 0 27.6 p4 add 100 50 0  p6 add 100 0 27.6 p5 add 0 50 27.6 p7 add 100 50 27.6 group 粗粒砂岩
generate zone brick size 25 25 5 p0 0 0 9 p1 100 0 9 p2 0 50 9 p3 0 0 37.5 p4 100 50 9 p6 100 0 37.5 p5 0 50 37.5 p7 100 50 37.5  group 细粒砂岩
generate zone brick size 25 25 9 p0 0 0 0 p1 100 0 0 p2 0 50 0 p3 0 0 9 p4 100 50 0 p6 100 0 9 p5 0 50 9  p7 100 50 9 group 煤
generate zone brick size 25 25 2 p0 0 0 -2.5 p1 100 0 -2.5 p2 0 50 -2.5 p3 0 0 0 p4 100 50 -2.5  p6 100 0 0  p5 0 50 0  p7 100 50 0  group 泥岩
generate zone brick size 25 25 3 p0 0 0 -6.7 p1 100 0 -6.7 p2 0 50 -6.7 p3 0 0 -2.5 p4 100 50 -6.7 p6 100 0 -2.5 p5 0 50 -2.5 p7 100 50 -2.5  group 细粒砂岩
attach face
2.煤的抗剪强度0.49？不对吧，我做实验也没这么小的
3.边界约束，这个真没看懂。是不是乱加的啊，一般都是地面全约束，四周用滚筒，也就是u*n=0,u为位移，n为边界法向：
fix x range x -.01 .01
fix y range y 49.99 50.01
fix y range y -.01 .01
fix x range x 99.9 100.1
fix x y z range z -6.6  -6.8
set mech ratio 1e-4

写错了 是模量。

感觉网格没对齐，看着云图显示都不连续

给你吧，flac3d真是让我伤心。

对应位置 改成我给你发的那个 就不伤心了