无网格方法.ppt
合集下载
无网格法介绍
无网格法(MFree)简介
1
目
录
1
无网格法概述 无网格法分类 构造无网格形函数 导出无网格法公式
2
3 4 5
无网格法研究主要进展及参考文献
2
无网格法概述
无网格法定义
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)
20
导出无网格法公式
基于局部弱式的无网格法
MFree局部弱式法核心是对控制方程在每一个局部子域采用局部加权 残值法，使一个需要在全域求解的问题简化为在各个子域上对局部 伽辽金方程的求解问题,从而避免了全域的数值积分 代表方法：无网格局部伽辽金法(MLPG)，局部点插值法(LPIM) MLPG特点。优点：避免全局数值积分，减小了对背景积分网格依赖； 缺点： “刚度阵”带状但不对称，增加了计算难度，尽管在大部分 的边界积分可通过选用适当的权函数消除掉 ，但在问题域的边界上 或附近,边界积分是不可避免的，使其难以应用于复杂边界。 LPIM特点。优点：更易满足本质边界条件且有较好 的精度和收敛性； 不需满足全域相容性，更有应用前景；缺点：插值力矩阵容易发生 奇异，需要特殊处理。
在移动最小二乘近似（MLS）中，系数a(x)的选取使近 似函数 u h ( x) 在计算点x的邻域内待求函数u(x)在某种最 小二乘意义下的最佳近似。近似函数在节点 xi 处的误差 数。
1
目
录
1
无网格法概述 无网格法分类 构造无网格形函数 导出无网格法公式
2
3 4 5
无网格法研究主要进展及参考文献
2
无网格法概述
无网格法定义
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)
20
导出无网格法公式
基于局部弱式的无网格法
MFree局部弱式法核心是对控制方程在每一个局部子域采用局部加权 残值法，使一个需要在全域求解的问题简化为在各个子域上对局部 伽辽金方程的求解问题,从而避免了全域的数值积分 代表方法：无网格局部伽辽金法(MLPG)，局部点插值法(LPIM) MLPG特点。优点：避免全局数值积分，减小了对背景积分网格依赖； 缺点： “刚度阵”带状但不对称，增加了计算难度，尽管在大部分 的边界积分可通过选用适当的权函数消除掉 ，但在问题域的边界上 或附近,边界积分是不可避免的，使其难以应用于复杂边界。 LPIM特点。优点：更易满足本质边界条件且有较好 的精度和收敛性； 不需满足全域相容性，更有应用前景；缺点：插值力矩阵容易发生 奇异，需要特殊处理。
在移动最小二乘近似（MLS）中，系数a(x)的选取使近 似函数 u h ( x) 在计算点x的邻域内待求函数u(x)在某种最 小二乘意义下的最佳近似。近似函数在节点 xi 处的误差 数。
无网格方法(刘欣著)PPT模板
5.6.1界面问 题的增强函 数
5.6.3数值 计算
07
ONE
第6章有限点方法
第6章有限点方法
6.1对流-扩散方程的有限点形式 6.2对流-扩散方程的有限点法求解 6.3Burgers方程的高阶时间格式有限点方法求解 6.4油藏数模的有限点法 6.5有限点方法在金融工程中的应用
第6章有限点 方法
10
ONE
第9章流体-结构相互作用的无网 格方法研究进展
第9章流体-结构相互作用的无网 格方法研究进展
9.1流体-结构相互作用的计算研究 概述 9.2流体-结构相互作用模型描述 9.3FSI问题的扩展有限元方法求解 9.4浸入粒子方法 9.5气动弹性计算中的径向基函数法
第9章流体-结构相互作用的无网格方法研究进展
5.4增强型单位分解有限元方法
5.4.1增强 型覆盖函数 的实现
5.4.2数值 计算
第5章单位分解 有限元方法
5.5单位分解有限元在断裂力学中 的应用
1
5.5.1裂纹尖端附近的渐近解
2
5.5.2平面裂纹的单位分解有限 元计算
第5章单位分 解有限元方法
5.6单位分解有限元在界面问题中 的应用
5.6.2界面问 题的增强方 式
7.6.3数值求解
09
ONE
第8章自适应无网格方法
第8章自适应无网 格方法
8.1自适应无网格Galerkin法 8.2结构动力问题的自适应无网 格计算 8.3hp自适应无网格方法
第8章自适应无网格方法
8.1自适应无网格Galerkin法
8.1.1后验误差估计
8.1.2背景网格重构 算法
8.1.3自适无网格静 力分析
无网格方法(刘欣著 )
无网格法的理论及应用
为了验证该方法的有效性和可行性，我们进行了一系列实验。实验过程中采 用了某稠油油田的实际数据集，包括地层压力、温度、渗透率等参数。同时，采 用了可视化评估指标，以便直观地评估计算结果的准确性。实验结果表明，该方 法在稠油热采数值模拟过程中具有较高的计算精度和计算效率，可为稠油热采技 术的优化提供有力支持。
1、算法开发：针对稠油热采的物理化学过程，开发相应的数值模拟算法， 如有限元法、有限差分法等。
2、软件架构：设计并实现数值模拟软件的架构，包括前后处理、求解器等 模块，以便用户进行快速高效的计算。
3、数据处理：针对稠油热采数值模拟过程中产生的大量数据，开发相应的 数据处理技术，如数据压缩、可视化等。
无网格法的数值积分采用移动最小二乘法（Moving Least Squares，MLS） 来实现。该方法通过对节点进行加权，构造一个局部近似函数来逼近真实的解。 数值积分通过在节点上建立局部近似函数，然后对该函数进行求导和积分来计算。 无网格法的数值积分具有高精度和高效性，同时避免了传统网格法中的网格生成 和数据处理问题。
1、结构分析
无网格法在结构分析中具有广泛的应用，可以处理各种复杂形状和材料属性 的结构。例如，桥梁、建筑物和飞机等结构分析中，无网格法能够适应复杂的几 何形状和非均匀的材料属性，同时提高计算效率和精度。此外，无网格法在疲劳 分析和振动分析中也得到了广泛应用。
2、流体分析
无网格法在流体分析中也有着广泛的应用，可以处理各种复杂的流体流动问 题。例如，无网格法可以应用于计算流体动力学（CFD）中的复杂流场模拟、燃 烧模拟以及噪声辐射模拟等。无网格法能够适应复杂的几何形状和流场特性，提 高计算精度和效率。
参考内容
稠油热采是一种重要的石油开采方法，具有提高采收率、降低开采成本等优 势。随着计算机技术的不断发展，数值模拟已成为稠油热采领域的重要工具。本 次演示旨在探讨稠油热采数值模拟自适应网格法计算软件的开发研究及实例应用。
无网格法的应用
无网格法的应用无网格方法的研究应用与进展引言有限元法（FEA）是随着电子计算机的发展而迅速发展起来的一种现代计算方法，但FEA 是基于网格的数值方法，在分析涉及特大变形（如加工成型、高速碰撞、流固耦合）、奇异性或裂纹动态扩展等问题时遇到了许多困难。
同时，复杂的三维结构的网格生成和重分也是相当困难和费时的。
近年来，无网格得到了迅速的发展，受到了国际力学界的高度重视。
与有限元的显著特点是无网格法不需要划分网格，只需要具体的节点信息，采用一种权函数（或核函数）有关的近似，用权函数表征节点信息。
克服了有限元对网格的依赖性，在涉及网格畸变、网格移动等问题中显示出明显的优势。
无网格方法的概述无网格方法（Meshless Method）是为有效解决有限元法在数值模拟分析时网格带来的重大问题而产生的，其基本思想是将有限元法中的网格结构去除，完全用一系列的节点排列来代之，摆脱了网格的初始化和网格重构对问题的束缚，保证了求解的精度[1]。
是一种很有发展的数值模拟分析方法。
目前发展的无网格方法有：光滑质点流体动力学法（SPH）、无网格枷辽金法（EFGM）、无网格局部枷辽金法（MLPGM）、扩散单元法（DEM）、Hp－clouds 无网格方法；有限点法（FPM）、无网格局部Petrov－Galerkin方法（MLPG）、多尺度重构核粒子方法（MRKP）、小波粒子方法（WPM）、径向基函数法（RBF）、无网格有限元法（MPFEM）、边界积分方程的无网格方法等。
这些方法的基本思想都是在问题域内布置一系列的离散节点，然后采用一种与权函数或核函数有关的近似，使得某个域上的节点可以影响研究对象上的任何一点的力学特性，进而求得问题的解。
无网格方法国内外研究的进展无网格法起源于20 世纪70 年代。
Perrone，Kao 最早采用任意网格技术将传统有限差分进行扩展，提出了有限差分法，这可看作无网格技术的最初萌芽。
1977年Lucy 和Monaghan 首次提出了基于拉格朗日公式的光滑质点流体动力法（Smoothed Particle Hydrocynamics：SPH），这是一种纯拉格朗日法，无需网格。
无网格法介绍
无网格法是在建立问题域的系统代数方程时，不需要利用预定义的 无网格法是在建立问题域的系统代数方程时， 网格信息，或者只利用更容易生成的更灵活、 网格信息，或者只利用更容易生成的更灵活、更自由的网格进行域 离散的方法。（刘桂荣，2009） 离散的方法。（刘桂荣，2009） 。（刘桂荣
无网格法概述
无网格法求解过程 FEM对比 对比） （与FEM对比）
导出无网格法公式
基于弱强式的无网格法
• MFree弱-强式法 弱 强式法 强式法(NWS)的核心思想是针对某一问题同时采用强式和 的核心思想是针对某一问题同时采用强式和 局部弱式建立起离散系统方程式，即对不同组别的节点根据其不同 局部弱式建立起离散系统方程式， 条件分别形成不同类型的方程，其中局部弱式被用于位于或接近导 条件分别形成不同类型的方程， 数边界条件的所有节点，强式被用于除此之外的其他节点。 数边界条件的所有节点，强式被用于除此之外的其他节点。 • 代表方法：MWS 代表方法： • MWS特点。MWS法使用最少数量的背景网格用于积分，对各类力学 特点。 法使用最少数量的背景网格用于积分， 特点 法使用最少数量的背景网格用于积分 问题均可得到稳定而精确的解，是目前近乎理想的无网格法。 问题均可得到稳定而精确的解，是目前近乎理想的无网格法。
构造无网格形函数
PIM形函数性质
• 一致性 如果单项式的完备阶数是p，则该形函数具有 C p 一致性 如果单项式的完备阶数是 ， • 再生性 PIM基函数可再生包含在其基函数当中的任意函数。 基函数可再生包含在其基函数当中的任意函数。 基函数可再生包含在其基函数当中的任意函数 • 线形独立性 PIM基函数在支持域上是线性独立的 基函数在支持域上是线性独立的 • δ 函数性
目
第六章 LS-DYNA无网格法
登录协同仿真时代
登录协同仿真时代
谢 谢！
Welcome to
ANSYS
China@
登录协同仿真时代
登录协同仿真时代
登录协同仿真时代
登录协同仿真时代
创建了材料之后，选择 done 退出当前窗口，再选择 done 结束 材料的创建。 在sph的主窗口下，选择apply，可以看到我们创建的box形的 sph模型出现在窗口中。
ቤተ መጻሕፍቲ ባይዱ
登录协同仿真时代
演示
登录协同仿真时代
登录协同仿真时代
登录协同仿真时代
登录协同仿真时代
从计算角度看，该方法是将连续的物质表示为带有速 度的可运动的粒子的集合。每一个粒子均可以表示为其它 粒子的插值点。如粒子“I”与其相距一设定距离(通常为 2h)范围内的所有其它粒子“J”发生相互作用。它们间的 相互作用是由称为光滑函数(或核函数)的近似函数W(xx′,h )来衡量的，h称为光滑长度。因此，任意粒子“I”的 连续函数的值或其导数都可以利用周围粒子“J”的已知值 估计出来。这样，整个问题的解就转化为采用规则的插值 函数，对所有粒子进行插值计算的问题，守恒方程可以转 化为采用流量或内力表示的形式。
登录协同仿真时代
登录协同仿真时代
临近搜索
登录协同仿真时代
SPH计算步骤
登录协同仿真时代
建模过程
在ls-prepost中第7页有sphgen 的功能，利用这个功能可以 建立一些简单的无网格模型。 选择其中的new，表示新建一 个无网格的命令。可以看到 下面有几个选项：box， sphere，cylinder，sketch
登录协同仿真时代
选择box，出现如图2所示的菜单，即为创建一个具有长方体形 的无网格模型。
无网格讲义
令
和 是
得
上完全单调的。
CPD径向函数和多次单调函数
定理9.3
是正整数，设 不是至多m次的多项式，若 在 上 次单调且不是常数，那么 是 中m阶严格正定的且径向的。任意 s，使 得
类似前面，我们知道具有紧支集的径向 函数不是 (对所有s)中的严格正定函数, 同样，没有具有紧支集的条件正定径向 函数。 定理9.4： 若函数 有紧支集， 是m阶严格条件正定的，则m=0；也就 是， 是严格正定的。
其广义Fourier变换是 阶m=max(0, ) 其中 ：大于等于 的最小整数 ：修正的 阶第二类Bessel函数 注：应排除 ，因为这样的 会使多项式的 次数是偶数
由广义Fourier函数是正定的且在原点具m阶奇 异性，据TH7.3知 是 (or> ) 阶严格条件正定函数。 注：若 ，Fourier变换是一种古典的变换， 由第四章的广义逆multiquadrics，可知函数是 0阶严格条件正定的，也就是严格正定的。
特别的，可用薄板样条函数
和约束条件
求解 中的散乱数据插值问题，这样由线性 函数得到的插值结果是精确的
注：古典薄板样条（ ）是2阶严格条件正定的； 是 是3阶严格条件正定的 这些函数不是单调的 这两个图都包含负函数值的部分
径向幂和薄板样条通常统称为调和样条 注： 还没有理论可表明使用薄板样条插值的 解是适定的
定且迹非正，则 A有一个负的特征值及N-1个正的特征值。 证明：令A的特征值 由TH9.5知 所以A至少有N-1个正特征值 而 因此至少有一个负的特征值
定理9.7
假设 是一阶严格条件正定的且 (0) 0 对于任意不同的点 ，矩阵 A（ ）有N-1个正特征值和1个负特 征值，进一步，A是非奇异的。 证明： 显然A是一阶严格条件正定的； A的迹 ，由TH9.6立得结论。
计算固体力学(有限元以及无网格方法)全套教学【121P】PPT课件
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x O
三角形单元
将位移试函数代入上式，并求偏导数，得
xxyy222111 (((bcciiiuuviii
bjuj cjvj cjuj
bmum) cmvm) cmum)(bivi
bjvj
bmvm)
第二章 平面弹性力学的有限元法
反映了单元的位移形态，称为形函数
vm
m (xm, ym)
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x
三角形单元
同理有 vN iv i N jvj N m v m N kv k
则位移向量可表示为
i,j,m
{ } e 单元节点位移向量
ui
vi
{f
}
u v
Ni
0
0 Ni
求
L(u)0
解 域
u aiui
离 散
i
L'(ui) 0
AXB
各种数值方法
ui u(xi)离散节点的变量值
第一章 科学和工程中的数值方法
1.3 几个简单示例
(a) 开孔板力学模型
(b) 力学模型离散化
平面问题有限元法
第一章 科学和工程中的数值方法
BEM的变形
起重机吊钩
FEM的变形
第一章 科学和工程中的数值方法
2.2 三角形常应变单元
y
3 单元中的应变和应力
{}[B]{}e
由于[B]是常量，单元内各点应变分
量也都是常量，这是由于采用了线性位移 O 函数的缘故，这种单元称为三角形常应变 单元。
i(xi , yi ) u i
um vj
uj j(xj , yj )
x O
三角形单元
将位移试函数代入上式，并求偏导数，得
xxyy222111 (((bcciiiuuviii
bjuj cjvj cjuj
bmum) cmvm) cmum)(bivi
bjvj
bmvm)
第二章 平面弹性力学的有限元法
反映了单元的位移形态，称为形函数
vm
m (xm, ym)
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x
三角形单元
同理有 vN iv i N jvj N m v m N kv k
则位移向量可表示为
i,j,m
{ } e 单元节点位移向量
ui
vi
{f
}
u v
Ni
0
0 Ni
求
L(u)0
解 域
u aiui
离 散
i
L'(ui) 0
AXB
各种数值方法
ui u(xi)离散节点的变量值
第一章 科学和工程中的数值方法
1.3 几个简单示例
(a) 开孔板力学模型
(b) 力学模型离散化
平面问题有限元法
第一章 科学和工程中的数值方法
BEM的变形
起重机吊钩
FEM的变形
第一章 科学和工程中的数值方法
2.2 三角形常应变单元
y
3 单元中的应变和应力
{}[B]{}e
由于[B]是常量，单元内各点应变分
量也都是常量，这是由于采用了线性位移 O 函数的缘故，这种单元称为三角形常应变 单元。