第四章 有限元法

第四章有限元法第1页，共92页，2023年，2月20日，星期四主要内容有限元法的基本概念有限元法的分类有限元法的求解步骤（重点）常用有限元软件简介第2页，共92页，2023年，2月20日，星期四随着生产的发展，不断要求设计高质量、高水平的大型、复杂和精密的机械和工程结构。在产品加工中，在分析产品性能过程中，往往需要了解产品在承受载荷的情况下的应变与位移、温度场、流场、磁场的分布情况等，从而为产品的性能设计提供依据。

第3页，共92页，2023年，2月20日，星期四蓄水后大坝的位移与应变情况、地震时大坝的位移与应变情况等三峡大坝的受力情况

第4页，共92页，2023年，2月20日，星期四航天飞机飞行中的受热分析温度场分布

第5页，共92页，2023年，2月20日，星期四导弹、飞机飞行的流体动力学分析流场分布

第6页，共92页，2023年，2月20日，星期四磁场分布

分析卫星、飞船在轨运行时磁场的影响

第7页，共92页，2023年，2月20日，星期四传统方法在处理这类问题时，往往要对一个实际的物理系统作出多种假设，比如形状假设、连续性假设、物体的各项同性假设，然后通过经典理论方法得出问题的解析解，这种解析解从形式上看，可以得出关于实际问题的连续解，比如用方程描述三峡大坝每一点的位移和应变，但这样的解析解往往和实际情况有比较大的偏差。这对于精度要求不高的领域是可以的，但对于有些领域，就不能满足实际的需要了。

第8页，共92页，2023年，2月20日，星期四同时，实际中常常要遇到一些几何上复杂、不规则边界、有裂缝或厚度突变以及几何非线性、材料非线性的物理系统，对这些系统经典理论解决起来相当困难，有时甚至无法解决，也就是无法求得解析解。因此，寻求离散数值分析法就成了必由之路。常用的数值分析法有两种：差分法和有限元法。

第9页，共92页，2023年，2月20日，星期四

差分法是在传统方法的基础上，将传统方法建立的微分方程中的微分dx、dy、dz变成差分Δx，Δy，Δz，从而把微分方程变成代数方程，用一步步迭代的方法，逐步求出物理系统中各个离散点的物理量，用差分离散解代替连续解。这种方法要求能建立微分方程，并能给出边界条件的数学表达式，因此，对于一些不规则的几何形状和不规则的特殊边界条件难以应用。第10页，共92页，2023年，2月20日，星期四一、有限元法的基本概念1.什么是有限元法我们实际要处理的对象都是连续体，在传统设计思维和方法中，是通过一些理想化的假定后，建立一组偏微分方程及其相应的边界条件，从而求出在连续体上任一点上未知量的值。因为点是无限多的，存在无限自由度的问题，很难直接求解这种偏微分方程用来解决实际工程问题，因此需要采用近似方法来处理。

第11页，共92页，2023年，2月20日，星期四其中最主要的是离散化方法，把问题归结为只求有限个离散点的数值，把无限自由度问题变成有限个自由度。

把一个连续体分割成有限个单元，即把一个复杂的结构看成由有限个通过节点相连的单元组成的整体，先进行单元分析，然后再把这些单元组合起来代表原来的结构，以得到复杂问题的近似数值解。这种方法称为有限元法（TheFiniteElementMethod）。

第12页，共92页，2023年，2月20日，星期四有限元法是一种以计算机为手段，通过离散化将研究对象变换成一个与原始结构近似的数学模型，再经过一系列规范化的步骤以求解应力、应变、位移等参数的数值计算方法。所谓离散化就是将一个连续体分割成若干个通过节点相连的单元，这样一个有无限个自由度的结构就变换成一个具有有限个自由度的近似结构。该过程还包括对单元和节点进行编码以及局部坐标系和整体坐标系的确定。第13页，共92页，2023年，2月20日，星期四2.几个基本概念1）单元（element）将求解的工程结构看成是由许多小的、彼此用点联结的基本构件如杆、梁、板和壳组成的，这些基本构件称为单元。在有限元法中，单元用一组节点间相互作用的数值和矩阵（刚度系数矩阵）来描述。第14页，共92页，2023年，2月20日，星期四单元具有以下特征：每一个单元都有确定的方程来描述在一定载荷下的响应；模型中所有单元响应的“和”给出了设计的总体响应；单元中未知量的个数是有限的，因此称为“有限单元”。第15页，共92页，2023年，2月20日，星期四2）节点（node）单元与单元之间的联结点，称为节点。在有限元法中，节点就是空间中的坐标位置，它具有物理特性，且存在相互物理作用。载荷载荷节点:

