有限差分法FLAC2Dppt课件

1、有限差分法（1）FLAC2D,1. FLAC简介,1.1 FLAC特点,FLACFast Lagrangian Analysis of Continua FLAC建立在拉格朗日算法基础上，采用有限差分显式算法来获得模型全部运动方程(包括内变量)的时间步长解，从而可以追踪材料的渐进破坏和垮落，这对研究采矿工程设计是非常重要的。 FLAC适用模拟计算岩土材料力学行为，特别适合模拟大变形和扭曲，包括材料的高度非线性(应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的应力渗流耦合、热力耦合以及动力学问题等,1.2 岩体力学本构模型,1) 各向同性弹性材料模型； (2) 横观各向同性弹性材

2、料模型； (3) 莫尔-库仑弹塑材料模型； (4) 应变软化/硬化塑性材料模型； (5) 双屈服塑性材料模型； (6) 遍布节理材料模型； (7) 空单元模型，可用来模拟岩体开挖和开采,岩石,充填体,岩石应变软化/硬化模型,断层,断层模型,1.3 支护构件与本构模型,1) 杆单元模型； (2) 梁模型； (3) 桩模型； (4) 群支架模型,预应力砂浆锚杆的模拟体系,FLAC举例,a)巴西试验 (b) 边坡滑移 （间接拉伸,P,2. FLAC理论基础,2.1 有限差分法,2.2 简单力学问题分析,m,F(t,牛顿运动定律,对于一个连续体，可写成广义形式,式中, 体密度, xj 坐标矢量 (x,

3、y) ij 应力张量分量 gi 重力加速度,2.3 应力应变关系,除了运动定律，连续介质必须服从本构关系应力与应变之间的关系 对于弹性材料,一般地，本构关系可表述为,式中,式中，ijKronecker记号； t时间步； G，K分别是剪切模量和体积模量,2.4 有限差分基本格式,在有限差分法中, 每个在前一个公式中导出的量（运动量或应力应变）由所在网格特定位置处的相关变量代数式所代替。 代数式是完全显式的；在代数式右侧的所有量都是已知的。在一个计算时步，FLAC网格每个元素（域或结点）与相邻单元表现出在物理上是隔绝的,选取的时步非常小，乃至在此时步间隔内实际信息不能从一个单元传递到另一个单元（事

4、实上，所有材料都有传播信息的某种最大速度）,计算循环的基本要求,2.5 基本显式计算循环,平衡方程 (运动方程,应力-应变关系 (本构方程,对于所有结点,对于所有单元,新的应力,结点力,高斯定理,应变速率,速度,例如, 弹性,FLAC的网格内部由三角形构成，三角形组合成四边形单元。推导差分方程的方式如下,单元叠加 三角形单元 结点力向量 速度向量,高斯定理,用于推导任意形状单元的有限差分公式,对于多边形，公式为,下式用于计算应变增量, eij,一旦计算出全部应力，可以从作用每个三角形边界上产生的牵引力计算得到结点力。例如,然后，用“传统”的中间差分公式获得新的速度和位移,大变形模式,Fi,ni

5、(1,ni(2,s(1,s(2,在时间域上的解算方法,所有单元,非线性定律,所有结点,重复n 个时步 时步内不作迭代,信息在每个时步内 不在单元之间传播,假设 (u)固定,假设 (F)固定,改正，若,p-波速度,位移 u,力 F,x,F,应力,u,数值结点,显式,隐式,单元,结点,每个时步内解算整套方程 如果存在非线性，在时步 内迭,方法比较,显式, 时间-进程,隐式, 静态,1. 可以遵循非线性定律而无需内部迭代, 因为本构计算中位移被“冻结”。 2. 对同样问题，计算时间增加N3/2 。 3. 物理不稳定不会引起数值不稳定。 4. 因为无需存储矩阵，用少量内存 可以模拟大型问题。 5. 大

6、应变、大位移和转动模拟无需 额外机时,1. 整个迭代过程需要遵循非线性定律。 2. 解算时间增加 N2 甚至 N3。 3. 模拟物理不稳定性困难。 4. 需要大内存，或大容量硬盘存储。 5. 模拟大变形明显需要更多的机时,动态衰减,在动态衰减过程中，结点按照牛顿定律运动。结点的加速度与非平衡力成正比。该算法确定这组位移将导致系统趋于平衡，或表征失稳状态。 关于动态衰减的两个重点： 选择时步 阻尼作用,时步,为满足数值稳定，时步必须满足,式中， Cp 与 1 /mgp成正比。对于静态分析，结点质量按比例度量 ，致使局部临界时步等于（ t=1），这样提供了收敛的最优时步。然后，对结点惯性质量进行调

7、整来实现稳定性条件,注意：重力质量不受影响,阻尼,速度-比例阻尼导致影响计算的体积力。 FLAC中采用了局部阻尼结点处的阻尼力与不平衡力大小成正比，正负号的规定确保振动状态被阻滞,局部阻尼,阻尼力 Fd 为,在 FLAC 中，通过监测不平衡力比率（不平衡力 Fi 与施加力 Fm之比）来确定静力状态 当 Fi / Fm 0.001时，则（计算机缺省）认为模型处于平衡状态,式中，Fi是不平衡力,在运动方程中引进阻尼力,如图，一维杆件用数个等尺寸的有限差分网格划分，杆件的密度为，杨氏模量为E,1.6 FLAC计算流程举例,对于固体材料，微分形式的本构方程为,1,运动方程(或平衡方程)为,假设杆件无侧

8、向约束。对于单元i，与(1)式对应的中间有限差分公式为,2,3,式中，(t)表示其变量是在时刻t确定的，上标i表示单元或结点编号,4,对结点i的有限差分形式运动方程为,或写成,5,在显式算法中，所有有限差分方程右端的值均是已知的。因此，必须先用(3)式算出所有单元的应力，然后再由(5)式和(6)式计算所有结点的速度和位移。 在概念上，这个过程等价于对变量值“同时”更新，而不是像其它方法那样，方程右端混存有 “新”、“旧”值，对变量值“依次”更新,积分后得出位移,6,3. FLAC软件应用,1）设计计算模型的尺寸 （2）规划计算网格数目和分布 （3）安排工程对象（开挖、支护等） （4）给出材料的

9、力学参数 （5）确定边界条件 （6）计算模拟,3.1 建模,FLAC 基本称谓,单元编号,结点编号,3.2 网格生成 Grid i,j 例如：grid 30,20,3.3 网格规划：Gen x1,y1 x2,y2 x3,y3 x4,y4 例如：Gen 0,0 0,10 10,20 20,0,3.4 分区规划网格。例如：Gen xI1,yI1 xI2,yI2 xI3,yI3 xI4,yI4 i=1,10 j=1,21 （I区）Gen xII1,yII1 xII2,yII2 xII3,yII3 xII4,yII4 i=10,20 j=1,21 (II区,3.5 特殊形状的网格（1）圆形 gen c

10、ircle xc,yc rad,xc,yc,rad,3.5 特殊形状的网格（2）弧线 gen arc xc,yc xb,yb theta,xc,yc,xb,yb,3.5 特殊形状的网格（3）直线 gen line x1,y1 x2,y2,x1,y1,x2,y2,4）任意形状 tab 1 x1,y1, x2,y2, ,xn,yn, x1,y1 gen tab 1,x1,y1,x2,y2,x3,y3,x4,y4,x5,y5,x6,y6,3.6 赋给单元材料性质,mod e (弹性) prop d 1800e-6 bu 12.5 sh 5.77 i=1,20 j=1,10 prop d 2400e-

11、6 bu 1250 sh 577 i=1,20 j=11,20 mod m （弹塑性Mohr-Coulumb准则） prop d 1800e-6 bu 12.5 sh 5.77 c 0 fri 20 ten 0.015 reg i,j,特别提示,1）量纲统一,2）参数换算,3.7 赋给模型边界条件,1）固定边界 （结点） Fix x i=1, j=1,21 Fix y i=1,21 j=1,Fix x,i=1,j=1,Fix y,3.7 赋给模型边界条件,2）施加边界力 (结点) apply yf -10 i=1,21 j=21 或 apply syy -10 i=1,21 j=21 appl

12、y xf -5 i=21, j=1,21 或 apply sxx -5 i=21, j=1,21,i=21,j=21,3.7 赋给模型边界条件,3）赋单元内应力 （单元） ini sxx -10 i=1,20 j=1,20 ini syy -5 var 0 4 i=1,21 j=1,21,i=21,j=21,3.8 计算,Set grav 9.81 Set large Step 1000 Save test.sav,3.9 结果显示,Plot grid 显示网格 Plot bo 显示边界 Plot plas 显示塑性区 Plot sig1 fi 显示最大主应力1 Plot sig2 fi 显示

13、最小主应力2 Plot sdif fi 显示主应力差(1- 2) Plot str 显示主应力矢量场 Plot xdis fi 显示X方向位移 Plot ydis fi 显示Y方向位移 Plot disp 显示位移矢量场,3.10 保存与调用结果命令,Call test.txt (或ca test.dat) 调用数据 Save test.sav 保存结果 New 重新开始 Rest test.sav 调用结果 Quit 退出程序,4. FLAC解题技巧,4.1 模型尺寸 4.2 模拟开挖 Mod nu i=6,15 j=5,12 (或 region i,j,i=6,15 j=5,12,4. F

14、LAC解题技巧,4.3 模拟锚杆支护（端锚） struct cable begin grid i1,j1 end grid i2,j2 seg n prop 1 stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0 kbond 0.0001 stru prop 1 sfri 30 peri 0.2723 den 7.5e-3,端锚,i2,j2,i1,j1,4.3 模拟锚杆支护（全长锚固,struct cable begin grid i,j end x,y seg n prop 2 stru prop 1 e 2e5 yield 0.5 a 0.235e

15、-3 sbond 0.42 kbond 5.37 stru prop 1 sfri 30 perimeter 0.2723 den 7.5e-3,全长锚,x,y,i,j,4.3 模拟锚杆支护（预应力锚固） struct cable begin grid i1,j1 end i2,j2 seg n ten 10 prop 1 struct cable begin node i2,j2 end x,y seg n prop 2 stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0. kbond 0.1 stru prop 1 sfri 0 perimeter 0.2723 den 7.5e-3 stru prop 2 e 2e5 yield 0.5 a 0.235e-3 sb 0.42 kb 5.37 stru prop 2 sfri 30 perimeter 0.2723 den 7.5e-3,锚固段,x,y,i1,j1,自由段,1,2,i2,j2 node n,4.4 模拟断层,Mod nu j=37 int 1 as from i1,j1 to ,i1,j2 bs fr