第三章、有限元基本理论.doc

有限元基本理论 摘要：从一般的边值问题数值解理论出发，讲解了有限元法的基本过程和基本理论。有限元法基本过程包括问题几何区域的离散、近似解待定参数的确定、方程的建立等；基本理论包括单元的分类、单元形函数的性质、等参单元、单元积分和节点等。本章讲述的内容不受应用领域的限制。 有限元法是为了解决结构分析而发展起来的一种新的数值方法。经过近50年发展，它不但是结构分析强有力的工具，而且，在结构分析获得重大成功后，其理论也已日趋成熟，商务化软件系统也已有一定规模和数量，在其它领域边值问题的数值计算方面同样获得巨大成功。 设由边界围成区域，其基本解为未知函数的某一连续介质边值问题。在第一章中我们将此问题转化成等效积分形式，并用加权残数法进行数值解；第二章中对具有泛函极值形式的问题采用Litz进行数值解。但是以上两章并没有解决数值解中的试探函数（有限元中称形函数）的选取问题。 有限元方法的关键是待定参数和形函数的选取及计算，那么采用有限元数值解法，需要经过哪些基本理论和过程呢？ 有限元法概述 区域的离散化 将区域近似地离散成有限数量的，基本形状有一定限制的，尺寸远小于和的子区域集，称为有限单元（Element）集，它的元素称单元，记为或，对每个单元给予编号，即 (3.1.1) 单元的基本形状可根据的几何维数选择，例如一维几何区域为线单元；二维区域可选择三角形或四边形单元；而三维区域选择四面体、五面体和六面体单元等。图1.1的平面区域被离散成有限个三角形单元，详细的单元分类和性质请见3.3的讨论。控制单元形状和位置的点称为单元节点（element node,也有称结点或接点），简称节点（Node），例如图1.2。节点的集合记为，称节点集，并给予编号，即 (3.1.2) 围成单元的几何元素称为单元边界，例如图1.2中四边形单元的四条边（edge）、四个顶点节点和四个中间节点都属于单元边界。单元边界比之单元在几何维数上要低，根据几何维数不同，单元边界又可以分单元面、单元边、单元节点。在离散区域时，为了保证问题解的唯一性和连续性，两相邻单元的边界必须保持完全重合，即单元边界的节点被相邻单元完全共享。例如图1.3中的节点1被单元与共享，而节点2被与四个单元共享。如何保证单元之间的问题解的连续性将在3.3、3.4节中讨论。 确定待定参数集 在第一章中已经指出，边值问题的数值解是待定参数矢量集的线性组合 (3.1.3) 设节点的问题解的值为，组成的集合记为，即 (3.1.4) 虽然还不能完全等同近似解的待定参数集，但如果试探函数看成是对插值函数，从矢量运算角度考虑，(3.2.3)可以改写成和 (3.1.5) 其中试探函数（插值函数）在有限元中称为形函数（shape function），所以在得到后，就获得了问题的近似解，只是选定合适的形函数。 例如，图1.4由四个四边形单元组成固体力学平面应力应变问题，则由所有单元节点位移矢量所组成，简称位移矢量。所以问题的单元集、节点集和位移矢量分别为 在中，并不是所有参数是待定的。在本质边界上，节点的值是确定，在混合边界上，节点的受到边界条件方程的约束。例如固体力学问题位移解法中，位移边界上节点的位移值属于已知，混合边界上节点的位移受混合边界条件方程约束。但是不管节点的值如何获得，(3.1.5)的近似式仍然成立。所以在有限元方法对单元讨论，暂时把看成待定参数集，只是在后面求界待定参数方程组时，把已知的参数和约束方程代入方程组，从而减少方程组的数量，详细讨论见下章讨论。 单元形函数的基本要求 在单元中，设有个节点。为了分析方便，节点的编号仍然采用1至，称之为局部编号以区别节点的整体编号。记第个单元局部节点的问题解在有限元法中采用以下假定： 单元内的问题解近似值只是该单元节点问题解的值所决定，与其他单元问题解的值无关。 问题解的每个分量都采用相同的形函数。 所以单元内近似解的插值的矢量形式和分量形式为 (3.1.6) 其中为单元内问题解的第个分量，为单元的第个节点的形函数，是第个节点的的第个分量。 以上插值显然是Langrange插值法，只保证了近似解的阶连续。如果要提高问题解连续性阶数，则需采用Hermite插值法，这时以上第一条假定得取消。 为了保证问题解的唯一性和单元之间问题解的连续，（3.1.6）式形函数必须满足以下性质： 1）唯一性：在每个节点上插值函数的值有 (3.1.7) 2）连续性：单元边界（或是单元面，或是单元边，或是单元节点）上的形函数值，除了此边界上节点的形函数外，其他节点的形函数必须为0，即 (3.1.8) 单元边界可以是。例如图1.5三维20节点的六