陈育民对FLAC3D常见问题的解答

目 录 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的固流耦合计算模式 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对孔压固定或释放。流体涌入或渗漏或可由命令apply施加。 渗流计算由命令set fluid和solve控制。如，set fluid on或off命令开启或关闭渗流计算模式。具体使用开启或关闭模式取决于渗流分析的耦合程度。 渗流分析结果以下面这些命令给出。命令print gp pp和print zone pp分别给出节点和单元孔压。节点和单元孔压历史可由命令histroy gp pp和命令history zone pp进行监测。对于瞬时计算，孔压与时间的关系可由命令history fltime监测。命令plot contour pp绘出节点孔压云图。命令plot contour saturation绘出饱和度云图。命令plot fluid绘出流量矢量图。渗流模式的所有信息由命令print fluid命令给出。FISH还提供了一些渗流变量。其中一个与节点有关的变量gp_flow，只能通过FISH函数使用。该变量描述了通过节点的净流入或流出量。因为可以提供一个系统总的流入或流出量，这些流量的统计在孔压固定的边界是很有用的。 渗流边界条件，初始条件 FLAC3D默认为不透水边界，即认为所有节点上的孔压随着从邻近单元流入或流出的量发生自由变化。可以使用命令fix pp将节点上的孔压设为“自由”，也可使用free pp使节点上的孔压“固定”。如果孔压固定，流体可以在外边界上流入或流出节点。下面总结这两种边界条件的影响： 1，孔压自由 这是默认的不透水边界条件。节点与外界之间不发生流量交换。系统根据当前饱和度值和流体是否形成涡凹现象来计算压力和饱和度变化。 2，孔压固定 这是一种流体通过外界流入或流出的边界条件。如果设定孔压为0，饱和度才可能变化。否则，饱和度被设为1（FLAC3D假设孔压只在完全饱和材料中存在）。孔压不能被固定在低于拉力极限的值，如果出现这种情况，FLAC3D会将其设定到拉力极限值。 如前所述，边界条件不是任意的。FLAC3D在进行计算前会“检查”并“修正”这些条件。可使用fix pp命令将孔压固定在某个值，也可在外边界或内边界上使用命令apply pp。如果边界条件被用于一个非表面节点，则必须加关键字interior。apply命令具有可以用“历史”命令进行监测的优点。 渗流边界条件可以通过apply命令用在单个或部分节点、单元面或单元上。命令apply pwell为边界节点指定了一个流入或流出井。如果加上interior关键字，则该条件用于内部节点。命令apply discharge和apply leakage为边界单元的表面分别指定了涌出和渗漏边界条件。命令apply vwell为指定区域内的单元提供一个流速。这些边界条件除了apply leakage外，均可使用history监测命令。 具有固定孔压节点就好像是流入源或流出源。没有直接的命令显示这些节点的流入或流出量。但可通过FISH变量gp_flow来记录。 孔压的初始分布，孔隙率，饱和度和流体属性可通过命令initial或property施加。如果还加了重力，则孔压初始分布应与重力梯度，水的重度和节点饱和度和孔隙率相容。如果这些初始分布不相容，则计算开始时所有单元中将出现流体流动。因此，应在模拟开始时设一定的计算步来检验初始条件是否相容。 如果模型中含有接触面，有效应力将沿着这些接触面进行初始化（即：在节点应力初始化时，认为接触面应力包含孔压）。water table命令将包含沿着接触面的孔压，这是因为定义在单元节点上的孔压也在接触面节点上。如果接触面的上下两面连在一起，在没有阻力时，将发生穿越接触面的流体流动。但程序不对沿着接触面的流体流动（裂隙流）进行计算。 单渗流与渗流耦合问题 FLAC3D既能进行单渗流分析，也能进行固流耦合分析。耦合分析可由FLAC3D内置力学模型完成。但要注意，渗流模型中的空单元并不是力学空单元。必须用命令model fl_null给单元赋予流体空属性。 对于耦合过程，FLAC3D提供了几种计算模式。其中之一是假设孔压一旦被赋予便不再改变。该方法并不要求任何额外空间存储计算过程。除此之外涉及到渗流的计算模式都要求使用命令config fluid。命令model fl_iso使所有单元中都能发生渗流。