无压渗流问题的混合自适应法

针对移动网格法和固定网格法求解无压渗流问题的不足,基于无网格法,提出了一种自由面调整的混合自适应法.该法能统一处理接近水平或垂直的自由面,且不依赖于初值,收敛性好,适用性强,结果比较精确.     

基于Kriging插值无网格法求解轴对称弹性力学问题

基于Kriging插值无网格法,提出了实际应用中复杂轴对称弹性力学问题求解的一条新途径.Kriging插值无网格法是一种新型的无网格法,该方法构造的形函数满足Kronecker delta函数性质,可以直接施加本质边界条件.采用Kriging插值无网格法分析轴对称问题,得到了轴对称问题的无网格离散方程,并编制了相应的计算程序.通过厚壁圆筒的静力学和动力学分析,对所提方法进行了检验.数值算例结果表明,提出的方法对求解轴对称弹性力学问题是行之有效的.     

自适应无网格法在生物涂层接触问题中的应用

自适应无网格法是针对有限元法无法求解或者不易求解的复杂问题,利用无网格法节点排布灵活、易于增删节点、便于自适应分析等优点发展起来的. 在对自适应无网格法理论基础和发展进行总结基础上,采用基于应变能梯度的自适应无网格—— 有限元耦合方法,对等离子喷涂制备的HA 生物涂层材料的无摩擦接触问题进行分析,对制备的两种不同厚度的生物涂层材料进行求解,分别给出了von Mises 应力分布云图. 结果表明,自适应无网格法能较好地应用于生物涂层接触问题中.     

利用无网格法进行几何非线性热固耦合柔性机构拓扑优化设计

将无网格伽辽金法引入到拓扑优化中,利用其进行了几何非线性热固耦合柔性机构的优化设计研究.利用无网格法离散和求解了热固耦合场的控制方程.基于SIMP模型和无网格法,建立了柔性机构的优化模型,利用MMA方法求解.研究了基于无网格法的伴随敏度分析方法,并提出了解决拓扑结构出现不连续散乱点同题的敏度过滤方法.求解经典算例,表明本文方法的正确性和有效性.     

无网格稳定配点法及其在弹性力学中的应用

伽辽金型无网格法具有精度高、稳定性好的优点,但是实现高阶准确积分过程复杂,计算效率低.配点型无网格法的计算效率高,但是其在求解复杂问题时往往会出现精度和稳定性较差的结果.本文介绍一种新的无网格法-无网格稳定配点法,采用重构核近似作为近似函数,在规则子域内非常容易实现高阶准确积分,既保留了配点型无网格法效率高的特点,又具备伽辽金型无网格法精度高和稳定性好的特点,而且还兼具有限体积法满足局域离散方程守恒的特点.通过弹性力学算例验证了该算法的优越性,未来可将其进一步应用于流体和流固耦合问题分析.     

无网格法的理论及应用

详细论述了近年来迅速发展的无网格法的理论基础及其在各个领域内的应 用. 无网格法网格依赖性弱, 避免了传统的有限元、边界元等基于网格的数值方法 中可能出现的网格畸变和扭曲, 在一些有限元、边界元等方法难以较好处理的领域体现 出独特的优势. 以加权余量法为主线归纳了已有的30多种无网格法, 各类 无网格法的主要区别在于使用了不同的加权余量法和近似函数. 详尽介绍 了各种无网格近似方案(包括移动最小二乘近似、核近似和重构核近似、单位分 解近似、径向基函数近似、点插值近似、自然邻接点插值近似等)和无网格法 中常用的各类加权余量法(伽辽金格式、配点格式、局部弱形式、加权最小二乘 格式和边界积分格式等), 并讨论了数值积分方法和边界条件的处理等问题. 在 此基础上较系统地总结了无网格法在冲击爆炸、裂纹传播、超大变形、结 构优化、流固耦合、生物力学和微纳米力学等领域的应用, 展示了无网格法相 对于传统数值方法的优势.     

稳定节点积分伽辽金无网格法的应力计算方法研究

