Flac3D中文流体计算.docx

Flac3D中文流体计算.docx

《Flac3D中文流体计算.docx》由会员分享，可在线阅读，更多相关《Flac3D中文流体计算.docx（7页珍藏版）》请在冰豆网上搜索。

Flac3D中文流体计算

Flac3D中文手册

FLAC3D的计算模式中是否需要做孔压分析取决于是否采用configfluid命令。

1无渗流模式〔不使用configfluid〕

即使不使用命令configfluid，仍然可以在节点上施加孔压。

这种模式下，孔压将保持为常量。

如果采用塑性本构模型的话，材料的破坏将由有效应力状态来控制。

节点上的孔压分布可由initialpp命令或watertable命令来设定。

如果采用watertable命令，由程序自动计算水位线以下的静水孔压分布。

此时，必须施加流体密度〔waterdensity〕和重力〔setgravity〕。

流体密度值和水位位置可以用命令printwater显示。

如果水位线是由face关键字来定义的，那么可用命令plotwater命令显示水位。

这两种情况，单元的孔压都由节点孔压值平均求出，并在本构模型计算中用作有效应力。

这种计算模式下，体积力中不反映流体的出现：

用户必须根据水位线以上或以下相应地指定干密度和湿密度。

使用命令printgppp和priintzonepp可分别得到节点或单元孔压。

plotcontourpp命令可绘出节点孔压云图。

2渗流模式〔使用configfluid〕

如果使用命令configfluid，那么可进展瞬时渗流分析，孔压改变和潜水面的改变都可能出现。

在configfluid模式下，有效应力计算〔静态孔压分布〕和非排水计算均被执行。

除此之外，还可进展全耦合分析，这种情况下，孔压改变将使固体产生变形，同时体积应变反过来影响孔压的变化。

如果采用渗流模式，单元孔压仍由节点孔压平均求出。

但这种模式，用户只能指定干密度〔不管是水位以上还是以下〕，因为FLAC3D将流体的影响考虑到了体积力的计算中。

采用渗流模式时，渗流模型必须施加到单元上，使用命令modelfl_isotropic模拟各向同性渗流，modelfl_anisotropic模拟各向异性渗流，modelfl_null模拟非渗透物质。

注意，力学模型为空的单元并不代表渗流模型为空。

流体性质〔参数〕可施加到单元或节点上。

各向同性渗透率、孔隙率、比奥系数和非排水热系数等单元流体性质由命令property施加。

对于各向同性渗流，渗透率通过perm关键字赋予。

对各向异性渗流，渗透率的3个主值采用关键字k1,k2,k3赋予，主方向由关键字fdip,fdd,frot确定。

渗透率的主方向服从右手系统。

fdip和fdd分别为k1和k2确定的平面的倾向和倾角。

frot为k1轴和倾角矢量的旋转角。

如果不特别指定，比奥系数默认为1，孔隙率默认为0.5。

节点的渗流性质由命令initial指定。

这些性质包括流体重度、流体体积模量、比奥模量、流体抗拉强度和饱和度。

每种性质在空间上都可以变化。

流体重度也可以用water命令给出。

在渗流模式里，有必要知道可压缩性被定义在以下两种参数中：

〔1〕比奥系数和比奥模量；〔2〕流体体积模量和孔隙率。

第一种参数表征的是固体颗粒的可压缩性〔对不可压缩颗粒，比奥系数设为1〕。

对第二种参数，固体颗粒被认为是不可压缩的。

单元属性可由命令printzoneproperty显示，节点属性由printgp命令显示。

流体重度，如果随着水位位置被确定，那么可由printwater命令显示。

渗流性质可由命令plotbcontourproperty显示。

对于各向异性渗流，渗透率的各球形分量可通过使用单元的属性关键字kxx,kyy,kzz,kxy,kxz,kyz来显示〔注意，这些球形分量不可被直接初始化〕。

初始节点孔压分布的施加对于渗流模式和非渗流模式都是一样的〔如，要么用initialpp命令或用watertable命令〕。

在指定节点可用命令fixpp或freepp对孔压固定或释放。

流体涌入或渗漏或可由命令apply施加。

