应力波数值计算中的SPH方法  

Application of SPH in stress wave simulation

作　　者：孙晓旺[1] 章杰[1] 王肖钧[1] 李永池[1] 赵凯[1,2] 

机构地区：[1]中国科学技术大学近代力学系,安徽合肥230026 [2]北京理工大学科学与技术国家重点实验室,北京100081

出　　处：《爆炸与冲击》2017年第1期10-14,共5页Explosion and Shock Waves

基　　金：国家自然科学基金项目(11402266;11202206;11472008);爆炸科学与技术国家重点实验室开放基金项目(KFJJ13-9M);爆炸冲击防灾减灾国家重点实验室开放基金项目(DPMEIKF201401)

摘　　要：对一维波动方程的SPH(smoothed particle hydrodynamics)格式和有限差分格式进行比较,并采用SPH法模拟了一维应力/应变波,获得1个可衡量SPH法模拟应力波准确性的重要指标。结果表明,SPH法模拟应力波传播中采用的光滑长度必须不小于粒子间距;采用B-样条核函数和高斯型核函数能够获得良好的应力波图像,而二次型核函数不能,因此二次型核函数不适用于冲击动力学的数值计算。Obtaining accurate waveforms is significant in impact mechanics numerical calculation. This paper is to analyze how the kernel functions and smooth length affect the result of stress wave simulation. The SPH （smoothed particle hydrodynamics） formulations with different kernel functions and smooth lengths of one dimensional wave equation was compared with the finite difference formulation, which was derived in this paper. One dimensional stress and strain waves were simulated using the SPH method with different kernel functions and smooth lengths, and waveforms were gained accurately by B-spline and Gaussian kernels when the smooth length was equal to or greater than the parti- cle interval. The wave velocity obtained by the quadratic kernel is below the theoretical value, no matter what the smooth length is. A parameter was deduced in this paper as roughly equal to the di- mensionless wave velocity. Several conclusions were drawn. Firstly, the smooth length is equal to or greater than the particle interval, which is the necessary prerequisite for accurate stress wave simulation with SPH. Then, the quadratic kernel is not suitable in impact mechanics numerical calculation. Finally, the parameter deduced in this paper is a significant index to evaluate the stress wave simula- tion result of SPH.

关 键 词：爆炸力学 核函数 光滑长度 SPH方法 应力波 

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