爆炸与冲击动力学若干问题研究进展

介绍了爆炸与冲击动力学数值模拟领域近几年的研究重点及主要进展.包括材料本构关系、损伤破坏、冲击相变、射流和微喷等,其中材料本构关系的发展比较系统,有几个典型的模型可以参考,不仅应用广泛,而且有实验数据进行检验和对比.损伤模型只有唯象模型,虽然有不同的物理机制,但缺少直接的实验依据.射流的工程应用比较广泛,理论研究比较缺乏.微喷研究基本停留在定性阶段.随着实验技术及测量手段的发展,相关领域取得了一些令人欣慰的进展.     

爆轰波对碰驱动平面锡飞层的动力学及动载行为特征研究

强冲击下的物质变形、破坏及诱发的轻重介质混合问题,是内爆压缩科学和工程应用领域的研究重点.本文针对爆轰波对碰条件下的复杂加载动力学过程及其动载破坏形态特征,开展数值模拟研究与极曲线理论分析.设计了爆轰波对碰驱动平面锡飞层的计算模型,获得了爆轰加载动力学过程及波系相互作用物理图像,分析了锡飞层对碰区自由表面速度历史的典型特征.给出了锡飞层中折射激波对碰发生马赫反射的临界条件,解读了三波结构的传播行为,阐明了对碰区内存在"一维正冲击"区域,一维区外存在单次斜冲击向两次斜冲击过渡的复杂加载动力学过程,提出了对碰区冲击动力学模型,揭示了影响对碰区动载行为特征的机理.数值模拟结果与极曲线理论分析结果相互印证,符合较好.本文的研究成果,将为深入理解和解读对碰区特殊的物质破坏及混合现象提供重要的理论支撑.     

二维拉格朗日坐标系下气粒混合双向耦合对激波流场影响的计算

喷射颗粒与气体混合是内爆压缩领域的热点和难点. 针对喷射混合中的气粒双向耦合问题, 开展了理论建模、离散算法以及颗粒反馈对激波流场的影响研究. 建立了拉格朗日计算框架下的数学模型; 给出了耦合源项的离散算法; 开展了平面及汇聚构型条件下, 气粒双向耦合的数值模拟研究; 发现了颗粒反馈导致气体激波提速现象以及气区流场物理量分布形态的改变, 初步获得了量化分析结果. 本文建立的数学模型、计算方法和获得的新的物理认识, 为深入理解喷射混合现象、解决相关工程应用问题提供了重要理论支撑.     

FP-1装置铝套筒内爆动力学过程的一维磁流体力学模拟

作为一种重要的柱面会聚冲击和准等熵压缩加载源,磁驱动固体套筒内爆技术已广泛应用于高能量密度物理实验研究.针对FP-1装置驱动的固体套筒内爆动力学过程,建立了含强度的一维磁流体力学模型,并对典型实验进行了模拟.计算获得的套筒内爆速度同实验结果较为相符.模拟结果显示,该装置在40 kV充压条件下,可以将直径3 cm,厚0.5 mm的铝套筒加速至1.1 km/s,内壁速度超过1.5 km/s,同时保持大部分材料为固体状态.内爆套筒与相同材料靶筒碰撞产生的冲击压力约9 GPa.改变靶筒内部填充气体的压力,可以获得不同的靶筒运动速度、轨迹以及反弹半径,以满足不同类型实验的研究需要.     

Ejecta from periodic grooved Sn surface under unsupported shocks

Dynamic failure and ejection characteristics of a periodic grooved Sn surface under unsupported shock loading are studied using a smoothed particle hydrodynamics method. An "Eiffel Tower" spatial structure is observed, which is composed of high-speed jet tip, high-density jet slug, longitudinal tensile sparse zone, and complex broken zone between grooves. It is very different from the spike-bubble structure under supported shocks, and has been validated by detonation loading experiments. In comparison with that under supported shocks at the same peak pressure, the high-speed ejecta decreases obviously, whereas the truncated location of ejecta moves towards the interior of the sample and the total mass of ejecta increases due to the vast existence of low-speed broken materials. The shock wave profile determines mainly the total ejection amount, while the variation of V-groove angle will significantly alter the distribution of middle-and high-speed ejecta, and the maximum ejecta velocity has a linear correlation with the groove angle.     

