有限元分析基础教程
Fundamentals of Finite Element A nalysis
曾 攀
(清华大学)
全教程包括两大部分,共分 9 章;第一部分为有限元分析基本原理,包括第 1 章至第 5章,内容有:绪论、有限元分析过程的概要、杆梁结构分析的有限元方法、连续体结构分析的有限元方法、有限元分析中的若干问题讨论; 第二部分为有限元分析的典型应用领域,包括第 6 章至第 9 章,内容有:静力结构的有限元分析、结构振动的有限元分析、传热过程的有限元分析、弹塑性材料的有限元分析。本书以基本变量、基本方程、求解原理、单元构建、典型例题、MATLAB 程序及算例、ANSYS 算例等一系列规范性方式来描述有限元分析的力学原理、程序编制以及实例应用;给出的典型实例都详细提供有完整的数学推演过程以及ANSYS 实现过程。本教程的基本理论阐述简明扼要,重点突出,实例丰富,教程中的二部分内容相互衔接,也可独立使用,适合于具有大学高年级学生程度的人员作为培训教材,也适合于不同程度的读者进行自学;对于希望在 MATLAB 程序以及 ANSYS 平台进行建模分析的读者,本教程更值得参考。
本基础教程的读者对象:机械、力学、土木、水利、航空航天等专业的工程技术人员、科研工作者。
曾攀,男,1963 年生,海南省海口市人;1988 年在清华大学获博士学位,1988-1992年先后在大连理工大学和西南交通大学从事两站博士后研究(领域为计算力学),为国家杰出青年科学基金获得者(1998)、长江学者(2000)、德国“洪堡”学者(1994-1995) 、“新世纪百千万人才工程” 国家级人选。现为清华大学机械工程系主任、教授、博士生导师,为《机械工程学报》、《工程力学》、《塑性工程学报》等五个学术期刊的编委, 为上海交通大学振动冲击噪声国家重点实验室、华中科技大学塑性成形模拟与模具国家重点实验室学术委员会委员, 先后主持包括国家级重点基金项目、863 项目、霍英东基金项目等科研项目 30 多个,获教委科技进步二等奖、机械部一等奖、北京高等教育教学成果二等奖各一项,获国家发明专利授权两项;已出版及翻译学术著作有:《材料的概率疲劳损伤力学及现代结构分析原理》(曾攀, 1993) 、《有限元分析及应用》(曾攀,2004) 、《工程中的有限元方法》(T.R.Chandrupatla,曾攀译)(2006)、《有限元方法第 1 卷:基本原理》(O.C.Zienkiewicz(第 5 版),曾攀译)(2008) ,发表论文 100 多篇。主要从事计算力学、结构设计与分析、材料加工中的数值模拟等方面的研究。
有限元分析已经在教学、科研以及工程应用中成为重要而又普及的数值分析方法和工具;该基础教程力求提供具备现代特色的实用教程。在教材的内容体系上综合考虑有限元方法的力学分析原理、建模技巧、应用领域、软件平台、实例分析这几个方面, 按照教科书的方式深入浅出地叙述有限元方法,并体现出有限元原理“在使用中学习,在学习中使用”的交互式特点,在介绍每一种单元的同时,提供完整的典型推导实例、MATLAB 实际编程以及 ANSYS 应用数值算例,并且给出的各种类型的算例都具有较好的前后对应性,使学员在学习分析原理的同时,也进行实际编程和有限元分析软件的操作,经历实例建模、求解、分析和结果评判的全过程,在实践的基础上深刻理解和掌握有限元分析方法。
一本基础教材应该在培养学员掌握坚实的基础理论、系统的专业知识方面发挥作用, 因此,教材不但要提供系统的、具有一定深度的基础理论,还要介绍相关的应用领域,以给学员进一步学习提供扩展空间,本教程正是按照这一思路进行设计的;全书的内容包括两个部分,共分 9 章;第一部分为有限元分析基本原理,包括第 1 章至第 5 章,内容有:绪论、有限元分析过程的概要、杆梁结构分析的有限元方法、连续体结构分析的有限元方法、有限元分析中的若干问题讨论;第二部分为有限元分析的典型应用领域,包括第 6 章至第 9 章,内容有:静力结构的有限元分析、结构振动的有限元分析、传热过程的有限元分析、弹塑性材料的有限元分析。在基本原理方面, 以基本变量、基本方程、求解原理、单元构建等一系列规范的方式进行介绍;在阐述有限元分析与应用方面,采用典型例题、MATLAB 程序及算例、ANSYS 算例的方式,以体现出分析建模的不同阶段和层次,引导学员领会有限元方法的实质,还提供有大量的练习题。
本教程的重点是强调有限元方法的实质理解和融会贯通,力求精而透,强调学员综合能力(掌握和应用有限元方法)的培养,为学员亲自参与建模、以及使用先进的有限元软件平台提供较好的素材;同时,给学员进一步学习提供新的空间。
本教程力求体现以下特点。
(1)考虑教学适应性:强调对学员在数学原理、分析建模、软件应用几个方面的培养目标要求,注重学员在工程数值方面的基础训练,培养学员“使用先进软件+分析实际问题”的初步能力。
(2)考虑认知规律性:力求按照有限元分析方法的教学规律和认知规律,在教材中设计了“基本变量、基本方程、求解原理、单元构建”这样的模块;并体现出有限元原理“在使用中学习,在学习中使用”的交互式特点,在介绍每一种单元的同时,提供实用的 MATLAB实际编程和数值实例;在每一章还进行要点总结,给出典型例题,以引导学员领会有限元方法的实质,体现教材的启发性,有利于激发学员学习兴趣和便于自学。
(3)考虑结构完整性:本教程提供完整的教材结构:绪论、正文、典型例题、基于 MATLAB的编程算例与数值算例、具有一定深度的 ANSYS 算例、各章要点、习题、专业术语的英文标注、关键词中文和英文索引、参考文献,便于学员查阅。
(4)内容上的拓展性:除基本内容外,还介绍了较广泛的应用领域,包括:静力结构分析、结构振动分析、传热过程分析、弹塑性材料分析; 提供了有关的典型问题的建模详细分析过程,基本上反映了有限元分析在一些主要领域的应用状况及建模方法。
(5)编排上的逻辑性:本教程力求做到具有分明的层次和清楚的条理,在每一章中重点突出有限元方法的思想、数理逻辑及建模过程,强调相应的工程概念,提供典型例题及详解,许多例题可作为读者进行编程校验的标准考题(Benchmark),还提供了对应的 MATLAB 编程算例与 ANSYS 算例,特别是介绍了基于 APDL 参数化的 ANSYS 建模方法,并给出具体的实例,力求反映有限元分析的内在联系及特有思维方式。
由于作者的水平有限,所写的教程一定会出现不少错误,敬请读者提出批评。
曾 攀
于清华大学
2008-12-11
第 1 章 绪论
1.1 概况
有限元方法(finite element method)或有限元分析(finite element an alysis)[ 1][2],是求取复杂微分方程近似解的一种非常有效的工具,是现代数字化科技的一种重要基础性原理。将它用于在科学研究中,可成为探究物质客观规律的先进手段。将它应用于工程技术中, 可成为工程设计和分析的可靠工具。严格来说,有限元分析必须包含三个方面:(1)有限元方法的基本数学力学原理,(2)基于原理所形成的实用软件,(3)使用时的计算机硬件。随着现代计算机技术的发展,一般的个人计算机就能满足第(3)方面的要求;因此,本书的重点将在以上的第(1)和第(2)方面,将通过一些典型的实例来深入浅出地系统阐述有限元分析的基本原理,并强调原理的工程背景和物理概念;基于 MATLAB 平台来系统演示基于有限元原理的编程方法和过程;通过 ANSYS 分析平台来展示具体应用有限元方法的建模过程。
1.2 有限元方法的历史
有限元方法的思想最早可以追溯到古人的“化整为零 ”、“化圆为直”的作法,如“曹冲称象”的典故,我国古代数学家刘徽采用割圆法来对圆周长进行计算;这些实际上都体现了离散逼近的思想,即采用大量的简单小物体来“冲填”出复杂的大物体。
早在 1870 年,英国科学家 Rayleigh 就采用假想的“试函数”来求解复杂的微分方程, 1909 年 Ritz 将其发展成为完善的数值近似方法,为现代有限元方法打下坚实基础。
20 世纪 40 年代,由于航空事业的飞速发展,设计师需要对飞机结构进行精确的设计和计算,便逐渐在工程中产生了的矩阵力学分析方法;1943 年,Courant 发表了第一篇使用三角形区域的多项式函数来求解扭转问题的论文[3];1956 年波音公司的 Turner,Clough,Martin和 Topp 在分析飞机结构时系统研究了离散杆、梁、三角形的单元刚度表达式[4] ;1960 年Clough 在处理平面弹性问题,第一次提出并使用“有限元方法”(finite element method)的名称[5][6] ;1955 年德国的Argyris 出版了第一本关于结构分析中的能量原理和矩阵方法的书[7],为后续的有限元研究奠定了重要的基础,1967 年 Zienkiewicz 和 Cheung 出版了第一本有关有限元分析的专著;1970 年以后,有限元方法开始应用于处理非线性和大变形问题;我国的一些学者也在有限元领域做出了重要的贡献,如胡海昌于 1954 提出了广义变分原理[8] ,钱伟长最先研究了拉格朗日乘子法与广义变分原理之间关系[9],钱令希在 20 世纪五十年代就研究了力学分析的余能原理,冯康在 20 世纪六十年代就独立地、并先于西方奠定了有限元分析收敛性的理论基础。图 1-1 展示了有限元方法的发展过程。
随着计算机技术的飞速发展,基于有限元方法原理的软件大量出现,并在实际工程中发挥了愈来愈重要的作用;目前,专业的著名有限元分析软件公司有几十家,国际上著名的通用有限元分析软件有 ANSYS,ABAQUS,MSC/NASTRAN,MSC/MARC,ADINA,ALGOR, PRO/MECHANICA ,IDEAS,还有一些专门的有限元分析软件,如 LS-DYNA ,DEFORM, PAM-STAMP, AUTOFORM ,SUPER-FORGE 等;国际上著名的主要有限元分析软件状况见表 1-1。有关有限元分析的学术论文,每年也不计其数,学术活动非常活跃,表 1-2 列出的是刊登有限元分析论文的常见学术期刊。




