FLAC3D原理Word文档下载推荐.docx

FLAC3D原理Word文档下载推荐.docx

《FLAC3D原理Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《FLAC3D原理Word文档下载推荐.docx（18页珍藏版）》请在冰豆网上搜索。

〔9〕修正的Cam粘土材料模型〔ModifiedCam-clayModel〕

除上述本构模型之外，FLAC3D还可进展动力学问题、水力学问题、热力学问题等的数值模拟。

在边界条件及初始条件的考虑上，FLAC3D软件十分灵活方便，可在数值计算过程中随时调整边界条件和初始条件。

FLAC3D具有强大的后处理功能，用户可以直接在屏幕上绘制或以文件形式创立或输出打印多种形式的图形、文字，用户还可根据各自的需要，将假设干个变量合并在同一幅图形中进展研究分析。

FLAC3D软件还可对各种开挖工程或施加支护工程等进展数值仿真模拟，软件自身设计有锚杆、锚索、衬砌、支架等构造元素，可以直接模拟这些支护于围岩〔土〕体的相互作用。

FLAC3D拥有可以自行设计的FISH语言，用户可根据自身需求，自己设计材料的本构模型、屈服准那么、支护方案、复杂形状的开挖方式等工作。

特别注意的是，岩石是一种脆性材料，当外荷载到达岩石强度后，材料发生断裂破坏，产生弱化现象，应属于弹塑性体。

在FLAC3D中，一般对于弹塑性材料，判断其破坏与否的根本准那么有两个，即Drucker-Prager准那么和Mohr-Coulomb准那么。

根据室岩石力学性质试验结果，其典型应力应变曲线反映出岩体破坏包络线符合莫尔—库伦屈服准那么，故本次建立的本构力学模型选择莫尔—库伦弹塑性材料模型为宜。

2.2.2FLAC3D分析计算原理

计算所采用的数学模型是根据弹塑性理论的根本原理〔应变定义、运动定律、能量守衡定律、平衡方程及理想材料的连续性方程等〕而建立的。

2.2.2.1根本约定

在数学及数值模型的表达式中，符号有一定的约定含义，一般[A]表示量，Aij表示量[A]的〔i,j〕分量，[a]表示矢量，ai表示矢量[a]的i分量，α,i表示α对xi的偏导数。

xi，ui，vi和dvi/dt，〔i=1,3〕分别表示一点的位置矢量分量、位移矢量分量、速度矢量分量和加速度矢量分量。

2.2.2.2数学模型

〔一〕柯西〔Cauchy〕应力量与柯西公式

对于一个具有体积V的封闭曲面s的物体，在其上取一外表元素∆s，这个外表元素的单位外法向矢量为n，在某一时刻t，在外表元素

对于连续介质中一点，作用着对称的应力量σij，根据∆s上作用有力∆P，那么极限

称为外表力。

假设用ti表示T的分量，那么在三维直角坐标系中可有关系式

〔1〕

这个关系式称为柯西公式，其中，σij称为柯西应力量。

〔二〕应变速率和旋转速率

如果介质质点具有运动速度矢量[v]，那么在一个无限小的时间dt，介质会产生一个由vidt决定的无限小应变，对应的应变速率分量ξij为

〔2〕

而其旋转速率分量ωij为

〔3〕

〔三〕运动及平衡方程

根据牛顿运动定律与柯西应力原理，如果质点作用着应力σij与体力bi，且具有速度vi，那么在无限小时间段dt，它们之间的关系为

〔4〕

式中，ρ为质点密度。

〔4〕式称为柯西运动方程。

当质点的加速度为零时，上式变为静力平衡方程

〔5〕

〔四〕本构方程

上述〔4〕式与〔5〕式组成的方程组中含有9个方程，15个未知量，其中12个是应力与应变速率分量，3个是速度分量。

其余6个关系式那么由本构方程提供，本构方程一般具有如下形式

〔6〕

式中，

为应力变化速率，H表示一个特定的函数关系，κ为与荷载历史有关的参数。

2.2.2.3数值模型

FLAC3D的数值剖分网格在计算中是按照四面体进展的，四面体的节点也既是网格剖分的节点，因此，每个计算单元有4个面和4个节点〔见图2-1〕。

〔一〕空间微分的有限差分逼近

对于一个计算单元，假设部各质点速度为一连续的矢量场[v]，那么根据高斯〔Gauss〕积分原理有

〔7〕

式中，[n]为外法向单位矢量场。

由于单元体的应变速率是连续的，因此可以近似认为速度是线性变化的，那么〔7〕式可用下面的求和公式近似逼近

〔8〕

式中，V为单元体体积，

为单元体某一面的面积，f为单元体面数，f=1,4，

为面平均速度的i分量。

由于速度场是线性的，那么有

〔9〕

式中，l为单元体节点数，l=1,4；

vil为l节点的i速度分量。