应力计算是基于稳定节点积分的伽辽金无网格法的重要组成部分.该文着重研究稳定节点积分伽辽金无网格法的应力计算方法,对稳定节点积分方法的变分一致条件进行了讨论.证明当节点代表域内的应变采用非局郎光滑应变时,相应的应力在节点代表域内为常数,稳定节点积分伽辽金无网格离散方程是变分一致的.文中提出了三种节点应力计算方法,研究表明,基于位移梯度的节点应力计算方法不满足变分一致性要求,而采用光滑应变的节点应力计算方法和一致形心应力计算方法满足变分一致性要求.典型数值算例的误差分析表明,满足变分一致性不一定确保得到更为精确的结果.而基于光滑应变的一致形心应力计算方法总是较其它两种方法更为精确.     

不可压缩流场多体运动问题的两种数值解法

不可压缩流场中多体运动的数值模拟是计算流体动力学的前沿课题,目前可以求解此问题的两种高精度方法是重叠网格法和移动网格法。本文详细讨论了两种方法的理论基础,典型算法的实现过程,并采用两种方法对单圆柱和多圆柱运动问题进行了数值试验。通过分析比较,本文认为重叠网格法计算效率高于移动网格法,重叠网格法比移动网格法更适合于求解物体大位移运动问题,尤其当计算区域几何形状复杂时重叠网格法的优势更加明显。     

基于S-R和分解定理的三维几何非线性无网格法

S-R(strain-rotation)和分解定理克服了经典有限变形理论的一些缺点, 使其可以为几何非线性数值分析提供可靠的理论基础. 对于大变形问题, 由于无网格法(element-free method)避免了对单元网格的依赖, 从而从根本上避免了有限单元法(finite element method, FEM)的单元畸变问题, 保证了求解精度. 因此, 将无网格法和S-R和分解定理结合起来势必能建立一套更加合理可靠的几何非线性数值计算方法. 目前基于S-R 定理的无网格数值方法研究较少并且只能用于二维平面问题的求解, 但实际上绝大多数问题都必须以三维模型来进行处理, 因此建立适用于三维情况的S-R无网格法是非常有必要的. 本文给出了适用于三维情况的S-R 无网格法: 采用由更新拖带坐标法和势能率原理推导出来的增量变分方程, 利用基于全局弱式的无网格Galerkin 法(EFG)得到了用于求解三维空间问题的离散格式. 利用MATLAB编制三维S-R 无网格法程序, 对受均布载荷的三维悬臂梁和四边简支矩形板结构的非线性弯曲问题进行了计算. 最后将所得的数值结果与已有文献进行了比较, 验证了本文的三维S-R无网格数值算法的合理性、有效性和准确性. 本文的三维S-R无网格数值算法可以作为一种可靠的三维几何非线性数值分析方法.      

裂纹问题的一致性高阶无网格法

一致性高阶无网格法能高效精确地求解连续体问题,尤其是能得到高精度的应力场。本文将该方法拓展到应力解析精度至关重要的裂纹问题（即非连续体问题）的数值分析。采用背景积分网格描述裂纹几何,基于无需增加节点额外自由度的虚拟节点法描述裂纹处位移场的间断,提出了虚拟节点的引入算法和断裂单元的数值积分方法。为进一步模拟裂纹扩展,采用相互作用积分方法计算应力强度因子,裂纹的扩展方向由最大周向应力准则确定。数值结果表明,本文发展方法能够精确地通过间断分片试验;相较于标准的高阶无网格法和低阶一致性无网格法,本文的一致性高阶无网格法显著改善了应力强度因子的计算精度,能够准确预测裂纹扩展路径。     

无单元法研究和应用现状及动态

无单元法(也称无网格方法)是一种新兴的数值计算方法，它是 有限元等传统数值分析方法的重要补充和发展,极大地简化了前处理工 作与裂纹扩展等问题的计算分析.近年来,无单元法得到了迅速发展,受 到了国际计算力学界的高度重视.简要介绍了国内、外无单元法的发展 动态和应用现状，评述了无单元法的最新研究成果，归结出无单元法的 一些优点及目前尚有待解决的一些问题.最后指出了无单元法在工程应 用中将有着广阔的发展前景.     

断裂力学问题的杂交边界点方法

提出了一种求解断裂力学的新的边界类型无网格方法-杂交边界点法.以修正变分原理和移动最小二乘近似为基础,同时具有边界元法和无网格法的优良特性,求解时仅仅需要边界上离散点的信息.该文将杂交边界点方法应用到弹性断裂问题中,将移动最小二乘方法中的基函数扩充,能更好的模拟裂纹尖端应力场的奇异性,推导了求解断裂力学的杂交边界点法方程,与传统的元网格方法相比,文中方法具有后处理简单,计算精度高的优点.数值算例表明了该方法的稳定性和有效性.     

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

