FLAC第一讲(很好).docVIP

第一讲 FLAC技术的基本原理和应用范围 1、FLAC基本简介与本构关系 1.1 FLAC程序简介 FLAC(Fast Lagrangian Analysis of Continua，连续介质快速拉格朗日分析)Cundall和ITASCA公司开发出的有限差分数值计算程序主要适用地质和岩土工程的力学。该程序自1986年后，经不断改版，已经日趋完善。前国际岩石力学学会主席C. Fairhurst 评价它：“现在它是国际上广泛应用的可靠程序”(1994)。 根据计算对象的形状单元和区域构成相应的网格。每个单元在外载和边界约束条件下，按照约定的线性或非线性应力—应变关系产生力学响应特别材料达到屈服极限后产生的塑性流动。。由于FLAC程序主要是为岩土工程应用而开发的岩石力学计算程序，程序中包括了反映岩土材料力学效应的特殊计算功能，可解算岩土类材料的高度非线性(包括应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的固—流耦合、热—力耦合以及动力学行为等 另外，程序设有界面单元，可以模拟断层、节理和摩擦边界的滑动、张开和闭合行为。支护结构，如砌衬、锚杆、可缩性支架或板壳等与围岩的相互作用也可以在FLAC中进行模拟。此外，程序允许输入多种材料类型，亦可在计算过程中改变某个局部的材料参数，增强了程序使用的灵活性，极大地方便了在计算上的处理。同时，用户可根据需要在FLAC中创建自己的本构模型，进行各种特殊修正和补充。 FLAC程序建立在拉格朗日算法基础上，特别适合模拟大变形和扭曲。FLAC采用显式算法来获得模型全部运动方程(包括内变量)的时间步长解，从而可以追踪材料的渐进破坏和垮落，这对研究工程地质问题非常重要。FLAC程序具有强大的后处理功能，用户可以直接在屏幕上绘制或以文件形式创建和输出打印多种形式的图形。使用者还可根据需要，将若干个变量合并在同一副图形中进行研究分析。 FLAC程序中提供了由空模型、弹性模型和塑性模型组成的十种基本的本构关系模型，所有模型都能通过相同的迭代数值计算得到解决：给定前一步的应力条件和当前步的整体应变增量，能够计算出对应的应变增量和新的应力条件。注意，所有的模型都是在有效应力的基础上进行计算的，在本构关系调入程序之前，孔隙压力把整体应力转化成有效应力。下面将简要的介绍种模型的理论。 1、空单元模型 空单元用来描述被剥落或开挖的材料，其中应力为0，这些单元上没有质量力(重力)的作用。在模拟过程中，空单元可以阶段转化成不同材料，例如开挖后回填。 2、各向同性弹性模型 弹性模型的特点是卸载时材料的变形可逆，应力—应变之间保持线性的变化规律，和加载路径无关。 各向同性弹性模型对材料的力学行为提供了最简单的描述。它适用于只产生线性应力—应变关系，在卸载条件下没有产生滞后效应的均质材料、各向同性材料和连续介质。 在平面应变条件下，模型的应力—应变增量表达式由胡克定律得： (6-44) () 式中：； ； K——体积模量； G——剪切模量。 (6-45) 式中：——应变张量增量； ——位移量； ——时步。 在平面应力条件下，模型的应力—应变增量表达式为： (6-46) () 式中：； 。 3、横观各向同性弹性模型 横观同性模型能够模拟层状的弹性介质，模型中垂直于层状介质的弹性模量和平行于介质的弹性模量明显不同。弹性模型的各向同性平面位于X-Z平面内，弹性模量的相关定义如下： E1 ( Ex )——各向同性平面内的弹性模量； E2 ( Ey ) ——垂直于各向同性平面内的弹性模量； G12 ( Gxy) ——X-Y平面或者Y-Z平面内的剪切模量； G13 ( Gxz)——各向同性平面内的剪切模量； ()——由单轴Y方向上的应力引起的各向同性平面内X方向上的法向应力和垂直于各向同性平面上的Y方向上的法向应力组成平面内的泊松比； ()——由于单轴Z方向上的应力引起的各向同性平面内X方向上的法向应力和垂直于各向同性平面上的Z方向上的法向应力组成平面内的泊松比。 横观各向同性弹性材料的特征是包含五个独立的常量。弹性模量参数之间的具体关系： (或) (或) (或) (或) (或) (或) 图6-6 横观各向同性坐标系定义规则(X-Z方向为各向同性面) 模型中关于材料的弹性模量变化有一定的限制(Amadei,1982)，其限制条件为： (6-47) 4、Drucker-Prager(德鲁克-布拉格)塑性模型 Drucker-Prager塑性模型适用