有限元法发展综述

文章描述:-2022年3月30日发(作者:文德恩)有限元法发展综述 随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规 模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。 这 一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能, 需要对 结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。 例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响

-

有限元法发展综述2022年3月30日发(作者:文德恩)


有限元法发展综述
随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规 模
的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。 这
一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能, 需要对
结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。 例
如分析计算高层建筑和大跨度桥梁在地震时所受到的影响, 看看是否会发生破
坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡 轮机
叶片内的流体动力学参数,以提高其运转效率。这些都可归结为求解物理问 题的控制
偏微分方程式往往是不可能的。近年来在计算机技术和数值分析方法支 持下发展起来
的有限元分析(FEA Finite Element Analysis )方法则为解决这 些复杂的工程分析计算
问题提供了有效的途径。
有限元法是一种高效能、常用的计算方法.有限元法在早期是以变分原理为 基础
发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类 物理场中
(这类场与泛函的极值问题有着紧密的联系)。自从1969年以来,某些 学者在流体力
学中应用加权余数法中的迦辽金法 (Galerki n)或最小二乘法等同
样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理 场
中,而不再要求这类物理场和泛函的极值问题有所联系.
、有限元法的孕育过程及诞生和发展
大约在300年前,牛顿和莱布尼茨发明了积分法,证明了该运算具有整体对 局
部的可加性。虽然,积分运算与有限元技术对定义域的划分是不同的, 前者进
行无限划分而后者进行有限划分,但积分运算为实现有限元技术准备好了一个理 论基
础。
在牛顿之后约一百年,著名数学家高斯提出了加权余值法及线性代数方程组 的
解法。这两项成果的前者被用来将微分方程改写为积分表达式, 后者被用来求
解有限元法所得出的代数方程组。在 18世纪,另一位数学家拉格郎日提出泛函 分
析。泛函分析是将偏微分方程改写为积分表达式的另一途经。
在19世纪末及20世纪初,数学家瑞雷和里兹首先提出可对全定义域运用展 开
函数来表达其上的未知函数。1915年,数学家伽辽金提出了选择展开函数中 形函数
的伽辽金法,该方法被广泛地用于有限元。 1943年,数学家库朗德第一
次提出了可在定义域内分片地使用展开函数来表达其上的未知函数。 这实际上就
是有限兀的做法。
所以,到这时为止,实现有限元技术的第二个理论基础也已确立。
20世纪50年代,飞机设计师们发现无法用传统的力学方法分析飞机的应力、 应
变等问题。波音公司的一个技术小组,首先将连续体的机翼离散为三角形板块 的集合


来进行应力分析,经过一番波折后获得前述的两个离散的成功。 20世纪
50年代,大型电子计算机投入了解算大型代数方程组的工作,这为实现有限元 技术
准备好了物质条件。I960年前后,美国的教授及我国的冯康教 授分别独
立地在论文中提出了 “有限单元”,这样的名词。此后,这样的叫法被 大家接受,
有限元技术从此正式诞生。
1990年10月美国波音公司开始在计算机上对新型客机 B- 777进行“无纸
设计”,仅用了三年半时间,于1994年4月第一架B- 777就试飞成功,这是制 造技
术史上划时代的成就,其中在结构设计和评判中就大量采用有限元分析这一 手段。
在有限元分析的发展初期,由于其基本思想和原理的“简单”和“朴素” ,
以至于许多学术权威都对其学术价值有所鄙视,国际著名刊物 Journal of
Applied Mecha nics 许多年来都拒绝刊登有关于有限元分析的文章。然而现在, 有限
元分析已经成为数值计算的主流,不但国际上存在如ASY箋数种通用有限 元分析软
件,而且涉及到有限元分析的杂志也有几十种之多。
二、有限元法的基本思想
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限 于相
对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉 夫
(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”, 即
有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的) 满足
整个定义域边界条件的允许函数的 Rayleigh Ritz法,有限元法将函数定义 在简单几
何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数), 且不考虑
整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之