加权最小二乘无网格法

在最小二乘法和移动最小二乘近似的基础上提出了加权最小二乘无网格法．该方法除节点外又引入了一些辅助点，控制方程在所有节点和辅助点处的残差用最小二乘法予以消除，边界条件用罚函数法引入．另外对移动最小二乘近似进行了改进，并给出了最小二乘法中泛函的简化格式，因而提高了计算效率．与配点法相比，新方法精度高，稳定性好，并且系数矩阵是对称正定矩阵．与Galerkin法相比，该方法不需要进行高斯积分，因而计算量小．算例表明该方法具有效率高、精度高和稳定性好等优点，并且易于实现．     

基于核重构思想的最小二乘配点型无网格方法

介绍重构核点法的基本原理和近似函数的构造方法，并基于核重构思想，应用配点法和最小二乘原理，离散微分方程，建立求解的代数方程，提出了一种基于核重构思想的最小二乘配点型无网格方法．与一般配点法相比，该方法的系数矩阵是有对称正定的，计算精度高，稳定性好．该方法的实施不需要背景网格，不需要进行高斯积分，与Galerkin法相比，具有计算量小、边界条件处理简单的特点，是一种真正的无网格法．对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨．文中结合若干典型算例，检验了该方法的有效性．     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

基于局部Petrov-Galerkin离散方案的无网格法

基于局部Petrov-Galerkin离散方案,选用自然邻近插值构造试函数,用Shepard函数作为权函数,提出了一种无网格方法(MNNPG),这种方法充分发挥了局部Petrov-Galerkin法的优势,并且结合了自然邻近插值的特点,方便引入边界条件,由于以Shepard函数的圆形支集作为积分子域,用分片中点插值来完成区域积分,无需额外背景网格,是一种真正的无网格法。本文将该无网格方法用于求解二维弹性力学边值问题,算例结果很好地吻合了精确解,表明该方法具有良好的数值精度和稳定性。     

The analysis of composite laminated beams using a 2D interpolating meshless technique

Laminated composite materials are widely implemented in several engineering constructions. For its relative light weight, these materials are suitable for aerospace, military, marine, and automotive structural applications. To obtain safe and economical structures, the modelling analysis accuracy is highly relevant. Since meshless methods in the recent years achieved a remarkable progress in computational mechanics, the present work uses one of the most flexible and stable interpolation meshless technique available in the literature—the Radial Point Interpolation Method(RPIM).Here, a 2 D approach is considered to numerically analyse composite laminated beams. Both the meshless formulation and the equilibrium equations ruling the studied physical phenomenon are presented with detail. Several benchmark beam examples are studied and the results are compared with exact solutions available in the literature and the results obtained from a commercial finite element software. The results show the efficiency and accuracy of the proposed numeric technique.     

动力弹塑性分析的无网格自然单元法

基于无网格自然单元法,提出了结构动力弹塑性响应分析的一条新途径．自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然邻近插值的伽辽金法．自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势．在空间域上采用自然单元法离散,并运用加权余量法推导了动力弹塑性分析的离散控制方程．然后,采用预校正形式的Newmark法在时间域上进行求解．最后给出了数值算例,并验证了所提方法的有效性和正确性．     

Steady‐state solution of incompressible Navier–Stokes equations using discrete least‐squares meshless method

A fractional step method for the solution of the steady state incompressible Navier–Stokes equations is proposed in this paper in conjunction with a meshless method, named discrete least‐squares meshless (DLSM). The proposed fractional step method is a first‐order accurate scheme, named semi‐incremental fractional step method, which is a general form of the previous first‐order fractional step methods, i.e. non‐incremental and incremental schemes. One of the most important advantages of the proposed scheme is its capability to use large time step sizes for the solution of incompressible Navier–Stokes equations. DLSM method uses moving least‐squares shape functions for function approximation and discrete least‐squares technique for discretization of the governing differential equations and their boundary conditions. As there is no need for a background mesh, the DLSM method can be called a truly meshless method and enjoys symmetric and positive‐definite properties. Several numerical examples are used to demonstrate the ability and the efficiency of the proposed scheme and the discrete least‐squares meshless method. The results are shown to compare favorably with those of the previously published works. Copyright © 2010 John Wiley & Sons, Ltd.