理学有限单元法简介

2024/7/5有限单元法简介1课时安排共36

学时有限元法：30学时软件应用：6

学时2024/7/5有限单元法简介2考试考核成绩构成InstituteofMechanicalEngineeringandAutomation听课表现(30分)

有限元软件考题（70分）

2024/7/5有限单元法简介3课程学习基本要求

教材及参考书InstituteofMechanicalEngineeringandAutomation1、李人宪.有限单元法基础(第二版)[M].北京:国防工业出版社,20042、傅永华.有限单元法基础[M].武汉:武汉大学出版社,20033、王生洪等.有限单元法基础及应用[M].北京:国防工业出版社,19904、黄国权.有限单元法基础及ANSYS应用[M].北京:机械工业出版社,20045、王勖成、邵敏.有限单元法[M].北京:清华大学出版社,1997.2024/7/5有限单元法简介4学习目的及方法InstituteofMechanicalEngineeringandAutomation目的

了解FEM数学力学基础；把握FEM求解具体问题的基本过程；应用FEM，特别是运用已有的通用或专用软件求解实际工程技术问题；方法

注重深入理解FEM思想的建立，数学力学基础;

运用通用或专用有限元程序进行工程问题分析；一、数值模拟方法概述二、有限单元法简介三、有限单元法分析步骤四、利用有限元软件进行工程分析五、结后语第1章有限单元法概述2024/7/5有限单元法简介6一、数值模拟方法概述工程技术领域中的许多力学问题和场问题，如固体力学中的位移场、应力场分析；电磁学中的电磁场分析；热力学中的温度场分析；流体力学中的流场分析等，都可以归结为在给定边界条件下求解其控制方程的问题。虽然人们能够得到它们的基本方程和边界条件，但是能够用解析法求解的只是少数性质比较简单和边界比较规则的问题，实际结构的形状和所受到的载荷往往比较复杂，按解析法求解是非常困难的。2024/7/5有限单元法简介7一、数值模拟方法概述

解决这类复杂问题主要有两种方法：1、引入简化假设，使其达到能用解析法求解的状态，然后求其近似解（未必可行，容易导致不正确的解答）2、保留问题的复杂性，利用数值模拟方法求得问题的近似解（较多采用）数值模拟技术（即CAE技术，Computer-aidedEngineering）是人们在现代数学、力学理论的基础上，借助于计算机技术来获得满足工程要求的数值近似解，是现代工程仿真学发展的重要推动力之一。2024/7/5有限单元法简介8一、数值模拟方法概述

目前在工程技术领域内常用的数值模拟方法有：

1、有限单元法FEM（FiniteElementMethod）

2、边界元法BEM（BoundaryElementMethod）

3、有限差分法FDM（FiniteDifferenceMethod）

4、离散单元法DEM（DiscreteElementMethod）其中有限单元法是最具实用性和应用最广泛的。2024/7/5有限单元法简介9一、数值模拟方法概述

由数值模拟方法结合计算机技术形成的应用软件在工程中得到广泛的应用，国际上著名的有限元通用软件有：ANSYS,ABAQUS,MCS.PATRAN,MCS.NASTRAN,MCS.MARC,ADINA,FLAC等它们大多采用FORTRAN语言编写，不仅包含多种条件下的有限元分析程序，而且带有强大的前处理和后处理程序。大多数有限元通用软件拥有良好的用户界面、使用方便，功能强大。2024/7/5有限单元法简介10[FEM求解工程问题思路]

工程中的问题本构关系线性的、边界规则的问题数值分析法精确解近似解非线性的、边界不规则的问题解析法图1工程问题的求解思路

FEM犹如万能钥匙，是解决工程问题的主流数值分析方法只能解决极少的方程求解，而且解答复杂。工程问题求解的基本过程InstituteofMechanicalEngineeringandAutomation二、有限单元法简介2024/7/5有限单元法简介111.有限元法定义

是一种工程物理问题的数值分析方法，根据近似分割和能量极值原理，把求解区域离散为有限个单元的组合，研究每个单元的特性，组装各单元，通过变分原理，把问题化成代数方程组求解。分析指导思想

化整为零，裁弯取直，以简驭繁，变难为易二、有限单元法简介2024/7/5有限单元法简介122.发展简史1943年，Courant提出有限元法概念1956年，Turner和Clough第一次用三角形单元离散飞机机翼，借助有限元法概念研究机翼的强度及刚度1960年，Clough正式提出有限元法（FEM）20世纪60年代，我国数学家冯康把FEM总结成凡是椭圆形偏微分方程都可用FEM求解20世纪60年代以后，由于数学界的参与，FEM得到蓬勃发展，并且扩大了应用二、有限单元法简介2024/7/5有限单元法简介133.发展方向新型单元的研究有限元的数学理论向新领域扩展应用大型通用程序的编制和设计

