弹性接触问题的对偶混合有限元分析.pdfVIP

维普资讯 1999年 l1月 计 算 数 学 第 21卷第 4期 NOV．．1999 MATHEM ATICA NUMERJCA SINICA V01．21．NO4 弹性接触问题的对偶混合有限元分析 D 弓 一 1L 王烈衡王光辉 0z z／ ． (中科院计算教学与科学工程计算研究所．科学与工程计算国家重点试验室) DU AL M IXED FIN ITE ELEM NT M ETHoD FoR CoNTACT PRoBLEM IN ELASTICITY W angLieheng WangGuanghui (StateKeyLaboratoryofScientific~EngineeringComputing, InstituteofComputationalMathematicsand$cieutifci／Emgineeri~ Computing, ChineseAcademyofSciences，Be{ling) Abstract Inthispaper，basedonthemixedvariationalformulationin[9】1adualmixed variationalformulation forcontactpro blem in elasticity ispresented．The exis- tenac anduniquenessofthesolutionofthedual variational problem arediscussed， andtheerrorboundo(h3／)isobtainedofrRaviart—Thomas(k=1)element approximation． Keywords：Contactproblem in elasticity,Dualmixedvariational form ulation，Ravlart-Thom aselem ent 关键词 弹性接触 问囊．对偶混台变分形式． Raviart-Thomas元 胜刀芋 R—T无 既言彳i腑 1．引 言 用混合有限元方法求解弹性力学问题，其优点在于可同时求解位移和应力．力学问题的 混合变分形式是混舍有限元方法的基础．对于弹性接触问题，文献 [6]给出了一种混合变分 形式，以及相应的混合有限元分析 (也可见 [1o])．其混合变分形式是直接从位移变分方程和 Hook方程导出的，获得了应力cr在工(Q)而位移u在 日 (Q)的一个闭凸子集上求解的 混合变分 问题． 本文在 9『1中提出的混合变分形式的基础上，再引入另一个 Lagrange乘子，获得了三 重组混合变分形式．它能同时求解物体 内点的应力，位移和接触边界上的位移．而应力 在 H(div，Q)，位移矗在 (Q)中．值得注意的是，这里的位移是物体实际运动的位移，它包 括除小弹性范围以外的刚体位移，当部分边界固定，这种位移即为弹性位移． +1998年 11月 10 日收到． 1)国家自然科学基金资助项目 维普资讯 484 计 算 数 学 1999芷 本文的安排如下 ：第 2节给出一些记号，叙述无摩擦问题