2015年 35卷 第3期
2015, 35(3): 289-295.
doi: 10.11883/1001-1455-(2015)03-0289-07
在线阅读
PDF 
本文被引
摘要:
基于强爆炸火球光辐射的多群辐射流体力学方法, 采用算子分裂方法将方程组分裂为对流项和刚性源项, 其中源项部分根据方程形式, 进一步分裂为各群内的单独求解。数值计算表明:该方法克服了直接求解过程中辐射与流体耦合所带来的强不稳定性, 时间步长大幅提高, 给出的火球光辐射能谱特征与已有规律一致。可为定量分析光辐射能谱特征提供有效手段。
基于强爆炸火球光辐射的多群辐射流体力学方法, 采用算子分裂方法将方程组分裂为对流项和刚性源项, 其中源项部分根据方程形式, 进一步分裂为各群内的单独求解。数值计算表明:该方法克服了直接求解过程中辐射与流体耦合所带来的强不稳定性, 时间步长大幅提高, 给出的火球光辐射能谱特征与已有规律一致。可为定量分析光辐射能谱特征提供有效手段。
2015, 35(3): 296-303.
doi: 10.11883/1001-1455-(2015)03-0296-08
摘要:
为了解TiNi柱壳横向压缩力学性能以制造可重复使用抗冲吸能装置, 对有、无侧向约束的TiNi柱壳进行了横向冲击实验。利用改进的霍普金森压杆装置(SHPB), 配套波形分离方法, 实现了较长时间(~3 ms)的波形测量, 获得了TiNi柱壳在动态加载下的载荷压缩量曲线。通过高速摄影, 捕捉了柱壳的动态变形过程。结果表明, 无约束试件具有优良的可恢复变形能力, 承载力平台段特征明显。侧向约束的引入, 可以有效提高柱壳的承载力和耗能能力, 可以承受更高速度的冲击。选择合适的约束组合, 可望同时实现较大压缩行程和高耗能, 制造实用的抗冲吸能装置。
为了解TiNi柱壳横向压缩力学性能以制造可重复使用抗冲吸能装置, 对有、无侧向约束的TiNi柱壳进行了横向冲击实验。利用改进的霍普金森压杆装置(SHPB), 配套波形分离方法, 实现了较长时间(~3 ms)的波形测量, 获得了TiNi柱壳在动态加载下的载荷压缩量曲线。通过高速摄影, 捕捉了柱壳的动态变形过程。结果表明, 无约束试件具有优良的可恢复变形能力, 承载力平台段特征明显。侧向约束的引入, 可以有效提高柱壳的承载力和耗能能力, 可以承受更高速度的冲击。选择合适的约束组合, 可望同时实现较大压缩行程和高耗能, 制造实用的抗冲吸能装置。
2015, 35(3): 304-311.
doi: 10.11883/1001-1455-(2015)03-0304-08
摘要:
采用新型数字激光动态焦散线实验系统, 对爆炸荷载作用下空孔周围的动应力场分布及空孔对爆生主裂纹扩展行为的影响进行了研究。研究结果表明, 在空孔周围强应力场的影响下, 2条相向扩展的爆生主裂纹逐渐向空孔处偏转, 并在空孔处贯通; 空孔附近的主应力方向与炮孔连心线夹角基本稳定在约12°, 增大空孔尺寸对空孔附近的主应力方向影响不明显; 爆炸应力波与空孔相互作用, 产生反射拉伸波, 改变了主裂纹尖端的应力场, 降低了主裂纹的扩展速度, 且空孔直径越大, 主裂纹的扩展速度越低; 当爆生主裂纹扩展到空孔附近时, 主裂纹尖端动态应力强度因子再次出现上升的趋势。
采用新型数字激光动态焦散线实验系统, 对爆炸荷载作用下空孔周围的动应力场分布及空孔对爆生主裂纹扩展行为的影响进行了研究。研究结果表明, 在空孔周围强应力场的影响下, 2条相向扩展的爆生主裂纹逐渐向空孔处偏转, 并在空孔处贯通; 空孔附近的主应力方向与炮孔连心线夹角基本稳定在约12°, 增大空孔尺寸对空孔附近的主应力方向影响不明显; 爆炸应力波与空孔相互作用, 产生反射拉伸波, 改变了主裂纹尖端的应力场, 降低了主裂纹的扩展速度, 且空孔直径越大, 主裂纹的扩展速度越低; 当爆生主裂纹扩展到空孔附近时, 主裂纹尖端动态应力强度因子再次出现上升的趋势。
2015, 35(3): 312-319.
doi: 10.11883/1001-1455-(2015)03-0312-08
摘要:
对MDYB-3有机玻璃进行了多组不同应变率(10-3~3 000 s-1)下的压缩实验, 得到准静态下的屈服应力与动态下的峰值应力。沿其增强与面内2个方向进行准静态压缩实验, 以分析定向拉伸对屈服应力的影响, 修正了Ree-Eyring模型与Cooperative模型以描述定向有机玻璃的屈服行为。采用Johnson-Cook模型描述屈服后的黏塑性行为。结果表明Cooperative屈服模型比Ree-Eyring屈服模型更接近实验结果, 且能准确描述准静态屈服应力。动态压缩下的峰值应力为失效应力, 说明试样在1 500 s-1以上应变率下未达到屈服应力时已经发生破坏。Johnson-Cook模型对于单条曲线拟合良好, 但无法准确描述材料的应变率相关性。
对MDYB-3有机玻璃进行了多组不同应变率(10-3~3 000 s-1)下的压缩实验, 得到准静态下的屈服应力与动态下的峰值应力。沿其增强与面内2个方向进行准静态压缩实验, 以分析定向拉伸对屈服应力的影响, 修正了Ree-Eyring模型与Cooperative模型以描述定向有机玻璃的屈服行为。采用Johnson-Cook模型描述屈服后的黏塑性行为。结果表明Cooperative屈服模型比Ree-Eyring屈服模型更接近实验结果, 且能准确描述准静态屈服应力。动态压缩下的峰值应力为失效应力, 说明试样在1 500 s-1以上应变率下未达到屈服应力时已经发生破坏。Johnson-Cook模型对于单条曲线拟合良好, 但无法准确描述材料的应变率相关性。
2015, 35(3): 320-325.
doi: 10.11883/1001-1455-(2015)03-0320-06
摘要:
为在实验室内开展深水爆炸气泡脉动规律研究, 通过增加水面大气压强来模拟水中静水压的方法, 建立可模拟深水环境的爆炸容器。开展不同模拟水深条件下的3种当量炸药的水下爆炸实验, 得到了气泡脉动过程图像, 验证小当量深水爆炸模拟实验与自由场实验的等效性, 分析气泡脉动周期与最大半径同模拟水深的关系。实验结果表明:容器壁面反射效应对气泡脉动过程的影响可以忽略不计, 模拟实验可等效为自由场实验; 深水爆炸气泡脉动周期及最大半径随流体静力深度增加的衰减系数分别为-0.83和-0.364。
为在实验室内开展深水爆炸气泡脉动规律研究, 通过增加水面大气压强来模拟水中静水压的方法, 建立可模拟深水环境的爆炸容器。开展不同模拟水深条件下的3种当量炸药的水下爆炸实验, 得到了气泡脉动过程图像, 验证小当量深水爆炸模拟实验与自由场实验的等效性, 分析气泡脉动周期与最大半径同模拟水深的关系。实验结果表明:容器壁面反射效应对气泡脉动过程的影响可以忽略不计, 模拟实验可等效为自由场实验; 深水爆炸气泡脉动周期及最大半径随流体静力深度增加的衰减系数分别为-0.83和-0.364。
2015, 35(3): 326-334.
doi: 10.11883/1001-1455-(2015)03-0326-09
摘要:
针对不可压缩可压缩水/气多介质问题, 提出一种新的界面处理方法。在可压缩水/气界面处构造Riemann问题, 在水中设音速趋于无穷大, 求解Riemann问题得到不可压缩可压缩水/气界面处流体的准确流动状态; 然后以此状态结合GFM(ghost fluid method)方法分别为2种流体定义界面边界条件, 将两相流问题转化为单相流问题计算, 通过求解level set方程来跟踪界面的位置。对各种不同的界面边界条件定义方法进行了比较, 数值模拟结果表明算法能准确地捕捉各类间断的位置, 证明了算法的有效性和稳健性。
针对不可压缩可压缩水/气多介质问题, 提出一种新的界面处理方法。在可压缩水/气界面处构造Riemann问题, 在水中设音速趋于无穷大, 求解Riemann问题得到不可压缩可压缩水/气界面处流体的准确流动状态; 然后以此状态结合GFM(ghost fluid method)方法分别为2种流体定义界面边界条件, 将两相流问题转化为单相流问题计算, 通过求解level set方程来跟踪界面的位置。对各种不同的界面边界条件定义方法进行了比较, 数值模拟结果表明算法能准确地捕捉各类间断的位置, 证明了算法的有效性和稳健性。
2015, 35(3): 335-342.
doi: 10.11883/1001-1455-(2015)03-0335-08
摘要:
为研究起爆不对称性对EFP战斗部飞行特性的影响, 对不同偏心量下⌀60 mm弧锥结合罩EFP战斗部进行飞行弹道实验。实验结果表明:偏心起爆条件下, 当相对偏心量小于3.3%时, EFP在网靶穿孔接近圆形, 弹丸飞行稳定; 起爆相对偏心量达到6.7%时, 弹丸飞行过程中摆动幅值增大, 降低了对目标的打击精度和毁伤效果。利用LS-DYNA及CFX非线性动力学有限元程序对不同起爆偏心量下成型EFP的空气动力学特性进行数值模拟, 描述了偏心起爆影响EFP成型对称性, 改变弹丸在飞行过程中流场的分布特征, 从而导致弹丸飞行过程中无规则运动的全过程。
为研究起爆不对称性对EFP战斗部飞行特性的影响, 对不同偏心量下⌀60 mm弧锥结合罩EFP战斗部进行飞行弹道实验。实验结果表明:偏心起爆条件下, 当相对偏心量小于3.3%时, EFP在网靶穿孔接近圆形, 弹丸飞行稳定; 起爆相对偏心量达到6.7%时, 弹丸飞行过程中摆动幅值增大, 降低了对目标的打击精度和毁伤效果。利用LS-DYNA及CFX非线性动力学有限元程序对不同起爆偏心量下成型EFP的空气动力学特性进行数值模拟, 描述了偏心起爆影响EFP成型对称性, 改变弹丸在飞行过程中流场的分布特征, 从而导致弹丸飞行过程中无规则运动的全过程。
2015, 35(3): 343-349.
doi: 10.11883/1001-1455-(2015)03-0343-07
摘要:
针对深部岩体中由断层、节理等不连续性结构面引发的岩爆地质灾害, 根据深埋地下隧洞中潜在发震断裂的分布特征和几何形态建立数值分析模型, 采用离散元单元法模拟存在刚性平直断裂的深部围岩的开挖响应, 并分别考察开挖接近并通过断裂附近时围岩应力状态的变化特征。通过探讨断裂的存在对围岩应力状态改变的作用机理, 揭示出断裂型岩爆是开挖面附近一定范围内存在的断裂构造在高应力作用下发生错动, 导致能量突然释放, 对围岩造成强烈冲击作用的结果, 基本与地震的断层“粘滑”机制相类似。
针对深部岩体中由断层、节理等不连续性结构面引发的岩爆地质灾害, 根据深埋地下隧洞中潜在发震断裂的分布特征和几何形态建立数值分析模型, 采用离散元单元法模拟存在刚性平直断裂的深部围岩的开挖响应, 并分别考察开挖接近并通过断裂附近时围岩应力状态的变化特征。通过探讨断裂的存在对围岩应力状态改变的作用机理, 揭示出断裂型岩爆是开挖面附近一定范围内存在的断裂构造在高应力作用下发生错动, 导致能量突然释放, 对围岩造成强烈冲击作用的结果, 基本与地震的断层“粘滑”机制相类似。
2015, 35(3): 350-358.
doi: 10.11883/1001-1455(2015)03-0350-09
摘要:
为适应城市密集建筑物下隧道爆破对振动的高安全要求, 研究用普通爆破器材进行振速精确控制的爆破技术和参数确定方法。以渝中隧道为研究背景, 在开发准确延时非电雷管的基础上, 利用傅立叶函数和MATLAB软件拟合了不同药量单孔爆破振动波形, 分析了1~50 ms不同间隔下振动叠加的量化数据; 讨论了各微差间隔时间的降振效果; 在指定振速的情况下, 确定单孔药量和微差起爆时间; 实测并分析了现场使用雷管各段微差间隔特点, 据此进行针对性的爆破设计和采用逐孔掏槽爆破进行振速控制。现场应用表明:隧道爆破振速始终小于1.00 cm/s, 在此振速下避免振动叠加的最优单孔药量为1.2 kg, 爆破振速峰值位于主掏槽的第1段或第2段雷管起爆后, 且与理论分析结果吻合较好; 逐孔起爆60 ms后振速下降50%以上。研究表明:在高安全指标下, 以非电雷管实施精确控制爆破是可以实现的。
为适应城市密集建筑物下隧道爆破对振动的高安全要求, 研究用普通爆破器材进行振速精确控制的爆破技术和参数确定方法。以渝中隧道为研究背景, 在开发准确延时非电雷管的基础上, 利用傅立叶函数和MATLAB软件拟合了不同药量单孔爆破振动波形, 分析了1~50 ms不同间隔下振动叠加的量化数据; 讨论了各微差间隔时间的降振效果; 在指定振速的情况下, 确定单孔药量和微差起爆时间; 实测并分析了现场使用雷管各段微差间隔特点, 据此进行针对性的爆破设计和采用逐孔掏槽爆破进行振速控制。现场应用表明:隧道爆破振速始终小于1.00 cm/s, 在此振速下避免振动叠加的最优单孔药量为1.2 kg, 爆破振速峰值位于主掏槽的第1段或第2段雷管起爆后, 且与理论分析结果吻合较好; 逐孔起爆60 ms后振速下降50%以上。研究表明:在高安全指标下, 以非电雷管实施精确控制爆破是可以实现的。
2015, 35(3): 359-365.
doi: 10.11883/1001-1455(2015)03-0359-07
摘要:
针对混凝土两点爆炸起爆参数选择问题, 提出了一种基于灰色理论的参数优化方法。通过正交试验方法设计试验方案, 运用AUTODYN软件进行了不同起爆参数组合条件下的静爆试验, 计算了起爆参数与爆坑直径、爆坑深度的关联系数和关联度, 进行了单目标因素优化和多目标因素优化, 确定了一组各因素的优选组合, 并进行了试验验证。验证结果表明:采用优化的起爆参数时, 爆坑直径增大(4~42)%, 左爆坑深度增大(0~29)%, 右爆坑深度增大(0~32)%, 两点爆炸混凝土靶体的毁伤效果得到明显改善。
针对混凝土两点爆炸起爆参数选择问题, 提出了一种基于灰色理论的参数优化方法。通过正交试验方法设计试验方案, 运用AUTODYN软件进行了不同起爆参数组合条件下的静爆试验, 计算了起爆参数与爆坑直径、爆坑深度的关联系数和关联度, 进行了单目标因素优化和多目标因素优化, 确定了一组各因素的优选组合, 并进行了试验验证。验证结果表明:采用优化的起爆参数时, 爆坑直径增大(4~42)%, 左爆坑深度增大(0~29)%, 右爆坑深度增大(0~32)%, 两点爆炸混凝土靶体的毁伤效果得到明显改善。
2015, 35(3): 366-371.
doi: 10.11883/1001-1455(2015)03-0366-07
摘要:
以高压水射流冲击HTPB推进剂的动态加载过程和准静态加载过程在作用压力和持续时间上的巨大差异为基础, 在水锤压力和滞止压力计算的基础上分别进行了点火模式预判, 然后以模型类比和实验方法分析了动态和准静态加载过程的安全性。结果表明, 使用出口压力在300 MPa以内的高压水射流冲击HTPB推进剂装药在动态加载过程中不会有点火起爆危险性, 但使用100M Pa以上的高压水射流冲击HTPB推进剂装药在准静态加载过程中其内部可能会发生温度突跃情况, 这可能会引起热点火、甚至热起爆。
以高压水射流冲击HTPB推进剂的动态加载过程和准静态加载过程在作用压力和持续时间上的巨大差异为基础, 在水锤压力和滞止压力计算的基础上分别进行了点火模式预判, 然后以模型类比和实验方法分析了动态和准静态加载过程的安全性。结果表明, 使用出口压力在300 MPa以内的高压水射流冲击HTPB推进剂装药在动态加载过程中不会有点火起爆危险性, 但使用100M Pa以上的高压水射流冲击HTPB推进剂装药在准静态加载过程中其内部可能会发生温度突跃情况, 这可能会引起热点火、甚至热起爆。
2015, 35(3): 372-379.
doi: 10.11883/1001-1455-(2015)03-0372-08
摘要:
为探索气流特征对水平长管内粉尘爆炸火焰结构的影响, 对采用加压送气传输方式形成的石松子粉尘云经静电引燃后其火焰在水平长管内的传播特性进行实验。利用热线风速仪测量不同气流条件下沿管径方向的速度分布和湍流强度分布, 采用高速摄像系统记录了火焰在水平管道内的传播过程。实验观察到, 即使管内石松子粉尘质量分数相同, 仍然会出现2种不同类型的火焰结构:一种类型火焰轮廓规则、清晰, 火焰中心为连续的黄色发光区并由红色边缘火焰包裹; 另一种类型火焰空间离散, 火焰发光区局部存在, 散乱地呈现不规则状态。详细分析不同气流条件对火焰结构的影响。
为探索气流特征对水平长管内粉尘爆炸火焰结构的影响, 对采用加压送气传输方式形成的石松子粉尘云经静电引燃后其火焰在水平长管内的传播特性进行实验。利用热线风速仪测量不同气流条件下沿管径方向的速度分布和湍流强度分布, 采用高速摄像系统记录了火焰在水平管道内的传播过程。实验观察到, 即使管内石松子粉尘质量分数相同, 仍然会出现2种不同类型的火焰结构:一种类型火焰轮廓规则、清晰, 火焰中心为连续的黄色发光区并由红色边缘火焰包裹; 另一种类型火焰空间离散, 火焰发光区局部存在, 散乱地呈现不规则状态。详细分析不同气流条件对火焰结构的影响。
2015, 35(3): 380-385.
doi: 10.11883/1001-1455-(2015)03-0380-06
摘要:
在方管的基础上引入折纹结构, 利用几何关系建立折纹管的折角公式。采用LS-DYNA软件研究了6种折纹管在轴向冲击下的屈曲模态与能量吸收性能, 并与方管进行对比分析。结果表明, 折纹管在冲击载荷作用下屈曲变形过程可分为3个阶段, 初始峰值阶段、稳定渐进屈曲阶段和密实化阶段。折角是影响初始峰值载荷和平均载荷的重要因素之一, 折纹结构的引入有效的降低了初始峰值载荷, 减小了冲击力的波动幅度; 折纹管的比吸能低于方管, 但是在特定折角下, 折纹管的压缩力效率和比总体效率高于方管。
在方管的基础上引入折纹结构, 利用几何关系建立折纹管的折角公式。采用LS-DYNA软件研究了6种折纹管在轴向冲击下的屈曲模态与能量吸收性能, 并与方管进行对比分析。结果表明, 折纹管在冲击载荷作用下屈曲变形过程可分为3个阶段, 初始峰值阶段、稳定渐进屈曲阶段和密实化阶段。折角是影响初始峰值载荷和平均载荷的重要因素之一, 折纹结构的引入有效的降低了初始峰值载荷, 减小了冲击力的波动幅度; 折纹管的比吸能低于方管, 但是在特定折角下, 折纹管的压缩力效率和比总体效率高于方管。
2015, 35(3): 386-392.
doi: 10.11883/1001-1455(2015)03-0386-07
摘要:
利用轻气炮进行了半球形头杆弹正撞击单层板和等厚接触式三层板的实验, 得到了这两种结构靶体的初始-剩余速度曲线以及弹道极限。采用ABAQUS/EXPLICIT数值模拟软件对杆弹撞击金属板的过程进行了数值模拟研究, 通过对比数值模拟和实验结果, 验证了数值模拟材料模型和参数的有效性。研究了靶体结构对抗侵彻特性的影响, 并分析了弹体对靶体的撞击过程。研究结果表明:多层板的弹道极限高于等厚单层板。单层板主要失效模式为剪切, 而多层板的主要失效模式为整体的蝶形变形和局部的盘式隆起。对于多层板, 靶板具体的失效模式与其在靶中位置相关。
利用轻气炮进行了半球形头杆弹正撞击单层板和等厚接触式三层板的实验, 得到了这两种结构靶体的初始-剩余速度曲线以及弹道极限。采用ABAQUS/EXPLICIT数值模拟软件对杆弹撞击金属板的过程进行了数值模拟研究, 通过对比数值模拟和实验结果, 验证了数值模拟材料模型和参数的有效性。研究了靶体结构对抗侵彻特性的影响, 并分析了弹体对靶体的撞击过程。研究结果表明:多层板的弹道极限高于等厚单层板。单层板主要失效模式为剪切, 而多层板的主要失效模式为整体的蝶形变形和局部的盘式隆起。对于多层板, 靶板具体的失效模式与其在靶中位置相关。
2015, 35(3): 393-400.
doi: 10.11883/1001-1455-(2015)03-0393-08
摘要:
为了有效防治矿井瓦斯爆炸事故, 以瓦斯的主要成分甲烷作为模拟气体, 运用自主设计改装的XKWB-S型小尺寸石英玻璃管道实验系统, 结合高速摄影仪, 并采用FLACS数值模拟软件, 研究惰性气体抑爆条件下甲烷燃烧爆炸特性, 进行体积分数为6%~27%的CO2抑制体积分数为9%CH4爆炸的实验及数值模拟, 结果表明:各组分混合气体在爆炸传播过程中, 爆炸压力、火焰锋面速度和气体运动速度均呈现一定程度的波动, 且压力和速度没有同时达到最大值; CO2的加入有效抑制了甲烷/空气反应, 且添加CO2体积分数越大, 抑爆效果越明显, 模拟结果与实验结果基本吻合。
为了有效防治矿井瓦斯爆炸事故, 以瓦斯的主要成分甲烷作为模拟气体, 运用自主设计改装的XKWB-S型小尺寸石英玻璃管道实验系统, 结合高速摄影仪, 并采用FLACS数值模拟软件, 研究惰性气体抑爆条件下甲烷燃烧爆炸特性, 进行体积分数为6%~27%的CO2抑制体积分数为9%CH4爆炸的实验及数值模拟, 结果表明:各组分混合气体在爆炸传播过程中, 爆炸压力、火焰锋面速度和气体运动速度均呈现一定程度的波动, 且压力和速度没有同时达到最大值; CO2的加入有效抑制了甲烷/空气反应, 且添加CO2体积分数越大, 抑爆效果越明显, 模拟结果与实验结果基本吻合。
2015, 35(3): 401-408.
doi: 10.11883/1001-1455-(2015)03-0401-08
摘要:
研究多孔材料细观结构与宏观力学性能之间的关系, 建立具有固定相对密度的含随机固体填充孔的圆形蜂窝结构模型。在此模型的基础上具体讨论了不同孔洞填充比和冲击速度对圆形蜂窝结构变形模式、动态冲击平台应力以及能量吸收性能的影响。研究结果表明:填充孔在蜂窝变形过程中有局部牵制作用, 蜂窝材料变形模式仍为准静态模式、过渡模式、动态模式; 当变形模式为过渡模式或动态模式时, 结构的平台应力与速度的平方成线性关系, 存在明显的速度效应; 高速冲击下, 含固体填充孔的蜂窝结构单位质量吸收的能量高于规则蜂窝结构。研究结果可为蜂窝材料的研究和设计提供参考。
研究多孔材料细观结构与宏观力学性能之间的关系, 建立具有固定相对密度的含随机固体填充孔的圆形蜂窝结构模型。在此模型的基础上具体讨论了不同孔洞填充比和冲击速度对圆形蜂窝结构变形模式、动态冲击平台应力以及能量吸收性能的影响。研究结果表明:填充孔在蜂窝变形过程中有局部牵制作用, 蜂窝材料变形模式仍为准静态模式、过渡模式、动态模式; 当变形模式为过渡模式或动态模式时, 结构的平台应力与速度的平方成线性关系, 存在明显的速度效应; 高速冲击下, 含固体填充孔的蜂窝结构单位质量吸收的能量高于规则蜂窝结构。研究结果可为蜂窝材料的研究和设计提供参考。
2015, 35(3): 409-415.
doi: 10.11883/1001-1455-(2015)03-0409-07
摘要:
为了分析多级渐扩型燃烧室结构对燃烧轻气炮氢氧燃烧特性的影响, 通过计算流体力学方法, 分别对采用传统圆柱型燃烧室和多级渐扩型燃烧室的燃烧轻气炮氢氧燃烧发射过程进行数值模拟。对比结果表明, 多级渐扩型燃烧室结构能够明显地减小燃烧室压力波动幅度, 提高氢氧燃烧稳定性; 多级渐扩型燃烧室内形成回流区, 可以减小气流轴向运动速度; 火焰扩展形态与渐扩型结构相吻合, 燃烧反应区表面变化平稳; 多级渐扩型燃烧室结构对氢氧火焰传播过程和压力波动现象有着重要影响。
为了分析多级渐扩型燃烧室结构对燃烧轻气炮氢氧燃烧特性的影响, 通过计算流体力学方法, 分别对采用传统圆柱型燃烧室和多级渐扩型燃烧室的燃烧轻气炮氢氧燃烧发射过程进行数值模拟。对比结果表明, 多级渐扩型燃烧室结构能够明显地减小燃烧室压力波动幅度, 提高氢氧燃烧稳定性; 多级渐扩型燃烧室内形成回流区, 可以减小气流轴向运动速度; 火焰扩展形态与渐扩型结构相吻合, 燃烧反应区表面变化平稳; 多级渐扩型燃烧室结构对氢氧火焰传播过程和压力波动现象有着重要影响。
2015, 35(3): 416-422.
doi: 10.11883/1001-1455-(2015)03-0416-07
摘要:
鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。
鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。
2015, 35(3): 423-427.
doi: 10.11883/1001-1455(2015)03-0423-05
摘要:
采用浮阻力模型对激波管低压缩和激光加载高压缩情况下的Richtmyer-Meshkov不稳定性诱导混合现象进行了研究。通过与实验和理论分析结果进行比较发现:为了达到好的吻合, Richtmyer-Meshkov不稳定性情况下阻力系数的取值范围(2.0~5.36)比Rayleigh-Taylor不稳定性情况下的值(3.3~4.0)宽得多; 而在Richtmyer-Meshkov不稳定性情况下, 高压缩时阻力系数的不确定度(约为3.36)明显高于低压缩时的值(约为1.46), 模型的进一步完善还有待于更精确实验的验证。研究显示:指数律经验公式中指数随工况的不同而显著变化, 目前工程设计中采用指数律经验公式是粗糙的。
采用浮阻力模型对激波管低压缩和激光加载高压缩情况下的Richtmyer-Meshkov不稳定性诱导混合现象进行了研究。通过与实验和理论分析结果进行比较发现:为了达到好的吻合, Richtmyer-Meshkov不稳定性情况下阻力系数的取值范围(2.0~5.36)比Rayleigh-Taylor不稳定性情况下的值(3.3~4.0)宽得多; 而在Richtmyer-Meshkov不稳定性情况下, 高压缩时阻力系数的不确定度(约为3.36)明显高于低压缩时的值(约为1.46), 模型的进一步完善还有待于更精确实验的验证。研究显示:指数律经验公式中指数随工况的不同而显著变化, 目前工程设计中采用指数律经验公式是粗糙的。
2015, 35(3): 428-436.
doi: 10.11883/1001-1455-(2015)03-0428-09
摘要:
针对节理岩体同时含有节理、裂隙等宏观缺陷及微裂隙、微孔洞等细观缺陷的客观事实, 提出了在节理岩体动态损伤本构模型中应同时考虑宏细观缺陷的观点。为此, 首先对基于细观动态断裂机理的经典岩石动态损伤本构模型—TCK(Taylor-Chen-Kuszmaul)模型进行了阐述, 其次基于Lemaitre等效应变假设推导了综合考虑宏细观缺陷的复合损伤变量(张量), 进而在此基础上建立了相应的节理岩体动态损伤本构模型, 并利用该模型讨论了载荷应变率及节理条数对岩体动态力学特性的影响规律。结果表明, 在不同载荷应变率下试件在变形初始阶段是重合的, 而后随着应变的增加, 试件峰值强度、峰值应变及总应变均随载荷应变率的增加而增加; 随着节理条数的增加, 试件峰值强度逐渐降低, 但降低趋势逐渐变缓并趋于某一定值。上述研究结论与目前的理论及实验研究结果的基本规律是一致的, 说明了本模型的合理性。
针对节理岩体同时含有节理、裂隙等宏观缺陷及微裂隙、微孔洞等细观缺陷的客观事实, 提出了在节理岩体动态损伤本构模型中应同时考虑宏细观缺陷的观点。为此, 首先对基于细观动态断裂机理的经典岩石动态损伤本构模型—TCK(Taylor-Chen-Kuszmaul)模型进行了阐述, 其次基于Lemaitre等效应变假设推导了综合考虑宏细观缺陷的复合损伤变量(张量), 进而在此基础上建立了相应的节理岩体动态损伤本构模型, 并利用该模型讨论了载荷应变率及节理条数对岩体动态力学特性的影响规律。结果表明, 在不同载荷应变率下试件在变形初始阶段是重合的, 而后随着应变的增加, 试件峰值强度、峰值应变及总应变均随载荷应变率的增加而增加; 随着节理条数的增加, 试件峰值强度逐渐降低, 但降低趋势逐渐变缓并趋于某一定值。上述研究结论与目前的理论及实验研究结果的基本规律是一致的, 说明了本模型的合理性。
2015, 35(3): 437-441.
doi: 10.11883/1001-1455(2015)03-0437-05
摘要:
在爆炸容器中进行小药量空中爆炸实验, 利用传感器序列测量冲击波速度, 根据冲击波Rankine-Hugoniot关系获得测点近似理论峰值压力, 从而实现压力传感器的标定, 获得的灵敏度相对误差较小。同时测量了相应的冲击波参数, 并利用Modified-Friedlander公式进行数据后处理, 结果表明固定超压拟合更接近物理事实, 固定正相时间拟合也具有较高精度。最后进行了误差分析, 发现不同传感器特性及数据后处理方法都会带来一定误差。实验结果表明这种测量和后处理方法具有较高的精度, 可以同时标定传感器和测量冲击波参数。
在爆炸容器中进行小药量空中爆炸实验, 利用传感器序列测量冲击波速度, 根据冲击波Rankine-Hugoniot关系获得测点近似理论峰值压力, 从而实现压力传感器的标定, 获得的灵敏度相对误差较小。同时测量了相应的冲击波参数, 并利用Modified-Friedlander公式进行数据后处理, 结果表明固定超压拟合更接近物理事实, 固定正相时间拟合也具有较高精度。最后进行了误差分析, 发现不同传感器特性及数据后处理方法都会带来一定误差。实验结果表明这种测量和后处理方法具有较高的精度, 可以同时标定传感器和测量冲击波参数。
2015, 35(3): 442-448.
doi: 10.11883/1001-1455(2015)03-0442-07
摘要:
以质量守恒与动量守恒定律为基础, 建立了高压水射流冲击煤体的力学模型。运用此模型分析了高压水射流在冲击煤体的过程中, 未破水体、破碎水体、煤体的破碎区与扩孔区的力学特征, 利用严格的力学守恒关系得出高压水射流冲击煤体的简化常微分方程组。将理论计算结果与现场实验和数值模拟结果进行对比, 结果表明:理论计算结果与数值模拟结果和实验结果基本一致。此模型具有明确的力学意义, 且能够反映真实的冲击过程。
以质量守恒与动量守恒定律为基础, 建立了高压水射流冲击煤体的力学模型。运用此模型分析了高压水射流在冲击煤体的过程中, 未破水体、破碎水体、煤体的破碎区与扩孔区的力学特征, 利用严格的力学守恒关系得出高压水射流冲击煤体的简化常微分方程组。将理论计算结果与现场实验和数值模拟结果进行对比, 结果表明:理论计算结果与数值模拟结果和实验结果基本一致。此模型具有明确的力学意义, 且能够反映真实的冲击过程。