渗流计算由命令setfluid和solve控制。

如，setfluidon或off命令开启或关闭渗流计算模式。

具体使用开启或关闭模式取决于渗流分析的耦合程度。

渗流分析结果以下面这些命令给出。

命令printgppp和printzonepp分别给出节点和单元孔压。

节点和单元孔压历史可由命令histroygppp和命令historyzonepp进展监测。

对于瞬时计算，孔压与时间的关系可由命令historyfltime监测。

命令plotcontourpp绘出节点孔压云图。

命令plotcontoursaturation绘出饱和度云图。

命令plotfluid绘出流量矢量图。

渗流模式的所有信息由命令printfluid命令给出。

FISH还提供了一些渗流变量。

其中一个与节点有关的变量gp_flow，只能通过FISH函数使用。

该变量描述了通过节点的净流入或流出量。

因为可以提供一个系统总的流入或流出量，这些流量的统计在孔压固定的边界是很有用的。

渗流边界条件，初始条件

FLAC3D默认为不透水边界，即认为所有节点上的孔压随着从邻近单元流入或流出的量发生自由变化。

可以使用命令fixpp将节点上的孔压设为“自由〞，也可使用freepp使节点上的孔压“固定〞。

如果孔压固定，流体可以在外边界上流入或流出节点。

下面总结这两种边界条件

的影响：

1，孔压自由

这是默认的不透水边界条件。

节点与外界之间不发生流量交换。

系统根据当前饱和度值和流体是否形成涡凹现象来计算压力和饱和度变

化。

2，孔压固定

这是一种流体通过外界流入或流出的边界条件。

如果设定孔压为0，饱和度才可能变化。

否那么，饱和度被设为1〔FLAC3D假设孔压只在完全饱和材料中存在〕。

孔压不能被固定在低于拉力极限的值，如果

出现这种情况，FLAC3D会将其设定到拉力极限值。

如前所述，边界条件不是任意的。

FLAC3D在进展计算前会“检查〞并“修正〞这些条件。

可使用fixpp命令将孔压固定在某个值，也可在外边界或边界上使用命令applypp。

如果边界条件被用于一个非外表节点，那么必须加关键字interior。

apply命令具有可以用“历史〞命令进

行监测的优点。

渗流边界条件可以通过apply命令用在单个或局部节点、单元面或单元上。

命令applypwell为边界节点指定了一个流入或流出井。

如果加上interior关键字，那么该条件用于部节点。

命令applydischarge和applyleakage为边界单元的外表分别指定了涌出和渗漏边界条件。

命令applyvwell为指定区域的单元提供一个流速。

这些边界

条件除了applyleakage外，均可使用history监测命令。

具有固定孔压节点就好似是流入源或流出源。

没有直接的命令显示这些节点的流入或流出量。

但可通过FISH变量gp_flow来记录。

孔压的初始分布，孔隙率，饱和度和流体属性可通过命令initial或property施加。

如果还加了重力，那么孔压初始分布应与重力梯度，水的重度和节点饱和度和孔隙率相容。

如果这些初始分布不相容，那么计算开场时所有单元中将出现流体流动。

因此，应在模拟开场时设一定

的计算步来检验初始条件是否相容。

如果模型中含有接触面，有效应力将沿着这些接触面进展初始化〔即：

在节点应力初始化时，认为接触面应力包含孔压〕。

waterlable命令将包含沿着接触面的孔压，这是因为定义在单元节点上的孔压也在接触面节点上。

如果接触面的上下两面连在一起，在没有阻力时，将发生穿越接触面的流体流动。

但程序不对沿着接触面的流体流动

〔裂隙流〕进展计算。

单渗流与渗流耦合问题

FLAC3D既能进展单渗流分析，也能进展固流耦合分析。

耦合分析可由FLAC3D置力学模型完成。

但要注意，渗流模型中的空单元并不是力学空单元。

必须用命令modelfl_null给单元赋予流体空属性。

对于耦合过程，FLAC3D提供了几种计算模式。

其中之一是假设孔压一旦被赋予便不再改变。

该方法并不要求任何额外空间存储计算过程。

除此之外涉及到渗流的计算模式都要求使用命令configfluid。

