有限单元法
有限单元法名词解释
有限单元法名词解释
有限单元法（Finite Element Method）是一种数值计算方法，常用于工程领域，用于求解复杂的物理问题。
该方法将连续体分割为有限个小区域，即“单元”，并在每个单元内近似求解。
在有限单元法中，首先将待解问题建模为数学上的形式，选择适当的数学模型
和边界条件。
然后，将物理区域分割为有限个单元，每个单元内的数学形式由逼近函数表示。
每个单元的近似解通过如三角形和四边形等简单形状来表示。
通过解决每个单
元内的数学形式，得到整个物理区域的近似解。
这些单元共同构成了一个有限元模型。
有限单元法的优势在于可以处理各种形状、复杂的物理特性和非线性问题。
它
能够准确地描述材料、结构、流体等领域的行为，并能够提供与实际现象相匹配的数值结果。
此外，有限单元法还能够提供对问题的优化和灵活性，通过改变单元的大小和
形状，可以在所需精度和计算效率之间进行权衡。
总之，有限单元法是一种强大的数值计算方法，应用广泛于各个领域，因其可
靠性和灵活性而受到广泛的青睐。
它是工程分析和设计中不可或缺的工具，为我们解决复杂问题提供了有效的数值模拟手段。
有限单元法原理及应用
有限单元法原理及应用有限单元法（Finite Element Method，FEM）是一种数值分析方法，广泛应用于工程领域中结构力学、流体力学、热传导等问题的数值求解。
它的基本思想是将一个复杂的结构或物理现象分割成有限数量的简单单元，通过对单元的力学行为进行建模，最终得到整个系统的数值解。
本文将围绕有限单元法的原理及其在工程领域中的应用进行详细介绍。
有限单元法的原理。
有限单元法的原理基于力学原理和数学方法，其基本步骤包括，建立数学模型、离散化、单元划分、建立单元刚度矩阵和载荷向量、组装和求解方程、计算结果后处理等。
在建立数学模型时，需要根据实际问题选择合适的数学方程和边界条件，将问题转化为求解一组代数方程。
离散化是指将连续的物理问题划分成若干个小单元，每个单元内的物理行为可以用简单的数学方程描述。
单元划分是将整个结构或领域划分成若干个有限单元，通常采用三角形、四边形、四面体、六面体等几何形状。
建立单元刚度矩阵和载荷向量是对每个单元进行力学行为的建模，根据材料性质和几何形状计算单元的刚度矩阵和载荷向量。
组装和求解方程是将所有单元的刚度矩阵和载荷向量组装成整个系统的刚度矩阵和载荷向量，然后通过数值方法求解代数方程组。
最后，计算结果后处理是对数值解进行分析和可视化，评估结构的性能和稳定性。
有限单元法的应用。
有限单元法在工程领域中有着广泛的应用，包括结构力学、流体力学、热传导等方面。
在结构力学中，有限单元法可以用于分析和设计各种结构，如桥梁、建筑、机械零件等。
通过对结构的受力分析，可以评估结构的安全性和稳定性，指导工程设计和施工。
在流体力学中，有限单元法可以用于模拟流体的流动行为，如水力学、空气动力学等问题的数值模拟。
在热传导中，有限单元法可以用于分析材料的热传导性能，评估材料的热稳定性和散热效果。
总结。
有限单元法作为一种数值分析方法，在工程领域中有着重要的应用价值。
通过对结构、流体、热传导等问题的数值模拟，可以为工程设计和科学研究提供重要的参考和支持。
有限单元法基础
性体在各节点处的位移解。
3、单元分析---三角形单元
y
3.1 单元的结点位移和结点力向量
从离散化的网格中任取一个单元。三个结点 按反时针方向的顺序编号为：i, j, m。
结点坐标: (xi,yi) , (xj,yj) , (xm,ym) 结点位移: (ui,vi) , (uj,yj) , (um,vm) 共有6个自由度
单元位移插值函数： u(x, y) a1 a2 x a3 y
(3.1)
v(x, y) a4 a5x a6 y
插值函数的系数： a1 aiui a ju j amum / 2 A, a4 aivi a jv j amvm / 2 A,
a2 biui bju j bmum / 2 A, a5 bivi bjv j bmvm / 2 A,
um a1 a2 xm a3 ym , vm a4 a5 xm a6 ym ,
求解以上方程组得到以节点位移和节点坐标表示的6个参数：
a1 aiui a ju j amum / 2 A, a4 aivi a jv j amvm / 2 A, a2 biui bju j bmum / 2 A, a5 bivi bjv j bmvm / 2 A, a3 ciui c ju j cmum / 2 A, a6 civi c jv j cmvm / 2 A,
研究方法
从数学上讲它是微分方程边值问题(椭圆型微分方程、抛物型微分方程和双曲型微 分方程)的一种的数值解法，是一种将数学物理问题化为等价的变分问题的解法，并作 为一种通用的数值解法成为应用数学的一个重要分支。从物理上讲是将连续介质物理 场进行离散化，将无限自由度问题化为有限自由度问题的一种解方法。从固体力学上 认识，是瑞利-里兹法的推广。
有限单元法ppt课件
06
有限单元法的发展趋势和展 望
发展趋势
工程应用领域拓展
随着科技的发展，有限单元法在解决 复杂工程问题上的应用越来越广泛， 不仅局限于结构分析，还涉及到流体 动力学、热传导等领域。
与其他方法的结合
有限单元法正与其他数值方法（如有 限差分法、边界元法等）进行交叉融 合，形成更为强大的数值分析工具。
05
有限单元法的优缺点
优点
灵活性
有限单元法允许对复杂的几何形状进 行离散化，适用于解决各种形状和大 小的问题。
高效性
有限单元法能够处理大规模问题，通 过使用计算机技术，可以快速求解。
广泛的应用领域
有限单元法被广泛应用于工程、物理 、生物等领域，是一种通用的数值分 析方法。
易于理解和实现
有限单元法的基本概念直观易懂，且 实现起来相对简单。
01
利用线性代数方法，将 各个单元的数学模型和 节点信息组合成整体方
程组。
03
将节点的未知量返回到 原问题中，得到问题的
解。
05
根据问题的物理性质和 边界条件，建立单元的 数学模型和节点信息。
02
解整体方程组，得到节 点的未知量。
04
有限单元法的特点
适用范围广
可以用于解决各种类型的问题，如弹性力学 、流体力学、传热学等。
高精度与高效率
研究者们致力于开发更高效、精确的 算法，以解决大规模、非线性、动态 等复杂问题。
并行化与云计算应用
随着计算资源的丰富，有限单元法的 计算过程正逐步实现并行化，利用云 计算平台进行大规模计算已成为趋势 。
展望
理论完善与创新
随着工程实践的深入，有限单元法的理论体系将进一步完善，同时会 有更多创新性的算法和模型出现。
有限单元法
军事工业
子弹穿甲模拟
生物力学 人体颈椎后仰应力分析
3.
怎样学习有限元 怎样学习有限元 ？
两条原则
一、掌握有限元的基本理论和概念，这是灵活 掌握有限元的基本理论和概念， 使用有限元方法的基础。 使用有限元方法的基础。 二、掌握一种有限元分析软件（ANSYS）， 掌握一种有限元分析软件（ ）， 做工程实际课题，在实践中提高。 做工程实际课题，在实践中提高。
有限元分析 是利用数学近似的方法对真实物理系统
（几何和载荷工况）进行模拟。利用简单而又相互作 几何和载荷工况）进行模拟。 用的元素，即单元， 用的元素，即单元，用有限数量的未知量去逼近无限 未知量的真实系统。 未知量的真实系统。
1.2有限单元法的基本思想
1. 结构的离散化：假想把连续系统(杆件、连续体、 连续介质)分割成数目有限的单元，单元之间只在 数目有限的节点处相互连接，构成一个单元集合 体来代替原来的连续系统。在节点上引进等效载 荷（或边界条件），代替实际作用于系统上的外 载荷（或边界条件）。
J J
铰接) 三维杆单元 (铰接 铰接 UX, UY, UZ
I L
I K
三维梁单元 UX, UY, UZ, ROTX, ROTY, ROTZ
K J
二维或轴对称实体单元 UX, UY I
J P M L I N K J O P
有限单元法原理与应用
有限单元法原理与应用
有限单元法（Finite Element Method，FEM）是一种数值分析
方法，常用于求解复杂的物理问题。
它将连续物体的区域划分为许多小的离散单元，然后在每个单元内建立局部的数学模型和方程。
通过求解这些局部模型和方程，可以得到整个物体的行为和性能。
有限单元法的基本原理是将连续问题离散化为有限数目的独立子问题。
在每个小单元内，选择一个数学函数作为近似解，并通过将近似解与原问题的偏微分方程进行数值积分和数值迭代，得到近似解的解析解。
将每个小单元的解汇总起来，可以得到整个物体的解。
有限单元法的应用非常广泛，可以用于解决各种工程和科学领域的问题。
例如，它可以用来模拟结构的强度和刚度特性，预测材料的疲劳寿命，优化产品的设计，以及研究流体和热传导等问题。
在建筑工程中，有限单元法可以用来分析建筑结构的荷载和变形，评估结构的安全性。
在汽车制造业中，它可以用来模拟车辆的碰撞和破碎行为，提高车辆的安全性。
在航空航天领域，有限单元法可以用来优化飞机的结构和翼型，提高飞机的性能。
此外，有限单元法还可以应用于地震工程、地下水流动、电磁场分析等领域。
总之，有限单元法通过离散化连续问题，将其转化为独立的子问题，然后通过求解局部模型和方程，得到整体解。
它具有广泛的应用领域，为解决多种复杂问题提供了有效的数值分析方法。
有限单元法 数学术语
有限单元法有限单元法，是一种有效解决数学问题的解题方法。
其基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。
采用不同的权函数和插值函数形式，便构成不同的有限元方法。
有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。
内容简述在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。
在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。
从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。
令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。
插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。
有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。
有限单元法基础
有限单元法基础
有限单元法（Finite Element Method，FEM）是一种数值计算
方法，常用于求解连续介质力学问题。
它将连续的物理域划分为有限数量的离散单元（finite elements），通过在每个单元内构建近似函数来描述物理场，再根据物理方程建立离散方程组，通过求解离散方程组来得到物理场的近似解。
有限单元法的基本思路是将连续域离散化为有限数量的小单元，每个小单元内使用适当的数学函数进行插值，将大问题分解为很多个小问题，并利用变量之间的连续性建立全局的离散方程组。
然后通过求解离散方程组得到近似解。
有限单元法的基本步骤包括：
1. 网格划分：将要求解的区域划分为多个离散单元，并在每个单元内选择适当的形状函数。
2. 形函数构造：在每个单元内选择适当的形状函数，用于描述物理场的分布。
3. 整体方程组：根据物理方程在每个单元上的积分，建立整个问题的离散方程组。
4. 边界条件：根据边界条件，将边界上的节点处的值固定为已知值。
5. 求解方程组：利用数值方法求解离散方程组，得到物理场的
近似解。
6. 后处理：根据求解结果，计算所需的物理量并进行分析和验证。
有限单元法具有广泛的应用，适用于各种连续介质力学问题的数值求解，如结构力学、固体力学、流体力学、热传导等。
它可以处理复杂的几何形状和边界条件，且精度和收敛性能较高。
有限单元法
有限单元法
答：一、定义
有限单元法的基本前提是：将连续的求解域离散为一组有限个单元的组合体，这样的组合体能解析地模拟或逼近求解区域。
由于单元能按各种不同的连接方式组合在一起，且单元本身又可以有不同的几何形状，因此可以模型化几何形状复杂的求解域，有限元法作为一种数值分析方法的另一重要步骤是利用在每一个单元内假设的近似函数来表示全求解区域上待求的未知场函数。
单元内的近似函数通常由未知场函数在各个单元节点上的数值以及插值函数表达。
这样一来，一个问题的有限单元分析中，未知场函数的节点值就成为新的未知量，从而使一个连续的无限自由度问题变成离散的有限自由度问题。
一旦求解出这些未知量，就可以利用插值函数确定单元组合体上的场函数。
显然，随着单元数日的增加，单元尺寸的缩小，解的近似程度将不断改进，如果单元是满足收敛要求的，近似解最后将收敛于精确解。
二、有限单元法主要学什么
1、有限单元法主要讲述线弹性有限元法的基本理论、matlab编程实现及相应商业有限元软件的应用，对线弹性动力有限元法及材料、几何和接触三类非线性有限元法的基本概念和程序应用也进行了介绍。
2、主要内容是：matlab编程及符号运算、分部积分、泛函极值与变分法、直接刚度法、有限元求解方法、杆单元力学基础、单元组
装、弹性固体结构、板壳结构。
有限单元法简介
3．非线性边界(接触问题) 在加工、密封、撞击等问题中，接触和摩擦 接触和摩擦的作用不可忽 接触和摩擦 视，接触边界属于高度非线性边界。 平时遇到一些接触问题，如： • 齿轮传动； • 冲压成型； • 轧制成型； • 橡胶减振器； • 紧配合装配等 当一个结构与另一个结构或外部边界相接触时通常要考虑 非线性边界条件。实际的非线性可能同时出现上述两种或三种 非线性问题。
(2)用每个单元内所假设的近似函数来分片地表示全求解 域内待求的未知场变量。 • 每个单元内的近似函数由未知场函数(或其导数)在单元各 个节点上的数值和与其对应的插值函数来表达(此表达式 通常表示为矩阵形式)。 • 由于在联结相邻单元的节点上，场函数应具有相同的数 值，因而将它们用作数值求解的基本未知量。
2．几何非线性问题 当物体的位移较大时，应变与位移的关系是非线性关系 应变与位移的关系是非线性关系，这意味 应变与位移的关系是非线性关系 着结构本身会产生大位移或大转动，而单元中的应变却可大可小。 研究这类问题时一般都假定材料的应力与应变呈线性关系 假定材料的应力与应变呈线性关系。 假定材料的应力与应变呈线性关系 这类问题包括： • 大位移大应变问题 如：橡胶部件成形过程 • 大位移小应变问题 如：如结构的弹性屈曲问题
6 有限元法的发展、现状和未来 有限元法的发展、
有限元法的早期工作
•从应用数学的角度考虑，有限元法的基本思想可以追溯到Courant在1943年的工作。 他首先尝试应用在一系列三角形区域上定义的分片连续函数和最小位能原理相结合， 来求解St.Venant扭转问题。 •此后，不少应用数学家、物理学家和工程师分别从不同角度对有限元法的离散理论、 方法及应用进行了研究。 •有限元法的实际应用是随着电子计算机的出现而开始的。首先是Turner，Clough等 人于1956年将刚架分析中的位移法推广到弹性力学平面问题，并用于飞机结构的分 析。他们首次给出了用三角形单元求解平面应力问题的正确解答。三角形单元的特 性矩阵和结构的求解方程是由弹性理论的方程通过直接刚度法确定的。他们的研究 工作开始了利用电子计算机求解复杂弹性力学问题的新阶段。 •1960年Clough进一步求解了平面弹性问题，并第一次提出了“有限单元法”的名称， 使人们更清楚地认识到有限单元法的特性和功效。
有限单元法的解题思路
解有限元方程
选择合适的求解器
根据有限元方程的特点和求解规模，选择合适的求解器，如直接法、迭代法等。
求解有限元方程
利用选择的求解器，求解有限元方程，得到节点自由度的解。
结果后处理与验证
结果后处理
对求解结果进行后处理，提取有用的 信息，如位移分布、应力分布等。
结果验证
将求解结果与实验结果或已知解进行 对比，验证求解的正确性和精度。
边界条件可以分为两类：本质边界条件和自然边界条件。本质边界条件是指那些必须满足的 约束条件，如固定位移、固定载荷等；自然边界条件是指在某些特定条件下系统自动满足的 约束条件，如无滑动、无渗透等。
在有限元分析中，需要对每个单元的边界进行处理，将边界条件转化为对每个单元的约束， 以保证整个系统的能量平衡。
发展
随着计算机技术的进步，有限单元法 在20世纪60年代得到迅速发展，广泛 应用于各种工程领域。
02
有限单元法的基本原理
离散化与有限元
离散化
将连续的物理问题离散化，将连续域 划分为有限个小的单元，每个单元具 有特定的形状和大小。
有限元
在离散化的基础上，选取每个单元的 中心点或节点作为代表点，通过这些 代表点将各个单元连接起来，形成一 个整体的有限元模型。
建立数学模型
01
确定问题类型
明确问题是静态、动态还是流体 问题，以及问题的边界条件和初 始条件。
02
确定物理模型
03
建立数学方程
根据问题类型，建立相应的物理 模型，包括受力分析、位移分析 等。
根据物理模型，建立相应的数学 方程，如平衡方程、运动方程等。
离散化处理
选择合适的单元类型
根据问题特点和求解精度要求，选择合适的单元类型，如一维、 二维或三维单元。
有限单元法
M
e
AL
1
210
0
0
1
0
6
0
11 l 0 210
9 70
13 420
l
1 l2 0
13 l
1
l2
105
420 140
0
1
0
0
6
3
0
9 70
13 l 0 420
13 35
11 210
l
0
13 l 420
1 l2 140
0
11 l 210
1 l2 105
一、基本思想
2、单元矩阵
c0 w1, c1 l1, c2 3w1 3w2 2l1 l2
c3 2w1 2w2 l1 l2
w(x,t) w1(t)w1(x) 1w2 (x) w2 (t)w3 (x) 2 w4 (x)
w1(x) 1 3x l2 2x
w3(x) 3 x l 2 2 x
l
l
3 ,w2 (x) l
u(x,t)
w( x, t )
N
(
x)
q
e
N
u1
0
0
w1
0
w2
u 2
0
0
w3
0
w4
一、基本思想
2、单元矩阵
T 1 2
l 0
A
u 2
w2
Hale Waihona Puke dx1 2l
0
Au
w
u w
dx
1 2
l
A0
qeT
N
T
N
qe
dx
1 2
qe
T M e
有限单元法
1.1 概述
工程和科学中典型问题
工程中的问题 (力学、物理) 各种方程及相应的定解条件 (边界条件及初始条件) 线性的、边界规则的问题 解析法 精确解 非线性的、边界不规则的问题 数值分析法 近似解
目录
上页 下页 退出
1.1 概述
数学物理学
2
有限 单元法
3 计算方法
目录
力学
1
4 计算机技术
上页 下页 退出
有限单元法实质
连续体 单元体
目录
上页 下页 退出
1.1 概述
有限单元法实质
无限自由度 偏微分方程 有限自由度 代数方程
目录
上页 下页 退出
1.1 概述
基本思想
化整为零
目录
集零为整
上页 下页 退出
1.1 概述
目录
上页 下页 退出
1.1 概述
科学研究方法
理论分析
指用现代科学理论进行实际问题分析的方法。
科学实验
人工控制条件下，通过干预和控制科研对象而观察 和探索科研对象有关规律和机制。
目录
上页 下页 退出
科学计算
用计算机处理科学研究和工程技术中的数学计算。 包括建立模型、建立求解的计算方法和计算机实现。
1.1 概述
问题求解方法
目录
连续系统：可以被无限分割，其中的问题只有利用无 上页 穷小的数学观念才能定义，意味着由无限个单元组成。 如一块受力平板，一个活塞，一根轴等。 下页
退出
1.1 概述
工程和科学中典型问题
水坝
已建立了连续系统的基本方 程，由于边界条件的限制， 无法给出精确的解答，要用 近似方法。
目录
上页 下页 退出
有限单元法简介课案课件
06
结论与展望
总结有限单元法的主要内容与特点
总结内容
有限单元法是一种广泛应用于工程和科 学计算中的数值分析方法，其主要思想 是将连续的求解域离散化为一组单元的 组合体，并在每个单元内假设一个近似 函数，然后通过单元组合体的方式求解 整个域的解。其主要特点包括离散化、 单元划分、近似函数和整体组装四个方 面。
有限单元法的物理原理
物理问题的离散化
将连续的物理问题离散化为有限个离 散的单元，每个单元内的物理量（例 如，位移、温度等）可以近似为常数 。
单元之间的相互作用
考虑单元之间的相互作用和边界条件 （例如，位移边界条件、温度边界条 件等），将各个单元连接起来形成一 个整体的求解对象。
有限单元法的应用范围与限制
求解方程
1 2
选择求解器
根据方程的特点和需要，选择合适的求解器进行 求解。
导入求解器
将方程导入到求解器中，进行求解。
3
分析求解结果
根据求解结果，分析方程的解是否符合要求，如 果不符合要求，需要重新进行求解。
结果分析
结果可视化
将求解结果进行可视化处理，生成模型在不同时刻的 状态图。
结果评估
对求解结果进行评估，分析模型的位移、应力、应变 等参数是否符合实际情况。
结果优化
根据结果评估的结果，对模型进行优化设计，提高模 型的性能和稳定性。
04
有限单元法的应用实例
结构分析
总结词
有限单元法在结构分析中得到广泛应用，能够解决各种复杂结构问题。
详细描述
通过将结构离散化为有限个单元，并对每个单元进行受力分析，可以得出结构的整体受力情况和变形，广泛应用 于桥梁、建筑、机械等领域。
分。
有限单元法
线性、二次和三次单元 Plate element
73
74
75
76
77
78
单元划分原则 •误差与单元尺寸成正比
79
•误差与单元单元最小内角正弦成反比
80
•分界线处单元划分
81
•载荷分布
82
网格疏密
83
•凹槽，孔
84
逐步加密
85
86
87
§4.2.3 单元分析
单元分析的主要内容： ➢ 由节点位移求单元内部任一点的位移。 ➢ 根据几何方程，任一点的位移对坐标求偏导得到该
38
例如当气流流过一个很高的铁塔时就会使 铁塔产生变形,而塔的变形又反过来影响到 气流的流动……这就需要用固体力学和流 体动力学的有限元分析结果交叉迭代求解, 即所谓"流固耦合"的问题。
39
•增强可视化的前置建模和后置数据处理功能 早期有限元分析软件的研究重点在于推导新的高
效率求解方法和高精度的单元。随着数值分析方法的 逐步完善,尤其是计算机运算速度的飞速发展,整个计 算系统用于求解运算的时间越来越少,而数据准备和 运算结果的表现问题却日益突出。在现在的工程工作 站上,求解一个包含10万个方程的有限元模型只需要 用几十分钟。
31
在国内，我国数学家冯康独立于西方提出了有限 元法。1965年，他发表论文《基于变分原理的差 分格式》，标志着有限元法在我国的诞生。冯康 的这篇文章不但提出了有限元法，而且初步发展 了有限元法。他得出了有限元法在特定条件下的 表达式，独创了“冯氏大定理”并且初步证明了 有限元法解的收敛性。虽然冯康创造的有限元法 不成熟，但他能在当时的条件下独立提出有限元 法已十分不易。对于他的这项成就，国内外专家 学者和国家领导人都有很高的评价。
有限单元法基本原理和数值方法 (2)
有限单元法基本原理和数值方法1. 引言有限单元法（Finite Element Method，简称FEM）是一种用于求解工程问题的数值计算方法。
它的基本原理是将连续体分割为离散的有限单元，通过建立有限单元间的关系，近似求解连续体的行为。
本文将介绍有限单元法的基本原理和数值方法。
2. 有限单元法基本原理有限单元法基于两个基本假设：一是一个连续物体可以用小的有限单元来近似表示；二是连续物体在每个有限单元内有近似均匀的力和位移。
有限单元法的基本原理可以概括为以下几个步骤：2.1 离散化将连续物体划分为有限个离散的单元，每个单元都有自己的性质和参数。
通常采用三角形、四边形、四面体等简单形状的单元。
2.2 建立单元间的关系通过节点和单元之间的连接关系来构建整个有限元模型。
每个单元都与相邻的单元共享一些节点，通过共享的节点建立单元间的关系。
2.3 定义单元的属性为每个单元定义材料性质、几何属性和荷载条件等参数，这些参数将用于描述单元的行为。
2.4 定义求解问题的边界条件为有限元模型定义相应的边界条件，如位移边界条件、力边界条件等。
2.5 利用单元间的关系建立方程通过应变能最小原理，利用单元间的关系建立求解整个结构的方程。
2.6 求解方程将建立的方程离散化，采用数值方法求解得到解。
3. 有限单元法数值方法有限单元法中常用的数值方法有直接法和迭代法。
3.1 直接法直接法是指直接求解线性方程组的方法，通常使用高斯消元法、LU分解法等。
直接法的优点是计算简单，稳定性好。
但是当方程组规模较大时，计算量会很大。
3.2 迭代法迭代法是指通过迭代逼近求解方程组的方法，常用的迭代法有Jacobi迭代法、Gauss-Seidel迭代法等。
迭代法的优点是计算量相对较小，适用于大规模方程组。
但是迭代法的收敛性需要保证，且需要选择合适的迭代停止准则。
4. 有限单元法应用有限单元法广泛应用于工程领域的结构分析、流体力学、电磁场分析等。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
有限单元法
有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。
采用不同的权函数和插值函数形式，便构成不同的有限元方法。
有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。
在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。
从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。
令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。
插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。
有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。
单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。
常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。
在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。
对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。
对于有限元方法，其基本思路和解题步骤可归纳为
(1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。
(2)区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。
区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
(3)确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。
有限元方法中的基函数是在单元中选取的，由于各单元具有规则的几何形状，在选取基函数时可遵循一定的法则。
(4)单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点的参数值)的代数方程组，称为单元有限元方程。
(5)总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。
(6)边界条件的处理：一般边界条件有三种形式，分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。
对于自然边界条件，一般在积分表达式中可自动得到满足。
对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。
(7)解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。