一种周期库仑作用势优化法的改进

利用简单多项式基组拟合变化平缓的长程作用部分,并增加多个限定条件以提高优化势与库仑势的局部重合度,对Natoli优化法进行改进.在不增加计算复杂度的前提下,得到平均偏差明显降低的优化势;建议优化时划分r空间与k空间的准则取k c r c≥15.在稠密液氢分子模拟中,优化法获得的能量及压力均小于Natoli法与Ewald法的结果.     

修正金属本构模型在超高速撞击模拟中的应用

为更加准确地计算93钨合金弹超高速撞击Q345钢板问题,构建了修正的金属本构模型。引入GRAY三相物态方程描述材料相态变化,采用Johnson-Cook强度模型描述撞击后期材料的力学行为。结合封加波损伤演化模型以及Johnson-Cook失效模型描述不同应力三轴度下材料的拉伸、剪切失效行为;引入曹祥提出的断裂演化模型,描述材料失效后应力归零的过程。通过对比超高速撞击数值模拟结果与实验结果,验证了本构模型的适用性,并进一步分析了典型弹靶撞击条件下破片群的空间分布特征。研究结果表明:基于修正金属本构模型获得的超高速撞击靶板穿孔直径、弹体侵蚀长度、破片群扩展速度结果与实验结果一致;GRAY三相物态方程能够相对准确地给出弹体撞击首层靶板以及剩余弹体、破片群撞击第2层靶板时弹靶材料的熔化情况;封加波损伤演化模型能够准确判断超高速撞击过程中靶板是否产生层裂破坏;综合封加波损伤演化模型、Johnson-Cook失效模型以及曹祥提出的断裂演化模型后,数值模拟获得的破片群撞击后效靶板的穿孔面积与累积数量的统计曲线结果与实验结果一致;获得了典型条件下的柱形93钨弹体超高速撞击Q345靶板破片群空间分布结果,破片群的前端具有较高的质量、轴向动量以及横向动量（绝对值）。     

防雷舱结构在聚能装药水下爆炸作用下的毁伤研究

通过实验以及数值模拟，开展了防雷舱结构在聚能装药水下爆炸作用下的毁伤研究，探讨冲击波在多介质结构中的传播规律及其对结构的毁伤机理。结果表明:相同爆炸当量条件下，爆破型装药仅对空舱产生破坏，防雷舱的多介质复合结构对冲击波具有较强的衰减作用，聚能型装药形成的侵彻体可造成液舱前板及后板的穿孔，孔径约为装药直径的1/3~1/2。加宽液舱以及加厚液舱后板可以增强防雷舱的抗爆能力。     

共振价键波函数在高压液氢量子蒙卡模拟中的适用性研究

在共振价键理论基础上, 选取高压液氢电子主要占据轨道的线性组合作为基组, 构建由Jastrow项和反对称孪生函数乘积项 (AGP) 组成的波函数. 考虑电子关联作用的共振价键 (RVB) 波函数得出的能量值低于LDA能量值; 当满足r_s<1.75或T >15000 K时引入backflow项以改善波函数结点面, 改善后的能量值下降约1 mHa/atom, 能量方差值变小. 将构建的RVB波函数与电子-离子耦合的蒙特卡罗法 (CEIMC) 相结合, 计算结果与实验及其他ab-initio结果相符合, 获得的液氘单次冲击Hugoniot曲线基本通过所有加载类型实验误差棒, 液氘在50.3 GPa处具有最大压缩率4.48, 在100–120 GPa内未发现压缩率有急剧增大的现象. 构建的RVB 波函数能够适用于较宽密度与温度范围内(1.0< r_s<2.2, 2800 K< T<60000 K)液氢的模拟, 与CEIMC法相结合可提高液氢冲击特性的模拟精度. 关键词： 共振价键理论 波函数 量子蒙卡法 液氢     

钨合金弹体对混凝土靶的超高速侵彻机理

