三角波强加载下延性金属多次层裂破坏问题

基于应力瞬时断裂判据和Tuler-Butcher损伤累积判据，分析了没有升压的简单三角波强加载下延性金属的多次层裂破坏问题。分析结果显示:层裂片的厚度随着冲击波宽度与强度的比值的增大而增大，材料的破坏深度小于冲击波宽度的一半；在强冲击加载下，材料的破坏深度约为冲击波宽度的一半。     

三维光滑粒子流体动力学并行计算程序

介绍了基于消息传递并行程序设计平台研制的三维光滑粒子流体动力学并行程序CSPH3D.包括计算格式、并行方案、并行程序逻辑,以及加快邻域粒子搜索的处理方法.对三维微喷射和斜侵彻的计算表明:CSPH3D程序可以较好地计算这类问题.并且程序具有较高的并行效率.对于粒子总数为1527402的微喷射算例和粒子总数为1454225的斜侵彻算例,使用100个处理器时,并行效率可以达到80%.     

圆柱形含能材料热爆炸临界参数计算

根据变分原理,计算了各表面散热条件不同的有限长圆柱含能材料的热爆炸临界参数和临界温度,研究了长径比对热爆炸临界参数和临界温度的影响;计算了各表面散热条件不同的有限长圆柱材料的热爆炸临界参数,得出了Biot数和对热爆炸临界参数和临界温度的作用。计算表明,用变分的方法将求解热爆炸导热方程问题转化为对一个特征值方程求解,求解过程要简单得多,且精度高、速度快。     

Molecular dynamics simulations of jet breakup and ejecta production from a grooved Cu surface under shock loading

Large-scale non-equilibrium molecular dynamics simulations are performed to explore the jet breakup and ejecta production of single crystal Cu with a triangular grooved surface defect under shock loading. The morphology of the jet breakup and ejecta formation is obtained where the ejecta clusters remain spherical after a long simulation time. The effects of shock strength as well as groove size on the steady size distribution of ejecta clusters are investigated. It is shown that the size distribution of ejecta exhibits a scaling power law independent of the simulated shock strengths and groove sizes. This distribution, which has been observed in many fragmentation processes, can be well described by percolation theory.     

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

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

缺陷形状及角度对于铝材料微喷射的影响

目前，一般认为引起微喷射的主要物理机制有：(1)试件的表面缺陷(如凹坑、划痕等)形成的微射流机制，很可能是微喷射的主要来源。(2)在冲击波作用下，晶界处产生局部能量沉积，而形成“热点”，致使物质局部熔化甚至汽化。(3)材料表面的杂质可能产生微层裂或剥落而形成微喷射。在材料发生微喷射的过程中，这些机理经常同时存在，对于不同的材料，在不同的条件下，可能是不同的机制起主导作用。在实验中发现，铝试件在冲击压为40GPa时，其表面温度大约为500K，这大大低于931K的熔化温度。因此，可以认为铝材料在冲击压力低于40GPa时，没有熔化现象发生。故而对于纯铝试件，当冲击压力达到约30GPa时，铝材料表面缺陷形成的微射流是其微喷射的主要机制。     

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

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

金属射流发散性理论分析及数值模拟

分析经典射流理论和相关文献,给出了在汇聚点坐标系中金属Cu对碰形成射流的汇聚射流区、无射流激波封锁区、无射流强度封锁区和形成发散射流区域。对金属超声速射流形成的发散性问题进行了理论分析,诠释了Walker基于实验提出射流发散理论,证明了金属射流超声速部分可形成发散射流也可形成汇聚射流,且超过1.23倍体声速的金属射流必定是发散的。最后,应用自编的欧拉计算程序MEPH对金属Cu以不同速度、倾角对碰射流形成过程进行数值模拟,得到了分叉射流、空洞射流和密实稀疏射流等的典型射流发散模式图像,印证了理论分析的结果。     

微喷混合问题的初步研究

基于颗粒轨道模型,采用二元碰撞假设,在二维拉氏流体动力学程序上实现了气-粒两相混合程序的研制.使用气-粒两相混合程序对不同充气压力条件下的微喷混合问题进行了数值模拟,给出了喷射场随时间演化过程,计算结果与实验定性一致.     

一维应变加载下单晶铁结构转变的微观研究

采用嵌入原子势和分子动力学方法,模拟了单晶铁在一维应变条件下由体心立方（bcc）转变为六角密排（hcp）结构的微观过程. 当应变加载至相变临界值时,hcp相开始均匀形核并沿{011}晶面长大为薄片状体系.弹性常数C₃₁和C₃₂在相变前被逐渐硬化,C₃₃则在相变前出现软化行为;当体系完全相变后,上述各弹性常数显示开始随体积压缩而迅速硬化,温度效应对晶格具有软化作用,可削弱C₃₃的硬化和软化过程;样品在压缩过程可出现孪晶结构,孪晶结构使晶格发生剪切变形.混合相中,hcp相势能比bcc相高,最大剪应力方向与bcc相反向;系统的偏应力与hcp相质量分数近似呈线性关系. 关键词： 结构转变 分子动力学 一维应变