【摘 要】
：
由于对流-扩散-反应方程中对流项、扩散项和反应项的时间尺度不一致，特别是反应项尺度与前两者存在很大差异，因此常微分方程存在很强的刚性。采用L-稳定的IMEX (implicit-explicit )型Additive Runge-Kutta (ARK )方法处理CDR 方程，对对流-扩散项进行显式处理，反应项则用隐式格式处理。着重于对ARK 格式的稳定性进行分析，同时对多个算例包括一维气相爆轰波进
【机 构】
：
北京理工大学爆炸科学与技术国家重点实验室,北京,100081
【出 处】
：
中国力学大会2011暨钱学森诞辰100周年纪念大会
论文部分内容阅读
由于对流-扩散-反应方程中对流项、扩散项和反应项的时间尺度不一致，特别是反应项尺度与前两者存在很大差异，因此常微分方程存在很强的刚性。采用L-稳定的IMEX (implicit-explicit )型Additive Runge-Kutta (ARK )方法处理CDR 方程，对对流-扩散项进行显式处理，反应项则用隐式格式处理。着重于对ARK 格式的稳定性进行分析，同时对多个算例包括一维气相爆轰波进行数值模拟，通过对数值试验结果进行对比分析，该方法能够很好的处理CDR 方程的刚性问题。
其他文献
本文依据爆炸焊接理论，通过建立临界参数的计算模型，对比分析其适用条件后，编程设计、计算并绘制出了铜/钼金属爆炸可焊性窗口。用计算结果指导试验参数的选择，调整选取最佳爆炸复合工艺参数。经实验验证表明，该铜/钼可焊性窗口和最佳工艺参数具有较高的可靠性，此方法用于研究爆炸焊接工艺，能够减少重复实验，提高研究效率。
动能弹侵彻机场跑道的动力学响应与倾角、着速、攻角和靶体材料等因素密切相关。为研究这些因素对动能弹侵彻响应的影响，本文采用正交试验方法设计了数值试验和4 参数3 水平的正交表，进行了侵彻过程中的参数敏感性分析。经过对参数的研究，得到了影响动能弹侵彻动力学响应的主要因素，研究表明：着速和倾角对弹体减速度影响较大，影响剩余速度的主要因素是着速、倾角和基层强度，倾角和攻角对弹体的偏转影响较大。
利用320mm 平衡炮，采用次口径加载技术，将50kg 弹体加速到1100m/s .设计了适合高速侵彻的金属弹体，采用1.5m ×1.5m ×6m 钢筋混凝土作为靶标，进行高速侵彻试验.试验结果表明：次口径弹托与弹丸完全分离，弹体飞行姿态稳定，弹体飞行攻角小于2 °.侵彻后回收的金属弹体结构完整，无明显弯曲，仅弹体头部存在一定塑性变形，弹体质量损失约1.2％，弹体侵蚀现象不明显.
基于我们近年来对金属空心球材料及圆胞蜂窝材料的实验研究、数值模拟和理论分析，本文将对这两类多胞材料的动态力学行为进行总结与比较。由于基体材料行为的应变率相关性，材料的惯性，以及应力波效应，这些多胞材料在动态变形中不但承载能力发生动力增强，而且显现非均匀的变形模式和变形局部化的现象，在冲击速度很高市还出现变形以冲击波的形式传播。
采用基于冲击波衰减动力学的实验技术，利用电炮作为一维平面加载源，驱动Mylar 膜飞片碰撞不同厚度的Zr51Ti5Ni10Cu25Al9 块体非晶合金样品，产生三角形冲击脉冲。采用速度干涉(VISAR )测试技术记录样品自由表面历史，获得了不同的层裂信号。通过对层裂的信号的分析，获得了Zr51Ti5Ni10Cu25Al9 块体非晶合金在拉伸应变率(～105s-1 )下的层裂强度(～5Gpa )，对
为研究柱形装药水中爆炸气泡的坍塌机制和射流特性，采用水中爆炸实验和数值模拟相结合的手段研究了典型的柱形装药自由场水中爆炸气泡运动过程。对比实验和数值模拟结果发现，气泡坍塌所形成的射流与以往的球形气泡坍塌射流行为存在差异。研究结果表明：非球形气泡的坍塌过程与气泡初始形态的局部曲率半径有关，初始气泡表面曲率半径大的区域的初始膨胀速度较曲率半径小的区域的小，坍塌时刻也较早，导致气泡产生非球形形式的坍塌。
本文用矩形可变形平板来简化蝙蝠翼，并将蝙蝠拍翼分解为拍动和变形两部分，其中变形简化为弦向弯曲变形(cambering)、展向弯曲变形(bending)以及扭转变形(twisting)三种基元变形，然后利用面元法求解了基元变形下三维大变形翼的非定常气动力，由此探讨了三种变形下的附加惯性效应、尾涡及侧缘涡对气动力的影响。通过对比计算结果，可以发现：1)附加惯性效应产生高的气动力的瞬时峰值，但气动力的平
本文将混凝土材料看成水泥砂浆基体、粗骨料颗粒和微裂纹组成的三相复合材料，假设水泥砂浆基体为弹性的、均匀的、各向同性的，粗骨料颗粒为均匀的、各向同性的，微裂纹为钱币状。基于Mori-Tanaka 理论和Eshelby等效夹杂理论，建立了混凝土材料在冲击荷载作用下动态本构模型，并给出了微裂纹损伤演化方程。通过模拟曲线与实验曲线进行比较，发现两者拟合良好，因而可以用该模型模拟混凝土材料在冲击荷载下的动态
对纵向导波在水泥砂浆锚固锚杆中的传播特性进行了分析。根据弹性动力学理论,首先建立了三维柱坐标情况下,水泥砂浆锚固锚杆中纵向导波的频散方程,并数值求解了纵向导波的相速度、能量速度和衰减频散曲线。由频散曲线可知,500kHz 范围内水泥砂浆锚固锚杆中存在5 个纵向导波模态；然后分析了特定频率下,纵向导波在水泥砂浆锚固锚杆横截面的位移,结果表明：锚杆截面上轴向位移越大,代表纵向导波对锚杆体的周向缺陷越敏
针对三维流体弹塑性计算程序MMIC-3D 计算网格为直六面体网格的特点，开发出了三维前处理程序MESH-3D .MESH-3D 采用CSG 和STL 模型两种建模方式进行复杂实体建模，并采用基于边的整体切片法的剖分算法，借鉴计算机图形学中的扫描线填充算法完成有限差分网格剖分.在绘制网格时，将网格单元的公共面予以删除，减少大量计算时间和存储空间，实现对网格的快速消隐绘制.MESH-3D 可实现亿级网