(精)第三章 有限元基本——最新.ppt

有限元分析岩土工程数值计算 主讲：翁其能 2010年10月 第 三 章 有限元基本 3.1 概述 1.有限元法（Finite Element Method） 简称FEM，是弹性力学的一种近似解 首先将连续体变换为离散化结构，然后再利用分片插值技术与虚功原理或变分方法进行求解。 2. FEM特点 （1）具有通用性和灵活性。 （2）对同一类问题，可以编制出通用程序，应用计算机进行计算。 （3）只要适当加密网格，就可以达到工程要求的精度。 3.1 概述 FEM简史 FEM是上世纪中期才出现，并得到迅速发展和广泛应用的一种数值解法。 1943年柯朗第一次提出了FEM的概念。 1956年，特纳等人提出了FEM。 20世纪50年代，平面问题的FEM建立，并应用于工程问题。 1960年提出了FEM的名称。 3.1 概述 FEM简史 20世纪60年代后，FEM应用于各种力学问题和非线性问题，并得到迅速发展。 1970年后，FEM被引入我国，并很快地得到应用和发展。 有限单元法的物理概念清晰，易于掌握和应用，计算速度快，精确程度高，具有灵活性和通用性，可以解决一些复杂的特殊问题，例如复杂的几何形状，任意的边界条件，不均匀的材料特性，结构中包含杆件、板、壳等不同类型的构件等。近二、三十年来，广泛应用于航空、造船、土木、水利、机械工业中。 3.2有限元法基本原理 基本思想是用有限个离散单元的集合体代替原连续体，采用能量原理研究单元及其离散集合体的平衡，以计算机为工具进行结构数值分析。 有限元模型是真实系统理想化的数学抽象。 材料的响应可以用状态变量描述。 位移（场） 应力（场） 应变（场） 一般地，状态变量是连续函数，求得状态变量解析解需要求解微分方程，这对于复杂问题是不可能的。 3.2有限元法基本原理 单元分析： 对于弹性力学问题，单元分析，就是建立各个单元的节点位移和节点力之间的关系式。 每个三角形单元仍然假定为连续的、均匀的、各向同性的完全弹性体。因单元内部仍是连续体，应按弹性力学方法进行分析。 取各结点位移 为基本未知量。然后对每个单元,分别求出各物理量,并均用 来表示。 结点位移 结点力 整体分析 对由各个单元组成的整体进行分析，建立节点外载荷与结点位移的关系，以解出结点位移，这个过程为整体分析。 在位移法中，主要的任务是求出基本未知量---结点位移。为此需要建立结点的平衡方程。 1.7.2 整体刚度矩阵 假设整体结构被划分为ne个单元和n个节点，在整体坐标系下，对于每个单元均有： 式(1-49）表明：ij边受线性分布面力： i点为[qsx, qsy]T，j点为0 时，其等价节点力可将总载荷的2/3分配给i点，1/3分配给j点，m点为零得出。 x y i j m qsi qs 体积力和表面力向节点的移置符合静力等效原理的前提条件是：线性位移模式。 现在，只考虑弹性体边界上的表面力，它只在部分单元上形成表面力（右下图）。设边界面上单位面积受到的表面力如下式： l—单元边界长度 h—单元厚度 A—表面力作用面积 ① ② ③ ④ ?qs? ?qs?沿厚度均匀分布，则单元表面力的势能Vs为 （3） 集中力势能 当结构受到集中力时，通常在划分单元网格时就把集中力的作用点设置为节点。于是单元集中力?Pc?的势能Vc为 p ① ② ③ ④ ③ ??? ③ p/2 ??C? （4）总势能 把(3-35）式中原括符内的部分用列阵?Fd?代替， 综合以上诸式，单元外力的总势能V为 (3-35） ?Fd?具有和???相同的行、列数。则 (3-36） 由单元的应变能U(3-34）和外力势能V(3-36），可得单元的总势能? (3-37） 将式(3-37）代入， 根据弹性力学最小势能原理：结构处于稳定平衡的必要和充分条件是总势能有极小值。 3、单元平衡方程 于是有， (1-34） (1-36） 式(3-38）是从能量原理导出的单元平衡方程。这个方程表达了单元力与单元位移之间的关系。其中，?Fd?和单元节点力?F?具有相同的意义。 (3-38） 即得单元平衡方程 1.6 单元刚度矩阵 平衡方程(3-38）中的矩阵[k]是单元力和单元位移关系间的系数矩阵，代表了单元的刚度特性，称为单元刚度矩阵。单元刚度矩阵的体积为nj ×nj， nj 是单元位移总数。其一般计算公式为： 1、一般计算公式 它与单元应变矩阵[B]和弹性矩阵[D]有关。 (3-33）