有限单元法概述.ppt

(3)有限元方程的解法 现在用于大型复杂工程问题的有限元分析，自由度达几十万个甚至上百万个已是经常的情况，这是与计算机软、硬件发展相配合的大型方程组解法的研究进展密不可分。有限元求解的问题从性质上可以归结为三类，即: ①独立于时间的平衡问题(或稳态问题) 最后归结为求解系数矩阵元素在对角线附近稀疏分布的线性代数方程组。此类问题至今主要是采用直接解法，先后发展了循序消去法、三角分解法、波前法等。近年来，为了适应求解大型、特大型方程时减少计算机存储和提高计算速度的需要，迭代解法特别是预条件共轭梯度法受到更多的重视，并已成功地应用。 ②特征值问题 它对应求解的是齐次方程。解答是使方程存在非零解的特征值和与之对应的特征模态；在实际应用中，它们代表的可能是振动的固有频率和振型，或是结构屈曲的临界载荷和屈曲模态等。针对求解经数值离散所导致的大型矩阵特征值问题，先后发展了幂迭代法、同步迭代法、子空间迭代法等。近10多年来，里兹(Ritz)向量直接叠加法和Laczos向量直接叠加法由于具有更高的计算效率而受到广泛的重视和应用。 * 有限单元法概述 是在当今工程分析中获得最广泛应用的数值计算方法。 由于它的通用性和有效性，受到工程技术界的高度重视。 伴随着计算机科学和技术的快速发展，现已成为： 计算机辅助设计(CAD)—— Computer Aided Design 计算机辅助制造(CAM)—— Computer Aided Manufacture 的重要组成部分。 有限单元法(或称有限元法，FEM )—— Finite Element Method 在工程或物理问题的数学模型: 基本变量; 基本方程; 求解域和边界条件等 确定以后，有限元法作为对其进行分析的数值计算方法的要点可归纳如下： 1 有限元法的要点 (1)将一个表示结构或连续体的求解域离散为若干个子域(单元)，并通过它们边界上的节点相互联结成为组合体。 下图表示将一个二维多连通求解域离散为若干个单元的组合体。图(a)和(b)分别表示采用四边形和三角形单元离散的图形。各个单元通过它们的角节点相互联结。 (2)用每个单元内所假设的近似函数来分片地表示全求解 域内待求的未知场变量。 每个单元内的近似函数由未知场函数(或其导数)在单元各 个节点上的数值和与其对应的插值函数来表达(此表达式 通常表示为矩阵形式)。 由于在联结相邻单元的节点上，场函数应具有相同的数 值，因而将它们用作数值求解的基本未知量。 转换为 求解原来待求场函数的无穷多自由度问题 求解场函数节点值的有限自由度问题 (3)通过和原问题数学模型(基本方程、边界条件)等效的变分原理或加权余量法，建立求解基本未知量(场函数的节点值)的代数方程组或常微分方程组。此方程组称为有限元求解方程，并表示成规范化的矩阵形式。接着用数值方法求解此方程，从而得到问题的解答。 (1)对于复杂几何构形的适应性 由于单元在空间可以是一维、二维或三维的，而且每一种单元可以有不同的形状，例如三维单元可以是四面体、五面体或六面体，同时各种单元之间可以采用不同的联结方式，例如两个面之间可以是场函数保持连续，可以是场函数的导数也保持连续，还可以仅是场函数的法向分量保持连续。这样一来，工程实际中遇到的非常复杂的结构或构造都可能离散为由单元组合体表示的有限元模型。 2 有限元法特性 下图所示是一水轮机转轮的有限元模型: 转轮由上冠、下环和13个叶片组成; 分别用三维块体单元和壳体单元离散; 叶片之间的水用三维流体单元离散。 (2)对于各种物理问题的可应用性 由于用单元内近似函数分片地表示全求解域的未知场函数，并未限制场函数所满足的方程形式，也未限制各个单元所对应的方程必须是相同的形式，所以尽管有限元法开始是对线弹性的应力分析问题提出的，很快就发展到弹塑性问题、粘弹塑性问题、动力问题、屈曲问题等。并进一步应用于流体力学问题、热传导问题等。而且可以利用有限元法对不同物理现象相互耦合的问题进行有效的分析。 下图表示金属板料成形过程的有限元模拟 其中图(a)表示冲头、模具和板料的图形； 图(b)是它们的有限元模型； 图(c)，(d)，(e)是冲头向下移动20mm、30mm、40mm时板料有限元模型的变形图。 图(a)是整车和假人的有限元模型。它由16,000个壳体单 元、刚体、弹簧、阻尼器以及特殊联结件组成; 图(b)和(c)分别是40m/s和70m/s时汽车和假人的变形图。 下图是一载有假人的整个汽车撞击刚性墙壁动态响应过程的有限元模拟: (3)建立于严格理论基础上的可靠性 因为用于建立有