流体力学(FLAC3D).docVIP

流体力学 — 单向流动 简介 FLAC3D通过具有渗流性的实体比如土来模拟流体的流动。流动模型的建立可以独立于力学计算而自动完成，或者说可以与力学模型同时建立，这样就可以考虑流体与土体之间的相互作用。流固耦合的一种类型是“固结”，即:空隙水压力逐渐消散而导致土体的沉降。这个过程包括两种力学反映：一，空隙水压的改变导致有效应力的变化，这将影响到土体的力学反映（如：有效应力的减小可能导致塑性区的产生）；二，力学实体中某一区域的流动会随着空隙水压的改变而改变。 该程序可以计算完全饱和情况下的流动，也可以模拟具有自由水面的流动。模拟具有自由水面的流动时，自由水面以上的部分空隙水压等于0，气相将不参与计算。对于不考虑毛细水压力颗粒较粗的材料可以采用这种模拟方法。 流体计算就有以下特点： 1 根据各项同姓和各项异性的渗流计算，相应采用两种流体运动定律。流动中的null材料用来模拟流动范围内的非渗流材料。 2 不同区域可以拥有不同的流动模型（isotropic, anisotropic or null）和模型参数。 3 可以事先指定流体的压力、流量、非渗流区边界条件。 4 流体源可以以电源，也可以以体源的形式插入到材料中，这些源对应于流体的流入或流出，可以随着时间而变化。 5 对于完全饱和流动，可以采用显式和隐式两种算法，但对于非饱和流动则只能采用显示计算。 6 任何力学和温度计算模型都可以与流体模型一起使用，在耦合计算中，可以考虑饱和体的压缩性和热膨胀性。 7.流体与力学计算的耦合通过提供比奥系数来实现。 8．与温度的耦合计算可以通过提供线性热膨胀系数αt和不排水温度系数β（undrained thermal coefficient,可能翻译的不对）来实现。 9.热-流体计算以线性理论为基础，假定材料参数为常数，不考虑对流。流体与实体的温度保持局部平衡。非线性行为可以采用FISH语言改变孔隙压力、材料特性来实现。。 动力水压力的产生以及循环荷载引起的液化也可以用FLAC3D模拟（3.4.4节）。FLAC3D不考虑部分饱和体中的毛细水压、粒子之间的电、化学力。但是可以根据材料的局部饱和情况，孔隙率或其他相关变量利用FISH语言施加内部应力来实现上述几种作用力。与此类似，由于气体溶解引起的流体刚度变化的影响也是不能直接模拟的，但是FISH编写一个压力，时间的函数就可以用来改变局部流体的模量 本章主要包括7个方面的内容： 1.2节和1.3节中分别介绍流体力学，流-力耦合以及流-热耦合的数学理论模型和数学表达式。 1.4中介绍流体力学模拟采用的计算模型和相关的命令。 1.5节介绍流体力学模拟采用的材料参数，各参数的单位。 1.6介绍流体力学模型采用的不同的边界条件，初始条件，流源和汇集。 1.7介绍流体分析的一般过程，该部分还包括几个例子，建议读者在做自己的流体分析以前首先研究清楚我们提供的例子。 1.8是几个考核题，用来验证flac3d流体分析的准确性。 1.9总结了所有与流体计算相关的命令，以供参考。 读者应该首先熟悉力学计算，然后才能将流体，力一起来考虑。流-力耦合是一个复杂的问题，需要读者有很好的洞察力来判断结果的正确性。 （1.2 1.3 我看就没有必要翻译了，有兴趣看看朱伯芳院士的书就可以了，理论是一样的） 1.4 流固耦合的计算模式 孔隙水压的计算模型和命令取决于节点是否已经设置了流动（也就是说，CONFIG fluid命令是否已经采用）1.4.1和1.4.2将分别介绍这两种情况。方便起见：下面所有的命令在文末1.9都作了汇总。 1.4.1 Grid Not Configured for Fluid Flow 如果不给出CONFIG fluid命令，仍旧可以在节点上设置孔隙水压力。这种计算模式下孔隙水压力不变，但是当材料为弹塑性体时有可能不收敛。 使用initial pp 命令加a gradient, 或用 WATER table 命令可以指定节点上孔隙水压力的分布。 如果使用WATER table 命令，水面一下部分的孔隙静水压力的分布程序将自动进行计算。这种计算模式，必须用WATER density设置流体密度，SET gravity 设置重力加速度。当用水面定义好以后，采用PRINT water命令可以输出流体密度，水面的位置，PLOT water命令可以显示出水面来 两种情况下微元体的孔隙压力计算是将个节点的值平均，从而计算有效应力。这种计算模式，体力的计算不考虑流体压力：水面以上和以下材料的天然重度和饱和重度要有用户自己指定。PRINT gp pp 和PRINT zone pp可以输出节点和体的孔隙水压力。PLOT contour pp 命令用来绘制节点孔隙水压的云图。 在《用户手册》中的3.6.3有这种计算模式的一个简单