光滑粒子动力学SPH方法应力不稳定性的一种改进方案_杨秀峰.pdfVIP

物理学报 Acta Phys. Sin. Vol. 61, No. 22 (2012) 224701 光滑粒子动力学SPH 方法应力不稳定性的 一种改进方案 杨秀峰 刘谋斌 ( 中国科学院力学研究所, 水动力学与海洋工程重点实验室, 北京 100190 ) ( 2012 年3 月16 日收到; 2012 年6 月12 日收到修改稿) 光滑粒子动力学方法是一种拉格朗日型无网格粒子方法, 在模拟大变形和自由表面流方面具有特殊的优势, 已 经在工程和科学领域得到了广泛的应用. 然而, 长期以来, 传统光滑粒子动力学方法一直受到应力不稳定性的困扰, 从而限制了它的进一步发展和应用. 应力不稳定性的根本原因在于应力状态与核函数的不匹配: 负压状态下粒子间 产生吸引力, 吸引力随着粒子间距的减小而增大, 导致拉伸不稳定性; 正压状态下粒子间产生排斥力, 排斥力随着粒 子间距的减小而先增大后减小, 导致压缩不稳定性. 本文通过改进光滑粒子动力学方法的核函数和离散格式, 使得 无论在正压还是负压状态下粒子间的作用力恒为