流体力学FLAC3D

1、1.流体力学一单向流动1.1.简介FLAC3D通过具有渗流性的实体比如土来模拟流体的流动。流动模型的建立可 以独立于力学计算而自动完成，或者说可以与力学模型同时建立，这样就可以考虑 流体与土体之间的相互作用。流固耦合的一种类型是“固结”，即：空隙水压力逐渐 消散而导致土体的沉降。这个过程包括两种力学反映：一，空隙水压的改变导致有 效应力的变化，这将影响到土体的力学反映（如：有效应力的减小可能导致塑性区 的产生）；二，力学实体中某一区域的流动会随着空隙水压的改变而改变。该程序可以计算完全饱和情况下的流动，也可以模拟具有自由水面的流动。 模拟具有自由水面的流动时，自由水面以上的部分空隙水压等于0,

2、气相将不参与计算。 对于不考虑毛细水压力颗粒较粗的材料可以采用这种模拟方法。流体计算就有以下特点：1根据各项同姓和各项异性的渗流计算，相应采用两种流体运动定律。流动中的null材料用来模拟流动范围内的非渗流材料。2不同区域可以拥有不同的流动模型（isotropic, anisotropic or null）和模型参 数。3可以事先指定流体的压力、流量、非渗流区边界条件。4流体源可以以电源，也可以以体源的形式插入到材料中，这些源对应于流体的流入或流出，可以随着时间而变化。5对于完全饱和流动，可以采用显式和隐式两种算法， 但对于非饱和流动则只 能米用显示计算。6任何力学和温度计算模型都可以与流体模

3、型一起使用，在耦合计算中，可以考虑饱和体的压缩性和热膨胀性。7.流体与力学计算的耦合通过提供比奥系数来实现。8 .与温度的耦合计算可以通过提供线性热膨胀系数at和不排水温度系数6（undrained thermal coefficient可能翻译的不对）来实现。9.热-流体计算以线性理论为基础，假定材料参数为常数，不考虑对流。流体与实体的温度保持局部平衡。非线性行为可以采用 FISH语言改变孔隙压力、材料特 性来实现。动力水压力的产生以及循环荷载引起的液化也可以用FLAC3D模拟（3.4.4节）cFLAC3D不考虑部分饱和体中的毛细水压、粒子之间的电、化学力。但是可以根据材 料的局部饱和情况，

4、孔隙率或其他相关变量利用FISH语言施加内部应力来实现上述 几种作用力。与此类似，由于气体溶解引起的流体刚度变化的影响也是不能直接模 拟的，但是FISH编写一个压力，时间的函数就可以用来改变局部流体的模量本意主要包括7个方面的内容：1. 1.2节和1.3节中分别介绍流体力学，流-力耦合以及流-热耦合的数学 理论模型和数学表达式。2. 1.4中介绍流体力学模拟采用的计算模型和相关的命令。3. 1.5节介绍流体力学模拟采用的材料参数，各参数的单位。4. 1.6介绍流体力学模型采用的不同的边界条件，初始条件，流源和汇集。5. 1.7介绍流体分析的一般过程，该部分还包括几个例子，建议读者在做 自己的流

5、体分析以前首先研究活楚我们提供的例子。6. 1.8是几个考核题，用来验证flac3d流体分析的准确性。7. 1.9总结了所有与流体计算相关的命令，以供参考。读者应该首先熟悉力学计算，然后才能将流体，力一起来考虑。流 -力耦合是 一个复杂的问题，需要读者有很好的洞察力来判断结果的正确性。（1.2 1.3我看就没有必要翻译了，有兴趣看看朱伯芳院士的书就可以了，理论是一样的）1.4流固耦合的计算模式孔隙水压的计算模型和命令取决于节点是否已经设置了流动（也就是说，CONFIGfluid命令是否已经采用）1.4.1和1.4.2将分别介绍这两种情况。方便起见: 下面所有的命令在文末1.9都作了汇总。1.4

6、.1 Grid Not Configured for Fluid Flow如果不给出CONFIG fluid命令，仍旧可以在节点上设置孔隙水压力。这种计算 模式下孔隙水压力不变，但是当材料为弹塑性体时有可能不收敛。使用initial pp命令加a gradient,或用 WATER table 命令可以指定节点上孔 隙水压力的分布。如果使用WATER table 命令，水面一下部分的孔隙静水压力的分布程序将自 动进行计算。这种计算模式，必须用 WATERensity设置流体密度,SETgravity 设 置重力加速度。当用水面定义好以后，采用 PRINT water命令可以输出流体密度， 水面

7、的位置，PLOT water命令可以显示出水面来两种情况下微元体的孔隙压力计算是将个节点的值平均，从而计算有效应力。 这种计算模式，体力的计算不考虑流体压力：水面以上和以下材料的天然重度和饱 和重度要有用户自己指定。PRINT gp pp和PRINT zone pp可以输出节点和体的孔 隙水压力。PLOT contour pp 命令用来绘制节点孔隙水压的云图。在用户手册中的3.6.3有这种计算模式的一个简单例子供读者学习。1.4.2 Grid Configured for Fluid Flow如果使用了 CONFIG fluid命令，就可以进行流体的瞬态分析。地下水面和孔隙 压力都会发生改变。

