有限元法基本原理及应用第3章重庆大学龙雪峰.pptVIP

有限元法作为一种数值分析的方法现已被广泛应用于固体力学、流体力学、传热学、电磁学等多个领域。 这种方法最早萌芽于弹性力学问题的求解。以位移、力或者力和位移的组合为未知量分别把有限元方法分为有限元力法、有限元位移法和混合有限元法。 有限元位移法由于计算过程的系统性和通用性而被广泛应用，目前大多数商业有限元软件在对弹性力学问题的分析中大多采用位移法。 本章介绍求解弹性力学问题的有限元位移法。 第3章 弹性力学有限元法 3.1 有限元法求解问题的基本步骤 有限元法作为一种分析方法，即可以分析比较简单结构的问题如一根轴、一块板等，也可以分析比较复杂结构的如潜艇、体育场馆等，但其分析问题的基本步骤大同小异，可归纳如下： 1．建立力学模型 首先，在尽可能反映工程实际的基础上对工程结构进行几何简化，一般分为一维问题、二维问题和三维问题。其中二维问题要分清是平面应变问题还是平面应力问题，如果结构具有对称性可利用来简化计算。 其次，要把结构上作用的载荷根据作用区域的大小简化为集中力、线载荷、面载荷或体力等。 第三，根据边界约束和运动情况提出边界条件。 2．连续体离散化 根据所建立的力学模型选用合适的单元将连续体划分为有限个具有规则形状的单元集合。 单元的选取应视所分析问题的性质、规模和精度要求而定，如杆单元有二维单元、三维单元，二维问题可选三角形单元、四边形单元，三维问题有四面体单元、六面体单元等。 3．单元分析 单元分析包括位移模式选择，单元力学分析两个内容。 位移模式也称位移函数或插值函数，在有限元位移法中是以节点位移为基本未知量，再由这些节点位移插值得到单元内任意一点的位移值。单元的位移模式一般采用多项式，因为多项式计算简便，并且随着项数的增加，可以逼近任何一段光滑的函数曲线。 单元力学分析 根据所选单元的节点数和单元材料性质，应用弹性力学几何方程和物理方程得到单元刚度矩阵。由于连续体离散化后假定力是通过节点在单元间传递的，因此要利用插值函数把作用在单元上的体积力、面积力和集中力按静力等效原则移到节点上。 4．整体分析和有限元方程求解 由已知的单元刚度矩阵和等效节点载荷列阵组装成整个结构的整体刚度矩阵和载荷列阵，得到一个由总体刚度矩阵[K]、总载荷向量{F}和整体节点位移向量{δ}表示的平衡方程式：[K] {δ}＝{F}。引进位移边界条件后求解得到整体节点位移向量。 有限元离散方程是一个代数方程组，代入边界条件处理以后的刚度矩阵是一个正定的对称稀疏方阵，这样一个代数方程组可以用高斯消元法、三角分解法、波前法和雅可比迭代法等多种方法求解。 5．结果后处理和分析 求解线性方程组得到位移矢量后，由几何和物理关系可以得到应变和应力。 由于应变（应力）来自位移的微分可能导致单元间应力不连续，这会使应力计算误差较大，要在节点附近进行平均化处理。 通过后处理还可得到位移、应变和应力的最大最小值及其所在位置以及主应力、主应变或其它定义的等效应力。 结果的输出可以应用图表、动画等各种方式。最后还要对这些结果进行分析以指导工程设计、产品开发等等。 3.2 连续体离散化 一个结构经过力学简化以后通常变为梁杆、板壳、实体或它们的组合，要进行有限元分析必须进行离散化，即建立有限元模型，就是把结构用不同的单元划分为只在结点上彼此连接的有限个单元组成的离散体。 为了使有限元模型能够准确地体现实际结构，必须选择适当的单元类型，下面简要介绍几种常用的单元类型。 3.2.1 杆状单元 平面杆单元 b) 空间杆单元 图3.1 杆单元 平面梁单元 b) 空间梁单元 图3.2 梁单元 平面桁架、空间网架是工程中常见的结构，这种结构的主要受力构件是杆件，其截面尺寸远小于轴向尺寸。这类单元主要有平面杆单元和空间杆单元，如图3.1所示。 杆单元有两个节点，只能承受轴向的拉压载荷，平面杆单元每个节点有两个自由度，空间杆单元每个节点有三个自由度。杆状单元还有用于求解平面刚架和空间刚架问题的梁单元，如图3.2所示。 常用梁单元有两个节点， 平面梁单元每个节点有三个自由度，两个线位移和一个角位移，可承受平面内的体力、集中力、分布力和垂直平面的弯矩的作用，只考虑单元的拉压和一个方向的弯曲变形。 空间梁单元每个节点有六个自由度，三个线位移和三个角位移，一般可承受各个方向的集中力、分布力和弯矩，能够考虑单元的弯曲、拉压、扭转变形。当切应变对挠度的影响比较大时，必须修正单元刚度矩阵以考虑切应变。 3.2.2 平面单元 平面问题是三维问题根据结构和外力特点，近似简化的结果。这种简化不仅为结构的理论分析带来了很大的便利，也为结构的有限元分析带来了便利。 弹性力学平面问题根据其应力与应变特点分为平面应力问题和平面应变问题，在进行有限元分析