空间中的坐标位置，具有一定响应，相互之间存在物理作用。单元:节点间相互作用的媒介，用一组节点相互作用的数值矩阵描述（称为刚度或系数矩阵)。第16页，共92页，2023年，2月20日，星期四.1node..信息是通过单元之间的公共节点传递的。分离但节点重叠的单元A和B之间没有信息传递具有公共节点的单元之间存在信息传递

..AB......AB...2nodes第17页，共92页，2023年，2月20日，星期四3）有限元模型（node）

有限元模型真实系统理想化的数学抽象。由一些形状简单的单元组成，单元之间通过节点连接，并承受一定载荷。

每个单元的特性是通过一些线性方程式来描述的。作为一个整体，所有单元的组合就形成了整体结构的数学模型。

有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。第18页，共92页，2023年，2月20日，星期四膨化饲料床层进风口出风口膨化饲料带式干燥机有限元模型XZ网格节点步长为40mm,共1113040个单元第19页，共92页，2023年，2月20日，星期四对于一个具体的工程结构，单元的划分越小，求解的结果就越精确，同时，其计算工作量也就越大。梯子的有限元模型不到100个方程；在ANSYS分析中，一个小的有限元模型可能有几千个未知量，涉及到的单元刚度系数几百万个。单元划分的精细程度，取决于工程实际对计算结果精确性的要求。第20页，共92页，2023年，2月20日，星期四4）单元形函数（node）

有限元法仅仅求解节点处的响应值。单元形函数是一种数学函数，规定了从节点响应值到单元内所有点处响应值的计算方法,因此，单元形函数提供一种描述单元内部结果的“形状”。单元形函数描述的是给定单元的一种假定的特性。单元形函数与真实工作特性吻合好坏程度直接影响求解精度。第21页，共92页，2023年，2月20日，星期四真实的二次曲线.节点单元二次曲线的线性近似

(不理想结果).2节点单元响应值二次分布..1节点单元线性近似(更理想的结果)真实的二次曲线.....3节点单元二次近似(接近于真实的二次近似拟合)(最理想结果)..4第22页，共92页，2023年，2月20日，星期四如果单元形函数不能精确描述单元内部的响应，就不能很好地得到导出数据，因为这些导出数据是通过单元形函数推导出来的。当选择了某种单元类型时，也就十分确定地选择并接受该种单元类型所假定的单元形函数。在选定单元类型并随之确定了形函数的情况下，必须确保分析时有足够数量的单元和节点来精确描述所要求解的问题。第23页，共92页，2023年，2月20日，星期四5）有限元分析有限元分析（FEA）就是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素（即单元），用有限数量的未知量去逼近无限未知量的真实系统。有限元分析是一种模拟设计载荷条件，并且确定在载荷条件下的设计响应的方法。它是用被称之为“单元”的离散的块体来模拟设计的。第24页，共92页，2023年，2月20日，星期四二、有限元法的特点与作用1.有限元法的特点1）离散化：把连续体划分成有限个单元，把单元间的连接点（节点）作为离散点；2）不考虑微分方程，而从单元本身特点进行研究；(研究未知量在单元内部及在单元节点上值的关系，从而导出单元节点响应和对应的载荷之间的关系，然后把它们组集起来,以求解一个以各节点响应为未知量的代数方程组)

3）理论基础简明，物理概念清晰，且可在不同的水平上建立起对该法的理解；

第25页，共92页，2023年，2月20日，星期四4）具有灵活性和适用性，适应性强。它可以把形状不同、性质不同的单元组集起来求解，故特别适用于求解由不同构件组合的结构，应用范围极为广泛。它不仅能成功地处理如应力分析中的非均匀材料、各向异性材料、非线性应力应变以及复杂的边界条件等问题，且随着其理论基础和方法的逐步完善，还能成功地用来求解如热传导、流体力学及电磁场领域的许多问题。5）在具体推导运算过程中，广泛采用了矩阵方法。

第26页，共92页，2023年，2月20日，星期四2.有限元法的作用1）减少模型试验的数量（计算机模拟允许对大量的假设情况进行快速而有效的试验）；2）模拟不适合在原型上试验的设计（例如：器官移植、人造膝盖）；3）节省费用，降低设计与制造、开发的成本；4）节省时间，缩短产品开发时间和周期；

