有限元分析与应用——第一章 PPT课件
合集下载
有限元分析课程 第一章 绪论PPT
1 T I [ y ( x)] = ∫ [ y L( y ) + yT f ]d Ω + b.t.( y, g ) Ω 2
其中: b.t.( y, g ) 与边界条件有关。）
14
若假设试探函数只选取一项，即
ϕ ( x ) = α1 ( x − x 2 )
5 易得 α1 = 9 ，则问题的近似解为 5 ϕ ( x) = ( x − x 2 ) 9 变分法的试探函数定义于整个求解域，且必须满足
23
转向机构支架的强度分析
24
动力分析
模态分析—计算线性结构的自振频率及振形. 谱分析—是模态分析的扩展，用于计算由于随机振动引起 的结构应力和应变 (也叫作响应谱).
整机的模态分析
25
谐响应分析—确定线性结构对随时间按正弦曲线变化的载 荷的响应. 旋转设备（如压缩机、发动机、泵、涡轮机械等）的支 座、固定装置和部件； 受涡流（流体的漩涡运动）影响的结构，例如涡轮叶片、 飞机机翼、桥和塔等。 瞬态动力学分析—确定结构对随时间任意变化的载荷的响 应. 可以考虑与静力分析相同的结构非线性行为. 显式动力分析—计算高度非线性动力学和复杂的接触问题。 用于模拟非常大的变形，惯性力占支配地位，并考虑所 有的非线性行为.
L=∫
b a
{ y( x)}
dy 1 + dx dx
2
L依赖于函数y(x)的形式，L随着曲线的形状而变化。L就是函 数y(x)的泛函。 12
假设试探函数为多项式： ϕ ( x) = α1 ( x − x 2 )+α 2 ( x − x 3 )+L +α n ( x − x n +1 )
P
meshing
P
其中: b.t.( y, g ) 与边界条件有关。）
14
若假设试探函数只选取一项，即
ϕ ( x ) = α1 ( x − x 2 )
5 易得 α1 = 9 ，则问题的近似解为 5 ϕ ( x) = ( x − x 2 ) 9 变分法的试探函数定义于整个求解域，且必须满足
23
转向机构支架的强度分析
24
动力分析
模态分析—计算线性结构的自振频率及振形. 谱分析—是模态分析的扩展，用于计算由于随机振动引起 的结构应力和应变 (也叫作响应谱).
整机的模态分析
25
谐响应分析—确定线性结构对随时间按正弦曲线变化的载 荷的响应. 旋转设备（如压缩机、发动机、泵、涡轮机械等）的支 座、固定装置和部件； 受涡流（流体的漩涡运动）影响的结构，例如涡轮叶片、 飞机机翼、桥和塔等。 瞬态动力学分析—确定结构对随时间任意变化的载荷的响 应. 可以考虑与静力分析相同的结构非线性行为. 显式动力分析—计算高度非线性动力学和复杂的接触问题。 用于模拟非常大的变形，惯性力占支配地位，并考虑所 有的非线性行为.
L=∫
b a
{ y( x)}
dy 1 + dx dx
2
L依赖于函数y(x)的形式，L随着曲线的形状而变化。L就是函 数y(x)的泛函。 12
假设试探函数为多项式： ϕ ( x) = α1 ( x − x 2 )+α 2 ( x − x 3 )+L +α n ( x − x n +1 )
P
meshing
P
有限元分析及应用课件
参数设置
设置材料属性、单元类型等参数。
求解过程
刚度矩阵组装
根据每个小单元的刚度，组装成全局的刚度矩阵。
载荷向量构建
根据每个节点的外载荷，构建全局的载荷向量。
求解线性方程组
使用求解器（如雅可比法、高斯消元法等）求解线性方程组，得到节点的位移。
后处理
01
结果输出
将计算结果以图形、表格等形式输 出，便于观察和分析。
有限元分析广泛应用于工程领域，如结构力学、流体动力学、电磁场等领域，用于预测和优化结构的 性能。
有限元分析的基本原理
离散化
将连续的求解域离散化为有限 个小的单元，每个单元具有特
定的形状和属性。
数学建模
根据物理问题的性质，建立每 个单元的数学模型，包括节点 力和位移的关系、能量平衡等。
求解方程
通过建立和求解线性或非线性 方程组，得到每个节点的位移 和应力分布。
PART 05
有限元分析的工程应用实 例
桥梁结构分析
总结词
桥梁结构分析是有限元分析的重要应用之一，通过模拟桥梁在不同载荷下的响应，评估 其安全性和稳定性。
详细描述
桥梁结构分析主要关注桥梁在不同载荷（如车辆、风、地震等）下的应力、应变和位移 分布。通过有限元模型，工程师可以预测桥梁在不同工况下的行为，从而优化设计或进
刚性问题
刚性问题是有限元分析中的一种 特殊问题，主要表现在模型中某 些部分刚度过大，导致分析结果 失真
刚性问题通常出现在大变形或冲 击等动态分析中，由于模型中某 些部分刚度过高，导致变形量被 忽略或被放大。这可能导致分析 结果与实际情况严重不符。
解决方案：为避免刚性问题，可 以采用多种方法进行优化，如采 用更合适的材料模型、调整模型 中的参数设置、采用更精细的网 格等。同时，可以采用多种方法 对分析结果进行验证和校核，以 确保其准确性。
设置材料属性、单元类型等参数。
求解过程
刚度矩阵组装
根据每个小单元的刚度，组装成全局的刚度矩阵。
载荷向量构建
根据每个节点的外载荷，构建全局的载荷向量。
求解线性方程组
使用求解器（如雅可比法、高斯消元法等）求解线性方程组，得到节点的位移。
后处理
01
结果输出
将计算结果以图形、表格等形式输 出，便于观察和分析。
有限元分析广泛应用于工程领域，如结构力学、流体动力学、电磁场等领域，用于预测和优化结构的 性能。
有限元分析的基本原理
离散化
将连续的求解域离散化为有限 个小的单元，每个单元具有特
定的形状和属性。
数学建模
根据物理问题的性质，建立每 个单元的数学模型，包括节点 力和位移的关系、能量平衡等。
求解方程
通过建立和求解线性或非线性 方程组，得到每个节点的位移 和应力分布。
PART 05
有限元分析的工程应用实 例
桥梁结构分析
总结词
桥梁结构分析是有限元分析的重要应用之一，通过模拟桥梁在不同载荷下的响应，评估 其安全性和稳定性。
详细描述
桥梁结构分析主要关注桥梁在不同载荷（如车辆、风、地震等）下的应力、应变和位移 分布。通过有限元模型，工程师可以预测桥梁在不同工况下的行为，从而优化设计或进
刚性问题
刚性问题是有限元分析中的一种 特殊问题，主要表现在模型中某 些部分刚度过大，导致分析结果 失真
刚性问题通常出现在大变形或冲 击等动态分析中，由于模型中某 些部分刚度过高，导致变形量被 忽略或被放大。这可能导致分析 结果与实际情况严重不符。
解决方案：为避免刚性问题，可 以采用多种方法进行优化，如采 用更合适的材料模型、调整模型 中的参数设置、采用更精细的网 格等。同时，可以采用多种方法 对分析结果进行验证和校核，以 确保其准确性。
《有限元分析及应用》PPT课件
41
2.3 基本变量的指标表达
指标记法的约定：
自由指标：在每项中只有一个下标出现，如
，
i，j为自由指标，它们可以自由变化；在三维ij 问题
中，分别取为1，2，3；在直角坐标系中，可表示
三个坐标轴x, y, z。
哑指标：在每项中有重复下标出现，如：
，j为哑指标。在三维问题中其变化的范ai围j x为j 1,b2i ,3
有限元方法的思路及发展过程
思路：以计算机为工具，分析任意变形体以获得所有 力学信息，并使得该方法能够普及、简单、高效、方 便，一般人员可以使用。 实现办法：
20
技术路线：
21
发展过程： 如何处理
对象的离散化过程
22
常用单元的形状
.点 (质量)
面 (薄壳, 二维实体,
.. 轴..对称实体.).......
3
有限元法是最重要的工程分析技术之一。 它广泛应用于弹塑性力学、断裂力学、流 体力学、热传导等领域。有限元法是60年 代以来发展起来的新的数值计算方法，是 计算机时代的产物。虽然有限元的概念早 在40年代就有人提出，但由于当时计算机 尚未出现，它并未受到人们的重视。
4
随着计算机技术的发展，有限元法在各个 工程领域中不断得到深入应用，现已遍及 宇航工业、核工业、机电、化工、建筑、 海洋等工业，是机械产品动、静、热特性 分析的重要手段。早在70年代初期就有人 给出结论：有限元法在产品结构设计中的 应用，使机电产品设计产生革命性的变化， 理论设计代替了经验类比设计。
由此得到
考虑 X 0
xyl ym zy n Y xl yxm zxn X
考虑
Z 0 xzl yzm zn Z
应力边界条件
第1章有限元基本理论ppt课件
x dx
li
E i
i
E (ui1ui )
x
x
li
1.8 直杆受自重作用的拉伸问题(续)
❖ 外载荷与结点的平衡方程
EA(uiui1 ) li1
EA(ui1ui ) li
q(li1 li ) 2
q(li1li ) 为第i个结点上承受的外载荷
2
1.8 直杆受自重作用的拉伸问题(续)
❖ 假定将直杆分割成3个单元，每个单元长为a=L/3， 则对结点2，3，4列出的平衡方程为：
单元: 一组节点自由度间相互作用的 数值、矩阵描述（称为刚度或系数 矩阵)。单元有线、面或实体以及二 维或三维的单元等种类。
载荷
有限元模型由一些简单形状的单元组成，单 元之间通过节点连接，并承受一定载荷。
1.6 节点和单元 (续)
信息是通过单元之间的公共节点传递的。
. . 2 nodes ...
. . . 1 node
1.1 有限元分析 (FEA)
有限元分析 是利用数学近似的方法对真实物理
系统（几何和载荷工况）进行模拟。它利用简 单而又相互作用的元素，即单元，用有限数量 的未知量去逼近无限未知量的真实系统。
1.2 有限单元法的基本思想
❖ 将连续的结构离散成有限个单元，并在每一单元中 设定有限个节点，将连续体看作只在节点处相连接 的一组单元的集合体。
I
J
O
N
三维实体结构单元
K UX, UY, UZ
P
M L
J
I
J
K J
O N
K J
三维梁单元 UX, UY, UZ, ROTX, ROTY, ROTZ
三维四边形壳单元 UX, UY, UZ, ROTX, ROTY, ROTZ
有限元及其分析绪论PPT课件
以处理很复杂的连续介质问题，是一种普遍方法。
60年代后期，J.T.Oden 等学者进一步研究了加权残值法与有限元法之间的关系，建立有限元法的计算格式， 并指出有限元法所利用的主要是Galerkin加权残值法，它可以用于即使泛函无法构造或泛函根本不存在的 问题，从而进一步扩大了有限元法的应用领域。
1972年，J.T.Oden 出版了第一本处理非线性连续介质问题的专著 《非线性连续体的有限元法 》。
• 在此期间，O.C.Zienkiewicz、卞学璜、董平等人进一步推动有限元的发展，分别提出了等参单元、杂交 单元的概念。1967年，O.C.Zienk iewic e 和Y.K.Cheung( 张佑启) 出版了第一本有关有限元分析的专著 《连续体和结构的有限元法》，此书是有限元法的名著，后更名为《有限单元法》。
V
•
Galerkin 方法
在Galerkin方法中，选择的加权函数wi为试函数（如取为形函数N，wi=Ni ）
L(x) EIv' ' ' ' p
R(x) EIv' ' ' ' p
L
0 wi(EIv''''' p)dx 0
i 0~n
• 以三角函数为试探函数求ci • 以幂级数为试探函数求ci • 以形函数为试探函数求ci
近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似；利 用与原问题的等效的变分原理（如最小势能原理）建立有限元基本方程（刚度方程）又体现了其明确的物理背 景。
• 厚实的理论基础，数值计算稳定、高效
• 有限元法计算格式的建立既可基于物理概念推得，如刚度法、虚功原理，也可基于纯数学原理推 得，如泛函变分原理、加权残值法。通常直接刚度法、虚功原理用于杆系结构或结构问题的方程 建立；而变分原理设计泛函极值，既适用于简单的结构问题，也适应于更复杂的工程问题（如温 度场问题）。当给定的问题存在经典变分叙述时，则利用变分原理很容易建立这类问题的有限元 方程，如加权残值法。加权残值法由问题的基本微分方程出发而不依赖于泛函，可用于处理一般 问题的有限元方程建立，如流固耦合问题。所以，有限元法不仅具有明确的物理背景，更具有坚 实的数学基础，且数值计算的收敛性、稳定性均可从理论上得到证明，有关这方面的内容可参考 相关资料。
60年代后期，J.T.Oden 等学者进一步研究了加权残值法与有限元法之间的关系，建立有限元法的计算格式， 并指出有限元法所利用的主要是Galerkin加权残值法，它可以用于即使泛函无法构造或泛函根本不存在的 问题，从而进一步扩大了有限元法的应用领域。
1972年，J.T.Oden 出版了第一本处理非线性连续介质问题的专著 《非线性连续体的有限元法 》。
• 在此期间，O.C.Zienkiewicz、卞学璜、董平等人进一步推动有限元的发展，分别提出了等参单元、杂交 单元的概念。1967年，O.C.Zienk iewic e 和Y.K.Cheung( 张佑启) 出版了第一本有关有限元分析的专著 《连续体和结构的有限元法》，此书是有限元法的名著，后更名为《有限单元法》。
V
•
Galerkin 方法
在Galerkin方法中，选择的加权函数wi为试函数（如取为形函数N，wi=Ni ）
L(x) EIv' ' ' ' p
R(x) EIv' ' ' ' p
L
0 wi(EIv''''' p)dx 0
i 0~n
• 以三角函数为试探函数求ci • 以幂级数为试探函数求ci • 以形函数为试探函数求ci
近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似；利 用与原问题的等效的变分原理（如最小势能原理）建立有限元基本方程（刚度方程）又体现了其明确的物理背 景。
• 厚实的理论基础，数值计算稳定、高效
• 有限元法计算格式的建立既可基于物理概念推得，如刚度法、虚功原理，也可基于纯数学原理推 得，如泛函变分原理、加权残值法。通常直接刚度法、虚功原理用于杆系结构或结构问题的方程 建立；而变分原理设计泛函极值，既适用于简单的结构问题，也适应于更复杂的工程问题（如温 度场问题）。当给定的问题存在经典变分叙述时，则利用变分原理很容易建立这类问题的有限元 方程，如加权残值法。加权残值法由问题的基本微分方程出发而不依赖于泛函，可用于处理一般 问题的有限元方程建立，如流固耦合问题。所以，有限元法不仅具有明确的物理背景，更具有坚 实的数学基础，且数值计算的收敛性、稳定性均可从理论上得到证明，有关这方面的内容可参考 相关资料。
《有限元法及其应用》课件
实例
某型战斗机的机翼设计过程中，通过有限元分析，优化了机翼的结构和材料分布，提高了机翼的抗弯和 抗扭能力，同时减小了机翼的气动阻力，为飞机的高性能提供了保障。
汽车碰撞模拟
01
总结词
利用有限元法模拟汽车碰撞过程，评估汽车的安全性能和 改进设计方案。
02 03
详细描述
汽车碰撞是交通事故中最为严重的一种情况，有限元法能 够模拟汽车碰撞过程，对汽车的结构、材料和吸能设计等 进行评估，为汽车的安全性能提供科学依据。同时，通过 模拟不同碰撞条件下的结果，可以为汽车设计提供改进方 案。
通过离散化的方法，将连续的偏微分 方程转化为离散的代数方程组。
刚度矩阵与载荷向量
刚度矩阵
描述了每个单元的刚度关系，反 映了单元之间的相互作用。
载荷向量
描述了作用在每个节点上的外力 。
位移求解与应力分析
位移求解
通过求解离散化的代数方程组，得到每个节点的位移。
应力分析
根据位移求解的结果，通过计算得到每个单元的应力应变状态。
有限元法的应用领域
结构分析
有限元法在结构分析中应用最为广泛，可 以用于分析各种结构的应力、应变、位移
等。
电磁场分析
有限元法可以用于分析电磁场中的电场强 度、磁场强度、电流密度等，如电磁兼容
性分析、天线设计等。
流体动力学
有限元法可以用于模拟流体在各种复杂环 境下的流动行为，如航空航天、船舶、汽 车等领域的流体动力学问题。
应用领域
广泛应用于科学研究和工 程领域，如化学、生物医 学、电磁学等。
FE-SAFE
概述
FE-SAFE是一款用于结构疲劳分析的有限元软件 ，基于有限元方法进行疲劳寿命预测。
特点
某型战斗机的机翼设计过程中，通过有限元分析，优化了机翼的结构和材料分布，提高了机翼的抗弯和 抗扭能力，同时减小了机翼的气动阻力，为飞机的高性能提供了保障。
汽车碰撞模拟
01
总结词
利用有限元法模拟汽车碰撞过程，评估汽车的安全性能和 改进设计方案。
02 03
详细描述
汽车碰撞是交通事故中最为严重的一种情况，有限元法能 够模拟汽车碰撞过程，对汽车的结构、材料和吸能设计等 进行评估，为汽车的安全性能提供科学依据。同时，通过 模拟不同碰撞条件下的结果，可以为汽车设计提供改进方 案。
通过离散化的方法，将连续的偏微分 方程转化为离散的代数方程组。
刚度矩阵与载荷向量
刚度矩阵
描述了每个单元的刚度关系，反 映了单元之间的相互作用。
载荷向量
描述了作用在每个节点上的外力 。
位移求解与应力分析
位移求解
通过求解离散化的代数方程组，得到每个节点的位移。
应力分析
根据位移求解的结果，通过计算得到每个单元的应力应变状态。
有限元法的应用领域
结构分析
有限元法在结构分析中应用最为广泛，可 以用于分析各种结构的应力、应变、位移
等。
电磁场分析
有限元法可以用于分析电磁场中的电场强 度、磁场强度、电流密度等，如电磁兼容
性分析、天线设计等。
流体动力学
有限元法可以用于模拟流体在各种复杂环 境下的流动行为，如航空航天、船舶、汽 车等领域的流体动力学问题。
应用领域
广泛应用于科学研究和工 程领域，如化学、生物医 学、电磁学等。
FE-SAFE
概述
FE-SAFE是一款用于结构疲劳分析的有限元软件 ，基于有限元方法进行疲劳寿命预测。
特点
第一章概述 有限元法基本原理及应用课件
第一章 概述
第一章 概述
有限元法的基本思想 有限元法的特点 有限元法的发展及其应用领域
1.1有限元法的基本思想
2．有限元法是一种应用已知求解未知的思想
在弹性力学领域，已经能用数学偏微分方程将问 题加以表达，但是运用解析方法求解这些方程有时会 很难甚至无法求解。而有限元法是应用人们对事物规 律的已有认识并结合研究对象的各种约束条件，组织 一个运用已知的参量和规律来求解未知问题的有机过 程。
西班牙的Onate E和波兰的Rojek J将DEM 和FEM结合解决地质 力学中的动态分析问题；
瑞典的Birgersson F和英国的Finnveden S针对FEM在频域中的 应用提出了SFEM 。
FEM也从分析比较向优化设计方向发展。印度Mahanty博士用 ANSYS对拖拉机前桥进行优化设计
物体的几何形状可以用大大小小的多种单元进行拼装，所以 有限元法可以分析包括各种特殊结构的复杂结构体。
单元之间材料性质可以有跳跃性的变化，所以能处理许多物 体内部带有间断性的复杂问题，以适应不连续的边界条件和载荷 条件。
三维实体的四面体单元划分
平面问题的四边形单元划分
1.2 有限元法的特点
7．适合计算机的高效计算
20世纪90年代以来，大批FEA系统纷纷向微机移植， 出现了基于各种微机版FEA系统。有限元法向流体力学、 温度场、电传导、磁场、渗流和声场等问题的求解计算 方面发展，并发展到求解一些交叉学科的问题。
1.3.1 有限元法的发展
3．有限元法的研究现状
美国的HeoFanis Strouboulis等人提出用GFEM 解决 分析域内含有大量孔洞特征的问题；比利时的Nguyen Dang Hung 和越南的Tran Thanh Ngoc 提出用HSM解 决实际开裂问题
第一章 概述
有限元法的基本思想 有限元法的特点 有限元法的发展及其应用领域
1.1有限元法的基本思想
2．有限元法是一种应用已知求解未知的思想
在弹性力学领域，已经能用数学偏微分方程将问 题加以表达，但是运用解析方法求解这些方程有时会 很难甚至无法求解。而有限元法是应用人们对事物规 律的已有认识并结合研究对象的各种约束条件，组织 一个运用已知的参量和规律来求解未知问题的有机过 程。
西班牙的Onate E和波兰的Rojek J将DEM 和FEM结合解决地质 力学中的动态分析问题；
瑞典的Birgersson F和英国的Finnveden S针对FEM在频域中的 应用提出了SFEM 。
FEM也从分析比较向优化设计方向发展。印度Mahanty博士用 ANSYS对拖拉机前桥进行优化设计
物体的几何形状可以用大大小小的多种单元进行拼装，所以 有限元法可以分析包括各种特殊结构的复杂结构体。
单元之间材料性质可以有跳跃性的变化，所以能处理许多物 体内部带有间断性的复杂问题，以适应不连续的边界条件和载荷 条件。
三维实体的四面体单元划分
平面问题的四边形单元划分
1.2 有限元法的特点
7．适合计算机的高效计算
20世纪90年代以来，大批FEA系统纷纷向微机移植， 出现了基于各种微机版FEA系统。有限元法向流体力学、 温度场、电传导、磁场、渗流和声场等问题的求解计算 方面发展，并发展到求解一些交叉学科的问题。
1.3.1 有限元法的发展
3．有限元法的研究现状
美国的HeoFanis Strouboulis等人提出用GFEM 解决 分析域内含有大量孔洞特征的问题；比利时的Nguyen Dang Hung 和越南的Tran Thanh Ngoc 提出用HSM解 决实际开裂问题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
0
0
k2u2 k2u3 k3u3 k3u4
k3u3 k3u4 k4u4 k4u5 0
k4u4 k4u5 P
写成矩阵的形式为
k1
=
k1 k1 k2 k2 0 0
k1 k1 0 0 0
0 k2 k 2 k3 k3 0
k1 k1 k2 k2 0 0
有限元方法与ANSYS简介
有限元方法是用于求解工程中各类问题的数值方法，应 力分析中稳态的、瞬态的、线性的或非线性的问题以及热传导、 流体流动和电磁学中的问题都可以用有限元方法进行分析解决。 现代有限元方法的20世纪早期开始，20世纪50年代，boeing公司 采用三角元对机翼进行建模，推动了有限元方法的应用。到20 世纪60年代，人们接受了“有限元”这个词。 ANSYS是一个通用的有限元计算机程序，其代码长度超 过10万行。应用ANSYS可以进行静态、动态、热传导、流体流 动和电磁学等分析。在过去的20多年里，ANSYS是主要的有限 元分析程序。现在ANSYS被广泛应用在如航天、汽车、电子、 核科学等领域。
第一章 概述
有限元方法是广泛用于解决应力分析、热传 递、电磁场和流体力学等工程问题的数值方 法。
本章的内容
（1）工程问题 （2）数值方法 （3）有限元方法与ANSYS简介 （4）有限元方法的基本步骤 （5）直接公式法 （6）最小总势能公式 （7）加权余数法 （8）结果的验证 （9）理解问题
工程问题
0
R1 0 0 0 0
0 k2 k 2 k3 k3 0
0 k3 k3 k 4 k4
0 u1 0 0 u 2 0 0 u3 0 k4 u 4 0 k4 P u5
上式表示成
R K u F
上式中， R为反作用力，K为刚度矩阵，u为位移矩阵，F为负 荷矩阵。由于杆的顶端是固定的，其位移为零，
u1 0
上面的方程简化为
有限元方法的基本步骤
预处理阶段 1 建立求解域并将之离散化成有限元，即将问题分解成节 点和单元。 2 假设代表单元物理行为的形函数，即假设代表单元解的 近似连续函数。 3 对单元建立方程。 4 将单元组合成总体的问题，构造总体刚度矩阵。 5 应用边界条件、初值条件和负荷。 求解阶段 6 求解线性或非线性的方程组，以得到节点的值，例如得 到不同节点的位移量或热传递问题中不同节点的温度值。 后处理阶段 7 得到其他重要的信息。如主应力、热量值等。
1.5 直接公式法 例1.1 杆的一端固定，另一端承受负荷P，杆的厚度t，长 度L。杆的弹性模量用E表示。假设应用的负荷比杆的 重量大的多。
预处理阶段
将问题域离散成有限的单元， 首先将问题分解成节点和单元 。我们用5节点4单元的模型代 替杆，如图所示，分的越细越 精确。杆的模型中有四个独立 的分段，每个分段都有统一的 横截面积。
数值方法
在许多实际工程问题中，由于微分方程的复杂性或边界 条件和初始条件的难以确定性，得不到系统的精确解。为此我 们借助数值方法来近似。解析解表明了系统在任何点上的精确 行为，而数值解只在称为节点的离散点上近似于解析解。 数值法的第一步都是离散化。这一过程将系统分成一些 单元和节点，然后对每一单元或节点建立代数方程组。 这种方法假设代表每个单元的近似函数是连续的。假设 单元间的边界是连续的，通过组合各单元的解产生系统的完全 解。
工程问题一般是物理问题的数学模型。数学模型是带有边 界条件和初始问题的微分方程，微分方程是通过对系统应用自 然的基本定律和原理推导出来的，如波动方程等。 这些微分方程代表了某种物理规律或平衡。（数理方程） 在可能的情况下，由给定的条件求解微分方程可以得到系 统的精确行为。在任何给定的工程问题中，存在两种影响系统 行为的参数。一种表示给定系统自然行为的参数，例如弹性模 量、热传导因子和粘度等。另一种是系统存在产生扰动的参数， 如外力、力矩、介质的温度差和流体的压力差等。 在有限元建模中，理解参数在刚度或传导矩阵以及负荷矩 阵中的作用是非常重要的。系统特性总是在刚度矩阵、传导矩 阵中得到体现，而扰动参数总是出现在负荷矩阵中。
这相当于弹簧。我们认为杆由四个弹性系数不同的弹簧组成， 对每个弹簧来说，
keq
Ai 1 Ai E 2l
Ai 和 Ai 1 分别是i 和i 1 处的节点的横截面积，l是单元的长度。
静力学要求每个节点上的应力和为零，所以 节点1： 节点2： 节点3： 节点4： 节点5：
R1 k1 u2 u1 0
k1 u2 u1 k2 u3 u2 0
k2 u3 u2 k3 u4 u3 0
k3 u4 u3 k4 u5 u4 0
k4 u5 u4 P 0
重组方程组，得到
k1u1 k1u2
R1
k1u1 k1u2 k2u2 k2u3
假设近似单元行为的近似解
考虑一个带有统一横截面A的实体的偏转量，横截面的 长度为l，承受的外力为F，如图所示。 实体的平均应力为
F A
实体的平均应变为
l l
应变与应力的关系为
E
这里E是弹性模量。
这样有
AE keq l
AE F Biblioteka l l 令，则
F keq l
0 0 k3 k3 k 4 k4
k1 0 0 0
在负荷矩阵中，将反作用力和负荷区分开来是很重要的，于是
0 u1 R 1 0 u 2 0 0 u3 = 0 0 k4 u 4 P k4 u 5