用数值方法构造二阶收敛的平面弹性问题非闭锁三角形元.pdfVIP

第25卷 第5期 中州大学学报 2o08年 10月 V01．25 No．5 JOURNAL 0F ZH0NGZHOU UNIVERSITY Oct．2008 用数值方法构造二阶收敛的 平面弹性问题非闭锁三角形元 任 国彪 (郑州大学数学系，郑州450O52) 摘 要：对于纯位移平面弹性问题，在不完全3次多项式空间中，用数值方法基于div ∈P。构造了一个二阶收 1 敛的非协调三角形单元，该单元是非闭锁的，数值算例验证了该单元的收敛性结果。 关键词：平面弹性；非闭锁；有限元；数值方法 中图分类号：0242．21 文献标识码：A 文章编号：10o8—37l5(2008)05一Ol13—05 ． 1引言 学界一致关心的问题。I．Babuska和M．suri针对一般 c- 本文考虑下述各向同性、均匀介质的纯位移齐次平面弹 king问题利用协调元进行了系统研究，指出高阶有限元格式 性问题 ：设Qc 是凸多边形区域 ，位移函数 =(u，u) 几乎都是 Locking—free的，而低阶有限元格式大都是 c— ∈( (n))满足 king的。1992年，s．c．Brenner和L_Y．Sung利用C—R元逼 口一(：A)grad(di)= inn， (1．1) 近纯位移边界条件的平面线性弹性问题，由凸多边形上平面 I ： onan 弹性问题的一致先验估计和一致正则性，获得了一致最优误 其中 ∈( (力))，(，A)∈[。，2]×(0，*)为 Lam6 差估计。。王烈衡和齐禾通过构造一个分片 投影算子的方 常数，且(O 。 )，grad，div分别表示梯度算子和散度算 法 ，提出了一个八节点和一个十节点的非协调矩形元 ，它们 子 ： 都是 king—free的。明平兵研究了非协调旋转 Q元，设 计了一种非协调有限元格式对纯应力边界条件的平面线性 = + ． ， 弹性问题是 cking—free的。以往的这些研究工作主要集 (1．1)的弱形式为 ： 中在能量模误差达到1阶的有限元格式上 ，本用数值方法基 于div；∈R．，在不完全三次多项式空间构造一个 15自由度 l 2 的IJocking—free非协调三角形元，该格式的 一模和能量 o(矗，；)= ；) v =( ， ) ∈( (力))． 其中