下面介绍几位有限元方法的创始人的情况。

Richard Courant(美国数学家,1888.01.08-1972.01.27,出生地:德国Lublinitz), 1910 年在哥廷根大学获得博士学位,1928 年创办了应用数学研究所,并在 1928~1933 年担任所长,开始应用数学的研究,1936 年去纽约大学创立了数学研究所,1964 年该研究所被命名为Courant数学科学研究所,出版了数学物理方法教科书,成为最有影响的书籍之一,还出版了科普读物“什么是数学? ”,至今还再版发行,Courant的名字还与有限元方法紧密联系在一起,是他给这种数值求解偏微分方程以坚实的理论基础, “Courant–Friedrichs–Lewy条件”以及“Courant最小值原理”都是以他的名字命名的重要定理。

John Argyris(德国人, 1913.08.19-2004.04.02,出生地:希腊Volos)是公认的计算科学、航空航天、流体力学领域的杰出专家,被誉为有限元分析的发明人和创始者之一,所创立的力学分析的矩阵位移方法为有限元方法的前身,被称为世界上最先进的数学工具,先后因在有限元方法以及混沌理论方面的贡献而获得菲利普王子奖章、美国最高荣誉的爱因斯坦奖,获得包括 5 个G8 国家在内的由 16 发达国家授予的杰出科技与工程奖;Argyris于 1950 年代在斯图加特大学创立计算机应用研究所,并担任所长近 40年,立至于有限元分析的研究,开发了几个著名的在欧洲广泛使用的有限元分析商业化软件,在有限元方法的应用方面也作出了杰出的贡献。

