岩土数值2014绪论.ppt

岩土数值分析－有限元法在岩土工程中的应用 学习要求： （1）课堂认真听讲，把我强调的重点在书上标记出来； （2）课程结束前至少3个星期开始复习； （3）复习时，一定要认真看书及ppt，其它资料只能参考，不能作为重点； （4）希望每一位同学一定要高度重视，是严格的闭卷考试，学校教务处统一负责的，监考很严格。 §1.概述 朱伯芳《有限单元法原理与应用》：有限单元法是在六七十年代发展起来的强有力的数值分析方法，它使许多复杂的工程问题迎刃而解，而且由于前后处理技术的发展，计算效率非常高，实际应用越来越广泛。 有限元（Finite Element Method Finite Element Modelling --- FEM Finite Element Analysis FEA) 一、有限元简介 §1.概述 有限元方法可以用来求解多种问题：这里应力变形问题 有限元方法把系统（结构）看作是由无限多个单元组成的连续体，在解这一连续体时将连续体离散化，然后将物理方程、平衡方程、几何方程结合起来，变换为求解线性方程组问题。单元与单元之间只通过结点连接 有限元分析可概括为六个步： 有限元基本思想 §1.概述 物理方程 有限元基本思想 平衡方程 几何方程 最终建立外荷载与位移之间的关系 §1.概述 有限元分析可概括为六步： 1. 结构的离散化：将分析结构系统划分成有限个单元体，并在单元体的指定点设置节点，把相邻的单元体在节点处连接起来组成单元的集合体，以代替原来的结构。一般情况，单元划分越细则表述变形情况越精确，即越接近实际变形，但计算量越大。 2. 选择位移模式：位移模式，就可导出用节点位移表示单元内任意点位移 有限元基本思想 §1.概述 3. 建立单元结点力和结点位移之间的关系 有限元基本思想 §1.概述 4. 计算结点荷载(包括集中力、面力、体积力) 5. 集合所有单元的刚度方程，建立整个结构的平衡方程 有限元基本思想 6. 方程组求解 直接法：高斯消去法、三角分解法、波前法等 迭代法：计算快、存储小 （1）用离散结构代替原结构：离散结构由许多 单元体组成，各单元只在结点处有力的联系。 （2）荷载简化到结点上，{R} （3）由s～ε→[D]→[K] （4）解线性方程组 [K]{δ}={R} （5）由位移{δ}→应变{ε}→应力{s} （6）对非线性问题，需重复（3）-（6）步。 有限元基本步骤 非线性：材料非线性与几何非线性 §1.概述 （1）平面（应力、应变）有限元：三结点、四结点、6结点、8结点 （2）轴对称有限元(本质上是二维问题) （3）三维有限元 有限元类型 （1）总应力法 （2）有效应力法 §1.概述 有限元应力变形计算时，无论采用哪种有限元软件，建立有限元模型是第一步必需做的工作。 一般说来，有限元模型应包含以下数据信息：（1）建立合适坐标系，将结构离散后给出所有结点编号及结点坐标；（2）所有单元结点编号信息；（3）各单元对应的材料信息，选择本构模型及材料参数（包括重度、本构模型参数等）；（4）约束结点的信息；（5）结构上的全部荷载及荷载分级信息，以及给出方程求解方法。 §1.概述 （6）对考虑固结的有效应力分析，还应该把结构的排水边界、孔压边界、每级荷载所经历的时间、计算时间步长、渗透性参数等告诉程序。 因此，对初学者，进行有限元计算时应检查上述（6）方面的基本数据是否已经具备，缺少任何一种数据信息都不能正常运行或得不到合理的结果。 1、有限元程序 Abaqus、adina、MAC、sap2000、ansys、 COMSOL Multiphysics Madas、plaxis、 2、有限差分 Flac 3、离散元 PFC 数值分析商业程序 运行方式： 1、窗口人机交互窗口模式—abaqus/CAE 建模、提交任务、监控运算过程、计算结果后处理 2、利用有关软件生成inp文件，在DOS状态下运行，如x.inp 生成：x.dat，x.odb abaqus job=x user=e:\temp\umat3d.for interactive §6 abaqus应用 inp文件 §6 abaqus应用 ** 1个单元 *node,nset=dam 1, .000, 5.000, 0.000 2, .000, 5.000, 5.000 3, 0.000, 0.000, 5.000 4, .000, 0.000, 0.000 5, 5.000, 5.000, 0.000 6, 5.000,