5）创造出高可靠性、高品质的产品。

第27页，共92页，2023年，2月20日，星期四三、有限元法的发展1.有限元法的产生有限元法分析的概念可以追溯到20世纪40年代。

1943年，柯朗特（Courant）第一次在他的论文中，取定义在三角形域上的分片连续函数，利用最小势能原理研究了圣维南（St.Venant）的扭转问题。然而，此方法发展很慢，几乎过了十年才再次有人用这些离散化的概念。

第28页，共92页，2023年，2月20日，星期四1956年Turner,Clough,Martin和Topp等人，在他们的经典论文中第一次给出了用三角形单元求得的平面应力问题的真正解答，他们利用弹性理论的方程求出了三角形单元的特性，并第一次介绍了今天人们熟知的确定单元特性的直接刚度法，其研究工作随同当时出现的数值计算机一起打开了求解复杂平面弹性问题的新局面。第29页，共92页，2023年，2月20日，星期四1960年美国的克劳夫(W.Clough)采用此方法进行飞机结构分析时首次将这种方法起名为“有限单元法”，简称“有限元法”。此后有限元法在工程界获得了广泛的应用。到20世纪70年代以后，随着计算机和软件技术的发展，有限元法也随之迅速的发展起来，发表的论文犹如雨后春笋，学术交流频繁，期刊、专著不断出现，可以说进入了有限元法的鼎盛时期，对有限元法进行了全面而深入的研究。第30页，共92页，2023年，2月20日，星期四2.有限元法的应用目前，有限元法广泛应用于固体力学、流体力学、热传导、电磁学、声学、生物力学等各个领域。有限元法已经成为结构分析的有效方法和手段，它几乎被用于所有连续介质和场的问题。第31页，共92页，2023年，2月20日，星期四1）结构分析结构分析用于确定变形、应变、应力及反作用力。静力分析用于静态荷载，可以考虑结构的线性及非线性行为，例如：大变形、大应变、应力刚化、接触、塑性、超弹及蠕变等。超弹密封第32页，共92页，2023年，2月20日，星期四动力分析模态分析：用于计算固有频率和振型。谐响应分析：用于确定结构对正弦变化的已知幅值和频率载荷的响应。瞬态动力学分析：用于确定结构对随时间任意变化载荷的响应，可以考虑与静力分析相同的结构非线性行为。特征屈曲分析：用于计算线性屈曲载荷并确定屈曲模态形状（结合瞬态动力学分析可以实现非线性屈曲分析）专项分析:

断裂分析、复合材料分析、疲劳分析第33页，共92页，2023年，2月20日，星期四2）热分析热分析用于确定物体中的温度分布。考虑的物理量是：热量、热梯度、热通量。有限元可模拟三种热传递方式（热传导、热对流、热辐射），可进行稳态分析和瞬态分析，还可模拟相变（蒸发与冷凝、熔化及凝固）第34页，共92页，2023年，2月20日，星期四3）电磁分析电磁分析用于计算电磁装置中的磁场，考虑的物理量是：磁通量密度、磁场密度、磁力和磁力矩、阻抗、电感、涡流、能耗及磁通量泄漏等。有限元法可进行静态磁场及低频电磁场分析，如模拟由直流电源、低频交流电或低频瞬时信号引起的磁场（螺线管制动器、电动机、变压器）；高频电磁场分析如模拟电磁波的传播装置；计算由电压或电荷激发引起的电场；分析电磁装置与电路的耦合问题等。第35页，共92页，2023年，2月20日，星期四4）流体分析计算流体动力学（ComputationalFluidDynamics

，CFD

）是一种典型的基于有限元的流体分析方法。用于确定流体中的流动状态和温度。能模拟层流和湍流、可压缩和不可压缩流体以及多组份流。广泛应用于航空航天、气象学、汽车设计等领域。典型的物理量是：速度、压力、温度、对流换热系数。第36页，共92页，2023年，2月20日，星期四用于模拟流体介质和周围固体的相互作用。典型的物理量是：压力分布、位移和自振频率。5）声学分析第37页，共92页，2023年，2月20日，星期四6）耦合场分析耦合场分析考虑两个或多个物理场之间的相互作用。因为两个物理场之间相互影响，所以单独求解一个物理场是不可能的。例如：双金属片受热变形热-应力分析（温度场和结构）流体热力学分析（温度场和流场）声学分析（流体和结构）热-电分析（温度场与电场）感应加热（磁场和温度场）第38页，共92页，2023年，2月20日，星期四第二节有限元法的分类第39页，共92页，2023年，2月20日，星期四一、结构有限元法的分类

