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

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

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

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

     

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

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

滑移爆轰驱动钢管的层裂

将VG损伤模型推广到二维情况,考虑了最大主应力方向对损伤演化的影响,并使用显式断裂算法对20钢管在GI-920炸药滑移爆轰驱动下层裂的问题进行二维数值模拟。分析了一维内爆和滑移爆轰两种加载方式下作用于钢管外表面的压力及钢管内部受力状态的区别,考察了滑移爆轰加载方式下钢管外表面的受力随炸药厚度变化的规律,进而研究了钢管内损伤的分布和演化,以及裂纹的产生和扩展现象。计算得到的层裂片初始厚度随炸药厚度的变化规律与实验结果符合较好。     

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

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

金属圆管内爆轰波相互作用效应的数值模拟研究

采用动力学有限元程序,针对大变形带来的网格扭曲造成计算的锁死现象进行了ALE技术处理,结合高速摄影和X光照相实验结果,对金属圆管内爆轰波的相互作用效应进行数值模拟,给出了爆轰波对碰作用下金属圆管运动计算结果,复现了金属圆管对碰区鼓包现象,计算结果与实验吻合较好。可以预见,基于该方法能够初步分析影响金属圆管运动规律的诸多因素。     

爆炸加载下分层金属圆管膨胀破裂过程初探

设计了分层结构金属圆管，用硝基甲烷炸药爆炸产生的近似轴对称滑移爆轰波从内部加载分层管。采用高速分幅摄影技术对整个动态破裂发展过程进行拍摄，得到了清晰的照片。通过与等厚度、同材料单层管动态破裂过程的比较，取得了分层管破裂特性的初步认识。     

外部爆轰加载过程中金属圆管断裂实验研究

利用脉冲X光照相技术和高速摄影对HR 2抗氢钢圆管在外部爆轰加载下的压缩过程进行了动态诊断。X光照相观察到抗氢钢圆管在压缩过程中管壁内出现了宏观裂纹 ,回收破片具有显著的剪切断裂特征 ,金相分析观察到了绝热剪切带及微孔洞与微裂纹。对半圆管结构的压缩过程进行了高速摄影 ,直接观测到抗氢钢圆管内壁裂纹萌生与发展的动态过程。由两种动态诊断实验结果对比分析推断 ,抗氢钢圆管在外部爆轰加载压缩过程中 ,裂纹首先在圆管内壁生成。     

用二维欧拉方法模拟滑移爆轰作用下金属圆筒的收敛运动

本文用二维流体-弹塑性欧拉编码对金属圆筒在滑移爆轰作用下的收敛运动作了数值模拟。文中考虑了起爆端对金属圆筒运动的影响,给出了圆筒自由面不同时刻的运动姿态和自由面速度随时间变化规律。计算结果与实验结果符合较好。     

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

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

钢圆管在内爆加载下的层裂特性研究

利用完整回收技术,获得了20钢圆管在不同固体装药厚度、装药类型以及外壳条件的滑移内爆加载下形变和层裂的特征和规律,得到了钢圆管层裂的相对临界装药厚度等初步结果,分析了装药类型、装药厚度、外壳对钢圆管层裂的影响。     

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

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

对碰波作用下金属圆管材料参数对其运动特征影响的数值模拟

基于两点起爆实验,采用动力学有限元程序,在不考虑圆管破裂的前提下,详细分析了同一炸药爆轰下金属圆管各材料力学参数对其运动的影响规律。计算结果表明:弹性模量越大,泊松比越大,金属圆管的膨胀速度越大,但它们对圆管膨胀的位移和速度影响极小,在工程中可以忽略;屈服应力对金属圆管的运动有一定影响,屈服应力越大,对碰部位的鼓包范围越宽,鼓包峰值越小,但圆管各处的膨胀量均落后于屈服应力较小的情况;密度对圆管运动有较大影响,密度越大,受鼓包影响范围越宽,鼓包峰高越小。     

圆柱形爆轰波的二维数值模拟

基于带化学反应的二维Euler方程,对圆柱形爆轰波的直接起爆和传播过程进行了二维数值模拟研究,拟分析起爆条件和初始压强对圆柱形爆轰波形成和传播的影响。研究发现,圆柱形爆轰波起爆成功向外传播的过程中,新的三波结构的生成标志着爆轰波进入稳定传播阶段。在起爆能量足够的情况下,起爆半径(曲率)的大小决定着三波结构初始形成时的数目和传播半径,起爆压强对其基本不产生影响;起爆半径大(曲率小)时,三波结构初始形成时的传播半径大、数目多,圆柱形爆轰波进入稳定传播阶段的传播距离长;数值模拟中,初始压强的提高,有助于圆柱形爆轰在较短的传播距离内进入稳定传播阶段。     

爆轰等容终点的四组理论解

用Hamilton原理描述爆轰产物粒子运动并最终耗散为热的过程,证明爆轰过程反应区有四组理论解可以达到等容终点,即反应区的多变量解、等温过程解、等压过程解、等容过程解;修正了经典的ZND模型由于忽略输运效应而认为爆轰过程不可能达到等容终点的结论。本文结果与实验结果较为一致。     

计算钝感炸药爆轰过程的一种新反应率模型

根据钝感炸药爆轰过程中含有激化过程和慢反应,建立了一种新的反应率模型。与其他反应率模型相比,这种反应率模型可以应用于较粗的网格。在每厘米50个网格条件下,炸药驱动铝飞片和钽飞片的自由面计算结果与实验结果很接近。同时,应用此反应率模型计算了钝感炸药驱动LiF过程,在每厘米50个网格的条件下,炸药与LiF间速度的计算值同实验值接近,且误差随网格尺寸变小而变小。这些表明,此反应率模型能够在较粗的网格条件下,比较准确地描述钝感炸药驱动飞片过程,有利于在工程实际中应用。     

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

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

波后慢放能过程对钝感炸药爆轰过程的影响

采用唯象反应率模型与慢放能函数相结合方法,研究钝感炸药波后缓慢放能过程对爆轰过程的影响,计算结果与实验对比分析表明,慢放能反应率模型可更好地描述钝感炸药的爆轰物理过程.