岩土数值分析第五章().ppt

1. FLAC简介 1.1 FLAC特点 FLAC——Fast Lagrangian Analysis of Continua FLAC建立在拉格朗日算法基础上，采用有限差分显式算法来获得模型全部运动方程(包括内变量)的时间步长解，从而可以追踪材料的渐进破坏和垮落，这对研究采矿工程设计是非常重要的。 FLAC适用模拟计算岩土材料力学行为，特别适合模拟大变形和扭曲，包括材料的高度非线性(应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的应力—渗流耦合、热—力耦合以及动力学问题等。 1.2 岩体力学本构模型 岩石应变软化/硬化模型 断层模型 1.3 支护构件与本构模型 预应力砂浆锚杆的模拟体系 FLAC举例 (a)巴西试验 (b) 边坡滑移 （间接拉伸） 2. FLAC理论基础 2.1 有限差分法 1.6 FLAC计算流程举例 如图，一维杆件用数个等尺寸的有限差分网格划分，杆件的密度为?，杨氏模量为E。 运动方程(或平衡方程)为： 在显式算法中，所有有限差分方程右端的值均是已知的。因此，必须先用(3)式算出所有单元的应力，然后再由(5)式和(6)式计算所有结点的速度和位移。 在概念上，这个过程等价于对变量值“同时”更新，而不是像其它方法那样，方程右端混存有 “新”、“旧”值，对变量值“依次”更新。 3. FLAC软件应用 （1）设计计算模型的尺寸 （2）规划计算网格数目和分布 （3）安排工程对象（开挖、支护等） （4）给出材料的力学参数 （5）确定边界条件 （6）计算模拟 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 circle xc,yc rad 3.5 特殊形状的网格（2）弧线 gen arc xc,yc xb,yb theta 3.5 特殊形状的网格（3）直线 gen line x1,y1 x2,y2 （4）任意形状 tab 1 x1,y1, x2,y2, ?,xn,yn, x1,y1  gen tab 1 3.6 赋给单元材料性质 mod e (弹性) prop d 1800e-6 bu 12.5 sh 5.77 i=1,20 j=1,10 prop d 2400e-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 *特别提示 3.7 赋给模型边界条件 （1）固定边界 （结点） Fix x i=1, j=1,21 Fix y i=1,21 j=1 3.7 赋给模型边界条件 （2）施加边界力 (结点) apply yf -10 i=1,21 j=21 或 apply syy -10 i=1,21 j=21 apply xf -5 i=21, j=1,21 或 apply sxx -5 i=21, j=1,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 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 显示最小主应力?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