求解重调和方程的混合MLS方法.pdfVIP

维普资讯 第8卷 第 13期 2008年7月 科 学 技 术 与 工 程 Vo1．8 No．13 July 2008 1671—1819(2008)13—3436—05 ScienceTechnologyandEngineering ⑥ 2008 Sci．Tech．Engng． 力学 求解重调和方程的混合 MLS方法 雷俊丽 (湖北经济学院统计与数学系 ，武汉430205) 摘 要 以重调和方程的混合变分形式为基础，采用移动最小二乘方法建立插值形函数空间，给出了重调和方程的混合 MLS 数值解法。这种方法降低 了对解空间的光滑度要求，同时高精度地获得了解的高阶导数值，并且避免了求解过程中划分网格 所带来的计算复杂性等问题。最后给出了这一方法的误差估计。 关键词 重调和方程 混合变分形式 移动最小二乘 混合MLS 中图法分类号 035； O241； 文献标志码 A 无网格方法是一种完全基于点近似的新兴数 有限元法和MF法的优点，处理起来更方便，精度也 值方法，由于该方法不需要节点间的拓扑连接关 更高。 系，完全避免了传统方法在求解时不断重分网格所 本文首先对一些数学符号和意义做了简要说 带来的精度降低及计算复杂性增加等问题，近来引 明，然后介绍了移动最小二乘方法的基本原理及其 起了学术界的广泛关注。目前已提出了一系列的 插值误差估计。之后引入 了重调和方程的混合 无网格方法，这些 方法基本上可 以分为三种类 MLS方法并对该方法的误差进行了分析。 型…：基于强形式的 MF方法 (如配点法等)，基于 弱形式的MF方法 (如 EFG，MLPG等方法)和粒子 1 准备工作 法 (如MPS，MD，LBE等)。其中，目前提出的基于 弱形式的MF方法均是建立于古典的Galerkin变分 记号Ol=(Ol。，Ot：，…，Ol)称为多重指标，其中 原理或 Ritz变分原理之上的。本文针对重调和方 Ol(i=1，2，…，n)为非负整数 ； 程提出了一种基于混合变分形式的混合无网格方 定义 lOll=Ol1+Ol2+… +Ol； 法，其中混合MIS形函数空间采用移动最／b--乘方 设 =(。， ，…，)为 中的任一元素，记 法建立。这种方法较之传统的MF方法有许多显著 ·X 2 … X‰ n ； 的特点，例如这种方法对解空间的光滑度要求降 假设 “为任意光滑的函数，记D “--xIvax2…呓 = 低 ，从而解空间更易构造；不仅能直接求出方程 的 宣!些 解本身 (力学中的位移)，还能同时求得解的更高阶 a a … a 导数值(应力或应变)，较之传统的有限元法和MF 定义通 常 的 Lebesgue空 间和 Sobolev空 间 法等数值方法是在求出了位移之后，再通过数值算 如下： 法求应力或应变的方法，不仅简化了计算过程，而 定义1 Lebesgue空间 ()定义为： ()= 且减小了计算误差；由于这种方法同时兼备了混合 {／lll／ll(n)∞}，1≤P≤∞ 其中，ll川 ( =( l )1) ，1≤p∞； 2008年3月6日收到 作者简介：雷俊丽(1980一)，女，E—mail：leijunli@126．coln。 Ilfll I 川 l。