改进的泰勒公式及其数值验证.pdf

　第 26 卷 　第 4 期 爆 炸 与 冲 击 Vol . 26 , No . 4 　 　2006 年 7 月 EXPL O SION AND SHOC K WAV ES J ul . , 2006 　 文章编号 : 100 11455 (2006) 改进的泰勒公式及其数值验证 刘 　义 ,秦承森 ,杭义洪 ,王 　裴 (北京应用物理与计算数学研究所 ,北京 100088) ( ) 　　摘要 : 将光滑粒子法 smoot hed p article hydrodynamics met ho d , SP H 应用于金属平板爆轰驱动过程的 数值模拟 ,得到了炸药与金属相互作用的清晰的图像和平板抛射角的变化 。根据抛射角方程对泰勒公式进行 了修正 ,得到了改进的泰勒公式 。利用数值模拟得到的结果比较了平板驱动几种不同的抛射角公式 ,改进的 泰勒公式准确反映出了抛射角下降并趋于终态的规律 ,其对平板抛射角的描述更为准确 。 　　关键词 : 爆炸力学 ;抛射角公式 ; SP H ;泰勒公式 ;爆轰驱动 　　中图分类号 : O383 　　　国标学科代码 : 130 ·3520 　　　文献标志码 : A 1 　引　言 　　爆轰过程中产生的高压能够驱动聚能罩及飞片使其加速和变形 ,这在爆炸加工 、爆炸焊接 、成型装 药设计等许多领域内都有广泛应用 。爆轰波的传播过程是一个非常复杂的物理化学过程 ,爆轰波与介 质尤其是金属介质的相互作用一直以来都为人们所关注 。在爆轰波后压力极高 ,远大于金属的屈服强 度 ,此时金属呈现流体的特征 ,在计算中通常将其视为流体 , 因此爆轰产物与金属平板的作用可以处理 为两种不同流体介质间的相互作用 ,使用统一的流体动力学方程组来进行求解 。 爆轰驱动飞片的飞行速度和飞行角度及其变化规律是最被关注的问题[ 1～3 ] 。爆轰对于二维平板的 ( ) 驱动问题是其中最简单的情形 ,也是颇具代表性的一种情况 。本文中应用光滑粒子法 SP H 计算金属 平板在爆轰波驱动下的变形及飞行过程 ,得到作用过程中清晰的波系图像和抛射角的变化过程 。 由于冲击波的瞬时作用及平行板方向和垂直板方向的速度分量的变化 ,平板抛射角初始较大 ,然后 [4 ] 下降为终态稳定值 。理论公式尚不能描述这一实际情况 ,本文中通过抛射角方程 对泰勒公式进行修 正 ,得到改进的泰勒公式 ,可以描述抛射角的这一变化 。通过计算得到的结果对几种抛射角公式进行比 较 ,改进的泰勒公式与数值计算结果一致性较好 。 2 　爆轰驱动平板的数值计算 　　爆轰对平板的驱动是炸药对物质作用的一种 简单情形 。这个问题在成型装药 、爆炸加工 、爆炸 焊接中都有着广泛的应用 。 采用光滑粒子法计算如图 1 所示的二维平面 三明治结构的爆轰作用过程 。 图 1 平板三明治结构 光滑粒子法是一种拉格 朗 日方法 , 最早在 Fig . 1 The flat sandwich configuration 19 7 7 年 由L ucy 等提 出 。它最早被应用于天体物 理中非轴对称现象的计算中 ,是一种无网格粒子方法 ,可以方便地由低维扩展到高维[ 5～7 ] 。光滑粒子法 的基本思想是积分插值 ,通过积分插值流体动力学微分方程