ANSYS,NASTRAN,ABAQUS开发微机用版本设计自动化及优化设计（CAD,CAE,CAM）二、有限单元法简介2024/7/5有限单元法简介144.有限元法的分类

以方程中未知数代表的意义分类

有限元位移法：未知数为位移有限元力法：未知数为力有限元混合法：未知数为力和位移

以推导方法分类

直接法变分法加权余数法二、有限单元法简介2024/7/5有限单元法简介15有限单元法的常用术语：二、有限单元法简介真实系统有限元模型

有限元模型是真实系统理想化的数学抽象。定义2024/7/5有限单元法简介16自由度（DOFs－degreeoffreedoms）自由度(DOFs)用于描述一个物理场的响应特性。结构DOFs结构 位移热温度电 电位流体压力磁 磁位物理场 自由度ROTZUYROTYUXROTXUZ2024/7/5有限单元法简介17节点和单元节点:空间中的坐标位置，具有一定自由度和存在相互物理作用。单元: 一组节点自由度间相互作用的数值、矩阵描述（称为刚度或系数矩阵)。单元有线、面或实体以及二维单元或三维单元等。有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。载荷载荷2024/7/5有限单元法简介18节点和单元(续)节点自由度随连接该节点的单元类型变化。JIIJJKLILKIPOMNKJIL三维杆单元(铰接)UX,UY,UZ三维梁单元二维或轴对称实体单元UX,UY三维四边形壳单元UX,UY,UZ,三维实体热单元TEMPJPOMNKJIL三维实体结构单元ROTX,ROTY,ROTZROTX,ROTY,ROTZUX,UY,UZ,UX,UY,UZ2024/7/5有限单元法简介19单元形函数FEA（Finiteelementanalysis）仅仅求解节点处的DOF值。单元形函数是一种数学函数，规定了从节点DOF值到单元内所有点处DOF值的计算方法。因此，单元形函数提供出一种描述单元内部结果的“形状”。单元形函数描述的是给定单元的一种假定的特性。单元形函数与真实工作特性吻合好坏程度直接影响求解精度。2024/7/5有限单元法简介20真实的二次曲线.节点单元二次曲线的线性近似

(不理想结果).2单元形函数(续)节点单元DOF值二次分布..1节点

单元线性近似(更理想的结果)真实的二次曲线.....3节点单元二次近似(接近于真实的二次近似拟合)

(最理想结果)..42024/7/5有限单元法简介21单元形函数(续)遵循:

DOF值可以精确或不够精确地等于节点处的真实解，但单元内的平均值与实际情况吻合得很好。这些平均意义上的典型解是从单元DOFs推导出来的（如，结构应力、热梯度）。如果单元形函数不能精确描述单元内部的DOFs，就不能很好地得到导出数据，因为这些导出数据是通过单元形函数推导出来的。遵循原则:

当选择了某种单元类型时，也就十分确定地选择并接受该种单元类型所假定的单元形函数。在选定单元类型并随之确定了形函数的情况下，必须确保分析时有足够数量的单元和节点来精确描述所要求解的问题。2024/7/5有限单元法简介22二、有限单元法简介有限单元法的基本思想：1、将一个连续域离散化为有限个单元，并通过有限个节点相连接的等效集合体。由于单元能按照不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模型化几何形状复杂的求解域。2、有限元法利用在每一个单元内假设的的近似函数来分片地表示全求解域上待求的未知场函数。单元内的近似函数由未知场函数在单元的各个节点的数值和其插值函数来表达。3、一个问题的有限元分析中，未知场函数在各个节点上的数值就成为新的未知量，从而使一个连续的无限自由度问题变成离散的有限自由度问题。4、一经求解出这些未知量，就可以通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上的近似解。显然，随着单元数目的增加，也即单元尺寸的缩小，或者随着单元自由度的增加以及插值函数精度的提高，解的近似程度将不断改进，如果单元是满足收敛要求的，近似解最后将收敛于精确解。2024/7/5有限单元法简介23三、有限单元法分析步骤有限元法分析问题的基本步骤：1、结构的离散化

