滑移爆轰驱动下飞板运动姿态的连续电阻测试法

飞板运动姿态的测定是爆炸焊接机理研究的基础,针对传统电测方法存在干扰因素多、易产生弯曲波等缺陷,设计了一种适用于野外大当量下爆炸焊接飞板姿态实验的连续电阻测试方法。研制了3种不同结构的梯形支架型连续电阻探针元件,利用有限元程序分析了探针的导通压力和响应时间,在此基础上,对3种探针实施了爆炸焊接实验,实验结果表明:金属丝网型探针元件具有最优的导通效果,各段测试曲线光滑无毛刺。以该探针数据计算获得了待测飞板的运动姿态曲线,并与Richter简化模型下的近似计算公式结果进行了对比,两者基本一致。所述测试方法实现了炸药爆速和飞板变形曲线的连续、可靠和快速测量,为滑移爆轰驱动问题、爆轰产物状态方程等的研究提供了测试方法补充。     

滑移爆轰密封的初步数值模拟

利用有限元方法,通过数值模拟得到爆炸密封铝、铜封管密封所需炸药药量;在此基础上,利用光滑粒子法对爆炸密封模型进行优化,分析炸药、保护外壳等形状等对封管密封性能的影响。基于对计算结果的分析,确定实验模型。实验回收结果和数值模拟结果符合较好,表明数值模拟结果对实验设计有指导作用。

     

液体炸药滑移内爆加载下钢管的变形与层裂破坏研究

用液体炸药滑移内爆加载方式,研究了不同炸药厚度和外壳条件下20钢圆管的内聚变形和破坏特征。实验中,钢管发生了较大径向应变,应变值大致随炸药厚度的增加而线性增大。在炸药厚度不低于3mm时钢管变形均匀,轴对称较好;药厚小于3mm时其变形不再是基本轴对称收缩,管子表面出现扭转褶皱。各种装药厚度下,均未观察到预期的层裂破坏。对未层裂的原因进行了初步分析。     

金属圆管在炸药滑移爆轰作用下的运动姿态

采用二维Ｌａｇｒａｎｇｅ数值摸拟方法，计算在管内装填岩石硝铰炸药滑移爆轰作用下，２０＃无缝钢管的运动过程，计算结果与用斜电阻丝法测量得到的结果基本相符。     

45钢层裂研究

对４５钢进行了层裂实验，对软回收试件进行了显微金相观测，观察了材料中微损伤的形核、增和聚合机制。将文献［３］中的动态损伤模型编入一维Ｌａｇｒａｎｇｅ有限差分动力学程序，对层裂实验进行了数值模拟，计算结果与实验吻合较好。     

内聚滑移爆轰加载下金属圆管压缩过程的工程分析

对内聚滑移爆轰加载下圆管压缩运动过程进行工程分析,分析解与数值模拟结果吻合。利用工程模型得到了不同管道参数（包括内径、厚度）下管道封闭所需要的炸药量和封闭时间表达式,可为相关研究工作提供参考。

     

爆轰产物驱动飞片运动数值模拟研究

利用数值模拟方法,研究了主装药与飞片之间的空气隙厚度变化对飞片中载荷幅度和载荷上升时间、飞片速度和飞片变形量的影响。数值模拟与验证实验结果吻合得较好。结果表明,带空气隙的化爆加载装置可以在飞片中获得较小的载荷幅度和较长的载荷上升时间,并且可以降低飞片击靶速度,从而实现较低的冲击压力,但飞片变形量随着空气隙厚度增大而增大。     

混凝土材料的层裂特性

利用Hopkinson压杆作为实验设备,通过试件后面吸收杆上应变波形研究了混凝土材料的层裂特性。不同强度混凝土在不同加载率下的实验结果表明,混凝土层裂强度与其压缩强度以及加载速率有关,给出了层裂强度和压缩强度以及加载率之间关系的经验公式,并指出不同加载速度下混凝土裂纹扩展方式不同是层裂强度率效应的主要原因。加载压缩波损伤的影响、重复加载实验以及多次层裂的顺序都表明,损伤对混凝土的层裂过程有重要影响。     

在两层不同爆速炸药滑移爆轰作用下飞板运动的数值计算研究

本文研究用低爆速炸药复盖高爆速炸药在定滑移爆轰作用下飞板运动的计算方法,并对惯性复盖与炸药复盖的结果作出相互比较。     

爆轰波管中铝粉尘爆轰的数值模拟

用两相流模型对爆轰波管中的铝粉尘的爆轰波进行了研究。模型考虑了气体和颗粒两相间速度和温度的不同及由于管壁引起的对流热传导和粘性引起的耗散,考虑了由于铝颗粒表面粗糙使得表面积增加的因素。铝颗粒的点火判据使用了新的判据,即铝颗粒在激波后的气流中温度达到铝的熔点且铝全部熔化即被点火。数值模拟了内径为15.2cm的爆轰波管中铝粉尘中爆轰波的传播和发展,得到了爆轰波速度及铝颗粒点火距离,还得到了爆轰流场中物理量的分布。从前导激波面到CJ面处,两相间的速度和温度有明显的差别。还考虑了粒子由于粗糙引起的表面积增加对爆轰波的影响,这个因素对铝颗粒的点火距离的影响较大,对这里计算的铝粉尘爆轰波的速度基本没有影响。结果表明,两相流模型可以较好地描述铝粉尘的爆轰过程,得到具有很粗糙表面、平均粒子直径为3.4m的铝粉尘浓度为304g/m3时爆轰波的速度为1.63km/s,点火距离为3mm,与实验值符合较好。     

