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

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

多组份计算中封闭性模型研究

在使用单速度多组份方法计算多介质混合网格中的物理量时,需要给出一种封闭模型来使控制方程封闭。本文分析了压力增量相等封闭性模型存在基础遭到破坏的原因,之后在保证流场守恒性的前提下构造可以保证假设基础的可解方程组对压力增量相等封闭性模型进行修正,给出了一种改进方法。使用改进的模型进行了爆轰驱动飞片问题的计算,计算结果表明,改进的模型能很好的追踪物质界面和冲击波位置,同时较好的抑制界面处压力、内能震荡。     

单调非线性互补问题的不可行内点算法

为了克服内点算法初始点不易给出的缺陷，本文给出了一个求解单调非线性互补问题的不可行内点算法，并证明了算法的收敛性。     

材料特性对爆轰驱动飞层对碰凸起影响的数值模拟

应用多介质弹塑性L-R两步欧拉计算方法,对炸药两端起爆情况下柱面爆轰驱动飞层对碰凸起现象进行了数值模拟。在飞层为铅材料情况下,给出了铅飞层熔化区域及断裂图像,将对碰凸起计算结果与实验X光图像进行比较的基础上,进一步探索凸起现象的形成机制。从材料的熔点、弹塑性、密度、声速等方面对金属飞层对碰区凸起的影响进行了分析。综合分析几个计算模型的数值模拟结果,初步认为:材料强度大、密度大、熔点高能使飞层对碰凸起钝化,而材料声速较小会造成对碰区波系关系发生变化且正压作用时间增加,从而形成较严重的对碰凸起,材料声速大小对对碰凸起的影响较强度、密度等因素更明显。     

炸药爆轰产物Jones-Wilkins-Lee状态方程不确定参数

Jones-Wilkins-Lee (JWL)状态方程是一种不显含化学反应、由实验方法确定参数的半经验状态方程, 能比较精确地描述爆轰产物的膨胀驱动做功过程. 在JWL状态方程中有多个未知(不确定)参数需要确定. 传统的确定JWL状态方程参数的方法是“调参数”, 人为因素影响较大, 无法给出参数的不确定性信息. 本文利用贝叶斯分析方法研究了炸药的不确定参数, 该方法能够基于以往的认识、实验和模拟数据标定(calibration)不确定参数. 在本文结果中, 不确定参数的后验分布均值与文献结果相符合, 基于参数标定结果的数值模拟90%置信区间完全包含实验数据. 数值标定结果说明贝叶斯参数标定适用于确定样品炸药的JWL状态方程参数. 特别是, 在本文JWL状态方程参数标定过程中极大地减少了人为因素的影响.     

钝感炸药点火增长模型的欧拉数值模拟

在自主研发的二维多介质欧拉弹塑性流体力学程序中,通过引入点火增长的反应率模型以及炸药 减敏模型,借助网格自适应技术,研究钝感炸药的冲击点火、直径效应以及死区形成等爆轰现象。数值模拟结 果表明,该程序能够正确模拟平面爆轰波的爆速、CJ状态、vonNeumann尖点等爆轰参数;并能够较好模拟炸 药的直径效应。另外,通过引入考虑减敏效应的反应率模型,能较好地模拟钝感炸药的死区形成过程。     

碎片云问题的三维欧拉数值模拟

针对材料超高速碰撞产生碎片云问题进行了多介质弹塑性流体力学欧拉方法数值模拟,介绍了所采用的数学模型、材料模型及相应的计算方法,模拟了铝球超高速碰撞铝板和斜铜柱高速撞击铝板形成碎片云的实验现象。进行了不同碰撞速度下模拟的碎片云形状和实验结果比较,利用经验公式对数值模拟的板靶开孔尺寸进行了分析。数值结果与实验结果、经验公式结果比较,说明多介质弹塑性流体力学欧拉方法能够较好地模拟超高速碰撞产生的碎片云问题,验证了数值方法的有效性。     

求解凸二次规划问题的不可行内点算法

该文对一般的凸二次规划问题，给出了一个不可行内点算法，并证明了该算法经过犗（狀２犔）步迭代之后，要么得到问题的一个近似最优解，要么说明该问题在某个较大的区域内无解．     

线性约束凸规划问题的仿尺度算法

对一般线性约束凸规划问题，给出了相应的仿尺度算法，并证明由该算法所产生的迭代点要么直接到达问题的最优解，要么其极限点满足问题的最优性条件。     

含氦泡辐照老化材料层裂损伤计算方法分析

辐照条件下，一些材料内部产生大量的氦泡等微缺陷，氦泡的大小和数密度随着辐照年限的增长而增长。氦泡分布特征的变化不仅影响材料本身的物理、力学性质，而且直接影响材料层裂损伤演化后期材料破坏颗粒度的分布特征。延性材料的层裂损伤演化过程一般包括孔洞的成核、增长和汇合，但因已有孔洞对新成核孔洞存在抑制作用，当初始孔洞数密度达到一定临界值时，材料内部没有新的孔洞成核，因此，层裂损伤的计算可以不考虑新孔洞成核的影响。本文中基于损伤早期演化的特征，给出了这一临界值的计算方法，并进一步探讨了含氦泡辐照老化钚材料层裂损伤的计算方法。同时，在完善孔洞增长(void growth, VG)层裂损伤模型中参数的确定方法的基础上，借助含氦泡常规铝材料的层裂实验结果，对此问题进行了定性的分析：在氦泡尺寸变化不大的情况下，当氦泡浓度低于临界氦泡浓度时，需要考虑初始氦泡以及新增孔洞的综合影响；反之，可以采用简单的层裂损伤模型，不需要计算孔洞成核，但由于增长孔洞之间的相互影响，损伤模型的初始损伤参数需要重新确定。