有限元分析及应用第六章.pdfVIP

有限元分析与应用 霍战鹏 第六章 非协调单元 第二章至 五章的讨论以最小势能原理为基础，要求在单元内假设的位移场（试探函数） 满足协调条件（在不同的单元内可以假设不同的的位移场）。满足协调条件的单元，它们的 收敛性等问题已在 四章中做了研究。等参数单元就是目前处理二阶问题时应用最广的一种 协调单元。 此外，还有一些单元，它们不满足协调条件，但仍可以收敛到真实解，这类单元称为非 协调单元，可以看成是对等参数单元的一种改进，目的在于 在计算量增加不多的情况下， 使单元的实际精度有所改善。对于四阶问题（例如板、壳），协调条件要求单元之间位移和 位移的一阶导数（转角）连续。在 七章中将会看到，实现上述协调条件不是件容易的事， 而且为此要增加相当大的计算量，因而人们在自编程序中常常对非协调单元感兴趣。 本章只讨论二阶问题，主要包括 非协调元的构造和分析方法，非协调元的理论基础（显 然不能再利用最小势能原理），收敛判别方法。这些结论对四阶问题同样适用。 从关于非协调元的讨论中，读者可以看到，有限元方法有了坚实的数学基础以后，在构 造方法时思路可以开阔很多。 6-1Wilson 非协调元 Wilson 非协调元可以看成是由等参数单元演变来的单元，现以二维情况为例。 １、母体单元 形函数 η 母体单元 边长为２的正方形 ４ ３ 自然坐标 η （1, 1） （-1,1） 取四个角 为节 ，在单元内的序号为１~4 。 形函数 ξ 1 N (x ,h) (1 x x )(1 h h) （i 1~ 4） i i i 4 ２、实际单元 e １ ２ 可看成母体单元 经变换F 得到 （-1,-1 ） （1,-1 ） ˆ F e e 图６－１ 利用上面定义的形函数，坐标的变换可写成 4 x N x 3 i i (x , y ) 3 3 i1 4 y N y y,v i i 2 i1