爆轰波对碰驱动下金属圆管膨胀变形特性研究

采用高速分幅摄影和激光速度干涉仪（velocity interferometer system for any reflector, VISAR）联合测试技术,对HR-2钢管在爆轰波对碰驱动下的膨胀变形特性进行了研究。实验得到了圆管膨胀变形的过程图像以及圆管特定部位径向膨胀的位移-时间曲线和速度-时间曲线。采用LSDYNA-3D程序对HR-2钢管的膨胀变形过程进行了数值模拟,计算结果与实验结果吻合较好。     

爆轰波对碰区产物驱动金属圆管的研究

利用闪光X光照相、光学高速分幅照相实验,测得了金属圆管在两个爆轰波对碰后产物驱动外壳的膨胀、变形过程,该过程表明爆轰波对碰作用使得对碰处金属圆管相对贯穿断裂时间明显提前,对碰处金属圆管径向膨胀速度较非对碰处明显增加。利用拉格朗日二维流体动力学程序TTD2C,数值模拟了对碰驱动过程,计算结果与实验结果符合很好,且用Taylor断裂判据得到的金属圆管相对贯穿断裂时间也明显提前。     

环形方管道中爆轰胞格的三维研究

通过实验和三维数值模拟研究了爆轰波在环形管道内的传播。实验采用烟迹板记录了爆轰波的胞格结构。数值模拟基于带化学反应的三维Euler方程,采用五阶精度的WENO格式捕捉激波,采用具有TVD性质的三阶Runge-Kutta法处理时间项,并结合并行技术,对爆轰波的传播进行了数值研究。结果表明,环形管道外壁为收敛壁面,由于其对流场的压缩效应,外壁面及附近的胞格较小,且较均匀。而内壁为发散壁面,其对流场起稀疏效应,内壁面及附近的胞格较大,且呈周期性变化。同时, 不同壁面的胞格结构均出现了拍波(slapping wave),其形状呈弯曲的折线。     

一点起爆半球形装药装置中散心爆轰波传播的数值模拟研究

运用自行研制的二维流体动力学程序TDY2D,对一点起爆半球壳装药装置中散心爆轰波的传播过程进行数值模拟计算.对冲击波到达炸药和飞层各界面的波形的计算值与实验值进行了对比,两者符合较好,从而验证了该程序的准确性和有效性;对爆轰波的传播过程进行了分析研究,获得了一些散心爆轰波传播的规律性认识,为进一步研究散心爆轰波特性提供了数值分析基础.     

爆轰波拐角传播三维数值模拟

为了研究爆轰波在拐角中的传播特性以及拐角爆轰低压流场现象,用LSDYNA 3D程序对三种常见的拐角装药的爆轰波传播特性进行了数值模拟,清晰显示了爆轰波因拐角处爆轰波传播面积的变化而产生的衰减 增长过程,讨论了侧向三通中爆速的变化情况,并得到了拐角处低压爆轰区的尺寸。     

单次脉冲爆轰发动机工作过程的数值模拟

通过求解14组分和19个基元反应的CH4-O2-N2详细化学反应动力学机理的二维轴对称Navier-Stokes方程,对爆轰管内半球形高温火团引发的爆轰过程和爆轰波进入外流场后的全流场分布进行数值模拟。模拟结果显示了爆轰波在管内的成长、稳定传播、进入尾喷管后衰减为激波和进入外流场的全过程,以及爆轰管出口端附近区域的复杂涡与激波的相互作用。对轴线上的压力分布和封闭端的压力等进行了讨论,为脉冲爆轰发动机的开发研制提供参考信息。     

直管内胞格爆轰的基元反应数值研究

基于基元反应和二维欧拉方程,对直管内胞格爆轰进行了数值模拟。采用5阶WENO(weighted essentially nonoscillatory scheme)求解对流项,采用2阶附加半隐的龙格-库塔法处理化学反应源相引起的刚性。获得了密度、压力、温度和典型组元质量分数流场及数值胞格结构等。结果表明：网格精度的差异明显影响胞格的规则性和爆轰的平衡模数,随着网格尺度的减小,胞格由不规则变为规则。预混气组成、初压、初温及管道宽度给定,三波点数收敛为确定值。足够强度的初始扰动可再现胞格爆轰,最终形成的自持胞格爆轰模数与初始扰动的形状、大小、位置均无关。沿胞格中心线,爆轰波速度变化范围为0.88DCJ~1.5DCJ,爆轰波平均速度与CJ爆轰速度仅偏差0.88％。峰值压力与初压之比为14~50。计算爆轰波平均速度、胞格宽长比与实验值基本一致,但计算胞格宽度比实验值略小。数值模拟加深了对横波的产生和发展、未反应气囊、爆轰胞格的二次起爆等胞格爆轰特性的认识。