2010-11-12_校庆报告-FLACFLAC3D理论与实践-陈育民.pdf

河海大学校庆学术报告 岩土工程的有限差分数值模拟 ——FLAC/FLAC3D理论与实践 陈育民 博士 河海大学土木与交通学院 安全与防灾工程研究所 2010.11.12 提纲 有限差分数值模拟基本理论 静力分析与前后处理 接触面 FISH语言 流固耦合分析 动力分析 自定义本构 结构单元 FLAC(2D)基本应用 2 提纲 有限差分数值模拟基本理论 静力分析与前后处理 接触面 FISH语言 流固耦合分析 动力分析 自定义本构 结构单元 FLAC(2D)基本应用 3 基本原理 FLAC/FLAC3D利用有限差分，显示方案，动态松弛方法 模拟连续体的非线性力学行为： 即使对准静态问题, 程序仍然求解完整的动力学方程。这种方法 的好处在于可以为物理非稳定过程例如塌方提供稳定解； 在“松弛”方法中，使用阻尼来吸收动能以模拟系统的“静 态”反应。这种方法可以用比其它方案如解矩阵法更为真实有 效地模拟塌方问题。 4 Lagrangian法 源自流体力学中的拉格朗日法 跟踪流体质点的运动状态 跟踪固体力学中结点，按时步用 Lagrangian法研究网格节点的运动 节点和单元随材料移动，边界和接 触面与单元的边缘一致 固体力学大变形理论 法国数学家、物理学家拉格朗日 5 混和离散技术 FLAC混和离散 = + /2 每个 为常应力/应变: 体积应变由整个四边形算出 . 应变偏量则有两个三角形 和 分别算出 (混合离散 过程) 解题过程中网格坐标按照“拉格朗日方式更新” ( 网格随材料移动), 且为显式 (一个时步内局部变化不会影响邻域) 6 FLAC3D混和离散 = + /2 7 FLAC3D混和离散 结构域离散为可由四面体单元组合形成的五面体或六面体等 单元; 以 为基本单元(常应力、常应变); 体应变的计算： ; 偏应变的计算： . 8 三维快速Lagrangian原理 空间导数