FLAC3D简述与使用步骤.docxVIP

FLAC 3D简述与使用步骤 （附：软件+实例）非常高兴看到 Matlab中文论坛 里有FLAC 3D的版块，本人在地质院校学习。乐意与大家一起学习这个非常好的岩土软件，同时需要说明的是它的应用绝对并不局限与岩土工程。下面是我自己在学习过程，通过书籍、论坛资料总结的一些资料。希望对大家有帮助。FLAC是快速拉格朗日差分分析（Fast Lagrangian Analysisof Continua），相信学过数学都应该听说过（拉格朗日^_^）。它最早由Willkins用于固体力学领域。FLAC 3D程序是美国ITASCA公司的产品，目前已经成为岩土力学计算中的重要数值方法之一。它是对FLAC二维的扩展，广泛应用于边坡稳定性评价、地下洞室、隧道工程、矿山工程、甚至是核废料处理。优点：（1）采用混合离散方法模拟材料的屈服或塑型流动特性，比起有限元更合理。（2）它的求解过程是动态的，能够模拟动态问题。（3）FLAC 3D是基于显式方法进行求解的，方便求出应力增量。不足：1 对于线性问题的求解，FLAC3D比其他有限元程序运行得要慢；但是，当进行大变形非线性问题或模拟实际可能出现不稳定问题时，FLAC3D是最有效的工具。 2 用FLAC3D求解时间取决于最长的自然周期和最短的自然周期之比。1FLAC 3D 基本术语Zone: 有限差分单元体，最小的几何区域。它的形状有长方体、椎体、四面体等。Grid Point：网格节点，连接有限差分单元体的交点。一个单元体由每个节点的坐标确定。Finite Difference Grid：有限差分网格，由多个单元体组成。可覆盖问题所分析的物理区域。Model Boundary：模型边界，若网格内部出现孔，其边界称为内部边界。Boundary Condition：边界条件。它规定对模型边界的约束或者控制条件。Initial Conditions:初始值。Constitutive Model：本构模型，模拟单元体的变形或强度的特性。Null Zone：空单元。Sub-Grid：子网格。网格可由子网格组成。Attached Faces：绑定面：是子网格被绑定和连接在一起的面。Interface：分界面，各类网格的分界面。Range：范围，用于描述三维空间范围。Group：组，是唯一命名的一组单元体。用来限制命令的范围。ID Number：ID号，模型中的任一元素都用ID标识。包括：节点、单元体、表、绘图等。Structural Element：结构元素，软件中有两种：两节点和三节点。Step：步，用于控制计算的步数。Static Solution：静态解。Unbalance Force：不平衡力。Dynamic Solution：动态解。Large Strain/Small Strain:用于指定变形的大小。2.FLAC程式使用前准备步骤 步骤1：依比例画出所欲分析的资料在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。 步骤2：换算输入资料成同一单位将现有地层资料，如 Density, Bulk modulus, Young`s modulus, tension, cohesion, frictionAngle 等资料，换算成同一单位。　附注 ：需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会造成FLAC无法计算，而产生ERROR信息。　 　　　　步骤3：应用公式简略计算应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。步骤4：建立x, y坐标与node i, j之间的关系在图上距离及深度的关系，建立x, y坐标系统，再由x, y坐标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格的位置及其相对应的x, y 坐标。建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行程序时，不宜使用网格太大的网格数目，因尽量使网格总数少于1000，以节省时间。2.FLAC 3D 程序的编写步骤：1 Config　________2 Grid　　________3 Model 　________4求起始的应力平衡 (1)建立x, y坐标与网格的关系，建议使用Gen指示：Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 详细指令参见使用手册，FLAC程式可自动产生x, y 坐标与网格的关系，但由于产生的网格坐标不易控制，将对其它的工作产生负面影响，故依使用前步骤4所建立的关系，将网格依其疏密程度需要的不同，实际控制网格的坐标。 (2)设定材料性质：prop (3)设定外力：Set Grav, Apply Pressure, inisxx, Syy(4)设定边界条件：fix, free(5)求起始