6数值模拟简介3D概要.ppt

FLAC3D数值模拟技术 内容提要 FLAC3D简介 FLAC3D语法 FLAC3D基本模块 FLAC3D应用实例 一、FLAC3D简介1、概述 FLAC3D(Fast Lagrangian Analysis of Continua in 3 Dimentions)是美国明尼苏达大学和美国Itasca Consulting Group Inc.开发的三维有限差分计算机程序 。 该程序主要适用模拟计算地质材料和岩土工程的力学行为 特别是材料达到屈服极限后产生的塑性流动 可计算地质类材料的高度非线性(包括应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的固—流藕合、热—力藕合以及动力学行为等。 FLAC3D是在二维FLAC程序基础上建立起来的。 可以模拟土层，岩石或其它材料承受塑性流变的三维结构行为。 用户可以根据目标物的形状来划分三维网格，这些网格单元就构成FLAC3D计算的依据。 另外， FLAC3D程序设有界面单元，可以模拟断层、节理和摩擦边界的滑动、张开和闭合行为。 FLAC还可以模拟支护结构，如砌衬、锚杆、可缩性支架或板壳等与围岩的相互作用。 同时，用户可根据需要在FLAC中创建自己的本构模型，进行各种特殊修正和补充。 FLAC3D程序是建立在拉格朗日算法基础上，特别适合模拟大变形和扭曲。 显现的Lagrangian计算系统和混合离散化数据区域技术被引用到FLAC3D中，这使得可精细地模拟塑性变形和流变 FLAC采用显式算法来获得模型全部运动方程(包括内变量)的时间步长解，从而可以追踪材料的渐进破坏和垮落，这对研究采矿设计是非常重要的。 （1）机械加载能力和变形 —— 边坡稳定性和地基设计； （2）渐变和坡断的演变 —— 用于硬岩开采和隧道设计； （3）断裂构造影响评估 —— 矿井设计； （4）地质体锚杆支护控制 —— 岩石锚固、回填和土层固定； （5）充分饱和的流体和孔隙压力建立及排水或不排水条件下的孔隙压力扩散、地下水流、固土结构研究； （6）不同材料随时间的蠕变行为—— 盐矿和钾盐矿开采设计； （7）有滑动倾向的地质构造的动力加载 —— 地震工程和矿井岩石突出研究； （8）爆炸加载和振动的动力作用 —— 隧道开挖或矿井开采设计； （9）构筑物的地震激发 —— 土坝设计； （10）因热力影响产生的荷载引起变形和机械失稳 —— 地下存储高辐射废料仓库的影响评价； （11）分析强变形材料 —— 仓中和矿石开采中材料膨胀流变研究。 FLAC模型（FLAC MODEL ） FLAC模型是用户为了模拟实际问题而创建的。 用户可使用FLAC的一系列命令，用数据方法定义问题的条件。 单元（ZONE） 有些差分单元是最小的几何空间，评估这个小区域的现象的变化（如压力与张力）， 不同形状的多面体区域（如块状、楔状、尖塔形（角锥）和四面体形状）用来创建模型，并且可用绘图工具显出来。 每个多面体区域包含一组或两组重叠的5个四面体区域。 默认情况下，两组重叠的5个四面体对解决问题精度更高。 区域的另一个术语是单元（element） 网格点（GRINPOINT） 网格点是有限差分单元各角的连接点。 根据单元形状，有５个、６个、７个或８个点与每个多面体连接。 每个网格点都有一组x、y、z坐标，这样就可以确定有限差分单元的确切位置。 网格点的另外术语是节和节点。 有限差分网格（FINITE DIFFERENCE GRID） 有限差分网格是由通过欲分析的物体范围的一个或多个有限差分单元而组成的集合。 “grid”的另一术语是“mesh（网格）”。 有限差分网格和具体位置也是固定的。在FLAC程序运行中，网格点位置储存着所有向量（如力、速度、位移）。而标量和张量（如应力和材料特征）储存在单元中心位置。 模型边界（MODE BOUNDARY） 模型边界是有限差分单元的周边界，内部边界（如网格中的洞）也是模型边界。 边界条件（BOUNDARY CONDITION） 边界条件是沿模型边界的一个约束或控制条件的规定（如力学问题中的固定位移或压力，地下水流问题中的不能渗透边界，热传导问题中的热绝缘边界）。 初始条件（INITIAL CONDIONS） 基本（或材料）模型代表FLAC模型中单元变形和承受力行为的限制。 FLAC中提供了几种基本模型来模拟地质学材料的不同类型的行为表现。 FLAC-3D程序本构模型 1)空单元模型(null model) (可用来模拟地下硐室的开挖和煤层开采 2)弹性模型（elastic model） 3)各向同性弹性材料模型(isotropic elastic model)； 4)正交各向异性弹性模型（elastic orthotrpic model） 5)Drucker-pr