有限单元法简介论述.pptVIP

(2)有限元法的理论基础和离散格式 研究工作的进展包括: 将多场变量的变分原理用于有限元分析，发展了混合型(单元内包括多个场变量)、 杂交型(某些场变量仅在单元交界面定义)的有限元表达格式，并研究了各自的收 敛性条件； 将与微分方程等效的积分形式——加权余量法，用于建立有限元的表达格式，从 而将有限元的应用扩展到不存在泛函或泛函尚未建立的物理问题； 有限元解的后验误差估计和应力磨平方法的研究进展，不仅改进了有限元解的精 度，更重要的是为发展满足规定精度的要求，以细分单元网格或提高插值函数阶 次为手段的自适应分析方法提供了基础。 (3)有限元方程的解法 现在用于大型复杂工程问题的有限元分析，自由度达几十万个甚至上百万个已是经常的情况，这是与计算机软、硬件发展相配合的大型方程组解法的研究进展密不可分。有限元求解的问题从性质上可以归结为三类，即: ①独立于时间的平衡问题(或稳态问题) 最后归结为求解系数矩阵元素在对角线附近稀疏分布的线性代数方程组。此类问题至今主要是采用直接解法，先后发展了循序消去法、三角分解法、波前法等。近年来，为了适应求解大型、特大型方程时减少计算机存储和提高计算速度的需要，迭代解法特别是预条件共轭梯度法受到更多的重视，并已成功地应用。 * 有限单元法简介 是在当今工程分析中获得最广泛应用的数值计算方法。 由于它的通用性和有效性，受到工程技术界的高度重视。 伴随着计算机科学和技术的快速发展，现已成为： 计算机辅助设计(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)分别是40ms和70ms