flac3d 自我总结，和交流 [复制链接]

最近遇到瓶颈了，不知道什么学了，我现在也有几个问题，一个就是地应力的是，到底初始应力是什么来的，前后模型参数不一样对初始应力场有什么影响等

这两天在仿真论坛上逛，吧初始应力给看了一遍，自己也把不同老大的东西整整。
结论如下

1.设定初始地应力场是十分必要的 没有地应力场进行计算会使模型在重力作用下发生过大的变形 从而不能与实际情况相符

2。apply是施加的外力  ini是初始的内力。

3. 有人问如果改参数这样做的话地应力不是被人为改变了吗？为什么要人为加大它的参数呢，很不合理？
回答：不能改变bulk和shear模量，因为这两个参数是用来计算应力的。而con和fric是用来判断破坏状态的。如果将bulk和shear模量改变了，那你在重力条件下计算的应力就不真实了。咱们改大参数的目的是让模型在不破坏的条件下得到正确的应力，所以只有将con和fric变大才能达到目的。
        当然如果在弹性模型下计算初始应力就没有这种问题了。因为弹性模型是没有破坏准侧的当然也没有描述达到破坏状态的参数了。

4. 大家正好讨论地应力的问题，我也有个疑问，平衡地应力后，水平方向的地应力如何确定的？是由本构模型来确定的吗？
这个问题现在我也没答案也是我的问题

5. gen zone brick size 10,10,10 p1 20,0,0 p2 0,20,0 p3 0,0,20
model mohr
prop bulk 5e9 shear 3e9 fric 35
ini density 2500  
set gravity 0,0,-10
fix x y z range z -0.1 0.1
ini szz = -5.0e6 grad 0,0,2.5e4
ini sxx = -2.5e6 grad 0,0,1.25e4
ini syy = -2.5e6 grad 0,0,1.25e
apply szz = -4.5e6 range z 19.9 20.1
apply szz = -5.0e6 range z -0.1 0.17
apply sxx = -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1
apply sxx = -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1
apply syy = -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1
apply syy = -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1
  U# P! y- P$ _+ D
它的意思是通过ini szz来赋予初始地应力,但是为何又有set gravity 0,0,-10
这样不是重复了吗? ini是的确是模型的内力，这个内力通过set grav ,然后就像一些帖子上介绍的那样,让程序计算不就行了吗?干嘛还有手动在设定一次呢?
我感觉它这样的做法有点重复了.

回答：（1）它这样做其实只是减少程序的计算步骤.5 _这样设置以后 整个模型就已经是平衡了
（2）设置重力的原因不是为了平衡减少程序计算步骤，而是因为如果未初始化应力的话，在set grav的时候，模型会在自重应力下发生应力变形。所以，初始化应力的原因是防止模型在自重应力下发生变形（我们认为地层自重变形已经发生，当然欠固结等现象是另外一回事）。

（3）其实就让模型发生变形，然后把位移清零也可以的

6。ini命令来对地层进行初始应力的施加，而APPLY是~7 u(施加应力边界条件
7. 用到ini初始平衡的时候要和重力计算的相同才有用

8。
gen zone brick size 10 10 10 &
p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20
model elas
prop bulk 5e9 shear 3e9
ini density 1600 range y -10,0
ini density 2000 range y -15,-10
ini density 2200 range y -25,-15
set gravity 0,-10
fix x range x -.1 .1
fix x range x 19.9 20.1
fix z range z -.1 .1
fix z range z 19.9 20.1
fix y range y -25.1 -24.9
solve
 
如果这样是不是就完成了初应力的计算了？
当然在这个基础上进一步计算需初始位移为0:

上面传的时候有点乱码，大家自动跳过