有限元法分析过程

有限元法的分析过程有限元法是一种数值分析方法，用于求解实际问题的物理场或结构的数学模型。

它将连续的实体分割成离散的小单元，通过建立节点和单元之间的关系，对物理问题进行逼近和求解。

以下是一般的有限元法分析过程。

1.问题建模和离散化在有限元分析中，首先需要对实际问题进行建模，确定物理场或结构的几何形状和边界条件。

然后，将几何形状分割成一系列小单元，例如三角形、四边形或四面体等。

2.网格生成根据问题的几何形状和离散化方式，生成网格。

网格是由一系列节点和单元组成的结构，节点用于描述问题的几何形状，单元用于划分问题域。

通常，节点和单元的位置和数量会直接影响有限元法的精度和计算效率。

3.插值函数和基函数的选择有限元法中的节点通常表示问题域中的几何点，而节点之间的关系由插值函数或基函数来描述。

插值函数用于建立节点和单元之间的关系，基函数用于对物理场进行逼近。

选择适当的插值函数和基函数是有限元法分析的关键。

4.定义系统参数和边界条件确定相关物理参数和材料性质，并将其转化为数值形式。

在有限元分析中，还需要定义边界条件，包括约束条件和加载条件。

5.定义数学模型和方程根据问题的物理场或结构和所选择的基函数，建立数学模型和方程。

有限元方法可以用来建立线性方程、非线性方程、静态问题、动态问题等。

具体建立数学模型和方程的过程需要根据问题的特点进行。

6.组装刚度矩阵和力载荷向量根据离散化的节点和单元，组装刚度矩阵和力载荷向量。

刚度矩阵描述节点之间的刚度关系，力载荷向量描述外部加载的作用力。

7.求解代数方程通过求解代数方程，确定节点的位移或物理场的数值解。

通常，使用迭代方法或直接求解线性方程组的方法来求解。

8.后处理和分析得到数值解后，可以进行后处理和分析。

包括计算节点和单元的应变、应力等物理量，进行矫正和验证计算结果的正确性。

还可以通过有限元法的网格适应性来优化问题的计算效率和精度。

以上是一般的有限元法分析过程，具体的步骤和方法可能会因不同的问题而有所不同。

有限元分析的一般过程一、结构的离散化将结构或弹性体人为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。

这一步要解决以下几个方面的问题:1、选择一个适当的参考系，既要考虑到工程设计习惯，又要照顾到建立模型的方便。

2、根据结构的特点，选择不同类型的单元。

对复合结构可能同时用到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。

3、根据计算分析的精度、周期及费用等方面的要求，合理确定单元的尺寸和阶次。

4、根据工程需要，确定分析类型和计算工况。

要考虑参数区间及确定最危险工况等问题。

5、根据结构的实际支撑情况及受载状态，确定各工况的边界约束和有效计算载荷。

二、选择位移插值函数1、位移插值函数的要求在有限元法中通常选择多项式函数作为单元位移插值函数，并利用节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。

位移插值函数需要满足相容（协调）条件，采用多项式形式的位移插值函数，这一条件始终可以满足。

但近年来有人提出了一些新的位移插值函数，如：三角函数、样条函数及双曲函数等，此时需要检查是否满足相容条件。

2、位移插值函数的收敛性（完备性）要求：1）位移插值函数必须包含常应变状态。

2）位移插值函数必须包含刚体位移。

3、复杂单元形函数的构造对于高阶复杂单元，利用节点处的位移连续性条件求解形函数，实际上是不可行的。

因此在实际应用中更多的情况下是利用形函数的性质来构造形函数。

形函数的性质：1）相关节点处的值为 1，不相关节点处的值为 0。

2）形函数之和恒等于 1。

1、建立数学模型（特征消隐，理想化，清除）（（即从 CAD 几何体→FEA 几何体），共有下列三法：▲ 特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。

▲ 理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理▲ 清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。

线弹性有限元法分析的基本步骤有限元法是工程领域应用最为广泛的一种计算方法，它不但可以解决工程中的结构分析问题，而且已成功地解决了热力学、流体力学、电磁学和声学等领域的问题。

经过数十年的发展，有限元方法的理论已相当完善。

将有限元理论、计算机图形学以及优化技术相结合而开发的各类专用有限元软件。