Olgierd Cecil Zienkiewiz (英国人,1921.05. 18-,出生地:英国Caterham) 。 英国Swansea大学的荣誉退休教授,是该校工程数值方法研究所的原主任,现在仍然是西班牙巴塞罗那Calalunya技术大学工程数值方法的UNESCO主席。从 1961 至 1989 年,担任Swansea大学土木工程系的主任,使该系成为有限元研究的重要中心之一。在 1968 年,创办了International Journal for Numerical Methods in Engineering杂志并任主编,该杂志至今仍然是该领域的主要刊物。他被授予 24 个荣誉学位和多种奖励。Zienkiewiez教授还是 5 所科学院的院士,这是对他在有限元方法领域的奠基性发展和贡献的赞誉。1978 年,成为皇家科学院和皇家工程院的院士;并先后被选为美国工程院的外籍院(1981),波兰科学院院士(1985),中国科学院院士(1998)和意大利国家科学院院士(1999)。1967 年,他出版了有限元领域的第一本专著《有限元方法》,该书目前也出版到第 5 版,可参见《有限元方法第 1 卷:基本原理》(O.C.Zienkiewicz, R.L.Taylor(第 5 版),曾攀译)(清华大学出版社,2008)。
1.3 有限元分析的作用
基于功能完善的有限元分析软件和高性能的计算机硬件对设计的结构进行详细的力学分析,以获得尽可能真实的结构受力信息,就可以在设计阶段对可能出现的各种问题进行安全评判和设计参数修改,据有关资料,一个新产品的问题有 60%以上可以在设计阶段消除,甚至有的结构的施工过程也需要进行精细的设计,要做到这一点,就需要类似有限元分析这样的分析手段。
下面举出几个涉及土木工程、车辆工程、航空工程以及生物工程的实例。
北京奥运场馆的鸟巢由纵横交错的钢铁枝蔓组成,它是鸟巢设计中最华彩的部分,见图1-2,也是鸟巢建设中最艰难的。看似轻灵的枝蔓总重达 42000 吨,其中,顶盖以及周边悬空部位重量为 14000 吨,在施工时,采用了 78 根支柱进行支撑,也就是产生了 78 个受力区域,在钢结构焊接完成后,需要将其缓慢而又平稳地卸去,让鸟巢变成完全靠自身结构支撑;因而,支撑塔架的卸载,实际上就是对整个钢结构的加载,如何卸载?需要进行非常详细的数值化分析,以确定出最佳的卸载方案。2006 年 9 月 17 日成功地完成了整体钢结构施工的最后卸载。

