有限元法发展综述及其特点(共9页).doc

《有限元法发展综述及其特点(共9页).doc》由会员分享，可在线阅读，更多相关《有限元法发展综述及其特点(共9页).doc（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上数 值 分 析 结 课 论 文有限元的发展历程及其特点论文题目：有限元的发展历程及其特点学 院: 专 业：学 号：姓 名：有限元法发展综述及其特点摘要：1965年“有限元”这个名词第一次出现，到今天有限元在工程上得到广泛应用，经历了三十多年的发展历史，理论和算法都已经日趋完善。有限元的核心思想是结构的离散化，就是将实际结构假想地离散为有限数目的规则单元组合体，实际结构的物理性能可以通过对离散体进行分析，得出满足工程精度的近似结果来替代对实际结构的分析，这样可以解决很多实际工程需要解决而理论分析又无法解决的复杂问题。关键词：有限元，积分法，加权余值法，边值，伽辽金(Ga

2、lerkin)法。引言有限元法是一种高效能、常用的计算方法有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系 有限元法的孕育过程及诞生和发展大约在300年前，牛顿和莱布尼茨发明了积分法，证明了该运算具有整体对局部的可加性。虽然，积分运算与有限元技术对定义域的划分是不同的，前者进行

3、无限划分而后者进行有限划分，但积分运算为实现有限元技术准备好了一个理论基础。在牛顿之后约一百年，著名数学家高斯提出了加权余值法及线性代数方程组的解法。这两项成果的前者被用来将微分方程改写为积分表达式，后者被用来求解有限元法所得出的代数方程组。在18世纪，另一位数学家拉格郎日提出泛函分析。泛函分析是将偏微分方程改写为积分表达式的另一途经。在19世纪末及20世纪初，数学家瑞雷和里兹首先提出可对全定义域运用展开函数来表达其上的未知函数。1915年，数学家伽辽金提出了选择展开函数中形函数的伽辽金法，该方法被广泛地用于有限元。1943年，数学家库朗德第一次提出了可在定义域内分片地使用展开函数来表达其上的

4、未知函数。这实际上就是有限元的做法。所以，到这时为止，实现有限元技术的第二个理论基础也已确立。20世纪50年代，飞机设计师们发现无法用传统的力学方法分析飞机的应力、应变等问题。波音公司的一个技术小组，首先将连续体的机翼离散为三角形板块的集合来进行应力分析，经过一番波折后获得前述的两个离散的成功。20世纪50年代，大型电子计算机投入了解算大型代数方程组的工作，这为实现有限元技术准备好了物质条件。1960年前后，美国的R.W.Clough教授及我国的冯康教授分别独立地在论文中提出了“有限单元”，这样的名词。此后，这样的叫法被大家接受，有限元技术从此正式诞生。1990年10月美国波音公司开始在计算机

5、上对新型客机B777进行“无纸设计”，仅用了三年半时间，于1994年4月第一架B777就试飞成功，这是制造技术史上划时代的成就，其中在结构设计和评判中就大量采用有限元分析这一手段。 在有限元分析的发展初期，由于其基本思想和原理的“简单”和“朴素”，以至于许多学术权威都对其学术价值有所鄙视，国际著名刊物Journal of Applied Mechanics 许多年来都拒绝刊登有关于有限元分析的文章。然而现在，有限元分析已经成为数值计算的主流，不但国际上存在如ANSYS等数种通用有限元分析软件，而且涉及到有限元分析的杂志也有几十种之多。和特点有限元方法与其他求解边值问题近似方法的根本区别在于它的

6、近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫（Clough）教授形象地将其描绘为：“有限元法=Rayleigh Ritz法分片函数”，即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解（往往是困难的）满足整个定义域边界条件的允许函数的Rayleigh Ritz法，有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其他近似方法的原因之一。有限元方法（FEM）的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元

7、内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式 ，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分

8、法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形 网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合 同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域 内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置

9、点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三

10、角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。对于有限元方法，其解题步骤可归纳为 ： 1.建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。   2.区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时还

11、需要列出自然边界和本质边界的节点序号和相应的边界值。 3.确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的，由于各单元 具有规则的几何形状，在选取基函数时可遵循一定的法则。 4.单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将 近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点 的参数值)的代数方程组，称为单元有限元方程。5.总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。 6.边界条件的处理：一般

12、边界条件有三种形式，分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件，一般在积分表达式中可自动得到满足。对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。  7.解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。三、有限元的应用及其发展趋势有限元的应用范围也是相当的广的。它涉及到工程结构、传热、流体运动、电磁等连续介质的力学分析中，并在气象、地球物理、医学等领域得到应用和发展。电子计算机的出现和发展是有限元法在许多

13、实际问题中的应用变为现实，并具有广阔的前景。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局(NASA)在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本,是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在,世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件,主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、E

14、LAS、MARC和STARDYNE等公司的产品。当今国际上FEA方法和软件发展呈现出以下一些趋势特征: 1.从单纯的结构力学计算发展到求解许多物理场问题 有限元分析方法最早是从结构化矩阵分析发展而来,逐步推广到板、壳和实体等连续体固体力学分析,实践证明这是一种非常有效的数值分析方法。而且从理论上也已经证明,只要用于离散求解对象的单元足够小,所得的解就可足够逼近于精确值。所以近年来有限元方法已发展到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算,最近又发展到求解几个交叉学科的问题。例如当气流流过一个很高的铁塔时就会使铁塔产生变形,而塔的变形又反过来影响到气流的流动这就需要用固体力学和

