FLAC3d与有限元程序联合建模.pptVIP

FLAC3D与有限元程序 上海大学：王栋 1、前言——研究问题的必要性 2、FLAC3D算例介绍 3、FLAC3D与有限元程序的联合使用 4、展望 1、问题研究的必要性 FLAC3D简介 FLAC是快速拉格朗日差分分析（Fast Lagrangian Analysis of Continua）的简写。FLAC3D是由美国Itasca Consulting GroupInc开发的以三维快速拉格朗日法为分析程序,应用于岩土工程力学计算的应用软件. 与有限元程序比较： 优点：（1）FLAC3D采用混合离散方法来模拟材料的屈服或塑性流动特性，这种方法比有限元方法中通常采用的降阶积分更合理；（2）FLAC3D利用动态的运动方程进行求解，使得FLAC3D能模拟动态问题，如振动、失稳和大变形等；（3）采用显式有限差分格式来求解场的控制微分方程。 显示差分法与有限元法对比 缺点：（1）对于线性问题，FLAC3D要比有限元程序花费更多计算时间，因而，FLAC3D在模拟非线性问题、大变形问题或动态问题时更有效；（2）FLAC3D的收敛速度取决于系统的最大固有周期与最小固有周期的比值，这使得它对某些问题的模拟效率非常低，如单元尺寸或材料弹性模量相差很大的情况；（3）对于复杂模型的建立不仅工作量大而且较难掌握。 总结： FLAC3D有着非常强大的计算功能,能够较好地模拟材料的力学行为,尤其是在岩体大变形方面,比其他有限元软件更加准确快速.但目前,FLAC3D建立计算模型时采用命令输入式的方法,虽然FLAC3D提供了12种初等模型,但是对于建立复杂模型来说,这些远远不够. 研究如何将具有强大建模功能的有限元程序与强大计算功能FLAC3D相结合有很大的必要性。 2、FLAC3D实例 问题：在任意一土体中开挖一2m×4m×3m的沟渠，进行应力、应变场分析？ S1：创建初始几何模型、划分网格 在FLAC3D提示符下输入命令：generate zone brick size 6 8 8 (该命令将产生在x轴6格、y轴8格和z轴8格的三维长方体网格——在初始 模型建立的同时划分好了网格) S2：定义材料模型及参数 在FLAC3D提示符下输入命令：model mohr (MODEL命令定义材料模型——mohr参数指摩尔-库仑模型) 在FLAC3D提示符下输入命令： property bulk =1e8 shear =0.38e8 friction=35 property cohesion=1e10 tension=1e10 (PROPERTY命令定义本构模型材料参数) bulk——体积模量 ；shear——切变模量 ； friction——内摩擦角 ； cohesion——内聚力 ；tension——抗拉强度 。 S3：加载及边界条件 在FLAC3D提示符下输入命令：set gravity 0, 0, -0.98 (SET命令用于设置FLAC3D的模拟条件或控制条件 ; 0.98指重力加速度值) 为了开挖沟渠，还需要初始化模型中土的密度 输入如下命令：initial density = 1000 （INITIAL命令初始化网格相关值； density 为网格质量密度--kg/m3） 接下来，设置边界条件： 输入如下命令： fix x range x -0.1 -0.1 fix x range x 5.9 6.1 fix y range y -0.1 -0.1 fix y range y 7.9 8.1 fix z range z -0.1 0.1 FIX命令的功能是保持网格节点指定参数（速度、压力和温度）的值不变 range 为指定范围 。 第一行的命令意思：在范围x= -0,1和x=0.1 这两个平面内所有的网格节点 保持在x轴方向的速度不变，因默认初速度为0，故位移也为0. 以上命令作用：对几何模型6个边界面中的5个固定。 S4：求解 监控变量 本例监控两个参数，一个是点（4， 4， 8）在 z 方向的位移迭代变化，另 一个是模型中最大不平衡力。命令如下： history nstep = 5 history unbalance history gp zdisplacement 4, 4, 8 HISTORY命令是采样（或记录）迭代时模型中变量的值 ；