无网格法介绍

无网格法是在建立问题域的系统代数方程时，不需要利用预定义的 无网格法是在建立问题域的系统代数方程时， 网格信息，或者只利用更容易生成的更灵活、 网格信息，或者只利用更容易生成的更灵活、更自由的网格进行域 离散的方法。（刘桂荣，2009） 离散的方法。（刘桂荣，2009） 。（刘桂荣
无网格法概述
无网格法求解过程 FEM对比 对比） （与FEM对比）
导出无网格法公式
基于弱强式的无网格法
• MFree弱-强式法 弱 强式法 强式法(NWS)的核心思想是针对某一问题同时采用强式和 的核心思想是针对某一问题同时采用强式和 局部弱式建立起离散系统方程式，即对不同组别的节点根据其不同 局部弱式建立起离散系统方程式， 条件分别形成不同类型的方程，其中局部弱式被用于位于或接近导 条件分别形成不同类型的方程， 数边界条件的所有节点，强式被用于除此之外的其他节点。 数边界条件的所有节点，强式被用于除此之外的其他节点。 • 代表方法：MWS 代表方法： • MWS特点。MWS法使用最少数量的背景网格用于积分，对各类力学 特点。 法使用最少数量的背景网格用于积分， 特点 法使用最少数量的背景网格用于积分 问题均可得到稳定而精确的解，是目前近乎理想的无网格法。 问题均可得到稳定而精确的解，是目前近乎理想的无网格法。
构造无网格形函数
PIM形函数性质
• 一致性 如果单项式的完备阶数是p，则该形函数具有 C p 一致性 如果单项式的完备阶数是 ， • 再生性 PIM基函数可再生包含在其基函数当中的任意函数。 基函数可再生包含在其基函数当中的任意函数。 基函数可再生包含在其基函数当中的任意函数 • 线形独立性 PIM基函数在支持域上是线性独立的 基函数在支持域上是线性独立的 • δ 函数性
目
1 2 3 4 5 无网格法概述 无网格法分类
录
构造无网格形函数 导出无网格法公式 无网格法研究主要进展
构造无网格形函数
无网格插值技术分类 分类 无网格近似技术 再生核粒子法（RKPM） 再生核粒子法（ ） 移动最小二乘（ 级数表达式 移动最小二乘（MLS） ） 点插值法（ 点插值法（PIM,RPIM） ） 单位分解法（ ） 单位分解法（PU） 广义有限差分法（ 差分表达式 广义有限差分法（GFDM） ） 光滑粒子动力学法（ 积分表达式 光滑粒子动力学法（SPH） ）
构造无网格形函数
移动最小二乘(MLS)形函数 移动最小二乘(MLS)形函数 (MLS) 同样，在计算点 处 同样，在计算点x处u(x)可近似表示为 可近似表示为
在移动最小二乘近似（MLS） 系数a(x)的选取使近 在移动最小二乘近似（MLS）中，系数a(x)的选取使近 在计算点x的邻域内待求函数 的邻域内待求函数u(x)在某种最 似函数 u h ( x ) 在计算点 的邻域内待求函数 在某种最 小二乘意义下的最佳近似。 小二乘意义下的最佳近似。近似函数在节点 x i 处的误差 加权平方和为
取最小值， 即可得最小二乘形函数。 令J取最小值，解得待定系数 取最小值 解得待定系数a(x),即可得最小二乘形函数。 即可得最小二乘形函数
构造无网格形函数
MLS形函数性质 形函数性质
• 一致性 如果单项式的完备阶数是p， 如果单项式的完备阶数是 ，则该形函数具有 C p 一致性 • 再生性 PIM基函数可再生包含在其基函数当中的任意函数 基函数可再生包含在其基函数当中的任意函数 • 单位分解性
构造无网格形函数
多项式基点插值法(PIM)形函数 形函数 多项式基点插值法 使用多项式作为插值函数是应用最早的差值方法之一。 使用多项式作为插值函数是应用最早的差值方法之一。 设定义在问题域中的一连续函数u(x)可由一组场节点表 设定义在问题域中的一连续函数 可由一组场节点表 在计算点x处 可近似表示为： 示，在计算点 处u(x)可近似表示为： 可近似表示为
无网格法分类
根据函数近似方法分类
•基于移动最小二乘近似的无网格法 基于移动最小二乘近似的无网格法 扩散单元法（ ），无单元 扩散单元法（DEM），无单元 ），无单元Galerkin法（EFG） ，无网格局部 法 ） Petrov-Galerkin法（MLPG） 法 ） •基于积分形式近似无网格法 基于积分形式近似无网格法 光滑粒子流体动力学法（ ），再生核粒子法 光滑粒子流体动力学法（SPH），再生核粒子法（RKPM） ），再生核粒子法（ ） •基于点插值无网格法 基于点插值无网格法 多项式基点插值法（ 径向基点插值法（ 多项式基点插值法（PIM）,径向基点插值法（RPIM） ） 径向基点插值法 ） •基于其他近似方法的无网格法 基于其他近似方法的无网格法 基于自然邻节点插值的自然单元法（ ），hp云法 基于自然邻节点插值的自然单元法（NEM）， 云法，单位分解 ）， 云法， 法（PU） ）
目
1 2 3 4 5 无网格法概述 无网格法分类
录
Байду номын сангаас
构造无网格形函数 导出无网格法公式 无网格法研究主要进展
无网格法研究主要进展
• Lucy（1977） 提出光滑粒子法（Smooth Particle Hydrodynamics method, Lucy（1977） 提出光滑粒子法（ SPH） 即SPH） • Lancaster（1979，1981）较为系统地研究了移动最小二乘法 Lancaster（1979，1981） • Belytschko（ 1994 ）提出无单元Galerkin方法（Element-Free Galerkin Belytschko（ 提出无单元Galerkin方法（ElementGalerkin方法 method， method，即EFG ） • Liu（1995）等人提出了重构核粒子法（Reproducing Kernel Particle Liu（1995）等人提出了重构核粒子法（ Method， Method，即RKPM) • Oden（1995）提出了Hp-clouds方法 Oden（1995）提出了Hp-clouds方法 Hp • Onate（1996）提出了有限点法（Finite Point Method,即FPM ） Onate（1996）提出了有限点法（ Method,即 • Golberg和Chen（1994，1996）研究了径向基函数法（radial basis Golberg和Chen（1994，1996）研究了径向基函数法（ function， RBF） function，即RBF） • Atluri（2000）提出了Meshless Local Petrov-Galerkin Method（即MLPG） Atluri（2000）提出了Meshless PetrovMethod（ MLPG）
FEM和MFree法流程图
无网格法概述
形成逼近函数 光滑粒子法 移动最小二乘法 单位分解法 重构核粒子法 径向基函数法
与基于网格的方法不同 与基于网格的方法相同
形成求解方程 加权残数法 变分原理 边界积分方程
+
无网格法概述
为何采用无网格法？ 为何采用无网格法？ • • • • 形成FEM网格时的计算成本高 网格时的计算成本高 形成 应力精度低 自适应分析困难 对某些问题分析的局限性 大变形问题（如冲压变形） 大变形问题（如冲压变形） 裂纹扩展问题 流固耦合问题 爆炸问题
构造无网格形函数
• 单位分解性 • 紧支性 PIM形函数是紧支域中的节点构造的，紧支域外的任意点 形函数是紧支域中的节点构造的， 形函数是紧支域中的节点构造的 处的函数值被认为是0 处的函数值被认为是 • 相容性 在应用PIM形函数时，如使用局部支持域，则全局域上的 形函数时， 在应用 形函数时 如使用局部支持域， 相容性将不能保证，即将节点进出或移动支持域时， 相容性将不能保证，即将节点进出或移动支持域时，其场 函数的近似式是不连续变化的
无网格法(MFree) 简介
目
1 2 3 4 5 无网格法概述 无网格法分类
录
构造无网格形函数 导出无网格法公式 无网格法研究主要进展
无网格法概述
无网格法定义
The meshfree method is used to establish a system of algebraic equations for the whole problem domain without the use of a predefined mesh, or uses easily generable meshes in a much more flexible or ‘‘freer’’ manner.(G R Liu,2009)
• δ 函数性 MLS形函数不具备 δ 函数性，本质边界条件不易施加 形函数不具备 函数性，
目
1 2 3 4 5 无网格法概述 无网格法分类
录
构造无网格形函数 导出无网格法公式 无网格法研究主要进展
导出无网格法公式
基于全局弱式的无网格法
• MFree全局弱式法中，控制偏微分方程连同其边界条件通过各种技术 全局弱式法中， 全局弱式法中 转化为一组积分方程， 转化为一组积分方程，再利用建立在问题域上的全局背景网格进行 数值积分操作而将其弱式转化为一组代数系统方程。 数值积分操作而将其弱式转化为一组代数系统方程。 • 代表方法：无单元 代表方法：无单元Galerkin法(EFG)，无网格径向基点插值法 法 ， (MRPIM) • EFG特点。优点：具有很好的精度和收敛性，对点的分布不敏感； 特点。 特点 优点：具有很好的精度和收敛性，对点的分布不敏感； 缺点：需要背景网格，应用本质边界条件困难，计算效率低。 缺点：需要背景网格，应用本质边界条件困难，计算效率低。 • MRPIM特点。优点：本质边界条件容易施加，插值稳定；缺点：径 特点。 特点 优点：本质边界条件容易施加，插值稳定；缺点： 向基参数不好确定，计算效率较低， 向基参数不好确定，计算效率较低，形状函数难于满足全域的相容 性条件。 性条件。
无网格方法模拟裂纹扩展