计算力学无网格学习报告..docVIP

无网格数值求解学习报告 无网格法的概念 当前数值模拟欧拉方程和 N-S方程的空间离散方法主要有三种:有限差分法、有限体积法和有限元法。有限差分法的主要思想是将微分型控制方程应用到每个网格节点上，对某个节点的导数值用相关网格节点值进行差分求得。再利用各种时间推进格式进行迭代求得流场的稳定解。有限体积法的主要思想是把流体力学的积分型控制方程应用到每一个由网格剖分得到的控制体微元上，从而把积分方程的求解转化为代数方程或常微分方程的求解，再利用各种时间推进格式进行迭代，最后在每个控制体微元上得到流场的稳定解。有限元方法的数值基础是变分原理和加权余量法，其主要思想是在每个网格单元上选择若干合适的点作为解函数的插值节点，方程的近似解由各网格单元的近似函数逼近，而网格中的近似函数又由单元基函数的线性组合得到。以上三种空间离散方法都是在网格的基础上进行的。 无网格方法是一种只需要节点信息而不需将节点连接成网格单元的数值解法，这是无网格算法区别于网格算法之处。因此，无网格算法不能直接沿用网格算法中的有限差分和有限体积法等的离散格式。无网格算法的基本思想是在求解域内布置一系列的离散的点称之为“节点”，然后采用一种形函数(或称为核函数)近似拟合曲线，要求每个小域(称为“点云”)内节点的个数多于建立的方程的数目，在该点云上建立矛盾方程组，通过解矛盾方程组求得变量的梯度，进而可以求得问题的解。 典型无网格法 无网格法大致可分成两类：一类是以Lagrange方法为基础的粒子法（Particle method），如光滑粒子流体动力学（Smoothed particle hydrodynamics，简称SPH）法，和在其基础上发展的运动粒子半隐式（Moving－particle semi－implicit，简称MPS）法等；另一类是以Euler方法为基础的无格子法（Gridless methods），如无格子Euler/N—S算法（Gridless Euler/Navier-Stokes solution algorithm）和无单元Galerkin法（Element free Galerkin，简称EFG）等。 特点：计算量大、精度高，稳定性好、需要背景网格进行积分、系数矩阵对称、不易施加本质边界条件处理。 背景网格积分： 加权余量法： 线弹性力学的控制方程 加权余量法不要求余量在各店均为零，而要求余量的加权积分为零 平均意义上满足方程： 等效积分形式 加权余量法的物理意义：选取合适的待定参数强迫余量在某种意义下为零 无网格方法可以方便地利用坐标点计算模拟复杂形状流场计算，但不足之处是在高雷诺数流动时提高数值计算精度较困难。 无网格方法中比较常见的还有径向基函数方法（Radious Basis Function），主要使用某径向基函数（如（MQ）f(r)=r^5）的组合，来逼近原函数。吴忠敏院士在这方面有比较突出的工作。以上方法中，无网格伽辽金法成为目前影响最大，应用最广的无网格计算方法，现有的 LS-dyna，Abaqus，Radioss等商业软件都加入了该方法的计算模块。 %-------------------------------------------------------------------------- tic % 计时开始 clear; % 清屏，将要讨论的是悬臂梁的弹性动力学问题 % DEFINE BOUNDARIES/PARAMETERS 定义边界条件和几何参数 Lx = 100; % Lx指的是问题域（矩形域）x向的边长 Ly = 10; % Ly指的是矩形域y向的边长 young = 210e9; % young是材料的弹性模量 nu=0.3; % nu是材料的泊松比 midu=10000; % midu是材料的弹性模量 %------------------------------------------------- nx = 20; % x方向的划分网格 ny = 6; % y方向的划分网格 %------------------------------------------------- ndivl=10; % x方向的积分背景划分网格 ndivw=5;