浅析有限元方法的发展与应用
浅析有限元方法的发展与应用
摘要:1965年“有限元”这个名词第一次在我国出现,到今天有限元在工程上得到广泛应用,经历了三十多年的发展历史,理论和算法都已经日趋完善。有限元法(Finite Element Method,简写为FEM)是求解微分方程的一种非常有效的数值计算方法,用这种方法进行波动数值模拟受到越来越多的重视。
关键字:有限元法发展应用
Abstract:1965 the term "finite element" first appeared in our country, to this day the finite ele ment is widely used in engineering, has experienced more than 30 years of development history, the ory and algorithm have been becoming more complete.FEM (Finite Element Method, abbreviated a s FEM) is a very effective to solve the differential equation of numerical calculation Method of wav e numerical simulation by using this Method is more and more attention.
Keywords: finite element method development Application
绪论
有限元法是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。它是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。
一、有限元的发展历程
有限元法的发展历程可以分为提出(1943)、发展(1944-1960)和后期(1961-二十世纪九十年代)三个阶段。有限元法是受内外动力的综合作用而产生的。
1943年，柯朗在《美国数学学会公报》(Bulletin of The American Mathematical Society)上发表了《平衡和振动问题的变分解法》 (Variational Methods for The Solution of Problems of Equilibrium And Vibration)一文，这篇文章实际上是他1941年在美国数学学
会演讲的书面稿，在其中柯朗提出了有限元法的核心思想。大约与柯朗同时，工程师阿格瑞斯在另一个领域独立地提出了有限元法。柯朗和阿格瑞斯各自在数学和工程学领域独立提出了有限元法，他们分别开创了有限元法的数学传统和工程学传统。
有限元法被提出来以后，经过一段时间的沉寂期，在二十世纪五十年代和六十年代初有了很大的发展。主要表现为在代数表达形式、单元划分、单元类型选择和解的收敛性研究上取得的突破。1960年，克劳夫在《平面应力分析中的有限元》(The Finite Element in Plan Stress Analysis)的论文中，第一次从数学上说明了将定义域划分成有限的单元能够成功的原因:他表明对一些特定类型的单元来说，随着单元尺寸的减小，近似解将收敛到精确解，这就在某些情况下证明了有限元法的收敛性。并第一次提出了“有限元法”这个名称，这个名称一直沿用至今，标志着有限元法早期发展阶段的结束。
有限元法后期阶段的发展有国外和国内两条线索。在国外的发展表现为: 第一，建立了严格的数学和工程学基础；第二，应用范围扩展到了结构力学以外的领域；第三，收敛性得到了进一步研究，形成了系统的误差估计理论；第四，发展起了相应的商业软件包。
在国内，我国数学家冯康在特定的环境中独立于西方提出了有限元法。1965年，他发表论文《基于变分原理的差分格式》，标志着有限元法在我国的诞生。冯康的这篇文章不但提出了有限元法，而且初步发展了有限元法。他得出了有限元法在特定条件下的表达式，独创了“冯氏大定理”并且初步证明了有限元法解的收敛性。虽然冯康创造的有限元法不成熟，但他能在当时的条件下独立提出有限元法已十分不易。对于他的这项成就，国内外专家学者和国家领导人都有很高的评价。
半个世纪来，有限单元法蓬勃发展，不仅已经成为结构分析中必不可少的工具，而且成为现象分析的一种手段。其应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、渗流与固结理论、热传导与热应力问题、磁场问题以及建筑声学与噪音问题。不仅涉及稳态场问题，还涵盖材料非线性、几何非线性、时间维问题和断裂力学等。
已出现多种新单元（先后有等参元、高次元、不协调元、拟协调元、杂交元、样条元、边界元、罚单元，还有半解析的有限条等不同单元）和求解方法（如半带宽与变带宽消去法、超矩阵法、波前法、子结构法、子空间迭代法等）。能解决各种复杂耦合问题的软件和软件系统不断涌现。对网格自动剖分和网格自适应过程的研究，大大加强了有限元法的解题能力，使有限单元法逐渐趋于成熟。有限元法作为一种离散化的数值解法，也已成为应用数学的一
个新的分支。
二、有限元法解题思路
有限元法分析计算的思路和做法可归纳如下：1.物体离散化；2.单元特性分析；3.单元组集；4.求解未知结点位移。
1.物体离散化
将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。
2.单元特性分析
(1)选择位移模式，通常，有限元法我们就将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数。
(2)分析单元的力学性质，据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵。
(3)计算等效节点力，物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上的力。
3.单元组集
利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程。
4.求解未知结点位移
最后利用已求出结点位移计算各个单元的应力，并经后处理软件整理、显示计算结果。有限元法是应用局部的近似解来建立整个定义域的解的一种方法。先把注意力集中在单个单元上，进行上述所谓的单元分析。基本前提是每一单元要尽可能小，以致其边界值在整个边界上的变化也是小的。这样，边界条件就能取某一在结点间插值的光滑函数来近似，在单元内也容易建立简单的近似解。因此，比起经典的近似法，有限元法具有明显的优越性。
三、有限元方法的实施过程
有限元法在计算机软件中应用尤为突出，可以说有限元法与生俱来就是为了使用计算来进行复杂的计算。有限元方法的实施过程可以分为三个步骤：
1.前处理。将整体结构或其一部分简化为理想的数学模型，用离散化的网格代替连续的实体结构。
2.计算分析。分析计算结构的受力、变形及特性。
3.将计算结果进行整理和归纳。
对于有限元程序使用者而言，第一步和第三步的工作量最大，一个有限元程序的好坏，在很大程度上取决于第一步的前处理和第三步的后处理功能是否强大。
前处理：对于第一步的前处理而言，要根据计算的目的和所关心的区域，将结构模型化、离散化。需要给出下列信息：
(1)节点的空间位置。
(2)单元与节点的连接信息。
(3)结构的物质特性和材料参数。
(4)边界条件或约束。
(5)各类载荷。
在构成离散模型时，为了使模型较为合理，必须遵循以下的原则：
1.使计算模型尽量简化，以减少计算时间和容量，但又必须抓住主要因素以不影响计算精度。
2.在所关心的区域加密计算网格。
后处理：有限元计算是一种大规模的科学计算，其特点是除了要花费巨大的计算机处理能力外，在计算过程中还会产生巨大数量的数字信息。只有在计算输出信息进行仔细分析理解之后，才能洞察计算中发生的情况和问题，才能获得对被研究对象的认识和见解。
在大多数情况下，被研究的对象都是三维介质中的场分布问题（应力分布、位移分布、压力分布、电场分布等），即所谓的“四维”问题。鉴于其计算结果分析的复杂性，人们提出了科学计算可视性的要求，即把四维的数据进行图形处理或称为可视化处理，使人们能够看到场的分布图象，从图象上直接进行分析、判断来获得有用的结论。这大大加快和加深了人们对计算对象的物理变化过程的认识，发现通常通过数值信息发现不了的现象，甚至获得意料之外的启发和灵感，从而缩短了研究和设计周期，提高了效率，获得更多的结果。
四、有限元的应用是在于多物理场的耦合
随着计算机技术的迅速发展，在工程领域中，有限元分析（FEA）越来越多地用于仿真模拟，来求解真实的工程问题。这些年来，越来越多的工程师、应用数学家和物理学家已经证明这种采用求解偏微分方程（PDE）的方法可以求解许多物理现象，这些偏微分方程可以用来描述流动、电磁场以及结构力学等等。有限元方法用来将这些众所周知的数学方程转化为近似的数字式图象。
早期的有限元主要关注于某个专业领域，比如应力或疲劳，但是，一般来说，物理现象都不是单独存在的。例如，只要运动就会产生热，而热反过来又影响一些材料属性，如电导率、化学反应速率、流体的粘性等等。这种物理系统的耦合就是我们所说的多物理场，分析起来比我们单独去分析一个物理场要复杂得多。很明显，我们现在需要一个多物理场分析工具。在上个世纪90年代以前，由于计算机资源的缺乏，多物理场模拟仅仅停留在理论阶段，有限元建模也局限于对单个物理场的模拟，最常见的也就是对力学、传热、流体以及电磁场的模拟。看起来有限元仿真的命运好像也就是对单个物理场的模拟。
现在这种情况已经开始改变。经过数十年的努力，计算科学的发展为我们提供了更灵巧简洁而又快速的算法，更强劲的硬件配置，使得对多物理场的有限元模拟成为可能。新兴的有限元方法为多物理场分析提供了一个新的机遇，满足了工程师对真实物理系统的求解需要。
下面通过几个例子来展示多物理场的有限元分析在未来的一些潜在应用:
压电扩音器（Piezoacoustic transducer）可以将电流转换为声学压力场，或者反过来，将声场转换为电流场。这种装置一般用在空气或者液体中的声源装置上，比如相控阵麦克风，超声生物成像仪，声纳传感器，声学生物治疗仪等，也可用在一些机械装置比如喷墨机和压电马达等。压电扩音器涉及到三个不同的物理场：结构场，电场以及流体中的声场。只有具有多物理场分析能力的软件才能求解这个模型。压电材料选用PZT5-H晶体，这种材料在压电传感器中用得比较广泛。在空气和晶体的交界面处，将声场边界条件设置为压力等于结构场的法向加速度，这样可以将压力传到空气中去。另外，晶体域中又会因为空气压力对其的影响而产生变形。仿真研究了在施加一个幅值200V，震荡频率为300 KHz的电流后，晶体产生的声波传播。这个模型的描述及其完美的结果表明在任何复杂的模型下，我们都可以用一系列的数学模型进行表达，进而求解。
将基片的电磁、电阻以及传热行为耦合起来需要一个真正的多物理场分析工具。一个典型的应用是在半导体的加工和退火的工艺中，有一种利用感应加热的热壁熔炉，它用来让半导体晶圆生长，这是电子行业中的一项关键技术。
例如，金刚砂在2000℃的高温环境下可以取代石墨接收器，接收器由功率接近10KW的
射频装置加热。在如此高温下要保持炉內温度的均匀，炉腔的设计至关重要。经过多物理场分析工具的分析，发现热量主要是通过辐射的方式进行传播的。在模型內不仅可以看到晶圆表面温度的分布，还可以看到熔炉的石英管上的温度分布。在电路设计中，影响材料选择的重要方面是材料的耐久性和使用寿命。电器小型化的趋势使得可在电路板上安装的电子元件发展迅猛。众所周知，安装在电路板上的电阻以及其他一些元件会产生大量的热，进而可能使得元件的焊脚处产生裂缝，最后导致整个电路板报废。多物理场分析工具可以分析出整个电路板上热量的转移，结构的应力变化以及由于温度的上升导致的变形。这样做可以用来提升电路板设计的合理性以及材料选择的合理性。
计算机能力的提升使得有限元分析由单场分析到多场分析变成现实，未来的几年内，多物理场分析工具将会给学术界和工程界带来震惊。单调的“设计－校验”的设计方法将会慢慢被淘汰，虚拟造型技术将让你的思想走得更远，通过模拟仿真将会点燃创新的火花。
五、总结
本文阐述了有限元法的发展历史、解题思路以及对未来发展应用的展望，总结了有限元法在计算机程序中的实现过程。使大家对有限元法的理论有了进一步的了解，并使有限元法这高级计算方法更接近我们的生活，更容易地理解到有限元法的实施过程。随着科技的发展与软件的更新，有限元法将得到更广泛的应用。
参考文献
[1]高攀,黄放.有限元方法的发展状况和应用[J].电机技术.1999年02
[2]李昱君.语境视角下的有限元法发展史[J].山西大学，2008-06-01
[3]刘英魁.有限元分析的发展趋势作者[J].中国新技术新产品.2009年06期
[4]张永刚.有限元法发展及其应用[J]. 科技情报开发与经济.2007年11期
[5]董志高，吴继敏.岩土工程有限元软件应用与发展[J].工程地质计算机应用.2010年3期
[6]王崧.有限元分析-ansys理论与应用[J].电子出版社.2009-04
[7]张洪伟，张庆生.非线性有限元分析方法[J].水利水电出版社.2013-04
[8]梁醒培，王辉.应用有限元分析[J].清华大学出版社.2010-06
[9]武建华.有限元分析基础[J].重庆大学出版社.2007-03
[10]黄艾香.有限元理论与方法[J].科学出版社.2009-05
非线性有限元方法及实例分析
非线性有限元方法及实例分析 梁军 河海大学水利水电工程学院，南京（210098） 摘 要：对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究，讨论了非线性计算的迭代收敛准则，并利用非线性有限元方法分析了一个钢棒单轴拉伸的实例。 关键词：非线性有限元，方程组求解，实例分析 1引 言 有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题，其应用范围从固体到流体，从静力到动力，从力学问题到非力学问题。有限元的线性分析已经设计工具被广泛采用。但对于绝大多数水利工程中遇到的实际问题如地下洞室等，将其作为非线性问题加以考虑更符合实际情况。根据产生非线性的原因，非线性问题主要有3种类型[1]： 1．材料非线性问题（简称材料非线性或物理非线性） 2．几何非线性问题 3．接触非线性问题（简称接触非线性或边界非线性） 2 非线性方程组的求解 在非线性力学中，无论是哪一类非线性问题，经过有限元离散后，它们都归结为求解一个非线性代数方程组[2]： ()()()00 021212211=… …==n n n n δδδψδδδψδδδψΛΛΛ （1.1） 其中n δδδ,,,21Λ是未知量，n ψψψ,,,21Λ是n δδδ,,,21Λ的非线性函数，引用矢量记 号 []T n δδδδΛ21＝ （1.2） []T n ψψψψΛ21= （1.3） 上述方程组（1.1）可表示为 ()0=δψ （1.4） 可以将它改写为 ()()()0=?≡?≡R K R F δδδδψ （1.5） 其中()δK 是一个的矩阵，其元素 是矢量的函数，n n ×ij k R 为已知矢量。在位移有限 元中，δ代表未知的结点位移，()δF 是等效结点力，R 为等效结点荷载，方程()0=δψ表示结点平衡方程。 在线弹性有限元中，线性方程组
有限元分析理论基础
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
有限单元法与有限元分析
有限单元法与有限元分析 1.有限单元法 在数学中，有限元法（FEM，Finite Element Method）是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解，每个子区域都成为简单的部分，这种简单部分就称作有限元。它通过变分方法，使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想，有限元法包含了一切可能的方法，这些方法将许多被称为有限元的小区域上的简单方程联系起来，并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 随着电子计算机的发展，有限单元法是迅速发展成一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 1.1.有限元法分析本质 有限元法分析计算的本质是将物体离散化。即将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 1.2.特性分析 1）选择位移模式： 在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如
第9章 非线性问题的有限单元法
第9章非线性问题的有限单元法 9.1 非线性问题概述 前面章节讨论的都是线性问题，但在很多实际问题中，线弹性力学中的基本方程已不能满足，需要用非线性有限单元法。非线性问题的基本特征是变化的结构刚度，它可以分为三大类：材料非线性、几何非线性、状态非线性。 1. 材料非线性(塑性, 超弹性, 蠕变) 材料非线性指的是材料的物理定律是非线性的。它又可分为非线性弹性问题和非线性弹塑性问题两大类。例如在结构的形状有不连续变化（如缺口、裂纹等）的部位存在应力集中，当外载荷到达一定数值时该部位首先进入塑性，这时在该部位线弹性的应力应变关系不再适用，虽然结构的其他大部分区域仍保持弹性。 2. 几何非线性(大应变, 大挠度, 应力刚化) 几何非线性是有结构变形的大位移引起的。例如钓鱼杆，在轻微的垂向载荷作用下，会产生很大的变形。随着垂向载荷的增加，杆不断的弯曲，以至于动力臂明显减少，结构刚度增加。 3. 状态非线性(接触, 单元死活) 状态非线性是一种与状态相关的非线性行为。例如，只承受张力的电缆的松弛与张紧；轴承与轴承套的接触与脱开；冻土的冻结与融化。这些系统的刚度随着它们状态的变化而发生显著变化。 9.2 非线性有限元问题的求解方法 对于线性方程组，由于刚度方程是常数矩阵，可以直接求解，但对于非线性方程组，由于刚度方程是某个未知量的函数则不能直接求解。以下将简要介绍借助于重复求解线性方程组以得到非线性方程组解答的一些常用方法。 1.迭代法 迭代法与直接法不同，它不是求方程组的直接解，而是用某一近似值代人，逐步迭代，使近似值逐渐逼近，当达到允许的规定误差时，就取这些近似值为方程组的解。 与直接法相比，迭代法的计算程序较简单，但迭代法耗用的机时较直接法长。它不必存贮带宽以内的零元素，因此存贮量大大减少，且计算中舍入误差的积累也较小。以平面问题 为例，迭代法的存贮量一般只需直接法的14左右。在求解非线性方程组时，一般采用迭代 法。 2. 牛顿—拉斐逊方法 ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。 一种近似的非线性救求解是将载荷分成一系列的载荷增量，即逐步递增载荷和平衡迭代。它可以在几个载荷步内或者在一个载荷步的几个子步内施加载荷增量。在每一个增量的
有限元法及其在工程中的应用
机械与汽车学院 曹国强 主要内容： 1、有限元法的基本思想。 2、结构力学模型的简化和结构离散化。 3、有限元法的实施过程。 一、有限元法的基本思想 有限元法是随着计算机的发展而发展起来的一种有效的数值方法。其基本思想是：将连续的结构分割成数目有限的小单元体（称为单元），这些小单元体彼此之间只在数目有限的指定点（称为节点）上相互连接。用这些小单元体组成的集合体来代替原来的连续结构。再把每个小单元体上实际作用的外载荷按弹性力学中的虚功等效原理分配到单元的节点上，构成等效节点力，并按结构实际约束情况决定受约束节点的约束。这一过程称为结构的离散化。其次，对每个小单元体选择一个简单的函数来近似地表示其位移分量的分布规律，并按弹性力学中的变分原理建立起单元节点力和节点位移之间的关系（单元刚度方程），最后，把全部单元的节点力和节点位移之间的关系组集起来，就得到了一组以结构节点位移为未知量的代数方程组（总体刚度方程），同时考虑结构的约束情况，消去那些结构节点位移为零的方程，再由最后的代数方程组就可求得结构上有限个离散节点的各位移分量。求得了结构上各节点的位移分量之后，即可按单元的几何方程和物理方程求得各单元的应变和应力分量。 有限元法的实质就是把具有无限个自由度的连续体，理想化为有限个自由度的单元的集合体，使问题简化为适合于数值解法的结构型问题。 经典解法（解析法）与有限元法的区别 解析法 { } 建立一个描述连续体性质的偏微分方程组 有限元解法 连续体 数目增加到∞ 大小趋于0 微元 有限元 离散化 （单元分析）集合 总体分析 求得近似解
二、结构力学模型的简化和结构离散化 （一）结构力学模型的简化 用有限元法研究实际工程结构问题时，首先要从工程实际问题中抽象出力学模型，即要对实际问题的边界条件、约束条件和外载荷进行简化，这种简化应尽可能地反映实际情况，不至于使简化后的解答与实际差别过大，但也不要带来计算上的过分复杂，在力学模型的简化过程中，必须判断实际结构的问题类型，是二维问题还是三维问题。如果是平面问题，是平面应力问题，还是平面应变问题。同时还要搞清楚结构是否对称，外载荷大小和作用位置，结构的几何尺寸和力学参数（弹性模量E、波松比μ等）。 （二）结构的离散化 将已经简化好的结构力学模型划分成只在一些节点连续的有限个单元，把每个单元看成是一个连续的小单元体，各单元之间只在一些点上互相联结，这些点称作节点，每个单元体称为一个单元。用只在节点处连接的单元的集合体代替原来的连续结构，把外载荷按虚功等效原理移置到有关受载的节点上，构成节点载荷，把连续结构进行这样分割的过程称为结构的离散化。现举例说明。 设一平面薄板，中间有一个园孔，其左端固定，右端受面力载荷q，试对其进行有限元分割和力学模型简化。
第八章几何非线性问题的有限元法
第八章 几何非线性问题的有限元法 引言 前面各章所讨论的问题都是在小变形假设的前提下进行的，即假定物体所发生的位移远小于物体自身的几何尺寸，应变远小于1。在此前提下，建立物体或微元体的平衡条件时可以不考虑物体的位置和形状（简称位形）的变化，因此在分析中不必区别变形前后位形的差别，且应变可用一阶无穷小的线性应变表达。 实际上，上述假设有时是不成立的。即使实际应变可能是小的，且不超过材料的弹性极限，但如果需要精确地确定位移，就必须考虑几何非线性，即平衡方程应该相对于变形后的位置得出，而几何关系应该计及二次项。例如平板大挠度理论中，由于考虑了中面内的薄膜应力，求得的挠度比小挠度理论的结果有显著的减低。再如在结构稳定性问题中，当载荷达到一定数值后，挠度比线性解答予示的结果更剧烈地增加，并且确实存在承载能力随继续变形而减低的现象。在冷却塔、薄壁结构及其它比较细长的结构中，几何非线性分析都显得十分重要。 几何非线性问题可以分为以下几种类型： （1）大位移小应变问题。一般工程结构所遇到的几何非线性问题大多属于这一类。例如高层建筑或高耸构筑物以及大跨度网壳等结构的分析常需要考虑到结构大位移的影响。 （2）大位移大应变问题，如金属压力加工中所遇到的问题就属于这一类型。 （3）结构的变形引起外载荷大小、方向或边界支承条件的变化等。 结构的平衡实际上是在结构发生变形之后达到的，对于几何非线性问题来说，平衡方程必须建立在结构变形之后的状态上。为了描述结构的变形需要设置一定的参考系统。一种做法是让单元的局部坐标系始终固定在结构发生变形之前的位置，以结构变形前的原始位形作为基本的参考位形，这种分析方法称作总体的拉格朗日（Ｌagrange ）列式法；另一种做法是让单元的局部坐标系跟随结构一起发生变位，分析过程中参考位形是不断被更新的，这种分析方法称作更新的拉格朗日列式法。 本章首先对几何非线性问题作一般性讨论，从中导出经典的线性屈曲问题的公式；然后建立平板大挠度问题和壳体的大位移（及大转动）分析的有限方法公式；接着还给出了大应变及大位移的一般公式，最后还详细讨论了杆系结构几何非线性问题的有关公式。在讨论中我们采用总体的拉格朗日列式法，但对杆系结构，为应用方便我们给出了两种列式法的公式。 & 一般性讨论 理论基础 无论是对于何种几何非线性问题，虚功原理总是成立的。由虚功原理，单元的虚功方程可以写成如下的形式 {}{}{}{}0=-???**v e eT e eT F dv δσε （） 其中{}F 为单元节点力向量，{}e *ε为单元的虚应变，{}e *δ为节点虚位移向量。 增量形式的应变一位移关系可表示为 {}[] {}e e d B d δε= （）
有限元方法的发展及应用
有限元方法的发展及应用 摘要：有限元法是一种高效能、常用的计算方法。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描 述的各类物理场中。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法或最小二乘法等同样获得了有限元方程，因而有限元法可应用于 以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值 问题有所联系。基本思想：由解给定的泊松方程化为求解泛函的极值问题。 1有限元法介绍 1.1有限元法定义 有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后再求解。它是起源于20世纪50年代末60年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。 有限元法的基本思想是将求解域看成是由许多称为有限元的小的互连子域 组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总 的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而 是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得 到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行 之有效的工程分析手段。有限元法最初应用在工程科学技术中,用于模拟并且解 决工程力学、热学、电磁学等物理问题。 1.2有限元法优缺点 有限元方法是目前解决科学和工程问题最有效的数值方法，与其它数值方 法相比，它具有适用于任意几何形状和边界条件、材料和几何非线性问题、容 易编程、成熟的大型商用软件较多等优点。 （1）概念浅显，容易掌握，可以在不同理论层面上建立起对有限元法的理解，既可以通过非常直观的物理解释来理解，也可以建立基于严格的数学理论 分析。 （2）有很强的适用性，应用范围极其广泛。它不仅能成功地处理线性弹性
钢筋混凝土梁非线性有限元分析方法
第28卷第1期 V ol.28 No.1 工 程 力 学 2011年 1 月 Jan. 2011 ENGINEERING MECHANICS 82 ——————————————— 收稿日期：2009-06-19；修改日期：2010-03-11 基金项目：国家科技支撑计划项目(2006BA904B03) 作者简介：*周凌远(1968―)，男，四川成都人，副教授，工学博士，从事桥梁结构行为分析研究(E-mail: zhoulingyuan@https://www.360docs.net/doc/8317464806.html,)； 李 乔(1954―)，男，黑龙江铁力人，教授，工学博士，博导，西南交通大学土木工程学院院长，从事桥梁结构行为分析研究 文章编号：1000-4750(2011)01-0082-05 钢筋混凝土梁非线性有限元分析方法 * 周凌远，李 乔 (西南交通大学土木工程学院，成都 610031) 摘 要：针对钢筋混凝土结构有限元分析中，材料进入非线性阶段后，难以通过梁理论准确描述混凝土截面和钢筋应力状态的问题，提出了基于柔度法和分布式塑性理论的钢筋混凝土梁单元材料非线性方法——网格截面法。这种方法采用平面等参单元将梁单元网格化，由单元轴向积分点位置截面网格积分点的混凝土应力描述单元截面应力分布，同时考虑钢筋对刚度的贡献，并通过对截面网格材料的积分计算积分点位置的截面刚度矩阵，再利用力插值函数和能量原理得到梁单元的柔度矩阵，进而对柔度矩阵求逆计算单元刚度矩阵。通过算例验证该方法在钢筋混凝土承载力分析时的准确性。 关键词：有限元；钢筋混凝土梁；柔度法；网格截面；极限承载力 中图分类号：TU375.1; O241.82 文献标识码：A AN APPROACH OF NONLINEAR FINITE ELEMENT ANALYSIS OF REINFORCED CONCRETE BEAM * ZHOU Ling-yuan , LI Qiao (School of Civil Eng, Southwest Jiaotong University, Chengdu 610031, China) Abstract: A beam element with a meshed section based on distributed plasticity and flexibility theory is presented for the material nonlinear finite element analysis of a reinforced-concrete framed structure, the sections of a concrete beam element are discretized into the plane isotropic components in this formulation, the stress distribution on the sections is described with the stresses at quadrature points in the mesh, the stiffness matrices of the sections are calculated by integration of the stress-strain relations of the material on the meshes and the contribution of the stiffness by reinforcing steel is also counted, the flexibility matrix of the element is formed by integration of section flexibility matrices with force-interpolation functions, and then it is inverted to obtain the element stiffness matrix. Finally, a numerical example of the ultimate load capacity analysis of a reinforced concrete beam illustrates the accuracy of the formulation. Key words: finite element; reinforced concrete beam; flexibility method; meshed section; load capacity 钢筋混凝土结构的整体承载力问题一直为工程界所关注，材料非线性有限元方法是研究这类问题的有效手段，其分析模型主要包括集中塑性铰 法[1]和纤维模型法，1977年，Kang 提出了基于纤维模型的二维梁单元[2]，并运用于预应力混凝土框 架的分析，1993年Izzuddin B A 等提出了三次多项式插值的分布式塑性方法分析空间梁单元[3 ―4] ，通 过对沿梁轴方向两个积分点位置的截面划分监控区域，并假定每个监控区域内的法向应力均匀，得到单元的刚度矩阵和节点力，这样在同一个单元内
有限元方法理论及其应用
1 课程论文：弹性力学有限元位移法原理（30分） 撰写一篇论文，对有限元位移法的原理作一般性概括和论述。要求论文论及但不限于下列内容：1）弹性力学有限元位移法的基本思想和数学、力学基础；2）有限元法求解的原理和过程，推导计算列式；对基本概念和矩阵符号进行解释和讨论；3）等参单元的概念、原理和应用。 1.1 对一维杆单元有限元形式的理解 我对此提出了几点疑问： 1)为什么边界条件u1=0，就要划去刚度矩阵[K]中对应的行列再解方程？ 2)为什么刚度矩阵[K]会奇异？ 3)为什么平衡方程本身是矛盾的，而加上边界条件u1=0之后就能解出一 个唯一的近似解？ 4)为什么刚度矩阵[K]是对称的？ 下面我谈谈自己的理解:节点平衡方程是在u1不定的前提下，假设单元内位移都是线性变化推导出来的，由此u1相当于一个不确定的定值约束，再加上中间两个节点的连续性要求，系统实际上只有三个独立的自由度（广义坐标）。 对于第一个问题，其实刚度矩阵[K]中的元素不是一成不变的，相反它是伴随边界条件动态变化的。当u1=0时由刚度矩阵的推导过程可以知道，刚度矩阵的第一行和第一列都会变为0，所以此时第一行和第一列对于求解方程是没有作用的。 对于第二个问题，由于系统自由度（广义坐标）只有三个，而我们的方程却列出
了四个，显然
这四个方程不可能线性无关，所以刚度矩阵奇异。 对于第三个问题，首先我们应该明确方程区别于等式，虽然左右两边都是用“=”连接，但是方程只在特殊条件下取得定解。由于平衡方程是在没有约束的条件下推导出来的，显然它不可能满足等式要求。宏观上看，系统在没有外部约束，而又施加有外力，显然系统会产生加速度而绝不会平衡。所以平衡方程本身是矛盾的。而加上边界条件之后，不但满足了平衡的前提，还改变了矩阵的结构和性质，所以有解。但是，由于我们提前假设了位移线性变化，相当于人为对单元施加了额外约束，让位移按照我们假设的规律变化，所以得到的解是过刚的近似解。但对于方程本身而言是精确解。 对于第四个问题，其力学的作用机理类似于作用力与反作用力，由于刚度矩阵不表征方向，所以其大小是相等的。 1.2 有限元法的思想 有限元法是求解连续介质力学问题的数值方法，更一般意义是一种分析结构问题和连续场数学物理问题的数值方法。 有限元法的基本思想是离散化和分片插值。 即把连续的几何机构离散成有限个单元，并在每一个单元中设定有限个节点，从而将连续体看作仅在节点处相连接的一组单元的集合体，同时选定场函数的节点值作为基本未知量并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律，再建立用于求解节点未知量的有限元方程组，从而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题。 求解得到节点值后就可以通过设定的插值函数确定单元上以至个集合体上的场函数。对每个单元，选取适当的插值函数，使得该函数在子域内部、在子域分界面上以及子域与外界面上都满足一定的条件。单元组合体在已知外载荷作用下处于平衡状态时，列出一系列以节点、位移为未知量的线性方程组，利用计算机解出节点位移后，再用弹性力学的有关公式，计算出各单元的应力、应变，当各单元小到一定程度，那么它就代表连续体各处的真实情况。
浅析有限元方法的发展与应用
浅析有限元方法的发展与应用 摘要:1965年“有限元”这个名词第一次在我国出现,到今天有限元在工程上得到广泛应用,经历了三十多年的发展历史,理论和算法都已经日趋完善。有限元法(Finite Element Method,简写为FEM)是求解微分方程的一种非常有效的数值计算方法,用这种方法进行波动数值模拟受到越来越多的重视。 关键字:有限元法发展应用 Abstract:1965 the term "finite element" first appeared in our country, to this day the finite ele ment is widely used in engineering, has experienced more than 30 years of development history, the ory and algorithm have been becoming more complete.FEM (Finite Element Method, abbreviated a s FEM) is a very effective to solve the differential equation of numerical calculation Method of wav e numerical simulation by using this Method is more and more attention. Keywords: finite element method development Application 绪论 有限元法是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。它是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。 一、有限元的发展历程 有限元法的发展历程可以分为提出(1943)、发展(1944-1960)和后期(1961-二十世纪九十年代)三个阶段。有限元法是受内外动力的综合作用而产生的。 1943年，柯朗在《美国数学学会公报》(Bulletin of The American Mathematical Society)上发表了《平衡和振动问题的变分解法》 (Variational Methods for The Solution of Problems of Equilibrium And Vibration)一文，这篇文章实际上是他1941年在美国数学学
有限元分析的基本步骤
一个典型的ANSYS分析过程可分为以下6个步骤： 1定义参数 2创建几何模型 3划分网格 4加载数据 5求解 6结果分析 1定义参数 1.1指定工程名和分析标题 启动ANSYS软件，选择File→Change Jobname命令 选择File→Change Title菜单命令 1.2定义单位 (2) 设置计算类型 ANSYS Main Menu: Preference→Material Props →Material Models →Structural →OK (3) 定义分析类型 ANSYS Main Menu: Preprocessor →Loads →Analysis Type →New Analysis→STATIC →OK 1.3定义单元类型 选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令 单击[Options]按钮，在[Element behavior]下拉列表中选择[Plane strs w/thk]选项，单击确定 1.4定义单元常数 在ANSYS程序主界面中选择Main Menu→Preprocessor→Real Constants→Add/Edit/Delete命令 单击[Add]按钮，进行下一个[Choose Element Type]对话框 1.5定义材料参数 在ANSYS程序主界面，选择Main Menu→Preprocessor→Material Props→Material Models命令 (1)选择对话框右侧Structural→Linear→Elastic→Isotropic命令，并单击[Isotropic]选项，接着弹出如下所示[Linear Isotropic Properties for Material Number 1]对话框。 在[EX]文本框中输入弹性模量“200000”，在[PRXY]文本框中输入泊松比“0.3”，单击OK 2创建几何模型 在ANSYS程序主界面，选择Main Menu→Preprocessor→Modeling→Creat→Areas→Rectangle →By 2Corners命令 选择Main Menu→Preprocessor→Modeling→Creat→Areas→Circle→Solid Circle命令 3网格划分(之前一定要进行材料的定义和分配) 选择Main Menu→Preprocessor→Modeling→Operate→Booleans→Subtract→Arears Circle命令 选择Main Menu→Preprocessor→Meshing→Mesh→Areas→Free命令，弹出实体选择对话框，单击[Pick All]按钮，得到如下所示网格 4加载数据 (1)选择Main Menu→Preprocessor→Loads→Define Loads→Apply→Structural→Displacement→On Lines命令， 出现如下所示对话框，选择约束[ALL DOF]选项，并设置[Displacement value]为0，单击OK。
非线性有限元分析
非线性有限元分析 1 概述 在科学技术领域内，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。 已经发展的数值分析方法可以分为两大类。一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。 另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。 1960年，发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。有限单元法的基本思想是将连续的求解区域离散为一组有限个，且按一定方式相互联接在一起的单元的组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模型化几何形状复杂的求解域。并且可以利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数，从而使一个连续的无限自由度问题变成离散的有限自由度问题。 现已证明，有限单元法是基于变分原理的里兹法的另一种形式，从而使里兹法分析的所有理论基础都适用于有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法。利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，而且事先不要求满足任何边界条件，因此可以用来处理很复杂的连续介质问题。 在短短四十余年的时间里，有限单元的分析方法已经迅速地发展为适合于使用各种类型计算机解决复杂工程问题的一种相当普及的方法。如今，有限元广泛地应用于各个学科门类，已经成为工程师和科研人员用于解决实际工程问题，进行科学研究不可或缺的有力工具。有限单元法的应用范围已由弹性力学平面问题扩展到空间问题，板壳问题，由静力平衡问题扩展到稳定问题，动力问题和波动问题。分析的对象从弹性材料扩展到塑性，粘弹性，粘塑性和复合材料等，从固体力学扩展到流体力学，传热学等连续介质力学领域。在工程分析中的作用已从分析和校核扩展
国内外主要有限元分析软件比较
有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC四个比较知名比较大的公司。 常见软件 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC四个比较知名比较大的公司，其中ADINA、ABAQUS在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。 软件对比 ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：1、ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：1、ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：1、ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS软件与ANSYS软件的对比分析 1．在世界范围内的知名度 两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。 由于ANSYS产品进入中国市场早于ABAQUS，并且在五年前ANSYS的界面是当时最好的界面之一，所以在中国，ANSYS软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北
浅析有限元方法的发展与应用
浅析有限元方法的发展与应用 1965年“有限元”这个名词第一次在我国出现,到今天有限元在工程上得到广泛应用,经历了三十多年的发展历史,理论和算法都已经日趋完善。有限元法(Finite Element Method,简写为FEM)是求解微分方程的一种非常有效的数值计算方法,用这种方法进行波动数值模拟受到越来越多的重视。 绪论 有限元法是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。它是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。 一、有限元的发展历程 有限元法的发展历程可以分为提出(1943)、发展(1944-1960)和后期(1961-二十世纪九十年代)三个阶段。有限元法是受内外动力的综合作用而产生的。 1943年，柯朗在《美国数学学会公报》(Bulletin of The American Mathematical Society)上发表了《平衡和振动问题的变分解法》(Variational Methods for The Solution of Problems of Equilibrium And Vibration)一文，这篇文章实际上是他1941年在美国数学学会演讲的书面稿，在其中柯朗提出了有限元法的核心思想。大约与柯朗同时，工程师阿格瑞斯在另一个领域独立地提出了有限元法。柯朗和阿格瑞斯各自在数学和工程学领域独立提出了有限元法，他们分别开创了有限元法的数学传统和工程学传统。 有限元法被提出来以后，经过一段时间的沉寂期，在二十世纪五十年代和六十年代初有了很大的发展。主要表现为在代数表达形式、单元划分、单元类型选择和解的收敛性研究上取得的突破。1960年，克劳夫在《平面应力分析中的有限元》(The Finite Element in Plan Stress Analysis)的论文中，第一次从数学上说明了将定义域划分成有限的单元能够成功的原因:他表明对一些特定类型的单元来说，随着单元尺寸的减小，近似解将收敛到精确解，这就在某些情况下证明了有限元法的收敛性。并第一次提出了“有限元法”这个名称，这个名称一直沿用至今，标志着有限元法早期发展阶段的结束。 有限元法后期阶段的发展有国外和国内两条线索。在国外的发展表现为: 第一，建立了严格的数学和工程学基础；第二，应用范围扩展到了结构力学以外的领域；第三，收敛性得到了进一步研究，形成了系统的误差估计理论；第四，发展起了相应的商业软件包。 在国内，我国数学家冯康在特定的环境中独立于西方提出了有限元法。1965年，他发表论文《基于变分原理的差分格式》，标志着有限元法在我国的诞生。冯康的这篇文章不但提出了有限元法，而且初步发展了有限元法。他得出了有限元法在特定条件下的表达式，独创了“冯氏大定理”并且初步证明了有限元法解的收敛性。虽然冯康创造的有限元法不成熟，但他能在当时的条件下独立提出有限元法已十分不易。对于他的这项成就，国内外专家学者和国家领导人都有很高的评价。 半个世纪来，有限单元法蓬勃发展，不仅已经成为结构分析中必不可少的工具，而且成为现象分析的一种手段。其应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、渗流与固结理论、热传导与热应力问题、磁场问题以及建筑声学与噪音问题。不仅涉及稳态场问题，还涵盖材料非线性、几何非线性、时间维问题和断裂力学等。
有限单元法作业非线性分析+程序
几何非线性大作业荷载增量法 和弧长法程序设计 系（所）：建筑工程系 学号：1432055 姓名：焦联洪 培养层次：专业硕士 指导老师：吴明儿 2015年6月19日
一、几何非线性大作业（ Newton-Raphson法） 用荷载增量法（Newton-Raphson法）编写几何非线性程序： （1）用平面梁单元，可分析平面杆系 （2）算例：悬臂端作用弯矩。悬臂梁最终变形形成周长为悬臂梁长度的圆。 1.1 Newton-Raphson算法基本思想 图1.1 Newton-Raphson算法基本思想 1.2 悬臂梁参数 基本参数：L=2m, D=0.03m, A=7.069E-4m2, I=3.976E-08m4 ,E=2.0E11N/m2
图1.2 悬臂梁单元信息 将悬臂梁分成10个单元，如图1.2所示 2.1 MATLAB输入信息 材料信息单元信息 约束信息（0为约束，1为放松）荷载信息（FX,FY,M）
节点信息 2.2 求解过程 梁弯成圆形：理论弯矩M=EIY"=24981.944N.m ,直径为0.642m 运用ABAQUS和MATLAB进行求解对比： 图1.3 加载图 图1.4 ABAQUS变形图
图1.5 MATLAB变形曲线 ABAQUS和MATLAB变形对比，最终在理论荷载作用下都弯成了一个圆，其直径为0.64716m，与理论值相对比值为：（0.64716-0.642）/0.642=0.00804.非常接近。 2.3 加载点荷载位移曲线 图1.5 加载点Y方向的荷载位移曲线
加载点的最大竖向位移分别为1.4525m和1.45246m，相对比值（1.4525-1.45246）/1.45246=2.75395E-05。完全相同，说明MATLAB的计算结果很好。
有限元分析方法在工程中的应用
有限元分析方法在工程中的应用 Application of finite element analysis method in Engineering 一、引言 从20世纪50年代诞生到现在，有限元方法和技术经历了60年的发展历程，已经成为当今科学与工程领域中分析和求解微分方程的系统化数值计算方法。由于有限元分析方法适用性强、形式简单、理论可靠等众多优点，近年来已被推广应用到航空航天、土木建筑、机械等相关科学领域。本文以ANSYS软件为例，介绍其功能和应用，包括几何建模技术、网格划分与有限元建模技术、施加载荷与求解过程、结果后处理技术等。图1是用有限元方法分析工程问题时的具体步骤[1]。 本文以车轮钢的疲劳性能研究为例，介绍有限元分析方法在其中的应用。 图1. 有限元方法进行计算机辅助工程分析的步骤 二、ANSYS操作步骤 ANSYS的基本操作步骤包括建模、划分网格、加载求解和后处理等步骤。进入ANSYS系统后有六个系统，提供使用者和软件之间的交流凭借这六个窗口可以实现输入命令、检查模型的建立、观察分析结果及图形输出与打印。ANSYS
各窗口及工具条如图2所示。 图2. ANSYS的窗口及工具条 1、建立模型 首先必须指定作业名和分析标题，接着使用PREP7前处理器定义单元类型、单元实常数、材料特性，然后建立几何模型。需要注意的是，ANSYS的GUI界面下没有类似WORD中的后退操作按钮，所以就出现了一个常见问题：做错一步操作如何后退？这里可以采用三种方法：（1）建模阶段可以使用Delete(删除)图元命令，划分网格阶段可以使用Clear（清除）单元命令。（2）每完成一个模块的操作，都用SA VE AS保存数据到不同名的数据库文件中，出错后点击Resum Form恢复。（3）使用命令：UNDO，ON以便激活ANSYS内部的返回命令。 本文以车轮钢为例，建立好的模型与图2类似，只是未划分网格。 2、单元网格划分 一个实体模型进行网格划分（meshing）之前必须指定所产生的单元属性（element attribute）。ANSYS有限元网格划分是进行数值模拟分析至关重要的一步，它直接影响着后续数值计算分析结果的精确性。ANSYS软件平台提供了映射网格划分和自由网格划分的策略。映射网格划分用于曲线、曲面、实体的网格划分方法，自由网格划分方法用于空间自由曲面和复杂实体。