2003年 23卷 第6期
2003, 23(6): 481-487.
doi: 10.11883/1001-1455(2003)06-0481-7
PDF 
本文被引
摘要:
研究了一般物态方程形式下一维平面和柱、球面对称CJ爆轰波后产物的等熵流动,给出Riemann不变量的表达式、Euler和Lagrange形式的自相似流动方程组,得到一个高精度近似解析解,在实验测量合理的范围内本文理论与实验基本符合。为了具体应用上述理论,给出四种普遍公认的物态方程(律、JWL、HOM和W.C.Davis近来提出的形式)的等熵线和声速曲线,以及较一般情形中产物或惰性介质沿等熵线的比内能表达式。四种等熵线差别不大,表明爆轰数值模拟中不但要应用产物的等熵线,还应注意产物能量状态方程是否适当。本文的理论计算方法,可作为产物一维流动数值模拟正确性的有效校核手段。
研究了一般物态方程形式下一维平面和柱、球面对称CJ爆轰波后产物的等熵流动,给出Riemann不变量的表达式、Euler和Lagrange形式的自相似流动方程组,得到一个高精度近似解析解,在实验测量合理的范围内本文理论与实验基本符合。为了具体应用上述理论,给出四种普遍公认的物态方程(律、JWL、HOM和W.C.Davis近来提出的形式)的等熵线和声速曲线,以及较一般情形中产物或惰性介质沿等熵线的比内能表达式。四种等熵线差别不大,表明爆轰数值模拟中不但要应用产物的等熵线,还应注意产物能量状态方程是否适当。本文的理论计算方法,可作为产物一维流动数值模拟正确性的有效校核手段。
2003, 23(6): 488-492.
doi: 10.11883/1001-1455(2003)06-0488-5
摘要:
基于带化学反应的Navier Stokes方程和有关的热力学和反应动力学数据,利用改进的VLS格式,对甲烷 空气混合物中激波与火焰的相互作用进行了数值模拟。根据计算结果,讨论了激波掠过火焰时的变形、分叉和发展,以及激波作用下火焰的失稳、变形、破碎和相应的带旋涡的流场。
基于带化学反应的Navier Stokes方程和有关的热力学和反应动力学数据,利用改进的VLS格式,对甲烷 空气混合物中激波与火焰的相互作用进行了数值模拟。根据计算结果,讨论了激波掠过火焰时的变形、分叉和发展,以及激波作用下火焰的失稳、变形、破碎和相应的带旋涡的流场。
2003, 23(6): 493-500.
doi: 10.11883/1001-1455(2003)06-0493-8
摘要:
针对大变形流体动力学数值计算中经常需要应用的网格重构与物理量重映技术,提出了一种逻辑简单的质点积分守恒重映方法。将旧网格细分为众多有体积的质点,并将旧网格的物理量分配到各个质点,新网格各守恒量的积分直接由落在新网格内的所有质点的物理量累加。建立了收敛速度极快的计算格式,采用的控制体很好地解决了速度的重映计算问题。分析了此守恒重映方法的收敛性与守恒性,研究了积分控制体对速度计算的影响。
针对大变形流体动力学数值计算中经常需要应用的网格重构与物理量重映技术,提出了一种逻辑简单的质点积分守恒重映方法。将旧网格细分为众多有体积的质点,并将旧网格的物理量分配到各个质点,新网格各守恒量的积分直接由落在新网格内的所有质点的物理量累加。建立了收敛速度极快的计算格式,采用的控制体很好地解决了速度的重映计算问题。分析了此守恒重映方法的收敛性与守恒性,研究了积分控制体对速度计算的影响。
2003, 23(6): 501-508.
doi: 10.11883/1001-1455(2003)06-0501-8
摘要:
就大气呼吸模式的推进机理作了较为系统的分析和研究,导出了考虑反射面几何约束的冲量及冲量耦合系数的近似解析表达式。在数值模拟方面,采用前期点爆炸自模拟解和后期高分辨率PPM格式相结合的新方法,计算了固定的抛物面反射镜聚焦入射平行激光击穿空气形成的高温等离子体流场及其对反射面产生的推动作用。结果表明,推进效应(飞行器获得的推力、总冲量、冲量耦合系数等)与反射面的几何参数和入射激光强度密切相关。83J的激光单脉冲能量算例得到的冲量耦合系数为246Ns/J,略高于同种工况下W.O.Schall等实验得到的耦合系数。
就大气呼吸模式的推进机理作了较为系统的分析和研究,导出了考虑反射面几何约束的冲量及冲量耦合系数的近似解析表达式。在数值模拟方面,采用前期点爆炸自模拟解和后期高分辨率PPM格式相结合的新方法,计算了固定的抛物面反射镜聚焦入射平行激光击穿空气形成的高温等离子体流场及其对反射面产生的推动作用。结果表明,推进效应(飞行器获得的推力、总冲量、冲量耦合系数等)与反射面的几何参数和入射激光强度密切相关。83J的激光单脉冲能量算例得到的冲量耦合系数为246Ns/J,略高于同种工况下W.O.Schall等实验得到的耦合系数。
2003, 23(6): 509-515.
doi: 10.11883/1001-1455(2003)06-0509-7
摘要:
为了提高军盔等防护装备的防弹性能,采用非线性动态显式有限元分析软件ANSYS/LS DYNA5.6.1模拟了几种复合材料防弹靶板的抗弹道侵彻过程,弹头采用标准1.1g模拟弹片。结合v50弹击实验分析了靶板在此类弹击下的破坏机理和吸能方式,对靶板的选材、铺层顺序、层数和制造工艺等提出合理化建议,主要是沿靶板厚向应采用非均匀的三段式结构和工艺。
为了提高军盔等防护装备的防弹性能,采用非线性动态显式有限元分析软件ANSYS/LS DYNA5.6.1模拟了几种复合材料防弹靶板的抗弹道侵彻过程,弹头采用标准1.1g模拟弹片。结合v50弹击实验分析了靶板在此类弹击下的破坏机理和吸能方式,对靶板的选材、铺层顺序、层数和制造工艺等提出合理化建议,主要是沿靶板厚向应采用非均匀的三段式结构和工艺。
2003, 23(6): 516-522.
doi: 10.11883/1001-1455(2003)06-0516-7
摘要:
通过实验和数值模拟对泡沫铝中冲击波传播特性进行了研究,结果表明:冲击波在泡沫铝中传播时显示明显的衰减特性;与此同时波头升时逐渐增加。这种衰减耗散特性主要来源于泡沫铝本身的本构粘性效应,而追赶卸载效应又会进一步促进冲击波的衰减。这为泡沫铝作为新型抗冲击缓冲材料提供设计基础。
通过实验和数值模拟对泡沫铝中冲击波传播特性进行了研究,结果表明:冲击波在泡沫铝中传播时显示明显的衰减特性;与此同时波头升时逐渐增加。这种衰减耗散特性主要来源于泡沫铝本身的本构粘性效应,而追赶卸载效应又会进一步促进冲击波的衰减。这为泡沫铝作为新型抗冲击缓冲材料提供设计基础。
2003, 23(6): 523-528.
doi: 10.11883/1001-1455(2003)06-0523-6
摘要:
用液体炸药滑移内爆加载方式,研究了不同炸药厚度和外壳条件下20钢圆管的内聚变形和破坏特征。实验中,钢管发生了较大径向应变,应变值大致随炸药厚度的增加而线性增大。在炸药厚度不低于3mm时钢管变形均匀,轴对称较好;药厚小于3mm时其变形不再是基本轴对称收缩,管子表面出现扭转褶皱。各种装药厚度下,均未观察到预期的层裂破坏。对未层裂的原因进行了初步分析。
用液体炸药滑移内爆加载方式,研究了不同炸药厚度和外壳条件下20钢圆管的内聚变形和破坏特征。实验中,钢管发生了较大径向应变,应变值大致随炸药厚度的增加而线性增大。在炸药厚度不低于3mm时钢管变形均匀,轴对称较好;药厚小于3mm时其变形不再是基本轴对称收缩,管子表面出现扭转褶皱。各种装药厚度下,均未观察到预期的层裂破坏。对未层裂的原因进行了初步分析。
2003, 23(6): 529-533.
doi: 10.11883/1001-1455(2003)06-0529-5
摘要:
采用前照明的高速分幅照相对不同壁厚的45钢柱壳在爆轰加载下的膨胀破裂过程进行了研究,实验成功拍摄到45钢柱壳前表面裂纹生成、扩展及产物泄漏的过程图像。研究表明:膨胀应变率在104s-1附近时,随着应变率的增加,45钢柱壳的各特征应变均有所增加,断裂模式由拉剪混合向单一剪切转变;45钢柱壳的膨胀破裂有一个过程,从柱壳的外壁面可以观察到裂纹生成至产物开始泄漏经历时间历程的长短与柱壳壁厚相关,柱壳壁越厚,从外壁出现裂纹到产物开始泄漏的时间间隔越长;由回收破片厚度测量推算出的膨胀断裂应变与动态实验过程中产物开始泄漏时的断裂应变基本一致。
采用前照明的高速分幅照相对不同壁厚的45钢柱壳在爆轰加载下的膨胀破裂过程进行了研究,实验成功拍摄到45钢柱壳前表面裂纹生成、扩展及产物泄漏的过程图像。研究表明:膨胀应变率在104s-1附近时,随着应变率的增加,45钢柱壳的各特征应变均有所增加,断裂模式由拉剪混合向单一剪切转变;45钢柱壳的膨胀破裂有一个过程,从柱壳的外壁面可以观察到裂纹生成至产物开始泄漏经历时间历程的长短与柱壳壁厚相关,柱壳壁越厚,从外壁出现裂纹到产物开始泄漏的时间间隔越长;由回收破片厚度测量推算出的膨胀断裂应变与动态实验过程中产物开始泄漏时的断裂应变基本一致。
2003, 23(6): 534-538.
doi: 10.11883/1001-1455(2003)06-0534-5
摘要:
为了研究爆轰波在拐角中的传播特性以及拐角爆轰低压流场现象,用LSDYNA 3D程序对三种常见的拐角装药的爆轰波传播特性进行了数值模拟,清晰显示了爆轰波因拐角处爆轰波传播面积的变化而产生的衰减 增长过程,讨论了侧向三通中爆速的变化情况,并得到了拐角处低压爆轰区的尺寸。
为了研究爆轰波在拐角中的传播特性以及拐角爆轰低压流场现象,用LSDYNA 3D程序对三种常见的拐角装药的爆轰波传播特性进行了数值模拟,清晰显示了爆轰波因拐角处爆轰波传播面积的变化而产生的衰减 增长过程,讨论了侧向三通中爆速的变化情况,并得到了拐角处低压爆轰区的尺寸。
2003, 23(6): 539-544.
doi: 10.11883/1001-1455(2003)06-0539-6
摘要:
用二级轻气炮加载技术将W89Mo9Ni1Fe1飞片加速到约5.0km/s,撞击封装有接近常态的甲烷气体的LY12铝靶。利用六通道高温计和示波器记录到冲击压缩的甲烷气体的光辐射历史曲线,并获得甲烷气体的冲击波速度与温度。利用一维Euler方程与化学反应方程的全流场耦合,采用甲烷体系的13组元40反应步的化学反应模型,并用3阶WENO数值格式对空间进行离散。为避免刚性过强,化学反应源项用点隐方法处理。由此得到了甲烷的冲击反应流场和波后热力学参数。通过比较,数值结果与实验结果符合较好。最后,解释了甲烷用作电探针保护气体的物理机制。
用二级轻气炮加载技术将W89Mo9Ni1Fe1飞片加速到约5.0km/s,撞击封装有接近常态的甲烷气体的LY12铝靶。利用六通道高温计和示波器记录到冲击压缩的甲烷气体的光辐射历史曲线,并获得甲烷气体的冲击波速度与温度。利用一维Euler方程与化学反应方程的全流场耦合,采用甲烷体系的13组元40反应步的化学反应模型,并用3阶WENO数值格式对空间进行离散。为避免刚性过强,化学反应源项用点隐方法处理。由此得到了甲烷的冲击反应流场和波后热力学参数。通过比较,数值结果与实验结果符合较好。最后,解释了甲烷用作电探针保护气体的物理机制。
2003, 23(6): 545-550.
doi: 10.11883/1001-1455(2003)06-0545-6
摘要:
详细描述了依据损伤原理建立的连续损伤模型,并对该模型进行了改进。在LS DYNA程序用户自定义材料模型中加入改进的连续损伤模型,并对弹体侵彻钢筋混凝土的穿孔过程进行了数值模拟,其结果与实验结果相吻合,模型可以用于钢筋混凝土的动态破坏预报。
详细描述了依据损伤原理建立的连续损伤模型,并对该模型进行了改进。在LS DYNA程序用户自定义材料模型中加入改进的连续损伤模型,并对弹体侵彻钢筋混凝土的穿孔过程进行了数值模拟,其结果与实验结果相吻合,模型可以用于钢筋混凝土的动态破坏预报。
2003, 23(6): 551-560.
doi: 10.11883/1001-1455(2003)06-0551-10
摘要:
提出了一个钢坯爆炸焊接的计算模型。模型的数学表达式表明了振型函数、广义坐标及其频率与钢坯力学性能、土基力学性能和钢坯尺寸的具体关系。通过挠度的解析表达式,可获得这种爆炸加工的动态响应特性,为工程实践中优化参数和动态分析提供了有效工具。
提出了一个钢坯爆炸焊接的计算模型。模型的数学表达式表明了振型函数、广义坐标及其频率与钢坯力学性能、土基力学性能和钢坯尺寸的具体关系。通过挠度的解析表达式,可获得这种爆炸加工的动态响应特性,为工程实践中优化参数和动态分析提供了有效工具。
2003, 23(6): 561-564.
doi: 10.11883/1001-1455(2003)06-0561-4
摘要:
采用多点起爆方式,设计了带尾翼翻转型爆炸成形弹丸(EFP)试验装置。用X光和SVR数字相机拍摄了EFP的外形和速度;用多层纸靶测试了EFP在不同飞行距离时的飞行姿态;进行了EFP穿靶能力的检验;并利用泡沫和锯末进行了EFP的软回收。由试验结果知,EFP速度为1.56～1.72km/s,长径比最大达到了3.69,由药型罩转变为EFP的质量转换率达到了98%,EFP具有较好的尾翼结构和气动稳定外形。该EFP能穿透厚度为50mm的厚钢靶或厚度为6mm、间距各为1m的5层薄钢靶。
采用多点起爆方式,设计了带尾翼翻转型爆炸成形弹丸(EFP)试验装置。用X光和SVR数字相机拍摄了EFP的外形和速度;用多层纸靶测试了EFP在不同飞行距离时的飞行姿态;进行了EFP穿靶能力的检验;并利用泡沫和锯末进行了EFP的软回收。由试验结果知,EFP速度为1.56～1.72km/s,长径比最大达到了3.69,由药型罩转变为EFP的质量转换率达到了98%,EFP具有较好的尾翼结构和气动稳定外形。该EFP能穿透厚度为50mm的厚钢靶或厚度为6mm、间距各为1m的5层薄钢靶。
2003, 23(6): 565-568.
doi: 10.11883/1001-1455(2003)06-0565-4
摘要:
用动态光弹性试验模拟研究了地质断层对应力波的影响及相互作用,比较分析了冲击荷载作用下水平和30方向的裂缝的动态响应,得到指定点处的条纹级数 时间关系曲线图和应力 时间关系曲线图,研究了结构物顶部震塌的机理,试验结论可为防护结构的改进提供参考。
用动态光弹性试验模拟研究了地质断层对应力波的影响及相互作用,比较分析了冲击荷载作用下水平和30方向的裂缝的动态响应,得到指定点处的条纹级数 时间关系曲线图和应力 时间关系曲线图,研究了结构物顶部震塌的机理,试验结论可为防护结构的改进提供参考。
2003, 23(6): 569-572.
doi: 10.11883/1001-1455(2003)06-0569-4
摘要:
在进行炸药驱动铜板实验中,用激光速度干涉仪(VISAR)记录铜板自由面速度,观察到了铜板层裂现象,并用非线性有限元方法对爆炸载荷下铜板层裂过程进行数值模拟分析,得到了与实验结果吻合的计算结果。
在进行炸药驱动铜板实验中,用激光速度干涉仪(VISAR)记录铜板自由面速度,观察到了铜板层裂现象,并用非线性有限元方法对爆炸载荷下铜板层裂过程进行数值模拟分析,得到了与实验结果吻合的计算结果。
2003, 23(6): 573-576.
doi: 10.11883/1001-1455(2003)06-0573-4
摘要:
介绍了内爆炸试验的发展及现状,通过实验研究建立了可以对高能炸药在坑道内爆炸的空气冲击波到时进行预计的公式。该公式适用于爆点在固定横截面的直通道口外、口内及口部处爆炸的情况。利用该公式可以求出空气冲击波在坑道中传播速度的变化。
介绍了内爆炸试验的发展及现状,通过实验研究建立了可以对高能炸药在坑道内爆炸的空气冲击波到时进行预计的公式。该公式适用于爆点在固定横截面的直通道口外、口内及口部处爆炸的情况。利用该公式可以求出空气冲击波在坑道中传播速度的变化。
