高等岩石力学岩石力学有限元法讲解.ppt

高 等 岩 石 力 学 第四章 岩石力学问题有限元法 岩石力学有限元法 第一节 概论 第二节 施工建造过程的模拟 第三节 节理及不连续面的模拟 第四节 多节理岩体的模拟 第五节 无限域单元及其应用 第六节 岩体工程中的弹塑性问题 第七节 无拉力分析及节理非线性分析 第八节 岩土工程三维非线性有限元程序系列 第一节 概 论 一、有限元法在工程中的应用 有限元自50年代发展至今，以成为求解复杂的岩石力学及岩土工程问题的有力工具，并已为愈来愈多的工程科技人员所熟悉。在求解像弹塑性及流变,动力，非稳态渗流等时间相关性问题，以及温度场在求解像弹塑性及流变、动力、非稳定渗流等时间相关性问题，以及温度场、渗流、应力场的耦合问题等复杂的非线性问题中的效能已使它成为在岩石力学领域中应用最广泛的数值分析手段。特别是近十余年，在工程应用方面已有了较大的进展，并引起广大工程科技人员的兴趣。在岩土工程有关专业的大学生和研究生中，有限元已被列为独立的课程。在本章中将着重讨论有限元法在岩石力学中应用的有关问题。涉及有限元法基本原理，方法及求解技术方面的基础知识。 有限元法基本方程 就数学概念来说，有限元法是通过变分原理(或加权余量法)和分区插值的离散化处理把基本支配方程转化为线性代数方程，把求解待解域内的连续场函数转化为求解有限个离散点（节点）处的场函数值。显然，这种离散化的处理是一种近似，因面中有当单元划分得充分小时，才能保证满意的求解精度。 由于单元充分小，在一个微小的单元内，未知场函数可以采用十分简单的代数多项式近似地表述。通常可取为如下的插值形式： (4-1) 下标m表示单元的节点数目，下标e表示单元的序号。 有限元法即是以所有节点处的1值作为基本未知量。对于二维的工程物理问题，如热传导、渗流等，其基本支配方程为如下形式的准谐方程： (4-2) 边界条件为： ① ，边界 在上 ② 边界在 上 当 =常量时，上式即退化为波松方程。由变分原理可知，满足方程式（4-2）及相应边界条件的场函数，应使如下的泛函数驻值，即： (4-3) 在对被考察域A进行离散化的情况下，此泛函驻值转化为多元函数的驻值： （4-4） 对于一个单元体的子域内，把方程式（4-1） 带入方程式 （4-3)并进行变分后即可得到单元的基本方程： (4-5) 式中 即单元特性矩阵： 即为节点未知量组成的矢量： 是由在 上的边界条件给定的a（x,y）有关的矩阵，故仅对a≠0的二类边界上的单元才予考虑。对整个求解域A由变分驻值条件式（4-4）可得 或写为 (4-6) 式中，[k]为由各单元特性矩阵 及 按节点号组集得到的总体特性矩阵：{u}为所有节点的待求值组成的矢量：{p}为与Q（x,y）有关的矢量及q(x,y)有关的矢量。 解线性方程式（4-6）即可求得场函数 在各单元节点处的值。 对于弹性力学问题，可通过最小势能原理或虚功原理导出有限元法的基本方程。有限元法求解弹性力学问