SPH方法在聚能装药射流三维数值模拟中的应用  

Application of SPH method to numerical simulation of shaped charge jet

作　　者：李磊[1] 沈兆武[1] 李学岭 倪小军[1] 

机构地区：[1]中国科学技术大学近代力学系,安徽合肥230027 [2]安徽方圆机电股份有限公司,安徽蚌埠233010

出　　处：《爆炸与冲击》2012年第3期316-322,共7页Explosion and Shock Waves

基　　金：国家自然科学基金项目(51134012,51174183)~~

摘　　要：采用有限元软件LS-DYNA中的SPH算法实现聚能装药射流形成过程的三维数值模拟。将SPH方法和FEM方法与理论计算结果相比较。研究结果表明:SPH数值模拟计算过程十分稳定,避免了有限元数值模拟过程中的网格扭曲、缠绕和物质穿透等问题,而且计算精度比有限元方法更高;采用SPH方法计算的射流速度、射流长度与有限元计算结果和理论计算结果的误差均小于5%,因此该方法对于模拟诸如聚能装药等多介质、大变形问题十分有效;初始相邻粒子数N的最佳取值范围为600～1 800,在该范围内计算结果的误差小于3%,且计算效率高。Based on the smoothed particle hydrodynamics（SPH） method in LS-DYNA software,the shaped charge jet process was simulated.And the mechanism of the shaped charge jet was discussed according to the numerical simulation results.To analyze the accuracy,effectiveness and superiority about this method,as well as the error brought by different initial numbers N of touching neighbors per particle,the numerical simulation results were compared with two-dimensional（2D） finite element results and theoretical calculations.The investigated results show that some problems,such as grid distortion,winding and contact penetration which often happen in the finite element method,can be avoided in the SPH method,so the calculation process by the SPH method is stable,and the calculation precision is higher than that by the 2D finite element method.The relative error of the jet tip velocity and the jet length calculated by the SPH method compared with the finite element results and theoretical calculations are both lower than five percent,therefore the SPH method is adaptable in multi-material,large deformation problems such as shaped charges.The most ideal N is from 600 to 1 800,in which the relative error is lower than 3% and the computational efficiency is high.

关 键 词：流体力学 初始相邻粒子数 聚能射流 SPH 

分 类 号：O358[理学—流体力学]