《FLAC3D数值模拟基础（中国矿业大学）》.pptVIP

FLAC3D数值模拟基础 刘升贵 中国矿业大学力学系 liushg2002@163.com 主要内容 FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并（联结）问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析 主要内容 FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并（联结）问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析 FLAC3D简介 Fast Lagrangian Analysis of Continua 美国Itasca（依泰斯卡）咨询公司开发2D程序(1986) 1990年代初引入中国 有限差分法(FDM) DOS版→2.0 →2.1 →3.0 Itasca其他软件 FLAC3D简介 应用： 岩土力学/岩石力学分析，例矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等 岩土工程、采矿工程、水利工程、地质工程 特色： 大应变模拟 完全动态运动方程使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍 显示求解具有较快的非线性求解速度 FLAC3D简介 基本原理 基本原理 有限差分法 Lagrangian网格 空间混合离散技术 Lagrangian格式动量平衡方程 FLAC3D的求解过程 FLAC3D的本构模型 有限差分法 古老的方法(上世纪40年代) 用差分格式转化控制方程中的微商格式 流体力学；土工渗流问题；固结 FDM FEM的混合求解 FDM的新进展 Lagrangian网格 源自流体力学中的拉格朗日法 跟踪流体质点的运动状态 跟踪固体力学中结点，按时步用Lagrangian法研究网格节点的运动 节点和单元随材料移动，边界和接触面与单元的边缘一致 固体力学大变形理论 FLAC3D的求解过程 FLAC3D中的本构模型 开挖模型null 3个弹性模型 各向同性弹性 横观各向同性弹性 正交各向同性弹性 8个塑性模型（Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型、修正剑桥模型和胡克布朗模型） FLAC3D中的本构模型 FLAC3D的前后处理 术语 FLAC3D的前后处理 术语 FLAC3D的前后处理 术语 FLAC3D的前后处理 术语 FLAC3D的前后处理 术语 FLAC3D的前后处理 术语 FLAC3D的前后处理 命令驱动(推荐) 程序控制 图形界面接口 计算模型输出 指定本构模型及参数 指定初始条件及边界条件，指定结构单元 指定接触面 指定自定义变量及函数(FISH) 求解过程的变量跟踪 进行求解 模型输出 菜单驱动(计算模式) 菜单驱动(Plot) 一个最简单的例子 gen zon bri size 3 3 3 ;建立网格 model elas ;材料参数 prop bulk 3e8 shear 1e8 ini dens 2000 ;初始条件 fix z ran z -.1 .1 ;边界条件 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve ;求解 app nstr -10e4 ran z 3 x 1 2 y 1 2 solve 前后处理功能的优点 多种zone类型 后处理快捷、方便、丰富 计算过程中的hist变量动态显示 FISH可进行参数化模型设计 单元状态的可编程 计算暂停时的后处理与可保存 前后处理功能的缺点 复杂模型的建模功能不强 可以编程导入其他软件形成的网格(比如：Ansys、Adina、GeoCAD) 无等值线的后处理功能(3D) 可编程将.sav文件写入TecPlot等其他后处理软件 全命令操作，学习困难 鼠标功能单一(双击取击点坐标) 后处理 常用命令 主要内容 FLAC3D软件简介 1、基坑开挖 2、浅基础的稳定性分析 3、网格的合并联结 4、界面的生成 5、隧道的生成 6、模型材料问题 7、Interface 合并（联结）问题 8、初始条件问题 9、 破坏问题分析 10、综合实例-煤巷应力分析 Case-1： 土体中挖了一个长宽为2m×4m，深为4m的沟，挖土体的同