点引爆圆板装药驱动飞板的实验及二维数值模拟

本文报导点引爆圆板装药驱动飞板的实验及二维数值模拟研究结果。实验测定了有机玻璃、铝、钢和铜四种飞板在爆轰产物驱动下的碰靶时间及碰靶波形,并将二维数值模拟结果与测量结果进行了比较,二者符合较好。文中还就数值模拟中存在的一些问题进行了简要分析、讨论。     

点引爆圆板装药驱动飞板的实验及二维数值模拟

本文报导点引爆圆板装药驱动飞板的实验及二维数值模拟研究结果。实验测定了有机玻璃、铝、钢和铜四种飞板在爆轰产物驱动下的碰靶时间及碰靶波形,并将二维数值模拟结果与测量结果进行了比较,二者符合较好。文中还就数值模拟中存在的一些问题进行了简要分析、讨论。     

准一维拉氏分析方法及其在凝聚炸药冲击波起爆研究中的应用

本文提出了一种准一维拉氏分析方法,用二维拉氏实验方法测得的不同拉氏位置上的压力和径向位移的变化曲线作为输入量,就能求解二维冲击波起爆过程中对称轴附近的流场。作为初步应用,文中列出了在两种加载条件下,压装TNT的二维冲击波起爆过程的分析结果。     

柱壳约束对散心冲击波空间分布的影响

结合实验和数值模拟，研究了散心冲击波在金属柱壳约束下沿有机玻璃内部的空间分布。进行了点起爆柱状炸药驱动飞片加载实验，采用Polyvinylidence Fluoride (PVDF)测试方法对有机玻璃内部的压力进行测试。实验结果显示:在冲击波传播过程中，在特定传播距离处，离中心轴越近，冲击波第一幅值压力越小，这是因为散心冲击驱动飞片成前凸形状，在飞片飞行过程中与有机玻璃碰撞面积越来越大，在远离对称轴部位冲击压力叠加累积效应更强引起的；但在随后的冲击波传播过程中，由于受到柱壳约束影响，离对称轴越近，冲击波幅值越小，这是由散心冲击波在约束柱壳边界反射与冲击波波阵面叠加的结果。通过对炸药网格大变形溢出柱壳翻转进行合理处理，对实验进行了数值模拟。数值模拟结果所得的冲击压力沿径向分布规律计算结果与实验结果定性相符。最后探讨了不同约束程度对这一规律的影响程度，结果表明，后续的冲击波幅值随着约束的增加而急剧增加。     

侧向稀疏波对非均质凝聚炸药冲击波起爆过程的影响

以压装TNT为例,开展二维拉氏实验,利用二维锰铜-康铜组合拉氏量计测量了轴对称位置和径向位置的压力时程曲线和径向位移曲线.同时利用二维拉格朗日分析计算程序RFLA对压装TNT的冲击波起爆过程进行了计算,分别得到了在2个不同径向位置上各个Lagrange位置上的质点速度、相对比容、比内能等时程曲线.研究结果展现了侧向稀疏...     

PELE贯穿薄靶后外壳破片径向速度计算方法与影响因素分析

基于横向效应增强型弹丸(PELE)侵彻金属薄靶板过程分析,将弹体前端在撞击作用下的变形过程分解为轴向一维压缩和径向自由膨胀两个变形阶段;依据冲击波理论,给出了弹体前端的冲击波压缩势能,由功能转化原理,给出了PELE前端外壳在靶后形成破片的最大径向飞散速度计算公式。计算结果在多种工况下均与文献的实验结果较为一致。计算结果表明:PELE靶后外壳破片的最大径向飞散速度与外壳和内芯材料的体积模量和泊松比有关,且随二者的增大而增大;PELE外壳破片的最大径向飞散速度是壳体和内芯在冲击波压缩作用下共同径向膨胀的结果,且外壳膨胀能在弹体整体膨胀能中所占比例较大,计算中应当同时考虑弹体外壳和内芯材料的横向膨胀效应对弹体破片径向飞散速度的影响。     

缸内激波对锥顶型燃烧室的活塞破坏机理