将〔9〕式代入〔8〕式可得

〔10〕

根据正交原理有

〔11〕

那么将〔10〕式两边除以V，并将〔11〕式代入可得

〔12〕

因此有

〔13〕

〔二〕运动方程的节点公式

根据前面对质点运动方程的讨论，对于连续介质，当处于平衡状态时，其平衡方程为

〔14〕

〔15〕

介质可由假设干个作用着体力[B]，各自产生一定变形的四面体组成，设节点力为[f]n，n=1,4，利用虚功原理，假定单元体节点具有速度δ[v]n，部具有变形速率δ[ξ]，那么由节点力[f]n和体力[B]所做的外功率与σij所做的功率应当相等。

外功率E可用下式表示

〔16〕

而功率I为

〔17〕

由〔13〕式，对于恒定应变速率单元体可有

〔18〕

由于应力量是对称的，定义矢量Tl

〔19〕

可得

〔20〕

将〔15〕式代入〔16〕式得

〔21〕

式中，Eb和EI分别是体力ρbi和惯性力产生的外功率，对于单元体恒定的体力ρbi，Eb可写成

〔22〕

而EI可以写成

〔23〕

如前所述，单元体部速度场以线性变化，为方便描述，选取单元体质心为原点，

、

为坐标轴，可有

〔24〕

式中，Nn（n=1,4）是具有如下形式的线性函数

〔25〕

〔n=1,4〕那么由下面方程决定

〔26〕

δnj为克罗克尔〔Kronecker〕记号，根据质心定义可得

〔27〕

用克莱默〔Cramer’s〕法那么解得

，那么

〔28〕

而

〔29〕

那么

〔30〕

在稳定状态下，功率与外功率必定相等，那么有

〔31〕

上式最后一项，如果单元体ρ恒定，那么根据质心定义有

〔32〕

假设将ρV/4看作是假定的节点质量mn，那么上式变为

〔33〕

故

〔34〕

根据牛顿定律有

〔35〕

式中，Fi<

l>

为节点l所受力的i分量，M<

为节点l的质量，

为节点l在Fi<

作用下产生的加速度，nn为包含在全部连续介质中单元节点总数。

设质点质量

，不平衡力[F]<

为

〔36〕

〔三〕时间微分的显式有限差分

对于上面〔35〕式，可写成

〔37〕

FLAC3D计算中，对于上式的逼近是在一个足够小∆t时间段，采用中心差分格式，即上式变为

〔38〕

同理，节点的坐标差分公式与位移差分公式分别为

〔39〕

〔40〕

〔四〕本构方程的增量形式

FLAC3D中，假定在∆t时间段，速度是不变的，本构方程的增量形式可表示为

〔41〕

记

〔42〕

那么∆εij为∆t时间段应变增量。

应力增量由下式确定

〔43〕

由下式定义

〔44〕

〔45〕

2.2.2.4本构模型

〔一〕弹塑性理论的增量关系

一般来讲，弹塑性体破坏的判断准那么具有如下形式

〔46〕

式中，f为一特定的函数关系，

为n维应力矢量i分量。

弹塑性体的变形

为弹性变形

与塑性变形

的总和，即

〔47〕

弹性应变与应力之间的关系为

〔48〕

其中，Si为线性函数。

塑性变形可用下式表示

〔49〕

式中，λ为一常数，g为

的某一函数关系。

将〔47〕、〔49〕式代入〔48〕式得

〔50〕

对于新应力

，仍有

〔51〕

当

为线性函数时，〔51〕式可写成

〔52〕

其中，

。

由于

，再将〔50〕式代入〔52〕式可得

〔53〕

假设记

那么有

〔54〕

〔55〕

那么新应力为

〔56〕

〔二〕莫尔—库伦〔Mohr-Coulomb〕模型本构关系

本次计算的本构模型选择的是莫尔—库伦弹塑性材料模型，弹塑性体可产生弹性及塑性两局部变形，根据虎克定律〔Hooke’slaw〕，应力应变的关系为

〔57〕

其中，[E]为刚度矩阵。

对于弹性变形，应力增量可由下式确定

〔58〕

式中，α1、α2可由剪切模量与体积模量得出

〔58〕式还可写成

〔59〕

〔三〕破坏准那么与流动法那么

假设

莫尔—库伦破坏准那么具有形式fs=0，其中

〔60〕

对于拉破坏，有ft=0的形式，其中

〔61〕

上两式中，φ为摩擦角，c为聚力，σt为抗拉强度，

〔62〕

对于材料的抗拉强度σt，其值不可能超过σ3〔见图2-2〕，σt的最大值为

〔63〕

对于塑性流动状态，设gs与gt分别为剪切破坏和拉破坏所对应的函数关系，有

〔64〕

〔65〕

式中，ψ为膨胀角，

〔66〕

流动准那么可定义为函数关系