超高速碰撞成坑特性分子动力学模拟
|
摘要
基于开源分子动力学程序LAMMPS,对直径为4.86nm的球形铝弹丸以10km/s超高速撞击半无限厚铝靶进行模拟.弹坑形成的物理过程与超高速碰撞宏观现象相似,弹坑深度与宏观经验公式计算结果基本一致,获得了弹丸头部相对于碰撞点的位移随时间的变化规律;分析了靶板中冲击波传播特性,碰撞初期冲击波阵面传播速度达到12km/s,随后冲击波传播速度逐渐减小,接近于弹性波速;弹坑周围观测区发生熔化相变,熔化时间持续0.07ps,熔化层厚度为2.9nm;弹坑周围观测区冷却速率达到1015 K/s量级,抑制了原子重结晶,最终呈现为固相非晶结构.
Based on the open source molecular dynamics program LAMMPS,hypervelocity impact of a spherical aluminum projectile with a diameter of 4.86 nm on the semi-infinite thick aluminum target at the speed of 10 km/s was simulated.The physical process of crater formation was similar to that from a macroscopic impact.The crater depth was consistent with the data from the macro empirical formula and its variation with time was obtained.The propagation characteristics of the shock wave in the target were analyzed.The propagation speed of the shock wave front reached up to 12 km/s during the initial impact,and then the speed gradually decreased and was close to the elastic wave velocity.The observation region around the crater experienced melting state which lasted about 0.07 ps and had a layer thickness of 29 nm.The cooling rate of the observation region around the crater reached the order of 1015 K/s which prevented atom recrystallization.The observation region finally turned into a solid phase amorphous structure.
Based on the open source molecular dynamics program LAMMPS,hypervelocity impact of a spherical aluminum projectile with a diameter of 4.86 nm on the semi-infinite thick aluminum target at the speed of 10 km/s was simulated.The physical process of crater formation was similar to that from a macroscopic impact.The crater depth was consistent with the data from the macro empirical formula and its variation with time was obtained.The propagation characteristics of the shock wave in the target were analyzed.The propagation speed of the shock wave front reached up to 12 km/s during the initial impact,and then the speed gradually decreased and was close to the elastic wave velocity.The observation region around the crater experienced melting state which lasted about 0.07 ps and had a layer thickness of 29 nm.The cooling rate of the observation region around the crater reached the order of 1015 K/s which prevented atom recrystallization.The observation region finally turned into a solid phase amorphous structure.
引文
[1]张庆明,黄风雷.超高速碰撞动力学引论[M].北京:科学出版社,2000.Zhang Qingming,Huang Fenglei.Introduction on hypervelocity impact dynamics[M].Beijing:Science Press,2000.(in Chinese)
[2]于超,李平.钨合金弹丸超高速撞击的分子动力学研究[J].北京理工大学学报,2014,34(增刊1):64-66.Yu Chao,Li Ping.Molecular dynamic investigation of hypervelocity impact by a tungsten alloy projectile[J].Transactions of Beijing Institute of Technology,2014,34(suppl 1):64-66.(in Chinese)
[3]Anders C,Bringa E M,Urbassek H M.Crater production by energetic nanoparticle impact on Au nanofoams[J].Applied Physics Letters,2016,108:4417.
[4]巨圆圆,张庆明,龚良飞,等.球形弹丸超高速撞击铝靶的分子动力学模拟[J].航天器环境工程,2018,35(2):153-157.Ju Yuanyuan,Zhang Qingming,Gong Liangfei,et al.Molecular dynamics simulation for hypervelocity impact of spherical projectile to aluminum target[J].Spacecraft Environment Engineering,2018,35(2):153-157.(in Chinese)
[5]Mei J,Davenport J W,Fernando G W.Analytic embedded-atom potentials for fcc metals:application to liquid and solid copper[J].Physical Review B,1991,43:4653.
[6]Ju Y Y,Zhang Q M,Gong Z Z,et al.Molecular dynamics simulation of shock melting of aluminum single crystal[J].J Appl Phys,2013,114:093507.
[7]Branicio P S,Kalia R K,Nakano A,et al.Atomistic damage mechanisms during hypervelocity projectile impact on AlN:a large-scale parallel molecular dynamics simulation study[J].Journal of the Mechanics and Physics of Solids,2008,56:1955-1988.
[8]Frost V C.Meteoroid damage assessment[J].NASA Spec Publ,1970,8042:30.
[9]冉宪文,汤文辉,谭华,等.含电子影响的冲击温度计算[J].高压物理学报,2006,20:153-156.Ran Xianwen,Tang Wenhui,Tan Hua,et al.Calculation of shock temperature by considering electronic contribution[J].Chinese Journal of High Pressure Physics,2006,20:153-156.(in Chinese)
[10]Saiz F,Gamero-Casta1o M.Molecular dynamics of nanodroplet impact:the effect of the projectiles molecular mass on sputtering[J].AIP Advances,2016,6:065319.
[2]于超,李平.钨合金弹丸超高速撞击的分子动力学研究[J].北京理工大学学报,2014,34(增刊1):64-66.Yu Chao,Li Ping.Molecular dynamic investigation of hypervelocity impact by a tungsten alloy projectile[J].Transactions of Beijing Institute of Technology,2014,34(suppl 1):64-66.(in Chinese)
[3]Anders C,Bringa E M,Urbassek H M.Crater production by energetic nanoparticle impact on Au nanofoams[J].Applied Physics Letters,2016,108:4417.
[4]巨圆圆,张庆明,龚良飞,等.球形弹丸超高速撞击铝靶的分子动力学模拟[J].航天器环境工程,2018,35(2):153-157.Ju Yuanyuan,Zhang Qingming,Gong Liangfei,et al.Molecular dynamics simulation for hypervelocity impact of spherical projectile to aluminum target[J].Spacecraft Environment Engineering,2018,35(2):153-157.(in Chinese)
[5]Mei J,Davenport J W,Fernando G W.Analytic embedded-atom potentials for fcc metals:application to liquid and solid copper[J].Physical Review B,1991,43:4653.
[6]Ju Y Y,Zhang Q M,Gong Z Z,et al.Molecular dynamics simulation of shock melting of aluminum single crystal[J].J Appl Phys,2013,114:093507.
[7]Branicio P S,Kalia R K,Nakano A,et al.Atomistic damage mechanisms during hypervelocity projectile impact on AlN:a large-scale parallel molecular dynamics simulation study[J].Journal of the Mechanics and Physics of Solids,2008,56:1955-1988.
[8]Frost V C.Meteoroid damage assessment[J].NASA Spec Publ,1970,8042:30.
[9]冉宪文,汤文辉,谭华,等.含电子影响的冲击温度计算[J].高压物理学报,2006,20:153-156.Ran Xianwen,Tang Wenhui,Tan Hua,et al.Calculation of shock temperature by considering electronic contribution[J].Chinese Journal of High Pressure Physics,2006,20:153-156.(in Chinese)
[10]Saiz F,Gamero-Casta1o M.Molecular dynamics of nanodroplet impact:the effect of the projectiles molecular mass on sputtering[J].AIP Advances,2016,6:065319.