有限元

有限元

本词条由“科普中国”科学百科词条编写与应用工作项目 审核 。

在数学中，有限元法（FEM，Finite Element Method）是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解，每个子区域都成为简单的，这种简单部分就称作有限元。

它通过变分方法，使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想，有限元法包含了一切可能的方法，这些方法将许多被称为有限元的小区域上的简单方程联系起来，并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的（较简单的）近似解，然后推导求解这个域总的满足条件（如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

计算简介

编辑

英文：Finite Element

有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域—飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛地应用于求解热传导、电磁场、流体力学等连续性问题。

方法

编辑

有限元法分析计算的思路和做法可归纳如下：

物体离散化

将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。

选择位移模式

在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。

当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法我们就将位移表示为坐标变量的简单函数，这种函数称为位移模式或位移函数。

分析力学性质

根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。

等效节点力

物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上的力。

 

 

有限元中单元节点和积分点有什么区别

 

原文地址：http://blog.sina.com.cn/s/blog_9187667601013lsw.html

有限元中单元节点和积分点有什么区别

学过数值积分的应该知道，有限元中的积分点指高斯积分点，因为这些点的收敛性好，精度高。

1. 节点
    在单元内，采用形函数来表述单元内变量的分布规律。而节点值是在节点处的对应物理量。
以简单矩形单元的温度为例：四个节点i,j,m,n的温度分别为Ti,Tj,Tm,Tn.
则以单元内自然坐标(x,y),(-1,-1),(-1,1),(1,-1),(1,1)分别为四个节点，单元内温度分布为：

T={Si, Sj, Sm, Sn} {Ti, Tj, Tm, Tn}
其中，Si=1/4(1-x)(1-y)，Sj=1/4(1+x)(1-y)] ，Sm=1/4(1+x)(1+y)，Sn=1/4(1-x)(1+y)

(单元的形函数我们可以从手册中查到)
从而我们知道了温度在单元内的分布。
2. 积分点
    我们需要对温度在单元内的面积上进行积分时，因为节点的温度显然与x,y无关，我们只需要考虑对形函数积分。采用Gauss_Legendre多项式计算积分时，我们只需要计算根据特定积分点的值（在自然坐标系下是固定的，可以查手册，这些点也叫高斯点、积分点）并加以权重就可以。这就把复杂的积分问题变成了简单的代数问题。因为形函数只与单元有关，所以积分点也只与单元形状有关。
    应力一般采用多个积分点的相互插值或外延来计算节点应力。这只是为了减少误差。因为在积分点应力比节点具有更高阶的误差。
    从理论上说，形函数已知后，用Maple或者Mathematic等软件进行符号积分的话，是可以精确计算出刚度矩阵和质量矩阵，但是这样做的话，对于工程实际应用来说并不合适。
原因：1，费时；2，Mindlin中厚板有剪力锁死问题，有时候需要采用缩聚积分），所以有些书上会把2节点梁单元的刚度阵直接写出来，但是再复杂点的单元，就使用数值积分（Newton-Cotes积分和高斯积分）
牛顿-科斯的积分点就是节点，这样得到的质量矩阵是集中质量阵形式

 

个人理解：

1.节点作用：构造形函数，节点的多少描述规则形状单元内的应力的近似分布情况，并获取节点上的位移值

2.积分点作用：构造规则形状单元与曲边（曲面）单元的转化的变换函数，积分点的选取多少和选取的位置直接关系到这种“映射”的精确程度，刚度矩阵、边界条件的转化都用到了坐标变换的积分关系，一般取高斯积分点能使被积函数计算精度尽量高。对于newton-cote积分点的选取，这种“映射”看起来，节点和积分点是同一个位置或说是同一点，而对于高斯积分点位置与节点是不同的。

故有如下结果：

 1.由于高斯积分点的这种变换比较高，在方程求解结束，返回积分点上的应力解比较准确。
 2.至于Mindlin中厚板有剪力锁死问题，采用缩聚积分，也是应为这种坐标的变换关系（可见《有限单元法基本原理和数值方法p345页10.4.11式可知），力的边界条件只有剪切，采用缩聚积分可以较大降低剪切力的影响，但是也可能引起刚度矩阵的奇异，所以对于中厚板的积分点选取不同一般的方案。

1.ANSYS手册(Chapter 13)上列出各种单元的积分点位置。
2.王瑁成的《有限单元法》第五章，有解释为什么积分点应力更加精确。

3.因为积分点应力更精确，所以我们一般采用积分点的应力内插或外延确定节点应力。特殊情况除外。

 

单元节点和积分点是不同的两个概念！