有限元法的基本原理

第二章有限元法的基本原理有限元法吸取了有限差分法中的离散处理内核，又继承了变分计算中选择试探函数并对 区域积分的合理方法有限元法的理论基础是加权余量法和变分原理，因此这里首先介绍加 权余量法和变分原理2.1等效积分形式与加权余量法加权余量法的原理是基于微分方程等效积分的提法，同时它也是求解线性和非线性微分 方程近似解的一种有效方法在有限元分析中，加权余量法可以被用于建立有限元方程，但 加权余量法本身又是一种独立的数值求解方法2.1.1微分方程的等效积分形式工程或物理学中的许多问题，通常是以未知场函数应满足的微分方程和边界条件的形式 提出来的，可以一般地表示为未知函数u应满足微分方程组, A1（u）、A（u） = A2（u） = 0 （在内） （2-1）I : J域可以是体积域、面积域等，如图2-1所示同时未知函数u还应满足边界条件/ Bi（u）、B（u） = B2（u） = 0 （在 r 内） （2-2）I : J要求解的未知函数u可以是标量场（例如压力或温度），也可以是几个变量组成的向量 场（例如位移、应变、应力等）A，B是表示对于独立变量（例如空间坐标、时间坐标等） 的微分算子微分方程数目应和未知场函数的数目相对应，因此，上述微分方程可以是单个 的方程，也可以是一组方程。

所以在以上两式中采用了矩阵形式以二维稳态的热传导方程为例，其控制方程和定解条件如下：d ^6 QA（©） = f（k ） + f（k ） + q = 0 （在内） （2-3）ox ox o~y oyc ± J-八 = 0BS） = \ 80 - 口（在L上）0（在「上）q（2-4）这里0表示温度（在渗流问题中对应压力）;k是流度或热传导系数（在渗流问题中对K 虱-q=0应流度K /目）；矿和q是边界上温度和热流的给定值（在渗流问题中分别对应边界上的压 力和边界上的流速）；n是有关边界r的外法线方向；q是源密度（在渗流问题中对应井的 产量）在上述问题中，若k和q只是空间位置的函数时，问题是线性的若化和q是0及其导 数的函数时，问题则是非线性的由于微分方程组（2-1）在域中每一点都必须为零，因此就有（2-5）j VtA（u）d三 j （vA （u） + v A （u） + ・・・）d三0Q Q 1 1 2 2其中（2-6）其中V是函数向量，它是一组和微分方程个数相等的任意函数式（2-5）是与微分方程组（2-1 ）完全等效的积分形式我们可以说，若积分方程对于 任意的V都能成立，则微分方程（2-1）必然在域内任一点都得到满足。

同理，假如边界条 件（2-2 ）亦同时在边界上每一点都得到满足，对于一组任意函数，下式应当成立j VB（u）dr 三 j （vB （u） + vB （u） + ・・・）dr 三 0r r 1 1 2 2因此积分形式j VtA（u）dQ + j VtB（u）dr = 0对于所有的V和V都成立是等效于满足微分方程（2-1）和边界条件（2-2）我们把（2-7） 式称为微分方程的等效积分形式2.1.2等效积分的“弱”形式在一般情况下，对（2-7）式进行分部积分得到另一种形式：j Ct （v）D（u）dQ + j Et （v ）F（u）dr = 0 （2-8）Qr其中C，D，E，F是微分算子，它们中所包含的导数的阶数较（2-7）式的低，这样 对函数u只需要求较低阶的连续性就可以了在（2-8）式中降低连续性要求是以提高V和V 的连续性要求为代价的，由于原来对V和V （在（2-7）式中）并无连续性要求，但是适当提高对其连续性的要求并不困难，因为它们是可以选择的已知函数这种降低对函数u连续 性要求的作法在近似计算中，尤其是在有限单元法中是十分重要的2-8）式称为微分方程 （2-1）和边界条件（2-2）式的等效积分“弱”形式。

值得指出的是，从形式上看“弱”形 式对函数u的连续性要求降低了，但对实际的物理问题却常常较原始的微分方程更逼近真正 解，因为原始微分方程往往对解提出了过分“平滑”的要求2.1.3加权余量法在求解域中，若场函数u是精确解，则在域中任一点都满足微分方程（2-1）式， 同时在边界r上任一点都满足边界条件（2-2）式，此时等效积分形式（2-7）式或（2-8） 式必然严格地得到满足但是对于复杂的实际问题，这样的精确解往往是很难找到的，因此 人们需要设法找到具有一定精度的近似解对于微分方程（2-1）式和边界条件（2-2）式所表达的物理问题，未知场函数u可以采 用近似函数来表示近似函数是一族带有待定参数的已知函数，一般形式是u 牝 u = ^Ln a = Na （ 2-9）i=1其中，ai是待定参数；N,是试探函数（或称基函数、形函数），为已知函数，它取自 完全的函数序列，是线性独立的所谓完全的函数系列是指任一函数都可以用此序列表示 近似解通常选择使之满足强制边界条件和连续性的要求例如当未知函数u是压力时，可取 近似解u = Nu + Nu +•.. + N u =£nu1 1 2 2 n n i ii=1其中a,是待定参数，共有n个。

