再生核质点法的研究进展.pdfVIP

维普资讯 第32卷 第4期 力 学 进 展 Vo1．32 No．4 2002年 l1月25日 ADVANCESIN MECHANICS Nov．25， 2002 再生核质点法研究进展 周进雄 李梅娥 - 张红艳 张 陵 西安交通大学建筑工程与力学学院，西安 710049 西安交通大学机械工程学院，西安 710049 摘 要 概述了无网格再生核质点法 的近似方案和本质边界条件的处理；综述了再生核质点法在结 构动力学、大变形分析、计算流体力学、断裂及损伤力学、微机电结构分析中的应用；从工程应用和 改进再生核质点法的角度提出了未来的若干研究方向． 关键词 再生核质点法，无网格法，多分辩率分析，大变形，结构功力学，损伤， CFD 大量复杂的工程实际问题为计算力学提出了许多迫切需要解决的难题．例如在模拟材料挤 压或模塑成形过程中，需要处理非常显著的结构大变形问题；在铸造充型过程及金属微观组织模 拟过程中，需要精确地计算出自由表面及液固相边界的移动和变化；在结构破坏失效预测中， 模拟任意复杂裂纹的扩展过程至关重要．传统的依赖于网格的有限元法在处理大变形 问题时经 常 由于网格缠结 (meshentanglement)而导致求解失败，而且局部应力集 中等现象的精细分析必 须进行网格细化并反复迭代求解；在求解裂纹扩展等移动不连续问题时，每求解一步必须将网 格重新剖分．这使得通常的有限元法在处理这一类 问题时不仅要花费大量的时间，而且求解过 程非常繁琐且计算精度较差．基于上述原因，一种部分或完全抛弃传统有限元网格的近似方法 — — 无网格法 (meshlessmethod)近几年来引起 了国内外学者的广泛关注．无网格法无需计算 网格，可以避免大变形分析网格畸变而引起的计算困难；对感兴趣的区域进行精细分析只需在 相应区域加入点即可，易于进行 自适应分析；在计算移动不连续问题时，无需网格重构；选择光 滑连续的核函数或权函数，所得到的解及其各阶导数也是全局连续的．从而使无网格法在处理 移动不连续、大变形、高梯度 问题等方面 比基于网格的近似方法具有特殊的优越性． 尽管无网格法已经有 20多年的历史了，但直到最近几年才真正引起人们的关注并开展了深 入的研究，涌现 了各种各样的无网格方法 _lJ．其中主要有：辐射元法 (disuseelementmethod， DEM)，无网格 Galerkin法 (elementfreeGalerkin，EFG)，有限点法 (finitepointmethod，FPM)， Hp云团法 (Hpclouds)，无网格局部边界积分方程法 (meshlesslocalboundaryintegralequation， MLBIE)，无网格局部 PetrovGalerkin法 (meshfreelocalpetrovGalerkin，MLPG)，单位分解 法 (partitionofunitymethod，PUM)，光滑质点流体动力学法 (smoothpraticlehydrodynamics， SPH)和再生核质点法 (reproducingkernelparticlemethod，RKPM)． 收稿 日期 ：2001．O1．O4，修回日期 ：2002—05-13 国家自然科学基金 和教育部高等学校骨干教师计划资助项 目 · 535 · 维普资讯 RKPM 是在 SPH法 的基础上发展起来的． SPH法是最早 出现 的无网格方法，它是一种纯 Lagrangian方法，主要用来处理无边界的天体物理 问题．但离散形式的SPH法对于非均匀布置 的质点不满足线性相容性条件，即使对于均匀布置的质点在边界处亦不满足线性相容性条件 -lJ． 因此 SPH法在处理有限域问题或仅采用少量离散质点时计算精度较差，相容性条件的丧失使其 在边界处变得很不稳定 J．针对 SPH法的不足，美 国西北大学计算力学学者 W ．K．Liu通过 对 SPH法核函数乘 以一项修正函数以满