结构有限元法可以分为两类，即线弹性有限元法和非线性有限元法。其中线弹性有限元法是非线性有限元法的基础，二者不但在分析方法和研究步骤上有类似之处，而且后者常常要引用前者的某些结果。第40页，共92页，2023年，2月20日，星期四1.线弹性有限元线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；位移与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。第41页，共92页，2023年，2月20日，星期四2.非线性有限元非线性问题与线弹性问题的区别：非线性问题的方程是非线性的，一般需要迭代求解；非线性问题不能采用叠加原理；非线性问题不总有一致解，有时甚至没有解。以上三方面的因素使得非线性问题的求解过程比线弹性问题更加复杂、费用更高和更具有不可预知性。第42页，共92页，2023年，2月20日，星期四1）材料非线性问题材料的应力和应变是非线性的，但应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的应力与应变的函数关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。第43页，共92页，2023年，2月20日，星期四2）几何非线性问题几何非线性问题是由于应变与位移之间存在非线性关系引起的。当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。第44页，共92页，2023年，2月20日，星期四3）非线性边界在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。实际的非线性可能同时出现上述两种或三种非线性问题。第45页，共92页，2023年，2月20日，星期四二、有限单元的类型根据研究对象的不同，有限元法中采用的单元形式也不相同。

通常，按照单元结构，可将单元划分维一维单元（线单元）、二维单元（面单元）和三维单元JIJKLI一维单元二维单元POMNKJIL三维单元第46页，共92页，2023年，2月20日，星期四按照单元结构特点和受力特点，可将单元划分为：1）桁架杆单元：主要应用于受轴向力作用的杆和杆系，如桁架结构；2）刚架杆单元：用于梁及刚架结构分析；3）三角形平面单元：主要用于弹性力学中平面应力和平面应变问题的有限元分析；4）三棱圆环单元：用于轴对称问题的有限元分析；5）等参数单元：用于一些具有曲线轮廓的复杂结构。第47页，共92页，2023年，2月20日，星期四第三节有限元法的求解步骤第48页，共92页，2023年，2月20日，星期四一、有限元法分析过程

有限元法分析过程大体可分为：前处理、分析、后处理三大步骤。1.前处理

对实际的连续体离散化后就建立了有限元分析模型，这一过程是有限元法的前处理过程。在这一阶段，要构造计算对象的几何模型，要划分有限元网格，要生成有限元分析的输入数据，这一步是有限元分析的关键。

第49页，共92页，2023年，2月20日，星期四2.有限元分析

有限元分析过程主要包括：单元分析、整体分析、载荷移置、引入约束、求解约束方程等过程。这一过程是有限元分析的核心部分，有限元理论主要体现在这一过程中。

结构有限元法包括三类：有限元位移法、有限元力法、有限元混合法。第50页，共92页，2023年，2月20日，星期四

有限元位移法以节点位移作为基本未知量；有限元力法以选节点力作为未知量；有限元混合法的一部分基本未知量为节点位移，另一部分基本未知量为节点力。有限元位移法计算过程的系统性、规律性强，特别适宜于编程求解。一般除板壳问题的有限元应用一定量的混合法外，其余全部采用有限元位移法。

因此，一般不做特别声明时，有限元法指的是有限元位移法。

第51页，共92页，2023年，2月20日，星期四3.有限元分析的后处理有限元分析的后处理主要包括对计算结果的加工处理、编辑组织和图形表示三个方面。它可以把有限元分析得到的数据，进一步转换为设计人员直接需要的信息，如应力分布状态、结构变形状态等，并且绘成直观的图形，从而帮助设计人员迅速的评价和校核设计方案。第52页，共92页，2023年，2月20日，星期四一、有限元法求解的基本步骤1.结构离散化：对整个结构进行离散化，将其分割成若干个单元，单元间彼此通过节点相连；2.求出各单元的刚度矩阵：是由单元节点位移向量求单元节点力向量的转移矩阵，其关系式为；

3.集成总体刚度矩阵[K]并写出总体平衡方程：[K]是由整体节点位移向量求整体节点力向量的转移矩阵，其关系式为，这就是总体平衡方程；第53页，共92页，2023年，2月20日，星期四确定总体刚度矩阵[K]的办法1）直接利用总体刚度系数的定义

