分片实验与有限元法(2)

    可以增加位移函数中的广义参数，通过分片试验的条件消去这些多余的广义参数，这样得到的位移插值函数会得到改善或完全满足分片试验的要求。这种方法的实质是改善了位移函数的空间，但它的应用还非常少，其主要原因是计算中涉及求逆运算。目前技术及软件的高速发展，尤其是代数运算软件的出现，这种做法也许会有一些生命力。下面举一个通过这种方法改善单元性能的例子。

   在构造三角形单元时，人们呈为完全的三次式中十个基函数的取舍大费周折，面积坐标的应用解决了对称性的问题，但Zienkiewicz元（BCIZ元）的性能不佳也是人所共知的。今位移函数的基取完全的三次式，含十个基函数，采用面积坐标可写成如下形式：

                                        (20)

其中 为Zienkiewicz元的单元位移函数， (i=1,2,3)为三个面积坐标，C为待定参数。以下通过C的确定来改善单元的性质。因只有一个待定参数，方程(6)不可能完全得到满足，考虑到对称性将(6)中的前两式相加得到方程：

                                           (21)

对边长为0.5的方板做图示两种网格划分，坐标原点在1点，其中图二中5点坐标为(0.2，0.15)，边界结点的位移参数按任意的二次挠度场 给定，计算5点的挠度及转角，表1列出了Zienkiewicz元和改进的Zienkiewicz元结果。

表1 分片试验

可以看出改进Zienkiewicz元的性能有很大的改善，以下做一算例。

算例：方板中心受集中力，根据对称性，取板的四分之一，采用交叉网格的计算结果如表2。

表2  BCIZ元改进前后板中心挠度计算

    由算例可以看出改进Zienkiewicz元的收敛性能有了很大的改善，而且单元采用的位移函数不仅具有几何对称性，各结点的挠度和转角值也表达精确。在三次位移函数的单元中，这种单元的位移函数的插值空间得到了进一步改进。

5

通过前面的讨论可以看出，各有限元法与分片试验是密不可分的，它们自觉或不自觉得满足了分片试验的要求。这些有限元法合理的共同原因也许在于它们能通过分片试验。

满足了应变约束条件的有限元法，一般是以损失连续性方程的严格性为代价的，这一点对计算结果一般影响不大，而且往往会改善计算精度，这些有限元法对分片试验的满足十分自然，但有些时候会涉及秩的问题；

使用了位移约束条件的有限元法，以损失位移函数在单元结点的准确程度为代价，换取了单元总体性能的改进，或者改善了位移试函数的插值空间，这类有限元法对在保持位移函数的几何对称性上有些困难。以上两类有限元法都得出了很多属于自己特色的单元。

   本文得出的是常应变分片试验的要求，同样可以得出应变或位移在什么情况下，能够通过线性应变的分片试验。如果单元的位移参数较多，位移插值函数已含完全三次多项式，单元片在线性应变情况下也应计算准确，这样才更值得我们增加参数。

（转载自科教范文网http://fw.nseac.com）

参    考    文    献

[1]    O.C.Zienkiewicz and R.L.Taylor, The Finite Element Method, (Fourth Edition), Mcgraw-Hill Book Company, 1988.

[2]    唐立民，有限元分析的若干基本问题，大连院学报，1979，18(2), 1-15

[3]    唐立民，陈万吉，刘迎曦，有限元分析中的拟协调元，大连工学院学报，1980,.19(2), 19-35

[4]    T.H.H.Pian, Derivation of Element Stiffness Matrices by Assumed Stress Distributions, A.I.A.A.J., 1964, 2(7), 1333-1336

[5]    T.H.H.Pian, and Dapeng Chen, Alternative Ways for Formulation of Hybrid Stress Elements, Int. J. Num. Meth. Eng., 1982, 18, 1679-1684

[6]    龙驭球，辛克贵，广义协调元，土木工程学报，1987，1，1-14

[7]    陈绍春，石钟慈，构造单元刚度矩阵的双参数法，计算，1991，3，286-296

[8]    陈万吉，单变量有限元的新思考：精化直接刚度法，计算结构力学及其应用，1993，10(4)：263-268

[9]    Tang Limin, Chen Wanji and Liu Yingxi, String Net Function Approximation and Quasi-Conforming Technique, Hybrid and Mixed Finite Element Methods, S.N.Atluri, R.H.Gallagher and O.C.Zienkiewicz, John Wiley & Sons, 1983.

[10]石钟慈，陈绍春，九参数广义协调元的收敛性，计算数学，1991，2，193-203

Patch Test and Finite Element Method

Abstract: This paper realized that the Patch Test is very important to Finite Element Methods. Derive the requirement to strain and displacement of Patch Test. Give the relationship between weak form of continuity equation and the Patch Test, through which the hybrid element method and quasi-conform element method are analyzed. Refined direct stiffness method and generalized conforming elements are also analyzed about why they can pass patch test. At last, as an example of using the requirement of patch test for displacement, improved the BCIZ Element.

（科教范文网 lw.nseaC.Com编辑发布）

Key Words: patch test, weak form, string-net function, finite element method