爆炸驱动动能杆定向抛撒散布试验

动能反导技术作为新型的反导技术，原理是在拦截弹和目标的相对速度方向上，抛射出大量动能杆条毁伤目标。其核心技术之一就是动能元素的抛撒，特别是定向抛撒。国外研究一般集中针对既定抛撒装置结构的杆条抛撒结果验证性报道，缺乏对杆条抛撒影响因素的详细研究。本文提出了采用周向弧形状药选择区域爆炸驱动中心杆条束定向侧向飞散的动能扦定向抛撒装置方案，设计不同结构参数的杆条定向抛撒装置的试验，研究抛撒结构参数对动能杆的驱动特性影响规律，为新概念反导战斗部设计提供技术帮助。     

点起爆炸药驱动平板飞片运动的数值模拟研究

 为研究点起爆发散爆轰波驱动飞片的运动特性,自行研制了二维流体动力学程序TDY2D,采用由实验所标定炸药的状态方程参数,对3种长径比相差较大的点起爆炸药驱动平板飞片运动的模型进行了数值模拟计算,并对结果进行了分析。计算所得的飞片运动历史和击靶波形均与实验符合较好,具有较高的计算精度,从而验证了该程序的准确性和有效性。     

燃料爆炸抛撒过程的实验研究

 通过高速录像对3种量级（5、20、50 kg）FAE装置的静爆实验进行了现场记录，通过对燃料爆炸抛撒形成云雾的成长过程的实时记录结果，把燃料的爆炸抛撒过程分为3个阶段：喷出阶段、过渡阶段和膨胀阶段，其依据是爆炸驱动力和气动阻力的相互作用而导致的云雾径向膨胀速度的变化；对于结构相似的FAE装置，爆炸抛撒各阶段结束时刻及最终的云雾直径与装填燃料量的立方根成正比，说明小型装置可以模拟大型FAE武器；此外对应于不同阶段结束时刻的云雾直径与装置直径的比值分别为5.5、30和45。     

水下爆炸冲击波数值仿真研究

由于在水下爆炸冲击波的数值仿真研究中,水的状态方程、人工黏性系数和网格尺寸对数值计算结果影响很大,采用常规TNT炸药的水下爆炸为例,以冲击波的峰值压力和比冲量为衡量指标,研究了这3个主要影响因素对数值仿真结果的影响。首先,通过采用常用的5种水的状态方程进行系列仿真,给出了各种状态方程的适用范围;其次,讨论了人工黏性系数对计算结果的影响,并给出了一次与二次人工黏性系数的建议取值范围;最后,通过对不同炸药当量及不同网格尺寸开展系列运算,从而得到不同炸药当量在满足工程计算精度要求下所对应的建议网格尺寸,并得到了不同炸药当量所对应的建议网格尺寸的表达式。     

燃料抛撒成雾及其燃烧爆炸的光滑离散颗粒流体动力学方法数值模拟研究

燃料在炸药爆炸驱动下形成燃料空气爆炸云团, 进而引燃爆炸, 对目标造成毁伤. 本文在前期提出的光滑离散颗粒流体动力学方法(SDPH)的基础上, 引入描述炸药由爆轰到膨胀整个过程的Jones-Wilkins-Lee状态方程及描述气体快速燃烧过程的EBU-Arrhenius燃烧模型, 建立了求解战斗部起爆、燃料抛撒和燃料二次引燃爆炸问题的新型SDPH方法. 设计了圆环形燃料颗粒在炸药爆炸驱动下运动抛撒的算例进行数值验证, 结果与理论相符; 对燃料空气炸药(FAE)云雾的形成和发展过程进行了数值模拟, 分析了云雾的形态, 并与实验结果进行对比, 符合较好, 同时分析了不同起爆方式对云雾团成型的影响; 最后, 在云雾团成型的基础上, 引入蒸发燃烧模型对FAE的燃烧爆炸过程进行了模拟研究. 结果表明, 本文建立的数学模型和计算方法可以较好的模拟燃料空气炸药抛撒成雾及云雾燃烧爆炸过程, 为该类武器装备的设计研究提供了较好的数值方法.     