命

令modelfl_iso使所有单元中都能发生渗流。

不同的耦合计算模式在下面讨论。

一般情况下，在能跟所模拟问题的物理过程相似的情况下，应使用尽可能简单的模式。

计算模式的选择

根据以下几个方面确定。

时间比例

对所需模拟的渗流或耦合问题用FLAC3D估计与涉及的不同进程相关的时间比例是非常有用的。

对有关研究问题的时间度量和扩散性的认识有助于估计最大网格宽度、最小区域尺寸、时步大小和计算可行性。

如果不同进程的时间比例相差太大，那么很可能采用一种简单的〔非

耦合〕方法。

时间比例可用特征时间给出。

以下这些由量纲分析得出的定义，都是基于解析的连续源理论表达式。

它们可用于得出FLAC3D分析的大

致时间比例。

力学过程特征时间、流体扩散过程特征时间

流体扩散率

FLAC3D中使用了取决于控制过程的储水系数的几种形式：

流体存储系数、地下潜水相存储系数、弹性存储系数

以上定义，有几点特性值得注意：

〔1〕因为FLAC3D中显式的时步对应于最小区域息从一个节点传到下一节点所需要的时间，时步的大小可用计算特征时间公式中特征长度的最小区域来估计。

重要的是注意FLAC3D中在用流体扩散率〔即使是在耦合模拟中〕计算显式流体时步。

因此，时步的大小可

用特征长度的最小区域尺寸来估计。

〔2〕在饱和流体问题中，简化的体积模量不但导致时步的增加，同样导致到达稳定状态时间的增加，所以总步数增加，该总步数可用模

型和最小区域的特征长度来估计。

〔3〕在局部饱和流体流动问题中，可通过调整流体体积模量加速收敛以趋于稳定状态，但要注意不可将体积模量减小太多以至产生数值不稳定。

数值稳定条件能由流体储量在一个特征长度区域的高度上必

须保持低于地下潜水储量的要求推导而出。

〔4〕为防止扩散问题中的边界效应，模型的特征长度必须大于某个尺度。

同样，最小模拟时间由某个关系式控制。

〔5〕在耦合流体问题中，实际扩散率由流体刚度与岩土介质的刚度

比来控制。

完全耦合模拟方法的选择

用FLAC3D进展完全耦合的准静态固流耦合分析通常要消耗大量时间，且有时候并不必要。

很多情况下，可使用不同程度的非耦合方法简化分析并加快计算速度。

下面的例子给出了对应于流固耦合的不同水平的FLAC3D模拟方法。

选择计算方法时有3个主要的因素需要

考虑：

〔1〕模拟时间比例和扩散过程的特征时间；

〔2〕耦合过程中强制扰动特性；

〔3〕流固刚度比。

时间比例

首先通过从扰动的开场阶段计算时间来考虑时间比例因素。

定义分析所需要的时间〔模拟时间〕，对应于耦合扩散过程的特征时间。

短期行为〔不排水〕

如果对应于耦合扩散特征时间，分析所需时间非常短，在模拟结果中流体流动的影响几乎可以忽略不计，那么可采用不排水模拟〔configfluied,setfluidoff〕。

数值模拟中不涉及真实的时间，但如果给流体体积模量一个实际值，那么体积应变将导致孔压的变化。

长期行为〔排水〕

如果分析所需时间大于耦合扩散特征时间并在模拟时间到达时排水，那么孔压场可不耦合到立场中。

稳定状态的孔压场可用单纯流动模拟确定〔setfluidon,setmechoff〕〔不代表扩散率〕，然后力学场可通过在设置流体模量为0的力学模式中〔setmechon,setfluidoff〕计算到平衡状态获得。

〔严格说，这种方法仅对弹性材料有效，因为塑性材料力学行为是与路径有关的〕。

另外一种描述时间比例的方法是不排水行为和排水行为。

严格地说，不排水表示与模型外界无流体交换。

排水那么是与模型外界有完全的流体交换，这就意味着流体压力能在各处到达平衡。

由于不排水试验通常所需时间很短，而排水试验那么需要很长的时间以使多余的流体压力消散，因此，“排水〞和“不