流固耦合手册现场翻译.docVIP

1 流体－固体耦合与单相渗流 1.1介绍 FLAC3D模拟了流体流过可渗透的介质，例如土体。渗流模型可以独立于通常FLAC3D的固体力学计算，而只考虑渗透；或者为了描述流体和固体的耦合特性，与固体模型并行计算。例如，固结是一类流固耦合的现象，在固结过程中孔隙压力逐渐消散，从而导致了土体的位移。这种行为包含了两种力学效应。其一，孔隙水压力的改变导致了有效应力的改变，有效应力的改变影响了固体的力学性能，例如有效应力的降低可能引发塑性屈服；其二，土体中的流体对孔隙体积的变化产生反作用，表现为孔隙水压力的变化。 本程序可以不仅可以解决完全饱和土体中的渗流，也可以分析有浸润线定义的饱和与非饱和区的渗流计算。该条件下，浸软面以上的土体的孔隙水压力为零，气体是被动的（气体的压力考虑成负的）。这种方法用于颗粒比较粗的毛细现象可以忽略土体。 渗流分析中有如下的特征： 对应于渗流各向同性和各向异性渗透性采用两种不同的流体传输定律。渗流区域中的不可渗透的区域用流体的null材料定义。 不同的zone可以赋予不同的渗流模型(isotropic, anisotropic or null)和属性。 流体压力，涌入量，渗漏量和不可渗透边界都可以定义。 土体中可以加入抽水井，考虑成点源或者体积源。 计算完全饱和土体中的渗流问题，可以采用显式差分法或者隐式差分法；而非饱和渗流问题只能采用显式差分法。 渗流模型可以和固体力学模型和传热模型一起使用。在耦合问题中，可以考虑饱和材料的压缩和热膨胀。 流体和固体的耦合程度依赖于土体颗粒（骨架）的压缩程度，用Biot系数表示颗粒的可压缩程度。（即用biot系数确定颗粒变形，来模拟力学与流体的耦合） 用线性热膨胀系数和非排干导热系数来考虑热耦合计算。 对于热-流体-流动计算是基于线性理论，假设材料属性为常数，不考虑对流。流体和固体温度局部平衡。如果要模拟非线性行为，需要制定孔隙压力和用fish函数确定材料属性。 由于循环荷载引起的动水压力和液化问题也可以用FLAC3D模拟。FLAC3D不考虑毛细现象，土体颗粒间的电、化学作用力。然而，可以根据土体的局部饱和度，孔隙率，或者其他的变量，通过编写一段FISH语言来考虑这种力。类似的，由于液体中溶解了空气而引起的液体刚度变化，也不能显式的模拟，而通过FISH将液体刚度表示为压力，时间和其他变量的函数。 这以章节可以分为七个主要部分： 数学模型描述（章节1.2）和相应的数值方法（章节1.3）（单相渗流和流固耦合计算）。 计算模式和渗流分析相关的命令。（章节1.4） 渗流分析所需要的材料属性的单位。（章节1.5） 不同边界条件，初始条件，和流体的源和汇的的描述。（章节1.6） 求解单相渗流问题和流固耦合问题的推荐方法。（章节1.7）建议先练习这里的例子再做自己的流体分析。 （章节1.8）提供了几个验证算例，演示了用falc3d的流体计算方法的准确性。 （章节1.9）总结了与流体计算的相关命令。 用户在尝试解决流固耦合问题之间，强烈建议先熟悉一下FLAC3D求解固体力学问题的步骤。流固耦合的力学行为通常非常复杂，需要用户对计算结果是否真确需要远见和判断。在开始计算一个大的项目支前，有必要在一个网格比较少的例子上作试验，尝试不同的边界条件和建模策略。“浪费”在试验上的时间，必定可以通过计算时间来弥补回来。 1.2 FLAC3D流体－热-力学方程描述 大多数工程分析（包括变形和扩散同时发生机理），通过非耦合的技术来模拟（章节1.7）。FLAC3D除了提供这些非耦合的计算模式，也提供流体-热-力学耦合计算的选项。在该方法中，多孔介质的力学响应通过瞬态流和热的条件来研究。虽然本章节主要考虑变形－流体扩散的问题，流体－热－固体力学耦合作用的方程完整的列出。 FLAC3D中力学变形－流体消散的描述在准静态Biot理论的框架下完成，而且可以应用到多孔介质中遵循Darcy定律的单相渗流的问题。不同类型的流体，包括气体和水，可以用这个模型表述。 描述多孔介质中流体渗流的变量是孔隙水压力，饱和度和特定排水向量的三个分量。相关的变量遵循流体的质量守恒定律，达西定律，和流体响应孔隙水压力改变，饱和度改变，体积应变改变和温度改变的本构方程，和与在非饱和区域与孔隙压力和饱和度相关的状态方程。力学本构定律中考虑了孔隙压力和温度的影响 ，以实现热-流体流动-力学耦合。 假设体积应变和温度应变率，在流体本构关系中使用达西定律，根据孔隙压力和饱和度服从不同方程以解决特定的几何、属性、边界和初始条件的问题，取代了质量平衡方程。  1.3数值方程 流体连续性方程： 在FLAC3D的数值方法中，渗流区域离散为由8节点定义的砖块状的zone。孔隙压力和饱和度设定为节点变量（如果有力学和热的计算，也采用相同的离散方法）。从内