FAE装置参数对燃料抛撒与爆炸威力影响的实验研究

 采用高速运动分析系统对FAE实验装置爆炸抛撒过程进行观测,描述了燃料抛撒过程的不同阶段,实验研究了比药量、长径比、壳体材质等装置参数对燃料抛撒与爆炸威力的影响。结果表明,燃料抛撒过程可分为射流形成与扩散运动阶段、燃料两向膨胀运动阶段和气液融合运动阶段。不同阶段对应不同的流体动力学特征,对云团形成的贡献不同。在实验装置总体优化条件下,适当增大比药量可提高云团覆盖面积与体积;在实验范围内长径比不是影响云团状态的显著因素,但长径比较大时可使燃料抛撒均匀性更好;采用钢质壳体时云雾抛撒状态明显优于铝质壳体。实验证明,采用碳钢壳体、比药量3%左右、长径比为3~5且装置参数良好匹配时,可获得理想的云团状态和高威力爆炸波毁伤效应。     

水中障碍物爆炸毁伤效应的二维数值模拟

 建立了水下障碍物的爆炸毁伤效应的物理模型和计算模型。建立了守恒方程组、混凝土的动态弹塑性本构关系和破坏准则,并确定了其在大应变、高应变率及高压条件下的的屈服准则。基于多物质流体的Euler算法,用面向对象的C++语言自行编制了NM-MMIC通用多物质二维流体弹塑性程序,对不同炸点位置、不同炸药量对轨条砦爆炸破坏过程进行数值模拟,并采用自行开发的后处理软件VISC2D实现了可视化,定量地反应了障碍物迎爆面所承受的压力变化过程,确定了轨条砦毁伤的最佳药量,并对毁伤程度进行了评价,得到了爆炸破障的毁伤规律。研究结果对设计和评价反水中障碍物的弹药具有一定的参考价值。     

深水爆炸冲击波作用下圆柱壳动态响应影响因素的数值模拟研究

 运用通用有限元程序ABAQUS,对圆柱壳在深水爆炸冲击波作用下的动态响应进行了数值模拟,研究了圆柱壳所处深度、爆心方位和预应力对于圆柱壳动态响应的影响。研究结果表明：在同样的水下爆炸冲击波作用下,随着深度的增大,圆柱壳的毁伤逐渐加重,而且响应过程中的变形逐渐增大,低频响应阶段的速度曲线趋于平坦;爆距不变时,不同的爆心方位在不同深度上的毁伤效果存在一定差异;预应力的存在使圆柱壳的毁伤加重,并且预应力对圆柱壳动态响应的影响与爆心方位有关。研究结果对于潜体结构的生命力评估具有一定参考价值。     

内部爆炸载荷作用下混凝土动力响应的数值模拟

 目前,采用实验方法对爆炸载荷破坏效应进行全面、准确描述仍是一个难题。运用ANSYS/LS-DYNA非线性显式动力学有限元程序,采用流固耦合算法,对混凝土在内部爆炸载荷作用下的空腔形成和发展规律等动力响应问题进行了数值模拟。通过实际算例,比较研究了混凝土敷设钢板与否对成坑大小、损伤区、压力云图以及应力时程曲线等的影响,进一步研究了有钢板加固时混凝土在爆炸荷载作用下的破坏机理。研究表明：采用流固耦合算法可以避免由于物质的大变形引起的单元畸变和网格消蚀等问题,提高了计算效率;很好地模拟了混凝土内部爆炸时,中心爆炸空腔压缩面和自由面的拉伸层裂破坏过程,且与理论计算值吻合较好;通过表面敷设钢板,可以减小混凝土的拉应力峰值,抑制混凝土中裂纹的发展,从而提高其防震塌能力以及结构的完整性。     

爆炸冲击波对玻璃幕墙破坏作用的多物质ALE有限元模拟

应用LS-DYNA显式程序,基于多物质的ALE有限元法,利用罚函数实现了爆炸冲击波与玻璃幕墙结构间的耦合作用,重现了玻璃幕墙防爆炸实验中爆炸荷载的产生及其与玻璃幕墙相互作用的三维动态过程。研究了爆炸冲击波作用下幕墙结构的动态响应行为,分析了幕墙内外层玻璃的破坏情况。研究表明,数值模拟结果与实验结果一致,这为玻璃幕墙的抗爆设计及改进提供了重要的参考依据。     