单晶Cu(001)薄膜塑性变形的分子动力学模拟

使用分子动力学方法,模拟研究了单晶Cu(001)薄膜在双向等轴拉伸应变下的塑性变形行为.当应变超过一定值时,样品通过产生位错、层错及孪晶而发生塑性变形.当应变相对较低时,不全位错首先在薄膜表面形核并在密排面上滑移,留下堆积层错;当应变增加时,位错在表面与内部同时成核生长,层错数量也随之增加.分析了相邻滑移面上的位错之间相互作用形成孪晶的微观过程.材料内部形成大量堆积层错及孪晶后,较大孪晶的密排面上的原子也会发生滑移,形成孪晶内部的层错结构以释放残余应力.     

金属Al表面熔化各向异性的分子动力学模拟

采用嵌入原子势,使用分子动力学方法对金属Al不同低指数晶面的表面熔化现象分别进行了模拟.分析了熔化过程中样品结构组态的变化.模拟结果表明对于不同的自由表面,表面熔化呈现出明显的各向异性行为.Al（110）面在低于熔点的温度之下发生预熔化;（111）与（001）面都出现过热现象.与（111）面不同,（001）面发生过热现象时表面原子层为类液层,而（111）面仍然保持很好的晶格结构,即预熔化的Al（001）面在高于熔点的温度下,仍可以在很长的时间内处于相对稳定的亚稳态.由模拟得到Al的热力学熔点为950 K左右,与实验值基本符合. 关键词： 分子动力学 表面熔化 过热     

强动载作用下孔洞汇合对延性金属层裂损伤演化过程的影响

针对强动载作用下延性金属的层裂问题,在分析孔洞之间几何关联的基础上,定义了一个新的耦合损伤及孔洞几何信息的孔洞汇合判定方法,同时,基于能量守恒原理,解析了孔洞汇合对损伤快速增长影响的物理机理.通过分析数值计算结果和对比相关文献的实验可知:孔洞汇合后不仅引起损伤增长,而且导致了损伤材料内部微孔洞数目的减少、孔洞平均尺寸的增加。     

一种新的光子多普勒速度频谱分析方法

光子多普勒速度计可给出飞层表面某一速度带内颗粒群速度随时间演化的频谱数据, 在冲击动力学实验尤其是微喷射及其混合研究中得到广泛应用. 本文提出一种新的光子多普勒频谱数据分析方法, 可推断出混合区厚度变化和前端等效颗粒尺度. 利用该方法, 对一些典型状态下喷射混合速度频谱开展分析, 获得了不同冲击压力、气体条件下颗粒度数据, 证实了气体环境下喷射颗粒的气动破碎现象, 以及破碎后尺度与初始条件的依赖性, 为喷射混合物理规律研究提供了重要依据.     

L-R两步欧拉方法在铝材料微喷射中的应用

应用三维弹塑性流体力学Lagrangian-Remapping两步欧拉计算方法对铝材料微喷射现象进行了数值模拟研究。计算了Asay实验中表面刻有相同深度、不同夹角沟槽的金属铝微喷射模型,计算得到的微喷物总质量、最大射流速度和实验结果均符合较好。进一步展开了对相同深度、更大夹角范围沟槽微喷射的数值模拟。分析认为喷射最大速度随沟槽角度的增大呈线性下降趋势。同时给出了喷射系数随沟槽角度的变化的拟合关系曲线,看到由于材料强度及沟槽角度变化后造成的波系关系变化的影响,随着沟槽角度增加,喷射系数曲线呈明显非线性发展。