离散化就是将要分析的结构分割成有限个单元体，并在单元的指定位置设置节点，使相邻单元的有关参数具有一定的连续性，构成单元的集合体代替原来的结构。结构离散化时，划分的单元大小和数目应根据计算精度的要求和计算机的容量来决定

选取坐标（右手法则）选择合适的单元，离散结构物为有限个单元，并对单元、节点进行编号2024/7/5有限单元法简介24三、有限单元法分析步骤2、选择位移插值函数

为了能用节点位移表示单元体的位移、应变和应力，在分析连续体问题时，必须对单元中位移的分布做出一定的假设，一般假定位移是坐标的某种简单函数。选择适当的位移函数是有限单元法中的关键。

{f}—单元内任意点的位移列矩阵[N]—单元形函数矩阵—单元节点位移的列矩阵2024/7/5有限单元法简介25三、有限单元法分析步骤3、分析单元的力学特性

利用几何方程、本构方程和变分原理得到单元的刚度矩阵和载荷矩阵{R}e=[K]e{δ}e{R}e—单元节点力[K]e—单元刚度矩阵2024/7/5有限单元法简介26三、有限单元法分析步骤4、集合所有单元平衡方程，得到整体结构的平衡方程

先将各个单元刚度矩阵集合成整体刚度矩阵[K]，然后将各单元的等效节点力列阵集合成总的载荷列阵[K]=Σ[K]e5、由平衡方程求解未知节点位移{δ}

按照问题的边界条件修改总的平衡方程，并进行求解。2024/7/5有限单元法简介27三、有限单元法分析步骤6、单元应变和应力的计算

根据已知结点的位移利用弹性力学方程和位移插值函数算出单元的应变和应力。2024/7/5有限单元法简介28四、利用有限元软件进行工程分析

所有的通用有限元软件都包括：前处理、求解器、后处理三个有逻辑顺序的模块。在进行实际工程分析时，也该按照以上三个模块来进行。进入求解器进行求解（设定分析步骤，输出变量）前处理（建模、材料特性、单元选择及网格划分）进入后处理（变形图、等值线图，列表显示等等后处理）2024/7/5有限单元法简介29四、利用有限元软件进行工程分析利用有限元软件进行工程问题的分析，一般应按下列步骤进行：（一）、制订分析方案需考虑以下几个方面：important2024/7/5有限单元法简介301、分析领域

几何体

载荷

物理系统结构热电磁2024/7/5有限单元法简介312、分析目标力？位移？温度？还是其他？2024/7/5有限单元法简介323、线性/非线性分析“我的物理系统是在线性还是非线性状态下工作？线性求解能满足我的需要吗？如果不能，必须考虑哪种非线性特性？”许多情况和物理现象都要求进行非线性计算。(a)订书钉t0t1t2t3Fu

(b)木制书架b1b2

(c)气动带FuFu2024/7/5有限单元法简介334、静力/动力分析静力求解能否满足你的分析要求？如果不能，应当进行哪种动力分析？动力分析的所有载荷都是随时间变化的，但在许多情况下动力影响可以忽略不计。一般情况下，激励频率低于结构最小固有频率的1/3时静力求解就足够了。惯性力是动力问题不同于静力问题的关键之处。2024/7/5有限单元法简介345、分析细节的考虑在建立分析模型之前必须制订好建模方案：必须考虑哪些细节问题？对称/反对称/轴对称？模型中存在应力奇异？选用哪种类型的单元？线单元壳单元X－Y平面单元平面应力或应变单元轴对称单元谐单元实体单元专用单元线性单元/高阶单元/P单元四边形单元/三角形单元，块单元/四面体单元2024/7/5有限单元法简介356、充分利用结构的对称性PPPPP2024/7/5有限单元法简介367、网格密度相邻单元的尺寸尽可能接近应力变化大处单元应密集一些。结点的多少与疏密要考虑计算机的容量和计算精度结点所连接的单元个数尽可能一致宜不宜2024/7/5有限单元法简介378、单位制注意：ANSYS和ABAQUS大型有限元软件中，没有固定的单位制，大家在使用的过程中，可以自己选用前后一致的一套单位制，则最后所得结果的单位即为即为所选单位制对应的单位。建议：尽量采用国际单位制2024/7/5有限单元法简介389、材料特性材料特性是有限元分析必须提供的数据，其准确与否直接影响到计算的精度；必要的时候需通过试验提供；一个复杂分析中可能包含很多种性质截然不同的材料，建模的时候应以足够的关键字以识别；很多有限元软件中，都提供用户接口，若软件中没有现成的材料模型，用户