超高速撞击下碎片云的OTM分析  

OTM analysis of debris cloud under hypervelocity impact

作　　者：廖祜明 黎波 樊江[1] 焦立新 于帅超 林健宇 裴晓阳[4] LIAO Huming;LI Bo;FAN Jiang;JIAO Lixin;YU Shuaichao;LIN Jianyu;PEI Xiaoyang(School of Energy and Power Engineering,Beihang University,Beijng 102206,China;ESCAAS Co.Ltd.,Beijing 100027,China;College of Engineering,Peking University,Beijing 100871,China;Institute of Fluid Physics,China Academy of Engineering Physics,Sichuan Mianyang 621900,China)

机构地区：[1]北京航空航天大学能源与动力工程学院,北京102206 [2]云翼超算(北京)软件科技有限公司,北京100027 [3]北京大学工学院,北京100871 [4]中国工程物理研究院流体物理研究所,四川绵阳621900

出　　处：《爆炸与冲击》2022年第10期48-58,共11页Explosion and Shock Waves

摘　　要：空间碎片超高速撞击是典型的高温、高压、高应变率的极限力学问题,涉及材料复杂的动态响应,对传统的数值方法提出了巨大挑战。最优运输无网格(OTM)方法通过有机结合最优运输时间积分理论、局部最大熵无网格近似、物质点抽样、基于物理的裂纹扩展算法以及大规模并行计算策略,克服了传统数值方法瓶颈,在理论上保证了不同形式能量耗散的自主耦合分配,为超高速撞击仿真预测提供了高效的解决方案。采用基于OTM方法自主研发的极限力学仿真软件ESCAAS,对不同质量(3、10 g)的铜飞片以不同撞击角度(5.4°、11.7°)和不同撞击速度(5.55、5.12 km/s)撞击铝合金靶板的过程进行数值模拟,获得碎片云的形貌、靶板穿孔孔径等结果,与实验测量数据吻合良好,显示出OTM方法及ESCAAS软件可以作为超高速撞击的有力数值分析手段。The hypervelocity impact(HVI)of space debris is a typical extreme mechanics problem at high temperature,high pressure and high strain rate.The HVI involves the complex dynamic response of materials.Numerical methods have become very useful and important tools to predict the complex phenomena and look into the details of the entire process.However,the accurate simulation of the HVI is a grand challenge in scientific computing that places exacting demands on physics models,numerical solvers and computing resources.The optimal transportation meshfree(OTM)method is a meshfree updated-Lagrangian methodology for fluid and solid dynamic flows,possibly involving multiple phases,viscosity and general equations of state,general inelastic and history-dependent constitutive relations,arbitrary variable domains and boundary conditions and the interaction between fluid flows and highly deformable structures.The rationale behind the approach is combining concepts from the optimal transportation(OT)theory with material-point sampling,local maximum entropy(LME)approximation,the seizing contact,variational material point failure algorithm,and overcomes the essential difficulties in grid-based numerical methods like Lagrangian and Eulerian finite element method.Owing to those advantages,the OTM method provides an efficient and accurate solution for HVI simulation.In this paper,large scale three-dimensional numerical simulations on the HVI of the copper projectiles with different thicknesses(3.45,5.13 mm),different masses(3,10 g),different impact angles(5.4°,11.7°)and different impact velocities(5.55,5.12 km/s)impacting the Al6061-T6 plate with different thicknesses(2.87,4.39 mm)were performed within the software ESCAAS based on the OTM method using a dynamic load balancing MPI/Pthreads parallel implementation.The dynamic response of material including phase transition in the high strain rate,high pressure and high temperature regime expected in this paper was described by the use of a variational thermomechanical coupling constit

关 键 词：超高速撞击 OTM方法 裂纹扩展 无网格方法 

分 类 号：O385[理学—流体力学]