战斗部爆轰产物温度效应数值模拟

 利用BKW代码、JWL方程和LS-DYNA程序,研究了无约束、弱约束和强约束等三种战斗部约束条件下炸药爆轰产物密度随时间的变化关系,再结合爆轰产物HOM状态方程获得爆轰产物温度随时间的衰变函数。给出了战斗部爆轰产物温度的理论计算方法。研究表明,爆轰区内爆轰产物温度及其衰减速率并不均匀,战斗部炸药的约束条件对爆轰产物密度和温度有重要影响,强约束可以延缓爆轰产物密度和温度的衰减,相同情况下仅战斗部壳体厚度加倍并不会导致爆轰产物密度和温度的衰减情况同比例变化。     

面层混凝土对跑道内爆毁伤影响的仿真研究

 采用考虑混凝土拉伸破坏的TCK失效模型,运用AUTODYN软件对跑道内爆毁伤进行了数值模拟。研究了炸点置于面层混凝土底部时,不同装药量和不同面层混凝土厚度对爆破漏斗坑参数的影响规律。结果表明：数值模拟结果较好地反映了混凝土因拉伸损伤引起的破裂和抛掷漏斗现象。基于数值模拟得到的漏斗坑参数值,建立了漏斗坑半径、空腔半径与面层厚度的无量纲表达式,计算得到的漏斗坑参数值得到了实验验证。     

强爆炸早期火球现象的一维数值研究

采用一维球对称辐射流体力学方程组,构建直接求解辐射输运方程的一种数值方法,得到强爆炸早期火球发展过程中火球阵面、弹壳冲击波的形成发展过程.计算结果与Brode结果进行对比,二者一致,证明数值方法的合理性.数值模拟爆炸源区不同当量、不同密度以及不同爆炸高度等条件下强爆炸早期火球发展过程中火球阵面、弹壳冲击波的形成发展过程,对计算结果进行比较分析.     

高速杆式弹丸三维数值模拟

 应用LS-DYNA动力学分析软件对变球缺型罩装药结构下高速杆式弹丸的形成过程进行数值模拟研究,揭示了其成型机理,并研究了中心点起爆、四点同时起爆和环形起爆等三种起爆方式对高速杆式侵彻体的成型和性能参数的影响,得到了弹体的速度分布。数值计算结果显示：变球缺型罩装药结构能形成速度高、长径比大、质量分布合理的高速杆式侵彻体,其成型机理兼有聚能射流和翻转型爆炸成型弹丸的特点,即上半部罩微元直接向轴线压合,下半部罩微元向轴线压合时逐渐向后翻转变形。不同起爆方式下所形成的高速杆式弹丸的形状和性能参数有所不同,其中环形起爆所形成的弹体性能参数最优。实验结果所反映的起爆方式的影响规律与数值计算结果是一致的。     

不同方式地下爆炸地震耦合效应的数值模拟

通过一种Euler与Lagrange相结合的有限差分法,开展地下爆炸地震耦合效应的数值模拟.比较花岗岩介质中填实和空腔两种爆炸方式在地震耦合效应方面的差异,进一步分析初始空腔大小在填实以及空腔爆炸地震效应中的作用以及空腔爆炸解耦因子随频率的变化特性.计算结果表明,花岗岩介质中完全空腔解耦因子能达到50~60倍.     

爆炸荷载下不同壁厚圆柱壳动力学行为的研究

对爆炸荷载下圆柱壳的动力学行为进行了实验研究及数值模拟。将外径均为100mm的3种壁厚的Q235钢质圆柱壳置于由TNT药柱产生的爆炸场中进行冲击实验,系统分析了在不同装药高度及壳壁厚度参数条件下圆柱壳的冲击变形模式,即迎爆面局部凹陷变形模式(ModeⅠ)、迎爆面局部凹陷与壳整体弯曲变形耦合模式(ModeⅡ)、整体变形失效模式(ModeⅢ)及局部穿透与整体变形失效耦合模式(ModeⅣ)。采用LS-DYNA有限元程序及Lagrangian-Eulerian流固耦合算法,对圆柱壳的非线性动力响应过程进行了数值模拟,分析了圆柱壳的变形历程及最终残余变形的情况,计算结果与实验现象吻合较好。研究结论可为圆柱壳结构爆炸破坏分级及抗爆技术设计提供科学依据。