为研究钨合金弹体超高速侵彻混凝土靶的相关机理,构建了适用于超高速撞击的金属强度模型、失效模型和混凝土的本构模型,对93钨合金弹体超高速撞击混凝土靶问题进行了数值模拟。开展了钨合金弹体超高速撞击混凝土靶实验,分析了靶板成坑特性,研究了侵彻总深度和残余弹体长度随撞击速度的变化规律,理论分析了长杆钨弹超高速撞击混凝土的侵彻模型和混凝土靶内的应力波传播。得到以下主要结论:（1）利用金属及混凝土的新本构模型获得的超高速撞击混凝土靶的破坏形貌数值模拟结果与实验结果一致;（2）超高速撞击条件下混凝土靶成坑为“弹坑＋弹洞”形,成坑体积与弹体动能近似成正比;（3）超高速撞击条件下,侵彻深度随弹速提高呈现先增大后减小的现象,高速段侵深降低是弹体经历销蚀侵彻后“刚体侵彻阶段”减少造成的;（4）建立的钨合金超高速撞击混凝土侵彻分析模型,可用来预估侵彻深度、残余弹长、蘑菇头直径等参数;（5）采用建立的超高速撞击混凝土靶内应力波传播理论模型得到的计算结果与实验结果吻合较好。     

杆式钨合金弹超高速撞击薄靶的能量损耗

超高速撞击过程伴随着复杂的物理过程。为分析杆式圆柱形钨合金弹超高速撞击薄钢靶时的物理过程,采用AUTODYN/SPH数值仿真计算方法获得了撞击过程模型及每个光滑粒子流体动力学信息,并通过广度搜索破片识别程序识别每个破片所含粒子,利用MATLAB编程对破片粒子数据信息进行统计分析,获得弹靶撞击过程的变化特性、弹靶破片数量、相关能量随撞击时间的变化规律。通过分析发现:随着弹体撞击速度的增加,剩余弹体被严重侵蚀,且弹体能量损耗增加,弹体损失的能量主要转变为弹靶破片动能;计算得到了撞击20μs时的能量损耗直方图,同时分析了发生撞击时靶板的能量变化过程,并简要描述了该过程。     

液氘单次冲击压缩的QMC模拟研究

液氢在较宽压力与温度范围内具有复杂的物理特性,须从分子层面构建精确模型开展探索研究。利用电子-核耦合的CEIMC法模拟液氘单次冲击实验,分析了液氘冲击特性,当压力达50.3GPa时液氘具有最大压缩率4.48,在110GPa冲击压力附近未发现有压缩率急剧增大的迹象。选用合适的Al基板材料模型,建立了液氘单次压缩状态与实验条件间的关系,总结了单次冲击实验规律。得到的状态方程与现有动高压实验结果一致,也与经修正后的100GPa以上压力的Omega激光实验值吻合,说明采用基于共振价键理论的波函数后,CEIMC法可应用于液氘的冲击模拟。     

93W杆式弹超高速撞击多层Q345钢靶毁伤及微观分析

为了解杆式弹超高速撞击多层薄钢靶的破坏过程及毁伤机理,开展了克级93W杆式弹正撞击多层Q345钢靶实验及数值模拟研究,通过扫描电子显微镜（scanning electron microscope,SEM）及金相显微镜,分析了超高速撞击实验后靶板材料的微观组织及成分。结果表明,超高速撞击作用下,靶板呈现出“翻唇”穿孔变形、花瓣状塑性变形、撕裂、撞击成坑及鼓包等破坏模式。靶板前3层毁伤以超高速穿孔为主,孔洞数目多但面积小,后几层靶板毁伤孔洞数目少且孔径呈先增大后减小趋势。微观分析表明靶材在强冲击压力下发生晶粒碎化、熔化及再结晶,撞击过程中会形成微孔聚集与微裂纹,可见靶板失效主要是熔融混合物冷却过程中产生的热应力与切应力下的剪切撕裂综合作用的结果。     

薄壁弹体高速侵彻钢筋混凝土实验研究

设计了一种优化结构的薄壁高速侵彻弹体,采用YOUNG方程预估弹体高速侵彻钢筋混凝土的侵彻深度,通过建立三维有限元模型分析并计算了薄壁弹体高速侵彻2.4m和2.8m钢筋混凝土靶的能力和结构变形情况。在100mm口径的火炮上,开展了8kg弹体高速侵彻钢筋混凝土靶实验,考核了弹体的结构强度和侵彻能力。结果表明:薄壁弹体高速侵彻钢筋混凝土靶板后结构完整,在1020m/s和1200m/s速度下具备穿透2.4m和2.8m钢筋混凝土靶的能力。