以二维数值模拟为基础,研究了锥顶型燃烧室内的冲击波发展的震荡过程,得到作用于活塞不同位置处的超压分布。模拟结果表明:由于燃烧室结构的独特性,导致冲击波能在特定区域进行汇聚,致使该区域超压明显高于其他区域。将该模拟结果与实际破坏失效的活塞进行对比,发现冲击波汇聚区域往往就是活塞被破坏的地方。数值模拟结果与实际破坏结果吻合很好。这为设计燃烧室形状以避免冲击波对活塞造成破坏提供了理论基础。     

预制孔靶板在爆炸冲击波载荷作用下的动态响应

对四边固定带预制孔靶板在爆炸冲击波载荷作用下的动态响应问题进行了实验研究。通过设计 爆炸冲击波载荷对2种材料不同孔数、不同孔径的靶板的冲击工况,获得了各靶板的中心点挠度。通过分析, 得到了预制孔靶板中心点挠度与预制孔数、孔径的相关规律,结合无孔平板在爆炸冲击波载荷下的挠度公式, 建立了广义的靶板中心点挠度公式,可用于计算不同孔靶板的中心点挠度。研究结果可作为平板材料在破片 冲击波联合毁伤研究的基础,为杀爆战斗部威力评估提供一些依据。     

不同当量强爆炸早期火球现象的数值模拟

采用一维球对称辐射流体力学方程组,研究了不同当量强爆炸早期火球阵面、冲击波的形成发展过程,对计算结果进行了相似律分析。结果表明,随爆炸当量增加,火球阵面、冲击波扩张速度加快,火球扩张过渡阶段开始结束位置及持续范围增加,火球中心温度下降变慢;在火球阵面以冲击波扩张传播以后,火球阵面参量满足立方根相似律。     

点起爆炸药驱动圆板击靶波形的数值模拟研究

用自行研制的二维流体动力学程序TDY2D对点起爆驱动圆板运动的3个模型进行数值模拟计算,得到的飞片击靶波形与实验波形均符合较好,具有较高的计算精度,验证了该程序的准确性和有效性。通过对爆轰驱动飞片运动过程数值模拟结果的分析,得到了装置参数对飞片击靶波形影响规律的初步认识。     

CE/SE方法在内弹道瞬态两相流中的应用

根据CE/SE方法思想,构造在变换坐标系下内弹道两相流动数值计算的 格式. 采用此格式研究了膛内瞬态二维两相流动点火过程. 数值计算结果反映实际 规律并和试验结果符合. 结果表明,对中心点火方式,在点传火过程初期,膛内压力在 径向和轴向上存在较大的压力梯度,当火药全部着火后,药室内压力在径向上逐渐趋于 一致. 这表明CE/SE方法能够完全模拟瞬态二维两相流动的情况.     

内埋炸药下超高韧性水泥基复合材料的抗爆性能

为研究超高韧性水泥基复合材料(ultra-high toughness cementitious composites, UHTCC)在内埋炸药爆炸下的抗爆性能和损伤破坏规律,对不同炸药埋深下的UHTCC和高强混凝土(high-strength concrete, HSC)进行了内埋炸药抗爆实验。得到了两种材料靶体的破坏状态,并利用接触爆炸的实验结果计算出了两种材料的抗爆性能参数。结果表明,在相同条件下,UHTCC抗爆性能优于高强混凝土。为了进一步探究UHTCC的抗压强度、抗拉强度以及拉伸韧性对靶体在内埋炸药下抗爆性能的影响,首先,采用改进的K&C模型对炸药埋深为40 mm的超高韧性水泥基复合材料靶体进行数值模拟,模拟结果与实验结果基本吻合,并根据数值模拟的结果得到了爆炸冲击波沿靶体径向衰减速度大于轴向衰减速度这一规律,验证了数值模型的有效性;然后,通过调整改进K&C模型中与抗压强度、抗拉强度以及拉伸韧性相关的参数,数值预测了不同抗压强度、抗拉强度以及拉伸韧性下UHTCC靶体的破坏状态,发现增强UHTCC的韧性可以有效防止靶体发生整体性破坏,增大UHTCC的抗拉强度可以减小靶体迎爆面的开坑直径,增大UHTCC的抗压强度对减小开坑直径效果不明显。     

