基于Voronoi结构的无网格局部Petrov-Galerkin方法.pdfVIP

基于Voronoi结构的无网格局部 Petrov．Galerkin方法 ∞A 蔡永昌朱合华 王建华 f同济大学地下建筑与工程系，上海 200092)(上海交通大学土木建筑工程系．上海200030) 摘要基于自然邻结点近似位移函数提出了一种用于求解弹性／U学平面问题的无网格局部Petrov-Galerkin 方法．这种方法在结构的求解域仃内任意布置离散的结点，并且利用需求结点的自然邻结点和Voronoi结构 来构造整体求解的近似位移函数．对于构造好的近似位移函数，在局部的Delaunay三角形于域上采用局部 Petrov．Galerkin方法建立整体求解的平衡控制方程，这样平衡方程的积分可在背景三角形积分网格的形心上 解析计算得到，而采用标准Galerkin方法的自然单元法需要三个数值积分点．该方法能够准确地施加边界条 件，得到的系统矩阵是带状稀疏矩阵，对软件用户来说，它还是一种完全的、真正的无网格方法．所得计算结 果表明，该方法的计算精度与有限元法四边形单元相当，但计算和形成系统平衡方程的时间比有限元法四边形 单元提高了将近一倍，是一种理想的数值求解方法． 关键词无网格，自然单元，Voronoi图，局部Petrov-Galerkin方法，Delaunay三角化 移边界上难以准确地施加边界条件等，这些缺点的 存在极大地限制了EFG法和MLPG法的应用和发 无网格方法由于只需结点而不需单元或网格信 展，使得它们虽然名义上比有限元法简单、方便，实 息，前处理特别简单，能够方便地用于裂纹开裂、 际上却无法真正应用于实际工程问题的计算，目前 非线性大变形等问题的数值计算，具有许多传统的 仅停留在理论上进行探讨． 计算方法(如有限元法)无法比拟的优点，正吸引着 无阿格方法主要是为了克服有限元法在前处理 越来越多的力学工作者和工程技术人员来对其进行 网格划分方面的困难而提出的，但是它却在计算时 研究和使用．经过十多年的发展，目前无网格方法 间、边界条件和不连续面的处理等方面付出了极为 的种类已十分繁多【““，比较典型的有光滑粒子法 昂贵的代价．实际上，在有限元法里，如果对结构分 Particle (Smooth Hydrodynamics，简记为SPH)、 析域完全采用三角形网格来进行计算，无论对于多 of 单位分解法(PartitionUnity)、有限点法(F} 么复杂的几何形状，用Delamlay准则均能全自动地 nitePoint Free Method)、无单元伽辽金法(Element 生成形状良好的三角形或四面体，根本就不会存在 Galerkia Method，简记为EFG)和无网格局部伽辽 网格划分的困难．但遗憾的是有限元法里三角形单 LocalPetrov-Galerk[n 金法(Mesbless Method，简记 元和四面体单元是常应变单元，计算精度很低，因 为MLPG)等，这些方