弱不连续问题高阶有限元离散系统的斍斄斖斍法.PDFVIP

卷 期 棳 棾棿 棻 计算 力 学 学 报 斨旓旍棶棾棿 斘旓棶棻 年 月 旇旈旑斿旙斿斒旓旛旘旑斸旍旓旀 旓旐旛旚斸旚旈旓旑斸旍斖斿斻旇斸旑旈斻旙 棽棸棻椃 棽 斊斿斺旘旛斸旘棽棸棻椃 旔 旟 椇 肖映雄 斈斚斏棻棸棶椃椀棻棻 旙旍旞棽棸棻椃棸棻棸棸棿 棽棸棻椀棸椆棸棸棸棽 旉 弱不连续问题高阶有限元离散系统的斍斄斖斍法 灣 棳 棳 肖映雄 王 彪 李真有 棬 棳 棭 湘潭大学 土木工程与力学学院 湘潭 棿棻棻棻棸椀 椇 棬 棭 暎 摘 要 弱不连续问题 如含夹杂问题 是固体力学计算中的一类重要问题 高阶有限元方法由于其具有更好的 棳 暎 棳 逼近效果 是确保数值解在界面保持较高精度的计算方法之一 但与线性元相比 高阶单元需要更多的计算机存 棳 暎 棳 储单元 具有更高的计算复杂性 本文利用两水平算法的思想 将高阶有限元离散系统化归于线性元离散系统的 棳 棬 棭 棳 求解 为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格 斍斄斖斍 法 并 暎 棳 棳 应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解 数值试验结果表明 相比于常用 斍斄斖斍法 新 暍 棳 棳 方法的迭代次数基本不依赖于问题规模 单元阶次以及杨氏模量的间断性 斆斝斦计算时间得到明显改善 具有更 棳 暎 好的计算效率和鲁棒性 可大大提高弱不连续问题有限元分析的整体效率 椇 椈 椈 椈 椈 关键词 弱不连续问题 高阶单元 条件数 两水平方法 代数多重网格法 中图分类号椇 椈 文献标志码椇 文章编号椇 棬 棭 斚棾棿棶 斣斅棻棻椀 斄 棻棸棸椃灢棿椃棸椄棽棸棻椃棸棻灢棸棸椀灢棸椄 棳 目前有许多这方面的研究工作 如基于单元刚度矩 棻 引 言 椲棳椵