最初的SPH方法是Lucy、Gingold和Monaghan首次于1977年分别提出的[15,16],Johnson 和 Beissel提出了归一化的光滑函数算法,这一算法能够提高SPH的计算精度,并能通过分片试验。在冲击碰撞方面,Johnson 等[20,21]采用类似于SPH的广义粒子算法(generalized particle algorithm)进行圆杆的撞击以及弹体侵彻的数值计算,Parshikov等[56]用改进的SPH方法分析了子弹冲击靶体的临界穿透速度,结果与实验吻合得较好。上述冲击碰撞算例中,冲击体与靶体都是固体结构,其发生的变形相对液体-固体的冲击碰撞要小的多,并且SPH粒子均采用规则建模,对于复杂三维模型分配粒子位置和质量时存在局限性。