2012年 32卷 第1期
2012, 32(1): 1-7.
doi: 10.11883/1001-1455(2012)01-0001-07
PDF 
本文被引
摘要:
首先基于系统能量守恒条件,提出了一种计算蒸气等离子体压力的一维耦合计算模型。模型中不 仅考虑了蒸气等离子体界面压力与质点速度的非线性效应,同时也考虑了界面烧蚀所致的运动速度,将蒸气 等离子体的膨胀与约束介质的变形耦合。在耦合模型的基础上,采用显式差分计算程序与显式有限元计算程 序LS-DYNA互相迭代求解的方法,对不同激光功率密度分布下的蒸气等离子体压力进行了计算。结果表 明,计算结果与实验测量结果具有很好的一致性,证明了计算模型的合理性。
首先基于系统能量守恒条件,提出了一种计算蒸气等离子体压力的一维耦合计算模型。模型中不 仅考虑了蒸气等离子体界面压力与质点速度的非线性效应,同时也考虑了界面烧蚀所致的运动速度,将蒸气 等离子体的膨胀与约束介质的变形耦合。在耦合模型的基础上,采用显式差分计算程序与显式有限元计算程 序LS-DYNA互相迭代求解的方法,对不同激光功率密度分布下的蒸气等离子体压力进行了计算。结果表 明,计算结果与实验测量结果具有很好的一致性,证明了计算模型的合理性。
2012, 32(1): 8-14.
doi: 10.11883/1001-1455(2012)01-0008-07
摘要:
对某高强度钢制成的光滑圆棒和缺口圆棒进行了系列准静态拉伸实验,采用ABAQUS对每个试 件进行了数值模拟,得到了该材料的真实应力应变曲线,拟合出了J-C本构模型和失效模型的部分材料常数。 最后,对该高强度钢制成的平板进行了撞击实验,并用得到的J-C模型对平板撞击实验进行了数值模拟,计算 结果与实验结果吻合很好,证明利用数值模拟并修正应力三轴度的方法是可行的。
对某高强度钢制成的光滑圆棒和缺口圆棒进行了系列准静态拉伸实验,采用ABAQUS对每个试 件进行了数值模拟,得到了该材料的真实应力应变曲线,拟合出了J-C本构模型和失效模型的部分材料常数。 最后,对该高强度钢制成的平板进行了撞击实验,并用得到的J-C模型对平板撞击实验进行了数值模拟,计算 结果与实验结果吻合很好,证明利用数值模拟并修正应力三轴度的方法是可行的。
2012, 32(1): 15-22.
doi: 10.11883/1001-1455(2012)01-0015-08
摘要:
基于H.Meng的动摩擦实验结果,建立了一个端面动摩擦模型,以描述端面动摩擦因数随SHPB 塑料试样与钢压杆接触界面间最大径向相对滑动速度的关系。在SHPB实验的数值模拟中,分别考虑了该 端面动摩擦模型和以往被广泛采用的常摩擦因数模型,结果发现,选取不同的端面摩擦模型不影响通过塑料 SHPB实验的数值模拟获得的转折应变率(量级为102s-1),然而当应变率超出转折应变率以上时,不同的端 面摩擦模型会对塑料SHPB实验的数值模拟结果产生显著影响,使用提出的端面动摩擦模型可得到更准确 的塑料SHPB实验的数值模拟结果。
基于H.Meng的动摩擦实验结果,建立了一个端面动摩擦模型,以描述端面动摩擦因数随SHPB 塑料试样与钢压杆接触界面间最大径向相对滑动速度的关系。在SHPB实验的数值模拟中,分别考虑了该 端面动摩擦模型和以往被广泛采用的常摩擦因数模型,结果发现,选取不同的端面摩擦模型不影响通过塑料 SHPB实验的数值模拟获得的转折应变率(量级为102s-1),然而当应变率超出转折应变率以上时,不同的端 面摩擦模型会对塑料SHPB实验的数值模拟结果产生显著影响,使用提出的端面动摩擦模型可得到更准确 的塑料SHPB实验的数值模拟结果。
2012, 32(1): 23-28.
doi: 10.11883/1001-1455(2012)01-0023-06
摘要:
采用显式动力有限元方法研究了具有不同胞元结构的六角形铝蜂窝在冲击荷载下的力学性能,讨 论了铝蜂窝的变形模式和动态承载力以及影响因素。通过改变胞壁夹角得到5种不同的胞元结构,计算采用 了3种冲击速度。结果表明,在准静态变形模式下,胞元的几何因素对铝蜂窝的承载力起主导作用;一旦蜂窝 的变形呈现动态模式后,惯性效应显著,对蜂窝承载力起决定作用,胞元几何因素的影响不再明显;在过渡模 式下,惯性效应与几何因素共同主导蜂窝的动态承载力,并且冲击速度越高,惯性效应的影响越大。
采用显式动力有限元方法研究了具有不同胞元结构的六角形铝蜂窝在冲击荷载下的力学性能,讨 论了铝蜂窝的变形模式和动态承载力以及影响因素。通过改变胞壁夹角得到5种不同的胞元结构,计算采用 了3种冲击速度。结果表明,在准静态变形模式下,胞元的几何因素对铝蜂窝的承载力起主导作用;一旦蜂窝 的变形呈现动态模式后,惯性效应显著,对蜂窝承载力起决定作用,胞元几何因素的影响不再明显;在过渡模 式下,惯性效应与几何因素共同主导蜂窝的动态承载力,并且冲击速度越高,惯性效应的影响越大。
2012, 32(1): 29-33.
doi: 10.11883/1001-1455(2012)01-0029-05
摘要:
为了研究双层药型罩形成的串联EFP,设计了4种不同材料组合的双层药型罩,利用X光获得了 串联EFP形成过程的照片,并开展了相应的数值模拟。发现具有合适结构和药型罩材料的双层药型罩可以 形成串联聚能侵彻体,药型罩材料组合和起爆位置对串联聚能侵彻体的形成具有较大的影响。
为了研究双层药型罩形成的串联EFP,设计了4种不同材料组合的双层药型罩,利用X光获得了 串联EFP形成过程的照片,并开展了相应的数值模拟。发现具有合适结构和药型罩材料的双层药型罩可以 形成串联聚能侵彻体,药型罩材料组合和起爆位置对串联聚能侵彻体的形成具有较大的影响。
2012, 32(1): 34-42.
doi: 10.11883/1001-1455(2012)01-0034-09
摘要:
为了讨论率敏感材料软钢钢梁受矩形脉冲载荷作用下的动力响应问题,通过直接离散有限变形弹 塑性连续体最小加速度原理中的Lee泛函得到基本控制方程,并将包含应变率的Cowper-Symonds方程嵌入 应力-应变关系方程,使该计算模型计及材料的应变率效应,因而能够准确描述钢梁受爆炸和冲击载荷作用下 的动力响应问题。该计算模型的有效性通过与通用有限元程序ABAQUS比较而得到了验证,并进一步与已 有的刚塑性解做了对比。利用该计算模型进行数值计算,分析了均布和集中脉冲载荷作用下钢梁动力响应的 应变率效应,结果发现对于钢梁存在新的异常行为响应模式,应变率导致异常区域偏移和扩大,已有的刚塑性 解在一定载荷强度范围内不能准确预报钢梁的实际位移。
为了讨论率敏感材料软钢钢梁受矩形脉冲载荷作用下的动力响应问题,通过直接离散有限变形弹 塑性连续体最小加速度原理中的Lee泛函得到基本控制方程,并将包含应变率的Cowper-Symonds方程嵌入 应力-应变关系方程,使该计算模型计及材料的应变率效应,因而能够准确描述钢梁受爆炸和冲击载荷作用下 的动力响应问题。该计算模型的有效性通过与通用有限元程序ABAQUS比较而得到了验证,并进一步与已 有的刚塑性解做了对比。利用该计算模型进行数值计算,分析了均布和集中脉冲载荷作用下钢梁动力响应的 应变率效应,结果发现对于钢梁存在新的异常行为响应模式,应变率导致异常区域偏移和扩大,已有的刚塑性 解在一定载荷强度范围内不能准确预报钢梁的实际位移。
2012, 32(1): 43-46.
doi: 10.11883/1001-1455(2012)01-0043-04
摘要:
为了探索泡沫材料在空间防护上的应用,结合实验和数值模拟,绘制了弹速0.5~7.0km/s范围 内填充式泡沫铝防护结构的弹道极限曲线,并与Whipple结构进行了比较。结果表明:填充式泡沫铝防护结 构的防护性能优于同等面密度、总厚度的Whipple结构;孔隙率对填充式泡沫铝防护结构防护性能影响较小。
为了探索泡沫材料在空间防护上的应用,结合实验和数值模拟,绘制了弹速0.5~7.0km/s范围 内填充式泡沫铝防护结构的弹道极限曲线,并与Whipple结构进行了比较。结果表明:填充式泡沫铝防护结 构的防护性能优于同等面密度、总厚度的Whipple结构;孔隙率对填充式泡沫铝防护结构防护性能影响较小。
2012, 32(1): 47-50.
doi: 10.11883/1001-1455(2012)01-0047-04
摘要:
基于矿井生产环境,构建了相应的物理模型、数学模型,进行了数值模拟。模拟过程中,考虑了化 学反应,将爆炸热力学过程和基元反应耦合,采用TVD格式进行数值离散。结果揭示巷道中设置1~3辆矿 车时,瓦斯爆炸冲击波及其火焰波的结构和参数变化特征。数值模拟与实验结果对比分析表明,单一矿车和 多辆矿车在激励效应存在明显不同,矿车数量增加,激励效应更加显著。
基于矿井生产环境,构建了相应的物理模型、数学模型,进行了数值模拟。模拟过程中,考虑了化 学反应,将爆炸热力学过程和基元反应耦合,采用TVD格式进行数值离散。结果揭示巷道中设置1~3辆矿 车时,瓦斯爆炸冲击波及其火焰波的结构和参数变化特征。数值模拟与实验结果对比分析表明,单一矿车和 多辆矿车在激励效应存在明显不同,矿车数量增加,激励效应更加显著。
2012, 32(1): 51-54.
doi: 10.11883/1001-1455(2012)01-0051-04
摘要:
为了解决大面积钛钢板爆炸焊接工程应用中复合板边界区的局部脱焊现象,利用动力非线性数值 模拟,考察了大幅板爆炸焊接的动态过程,通过分析爆焊过程中不同时刻复板飞行姿态,发现了在大幅板爆炸 焊接过程中存在褶皱变形,发生位置与实际工程中大幅板的脱焊位置一致,指出应特别注意复板的褶皱变形, 并解释了褶皱变形出现的原因,为爆炸焊接工艺技术的改进和发展提供了参考。
为了解决大面积钛钢板爆炸焊接工程应用中复合板边界区的局部脱焊现象,利用动力非线性数值 模拟,考察了大幅板爆炸焊接的动态过程,通过分析爆焊过程中不同时刻复板飞行姿态,发现了在大幅板爆炸 焊接过程中存在褶皱变形,发生位置与实际工程中大幅板的脱焊位置一致,指出应特别注意复板的褶皱变形, 并解释了褶皱变形出现的原因,为爆炸焊接工艺技术的改进和发展提供了参考。
2012, 32(1): 55-60.
doi: 10.11883/1001-1455(2012)01-0055-06
摘要:
基于速度势理论,利用复数变量边界元法对二维楔形体常速入水冲击的时域解进行了数值研究。 以相似解作为数值计算的初始条件,采用时域解射流线性近似处理方法,利用复数变量边界元法进行求解,以 减少计算量和数值误差。深入讨论了扩展坐标系求时域解、射流处理、网格划分和网格更新等关键技术。最 后数值计算了不同斜升角楔形体入水时的自由液面隆起、射流飞溅和压力分布,经与相似解结果作比较,自由 液面隆起轮廓基本吻合,而压力分布更符合实际情况,从而证明了模型及分析方法的正确性。
基于速度势理论,利用复数变量边界元法对二维楔形体常速入水冲击的时域解进行了数值研究。 以相似解作为数值计算的初始条件,采用时域解射流线性近似处理方法,利用复数变量边界元法进行求解,以 减少计算量和数值误差。深入讨论了扩展坐标系求时域解、射流处理、网格划分和网格更新等关键技术。最 后数值计算了不同斜升角楔形体入水时的自由液面隆起、射流飞溅和压力分布,经与相似解结果作比较,自由 液面隆起轮廓基本吻合,而压力分布更符合实际情况,从而证明了模型及分析方法的正确性。
2012, 32(1): 61-66.
doi: 10.11883/1001-1455(2012)01-0061-06
摘要:
采用波动理论,提出了水下爆炸冲击波防护的设计原则。在此基础上针对水下爆破工程安全防护 的实际情况,对提出的空气隔层衰减冲击波设想进行了实验研究,发现空气隔层能有效衰减冲击波峰值压力。 为了对空气隔层衰减冲击波能量特性进行研究,提出了基于小波包分析水下爆炸冲击波能量分布规律的方 法,发现水下爆炸冲击波能量分布频率比较广,在空气隔层作用下,各频率段的冲击波能量都有很好的衰减作 用,空气隔层对绝大部分频率段的冲击波能量衰减都在50%以上。研究结果表明,空气隔层能在水下爆破工 程安全防护中发挥积极作用,可以有效降低水下爆炸冲击波的破坏效应,具有实用价值。
采用波动理论,提出了水下爆炸冲击波防护的设计原则。在此基础上针对水下爆破工程安全防护 的实际情况,对提出的空气隔层衰减冲击波设想进行了实验研究,发现空气隔层能有效衰减冲击波峰值压力。 为了对空气隔层衰减冲击波能量特性进行研究,提出了基于小波包分析水下爆炸冲击波能量分布规律的方 法,发现水下爆炸冲击波能量分布频率比较广,在空气隔层作用下,各频率段的冲击波能量都有很好的衰减作 用,空气隔层对绝大部分频率段的冲击波能量衰减都在50%以上。研究结果表明,空气隔层能在水下爆破工 程安全防护中发挥积极作用,可以有效降低水下爆炸冲击波的破坏效应,具有实用价值。
2012, 32(1): 67-72.
doi: 10.11883/1001-1455(2012)01-0067-06
摘要:
为了研究装药壳体厚度对水下爆炸冲击波特性的影响,对1kg柱形含铝炸药分别在厚6mm 的钢 壳或硬铝壳装药下进行了水下爆炸实验与数值模拟研究。实验结果表明,相对硬铝壳,钢壳装药的冲击波冲 量、衰减时间常数以及冲击波能都偏大,而冲击波峰值相差不大。针对不同厚度钢壳装药的数值模拟表明,随 着壳体厚度的增加,冲击波参数明显增强,当壳体厚度超过最优值时这种效应减弱;金属壳的存在导致冲击波 峰值爬升产生滞后效应;对一定质量的炸药,存在可有效提高冲击波压力峰值的最优壳厚,填装比可以作为衡 量效果的重要指标。
为了研究装药壳体厚度对水下爆炸冲击波特性的影响,对1kg柱形含铝炸药分别在厚6mm 的钢 壳或硬铝壳装药下进行了水下爆炸实验与数值模拟研究。实验结果表明,相对硬铝壳,钢壳装药的冲击波冲 量、衰减时间常数以及冲击波能都偏大,而冲击波峰值相差不大。针对不同厚度钢壳装药的数值模拟表明,随 着壳体厚度的增加,冲击波参数明显增强,当壳体厚度超过最优值时这种效应减弱;金属壳的存在导致冲击波 峰值爬升产生滞后效应;对一定质量的炸药,存在可有效提高冲击波压力峰值的最优壳厚,填装比可以作为衡 量效果的重要指标。
2012, 32(1): 73-78.
doi: 10.11883/1001-1455(2012)01-0073-06
摘要:
采用共节点分离式钢筋混凝土模型,通过比较分析塔体缺口中间有切缝与无切缝设计2种方案, 对典型高耸筒形结构冷却塔爆破拆除过程进行了三维数值模拟。对混凝土和钢筋单元的受力过程进行分析, 共节点分离式模型可以体现混凝土和钢筋材料的力学性能差异;采取切缝设计可减小结构倒塌过程中的下 坐与后坐倾向,促使塔体后部严重扭曲变形,使塔身解体比较完全,缩短整体倒塌时间,减小爆堆范围。
采用共节点分离式钢筋混凝土模型,通过比较分析塔体缺口中间有切缝与无切缝设计2种方案, 对典型高耸筒形结构冷却塔爆破拆除过程进行了三维数值模拟。对混凝土和钢筋单元的受力过程进行分析, 共节点分离式模型可以体现混凝土和钢筋材料的力学性能差异;采取切缝设计可减小结构倒塌过程中的下 坐与后坐倾向,促使塔体后部严重扭曲变形,使塔身解体比较完全,缩短整体倒塌时间,减小爆堆范围。
2012, 32(1): 79-83.
doi: 10.11883/1001-1455(2012)01-0079-05
摘要:
根据Held引爆理论,利用聚能射流使炸药只燃烧或发生不完全爆轰的特点,设计了微声聚能穿孔 装置,并对装药结构进行了正交优化。实验结果表明,装置的消声能力良好,用于销毁模拟爆炸药盒时,装药 仅发生低速爆轰,达到了安全销毁简易爆炸装置的目的。
根据Held引爆理论,利用聚能射流使炸药只燃烧或发生不完全爆轰的特点,设计了微声聚能穿孔 装置,并对装药结构进行了正交优化。实验结果表明,装置的消声能力良好,用于销毁模拟爆炸药盒时,装药 仅发生低速爆轰,达到了安全销毁简易爆炸装置的目的。
2012, 32(1): 84-90.
doi: 10.11883/1001-1455(2012)01-0084-07
摘要:
根据场地地震波的传播叠加原理,以实测单炮孔爆破振动波形为基础,考虑预测点位置与各炮孔 的相对位置关系,并按照实际起爆网路设计的各炮孔起爆时差和实测的地震波传播速度等参数,计算获得预 测点的爆破振动波形。不仅可以预测爆破振动速度峰值,而且可以预测完整的振动波形,并可获知爆破振动 持续时间及主振频率分布范围。根据现场应用数码电子雷管的深孔爆破实验,该方法计算的预测波形与实测 波形相当吻合,计算结果可靠性较好,可以在实际工程中推广使用。
根据场地地震波的传播叠加原理,以实测单炮孔爆破振动波形为基础,考虑预测点位置与各炮孔 的相对位置关系,并按照实际起爆网路设计的各炮孔起爆时差和实测的地震波传播速度等参数,计算获得预 测点的爆破振动波形。不仅可以预测爆破振动速度峰值,而且可以预测完整的振动波形,并可获知爆破振动 持续时间及主振频率分布范围。根据现场应用数码电子雷管的深孔爆破实验,该方法计算的预测波形与实测 波形相当吻合,计算结果可靠性较好,可以在实际工程中推广使用。
2012, 32(1): 91-96.
doi: 10.11883/1001-1455(2012)01-0091-06
摘要:
根据材料不可压缩假设,适当选取材料塑性本构关系,将圆柱壳塑性膨胀运动简化为以圆柱壳内 半径为因变量的常微分方程初值问题。对3种初始内外半径比纯铜厚壁圆柱壳膨胀过程的计算表明,应变率 硬化阻碍圆柱壳的塑性膨胀运动,应变硬化和温度效应的影响可以忽略,而与圆柱壳的厚度无关。
根据材料不可压缩假设,适当选取材料塑性本构关系,将圆柱壳塑性膨胀运动简化为以圆柱壳内 半径为因变量的常微分方程初值问题。对3种初始内外半径比纯铜厚壁圆柱壳膨胀过程的计算表明,应变率 硬化阻碍圆柱壳的塑性膨胀运动,应变硬化和温度效应的影响可以忽略,而与圆柱壳的厚度无关。
2012, 32(1): 97-102.
doi: 10.11883/1001-1455(2012)01-0097-06
摘要:
利用由轻气炮装置改装而成的生物激波管,开展不同冲击波超压值下兔的损伤效应研究。通过对 兔的多种生理特征和器官(如肺、肾脏等)的变化情况进行检测,综合分析了生物体生理系统对冲击波响应及 其损伤后的自身调节作用。建立了损伤肺湿干比与超压间的定量关系,发现了血压随冲击波增加而降低的现 象。除肺外,冲击波对生物体其他脏器同样存在显著影响。
利用由轻气炮装置改装而成的生物激波管,开展不同冲击波超压值下兔的损伤效应研究。通过对 兔的多种生理特征和器官(如肺、肾脏等)的变化情况进行检测,综合分析了生物体生理系统对冲击波响应及 其损伤后的自身调节作用。建立了损伤肺湿干比与超压间的定量关系,发现了血压随冲击波增加而降低的现 象。除肺外,冲击波对生物体其他脏器同样存在显著影响。
2012, 32(1): 103-107.
doi: 10.11883/1001-1455(2012)01-0103-05
摘要:
以3mm 厚的TA2钛板和26mm 厚的正火态Q345R为材料,通过爆炸焊接实验,对钛/钢复合板 爆炸焊接的动态参数进行了研究。结合复合板结合界面特征、复合板结合强度(剪切强度)以及界面波的金相 组织,讨论了钛/钢爆炸焊接时获得高强度结合和规则的界面结合波状形态的条件。对于3mm 厚TA2与 26mm厚正火态Q345R,该条件是动态碰撞角17,动态碰撞速度vp760m/s。根据界面波及基板轧制 金相组织形态,分析了形成界面波的机理,认为射流阻碍复板连续碰撞基板是形成界面波的一个主要原因。
以3mm 厚的TA2钛板和26mm 厚的正火态Q345R为材料,通过爆炸焊接实验,对钛/钢复合板 爆炸焊接的动态参数进行了研究。结合复合板结合界面特征、复合板结合强度(剪切强度)以及界面波的金相 组织,讨论了钛/钢爆炸焊接时获得高强度结合和规则的界面结合波状形态的条件。对于3mm 厚TA2与 26mm厚正火态Q345R,该条件是动态碰撞角17,动态碰撞速度vp760m/s。根据界面波及基板轧制 金相组织形态,分析了形成界面波的机理,认为射流阻碍复板连续碰撞基板是形成界面波的一个主要原因。
2012, 32(1): 108-112.
doi: 10.11883/1001-1455(2012)01-0108-05
摘要:
为深入了解液体燃料云雾爆轰机理,借助自行设计和建造的立式爆轰管,并采用升降法和烟迹技 术,对环氧丙烷等液体燃料云雾爆轰参数(爆速、爆压、临界起爆能和爆轰胞格尺寸)与当量比的关系进行了实 验研究。结果表明,环氧丙烷的爆速和爆压随当量比的增加先增大后平缓减小;碳氢液体燃料云雾爆轰的临 界起爆能与当量比呈U型关系,最佳值点偏向富燃料一侧;临界起爆能的大小与燃料的分子结构和挥发性 有密切关系,IPN和PO 临界起爆能相当,而烷烃类临界起爆能均较大。环氧丙烷在25和50 ℃时的爆轰胞 格宽度与当量比皆呈U型关系,最小值点偏向富燃料一侧;气相爆轰胞格宽度比云雾爆轰略小。常温下环 氧丙烷云雾爆轰主要由气相反应所控制。
为深入了解液体燃料云雾爆轰机理,借助自行设计和建造的立式爆轰管,并采用升降法和烟迹技 术,对环氧丙烷等液体燃料云雾爆轰参数(爆速、爆压、临界起爆能和爆轰胞格尺寸)与当量比的关系进行了实 验研究。结果表明,环氧丙烷的爆速和爆压随当量比的增加先增大后平缓减小;碳氢液体燃料云雾爆轰的临 界起爆能与当量比呈U型关系,最佳值点偏向富燃料一侧;临界起爆能的大小与燃料的分子结构和挥发性 有密切关系,IPN和PO 临界起爆能相当,而烷烃类临界起爆能均较大。环氧丙烷在25和50 ℃时的爆轰胞 格宽度与当量比皆呈U型关系,最小值点偏向富燃料一侧;气相爆轰胞格宽度比云雾爆轰略小。常温下环 氧丙烷云雾爆轰主要由气相反应所控制。