显然，在通常n取有限项数的情况下近似解是不能精确满足微分方程（2-1）式和边界 条件（2-2）的，它们将产生残差R及RA（ Na） = R; BA（ Na） = R残差R及R亦称为余量在（2-7）式中我们用个规定的函数来代替任意函数v及v，即u = Nu + N u +••. + N u = Nu1 1 2 2 n n i ii=1可以得到近似的等效积分形式j W A（Na）dQ + f WB（Na）dr = 0 （j = 1~ n） （2-10）j r j亦可以写成余量的形式j W RdQ + f WRdr = 0 （j = 1~ n） （2-11）j r j（2-10）式或（2-11）式的意义是通过选择待定系数at，强迫余量在某种平均意义下等 于零Wj和Wj称为权函数余量的加权积分为零就得到了一组求解方程，用以求解近似解的待定系数「，从而得到原问题的近似解答求解方程(2-10)的展开形式是\ WA(Na)dQ + J WB(Na)d「= 0 q 1 r 1J WA(Na)dQ + J WB(Na)d「= 0q 2 r 2J WA(Na)dQ + J WB(Na)d「= 0Q n r n其中若微分方程组A的个数为I，边界条件B的个数为气，则权函数W.(，• = L..顷 是m阶的函数列阵，W.(j = 1,...,n)是m2阶的函数列阵。

当近似函数所取试探函数的项数n越多，近似解的精度将越高当项数n趋于无穷时， 近似解将收敛于精确解对应于等效积分“弱”形式(2-8)式，同样可以得到它的近似形式J Ct(W.)D(Na)dQ + J Et(W.)F(Na)dr = 0 (j = 1,...,n) (2-12)采用使余量的加权积分为零来求得微分方程近似解的方法称为加权余量法加权余量法 是求微分方程近似解的一处种有效方法常用的权函数的选择有以下几种：(1) 配点法，这种方法相当于简单地强迫余量在域内n个点上等于零；(2) 子域法，该方法的实质是强迫余量在n个子域Q.的积分为零；(3) 最小二乘法，此方法实质是使得近似解和权函数组成的泛函取最小值；(4) 力矩法，该方法是强迫余量的各次矩等于零，通常又称此法为积分法；(5) 伽辽金法(Galerkin)加权余量法可以用于广泛的方程类型，选择不同的权函数，可以产生不同的加权余量法； 通过采用等效积分的“弱”形式，可以降低对近似函数连续性要求当近似函数满足连续性和 完备性要求、试探函数的项数不断增加时，近似解可趋近于精确解由于Galerkin具有广泛 的适用性，因此，下面简单介绍其基本原理：取W. = J，在边界上W. = -W. = -N.，即简单地利用近似解的试探函数序列作为权 函数。

近似积分形式可以写成J N辰N a )dQ + J N tB(£Na )dr = 0 (j = 1,…，n) (2-13)Q j . i i r j . i i由(2-9)式，可以定义近似解U的变分"为8 u = Nd a + N 8 a + • . • + N 8 a其中8a是完全任意的2-13)式可更简洁地表示为J 8 utA(u)dQ + J 8 utB(u)dr = 0q r对于近似积分的“弱”形式(2-12)式则有J Ct (8 U)D(U) dQ + J Et (8 U)F (U) dr = 0q r我们将会看到，在很多情况下，采用伽辽金法得到的求解方程的系数矩阵是对称的，这 是在用加权余量法建立有限元格式时几乎毫无例外地采用伽辽金法的主要原因，而且当存在 相应的泛函时，伽辽金法与变分法往往导致同样的结果2.2变分原理讨论一个连续介质问题的变分原理首先要建立一个标量泛函n，它由积分形式确定8u 、一n =J F u,丁，…dQ + J E u,丁，…dr （2-14）q k 办 J r \ dx J其中，u是未知函数，F和E是特定的算子，是求解域，r是n称为 未知函数的泛函，它随函数u的变化而变化。

连续介质问题的解u使泛函n对于微小的变 化5 u取驻值，即泛函的“变化”等于零5 n = 0 （2-15）这种求得连续介质问题解的方法称为变分原理或变分法如前所述，连续介质问题中经常存在着和微分方程及边界条件不同的，但却是等价的表 达形式，变分原理是另一种表达连续介质问题的积分表达形式在用微分公式表达时，问题 的求解过程是对具有已知边界条件的微分方程或微分方程组进行积分在经典的变分原理表 达中，问题的求解过程是寻求使得具有一定已知边界条件的泛函（或泛函系）取驻值的未知 函数（或函数系）这两种表达形式是等价的，一方面满足微分方程及边界条件的函数将使 泛函取极值或驻值，另一方面从变分的角度来看，使泛函取极值或驻值的函数正是满足问题 的控制微分方程和边界条件的解应注意到，经常有些物理问题可以直接用变分原理的形式来叙述，如表述力学体系平衡 问题的最小位能原理和最小余能原理等，但是并非所有以微分方程表达的连续介质问题都存 在这种变分原理研究表明，原问题等效积分的Galerkin提法等效于它的变分原理，即原问题的微分方程 和边界条件等效于泛函的变分等于零，亦即泛函取驻值反之，如果泛函取驻值则等效于满 足问题的微分方程和边界条件，而泛函可以通过原问题的等效积分的Galerkin提法而得到。

