有限元分析第六章.docVIP

非协调单元 第二章至第五章的讨论以最小势能原理为基础，要求在单元内假设的位移场（试探函数）满足协调条件（在不同的单元内可以假设不同的的位移场）。满足协调条件的单元，它们的收敛性等问题已在第四章中做了研究。等参数单元就是目前处理二阶问题时应用最广的一种协调单元。 此外，还有一些单元，它们不满足协调条件，但仍可以收敛到真实解，这类单元称为非协调单元，可以看成是对等参数单元的一种改进，目的在于：在计算量增加不多的情况下，使单元的实际精度有所改善。对于四阶问题（例如板、壳），协调条件要求单元之间位移和位移的一阶导数（转角）连续。在第七章中将会看到，实现上述协调条件不是件容易的事，而且为此要增加相当大的计算量，因而人们在自编程序中常常对非协调单元感兴趣。 本章只讨论二阶问题，主要包括：非协调元的构造和分析方法，非协调元的理论基础（显然不能再利用最小势能原理），收敛判别方法。这些结论对四阶问题同样适用。 从关于非协调元的讨论中，读者可以看到，有限元方法有了坚实的数学基础以后，在构造方法时思路可以开阔很多。 §6-1Wilson 非协调元 Wilson 非协调元可以看成是由等参数单元演变来的单元，现以二维情况为例。 母体单元　形函数 母体单元ê：边长为２的正方形 自然坐标：ξ、η 取四个角点为节点，在单元内的序号为１~4。 形函数 实际单元　e 可看成母体单元ê经变换F得到 利用上面定义的形函数，坐标的变换可写成 其中（xi, yi）为实际单元中节点的坐标。 至此，还看不出Wilson非协调单元与上一章介绍的等参数单元之间的差别。 单元内假设位移场 同四节点等参元相比，单元内假定的位移场多了四项： 它们有如下特性： 不影响节点处的位移值，故称αl为非节点自由度或单元的“内自由度”。在计算单元变形能和单元体积力做功时计入这些位移；但在计算边界外力做功（为了将边界力化为等效节点力）时不计这些位移。即在计算边界外力做功时只计Niui、Nivi各项。 补充这些项后，单元内的位移场是ξ、η的完全二次多项式。当实际单元e为矩形时，单元内位移场将是x、y的完全二次多项式。 协调性分析 　　沿单元的一边，例如节点１、２所在的边，η＝－１。u,v是ξ的二次函数，完全被u1, v1,α1, 和u2, v2,α3所决定。但由于不同单元的α1~α4彼此独立，故不能保证单元之间位移的协调性。对于不满足协调条件的单元，显然不能再用最小势能原理。尽管如此，人们仍然照搬协调单元的一些具体做法（这样，对协调元编制的程序基本上可以照搬）。至于能否保证收敛到真实解等问题，留待下面几节加以讨论。 现假如我们研究某个具体的平面应力问题。节点总数为n，单元总数为m。则总的自由度可区分为： 节点自由度 单元内自由度 系统的总势能定义为 其中　　　号单元的变形能；　　　号单元体积力做功；　　　为各边界外力在位移 　　　　　　　　　上做的功之和。 我们将把由 求得的ui, vi, αl(j)以及由此求得的应力做为非协调单元的的有限元解。 在（6-1-3）中共有2n+4m个未知量。比四节点等参元多了４m个未知量。但是α1(j)、α2(j)、α3(j)、α4(j)仅属于第j号单元，故有 这样，就可以在单元分析时先消去αl(j)（这一步骤称为静凝聚），只剩下ui, vi进入总体平衡方程。 ４、单元分析　　静凝聚 在单元分析中，节点仍取它在单元内的序号，并约定： 单元的外自由度为： 单元的内自由度为： （6-1-1）所定义的单元位移场可写成 几何矩阵 应变 其中[B]为几何矩阵。 单元刚度矩阵和体积力载荷向量 　单元刚度矩阵， 阶的方阵 单元变形能 由于[k]为对称阵，必有 以及 体积力　　　　 做功 其中体积力载荷向量 静凝聚 略去（6-1-4）中的单元编号j，以（6-1-9），（6-1-10）代入，（变形能对内部自由度取偏导）并注意到 则有 解得 将（6-1-12）代入（6-1-9）和（6-1-10）有 （6-1-13）右端第三项与{uΙ}无关，不影响πPh取驻值。第一项为{uΙ}的二次型，[k]为凝聚掉内自由度后的单元刚度矩阵。{r}为凝聚内自由度后的载荷向量。它们的计算式为 为了避免对一个 　　的方阵[kⅡ]求逆，静凝聚可以分四次进行，每次只凝聚一个内自由度。 （４）　边界外载荷　　　 　的等效节点力 积分沿单元受到载荷的边界SP进行。其中　　 为作用在单位面积上的边界载荷。 ５、组装及求解总体方程 由[k]组装总体刚度矩阵[k]。由{r}、{rs}