激光驱动液滴迁移的机理研究

液滴热毛细迁移是微重力流体科学中的典型科学问题, 微重力液滴动力学研究不仅具有流体力学的理论意义, 而且具有重要的实际应用价值. 建立了二维轴对称激光驱动液滴迁移模型, 通过仿真计算研究微重力环境下激光驱动液滴迁移的过程, 研究了液滴直径、母液参数等对液滴迁移速度及行为的影响. 首先研究了母液和液滴对激光系数均较小, 液滴初始位置不同时液滴的迁移行为; 然后研究了母液对激光吸收系数较小, 液滴对激光吸收系数较大时, 不同液滴直径与母液宽度比条件下液滴的迁移行为. 仿真结果表明: 当母液和液滴对激光的吸收系数都很小时, 液滴迁移的方向主要受到液滴初始位置的影响; 当母液对激光的吸收系数较小, 液滴对激光的吸收系数较大时, 液滴会朝激光方向迁移, 液滴初始位置对迁移方向影响较小, 但液滴直径与母液宽度之比会影响液滴迁移行为. 将模拟结果与YGB理论对比, 仿真结果与理论结果趋势一致. 研究激光驱动液滴迁移的物理机制, 探索界面张力作用机理, 得到激光驱动液滴迁移的规律, 探索对液滴的驱动控制方法.      

海拔高度对长直坑道内爆炸冲击波传播的影响

为有效表征不同海拔坑道内爆炸冲击波的传播特征,利用非线性显式动力学有限元软件AUTODYN,研究了海拔高度对长直坑道内爆炸冲击波传播的影响规律,探讨了高海拔环境对坑道内冲击波传播的影响,基于量纲分析,建立了适用于不同海拔高度典型坑道内冲击波峰值超压的计算模型,并通过数值计算进行了验证。结果表明：随着海拔高度升高,坑道内爆炸冲击波波阵面传播速度与径向的冲击波参数偏差增大,平面波形成距离增加,冲击波峰值超压降低;在0～4 000 m范围内,海拔高度每升高1 000 m,冲击波冲量降低约0.91%。结合Sachs无量纲修正方法和量纲分析,推导出不同海拔高度冲击波峰值超压的理论分析模型,模型计算结果与数值计算结果的相对偏差不大于10%,能够为高海拔环境下坑道内爆炸冲击波的传播提供理论依据。     

长杆弹撞击陶瓷靶的一种数值模拟方法

陶瓷材料具有高强度和低密度等特点，抗弹性能优越，被广泛用于各类装甲中。长杆弹撞击陶瓷靶时会发生径向流动、质量显著侵蚀而无明显侵彻的界面击溃现象，是陶瓷抗侵彻性能研究中具有重要研究价值的特殊现象。利用有限元软件AUTODYN建立了长杆弹撞击陶瓷靶的二维轴对称计算模型，采用Lagrange和光滑粒子流体动力学（smooth particle hydrodynamics, SPH）算法，模拟了柱形钨合金长杆弹撞击带盖板的碳化硅陶瓷，通过改变长杆弹的撞击速度，得到了界面击溃、驻留转侵彻和直接侵彻3个不同现象。讨论了不同建模算法、边界条件以及材料参数对模拟结果的影响。通过网格收敛性验证和与实验结果进行拟合，综合验证了计算模型中算法、边界条件和参数设定的可靠性。结果表明，在建模中若同时使用SPH算法和Lagrange算法，需要考虑粒子和网格大小对于模拟结果的影响。针对长杆弹撞击陶瓷靶的界面击溃模拟，不建议对陶瓷材料采用SPH粒子建模。相关建模和参数选择方法对后续陶瓷抗侵彻/界面击溃的数值模拟具有重要的指导意义。     

改进WENO格式及其在爆轰计算中的应用