能高速高效地解决各类有限元问题。

线弹性有限元法分析的基本步骤线弹性有限元法分析有限元一般分为6步。

第一步：定义形函数)(x N ，进而通过单元节点变量a 描述单元域内连续的变量)(x ue a x N x u )()(=第二步：定义单元材料的响应；如应力、应变和热流等。

e Ba x u L x ==)]([)(ε)()(x D x εσσ==第三步：形成单元矩阵，建立单元与外界的平衡关系:0=+e e e f a K式中e K ：单元刚度矩阵e a ：单元节点位移e f ：单元节点上的等效外力刚度矩阵按下式计算：⎰Ω=e DBdV B K T e节点等效外力为： F tdS x N bdV x N f ee T T e++=-⎰⎰ΓΩ)()( 单元刚度矩阵代表了结构单元的刚度或传热单元的传导性，节点外载荷单元内体力（或内热源）、面力（或热流）、及节点集中力的贡献。

第四步：集成。

将覆盖结构全域的所有单元的刚度矩阵和节点外力对平衡的贡献集成，建立整体结构的平衡方程。

∑=ee K K∑=ee f f第五步：求解平衡方程。

指定一些节点位移后，可将平衡方程按已知节点位移和未知的场变量分解为两部分。

⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡r a s u ss suus uu f f a a K K K K 式中 u a 为未知节点变量s a 为已知节点变量a f 为外加的节点力r f 为节点反作用力求解变成获得已知位移的节点反作用力和已知的载荷节点位移。

第六步：回代。

根据计算出的节点变量，代入第二步的表达式中，获得单元应变、应力或热流等量。

有限元分析过程有限元分析过程可以分为以下三个阶段：1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型——有限元模型，从而为有限元数值计算提供必要的输入数据。

有限元建模的中心任务是结构离散，即划分网格。

但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。

2.计算阶段: 计算阶段的任务是完成有限元方法有关的数值计算。

由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成。

3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理，并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。

注意：在上述三个阶段中，建立有限元模型是整个有限分析过程的关键。

首先，有限元模型为计算提供所以原始数据，这些输入数据的误差将直接决定计算结果的精度；其次，有限元模型的形式将对计算过程产生很大的影响，合理的模型既能保证计算结构的精度，又不致使计算量太大和对计算机存储容量的要求太高；再次，由于结构形状和工况条件的复杂性，要建立一个符合实际的有限元模型并非易事，它要考虑的综合因素很多，对分析人员提出了较高的要求；最后，建模所花费的时间在整个分析过程中占有相当大的比重，约占整个分析时间的70%，因此，把主要精力放在模型的建立上以及提高建模速度是缩短整个分析周期的关键。

原始数据的计算模型，模型中一般包括以下三类数据：1.节点数据: 包括每个节点的编号、坐标值等；2.单元数据：a.单元编号和组成单元的节点编号；b.单元材料特性，如弹性模量、泊松比、密度等；c.单元物理特征值，如弹簧单元的刚度系数、单元厚度、曲率半径等；d.一维单元的截面特征值，如截面面积、惯性矩等；e.相关几何数据3.边界条件数据：a.位移约束数据；b.载荷条件数据；c.热边界条件数据；d.其他边界数据.建立有限元模型的一般过程：1.分析问题定义在进行有限元分析之前，首先应对结果的形状、尺寸、工况条件等进行仔细分析，只有正确掌握了分析结构的具体特征才能建立合理的几何模型。

第十一章 有限元分析方法概述1、基本概念有限元分析方法是随着电子计算机的发展而迅速发展起来的一种现代没计计算方法。

它是20世纪50年代首先在连续体力学领域—飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快就广泛地应用于求解热传导、电磁场、流体力学等连续性问题。

在工程分析和科学研究中，常常会遇到大量的由常微分方程、偏微分方程及相应的边界条件描述的场问题，如位移场、应力场和温度场等问题。

求解这类场问题的方法主要有两种：用解析法求得精确解；用数值解法求其近似解。

应该指出，能用解析法求出精确解的只是方程性质比较简单且几何边界相当规则的少数问题。

而对于绝大多数问题，则很少能得出解析解。

这就需要研究它的数值解法，以求出近似解。

目前工程中实用的数值解法主要有三种：有限差分法、有限元法和边界元法。

其中，以有限元法通用性最好，解题效率高，目前在工程中的应用最为广泛。

下面通过一个具体例子，分别采用解析法和数值解法进行求解，从而体会一下有限元分析方法的含义及其相关的一些基本概念。

如下图所示为一变横截面杆，杆的一端固定，另一端承受负荷P ，试求杆沿长度方向任一截面的变形大小。

其中，杆的上边宽度为1w ，下边宽度为2w ，厚度为t ，长度为L ，杆的材料弹性模量为E 。

已知P ＝4450N ，1w ＝50mm ，2w =25mm ，t =3mm ，L =250mm ，E =72GPa 。

① 采用解析法精确求解假设杆任一横截面面积为)(y A ，其上平均应力为σ，应变为ε。

根据静力平衡条件有：0)(=-y A P σ根据虎克定律有：εσE =而任一横截面面积为：t y L w w w y A )()(121-+= 任一横截面产生的应变为：dydu=ε将上述方程代入静力平衡条件，进行变换后有：dy y EA Pdu )(=沿杆的长度方向对上式两边进行积分，可得：⎰⎰⎰-+==y yudy y Lw w w Et P dy y EA P du 01210)()(将)(y A 表达式代入上式，并对两边进行积分，得杆沿长度方向任一横截面的变形量：]ln )[ln()()(112112w y Lw w w w w Et PL y u --+-=当y 分别取0、62.5、125、187.5、250值时，变截面杆相应横截面处的沿杆长方向的变形量分别为：m u m u m u m u m u 6564636211080.142 ;1083.96 ;1027.59 ;1051.27 ;0----⨯=⨯=⨯=⨯==② 采用数值解法近似求解将变横截面杆沿长度方向分成独立的4小段，每一小段采用等截面直杆近似，等截面直杆的横截面面积为相应的变截面杆横截面面积的平均面积表示，每一小段称为一个单元，小段之间通过节点连接起来。

有限元分析过程：一，结构离散化1.选择单元类型2.单元划分；二，单元分析1.选择位移函数2.分析单元力学特性；三，整体分析1.集成整体结点载荷向量2.集成整体刚度方程3.引进边界约束条件，解总体刚度方程求出结点位移分量。

位移模式应满足下列收敛性条件：完备性 1.位移模式必须包含单元的常应变状态；2.位移模式必须包含单元的刚体位移；协调性 3.位移模式应尽可能反映位移的连续性。

单元刚度矩阵的性质：1.对称性；2.单元刚度矩阵与单元位置无关；3.奇异性。

总体刚度矩阵的性质：1.稀疏性；2.带状性；3.奇异性与对称性。

由单元刚度方程组集总纲时应满足的原则：1各单元在公共节点上协调地彼此连接，即在公共结点处具有相同的位移2结构的各节点离散出来后应满足平衡条件提高单元精度的方法:1增加结点数即提高位移模式的阶次2建立等参单元进行等参数变换等参数变换、等参数单元、等参单元具有哪些优越性？：1将局部坐标中几何形状的单元转换成总体坐标中几何形状复杂的单元且这种坐标变换和函数插值采用了相同数目的结点数参数和相同的插值函数2采用等参数变换的单元称为等参数单元3优点：可以很方便地用来离散具有复杂性体的结构。

由于等参变换的采用使等参单元特性矩阵的计算仍在单元的规则域内进行，边界条件：位移边界条件和应力边界条件引进位移边界条件的方法：对角元素改一和乘大数弹性力学中求解力学位移的方法：解析法或半解析法、数值法弹性力学的基本方程：平衡方程（静力平衡关系）、几何方程（应变分量与位移间的关系）、物理方程（应力分量与应变分量之间的关系）什么叫结点力和结点载荷？两者有什么不同？为什么应保留结点力的概念？：①结点力：结点对单元的作用力。

结点载荷：包括集中力和将体力、面力按静力等效原则移植到节点形成的等效载荷，原荷载和移植后的荷载在虚位移上的虚功相等②相对于整体结构来说，节点力是内力，结点载荷是外力③节点力的概念在建立单元刚度方程的时候需要用到在薄板弯曲理论中做了哪些假设？解：①板厚方向的挤压变形可忽略不计。