有限元方法(FEM的基础是变分原理和加权余量法,其基本求解思想是把 计算
域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为 求解函数
的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所 选用的插值函
数组成的线性表达式 ,借助于变分原理或加权余量法,将微分方
程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。有 限元
方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值 模拟。在
有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元, 在每个单元
内选择基函数,用单元基函数的线形组合来逼近单元中的真解, 整个
计算域上总体的基函数可以看为由每个单元基函数组成的, 则整个计算域内的解
可以看作是由所有单元上的近似解构成。 在河道数值模拟中,常见的有限元计算 方
法是由变分法和加权余量法发展而来的里兹法和伽辽金法、 最小二乘法等。根
据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函 数的
选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的 形状来划
分,有三角形网格、四边形网格和多边形 网格,从插值函数的精度来


划分,又分为线性插值函数和高次插值函数等。不同的组合 同样构成不同的有
限元计算格式。对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的
基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数 的平
方误差最小;在配置法中,先在计算域 内选取个配置点。令近似解在选
定的个配置点上严格满足微分方程,即在配置点上令方程余量为 0。插值函数
一般由不同次幕的多项式组成,但也有采用三角函数或指数函数组成的乘积表 示,但
最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值 多项式本身
在插值点取已知值,称为拉格朗日

Lagrange)多项式插值;另一种不 仅要求插值
多项式本身,还要求它的导数值在插值点取已知值,称为哈密特 (Hermite)多项式
插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对 称和不对称等。常采用
的无因次坐标是一种局部坐标系, 它的定义取决于单元的 几何形状,一维看作长度
比,二维看作面积比,三维看作体积比。在二维有限元 中,三角形单元应用的最早,
近来四边形等参元的应用也越来越广。 对于二维三
角形和四边形电源单元,常采用的插值函数为有 Lagrange插值直角坐标系中的
线性插值函数及二阶或更高阶插值函数、 面积坐标系中的线性插值函数、二阶或
更高阶插值函数等。
对于有限元方法,其解题步骤可归纳为 :
1. 建立积分方程,根据变分原理或方程余量与权函数正交化原理, 建立与微 分
方程初边值问题等价的积分表达式,这是有限元法的出发点。
2. 区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分 为
若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工 作,这
部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关
系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点 序号
和相应的边界值。
3. 确定单元基函数,根据单元中节点数目及对近似解精度的要求, 选择满足 一
定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选 取的,
由于各单元 具有规则的几何形状,在选取基函数时可遵循一定的法则。
4. 单元分析:将各个单元中的求解函数用单元基函数的线性组合表达式进行 逼
近;再将 近似函数代入积分方程,并对单元区域进行积分,可获得含有待定 系数
(即单元中各节点 的参数值

的代数方程组,称为单元有限元方程。
5. 总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按 一
定法则进行累加,形成总体有限元方程。
6. 边界条件的处理:一般边界条件有三种形式,分为本质边界条件

狄里克 雷
边界条件

、自然边界条件

黎曼边界条件

、混合边界条件

柯西边界条件


对于自然边界条件,一般在积分表达式中可自动得到满足。 对于本质边界条件和 混
合边界条件,需按一定法则对总体有限元方程进行修正满足。
7. 解有限元方程:根据边界条件修正的总体有限元方程组, 是含所有待定未 知
量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。
三、有限元的应用及其发展趋势
有限元的应用范围也是相当的广的。它涉及到工程结构、传热、流体运动、 电


磁等连续介质的力学分析中,并在气象、地球物理、医学等领域得到应用和发 展。电
子计算机的出现和发展是有限元法在许多实际问题中的应用变为现实, 并
具有广阔的前景。
国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具 有强
大功能的有限元分析程序。其中最为著名的是由美国国家宇航局
(
ASA)在
1965年委托美国计算科学公司和贝尔航空系统公司开发的 ASTRA有限元分析
系统。该系统发展至今已有几十个版本

是目前世界上规模最大、功能最强的有 限元
分析系统。从那时到现在

世界各地的研究机构和大学也发展了一批规模较 小但使用
灵活、价格较低的专用或通用有限元分析软件
,
主要有德国的ASKA英 国的 PAFEC
法国的 SYSTUS 美国的 ABQU
S
ADIA ASYS BERSAFEBOSOR COSMOSELAS
MARC口 STARDY等公司的产品。当今国际上 FEA方法和软件发 展呈现出以下
一些趋势特征

1. 从单纯的结构力学计算发展到求解许多物理场问题
有限元分析方法最早是从结构化矩阵分析发展而来

逐步推广到板、壳和实 体
等连续体固体力学分析

实践证明这是一种非常有效的数值分析方法。而且从 理论上
也已经证明

只要用于离散求解对象的单元足够小

所得的解就可足够逼 近于精确
值。所以近年来有限元方法已发展到流体力学、 温度场、电传导、磁场、 渗流和声
场等问题的求解计算

最近又发展到求解几个交叉学科的问题。例如当 气流流过一个
很高的铁塔时就会使铁塔产生变形

而塔的变形又反过来影响到气 流的流动,,这就需
要用固体力学和流体动力学的有限元分析结果交叉迭代求 解

即所谓流固耦合的问
题。
2. 由求解线性工程问题进展到分析非线性问题
随着科学技术的发展

线性理论已经远远不能满足设计的要求。例如建筑行 业
中的高层建筑和大跨度悬索桥的出现

就要求考虑结构的大位移和大应变等几 何非线
性问题

航天和动力工程的高温部件存在热变形和热应力

也要考虑材料
的非线性问题

诸如塑料、橡胶和复合材料等各种新材料的出现

仅靠线性计算理
论就不足以解决遇到的问题

只有采用非线性有限元算法才能解决。众所周知

非线性的数值计算是很复杂的

它涉及到很多专门的数学问题和运算技巧

很难 为
一般工程技术人员所掌握。为此近年来国外一些公司花费了大量的人力和投资 开发诸
如MARCABQU和 ADIA等专长于求解非线性问题的有限元分析软件
,
并 广泛应用
于工程实践。这些软件的共同特点是具有高效的非线性求解器以及丰富 和实用的非线
性材料库。
3. 增强可视化的前置建模和后置数据处理功能
早期有限元分析软件的研究重点在于推导新的高效率求解方法和高精度的 单
元。随着数值分析方法的逐步完善

尤其是计算机运算速度的飞速发展

整个计 算
系统用于求解运算的时间越来越少

而数据准备和运算结果的表现问题却日益 突出。
在现在的工程工作站上

求解一个包含10万个方程的有限元模型只需要用 几十分
钟。但是如果用手工方式来建立这个模型

然后再处理大量的计算结果则 需用几周的
时间。可以毫不夸张地说

工程师在分析计算一个工程问题时有 80% 以上的精力都
花在数据准备和结果分析上。因此目前几乎所有的商业化有限元程 序系统都有功能很
强的前置建模和后置数据处理模块。在强调 可视化的今天


很多程序都建立了对用户非常友好的 GUI(Graphics User In terface), 使用户能
以可视图形方式直观快速地进行网格自动划分

生成有限元分析所需数据

并按
要求将大量的计算结果整理成变形图、等值分布云图

便于极值搜索和所需数据
的列表输出。
4. 与CAD软件的无缝集成
当今有限元分析系统的另一个特点是与通用 CAD软件的集成使用

即在用
CAD软件完成部件和零件的造型设计后

自动生成有限元网格并进行计算

如果 分
析的结果不符合设计要求则重新进行造型和计算

直到满意为止

从而极大地 提高
了设计水平和效率。今天
,
工程师可以在集成的CAD和FEA软件环境中快捷 地解决
一个在以前无法应付的复杂工程分析问题。 所以当今所有的商业化有限元
系统商都开发了和著名的CAD软件
(
例如ProEGIEER Unigraphics、
SolidEdge、SolidWorks、IDEAS Bentley 和 AutoCAD等)的接口。
5. 在Win tel平台上的发展
早期的有限元分析软件基本上都是在大中型计算机
(
主要是Mainframe)
上开发和运行的

后来又发展到以工程工作站(EWS,Engineering WorkStation ) 为平


它们的共同特点都是采用 UIX操作系统。PC机的出现使计算机的应用 发生了
根本性的变化

工程师渴望在办公桌上完成复杂工程分析的梦想成为现 实。但是早期
的PC机采用16位CPU和 DOS操作系统

内存中的公共数据块受到 限制
,
因此当
时计算模型的规模不能超过 1万阶方程。Microsoft Windows 操作 系统和32位的
Intel Pentium 处理器的推出为将PC机用于有限元分析提供了
必需的软件和硬件支撑平台。因此当前国际上著名的有限元程序研究和发展机构 都纷
纷将他们的软件移植到 Win tel平台上。下表列出了用 ADIA V7.3版在PC 机的
WindowsT环境和SGI工作站上同时计算4个工程实例所需要的求解时间。 从中可
以看出最新高档PC机的求解能力已和中低挡的 EWSf相上下。为了将在 大中型计算
机和 EWS上开发的有限元程序移植到 PC机上
,
常常需要采用
Hummingbird公司的一个仿真软件Exceed。这样做的结果比较麻烦
,
而且不能充 分
利用PC机的软硬件资源。所以最近有些公司
,
例如IDEAS ADIA和 R&D开始 在
Windows平台上开发有限元程序

称作版本

同时还有在PC 机上
的Linux操作系统环境中开发的有限元程序包。
在大力推广CAD技术的今天,从自行车到航天飞机,所有的设计制造都离不 开
有限元分析计算,有限元法在工程设计和分析中将得到越来越广泛的重视。 目
前以分析、优化和仿真为特征的 CAE(Computer aided Engineering CAE)技术在 世
界范围内蓬勃发展。它通过先进的 CAE技术快速有效地分析产品的各种特性、 揭示
结构各类参数变化对产品性能的响, 进行设计方案的修改和调整,使产品达 到性能
和质量上的最优,原材料消耗最低。因此,基于计算机的分析、优化和仿 真的CAE
技术的研究和应用,是高质量、高水平、低成本产品设计与开发的保证。 有限元法也
必将在科技发展史上大放异彩。


有限元法发展综述
随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规 模
的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。 这
一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能, 需要对
结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。 例
如分析计算高层建筑和大跨度桥梁在地震时所受到的影响, 看看是否会发生破
坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡 轮机
叶片内的流体动力学参数,以提高其运转效率。这些都可归结为求解物理问 题的控制
偏微分方程式往往是不可能的。近年来在计算机技术和数值分析方法支 持下发展起来
的有限元分析(FEA Finite Element Analysis )方法则为解决这 些复杂的工程分析计算
问题提供了有效的途径。
有限元法是一种高效能、常用的计算方法.有限元法在早期是以变分原理为 基础
发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类 物理场中
(这类场与泛函的极值问题有着紧密的联系)。自从1969年以来,某些 学者在流体力
学中应用加权余数法中的迦辽金法 (Galerki n)或最小二乘法等同
样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理 场
中,而不再要求这类物理场和泛函的极值问题有所联系.
、有限元法的孕育过程及诞生和发展
大约在300年前,牛顿和莱布尼茨发明了积分法,证明了该运算具有整体对 局
部的可加性。虽然,积分运算与有限元技术对定义域的划分是不同的, 前者进
行无限划分而后者进行有限划分,但积分运算为实现有限元技术准备好了一个理 论基
础。
在牛顿之后约一百年,著名数学家高斯提出了加权余值法及线性代数方程组 的
解法。这两项成果的前者被用来将微分方程改写为积分表达式, 后者被用来求
解有限元法所得出的代数方程组。在 18世纪,另一位数学家拉格郎日提出泛函 分
析。泛函分析是将偏微分方程改写为积分表达式的另一途经。
在19世纪末及20世纪初,数学家瑞雷和里兹首先提出可对全定义域运用展 开
函数来表达其上的未知函数。1915年,数学家伽辽金提出了选择展开函数中 形函数
的伽辽金法,该方法被广泛地用于有限元。 1943年,数学家库朗德第一
次提出了可在定义域内分片地使用展开函数来表达其上的未知函数。 这实际上就
是有限兀的做法。
所以,到这时为止,实现有限元技术的第二个理论基础也已确立。
20世纪50年代,飞机设计师们发现无法用传统的力学方法分析飞机的应力、 应
变等问题。波音公司的一个技术小组,首先将连续体的机翼离散为三角形板块 的集合


来进行应力分析,经过一番波折后获得前述的两个离散的成功。 20世纪
50年代,大型电子计算机投入了解算大型代数方程组的工作,这为实现有限元 技术
准备好了物质条件。I960年前后,美国的教授及我国的冯康教 授分别独
立地在论文中提出了 “有限单元”,这样的名词。此后,这样的叫法被 大家接受,
有限元技术从此正式诞生。
1990年10月美国波音公司开始在计算机上对新型客机 B- 777进行“无纸
设计”,仅用了三年半时间,于1994年4月第一架B- 777就试飞成功,这是制 造技
术史上划时代的成就,其中在结构设计和评判中就大量采用有限元分析这一 手段。
在有限元分析的发展初期,由于其基本思想和原理的“简单”和“朴素” ,
以至于许多学术权威都对其学术价值有所鄙视,国际著名刊物 Journal of
Applied Mecha nics 许多年来都拒绝刊登有关于有限元分析的文章。然而现在, 有限
元分析已经成为数值计算的主流,不但国际上存在如ASY箋数种通用有限 元分析软
件,而且涉及到有限元分析的杂志也有几十种之多。
二、有限元法的基本思想
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限 于相
对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉 夫
(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”, 即
有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的) 满足
整个定义域边界条件的允许函数的 Rayleigh Ritz法,有限元法将函数定义 在简单几
何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数), 且不考虑
整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之