在求出整体结构中各节点力与节点位移关系的基础上获得总体刚度矩阵。此方法只在简单情况下才能采用。

2）集成法将整体坐标下的单元刚度矩阵进行迭加而得。这里所说的迭加不是简单的相加，而是将下角标相同的刚度系数相加，然后按总码的顺序对号入座。第54页，共92页，2023年，2月20日，星期四3）利用节点间的刚度系数直接写出总体刚度矩阵总体刚度矩阵对角线上的刚度系数等于在节点i汇交的几个单元的刚度系数之和；非对角线上的刚度系数等于联结节点i与节点j间几个单元的刚度系数之和。第55页，共92页，2023年，2月20日，星期四4.引入边界条件，求出各节点的位移节点的边界条件有两种：一种是节点n沿某个方向的位移为零，另一种是节点n沿某个方向的位移为一给定值。5.求出各单元内的应力和应变

第56页，共92页，2023年，2月20日，星期四第57页，共92页，2023年，2月20日，星期四三、有限元求解实例分析

【例1】一根由两段组成的阶梯轴，一端固定，另一端承受一个轴向载荷F3。这两段的横截面积分别为A(1)和A(2)，长度分别为L(1)和L(2)，弹性模量分别为E(1)和E(2)

，如图所示。求出这两段的应力和应变。已知数据分别为F3=100N，

第58页，共92页，2023年，2月20日，星期四1①2②3A(1)E(2)A(2)E(2)L(1)

L(2)

①②2Φ1F2F3Φ2Φ3F1F3第59页，共92页，2023年，2月20日，星期四【解】1.离散化

把这根阶梯轴看成是由两个单元组成的，节点选在截面积突变处，两个单元的连接处是一个节点，该阶梯轴的两端视为另外两个节点，所以整个结构共有三个节点。这根轴是一维结构，并只受轴向载荷，因此各单元内只有轴向位移。三个节点位置的位移量分别记为、、。在整个结构中节点载荷及节点位移均用大写字母标记，其角标为节点在总体结构中的编码，简称总码。

第60页，共92页，2023年，2月20日，星期四2.求单元刚度矩阵

该矩阵中任意一个元素都称为单元刚度系数，它表示：该单元内除节点j产生单位位移外，其余各节点的位移均为零时在节点i

处所引起的载荷。第61页，共92页，2023年，2月20日，星期四3.总体刚度矩阵的集成和总体平衡方程的写出该阶梯轴上三个节点位移、、和三个节点轴向力分别组成该整体结构的节点位移向量和节点轴向力向量。

两向量间的转换关系可表示为或

第62页，共92页，2023年，2月20日，星期四

上式中的转移矩阵称为总体刚度矩阵或总体特性矩阵，其阶数等于总体结构中的节点总数。

[K]中的元素称为总体刚度系数，它表示在整体结构中除了节点j

产生单位位移外，其余各节点的位移均为零时在节点i处所引起的载荷。

第63页，共92页，2023年，2月20日，星期四求出总体刚度矩阵是进行总体分析的主要任务，一旦获得总体刚度矩阵，可以很容易地写出总体平衡方程。求总体刚度矩阵[K]的方法主要由两种：一是直接法，即根据总体刚度系数的定义求解；另一种方法是集成法，即由各单元刚度矩阵求总体刚度矩阵。第64页，共92页，2023年，2月20日，星期四1)用集成法求总体刚度矩阵[K]

这种方法从单元刚度矩阵出发，根据迭加原理，利用刚度系数集成的方法获得总体刚度矩阵。这样，首先要写出各单元的刚度矩阵。

局部码1(1)2(1)1(2)2(2)总码1223局部码与总码对应关系第65页，共92页，2023年，2月20日，星期四约定：

第66页，共92页，2023年，2月20日，星期四第67页，共92页，2023年，2月20日，星期四2)用直接法求总体刚度矩阵[K]这种方法具有概念清晰的特点，但是在分析复杂结构时运算极其复杂，因而限制了它的应用。

直接法利用总体刚度系数的定义，直接求出每一个总体刚度系数第68页，共92页，2023年，2月20日，星期四一般而言，刚度矩阵具有如下特性：1）对称性

单元刚度矩阵和总体刚度矩阵都是对称方阵，即，即由第j个节点单位位移引起的第i个节点载荷和由第i个节点单位位移引起的第j个节点载荷是相等的。这是弹性结构一个共同的特点。这种对称性可减少矩阵存储运算时的内存量。

