flac学习总结（3）

1、  得到初始应力的方法：

方法1、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。如在手册中给的第一个示例中就是这样做的。下面是例子，These are only initial values that are used during the development of gravitational stresses within the body. In effect, we are forcing the body to behave elastically during the development of the initial in-situ stress state.* This prevents any plastic yield during the initial loading phase of the analysis.

Gen zone brick size 6 8 8

Mode mohr

Prop bulk 1e8 shear 0.3e8 fric 35

Prop cohesion 1e10 tens 1e10 ;注意在此这个值给的很大。

Init dens 1000

Set gravity 0 0 -10

Fix x range x -0.1 0.1

Fix x range x 5.9 6.1

Fix y range y -0.1 0.1

Fix y range y 7.9 8.1

Fix z range z -0.1 0.1

set mech force=50

solve

;---------------------- mode null ---------------------

Prop coh 1e3 tens 1e3 ;改为正常值(在此例中我们故意给小值)

Mode null range x 2 ,4 y 2 , 6 z 5, 10

Set large

Ini xdis 0 ydis 0 zdis 0 ;清零，不影响结果，为画图方便。

Step 2000

return

方法2：得到初始应力的一个不错办法是设为弹性模型，然后赋予地勘信息，so以后就ok了，得到初始应力了，再把相应的位移清零，下面就可以正式进入所分析的内容了。如上例还可以先设为弹性 (MODEL elastic)，之后在开挖之前再设为 Mohr-Coulomb model， 就行。

 

2、  创建Interface单元的方法

在FLAC3D中，有两种办法来创建INTERFACE单元，一种就是在实体单元的自由表面上通过指定范围来实现，一种就是通过先创建界面节点，再创建界面单元，然后在在节点和实体单元间创建联结关系，我想说的是第二种办法很受局限性，因为软件将通过这种方式创建的界面NODE AND ELEMENT都固定在创建处。

最大不平衡力：每个节点周围都绕着8个立方体所集中的应力，在平衡状态下，这些应力的代数和应该为0。并不一定在所有的情况下都为0，它只是评估是否达到平衡的一个方法，一般说来并不真正的为0，只是相比之下很小。最大不平衡力为一固定值时，说明会有两种情况发生，一种是达到了平衡，另一种是发生了塑性破坏或蠕动破坏。

3、  Tension-p shear-p tension-n shear-n表示的含义（在plot block state命令显示的图中）：在计算的循环里面,每个循环中,每个zone都依据failure criterion处于不同的状态,shear 和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态.

看了一个例子，算边坡的安全系数，但我不知道用什么命令查看我最终的计算结果呢，比如我想知道我算得的安全系数到底是多少

/网格省去

model mohr

prop bulk = 1.0E8 shear = 3.0E7 coh 12380.0 tens 1.0E10 &

friction = 20 dilation = 20

fix x y z range z -0.1 0.1

fix x range x 19.9 20.1

fix x range x -0.1 0.1

fix y

ini density=2000

set gravity = 10.0

set large

solve fos file slope3dfos.sav associated

ret

用什么命令查看最终的安全系数呢？

Plot fos