有限元方程的求解方法
有限元法介绍
通俗地说，有限元法就是一种计算机模拟技术，使人们能够在计算机上用软件模拟一个工程问题的发生过程而无需把东西真的做出来。
这项技术带来的好处就是，在图纸设计阶段就能够让人们在计算机上观察到设计出的产品将来在使用中可能会出现什么问题，不用把样机做出来在实验中检验会出现什么问题，可以有效降低产品开发的成本，缩短产品设计的周期。
有限元法也叫有限单元法（finite element m ethod, FEM），是随着电子计算机的发展而迅速发展起来的一种弹性力学问题的数值求解方法。
五十年代初，它首先应用于连续体力学领域—飞机结构静、动态特性分析中，用以求得结构的变形、应力、固有频率以及振型。
由于这种方法的有效性，有限单元法的应用已从线性问题扩展到非线性问题，分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料，从连续体扩展到非连续体。
有限元法最初的思想是把一个大的结构划分为有限个称为单元的小区域，在每一个小区域里，假定结构的变形和应力都是简单的，小区域内的变形和应力都容易通过计算机求解出来，进而可以获得整个结构的变形和应力。
事实上，当划分的区域足够小，每个区域内的变形和应力总是趋于简单，计算的结果也就越接近真实情况。
理论上可以证明，当单元数目足够多时，有限单元解将收敛于问题的精确解，但是计算量相应增大。
为此，实际工作中总是要在计算量和计算精度之间找到一个平衡点。
有限元法中的相邻的小区域通过边界上的结点联接起来，可以用一个简单的插值函数描述每个小区域内的变形和应力，求解过程只需要计算出结点处的应力或者变形，非结点处的应力或者变形是通过函数插值获得的，换句话说，有限元法并不求解区域内任意一点的变形或者应力。
大多数有限元程序都是以结点位移作为基本变量，求出结点位移后再计算单元内的应力，这种方法称为位移法。
有限元法本质上是一种微分方程的数值求解方法，认识到这一点以后，从70年代开始，有限元法的应用领域逐渐从固体力学领域扩展到其它需要求解微分方程的领域，如流体力学、传热学、电磁学、声学等。
有限元综述.(优选)
有限元综述蔡璟、吕丹丹、李川摘要：有限元法（Finite Element Method）是一种高效能、常用的数值计算方法。
1965年“有限元”这个名词第一次出现，经历了三十多年的发展历史，理论和算法都已经日趋完善。
如今，有限元在工程上得到广泛应用。
本文首先介绍了有限元的研究背景和意义，其次从它的诞生、主要特点以及解题步骤三方面阐述相关概念，再讨论传统有限元算法及优化算法、有限元与其他算法结合得到的混合算法两个方面来分类阐述各自的研究现状与特点，最后总结有限元算法的应用以及发展趋势。
关键词：有限元法，FEM，经典算法，优化算法，网格优化，Herrmann算法，时域有限元，混合算法，矩量法，时域有限差分，应用研究，边界元法，光滑粒子法，发展趋势前言有限元法（Finite Element Method）是一种高效能、常用的数值计算方法，其基本思想是由解给定的泊松方程化为求解泛函的极值问题。
有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。
自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，解决了物理场应用中的限制。
经历几十年的发展，有限元法已经被广泛用于各个领域。
1.研究背景和意义有限元法的思想首先由 R. Courant 在 1943 年提出，十九世纪六十年代数值分析科学家认识了有限元基本思想，建立了有限元方法的数学基础。
其中，我国数学家冯康独立地提出了有限元方法，将其命名为“基于变分原理的差分格式”，对有限元方法的创始及奠基工作做出了重要贡献。
以变分原理为基础建立起来的有限元法,因其理论依据的普遍性,不仅广泛地被应用于各种结构工程,而且作为一种声誉很高的数值分析方法已被普遍推广并成功地用来解决其他工程领域中的问题,例如热传导!渗流!流体力学、空气动力学、土壤力学、机械零件强度分析、电磁场工程问题等等。
第五章 有限元法-9-有限元方程组的求解
利用变分原理和离散化方法建立有限元矩阵方程后， 须求解以结点值为未知数的矩阵方程。 其方程写为： Ax b
（2.42）
式中系数矩阵A是一个n×n方阵，x是待求解的未知量，b表 示已知向量。
为精确描述实际问题，系数矩阵的维数（对应离散剖 分的结点值未知量个数）往往非常大，有庞大的计算 机内存需求和过长的计算时间。
除存储量降低外，有限元矩阵的特殊性质也能减少计算时间。大 量的零矩阵元素不需产生，加上适当设计算法，它们在解过程中 的运算也可避免。 正是这一为矩量法等积分方程方法所不具备的特殊性质，使得有 限元方法对分析电大尺寸问题时更有吸引力。
下面先介绍矩阵方程的解法，在此基础上然后 介绍Ansoft HFSS为在一定精度的要求上最大 限度的提高效率而设计的自适应迭代算法。
Lanczos法是有效的求解带状稀疏矩阵的本征值问题 的方法，Ansoft HFSS可以在solver中找到。
4 Ansoft HFSS的自适应迭代算法
矩阵方程的求解复杂度与有限元的剖分密度即未知数数目有很大 的关系，未知数数目越多，求解所需的时间越长。 然而，另一方面，有限元方法求解的精度与也随着未知数数目的 增加而更加准确。 因此，有限元方法的求解时间与准确度是一对矛盾。 为了在短的时间内取得越大的精度，Ansoft HFSS采取了自适应 迭代算法，如图2.5所示。 该算法一开始先选用较粗的剖分，采用前述的方法求解，然后看 其进度是否满足要求。 如不满足，进一步细化剖分，再次进行求解，直至达到给定的精 度。
如果矩阵可以分解为 A＝LU
（2.43）
其中，L是一个下三角矩阵，U是一个上三角矩阵。
有限元分析考试
有限元复习资料工程问题的研究对象：离散体：材料力学中的连续梁、建筑结构框架和桁架结构。
我们把这类问题，称为离散系统。
离散问题是可解的连续体：通常可以建立它们应遵循的基本方程，即微分方程和相应的边界条件。
通常只能得到少数简单问题的精确解答。
对于许多实际的工程问题，还无法给出精确的解答。
有限元方法的两大应用:科学计算或数值模拟数字设计或虚拟仿真(CAD CAM CAE)有限元法的实质：将复杂的连续体划分为有限个简单的单元体，化无限自由度问题为有限自由度问题。
将连续场函数的（偏）微分方程的求解问题转化为有限个参数的代数方程组的求解问题。
有限元法的基本思路：将连续系统分割成有限个分区或单元，对每个单元提出一个近似解，再将所有单元按标准方法组合成一个与原有系统近似的系统有限元分析的基本步骤：1）建立研究对象的近似模型。
2）将研究对象分割成有限数量的单元(结构离散化)3）用标准方法对每个单元提出一个近似解（单元分析）（1）选择位移模式（2）建立单元刚度矩阵（3）计算等效节点力4）将所有单元按标准方法组合成一个与原有系统近似的系统（单元集成）5）用数值方法求解这个近似系统。
（选择合适计算方法计算节点位移、应力、应变）6）计算结果处理与结果验证（如何显示、分析数据并找到有用的结论）弹性力学的基本假设：（1）连续性（2）均匀性（3）各向同性（4）完全弹性符合（1）-（4）假定的称为理想弹性体。
（5）小变形假定满足以上五个基本假设的弹性力学称为线弹性力学。
弹性力学中的基本量：1）外力（体力面力集中力）弹性体受外力以后，其内部将产生应力。
三种外力应力（正应力剪应力）应力成对出现，坐标面上的应力以正面正向，负面负向为正。
六个应力分量剪应力互等定律：作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是互等的。
(大小相等，正负号也相同)。
弹性体在受外力以后，还将发生变形。
物体的变形状态，一般有两种方式来描述：1、给出各点的位移；（运动）三个位移分量2、给出各微元体的变形（应变）六个应变分量在弹性体内部，三类基本方程：根据微分体上力的平衡条件（法向应力和剪应力即内力与其体力平衡），，建立平衡微分方程；根据微分线段上应变－位移的几何条件，建立几何方程；根据应力－应变间的物理条件，建立物理方程当沿X轴方向的两个对面受有均匀布的正应力时，弹性体在X方向的伸长还伴随有侧收缩，即在y和Z方向的单位缩短。
有限元方法
有限元方法求解微分方程，特别是椭圆型边值问题的一种离散化方法，其基础是变分原理和剖分逼近。
有限元方法是传统的里茨－加廖金方法的发展，并融会了差分法的优点，处理上统一，适应能力强，已广泛应用于科学与工程中庞大复杂的计算问题。
作为有限元方法出发点的变分原理，是表达物理基本定律的一种普遍形式。
其表述可概括如下：给出一个依赖物理状态v的变量J(v)（v是函数,J(v)在数学上称为泛函），同时给出J(v)的容许函数集V,即一切可能的物理状态,则真实的状态是V中使J(v)达到极小值的函数。
剖分逼近是有限元离散化的手段，把问题的整体（即求解域）剖分为有限个基本块，称为"单元"，然后通过单元上的插值逼近，得到一个结构简单的函数集，称为"有限元空间",它一般是容许函数集V的子集或有某种联系。
有限元方法就是在这个有限元空间中寻找J(v)的极小解作为近似解。
典型问题为具体说明有限元方法，讨论二维有界域Ω上的椭圆型方程, (1)变系数β表示介质不均匀。
物理学中许多平衡态或定常态问题都可归结为这个典型方程。
与方程(1)相配的有如下三类边界条件：第一类：；第二类：；第三类：。
这里的φ、g及α均为定义在边界дΩ上的已知函数，表示外法向导数，第二类边界条件是第三类当α=0时的特例。
为说明有限元方法能统一处理复杂的情况，假定讨论的问题是混合边值，并且介质有间断,即дΩ分成Г0和Г1两部分，分别有边界条件, (2),(3)β（x,y)有间断线，把Ω分为Ω-,Ω+两部分，在间断线上微分方程(1)无定义，而代之以接触条件, (4)及表示间断线上分别指向Ω+及Ω-的法向导数。
变分原理与微分方程(1)及附加条件(2)、(3)、(4)的边值问题相对应的是物理学中的极小能量原理。
构造"能量积分"并取J(v)的容许函数集V为一切满足边界条件(2)且一阶偏导数平方可积的函数,则使J(v)达到极小值的u,即,(6)也必满足方程(1)及(2)、(3)、(4)。
4.5.14.5平面问题有限元分析步骤及计算实例
K
88
K 12 11 K21 1
K 12 31
K41 2
K22 1 K32 1
K 12 33
K43 2
K
44
2
由于[Krs]=[Ksr]T，又单元1和单元2的节点号按1、2、
3对应3、4、1，则可得：
K11 1
K33 2
3E 16
3 0
0 1
K21 1 K43 2
K12 1
3E 8
3 1 0
0 0 1
3 1 1
1 3 1
0 0 1
013
q/E 0
q/E 0
3E 8
8q
0 /(3E) 0
0 q1
0
0
单元应力可看作是单元形心处的应力值。
7）引入约束条件，修改刚度方程并求解
根据约束条件：u1 =v1=0；v2=0；u4=0和等效节点力列
阵：F 0 0 0 0 0 q / 2 0 q / 2T
五. 边界条件的处理及整体刚度矩阵的修正 整体刚度矩阵的奇异性可以通过引入边界约束条件来排除弹性体的
刚体位移，以达到求解的目的。
（两种）方法 “化1置0法”
“乘大数法”
⑴修改后的总刚为非奇异，对应的总体平衡方程可求解； ⑵如果已知位移不等于0，采用第二种方法，固定约束用 第一种方法。 ※求解可以采用解方程组的任何一种方法。（高斯消去法 常用），可借用一些计算机软件：如Matlab，Excel等。
所以 q / E0 0 1/ 3 0 1/ 3 1 0 1T
习题和思考题
• 4.1三角形常应变单元的特点？ • 4.2平面问题有限元法的基本思想和解题步骤。 • 4.3简述形函数的概念和性质。 • 4.4平面问题整体刚度矩阵的推导过程。 • 4.5矩形单元的特点？ • 4.6有限元方法解的收敛准则。
c3d8有限元单元方程推导过程
有限元单元方程推导过程1.引言有限元分析是一种数值计算方法，用于求解结构力学、流体动力学等领域的物理问题。
在有限元分析中，有限元单元是构成整个有限元模型的基本单元，通过推导有限元单元的方程，可以实现对结构或系统的精确分析和计算。
本文将从有限元方法的基本原理出发，详细介绍有限元单元方程的推导过程。
2.有限元方法基本原理有限元方法是将连续的物理问题离散化，转化为有限个代表性元素的集合，通过对每个元素施加适当的边界条件和力学方程，最终得到整个系统的解。
有限元方法通过有限元单元之间的相互作用，从而模拟整个系统的行为。
3.有限元单元的概念有限元单元是有限元模型中最小的离散单元，它是对实际的结构或系统进行离散化的结果。
不同的物理问题和结构，可以采用不同类型的有限元单元进行离散化，如梁单元、壳单元、板单元等。
4.有限元单元方程的一般形式有限元单元方程的一般形式可以表示为：\[K_{e}U_{e}=F_{e}\]其中\(K_{e}\)为有限元单元的刚度矩阵，\(U_{e}\)为有限元单元的位移矢量，\(F_{e}\)为有限元单元的荷载矢量。
5.有限元单元方程推导的基本步骤有限元单元方程的推导主要包括以下几个基本步骤：5.1 单元刚度矩阵的推导首先需要根据有限元单元的几何形状和材料性质，推导出单元刚度矩阵。
单元刚度矩阵可以通过对单元内部的应变能量或者应力-应变关系进行积分得到。
5.2 单元位移矢量的表示在推导单元方程过程中，需要选择合适的位移矢量表示方式，可以采用基函数展开的方法，将位移矢量表示为一组未知系数乘以基函数的线性组合形式。
5.3 单元荷载矢量的求解单元荷载矢量是由外部施加的荷载和边界条件共同决定的，在推导单元方程的过程中需要将这些荷载转化为局部坐标系下的形式，并利用位移矢量的表示方式，将荷载矢量表达为位移矢量和未知系数的线性组合。
5.4 单元方程的组装需要将单元刚度矩阵、位移矢量和荷载矢量组装成完整的单元方程，可以通过坐标变换或者有限元单元之间的关系对单元方程进行组装。
偏微分方程的有限元法求解
16.901讲义笔记一维有限％首先，我们考虑•个比上一节稍微复杂点的问题; 豎二f(X),卫冲，V(O) = O.V(L)=O在这里，f(X)是)C的般函数，我们来看•个特别的情形：f(x)=x(L-x),此时，方程的梏确解如F：有限元方法利用加权残差的方法■其中:(1)设va)=£«Ma), v()()是我们对v(x)的近似，省为未知常数9 V|(x)是用户选择的歯数，即形状朗数：(2)定义N个加权残差LRj = p^(x)R(V)dx • j = l-> N to其中，RV)二器・f为绒差凹⑴足“用户”选择的加权函数，即权函数:(3)令加权残并为冬•町以确定⑷的值，即求耳使得对所fi 1=I->N, Rj=Oe令限元方法( )是加权残若法的一种，下血看看我们是如何用它来解决问题的。
一维有限元方法有限元方法（〉扌野个连续区域离散化-系列小单尤，这些单元与有限差分法（）或有限体积法（）产牛的网格完全相同，而佼之前两者主耍的优点在于：能够容易地把握单元的变化范囤。
对于我们讨论的一维问题,可以将区域（数轴〉离散化为如下图所示:这里，叫三单•元的个数。
我们还会用別下血i些定义：个三角划分；尽管令限元法对于一维，二维，三维甚至高细问题都是仃效的，们我们还是要谈及区域离散化的一种方浓，即三角划分。
4 T定义为第I个单元所在的区域。
对于_维问题,这表明，TS-个满足片心的X的集合。
接卜来耍确定的是毎个单兀该用什么样的函数，典型的函数形式就是用从一个单元到卜一个单兀保持解连续的多项式。
例如：一个线性有限元如卜團;i示:在毎个单元内的函数是线形的，在毎两个单元的交点处足连续的。
对于专门诜择的满足线件变化的形状函数,右估计残差时有一个很明显的问题：回忆前曲的内容，RV)二器一f,它在一个单冗里等于什么呢?因为函数是线性的,所以器=0,则有:R(V)=f ,即R(V)与无关。
冋时，满足线性变化的形状函数似乎也是一个好的近似，我们举-个例子来说明。
第3章 有限元分析的数学求解原理-三大步骤
U x x y y z z xy xy yz yz zx zx dV
X u Y v Z w dV X u Y v Z w d W
V V
用 * 表示；引起的虚 应变分量用 * 表示
j Vj
Ui
i Vi
0 X
y
¼ 1-9 Í
ui* * vi wi* * * u j ， v* j w*j
x* * y * z * * xy *yz * 18 zx
19
7.间接解法：最小势能原理
20
最小势能原理
W U 0
最小势能原理就是说当一个体系的势能最小时，系统会处于稳定 平衡状态。或者说在所有几何可能位移中，真实位移使得总势能取最小值
0 表明在满足位移边界条件的所有可能位移 最小势能原理： 中，实际发生的位移使弹性体的势能最小。即对于稳定平衡状态，实 际发生的位移使弹性体总势能取极小值。显然，最小势能原理与虚功 原理完全等价。 n m
虚功原理的矩阵表示
在虚位移发生时，外力在虚位移上的虚功是：
* 式中
U i u i* V i v i* W i w i* U j u *j V j v *j W j w *j
* 是 的转置矩阵。
T
*
F
T
同样，在虚位移发生时，在弹性体单位体积内，应力在虚应变上的虚 功是： * * * * * * * T x x y y z z xy xy yz yz zx zx
27
⑴解析法
有限元分析法概述
第十一章 有限元分析方法概述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小段，每一小段采用等截面直杆近似，等截面直杆的横截面面积为相应的变截面杆横截面面积的平均面积表示，每一小段称为一个单元，小段之间通过节点连接起来。
2+有限元法的基本概念与求解方法
a11x1+ a12x2 +…+ a1nxn =b1 a21x1+ a22x2 +…+ a2nxn =b2 …………
an1x1+ an2x2 +…+ annxn =bn 当其系数行列式 A 不等于零时
A
上述的方程组有唯一解：
x j
j
A
（a11 ann系数） （j＝１，２…n）
其中 A 是将Ａ中第j列元素替换为右端项而得到的行列式 j
综上所述，三角形常应变单元属于协调元
位移函数与形函数
面积坐标
面积坐标是利用三角形的面积关系表示三角形单元任一点位置 的一种方法。优点：简明，方便。
位移函数与形函数
面积坐标
对于图中三角形单元任一点P（x , y）可用下三个比值来确定：
l i ,l j ,l m
i
j
m
Li , Lj , Lm 称为P点的面积坐标，显然面积坐标具有以下性质： 性质１. Li + Lj + Lm =1 （ i+ j + m = ）
性质２. 平行于三角形jm边的直线上所有点其Li相同 AB变化时， hi’不变, 故i不变， Li 不变
位移函数与形函数
面积坐标
性质３. Li =1 Lj =0 Lm =0 （i） Li =0 Lj =1 Lm =0 （j） Li =0 Lj =0 Lm =1 （m）
性质４. Li = Lj = Lm=1/3 在三角形形心处 面积坐标与形函数的关系：
结构离散化： 1）网格划分 将结构划分为有限个单元； 2）载荷移置 将作用在结构上的非节点载荷等效地移置为节点载荷； 3）简化约束 把结构边界上的约束，用适当的节点约束代替。
有限元法
称为有限元刚度矩阵，但 不能直接求解，需要消去
1行、1列。
2. 一维有限元法
由边界条件对整个问题的代数方程组消元：
由问题的边界条件，第5 个节点电位为0.5V，已知，故消去该节点的方程：5 行5列。必有这一步，实际上原K矩阵行列式的值为0，本质上是找参考电位
5 5
1 0.1
1. 有限元法
上一讲，利用加权余数法和变分法将偏微分方程转化为代数方程组求解
KC f b
kij k ji j i d
f
＝
j
jq
d
bj
2
jh
d
通过尝试函数的 选取，近似解满 足1类边界条件，
该矩阵方程包括系数矩阵、激励源矩阵和边界矩阵，而计算这些矩阵的元素 时，常常用到分部积分法。如果为了计算精度而选取很多个尝试函数，那么 计算这些为数众多的分部积分既十分复杂又很费时间，并且很难用计算机进 行数值计算。
2. 一维有限元法
局部系数矩阵的计算
k
e ij
kkieie,1i,i
ke i ,i 1
ke i 1,i 1
f
e i
f
e i
f
e i 1
bie
bbieie1
kiej e j i de
2. 一维有限元法
本例，场域分割成4个单元，5个节点， 求场域内电势分布，转化为求5个节点的 电位即可。 场域内其它点（各单元内）的电位，由5 个节点电位来插值表示。(一阶插值、高 阶插值)
对一 维场域来说，单元就是一个线段； 对二维场域，有限元单元形状可为二角 形、矩形等，单元形状对有限元的简化 有影响，通常为三角形
第一章 有限元法基本原理
称为弹性体的虚功方程。
§1.5 单元刚度矩阵
假设杆单元L： ui Ui
l
其中：A, E , l 为参数。
u j U j
杆单元应力-应变关系为：
则：
U j E E u j ui
A
l
EA
U j l u j ui
由力的平衡条件： U i U j 0
则：
Ui
EA l
uj
§1.3 单元应力和应变
位移函数→几何方程（→应变）→物理方程（→应力） 杆单元的几何方程为：
x
du dx
x
d Ne
dx
d dx
1
x l
d dx
x l
e
11
l
1 e
简化为： B e
1.8
B 11 1
1.9
l
B 为几何矩阵。
杆单元的物理方程为：
或：
x E x
D
D 为弹性矩阵 。
参考书目：
❖ 《有限元法及其在锻压工程中的应用》吕丽萍主编， 西北工业大学出版社
❖ 《弹性和塑性力学中的有限单元法》丁皓江等主编， 机械工业出版社
❖ 《有限元分析的基本方法及工程应用》周昌玉、贺 小华 编著，化学工业出版社
§1.2 位移函数与形状函数
1、坐标系
以杆单元为例：
ui U i
Y
vi Vi
1960年， Clough处理 平面弹性问题， 第一次提出“ 有限单元法”。
1963—1964, Besseling, Melosh,Jones 证明有限元法是 基于变分原理 的里兹法， 确认了有限元 法是处理连续 介质问题的 一种普遍方法。
变分法建立 有限元方程 与经典里兹 法的主要 区别
有限元方法基本原理
有限元方法基本原理有限元方法（Finite Element Method, FEM）是一种数值计算方法，主要用于求解偏微分方程的数值解。
它最早由Courant、Bubnov和Galerkin等人在20世纪50年代提出，并在以后的几十年中得到了广泛的发展和应用。
有限元方法的基本原理是将要求解的区域分割成若干个小的子区域，通常称为有限元，每个有限元内部的物理量可以用一个简单的数学表达式来表示。
然后，通过在有限元之间建立连续性条件，将整个问题转化为一组代数方程，进而得到数值解。
有限元方法的基本步骤包括：建立有限元模型、离散化、建立代数方程、求解代数方程和后处理。
下面将详细介绍每个步骤的具体内容。
第一步，建立有限元模型。
该步骤主要是对要求解的问题进行数学建模，包括选择适当的坐标系、定义物理量和约束条件等。
通常，物理问题可以通过连续介质假设，将其离散化为一组小的有限元。
第二步，离散化。
将要求解的区域划分为有限个小的子区域，通常称为有限元。
常见的有限元形状包括三角形、四边形和六面体等。
有限元的选择通常是根据问题的几何形状和物理条件来确定的。
第三步，建立代数方程。
有限元方法的核心是建立代数方程，用于描述物理问题在离散点上的数值解。
代数方程通常是通过施加适当的数学形式和边界条件来建立的。
建立代数方程的基本思想是使用一组试验函数来近似描述有限元内部的解。
通常采用Galerkin方法，即在离散点上进行加权残差积分，使得残差的加权平均为零。
第四步，求解代数方程。
一旦代数方程建立完成，就可以使用数值方法求解这组代数方程。
常见的求解方法包括直接法和迭代法等。
直接法适用于方程较小的情况，而迭代法适用于方程较大的情况。
常见的迭代法有Jacobi迭代法、Gauss-Seidel迭代法和共轭梯度法等。
第五步，后处理。
求解代数方程后，需要对结果进行后处理和分析。
后处理包括计算和显示物理量、绘制图形以及进行误差估计等。
通过后处理，可以对模型进行验证，并对结果进行解释和解释。
有限元法的基本步骤
有限元法的基本步骤有限元法是一种用于求解较为复杂的实际工程问题的数值分析方法。
它将一个连续的物体或系统划分为许多小的单元，然后通过建立在这些单元上的数学方程来模拟和求解实际问题。
在这篇文章中，我们将探讨有限元法的基本步骤，并深入讨论其原理和应用。
1. 确定问题的边界和几何形状在使用有限元法求解实际问题之前，需要先确定问题的边界和几何形状。
通常情况下，问题的边界需要定义为固定边界或自由边界，以便在数学模型中进行处理。
问题的几何形状也需要被建模和描述，这样才能得到准确的计算结果。
2. 划分网格划分网格是有限元法中非常重要的一步。
网格划分是将问题的几何形状划分为一系列小的单元。
这些小单元称为有限元，它们可以是三角形、四边形或其他形状。
网格的划分需要根据问题的几何形状和求解精度来确定，并且需要保证各个有限元之间具有充分的连续性和相互联系，以确保模拟结果的准确性和可靠性。
3. 建立数学模型和方程在确定问题的边界和划分网格之后，下一步是建立与物理现象相关的数学模型和方程。
根据问题的具体情况，可以使用不同类型的方程，如静力学方程、热传导方程、流体力学方程等。
这些方程将物理现象转化为数学表达式，并可以通过有限元法进行求解。
4. 应用边界条件在建立数学模型和方程之后，需要应用边界条件。
边界条件可以是物体的固定边界条件，如固定端或自由端；也可以是物体的外部边界条件，如外力、温度等。
边界条件的正确应用对于求解实际问题非常重要，它们将影响模拟结果的准确性和可靠性。
5. 求解数学方程一旦建立了数学模型、划分网格并应用了边界条件，下一步就是使用数值方法求解数学方程。
有限元法将整个问题转化为一个求解代数方程组的问题，并通过迭代方法求解。
求解过程中需要根据初始条件和边界条件进行迭代计算，直到得到收敛的解。
通过以上的基本步骤，我们可以使用有限元法对复杂的实际工程问题进行数值求解。
有限元法的优点在于可以模拟各种不同的物理现象，并且可以对复杂的几何形状进行建模和求解。
有限元法的基本思想及计算步骤
用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1)连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见 的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单 元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单 元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约 束的结点上根据实际情况设置约束条件。 2)单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三 角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个 结点有两个位移分量u，v和两个结点力分量Fx，Fy。三个结点共六个结点位移分量可用列阵(δ)e 表示： ,δ-e=*ui vi uj vj um vm+T 同样，可把作用于结点处的六个结点力用列阵{F}e表示： {F}e=[Fix Fiy Fjx Fjy Fmx Fmy]T 应用弹性力学理论和虚功原理可得出结点位移与结点力之间的关系 ,F-e=*k+e,δ-e (1)式中 [k]e——单元刚度矩阵。
有限元语言及编译器finiteelementlanguagecompiler以下简称felac是中国科学院数学与系统科是具有国际独创性的有限元计算软件是pfepg系列软件三十年成果1983年2013年的总结与提升有限元语言语法比pfepg更加简练更加灵活功能更加强大
有限元法的基本思想及计算步骤
有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散 化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性 体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种 联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能 通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结 点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各 个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作 为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单 元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之 间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分 量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性 要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改 进，近似解最终将收敛于精确解。
有限元计算的流程
有限元计算的流程下载温馨提示:该文档是我店铺精心编制而成，希望大家下载以后，能够帮助大家解决实际的问题。
文档下载后可定制随意修改，请根据实际需要进行相应的调整和使用，谢谢!并且，本店铺为大家提供各种各样类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，如想了解不同资料格式和写法，敬请关注!Download tips: This document is carefully compiled by theeditor. I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classic articles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!有限元计算的流程一、问题定义阶段在进行有限元计算之前，首先需要明确要解决的问题。