第69页，共92页，2023年，2月20日，星期四2）奇异性

单元刚度矩阵和总体刚度矩阵都是奇异矩阵，即它们的行列式都等于0，这样，其逆阵就不存在。因此，对总体刚度矩阵要引入边界条件进行处理之后才能求解。3）稀疏性

总体刚度矩阵是零元素非常多的矩阵。结构越大，零元素越多。大型结构的总体刚度矩阵一般都是非常稀疏的矩阵。

第70页，共92页，2023年，2月20日，星期四总体平衡方程为：

第71页，共92页，2023年，2月20日，星期四4.引入支撑条件，计算节点位移上式中的未知量仍不能求出，因为[K]是一个奇异矩阵，必须引入支撑条件。在本例中支撑条件是节点1的位移为零，即这样总体平衡方程简化为

第72页，共92页，2023年，2月20日，星期四代入已知条件：

可求得：

第73页，共92页，2023年，2月20日，星期四5.求单元中的应力及应变单元1中的应变：单元2中的应变：单元1中的应力：单元2中的应力：

第74页，共92页，2023年，2月20日，星期四K2【例2】在光滑的水平面上有三个小车，他们彼此用四根弹簧相连，其连接方式如图1-2所示。小车1又通过弹簧k1与墙壁固连。每个小车上的作用力分别为F1、F2、F3。求每个小车的位移。

第75页，共92页，2023年，2月20日，星期四由此可求出三个节点即三个小车的位移。

第76页，共92页，2023年，2月20日，星期四【例3】用有限元法求解某系统的过程中，将系统划分成3个单元，这3个单元通过3个节点相连接，节点局部编码与总体编码的对照如下表所示，请写出该系统的总体刚度矩阵。

局部码1（1）1（2）2（2）1（3）2（3）总码11223已知：

第77页，共92页，2023年，2月20日，星期四第78页，共92页，2023年，2月20日，星期四第四节常用有限元软件简介

第79页，共92页，2023年，2月20日，星期四一、通用有限元软件的共同之处有限元的高度通用性与实用性导致了有限元通用程序的发展。四十多年来，有限元通用软件的发展在数量和规模上是惊人的。这些通用有限元软件的共同之处可归结为以下几点：

1.功能强大。一般都可以进行多种物理场分析，如结构分析、温度场分析、电磁场分析、流场分析、多场耦合分析等；

第80页，共92页，2023年，2月20日，星期四

2.具有丰富的材料库。可以处理多种材料，如金属、土壤、岩石、塑料、橡胶、木材、陶瓷、混凝土、复合材料等；

3.具有多种自动网格划分技术，自动进行单元形态、求解精度检查及修正；

4.具有强大的后处理及图像处理功能；

5.具有与多种CAD系统直接连接的接口；

6.具有良好的用户开发环境；

7.具有良好的维护和培训能力（ＧＵＩ）。第81页，共92页，2023年，2月20日，星期四二、几个著名的通用有限元软件简介ANSYSANSYS软件是美国ANSYS公司的产品，该公司成立于1970年，公司总部位于美国宾夕法尼亚的匹兹堡。ANSYS软件是融结构、流体、电磁场、声场和耦合场分析于一体的大型通用有限元软件，可广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等一般工业及科学研究。第82页，共92页，2023年，2月20日，星期四

ANSYS公司在北京、上海、成都相继成立了办事处，构成了ANSYS在中国完整的市场、销售及售后服务体系。它的显著特点是具有独一无二的多场耦合分析功能，可处理高速变形和高度非线性问题（如冲击、爆炸、碰撞、实体变形、板成形），边界元流体动力学问题（如水下结构振动、气弹颤振分析）。第83页，共92页，2023年，2月20日，星期四2.MSC.MarcMarc软件原为美国MARC公司的产品，该公司创建于1967年，他的创始人是美国著名的布朗大学教授、有限元分析的先驱者PedroMarcel。MARC公司致力于非线性有限元技术的研究、非线性有限元软件的开发、销售和售后服务。经过三十多年的不懈努力，Marc软件得到了学术界和工业界的大力推崇和广泛应用，建立了他在全球非线性有限元软件行业的领导地位。第84页，共92页，2023年，2月20日，星期四

1999年6月，美国MSC公司收购了MARC公司，相应地将该软件更名为MSC.Marc软