15、流体动力学的有限元分析结果交叉迭代求解,即所谓"流固耦合"的问题。 2. 由求解线性工程问题进展到分析非线性问题 随着科学技术的发展,线性理论已经远远不能满足设计的要求。例如建筑行业中的高层建筑和大跨度悬索桥的出现,就要求考虑结构的大位移和大应变等几何非线性问题;航天和动力工程的高温部件存在热变形和热应力,也要考虑材料的非线性问题;诸如塑料、橡胶和复合材料等各种新材料的出现,仅靠线性计算理论就不足以解决遇到的问题,只有采用非线性有限元算法才能解决。众所周知,非线性的数值计算是很复杂的,它涉及到很多专门的数学问题和运算技巧,很难为一般工程技术人员所掌握。为此近年来国外一些公司

16、花费了大量的人力和投资开发诸如MARC、ABQUS和ADINA 等专长于求解非线性问题的有限元分析软件,并广泛应用于工程实践。这些软件的共同特点是具有高效的非线性求解器以及丰富和实用的非线性材料库。3. 增强可视化的前置建模和后置数据处理功能 早期有限元分析软件的研究重点在于推导新的高效率求解方法和高精度的单元。随着数值分析方法的逐步完善,尤其是计算机运算速度的飞速发展,整个计算系统用于求解运算的时间越来越少,而数据准备和运算结果的表现问题却日益突出。在现在的工程工作站上,求解一个包含10万个方程的有限元模型只需要用几十分钟。但是如果用手工方式来建立这个模型,然后再处理大量的计算结果则需用几周

17、的时间。可以毫不夸张地说,工程师在分析计算一个工程问题时有80%以上的精力都花在数据准备和结果分析上。因此目前几乎所有的商业化有限元程序系统都有功能很强的前置建模和后置数据处理模块。在强调"可视化"的今天,很多程序都建立了对用户非常友好的GUI(Graphics User Interface),使用户能以可视图形方式直观快速地进行网格自动划分,生成有限元分析所需数据,并按要求将大量的计算结果整理成变形图、等值分布云图,便于极值搜索和所需数据的列表输出。4. 与CAD软件的无缝集成当今有限元分析系统的另一个特点是与通用CAD软件的集成使用,即在用CAD软件完成部件和零件的造型

18、设计后,自动生成有限元网格并进行计算,如果分析的结果不符合设计要求则重新进行造型和计算,直到满意为止,从而极大地提高了设计水平和效率。今天,工程师可以在集成的CAD和FEA软件环境中快捷地解决一个在以前无法应付的复杂工程分析问题。所以当今所有的商业化有限元系统商都开发了和著名的CAD软件(例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley和AutoCAD等)的接口。 5. 在Wintel平台上的发展早期的有限元分析软件基本上都是在大中型计算机 (主要是Mainframe) 上开发和运行的,后来又发展到以工程工作站(EWS,

19、Engineering WorkStation )为平台,它们的共同特点都是采用UNIX操作系统。PC机的出现使计算机的应用发生了根本性的变化,工程师渴望在办公桌上完成复杂工程分析的梦想成为现实。但是早期的PC机采用16位CPU和DOS操作系统,内存中的公共数据块受到限制,因此当时计算模型的规模不能超过1万阶方程。Microsoft Windows操作系统和32位的Intel Pentium 处理器的推出为将PC机用于有限元分析提供了必需的软件和硬件支撑平台。因此当前国际上著名的有限元程序研究和发展机构都纷纷将他们的软件移植到Wintel平台上。下表列出了用ADINA V7.3 版在PC机的W

20、indows NT环境和SGI工作站上同时计算4个工程实例所需要的求解时间。从中可以看出最新高档PC机的求解能力已和中低挡的EWS不相上下。为了将在大中型计算机和EWS上开发的有限元程序移植到PC机上,常常需要采用Hummingbird公司的一个仿真软件Exceed。这样做的结果比较麻烦,而且不能充分利用PC机的软硬件资源。所以最近有些公司,例如IDEAS、ADINA和R&D开始在Windows平台上开发有限元程序,称作"Native Windows"版本,同时还有在PC机上的Linux操作系统环境中开发的有限元程序包。 在大力推广CAD技术的今天，从自行车到航天飞

21、机，所有的设计制造都离不开有限元分析计算，有限元法在工程设计和分析中将得到越来越广泛的重视。目前以分析、优化和仿真为特征的CAE(Computer aided Engineering CAE)技术在世界范围内蓬勃发展。它通过先进的CAE技术快速有效地分析产品的各种特性、揭示结构各类参数变化对产品性能的响，进行设计方案的修改和调整，使产品达到性能和质量上的最优，原材料消耗最低。因此，基于计算机的分析、优化和仿真的CAE技术的研究和应用，是高质量、高水平、低成本产品设计与开发的保证，有限元法也必将在科技发展史上大放异彩。参考文献1AnnikaM,HakanC.TheDevelopmentofthe

22、FiniteElementMethod,Reporttvsm7040,LundUniversity,Sweden,19872ZienkiewiczOC.Thefiniteelementmethod.Thirdedition,Mcgraw-Hill,19773VichnevetskyR.ComputerMethodForPartialDiffrencialEquations,vol.1,Prenricehall,19814NielsO.Introductiontothefiniteelementmethod.PrenticeHall,19925CloughRW.Thefiniteelementmethodinplanestressanalysis.porceedings2ndasceconf.Onelfetorniccomputation,Pittsburgh,PA.,sep.8-9,19606林群.21世纪初科学发展趋势课题组.21世纪初科学发展趋势.北京:科学出版社。19967孔祥谦.有限单元法在传热学中的应用。第三版，科学出版社。1998专心-专注-专业