1-19 滑移爆轰作用下飞板运动的解析解谭多望 孙承纬 根据在一维化拟定态运动的参考系中的瞬时爆轰产物及其稀疏流动,建立滑移爆轰驱动飞板运动的理论模型,导出了任何k(炸药多方指数)值下滑移爆轰驱动飞板运动的计算公式,这种简单明了的理论模型可较好地再现目前有代表性的实验和模型计算结果 。考虑飞板中某一微元,t=0时爆轰波阵面掠过,使得微元表面炸药段发生瞬时爆轰。以该微元与炸药接触的表面为原点建立(x,t)坐标系,如图1所示。
在忽略空气阻力、飞板的强度和可压缩性、应力波在飞板中的传播以及假定爆轰产物作等熵流动、爆轰波到达飞板表面时反射弱冲击波的条件下,产物在(x,t)平面上的流动分为Ⅰ、Ⅱ、Ⅲ、Ⅳ四个区域,Ⅰ为产物静止区,Ⅲ为右行波C+波后区,Ⅱ、Ⅳ分别是微元运动同Ⅰ、Ⅲ相互作用区域。
由于假定微元内表面反射稀疏波,其前沿迹线SR必是特征线族C-之一支。当k=3时,产物运动的右行波C+不受微元运动引起的左行波C-的影响,黎曼不变量 因此Ⅲ区和Ⅳ区的流动可用以-1点为心的中心稀疏波描述,它们的差别在于Ⅲ区的左行波为均匀的平行直特征线,Ⅳ区的左行波为飞板运动发出的虚拟中心稀疏波。孙承纬已证明从飞板迹线 记飞板厚度为lm、密度为rm,炸药初始压力为ps=[2(k+1)]-1,初始声速为cs=[k/2(k+1)]1/2,并令h=ps/(lmpmcs3), a =2(1+4h cs+2h2 cs3)/(2+h cs)2, b=(2+4h cs+h 2 cs3)/[ h (2+h cs)]。根据黎曼不变量a 的函数及沿特征线的量不变等关系,可导出当k=3时飞板运动速度公式(1)。 当炸药多方指数
(1)
通 过与实验结果、数值模拟和其它方法的比较,在参数k和 |