8、在CONFIG fluid计算模式下，有效应力（静止孔隙压力的分布） 和不排水计算分析可以同时考虑，此外，可以进行完全耦合的计算分析，这时孔隙 水压的改变会导致变形同时体积应变乂会引起孔隙水压的变化。微元体的孔隙水压仍旧等于其各个节点上孔隙水压的平均值，但是这种计算模式用户必须指定材料的十密度（包括水面以上和以下的部分），因为在 FLAC计算 体力的时候考虑了流动的影响。模式下必须指定单元体的流动模式：MODEL fl_ isotropic命令指定为各向同性流动模式MODEL fl _anisotropic .命令指定为各向异性流动模式MODEL fl_null命令指定为不透水区域，注意力学分

9、析指定的null区域不能自 动成为流体中的不透水区域。流体的性质可以赋予单元体和节点，单元体的流体性质可以用命令赋值，各向 同性体包括：渗透系数、孔隙率、比奥系数。对于各向同性的流动，渗透系数采用关键词perm指定。各项异性流动采用关键词：Ki, K 2, K 3指定三个主方向的渗透系数，用关键词fdip , fdd , frot .指定三个主方向。三个主方向的渗透系数与K 1, K2, K 3相对应且服从右手定则。其中fdip和fdd分别是k1 k2两方向组成的平面的倾角和倾向角，其中，倾角是以从 xy-面向z- 轴负方向转动为正。倾向角是倾向在xy-面同y轴正向的火角（顺时针）。 是k1-

10、k2 面内倾向顺时针转到k1方向的角度。图1 1三个渗透主方向的定义默认条件下，比奥系数等丁 1,孔隙率为0.5。节点的流体性质用INITIAL命令付 值。其性质包括：流体密度，流体体积模量，比奥模数，饱和度。节点的性质可以 给定变化范围，流体密度可以采用 WATER命令一次性付值。下面表1-1是相关参数 及其关键字：表1-1相关的参数的确定性质关键字范围命令渗透系数（各向同性）permin zonesPROPERTY主方向渗透系数（各项异性）k1in zonesPROPERTY主方向渗透系数（各项异性）K2in zonesPROPERTY主方向渗透系数（各项异性）K3in zonesPROP

11、ERTY渗透主方向定义（各项异性）Fdipin zonesPROPERTY渗透主方向定义（各项异性）Fddin zonesPROPERTY渗透主方向定义（各项异性）frotin zonesPROPERTY孔隙率porosin zonesPROPERTY比奥系数Biot_cin zonesPROPERTY流体模量fmodat gridpointsINITIAL比奥模量biot modat gridpointsINITIAL饱和度satat gridpointsINITIAL流体拉极限ftensat gridpointsINITIAL流体密度fdensat gridpointsINITIALden

12、sgloballyWATER用 命令采用两种模式可以定义流体的压缩性：（1）定义比奥系数，比奥模 量；（2）定义流体模量和孔隙率。第一种情况考虑到土颗粒的压缩性，当比奥系数 等丁 1时，颗粒不具有压缩性，第二种情况假定土颗粒不可压缩。PRINT zone property 用命令输出实体单元的性质， PRINT gp命令输出节点 的性质。如给定了流体密度以及水面的位置， 可以用PRINT water命令输出，流体 流动的性质可以采用 PLOT bcontour property 命令或PLOT block property 命令 输出云图。对于各项异性流动，流体区域的整体渗透方向量可以使用实体

13、单元的关 键字 kxx, kyy, kzz, kxy,kxz, kzz 输出。节点初始孔隙压力分布，CONFIG fluid和non-CONFIG fluid的做法是相同的。 即：采用INITIAL pp和 WATER table 命令。FIX pp命令和FREE pp 命令可以为 选定的节点施加的接触约束。流源和会水电可以采用 APPL饰令。1.6节会介绍流 体边界和初始条件。流体计算采用的是SET fluid and SOLVE命令。不同的计算过程可以通过选择。 例如：SET fluid on or off命令为采用或取消流体力学的计算。 这些命令和关键字的 选择取决于流动分析的耦合程度。1.7节将介绍不同耦合程度应该采用的计算过程。 该部分还给出了几个例子，包括：仅作流动的分析和流固耦合的分析等。流体力学的计算可以采用好几种性是查看其结果。PRINT gp pp and PRINTzone pp命令可以输出节点和实体单元的孔隙水压力。HISTORY gp pp and HISTORY zone pp命令可以查看孔隙水压的变化。瞬态计算中，可以用HISTORY fltime命令记录流动时间的变化，从而可