FLAC3D深部工程建模S-B法的讨论 [复制链接]

在书上看到S-B法的例子
原来的命令
new
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 &
size 6 6 10
model elas
pro bulk 10e10 she 10e10
ini den 2500
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
apply sxx -1e9 grad 0 0 7.6666666e6 range x 59.9 60.1
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120
set grav 0 0 -10
step 30000
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
plot cont szz



<DIV id=read_2045051 mb10?>请问高手例子中的一些参数是如何得到的
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
apply sxx -1e9 grad 0 0 7.6666666e6 range x 59.9 60.1
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120
中的1.11111111e7 ,7.6666666e6和8.333333333333e5 <DIV mb10?><DIV id=read_2045055 mb10?>还有最后apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120，为什么是range z 0 120 ,不是和之前一样取两个面呢

同样困惑，顶起，希望高手解答。

应力梯度线性变化，x y方向施加的是约束！

这个例子有问题吧？z方向上只在顶面上施加向下的力，模型底部却没有固定，运算后 这个模型应该会移动吧？？？？？？？？？？？？

有什么问题啊

假定对山体内某一深层岩体进行某岩石工程的开挖，为计算其稳定性，现对其周围一定范围内的岩石进行模拟。其初始条件为：模型的长×宽为60m×60m；垂直方向延伸到地表，地表为一山坡；从模型地面到地表平均深度为120m，最深处150m。根据地应力勘测资料，地应力场远远大于自重应力场，模型底面处sxx = syy =1000MPa，szz=100MPa，从底面到地表呈线性分布。假定岩体密度2500 kg/m&sup3;；同时，为保证弹性模型中岩体的变形很小，假定岩体的体积模量K =1.0×106 MPa，剪切模量G =1.0×106 MPa。
参数说的很清楚
apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120   ;设置z=0 边界z 向竖向应力
根据原文的解释 这句命令写错了 应该是 ran z -0.1 0.1
之所以取得一个面 因为上面那个面是地表 不会受力的