基于比例移动最小二乘近似的误差分析-应用数学和力学.PDF

　 ， ３８ １１ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１７年１１月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３８，Ｎｏ．１１，Ｎｏｖ．１５，２０１７ 文章编号：１０００⁃０８８７（２０１７）１１⁃１２８９⁃１１ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 基于比例移动最小二乘近似的误差分析∗ 王青青，　 李小林 （重庆师范大学 数学科学学院，重庆４０１３３１） 摘要：　 相较于移动最小二乘近似方法，比例移动最小二乘近似法有效地克服了前者带来的矩阵 病态这一问题，展示出了更好的数值稳定性和更高的计算精度．给出了比例移动最小二乘近似对函 数及其任意阶导数的误差估计，并给出了数值算例来验证之前的理论分析结果，通过与移动最小 二乘近似的比较，表明比例移动最小二乘近似能得到更快的收敛性和更稳定的计算性． 关　 键　 词：　 无网格方法；　 比例移动最小二乘近似；　 稳定性；　 误差估计 中图分类号：　 Ｏ２４２．２　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．２１６５６／ １０００⁃０８８７．３７０２６０ 引　 　 言 近几十年来，有限元方法（ＦＥＭ）和边界元方法（ＢＥＭ）已成为计算力学中解决工程问题的 主要数值计算方法．但在近似时，基于网格，在处理一些碰撞问题和大变形问题时，需要进行网 格重构，比较费时．而无网格方法，采用基于点的近似，可以彻底或部分地消除网格，不需要网 格的初始划分和重构，不仅可以保证计算的精度，而且可以减小计算的难度． ［１］ 形函数是无网格方法的基石，移动最小二乘近似（ＭＬＳ） 是构造形函数最重要的方法．在 ＭＬＳ 中，都是基于离散的节点近似而不是基于网格的．因此，许多基于移动最小二乘近似的无 ［２］ ［３］ 网格方法也蓬勃发展起来，比如无单元Ｇａｌｅｒｋｉｎ（ＥＦＧ）方法 、无网格Ｐｅｔｒｏｖ⁃Ｇａｌｅｒｋｉｎ方法 、 ［４］ ［５⁃６］ 边界点方法（ＢＮＭ） 、Ｇａｌｅｒｋｉｎ ＢＮＭ 等． ［７⁃８］ 基于传统的最小二乘法，ＭＬＳ 的系数矩阵的条件数可能会变得很大 ．因此，对系数矩阵 ［９］ 取逆运算可能会导致在计算稳定性和计算精度方面的下降．近来，基于带权正交基函数 的改 进的移动最小二乘近似也发展起来，在改进的移动最小二乘近似方法里，尽管最后代数方程的 ［１０］ 求解不用通过矩阵求逆运算，但条件数增大这一问题仍然存在 ． 本文首先建立了比例移动最小二乘近似（ＳＭＬＳ）．通过选择比例基函数［８，１１］ ，ＳＭＬＳ 的系数 矩阵将不会变得病态．因此，通过ＳＭＬＳ可以避免由ＭＬＳ带来的坏条件这一缺点，同时，也期待 得到更高的计算稳定性和精确度． 形函数在无网格方法中占有非常重要的地位，ＭＬＳ 的误差是分析基于ＭＬＳ 的无网格方法 误差的基础．文献［４⁃８，１１⁃１５］在不同的假设下，通过不同的范数给出了ＭＬＳ 的误差估计．同 ∗ 收稿日期：　 ２０１６⁃０８⁃２６；修订日期：　 ２０１７⁃０９⁃２８ 基金项目：　 国家自然科学基金（面上项目）；重庆市基础科学与前沿技术研究重点项目 （ｃｓｔｃ２０１５ｊｃｙｊＢＸ００８３）；重庆市教委科学技术研究项目（ＫＪ１６００３３０） 作者简介：　 王青青（１９９２—