图 1-3 给出的是现代列车车厢整体结构的有限元分析模型[ 10];图 1-4 表明的是空客 A350后机身第 19 框的设计与有限元分析过程[ 11];图 1-5 为人体肩部区域的骨胳有限元分析模型以及计算结果[ 12] ,并与实际测试结果进行了比较。所有这些表明了有限元方法的广泛适应性以及在科学研究与重大工程中的重要性。

本章先通过一个简单的实例,采用直接的推导方法,逐步展示有限元分析的基本流程,从中可以了解有限元方法的思路形成过程,以及如何由具体的求解步骤归纳出一种通用的标准求解方法。
2.1 有限元分析的目的和概念
任何具有一定使用功能的构件(称为变形体(deformed body))都是由满足要求的材料所制造的,在设计阶段,就需要对该构件在可能的外力作用下的内部状态进行分析,以便核对所使用材料是否安全可靠,以避免造成重大安全事故。描述可承力构件的力学信息一般有三类:
(1) 构件中因承载在任意位置上所引起的移动(称为位移(displacement));
(2) 构件中因承载在任意位置上所引起的变形状态(称为应变(strain));
(3) 构件中因承载在任意位置上所引起的受力状态(称为应力(stress));
若该构件为简单形状,且外力分布也比较单一,如:杆、梁、柱、板就可以采用材料力学的方法,一般都可以给出解析公式,应用比较方便;但对于几何形状较为复杂的构件却很难得到准确的结果,甚至根本得不到结果。
有限元分析的目的:针对具有任意复杂几何形状变形体,完整获取在复杂外力作用下它内部的准确力学信息,即求取该变形体的三类力学信息(位移、应变、应力)。
在准确进行力学分析的基础上,设计师就可以对所设计对象进行强度(strength) 、刚度(stiffness)等方面的评判,以便对不合理的设计参数进行修改,以得到较优化的设计方案;然后,再次进行方案修改后的有限元分析,以进行最后的力学评判和校核,确定出最后的设计方案。
图 2-1 给出一个针对大型液压机机架的设计过程以及采用有限元分析的状况。

