Flac3d是有限差分数值分析软件，它可以模拟岩土和其它材料的三维力学特性。因此，在岩土工程力学中得到了广泛的应用，我们可以用该程序较好的模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学特性，特别适用于分析渐进破坏失稳以及模拟大变形。该软件中提供的强大丰富的本构模型可以进行静力、动力、蠕变、渗流、温度等多种计算模式，各种模式间可以互相耦合，也可以模拟多种结构形式，如梁、锚杆、土工格栅、衬砌、桩等。作为一位初学者，主要了解了下flac3d中的弹性模型和莫尔库仑模型，通过一段时间的练习主要有以下心得；在分析工程案列时，我们可以依照flac3d的一般求解流程进行剖析案列，将命令流归纳为三个部分，即分析模型、模拟及求解、计算结果输出。分析模型部分包括生成网格单元、设置初始条件和边界条件以及初始应力平衡；在模拟求解部分包括加载及连续建模、后续工程施工、场方程的有限差分求解；输出计算结果部分主要为图表的绘制、相关数据的输出等；通过这三个部分，我们可以写出一个完整的命令流，方便后续在运行过程中修改命令流以及查出具体在哪个环节出错了。这里我们应该注意在flac3d的建立分析模型部分，材料性质的定义、初始条件和边界条件的设置并无明显的先后顺序。

对于工程实例中，我们在分析模型中，首先看是否对称，如果对称我们可以简化模型，然后用gen zone reflect dip 关键字 dd 关键字 ori 关键字 这个命令来镜像模型。这里我们需要注意的就是dip和dd的具体含义才能把模型建好；dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线 ，dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．

gen zone reflect dip 90 dd 90是以YZ面.

gen zone reflect dip 90 dd 0 才是以xz为对称面的

当初始应力平衡后，我们可以求工程的影响，求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤： 

(1)调出起初的应力平衡：re_____ .sav 

(2)设定新的材料性质：model,prop 

(3)设定新的支撑性质：struct

(4)设定新的外力 apply

(5)设定边界条件 apply,fix

(6)求工程时的应力平衡 solve, step, cycle

(7)储存 save file.sav

在初始应力平衡过程中，我们要注意是否收敛，对于不收敛一般要检查模型是不是没有建好，相互之间是不是没有连接好；仔细检查参数 ，是不是已经发生了破坏；对于计算时间过长，我们可以适当修改模型的网格数以减少计算时间，或者使用step命令，这里的关键问题是要运行多少步。对于边坡稳定性的判断方法之一为最大不平衡力的收敛性，最大不平衡力的定义为：在每一个计算循环中，外力通过网格节点传递分配到体系各节点时，所有节点的外力与内力之差中的最大值。最大不平衡小于某一临界值作为边坡的收敛标准。命令流为：set mech force 。

在给材料参数初始化时，我们要注意initial、apply这两者的区别。initial初始化命令，如初始化计算体的应力状态等；" W. Z% C( ~- \\+ r' r7 ~1 U; `3 Papply边界条件命令，如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。

Flac3d提供了强大的后处理，我们可以用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了，也可以得到直观的图形，如plot block group网格划分图，plot vel red grid速度矢量云图，plot con disp 位移等值线图，plot con vel 速度等值线图。下载本文