陈育民对FLAC3D常见问题的解答重点.docVIP

目 录1．FLAC3D的固流耦合计算模式 1 2．FLAC3D固流耦合学习小结 5 3．关于流固耦合的问题 6 4．也谈采用FLAC3D对地下采矿的模拟 8 5．FLAC3D本构模型开发 8 6．FLAC3D自定义本构模型 11 7．数值计算中初始应力场的模拟 13 8．FLAC3D应变分析 13 9．FLAC3D的调参 14 10．开采沉陷垂直剖面等值线的生成 15 11．FLAC3D的应变硬化软化模型 16 12．FLAC3D的塑性流动格式 17 13．FLAC3D的动画制作 17 14．地下连续墙基坑开挖支护 18 15．一个汇的小例子 21 16．用3DEC生成岩体随机节理网络 23 17．固结小算例 24 FLAC3D的固流耦合计算模式英文原文 26 1．FLAC3D的固流耦合计算模式 /blog/static/323428402007102243049387/ FLAC3D的计算模式中是否需要做孔压分析取决于是否采用config fluid命令。 1 无渗流模式（不使用config fluid） 即使不使用命令config fluid，仍然可以在节点上施加孔压。这种模式下，孔压将保持为常量。如果采用塑性本构模型的话，材料的破坏将由有效应力状态来控制。 节点上的孔压分布可由initial pp命令或water table命令来设定。 如果采用water table命令，由程序自动计算水位线以下的静水孔压分布。此时，必须施加流体密度（water density）和重力（set gravity）。流体密度值和水位位置可以用命令print water显示。如果水位线是由face关键字来定义的，则可用命令plot water命令显示水位。 这两种情况，单元的孔压都由节点孔压值平均求出，并在本构模型计算中用作有效应力。这种计算模式下，体积力中不反映流体的出现：用户必须根据水位线以上或以下相应地指定干密度和湿密度。使用命令print gp pp和print zone pp可分别得到节点或单元孔压。plot contour pp命令可绘出节点孔压云图。 2 渗流模式（使用config fluid） 如果使用命令config fluid，则可进行瞬时渗流分析，孔压改变和潜水面的改变都可能出现。在config fluid模式下，有效应力计算（静态孔压分布）和非排水计算均被执行。除此之外，还可进行全耦合分析，这种情况下，孔压改变将使固体产生变形，同时体积应变反过来影响孔压的变化。 如果采用渗流模式，单元孔压仍由节点孔压平均求出。但这种模式，用户只能指定干密度（不论是水位以上还是以下），因为FLAC3D将流体的影响考虑到了体积力的计算中。 采用渗流模式时，渗流模型必须施加到单元上，使用命令model fl isotropic模拟各向同性渗流，model fl anisotropic模拟各向异性渗流，model fl null模拟非渗透物质。注意，力学模型为空的单元并不代表渗流模型为空。 流体性质（参数）可施加到单元或节点上。各向同性渗透率、孔隙率、比奥系数和非排水热系数等单元流体性质由命令property施加。 对于各向同性渗流，渗透率通过perm关键字赋予。对各向异性渗流，渗透率的3个主值采用关键字k1,k2,k3赋予，主方向由关键字fdip,fdd,frot确定。渗透率的主方向服从右手系统。fdip和fdd分别为k1和k2确定的平面的倾向和倾角。frot为k1轴和倾角矢量的旋转角。如果不特别指定，比奥系数默认为1，孔隙率默认为0.5。节点的渗流性质由命令initial指定。这些性质包括流体重度、流体体积模量、比奥模量、流体抗拉强度和饱和度。每种性质在空间上都可以变化。流体重度也可以用water命令给出。 在渗流模式里，有必要知道可压缩性被定义在以下两种参数中：（1）比奥系数和比奥模量；（2）流体体积模量和孔隙率。第一种参数表征的是固体颗粒的可压缩性（对不可压缩颗粒，比奥系数设为1）。对第二种参数，固体颗粒被认为是不可压缩的。 单元属性可由命令print zone property显示，节点属性由print gp 命令显示。流体重度，如果随着水位位置被确定，则可由print water命令显示。渗流性质可由命令plot bcontour property显示。对于各向异性渗流，渗透率的各球形分量可通过使用单元的属性关键字kxx,kyy,kzz,kxy,kxz,kyz来显示（注意，这些球形分量不可被直接初始化）。 初始节点孔压分布的施加对于渗流模式和非渗流模式都是一样的（如，要么用initial pp命令或用water table命令）。在指定节点可用命令fix pp或free pp对孔压固定