部分重叠覆盖流形法的覆盖加密方法-长江科学院院报.PDFVIP

　 第３０卷第７期 长　江　科　学　院　院　报 Ｖｏｌ．３０　Ｎｏ．７　　 　　２０１３年 ７月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｊｕｌ．２０ １３　 　 　　ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００１－５４８５．２０１３．０７．０１９ ２０１３，３０（０７）：９５－１００ 部分重叠覆盖流形法的覆盖加密方法 苏海东，祁勇峰 （长江科学院材料与结构研究所，武汉　４３００１０） 摘要：部分重叠覆盖的流形法是一种新型的数值流形方法，它以独立覆盖的分析为主，便于在各区域采用合适的覆 盖函数以适应物理场的局部特征。以往的研究表明，在实际物理场分布较复杂的区域采用较小的覆盖可以更好地 逼近真实解，这就涉及到与周边较大覆盖的连接问题。针对部分重叠的矩形覆盖提出覆盖加密方法，可以很容易 地将大覆盖加密成为较小的覆盖，并使其与周边的大覆盖保持协调过渡，相对于以往采用有限元网格的流形法而 言在大小网格过渡上要方便得多。应用重力坝应力分析和圆孔应力集中２个算例验证了该方法的有效性。 关　键　词：数值流形方法；部分重叠覆盖；覆盖加密；矩形覆盖 中图分类号：ＴＢ１１５；ＴＶ３１１　　　文献标志码：Ａ　　　文章编号：１００１－５４８５（２０１３）０７－００９５－０６ 的不匹配经常会造成计算精度的下降，流形法的计 ［５］ 算结果往往达不到常规有限元法的精度水平 。 １　研究背景和思路 作者在文献［４，６］中提出部分重叠覆盖的流形 ［１］ 数值流形方法 （以下简称流形法）是留美学 法，采用以独立覆盖为主的分析方式，便于在各区域 者石根华博士提出的一种新型数值计算方法。该方 用合适的覆盖函数以适应物理场的局部特征，独立 法采用了覆盖的概念和２套相互独立的网格体系。 覆盖之间仅用较小的重叠区域以保持覆盖的连续 覆盖，或称为数学覆盖，用以定义各个区域的局部函 性，从而将现有的基于完全重叠覆盖的流形法扩展 数。两套相互独立的网格体系，即反映数值解精度 到一般意义上的流形研究。作为部分重叠覆盖流形 的数学网格和表示几何边界和材料分区的物理网 法的首次尝试，该文基于矩形覆盖，给出部分重叠的 格，将整个研究区域划分成有限个相互重叠的区域。 覆盖形式、基本公式及其实现方式，通过单个矩形数 这些区域被称为物理覆盖，在各个覆盖上独立定义 学网格各结点之间的强制约束方式，很方便地将完 局部覆盖函数，通过联系函数（或称为权函数）加权 全重叠的覆盖形式及公式转化为部分重叠的覆盖形 平均得到整个求解域上的总体函数。２套网格的交 式及公式。 集，称为流形元，跟有限元一样，是基本计算单元，但 众所周知，孔周边等局部区域具有不同于一般 可具有复杂任意的形状。为保证任意形状的流形元 情况的复杂变形和应力分布。如果采用较大的数学 的积分精度，流形法多采用单纯形精确积分方 覆盖、用多项式函数去逼近这样复杂的分布，理论上 式［１－２］。 是可行的，但往往需要很高阶数的多项式级数，而且 当前的流形法研究一般采用有限单元网格来构 高阶情况下由于方程组的性态不佳，易引入数值误 ［１，３］ 差，因此效果很差。实践表明，在大的覆盖中单纯