为什么采用有限元方法就可以针对具有任意复杂几何形状的结构进行分析,并能够得到准确的结果呢?这时因为有限元方法是基于“离散逼近(discretized approximation)”的基本策略,可以采用较多数量的简单函数的组合来“近似”代替非常复杂的原函数。
一个复杂的函数,可以通过一系列的基底函数(base function)的组合来“近似”,也就是函数逼近,其中有两种典型的方法:(1)基于全域的展开(如采用傅立叶级数展开),以及(2)基于子域(sub-domain)的分段函数(pieces function)组合(如采用分段线性函数的连接);下面,仅以一个一维函数的展开为例说明全域逼近与分段逼近的特点。
【典型例题】2.1(1) 一个一维函数的两种展开方式的比较

比较以上两种方式的特点,可以看出,第一种方式所采用的基本函数φi (x∈[x0 , xL])非常复杂,而且是在全域上[x0 , xL]定义的, 但它是高次连续函数,一般情况下,仅采用几个基底函数就可以得到较高的逼近精度 ; 而第二种方式所采用的基本函数
ai + bix(x∈[xi , xi+1])非常简单,而且是在子域上[xi , xi+1]定义的,它通过各个子域组合出全域[x0 , xL],但它是线性函数,函数的连续性阶次较低,因此需要使用较多的分段才能得到较好的逼近效果,则计算工作量较大。
对于第一种的函数逼近方式, 就是力学分析中的经典瑞利-里兹方法(Rayleigh-Ritz principle)的思想,而针对第二种的函数逼近方式,就是现代力学分析中的有限元方法的思想,其中的分段就是“单元”的概念。
基于分段的函数描述具有非常明显的优势: (1)可以将原函数的复杂性“化繁为简 ”,使得描述和求解成为可能,(2)所采用的简单函数可以人工选取,因此,可取最简单的线性函数,或取从低阶到高阶的多项式函数,(3)可以将原始的微分求解变为线性代数方程。但分段的做法可能会带来的问题有:(1)因采用了“化繁为简”,所采用简单函数的描述的能力和效率都较低,(2)由于简单函数的描述能力较低,必然使用数量众多的分段来进行弥补,因此带来较多的工作量。
综合分段函数描述的优势和问题,只要采用功能完善的软件以及能够进行高速处理的计算机,就可以完全发挥“化繁为简”策略的优势,有限元分析的概念就在于此。
原资料见附件
免责声明:
本页面/内容部分素材来源于互联网公 开 信 息,旨在传递更多信息,不代表本平台立场。
版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。
本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险。