三维弹性接触问题边界元法的·种误差直接估计.pdfVIP

第 28卷 第 2期 燕 山大学学报 Vol1．28 NO．2 2004年 3月 JournalofYanshanUniversity M ar． 2004 三维弹性接触 问题边界元法的 ·种误差直接估计 姚振汉 刘永健 摘 要 将作者研究组关于二维边界元法误差直接估计的早期工作推广到三维弹性接触问题 ，提出了当域 内点 趋于边界时计算内点位移的极限的算法，并把接触面两边物体相应点的位移极限值之差的某种度量作为三维弹 性接触 问题边界元法误差的一种直接估计。给 出的三维弹性接触问题的算例证明了这种直接估计和相应的 自适 应计算方案的有效性。 关键词 边界元法，误差估计，自适应方案，三维弹性接触 问题。 A DirectErrorEstimationofBEM fOr3D ElasticContactProblem Abstract Anearlyworkondirecterrorestimmionof2D BEM forelasticiytproblem iSextendedto3D elasticcontact problem inthispaperAn algorithm forcomputingthelimitofinnerdisplacementasthepointapproachesboundaryfrom theinternalpointofthedomainispresented．Basedonthedifferencebetweenthecorrespondinglimitofdisplacement from bothsideofconatctedbodies，adirecterrorestimatorofBEM solutionofr3D elasticconatctproblem ispresented． Numericalexamplesaregiventoshow theeffectivenessofsuchdirectellOrestimationandthecorrespondingadaptive scheme Keywords boun daryelementmethod，errorestimation，adaptivescheme，3D elasticconatctproblem． 分方程，误差完全来源于边界条件不能精确满足。 0 引言 通常，我们虽然在求解时精确地代人了边界条件 的 有限元法的误差估计和相应的 自适应方法 已 给定量，但是精确满足域内微分方程的内点变量当 经从理论和方法上得到了很好的发展，但是边界元 趋于边界点时的极限与原来给定的边界给定量并 法的相应研究还有许多工作要做。边界元法误差估 不相同。只有对于没有离散误差的 “简单问题”， 计的近期成果可见于Paulino等的文章 …。一种误 在略去很小的计算误差的情况下，二者基本相等， 差的自然度量最早由Paulino提出乜，它基于同时 此时边界元法解的精度可以达~1]lO一，甚至更高。 利用标准边界积分方程和超奇异边界积分方程，称 基于上述思想，本文针对三维弹性力学问题， 为两级迭代法。该法不仅用于位势问题和线弹性问 提出了用于计算内点位移趋于边界的极限的算法， 题，而且成功用于对称伽辽金边界元法、无网格边 并对三维弹性接触问题，将接触界面两边这种极限 界点法和边界轮廓法等 口。 之差的某种度量作为误差的直接估计，给出了相应 但是本文的工作基于研究组早年对于二维边