Galerkin法的适用性比变分原理要强，原因是对于有的微分方程很难找到对应的泛函或根本找不到泛函，这时变分原理不适用，但Galerkin法仍然适用如前所述，无论是加权余量法还是变分原理，虽然可以得到微分程的近似解，但是由于 它是在全求解域中定义近似函数，因此实际应用中会遇到两方面的困难（1） 在求解域比较复杂的情况下，选取满足边界条件的试探函数，往往会产生难以克 制的困难，甚至有时做不到2） 为了提高近似解的精度，需要增加待定参数，即增加试探函数的项数，这就增加 了求解的繁杂性而且由于试探函数定义于全域，因此不可能根据问题的要求，在求解域的 不同部位对试探函数提出不同精度的要求，往往由于局部精度的要求使整个问题的求解增加 许多困难变分有限元法和加权余量有限元法就是分别以变分原理和加权余量法为理论基础，通过 对求解区域进行单元剖分，把整个的求解区域剖分成有限的小区域子域，然后在子域内定义 近似函数（近似解），因此称为变分有限元法和加权余量有限元法变分有限元法和加权余 量有限元法虽然在本质上与变分法和加权余量法是类似的，但由于近似函数在子域（单元） 上定义，因此可以克服上述两方面的困难，并由于和现代计算机技术的结合，使得有限元法 成为对物理、力学以及其它科学技术领域问题进行分析、求解的有效工具。

2.3有限元方法的一般步骤在有限元法中，把所研究的连续介质表示为一些小部分（称为有限元）的集合这些单 元可认为是一些称为结点的指定结合点处彼此连接的这些结点通常是置于单元的边界上， 并认为相邻单元就是在这些边界上与它相连接的由于不知道连续介质内部的场变量（在固 体力学中如位移、应力，在渗流问题中如压力、饱和度）真实的变化，因此，我们假设有限 元内场变量的变化可以用一种简单的函数来近似这些近似函数（也称为插值模式）可由场 变量在结点处的值确定当对整个连续介质写出场方程组（如平衡方程组）时，新的未知量 就是场变量的结点值求解场方程组（通常以矩阵方程形式表示），即得到场变量的结点值 一旦知道了这些结点值，则可由近似函数确定整个单元集合体的场变量有限元法求解一般的连续介质问题时，总是依次逐步进行的以与时间无关的物理问题 为例，说明有限元法的基本步骤见图2-21） 结构或求解域的离散化有限元法的第一步，是把求解域分割成许多小部分或称 为单元，因而对于一个具体的有限元分析问题，首先要用适当的有限元把结构进行剖分，并 确定单元的数量、类型、大小和布置2） 选择适当的插值模式由于在任意给定的约束作用下，问题的准确解为未知，因 此，我们假设用单元内的一些适当解来近似未知解。

从计算的观点看，假设的解必须简单， 而且应当满足一定的收敛性要求通常，把解的插值模式取为多项式形式3） 单元分析即进行单元刚度矩阵和载荷向量的推导根据假设的插值模式，利用 平衡条件或适当的变分原理，就可以推导出单元的刚度矩阵K^和载荷向量七，形成单元 平衡方程KePe = Fe（4） 总体合成集合各单元方程以得到总的平衡方程（组）由于结构是由若干个有限 元组成的，因此，应当把各个单元刚度矩阵和载荷向量按适当方式进行集合，从而建立如下 形式的总的平衡方程：KP = F其中，K称为集合刚度矩阵，或称总体刚度矩阵；P是整体结构的结点参数向量，F 是它的结点载荷向量在不同领域的问题中，P所代表的物理量含意不同，如在固体力学 问题中P代表结点处的位移，在渗流力学问题中P代表结点处的压力，在热学问题中P代 表结点处的温度图2-2有限元分析的一般过程（5） 引入约束条件在总体平衡方程的基础上，按问题的边界条件修改总的平衡方程 考虑了边界条件后，可以把平衡方程表示为KP F（6） 方程求解对上述方程进行求解，对于线性问题可以很容易地解出向量P，而对 非线性问题则要经过一系列的步骤才能求得解答，每一步都要对刚度矩阵K和载荷向量F 进行修正。

7） 计算其它参数在直接求得结点变量后，可由此计算其它参数，对于渗流问题可 求压力梯度和流量对于热学问题，可求温度梯度和传热量如对于固体力学问题，可求应变和 应力等。