• 1.摘要
  • 2.基本信息
  • 3.基础定义
  • 4.应用举例
  • 5.演绎过程
  • 6.派生

有限元法

数值计算方法

有限元法(finite element method)是一种高效能、常用的数值计算方法。科学计算领域,常常需要求解各类微分方程,而许多微分方程的解析解一般很难得到,使用有限元法将微分方程离散化后,可以编制程序,使用计算机辅助求解。有限元法在早期是以变分原理为基础发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中(这类场与泛函的极值问题有着紧密的联系)。自从1969年以来,某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理场中,而不再要求这类物理场和泛函的极值问题有所联系。基本思想:由解给定的泊松方程化为求解泛函的极值问题。

基本信息

  • 中文名

    有限元法

  • 外文名

    finite element method

  • 所属领域

    数学

  • 概述

    一种高效能常用的计算方法

  • 别名

    有限元分析方法

  • 应用范围

    弹性力学问题

基础定义

将连续的求解域离散为一组单元的组合体,用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数,近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。

有限元分析是使用有限元方法来分析静态或动态的物理物体或物理系统。在这种方法中一个物体或系统被分解为由多个相互联结的、简单、独立的点组成的几何模型。在这种方法中这些独立的点的数量是有限的,因此被称为有限元。由实际的物理模型中推导出来得平衡方程式被使用到每个点上,由此产生了一个方程组。这个方程组可以用线性代数的方法来求解。有限元分析的精确度无法无限提高。元的数目到达一定高度后解的精确度不再提高,只有计算时间不断提高。

有限元分析法(FEA)已应用得非常广泛,现已成为年创收达数十亿美元的相关产业的基础。即使是很复杂的应力问题的数值解,用有限元分析的常规方法就能得到。此方法是如此的重要,以至于即便像这些只对材料力学作入门性论述的模块,也应该略述其主要特点。 不管有限元法是如何的卓有成效,当你应用此法及类似的方法时,计算机解的缺点必须牢记在心头:这些解不一定能揭示诸如材料性能、几何特征等重要的变量是如何影响应力的。一旦输入数据有误,结果就会大相径庭,而分析者却难以觉察。所以理论建模最重要的作用可能是使设计者的直觉变得敏锐。有限元程序的用户应该为此目标部署设计策略,以尽可能多的封闭解和实验分析作为计算机仿真的补充。 与现代微机上许多字处理和电子制表软件包相比,有限元的程序不那么复杂。然而,这些程序的复杂程度依然使大部分用户无法有效地编写自己所需的程序。可以买到一些预先编好的商用程序1,其价格范围宽,从微机到超级计算机都可兼容。但有特定需求的用户也不必对程序的开发望而生畏,你会发现,从诸如齐凯维奇(Zienkiewicz2)等的教材中提供的程序资源可作为有用的起点。大部分有限元软件是用Fortran语言编写的,但诸如felt等某些更新的程序用的是C语言或其它更时新的程序语言。

在实践中,有限元分析法通常由三个主要步骤组成: 1、预处理:用户需建立物体待分析部分的模型,在此模型中,该部分的几何形状被分割成若干个离散的子区域——或称为“单元”。各单元在一些称为“结点”的离散点上相互连接。这些结点中有的有固定的位移,而其余的有给定的载荷。准备这样的模型可能极其耗费时间,所以商用程序之间的相互竞争就在于:如何用最友好的图形化界面的“预处理模块”,来帮助用户完成这项繁琐乏味的工作。有些预处理模块作为计算机化的画图和设计过程的组成部分,可在先前存在的CAD文件中覆盖网格,因而可以方便地完成有限元分析。 2、分析:把预处理模块准备好的数据输入到有限元程序中,从而构成并求解用线性或非线性代数方程表示的系统

u和f分别为各结点的位移和作用的外力。矩阵K的形式取决于求解问题的类型。 3、分析的早期,用户需仔细地研读程序运算后产生的大量数字,即 型,本模块将概述桁架与线弹性体应力分析的方法。商用程序可能带有非常大的单元库,不同类型的单元适用于范围广泛的各类问题。有限元法的主要优点之一就是:许多不同类型的问题都可用相同的程序来处理,区别仅在于从单元库中指定适合于不同问题的单元类型。

应用举例

1.弹性力学分析问题

2.平衡问题

3.固体力学

4.工程力学

演绎过程

步骤1:剖分:

将待解区域进行分割,离散成有限个元素的集合。元素(单元)的形状原则上是任意的。二维问题一般采用三角形单元或矩形单元,三维空间可采用四面体或多面体等。每个单元的顶点称为节点(或结点)。

步骤2:单元分析:

进行分片插值,即将分割单元中任意点的未知函数用该分割单元中形状函数及离散网格点上的函数值展开,即建立一个线性插值函数。

步骤3:求解近似变分方程

用有限个单元将连续体离散化,通过对有限个单元作分片插值求解各种力学、物理问题的一种数值方法。有限元法把连续体离散成有限个单元:杆系结构的单元是每一个杆件;连续体的单元是各种形状(如三角形、四边形、六面体等)的单元体。每个单元的场函数是只包含有限个待定节点参量的简单场函数,这些单元场函数的集合就能近似代表整个连续体的场函数。根据能量方程或加权残量方程可建立有限个待定参量的代数方程组,求解此离散方程组就得到有限元法的数值解。有限元法已被用于求解线性和非线性问题,并建立了各种有限元模型,如协调、不协调、混合、杂交、拟协调元等。有限元法十分有效、通用性强、应用广泛,已有许多大型或专用程序系统供工程设计使用。结合计算机辅助设计技术,有限元法也被用于计算机辅助制造中。

有限单元法最早可上溯到20世纪40年代。Courant第一次应用定义在三角区域上的分片连续函数和最小位能原理来求解St.Venant扭转问题。现代有限单元法的第一个成功的尝试是在 1956年,Turner、Clough等人在分析飞机结构时,将钢架位移法推广应用于弹性力学平面问题,给出了用三角形单元求得平面应力问题的正确答案。1960年,Clough进一步处理了平面弹性问题,并第一次提出了"有限单元法",使人们认识到它的功效。