一个平面弹性问题的虚位移原理和Locking-Free非协调有限元方法.pdfVIP

第42卷第5期 数学的实践与认识 V01．42，No．5 IN AND 2012年3月 MATHEMATICSPRACTICETHEORY MaLr．．2012 非协调有限元方法 赵中建·，谷留新t，陈绍春z (1．华北水利水电学院数学与信息科学学院，河南郑州450011) (2．郑州大学数学系，河南郑州450002) 摘要：证明了平面弹性问题的虚位移原理，提出了一个新的非协调有限元解决纯 敛性，并且证明了误差估计结果与参数A无关． 关键词：非协调有限元；Locking-Free；虚位移原理；最优误差估计 1引言 对于各向同性、均匀介质的平面弹性问题，当Lam6常数入—哼CX3时，即对于几乎不可 压介质，通常的协调有限元格式的解，往往不再收敛到原问题的解，或者达不到最优的收敛 与A有关，当A—一。。时，该系数将趋向无穷大．详细解释见Arnold[1|，Babu§kaSuril3_剖， 量较大． 本文证明了平面弹性体的虚位移原理Is]，所用单元是在Rannacher Robust和Optimal，并且能量模的收敛估计与参数入无关． 2位移边界条件下的平面弹性问题和虚位移原理 设QC (咒，B)T为体力． 考虑纯位移边界条件下的平面弹性问题 ：竺舍，矗一(p+A)grad(diV回=户，i扎。0’aQc·， 收稿日期：2010-11．07 资助项目：国家自然科学基金 5期 181 赵中建，等：一个平面弹性问题的虚位移原理和Locking-Free非协调有限元方法 其中(A，肛)∈(0，+。。)X[肛1，p2]，0p1p2． 设o(，)：[础(Q)]2 X[础(Q)]2一R是对称的双线型定义如下 (2) 其中矗∈旷=(硪(Q))2． (1)式的等价的变分形式为：求面∈矿，对于给定的F∈[L2(Q)]2 f o(矗，回=(F，回，V秽∈V， 1(户，叼：p谢zd可 (3) 定义1虚位移指任何无限小的位移，且它在物体结构内部必须是连续的，在边界上必须 满足运动学边界条件． 以平面应力问题为例，考虑单位厚度的平面弹性体，为简化起见，假定不存在体积力和初 应力，边界可分为两部分，在r。上位移等于零，在L上作用着边界力d，对应分量G。，G， 在体力作用下物体内产生的应力如，屯，亿Ⅳ，由弹性力学可知，这些应力分量满足下列平衡方 程 等+鲁-o，等+鲁=。 ㈥ 体力F与位移何的关系见(1)式，在r，这部分的边界条件为 如如+亿Ⅳ勺=只：， 如％+％掣如=玛 (5) 其中如，％为边界法线的方向余弦． 假定物体产生了虚位移矗=(“-，U2)T，体力在虚位移上所做的虚功为 aV=见Ul+日U2={田1{F} 对应于虚位移，相应的虚应变为 a乱】 a钍2 a札】 a札2 Ez2面，￡Ⅳ2瓦，％!，2瓦+瓦 由于上述虚应变，物体内部应力所产生的虚应变能为 ， 盯u=ii