有限元法概念意义与应用.docx
《有限元法概念意义与应用.docx》由会员分享,可在线阅读,更多相关《有限元法概念意义与应用.docx(13页珍藏版)》请在冰豆网上搜索。
![有限元法概念意义与应用.docx](https://file1.bdocx.com/fileroot1/2022-12/15/712281b0-ca1e-4844-9e8c-d1e121641804/712281b0-ca1e-4844-9e8c-d1e1216418041.gif)
有限元法概念意义与应用
有限元法概论、意义与应用
班级:
2013信息姓名:
张正
学号:
2013040692
指导老师:
曾伟梁
摘要:
有限元法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
关键词:
有限元法;变分原理;加权余量法;函数。
Abstract:
Finiteelementmethodisbasedonthevariationalprincipleandtheweightedresidualmethod,thebasicideaistosolvethecomputationaldomainisdividedintoafinitenumberofnon-overlappingunits,eachunit,selectsomeappropriatefunctionforsolvingtheinterpolationnodepointsas,thedifferentialvariablesrewrittenoritsderivativebythevariablevalueoftheselectednodeinterpolationfunctionsconsistingoflinearexpressions,bymeansofvariationalprincipleorweightedresidualmethod,thediscretedifferentialequationstosolve.Differentformsofweightfunctionsandinterpolationfunctions,itconstitutesadifferentfiniteelementmethod.
Keywords:
Finiteelementmethod;variationalprinciple;weightedresidualmethod;function。
引言
随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。
这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。
例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响,看看是否会发生破坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡轮机叶片内的流体动力学参数,以提高其运转效率。
这些都可归结为求解物理问题的控制偏微分方程式往往是不可能的。
近年来在计算机技术和数值分析方法支持下发展起来的有限元分析(FEA,FiniteElementAnalysis)方法则为解决这些复杂的工程分析计算问题提供了有效的途径。
有限元法是一种高效能、常用的计算方法.有限元法在早期是以变分原理为基础发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中(这类场与泛函的极值问题有着紧密的联系)。
自从1969年以来,某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理场中,而不再要求这类物理场和泛函的极值问题有所联系.
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。
从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。
令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。
插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。
有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。
单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。
常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。
在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。
对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。
1、有限元法的孕育过程及诞生和发展
大约在300年前,牛顿和莱布尼茨发明了积分法,证明了该运算具有整体对局部的可加性。
虽然,积分运算与有限元技术对定义域的划分是不同的,前者进行无限划分而后者进行有限划分,但积分运算为实现有限元技术准备好了一个理论基础。
在牛顿之后约一百年,著名数学家高斯提出了加权余值法及线性代数方程组的解法。
这两项成果的前者被用来将微分方程改写为积分表达式,后者被用来求解有限元法所得出的代数方程组。
在18世纪,另一位数学家拉格郎日提出泛函分析。
泛函分析是将偏微分方程改写为积分表达式的另一途经。
在19世纪末及20世纪初,数学家瑞雷和里兹首先提出可对全定义域运用展开函数来表达其上的未知函数。
1915年,数学家伽辽金提出了选择展开函数中形函数的伽辽金法,该方法被广泛地用于有限元。
1943年,数学家库朗德第一次提出了可在定义域内分片地使用展开函数来表达其上的未知函数。
这实际上就是有限元的做法。
所以,到这时为止,实现有限元技术的第二个理论基础也已确立。
20世纪50年代,飞机设计师们发现无法用传统的力学方法分析飞机的应力、应变等问题。
波音公司的一个技术小组,首先将连续体的机翼离散为三角形板块的集合来进行应力分析,经过一番波折后获得前述的两个离散的成功。
20世纪50年代,大型电子计算机投入了解算大型代数方程组的工作,这为实现有限元技术准备好了物质条件。
1960年前后,美国的R.W.Clough教授及我国的冯康教授分别独立地在论文中提出了“有限单元”,这样的名词。
此后,这样的叫法被大家接受,有限元技术从此正式诞生。
1990年10月美国波音公司开始在计算机上对新型客机B-777进行“无纸设计”,仅用了三年半时间,于1994年4月第一架B-777就试飞成功,这是制造技术史上划时代的成就,其中在结构设计和评判中就大量采用有限元分析这一手段。
在有限元分析的发展初期,由于其基本思想和原理的“简单”和“朴素”,以至于许多学术权威都对其学术价值有所鄙视,国际著名刊物JournalofAppliedMechanics许多年来都拒绝刊登有关于有限元分析的文章。
然而现在,有限元分析已经成为数值计算的主流,不但国际上存在如ANSYS等数种通用有限元分析软件,而且涉及到有限元分析的杂志也有几十种之多。
二、有限元法的基本思想
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。
20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:
“有限元法=RayleighRitz法+分片函数”,即有限元法是RayleighRitz法的一种局部化情况。
不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的RayleighRitz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
有限元方法(FEM)的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。
从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。
令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。
插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。
有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。
单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。
常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。
在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。
对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。
三有限元单元法的基本思路
弹性力学解法的问题弹性力学解法的问题在于:
不论是应力函数解法数解法、扭转函数解
法、挠曲函数解法、还是基于最小势能原还是基于最小势能原理的瑞利-李兹等方法,其困难在于如何给出一个在全求解区给出一个在全求解区域上均成立的试探函数。
在有限单元法里在有限单元法里,这个问题通过定义分片插值的位移或应力函数得到了巧妙的解决。
对于任意单元对于任意单元(i,j,m),,以结点位移以结点位移(u,u,u)为
待定系数,可以给出该单元的插值函数:
线性代数方程组的求解在数学上是极其容易的。
也就是说有限元法通过单元离散和最小势能原理小势能原理,避开了微分方程直接求避微分方程直接求解在数学上的困难,把定解条件下的微分方程组的求解巧妙地转化为线性方程组的运算,实现了任何复杂弹性力学问题轻易分析计算。
四有限元单元法求解问题的的基本步骤
1.
(1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。
(2)区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元。
区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插值条件的插值函数作为单元基函数。
有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则。
(4)单元分析:
将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程。
(5)总体合成:
在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程。
(6)边界条件的处理:
一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。
对于自然边界条件,一般在积分表达式中可自动得到满足。
对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足。
(7)解有限元方程:
根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。
五求解计算结果的整理和有限元法后处理
有限元方程是一个线性代数方程组,一般有两大类解法,一是直接解法,二是迭代法。
直接法有高斯消元法和三角分解法,如果方程规模比较大时,可用分块解法和波前解法。
迭代法有雅可比迭代法、高斯-赛德尔迭代法和超松弛迭代法等。
通过选用合适的的求解法求解经过位移边界条件小处理的公式后,得到整体节点位移列阵,然后根据单元节点位移由几何矩阵和应力矩阵得到单元节点的应变和应力,对于非节点处的位移通过形函数插值得到,再由几何矩阵和应力矩阵求得相应的应变和应力。
应变要通过位移求导得到,精度一般要比位移差一些,尤其对于一次单元,应变和应力在整个单元内是常数,应变和应力的误差会比较大,尤其单元数比较少时,误差更大,因此对于应力和应变要进行平均化处理:
(1)绕节点平均法,即依次把围绕节点所有单元的应力加起来平均,以此平均应力作为该节点的应力。
(2)二单元平均法,即吧相邻的两单元的应力加以平均并以此作为公共边的节点出的应力。
整理并对有限元法计算结果进行后处理,一是要得到结构中关键位置力学量得数值(如最大位移、最大主应力和主应变,等效应力等);二是得到整个结构的力学量得分布(根据计算结果直接绘制位移分布图,应力分布图等)。
三是后处理要得到输入量和输出量之间的响应关系。
梁的有限元建模与变形分析
计算分析模型如图1所示。
NOTE:
要求选择不同形状的截面分别进行计算。
图1梁的计算分析模型
梁截面分别采用以下三种截面(单位:
m):
矩形截面:
圆截面:
工字形截面:
B=0.1,H=0.2R=0.1w1=0.2,w2=0.2,w3=0.2,
t1=0.05,t2=0.05,t3=0.05
梁的弹性模量为:
200GPa,泊松比为:
0.3。
试用ANSYS求该梁的挠度及应力分布(包含轴向应力),画出梁的弯矩图和剪力图。
(1)单元类型定义后,设置梁类型与实常数。
(2)设置材料属性。
(3)建立几何模型。
(4)创建线并划分网格与单元。
(5)施加载荷与约束。
(6)求解。
三、有限元的应用及其发展趋势
有限元的应用范围也是相当的广的。
它涉及到工程结构、传热、流体运动、电磁等连续介质的力学分析中,并在气象、地球物理、医学等领域得到应用和发展。
电子计算机的出现和发展是有限元法在许多实际问题中的应用变为现实,并具有广阔的前景。
国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。
其中最为著名的是由美国国家宇航局(NASA)在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。
该系统发展至今已有几十个版本,是目前世界上规模最大、功能最强的有限元分析系统。
从那时到现在,世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件,主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。
当今国际上FEA方法和软件发展呈现出以下一些趋势特征:
1.从单纯的结构力学计算发展到求解许多物理场问题
有限元分析方法最早是从结构化矩阵分析发展而来,逐步推广到板、壳和实体等连续体固体力学分析,实践证明这是一种非常有效的数值分析方法。
而且从理论上也已经证明,只要用于离散求解对象的单元足够小,所得的解就可足够逼近于精确值。
所以近年来有限元方法已发展到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算,最近又发展到求解几个交叉学科的问题。
例如当气流流过一个很高的铁塔时就会使铁塔产生变形,而塔的变形又反过来影响到气流的流动……这就需要用固体力学和流体动力学的有限元分析结果交叉迭代求解,即所谓"流固耦合"的问题。
2.由求解线性工程问题进展到分析非线性问题
随着科学技术的发展,线性理论已经远远不能满足设计的要求。
例如建筑行业中的高层建筑和大跨度悬索桥的出现,就要求考虑结构的大位移和大应变等几何非线性问题;航天和动力工程的高温部件存在热变形和热应力,也要考虑材料的非线性问题;诸如塑料、橡胶和复合材料等各种新材料的出现,仅靠线性计算理论就不足以解决遇到的问题,只有采用非线性有限元算法才能解决。
众所周知,非线性的数值计算是很复杂的,它涉及到很多专门的数学问题和运算技巧,很难为一般工程技术人员所掌握。
为此近年来国外一些公司花费了大量的人力和投资开发诸如MARC、ABQUS和ADINA等专长于求解非线性问题的有限元分析软件,并广泛应用于工程实践。
这些软件的共同特点是具有高效的非线性求解器以及丰富和实用的非线性材料库。
3.增强可视化的前置建模和后置数据处理功能
早期有限元分析软件的研究重点在于推导新的高效率求解方法和高精度的单元。
随着数值分析方法的逐步完善,尤其是计算机运算速度的飞速发展,整个计算系统用于求解运算的时间越来越少,而数据准备和运算结果的表现问题却日益突出。
在现在的工程工作站上,求解一个包含10万个方程的有限元模型只需要用几十分钟。
但是如果用手工方式来建立这个模型,然后再处理大量的计算结果则需用几周的时间。
可以毫不夸张地说,工程师在分析计算一个工程问题时有80%以上的精力都花在数据准备和结果分析上。
因此目前几乎所有的商业化有限元程序系统都有功能很强的前置建模和后置数据处理模块。
在强调"可视化"的今天,很多程序都建立了对用户非常友好的GUI(GraphicsUserInterface),使用户能以可视图形方式直观快速地进行网格自动划分,生成有限元分析所需数据,并按要求将大量的计算结果整理成变形图、等值分布云图,便于极值搜索和所需数据的列表输出。
4.与CAD软件的无缝集成
当今有限元分析系统的另一个特点是与通用CAD软件的集成使用,即在用CAD软件完成部件和零件的造型设计后,自动生成有限元网格并进行计算,如果分析的结果不符合设计要求则重新进行造型和计算,直到满意为止,从而极大地提高了设计水平和效率。
今天,工程师可以在集成的CAD和FEA软件环境中快捷地解决一个在以前无法应付的复杂工程分析问题。
所以当今所有的商业化有限元系统商都开发了和著名的CAD软件(例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley和AutoCAD等)的接口。
5.在Wintel平台上的发展
早期的有限元分析软件基本上都是在大中型计算机(主要是Mainframe)上开发和运行的,后来又发展到以工程工作站(EWS,EngineeringWorkStation)为平台,它们的共同特点都是采用UNIX操作系统。
PC机的出现使计算机的应用发生了根本性的变化,工程师渴望在办公桌上完成复杂工程分析的梦想成为现实。
但是早期的PC机采用16位CPU和DOS操作系统,内存中的公共数据块受到限制,因此当时计算模型的规模不能超过1万阶方程。
MicrosoftWindows操作系统和32位的IntelPentium处理器的推出为将PC机用于有限元分析提供了必需的软件和硬件支撑平台。
因此当前国际上著名的有限元程序研究和发展机构都纷纷将他们的软件移植到Wintel平台上。
下表列出了用ADINAV7.3版在PC机的WindowsNT环境和SGI工作站上同时计算4个工程实例所需要的求解时间。
从中可以看出最新高档PC机的求解能力已和中低挡的EWS不相上下。
为了将在大中型计算机和EWS上开发的有限元程序移植到PC机上,常常需要采用Hummingbird公司的一个仿真软件Exceed。
这样做的结果比较麻烦,而且不能充分利用PC机的软硬件资源。
所以最近有些公司,例如IDEAS、ADINA和R&D开始在Windows平台上开发有限元程序,称作"NativeWindows"版本,同时还有在PC机上的Linux操作系统环境中开发的有限元程序包。
在大力推广CAD技术的今天,从自行车到航天飞机,所有的设计制造都离不开有限元分析计算,有限元法在工程设计和分析中将得到越来越广泛的重视。
目前以分析、优化和仿真为特征的CAE(ComputeraidedEngineeringCAE)技术在世界范围内蓬勃发展。
它通过先进的CAE技术快速有效地分析产品的各种特性、揭示结构各类参数变化对产品性能的响,进行设计方案的修改和调整,使产品达到性能和质量上的最优,原材料消耗最低。
因此,基于计算机的分析、优化和仿真的CAE技术的研究和应用,是高质量、高水平、低成本产品设计与开发的保证。
有限元法也必将在科技发展史上大放异彩。
参考文献
1张晋红,吴凤林有限元法及其应用现状【J】。
机械管理开发,2009,4
2樊杰,赵铁,徐清洁有限元法的一般思想【J】。
科技情报开发与经济,2008,3
3杨桂通弹性力学【M】北京:
高等教育出版社,1998
4邵敏有限元法基本原理和数值方法[M]。
北京:
清华大学出版社,1997
5傅永华有限元法分析基础【M】。
武汉:
武汉大学出版社,2003