FLAC-3D 初始地应力 [复制链接]

一下是我对初始地应力模拟的程序。模型是30*30*30，水平荷载22e6，垂直荷载20e6
程序如下：
new
gen zon brick p0 0 0 0 p1 300 0 p2 0 30 0 p3 0 0 30size 30 30 30
model elas
prop bulk 1e10 shear 1e10
ini dens 2000
fix x ran x -0.1 0.1
fix x ran x 29.9 30.1
fix x y z ran z -0.1 0.1
fix y ran y 29.9 30.1
fix y ran y -0.1 0.1
;
set large
set grav 0 0 -10

;
ini sxx -22e6
ini syy -22e6

ini szz -20e6
;
apply szz -20e6 ran z 29.9 30.1

solve

以上程序哪里有不妥之处希望各位前辈指教？
谢谢

谁牛X能告诉以上程序错在哪了？？？？？？？？？？？？

p1 300 0只发现一个低级错误，别的不清楚

体积量与弹性模量的比值一般是泊松比的函数有个取值范围，所以在设定这两个参数值要注意，一般是根据弹模和泊松比来换算这两个量的

请问 上面这个程序有问题吗
是对的吗 ？

做的深埋，没加重力。结果还行。思想是边界应力平衡内部应力。仅供参考。

ini sxx -27.69e6

ini syy -29.67e6

ini szz -42.57e6

apply sxx -27.69e6 range x -30.01 -29.99
apply sxx -27.69e6 range x 39.99 40.01
apply syy -29.67e6 range y 29.99 30.01
apply syy -29.67e6 range y -25.01 -24.99
apply szz -42.57e6 range z -0.01 0.01
apply szz -42.57e6 range z -60.01 -59.99

这个如果加上重力得到的效果怎么样，你做过吗？得到的SZZ SXX SYY  是不是必须是层状的

def ini_para              
rock_dens = 2.8e3  
x_pp_b = -0.134e6+(264+206)*0.0182e6*(-1)
y_pp_b = -0.445e6+(264+206)*0.0395e6*(-1)
z_pp_t = -0.13e6+264*0.0106e6*(-1)
z_pp_b = -0.13e6+(264+206)*0.0106e6*(-1)
rock_e_mod = 60e9
rock_p_ratio = 0.2
rock_b_mod = rock_e_mod / (3.0 * (1.0 - 2.0 * rock_p_ratio))
rock_s_mod = rock_e_mod / (2.0 * (1.0 + rock_p_ratio))                                    
end
ini_para    
model ela
prop bulk rock_b_mod shear rock_s_mod
ini den rock_dens

ini sxx x_pp_b grad 0,0,0.0182e6
ini sxx x_pp_b grad 0,0,0.0182e6
ini syy y_pp_b grad 0,0,0.0395e6
ini syy y_pp_b grad 0,0,0.0395e6
ini szz z_pp_b grad 0,0,0.0106e6

app sxx x_pp_b grad 0,0,0.0182e6 range x -49.9 -50.1
app sxx x_pp_b grad 0,0,0.0182e6 range x 138.9 139.1
app syy y_pp_b grad 0,0,0.0395e6 range y 39.9 40.1
app syy y_pp_b grad 0,0,0.0395e6 range y -69.9 -70.1
app szz z_pp_b grad 0,0,0.0106e6  range z  0
;set gravity 0,0,-10
solve
麻烦帮我看看这个行吗？

如果加重力的话，直接在竖向用应力梯度模拟，而不使用重力场。
我觉得这样好理解一点。

谢谢啦，那我这个程序就没有问题吧