对已有的一种改进型WENO-M格式进行了理论上的修正,得到了新的WENO权函数。与原WENO-M格式相比,新得到的WENO-M-P格式减少了约9%的CPU计算时间,同时也保证了在一阶极值点处不降低精度（经典的WENO格式在一阶极值点处精度下降）,仍然保持5阶精度。为了验证修正后的格式,采用二步反应模型数值模拟了几组一维和二维爆轰：在一维爆轰模拟中,对比了WENO-M格式和WENO-M-P格式在一定温度下的最低起爆压强以及一定压强下的最低起爆温度;在二维旋转爆轰的模拟中,对比了WENO-M格式和WENO-M-P格式的计算效果。结果表明：在同样的起爆温度下,WENO-M 格式和 WENO-M-P 格式的最低起爆压强均高于 WENO格式,但WENO-M-P格式比WENO-M格式的最低起爆压强低;三种格式在不同压强下的最低起爆温度相同;WENO-M-P格式具有和WENO-M格式一样的计算效果,优势是节省了计算时间。     

磁驱动等熵压缩实验构形的磁流体力学计算模拟

在一维弹塑性反应流体力学拉格朗日编码SSS 的基础上扩充编制成功的一维磁流体力学计算编码SSS/MHD,可以计算模拟磁驱动等熵压缩（ICE） 实验中多介质、多空腔的不同厚度平行样品构形（台阶靶） 运动,并且把磁流体计算与模拟脉冲功率驱动源的集中参数外电路方程组直接耦合,不依赖于实验当时测量的驱动电流（电压） 数据,直接根据原始数据对实验样品的力学和电磁学运动过程和各种参数剖面进行数值模拟. 介绍了对于纯铝、钽和塑料粘结炸药等3 种性能不同材料的台阶靶样品ICE 实验的计算模拟,得到了与激光速度干涉测量的不同厚度样品后自由面（或光学窗口界面） 速度历史十分符合的结果,对于实验结果的分析和实验设计的改进具有重要意义. 旨在为磁驱动实验系统的校核、设计和实验构形改进提供有一定预测能力、功能比较完善的计算模拟研究手段.      

MAGNETO-HYDRODYNAMICS MODELING OF CONFIGURATIONS IN MAGNETICALLY DRIVEN ISENTROPIC COMPRESSION EXPERIMENTS

在一维弹塑性反应流体力学拉格朗日编码SSS 的基础上扩充编制成功的一维磁流体力学计算编码SSS/MHD,可以计算模拟磁驱动等熵压缩（ICE） 实验中多介质、多空腔的不同厚度平行样品构形（台阶靶） 运动,并且把磁流体计算与模拟脉冲功率驱动源的集中参数外电路方程组直接耦合,不依赖于实验当时测量的驱动电流（电压） 数据,直接根据原始数据对实验样品的力学和电磁学运动过程和各种参数剖面进行数值模拟. 介绍了对于纯铝、钽和塑料粘结炸药等3 种性能不同材料的台阶靶样品ICE 实验的计算模拟,得到了与激光速度干涉测量的不同厚度样品后自由面（或光学窗口界面） 速度历史十分符合的结果,对于实验结果的分析和实验设计的改进具有重要意义. 旨在为磁驱动实验系统的校核、设计和实验构形改进提供有一定预测能力、功能比较完善的计算模拟研究手段.     

爆炸与冲击问题的可视化研究

针对二维多物质流体弹塑性程序SMMIC(自行开发 )的结果 ,采用VisualC+ + 语言 ,自行开发设计了一个切实可行的二维爆炸场可视化工具软件 ,以此来研究二维爆炸与冲击数值模拟结果的可视化问题。对混凝土中的爆炸问题做了详细分析 ,对数值模拟结果进行了可视化处理 ,可视化的效果比较令人满意。从可视化的结果中不仅可以展现爆炸与冲击过程的整体概貌 ,而且可以清晰的显示爆炸与冲击过程的细微过程 (如爆炸场中冲击波遇障碍物后的反射、绕流及马赫反射等 )。     

钝感炸药爆轰波传播和驱动的数值模拟研究

为实现对非理想爆轰波传播和爆轰驱动问题的精确计算,研究了与贴体坐标下Level Set方法相结合的爆轰冲击波动力学(DSD)计算方法;根据Hamilton-Jacobi方程的Godunov差分格式,提出了非正交的贴体坐标系中Level Set函数方程的差分格式及其相应的数值方法;并将该DSD方法与二维流体动力学程序TD...