有限元方法(FEM的基础是变分原理和加权余量法,其基本求解思想是把 计算
域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为 求解函数
的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所 选用的插值函
数组成的线性表达式 ,借助于变分原理或加权余量法,将微分方
程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。有 限元
方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值 模拟。在
有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元, 在每个单元
内选择基函数,用单元基函数的线形组合来逼近单元中的真解, 整个
计算域上总体的基函数可以看为由每个单元基函数组成的, 则整个计算域内的解
可以看作是由所有单元上的近似解构成。 在河道数值模拟中,常见的有限元计算 方
法是由变分法和加权余量法发展而来的里兹法和伽辽金法、 最小二乘法等。根
据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函 数的
选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的 形状来划
分,有三角形网格、四边形网格和多边形 网格,从插值函数的精度来


划分,又分为线性插值函数和高次插值函数等。不同的组合 同样构成不同的有
限元计算格式。对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的
基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数 的平
方误差最小;在配置法中,先在计算域 内选取个配置点。令近似解在选
定的个配置点上严格满足微分方程,即在配置点上令方程余量为 0。插值函数
一般由不同次幕的多项式组成,但也有采用三角函数或指数函数组成的乘积表 示,但
最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值 多项式本身
在插值点取已知值,称为拉格朗日

Lagrange)多项式插值;另一种不 仅要求插值
多项式本身,还要求它的导数值在插值点取已知值,称为哈密特 (Hermite)多项式
插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对 称和不对称等。常采用
的无因次坐标是一种局部坐标系, 它的定义取决于单元的 几何形状,一维看作长度
比,二维看作面积比,三维看作体积比。在二维有限元 中,三角形单元应用的最早,
近来四边形等参元的应用也越来越广。 对于二维三
角形和四边形电源单元,常采用的插值函数为有 Lagrange插值直角坐标系中的
线性插值函数及二阶或更高阶插值函数、 面积坐标系中的线性插值函数、二阶或
更高阶插值函数等。
对于有限元方法,其解题步骤可归纳为 :
1. 建立积分方程,根据变分原理或方程余量与权函数正交化原理, 建立与微 分
方程初边值问题等价的积分表达式,这是有限元法的出发点。
2. 区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分 为
若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工 作,这
部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关
系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点 序号
和相应的边界值。
3. 确定单元基函数,根据单元中节点数目及对近似解精度的要求, 选择满足 一
定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选 取的,
由于各单元 具有规则的几何形状,在选取基函数时可遵循一定的法则。
4. 单元分析:将各个单元中的求解函数用单元基函数的线性组合表达式进行 逼
近;再将 近似函数代入积分方程,并对单元区域进行积分,可获得含有待定 系数
(即单元中各节点 的参数值

的代数方程组,称为单元有限元方程。
5. 总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按 一
定法则进行累加,形成总体有限元方程。
6. 边界条件的处理:一般边界条件有三种形式,分为本质边界条件

狄里克 雷
边界条件

、自然边界条件

黎曼边界条件

、混合边界条件

柯西边界条件


对于自然边界条件,一般在积分表达式中可自动得到满足。 对于本质边界条件和 混
合边界条件,需按一定法则对总体有限元方程进行修正满足。
7. 解有限元方程:根据边界条件修正的总体有限元方程组, 是含所有待定未 知
量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。
三、有限元的应用及其发展趋势
有限元的应用范围也是相当的广的。它涉及到工程结构、传热、流体运动、 电


磁等连续介质的力学分析中,并在气象、地球物理、医学等领域得到应用和发 展。电
子计算机的出现和发展是有限元法在许多实际问题中的应用变为现实, 并
具有广阔的前景。
国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具 有强
大功能的有限元分析程序。其中最为著名的是由美国国家宇航局
(
ASA)在
1965年委托美国计算科学公司和贝尔航空系统公司开发的 ASTRA有限元分析
系统。该系统发展至今已有几十个版本

是目前世界上规模最大、功能最强的有 限元
分析系统。从那时到现在

世界各地的研究机构和大学也发展了一批规模较 小但使用
灵活、价格较低的专用或通用有限元分析软件
,
主要有德国的ASKA英 国的 PAFEC
法国的 SYSTUS 美国的 ABQU
S
ADIA ASYS BERSAFEBOSOR COSMOSELAS
MARC口 STARDY等公司的产品。当今国际上 FEA方法和软件发 展呈现出以下
一些趋势特征

1. 从单纯的结构力学计算发展到求解许多物理场问题
有限元分析方法最早是从结构化矩阵分析发展而来

逐步推广到板、壳和实 体
等连续体固体力学分析

实践证明这是一种非常有效的数值分析方法。而且从 理论上
也已经证明

只要用于离散求解对象的单元足够小

所得的解就可足够逼 近于精确
值。所以近年来有限元方法已发展到流体力学、 温度场、电传导、磁场、 渗流和声
场等问题的求解计算

最近又发展到求解几个交叉学科的问题。例如当 气流流过一个
很高的铁塔时就会使铁塔产生变形

而塔的变形又反过来影响到气 流的流动,,这就需
要用固体力学和流体动力学的有限元分析结果交叉迭代求 解

即所谓流固耦合的问
题。
2. 由求解线性工程问题进展到分析非线性问题
随着科学技术的发展

线性理论已经远远不能满足设计的要求。例如建筑行 业
中的高层建筑和大跨度悬索桥的出现

就要求考虑结构的大位移和大应变等几 何非线
性问题

航天和动力工程的高温部件存在热变形和热应力

也要考虑材料
的非线性问题

诸如塑料、橡胶和复合材料等各种新材料的出现

仅靠线性计算理
论就不足以解决遇到的问题

只有采用非线性有限元算法才能解决。众所周知

非线性的数值计算是很复杂的

它涉及到很多专门的数学问题和运算技巧

很难 为
一般工程技术人员所掌握。为此近年来国外一些公司花费了大量的人力和投资 开发诸
如MARCABQU和 ADIA等专长于求解非线性问题的有限元分析软件
,
并 广泛应用
于工程实践。这些软件的共同特点是具有高效的非线性求解器以及丰富 和实用的非线
性材料库。
3. 增强可视化的前置建模和后置数据处理功能
早期有限元分析软件的研究重点在于推导新的高效率求解方法和高精度的 单
元。随着数值分析方法的逐步完善

尤其是计算机运算速度的飞速发展

整个计 算
系统用于求解运算的时间越来越少

而数据准备和运算结果的表现问题却日益 突出。
在现在的工程工作站上

求解一个包含10万个方程的有限元模型只需要用 几十分
钟。但是如果用手工方式来建立这个模型

然后再处理大量的计算结果则 需用几周的
时间。可以毫不夸张地说

工程师在分析计算一个工程问题时有 80% 以上的精力都
花在数据准备和结果分析上。因此目前几乎所有的商业化有限元程 序系统都有功能很
强的前置建模和后置数据处理模块。在强调 可视化的今天


很多程序都建立了对用户非常友好的 GUI(Graphics User In terface), 使用户能
以可视图形方式直观快速地进行网格自动划分

生成有限元分析所需数据

并按
要求将大量的计算结果整理成变形图、等值分布云图

便于极值搜索和所需数据
的列表输出。
4. 与CAD软件的无缝集成
当今有限元分析系统的另一个特点是与通用 CAD软件的集成使用

即在用
CAD软件完成部件和零件的造型设计后

自动生成有限元网格并进行计算

如果 分
析的结果不符合设计要求则重新进行造型和计算

直到满意为止

从而极大地 提高
了设计水平和效率。今天
,
工程师可以在集成的CAD和FEA软件环境中快捷 地解决
一个在以前无法应付的复杂工程分析问题。 所以当今所有的商业化有限元
系统商都开发了和著名的CAD软件
(
例如ProEGIEER Unigraphics、
SolidEdge、SolidWorks、IDEAS Bentley 和 AutoCAD等)的接口。
5. 在Win tel平台上的发展
早期的有限元分析软件基本上都是在大中型计算机
(
主要是Mainframe)
上开发和运行的

后来又发展到以工程工作站(EWS,Engineering WorkStation ) 为平


它们的共同特点都是采用 UIX操作系统。PC机的出现使计算机的应用 发生了
根本性的变化

工程师渴望在办公桌上完成复杂工程分析的梦想成为现 实。但是早期
的PC机采用16位CPU和 DOS操作系统

内存中的公共数据块受到 限制
,
因此当
时计算模型的规模不能超过 1万阶方程。Microsoft Windows 操作 系统和32位的
Intel Pentium 处理器的推出为将PC机用于有限元分析提供了
必需的软件和硬件支撑平台。因此当前国际上著名的有限元程序研究和发展机构 都纷
纷将他们的软件移植到 Win tel平台上。下表列出了用 ADIA V7.3版在PC 机的
WindowsT环境和SGI工作站上同时计算4个工程实例所需要的求解时间。 从中可
以看出最新高档PC机的求解能力已和中低挡的 EWSf相上下。为了将在 大中型计算
机和 EWS上开发的有限元程序移植到 PC机上
,
常常需要采用
Hummingbird公司的一个仿真软件Exceed。这样做的结果比较麻烦
,
而且不能充 分
利用PC机的软硬件资源。所以最近有些公司
,
例如IDEAS ADIA和 R&D开始 在
Windows平台上开发有限元程序

称作版本

同时还有在PC 机上
的Linux操作系统环境中开发的有限元程序包。
在大力推广CAD技术的今天,从自行车到航天飞机,所有的设计制造都离不 开
有限元分析计算,有限元法在工程设计和分析中将得到越来越广泛的重视。 目
前以分析、优化和仿真为特征的 CAE(Computer aided Engineering CAE)技术在 世
界范围内蓬勃发展。它通过先进的 CAE技术快速有效地分析产品的各种特性、 揭示
结构各类参数变化对产品性能的响, 进行设计方案的修改和调整,使产品达 到性能
和质量上的最优,原材料消耗最低。因此,基于计算机的分析、优化和仿 真的CAE
技术的研究和应用,是高质量、高水平、低成本产品设计与开发的保证。 有限元法也
必将在科技发展史上大放异彩。

-

有限元法发展综述

发布时间:2022-03-30 03:32:22
文章版权声明:除非注明,否则均为IT技术网-学习WEB前端开发等IT技术的网络平台原创文章,转载或复制请以超链接形式并注明出处。

发表评论

评论列表 (有 5 条评论,641人围观)
宫颈癌症状V铁粉24 minutes ago Google Chrome 93.0.4577.82 Windows 10 x64
橡胶和复合材料等各种新材料的出现
两地生活V铁粉1 minute ago Google Chrome 93.0.4577.82 Windows 10 x64
2. 区域单元剖分
新汶矿务局V铁粉3 minutes ago Google Chrome 93.0.4577.82 Windows 10 x64
而且不能充 分利用PC机的软硬件资源
苦心人天不负V铁粉22 minutes ago Google Chrome 93.0.4577.82 Windows 10 x64
I960年前后