《FLAC3D理论基础屏幕版》由会员分享，可在线阅读，更多相关《FLAC3D理论基础屏幕版（49页珍藏版）》请在金锄头文库上搜索。
1、 .wd.FLAC3D理论根基这局部阐述的是FLAC3D的有关理论。FLAC3D很大一局部是二维FLAC的扩展，而显式有限差分法是FLAC和FLAC3D的共同的理论根基，有关这一局部，可参考FLAC用户手册。尽管如此，二维和三维的方程还是有一些明显的不同，特别是在数学模型的扩展上。这里主要讨论三维模型在FLAC3D中的实现方法。1.三维显示差分模型模型的构成FLAC3D是显式有限差分程序，可以模拟连续三维介质到达平衡状态或稳定塑性流动时的力学行为。这种力学行为，可以通过建设特定的数学模型和特定的数字模拟方法来实现。下面就来阐述这两方面的有关内容。1.1数学模型介质的力学特征可通过一般的力学关系
2、如应变的定义、运动方程等和理想介质的本构方程进展推导。所得到的数学表达式是一系列的偏微分方程及相关变量如：静力学中应力和动力学中的应变速率、速度等。对于特定的具有几何特征和特殊性质的介质，这些方程和变量在给定的边界条件和初始条件下，可以求解。尽管FLAC3D主要是研究处于极限平衡状态下的介质变形及应力状态，但它的模型里可以包含有运动方程是它的一大特色。在进展数字模拟过程中，由于惯性物体将到达稳定状态或平衡状态。1.1. 1符号约定在FLAC3D的拉格朗日公式中，用矢量，其中i=1，3来分别表示介质中点的空间位置、位移、速度和加速度。作为一种符号约定，据上下文的不同，斜体字可以矢量和张量。如：符
3、号表示笛卡儿坐标系下矢量的i分量；Aij表示张量A的i,j分量。还有，表示对xi的偏导数其中可以是标量，也可以矢量或张量的分量。规定：拉力和张力为正。爱因斯坦的求和约定只适用于i,j,ki,j,k=1，2，31.1.2应力给定点的应力状态可用一个对称的应力张量来表示。由柯西定理，假设一个面的单位法矢量为n，则它的拖曳矢量t： 11.1.3 应变速率与转动速率假定介质颗粒以速度v运动，则在无穷小的时间内，发生无穷小应变，相应的应变张量可写为： 2式中是对空间位置矢量的偏导数。在下面的论述中，第一应变速率变量是表征单元体积的膨胀率的量。介质的形变率除了张量外，还有刚体位移v以及转动速率： 3式中：
4、eijk为符号函数，w为转动速率张量，其分量定义为： 41.1.4 运动方程与平衡方程由运动定律的连续形式的柯西运动方程： 5式中：介质单位体积的质量，b作用于单位质量上的体力，dv/dt是速度对时间的偏导数。这些定理控制了介质单元体在力的作用下运动状态。当介质处于静力平衡时，加速度dv/dt为0，由上式得平衡状态下的偏微分方程： 61.1.5 边界条件及初始条件边界条件包括施加于边界上的牵引力和或速度由位移而产生的。另外，还应考虑体力和初始应力状态。1.1.5 本构方程运动方程式5和应变速率的定义式2共有9个方程和15个未知量，这15未知量包括6个应力速率分量、6个应变速率分量及3个速度矢量
5、分量。考虑特定的材料性质，有另外6个相联系的本构方程，通常用下式表示： 7式中：为共转co-rotational应力速率张量，H为一给定的函数，k是与加载路径有关的参数。共转co-rotational应力速率张量，等于给定参考系的介质内一点的应力的偏导数和以瞬时角速度的转动，数学表达式如下： 8式中，是应力对时间的偏导数，w为转动速率张量。1.2 数学模拟FLAC3D有以下三种逼近方法：（1） 有限差分法。假定在变量在空间和时间内线性变化，用变量对空间和时间的一阶导数来近似等于它的有限差分值。（2） 离散单元法。将连续介质离散为等效块体集合体，所有的力施加的作用力与相互作用力作用在三维网格的节
6、点上。（3） 动力学解法。运用运动方程求解所研究系统到达平衡状态时的参量。利用以上逼近方法，连续介质的运动定律可变为节点上的牛顿定律的离散形式。从而可通过显式有限差分法来求解一般的差分方程。等价介质空间偏导数在由速度定义的应变速率用到。因此，为了定义速度变量和相应的空间间距，介质需离散为常应变速率四面体单元，它的顶点即为网格的节点图1。图1 四面体单元的面和节点1.2.1 空间微分与有限差分的近似作为根基，下面由节点运动方程推导四面体的应变速率张量各分量的有限差分公式。四面体的节点1到4，面n指与节点n相对的面见图1。由四面体高斯离散理论，可知： 9上式，左边和右边分别为对四面体体积和面积积分
7、，n为面的单位外法矢量。 对常应变速率四面体，速度场是线性变化的，而每一个面的外法矢量n也为一常值。因此，式9积分后有： 10式中，上标f指面f的相关变量值，指i速度分量的均值。假设速度呈线性变化，则： 11上标l指节点l的值。 将上式代入式10，有： 12 在式9中，假设vi=1，则由离散理论可得： 13所以，式12两边同除以V，则有： 14而应变速率张量则可由下式表示： 151.2.2 运动定律的节点方程 对静力学问题，一般用虚功原理来推导运动定律的节点方程。节点习惯术语也与解相应的平衡方程的所用的相类似式6。一定的时间t内，我们认为静力平衡问题一般可通过以的平衡方程进展求解： 16式中，
8、体力的定义见式5，其中： 17 在这里所用的有限差分逼近格式中，假设介质为受体力B的常应变四面体的连续集合体。根据虚功原理，作用于单个四面体上的节点力fnn=1,4与所受到的四面体应力和等效体力到达一种静态平衡。假设节点的虚速度它在四面体中产生线性速度场和常应变速率，则由此产生的节点力和体力等外力功功率等于内部应力产生的内力功功率。根据1.1.1节的的符号约定上标表示变量的在该节点上的值和爱因斯坦的求和约定，外力功功率可表示为： 18而内力功功率： 19由式15，对常应变速率的四面体有： 20应力张量是对称张量，定义矢量Tl： 21则： 22式17代入式18，有： 23Eb和EI分别体力和惯性
9、力所作的外力功功率。假设四面体内体力为常数，则有： 24EI为： 25由前述有限差分逼近假定，四面体单元内的速度场符合线性变化。因此，我们引进一个参考坐标系它的坐标原点则四面体的中心上，则有： 26式中Nnn=1，4为一线性函数： 27其中，n=1，4为下述方程的解： 28式中，是克罗内克尔增量Kronecker delta。根据重心的定义，所有形如的积分为0，由式27，再将式26代入式23得： 29克莱姆Cramer解式28，再由重心的性质，有： 30将上式代入式29，有： 31同理，将式26代入式25得到： 32将式31和32代入式23： 33在等效问题中，处于静态平衡的四面体，在任意的虚
10、速度下，内力功功率式22等于外力功功率所以，由上述方程式，有： 34假设四面体内加速度在均值附近只有微小的空间变化，则有： 35四面体内，为一常数，由重心的性质式27和30，则上式可写为： 36考虑平衡方程数学上的对称性，用虚拟节点质量mn代替上式中的质量： 37则，式34可写为： 38由等效系统的平衡条件，各个节点上，等效静力-f与有何载荷集中力引起的作用在四面体内和节点上的外力P的和应等于0。为了说明这种情况，我们引入一些符号，上标表示在所有节点中标号为l的节点上的变量，.表示作用在所有与包含节点l的四面体上的变量的和。根据上述约定，可将节点牛顿定律表达为： 39式中，nn介质中的所有的节
11、点总数，节点质量定义为： 40不平衡力F定义为： 41当介质到达平衡时，不平衡力等于0。1.2.3 显式有限差分逼近的时间导数考虑本构方程式7和变形速率与节点速率之间的关系式15，式35可表示为一般的差分方程： 42式中，是指节点l速度的一个子集式39。在FLAC3D中，是利用显式差分方程进展求解的。在这种逼近方法中，假定节点速度随时间间隔线性变化，式42左边的速度对时间的偏导数用中心差分格式，因此，随与位移和力有关的半时间步长进展迭代。节点的速度利用下式进展计算： 43类似地，节点的位置也用中心差分进展迭代： 44由上述可以看出，中心差分可防止一阶误差，具有二阶准确性。因此，节点位移也有如下关系： 45其中，。1.2.4 本构方程的增量形式FLAC3D中，时间内为常数。本构方程式7增量表示为： 46在内，假设位移和位移
