超高速动能撞击密实小行星的动量传递规律.pdf

超高速动能撞击密实小行星的动量传递规律刘文近1，张庆明1，龙仁荣1，龚自正2，任健康1，任思远2，武 强2，宋光明2，陈 川2，张品亮2（1.北京理工大学 爆炸科学与技术国家重点试验室，北京 100081；2.北京卫星环境工程研究所，北京 100094）摘 要：动能撞击被认为是使潜在威胁小行星偏离与地球碰撞轨道的有效方法。基于撞击速度对动量传递系数的影响规律，利用二级轻气炮进行了6 mm铝弹丸以24 km/s超高速撞击玄武岩靶试验，利用AUTODYN光滑粒子流体动力学方法开展了数值模拟，将模拟结果和试验结果进行了比较，验证了计算模型和参数的准确性。模拟结果表明在不同撞击速度下抛射物的质量和速度分布基本一致，撞击速度越大累计抛射物质量越大。结合试验和数值模拟结果，获得了动能撞击密实小行星的动量传递系数相似律，密实小行星的动量增强系数随撞击速度的0.65次方的增大而增大。动量传递相似律可为动能撞击偏转小行星提供数据支撑。关键词：超高速碰撞；相似律；行星防御；数值模拟中图分类号：P185.7；O383；V412 文献标识码：A 文章编号：2096-9287（2023）04-0420-08DOI:10.15982/j.issn.2096-9287.2023.20230042引用格式：刘文近，张庆明，龙仁荣，等.超高速动能撞击密实小行星的动量传递规律J.深空探测学报（中英文），2023，10（4）：420-427.Reference format:LIU W J，ZHANG Q M，LONG R R，et al.Momentum transfer law of hypervelocity kineticimpacting dense asteroidsJ.Journal of Deep Space Exploration，2023，10（4）：420-427.引言虽然小行星和彗星撞击地球的概率低，但后果和危害极其严重1。动能撞击目前被认为是使对地球有潜在威胁小行星偏离与地球碰撞轨道的最有效方法2。2021年11月24日美国成功发射了“双小行星重定向测试”（Double Asteroid Redirection Test，DART）探测器，于2022年9月26日以超过6.2 km/s的速度撞击双小行星Didymos的次星Dimorphos，成功进行了人类历史上第1次全尺寸的小行星动能撞击偏转测试试验3。小行星是一个多样化的种群，在形状、密度、强度和结构方面有很大差异4。当通过动能撞击偏转技术改变小行星轨道时，小行星物质属性通常是未知的5。因此，需要通过理论分析、试验和数值模拟，分析撞击速度和小行星的物理特性对动量传递系数的影响规律，降低动能撞击偏转小行星的风险和不确定性。Gault等6通过超高速撞击试验获得了6.3 km/s撞击玄武岩的动量传递系数。Hoerth等7撞击不同孔隙的岩石结果表明，动量传递系数随着岩石孔隙率的增加而降低。Chourey等8对多孔干砂、玻璃微球、月球风化层模拟物的撞击试验结果表明动量传递系数均随着撞击速度的增加而增加。George等9在美国国家航空航天局（National Aeronautics Space and Administration，NASA）艾姆斯垂直靶场对各种陨石和陨石模拟物进行了一系列超高速碰撞试验，结果表明含水陨石模拟物动量传递系数均显著大于同等孔隙度无水陨石的动量传递系数，动量传递系数随着陨石孔隙率的增加而降低。受试验条件的限制，试验速度大都低于6.5 km/s，而且大部分超高速撞击试验缺乏对小行星模拟材料性能的完整描述。试验室模型试验通常不能直接重现动能撞击小行星表面所涉及的大撞击器尺寸和极高撞击速度10。数值模拟可有效弥补模型试验在尺寸和速度的不足，得出大尺度和高速条件下的碰撞结果11-12。Raducan等11使用iSALE以双小行星重定向测试任务为例进行了2D数值模拟，发现动量传递效率系数主要受材料初始孔隙度和内摩擦系数的影响。Rainey等12使用CTH模拟了DART撞击均质小行星模型，分析了模型参数对动量传递效率系数、陨石坑直径和深度的影 收稿日期：2023-04-07 修回日期：2023-05-07基金项目：国家自然科学基金（62227901，12202068）；民用航天预研项目（D020304）；国家国防科技工业局空间碎片与近地小行星防御科研项目（KJSP2020020101）第 10 卷 第 4 期深 空 探 测 学 报（中英文）Vol.10 No.42023 年 8 月Journal of Deep Space ExplorationAugust 2023响，结果表明动量传递效率系数、陨石坑直径和深度主要由材料的强度和孔隙率决定，但相关模拟未与试验结果比较，验证模型参数的准确性。因此，本文首先开展球形铝弹丸撞击玄武岩靶试验，获得不同撞击速度下成坑的结果和动量传递系数。利用AUTODYN软件光滑粒子流体动力学（Smoothed Particle Hydrodynamics，SPH）法进行数值模拟，将模拟结果与实验室规模试验数据进行了比较，验证了模型参数的准确性。在此基础上，模拟了最高10 km/s撞击速度对抛射物质量、速度分布和动量传递系数的影响。结合试验和数值模拟结果，获得动能撞击密实小行星的动量传递系数相似律，揭示撞击速度对动量传递系数的影响规律。1 超高速撞击分析 1.1 玄武岩靶123.3 mm123.350 mm10050 mm50玄武岩经水钻钻取成棒材，然后用磨床进行打磨得到如图1所示的圆柱靶，尺寸 mm，同时制备了测量岩石力学特性的试验试件，单轴压缩试验试件尺寸 mm，巴西圆盘劈裂试验试件尺寸 mm。在应变速率为1 104/s的条件下，采用单轴压缩和巴西圆盘试验测定了试件的抗压和抗拉强度。测得玄武岩单轴抗压和抗拉强度的平均值分别为146.0 MPa和9.8 MPa。通过测量抗拉试件质量和体积，得到玄武岩的密度2.876 g/cm3。靶板单轴压缩试件巴西圆盘劈裂试件 图 1 玄武岩靶和试件Fig.1 Basalt targets and specimens 为直接测量撞击玄武岩的动量传递系数，玄武岩靶被挂在2根42 cm的尼龙绳上，悬挂在靶室中央，形成弹道摆。为防止岩石材料在受高速撞击时崩落破坏，在玄武岩的侧面缠绕几圈玻璃纤维增强塑料胶带，如图2所示。1.2 弹丸和弹托试验采用如图3（a）所示的直径为6 mm球形2024铝合金弹丸，用聚碳酸酯弹托发射，弹托强制分离系统将弹丸与弹托分离，然后用磁测速测量弹丸速度13。为测量无弹丸情况下弹托的速度，将弹托切分为上下两瓣，弹托中心开槽，槽内放置一块铝片后，将分瓣弹托组合在一起，形成内部包含铝片的弹托。弹托内部铝片用于触发磁测速，测量无弹丸情况下弹托的速度，如图3（b）所示。1.3 试验系统和装置试验加载系统如图4所示，二级轻气炮加速弹丸和弹托，弹丸和弹托出炮口后，用弹托分离器将弹丸与弹托分离14，通过磁测速获得弹丸速度信号。为测量无弹丸情况下的弹托速度，将磁测速装置安装在弹托分离器前面，弹托首先通过磁测速装置，弹托内铝片触发磁测速，获得弹托速度信号。速度信号触发高速摄像，高速摄像记录弹道摆的整个摆动过程。高速摄像的拍摄帧率为1万帧/s，靶室内高速摄像机的对面设置了1个LED灯补光。靶室内环境的空气压力保持在100 Pa左右。当弹丸碰撞靶板后，轻气炮高压气室中的氮气和一级段泵管中的氢气会尾随弹丸进入靶室。为最大限度地减小轻气炮发射气体对目标的影响，在二级发射管与靶室间放置了一个圆形档板，档板中心开一个直径20 mm的圆孔，一方面可减小轻气炮发射气体对目标动量的影响；另一方面可减少进入靶室的黑烟，影响高速摄像的拍摄。1.4 动量传递系数的计算和误差分析 1.4.1 计算方法从高速摄像记录弹道摆运动过程图中找出弹道摆的初始位置和最大位置处如图5所示，将2幅图片导入 挡板中心孔靶观察窗吊绳 图 2 双吊绳弹道摆Fig.2 Double rope pendulum 组合弹托，内部含铝片铝片槽分瓣弹托弹丸弹丸和弹托弹托（a）弹丸和弹托（b）无弹丸弹托（b）图 3 弹丸和弹托Fig.3 Projectile and sabot第 4 期刘文近，等：超高速动能撞击密实小行星的动量传递规律421xmaxmaxhAdobe Photoshop软件，测量得到弹道摆最大水平位移。根据弹道摆的最大水平位移xmax和几何关系可算出弹道摆的最大摆角和弹道摆上升高度。hllxmax靶靶初始位置最大水平位移xmax=8.22 cm 图 5 撞击误差分析Fig.5 Impact error analysis max=arcsin(xmax/l)(1)lcos=lh h=l(1cos)(2)v根据能量守恒，计算弹丸撞击靶板后靶板获得的速度增量 为v=2gl1cos arcsin(xmax/l)(3)xbdbdmaxbd为排除发射气体对动量传递的影响，开展相同工况下不装载弹丸发射试验，标定发射气流引起弹道摆水平位移和最大摆角。弹丸引起弹道摆角度的实际变化 等于弹丸撞击弹道摆的最大摆角减去不装载弹丸发射试验最大摆角，实际动量传递系数的大小为=MmpU2gl1cosarcsin(xmax/l)arcsin(xbd/l)(4)mpUMg其中：为弹丸质量；为根据磁测速信号计算的弹丸速度；为撞击后靶和吊绳的质量；为重力加速度。1.4.2 误差分析x1在靶板质心建立Oxyz坐标系，如图6所示，坐标原点位于靶板质心，弹丸速度和 轴的夹角记为，撞击角度对动量传递系数带来的误差为1=1cos(5)0撞击角为10时，撞击角对动量传递系数带来的误差1.6%。实际试验中弹道摆垂直悬挂在靶室中，弹丸出炮口后均沿水平方向运动。因此，弹丸均垂直撞击靶板表面，因此撞击角度带来的误差约为0。zze超高速撞击动量传递测试系统如图6所示，当撞击点在靶板表面中轴线上下移动，如图6（b）所示，由于撞击点改变对动量系数无影响，当撞击点在靶板表面中轴线左右移动时，撞击点距 轴线距离为，撞击后弹道摆向后运动过程中同时发生旋转，因此试验前需要进行激光对中。由图5可知，弹道摆未发生偏转，试验带来误差可忽略不计。zxyOUe2zxyOU（a）撞击角度误差分析（b）撞击点误差分析吊绳l撞击速度l吊绳圆柱靶圆柱靶zyOO 图 6 超高速撞击动量传递测试系统示意图Fig.6 Schematic figures of the hypervelocity impact momentum transfertest system 泵管高压气室氢气活塞发射管高速摄像强制托磁测速示波器H2N2靶室靶弹托和弹丸观察窗吊绳锥段挡板膜片LED灯llxmax弹托和内部铝片电脑无弹丸情况 图 4 弹道摆最大水平位移Fig.4 Maximum horizontal displacement of pendulum422深空探测学报（中英文）2023年1.5 结果分析不同撞击速度下成坑直径、深度和动量传递系数如表1所示。撞击速度从2.3 km/s增加到3.9 km/s，从1.96增加到2.51。不同撞击速度下，弹丸撞击玄武岩靶成坑形貌基本一致。6 mm弹丸以3.9 km/s撞击玄武岩典型成坑结果如图7所示。玄武岩靶上形成的空腔由一个碗状中心小坑（图7中白色虚线圈内）和一片较大的浅层剥落区（图7中白色和红色虚线圈内）组成，与之前超高速撞击致密岩石成坑形貌一致7。碗状的中心坑是由冲击开挖流产生的，浅层剥落区是冲击波与稀疏波在自由表面的相互作用导致材料发生拉伸破坏形成的15。表 1 超高速撞击条件下的试验参数及结果Table 1 Experimental parameters in hypervelocity impactconditions and results序号mp/gU/（kms1）M/gmejectaxmaxxbd或/cmD/mmP/mm10.312.304 359.413.38.221.964210.020.313.474 355.633.313.602.395212.030.313.904 252.832.815.752.515513.540.002.554 388.21.5750.003.514 388.91.6260.003.914 385.61.70mejecta注：p为成坑深度；D为成坑直径；为抛射物质量。1 cm 图 7 撞击速度为3.9 km/s时，玄武岩的成坑结果Fig.7 Craters formed in basalt，with impact velocity of 3.9 km/s 2 数值模型 2.1 算法及模型AUTODYN软件与其它仿真软件在计算陨石坑大小、撞击过程中峰值压力和动量增强系数大小进行了基准比较16-18，广泛用于模拟动能侵彻、空间碎片防护、陨石成坑等超高速撞击问题19-22。动能撞击偏转小行星涉及超高速和大变形的成坑过程，相对于基于网格的方法，SPH模拟超高速撞击成坑问题不但可以避免网格畸变，而且可以追踪在陨石坑形成过程中抛射1出的粒子质量和速度。基于AUTODYN软件二次开发功能，可以导出所有SPH粒子在不同时间的空间位置、质量和速度信息。根据所有抛射物的质量和速度信息，将撞击方向上产生的所有抛射物粒子动量分量相加，动量传递系数等于抛射物在撞击方向上的总动量除以抛射物初始动量23。利用AUTODYN软件建立二维轴对称模型进行数值计算如图8所示，球形弹丸直径6 mm，采用航天器常用材料Al2024。靶厚度和直径均为123.3 mm的玄武岩，弹丸半径20个粒子，SPH粒子大小为0.15 mm16。Al2024玄武岩 图 8 简化计算模型Fig.8 Simplified calculation model 2.2 材料模型及参数Al2024采用Shock状态方程、Johnson-Cook强度模型和最大拉应力失效，失效应力取2.6 GPa24，材料参数来自AUTODYN材料库，具体材料参数见表2和表3。表 2 Al2024和玄武岩的Shock状态方程参数19，25Table 2 Parameters of shock for Al 2024-T4 and basalt材料0/（gcm3）c0/（ms1）s1Al-20242.7855 3281.3382玄武岩2.8762 4001.6001 表 3 Al-2024 的Johnson-Cook本构模型和失效参数26Table 3 The Johnson-Cook model and failuparameters for Al-2024材料A/MPaB/MPancmTm/KTr/KAl-20242654260.0190.340.4935300 YP玄武岩采用Shock状态方程，玄武岩采用Drucker-Prager强度模型，模型考虑了屈服强度与静水压力之间的线性关系。Y=Y0+f P(6)Y0f玄武岩初始屈服强度取146 MPa，摩擦系数取0.611，玄武岩采用最大拉应力失效，失效应力取50 MPa27。第 4 期刘文近，等：超高速动能撞击密实小行星的动量传递规律4232.3 计算结果与试验结果比较6 mm铝弹丸3.9 km/s速度撞击玄武岩的成坑结果如图9所示。从0.050.2 ms成坑深度不变，表明开坑完成。撞击0.2 ms后，层裂碎片与靶完全分离，速度不会再发生改变。层裂碎片速度较低，完全离开靶板表面需要较长，可根据导出的SPH粒子空间位置、质量和速度信息，计算反溅抛射物和层裂碎片的动量。因此，为节约计算时间，取不同撞击速度下撞击0.2 ms后，计算反溅抛射物和层裂碎片的动量。60402002040606040200204060 x/mm y/mm0.05 ms0.2 ms撞击速度=3.90 km/s层裂碎片 图 9 撞击后成坑结果，左0.05 ms，右0.2 msFig.9 Crater results，left 0.05 ms，right 0.2 ms PI撞击后0.2 ms不同撞击速度形成的抛射物的质量速度分布如图10所示。对抛射物质量速度分布曲线积分，得到抛射物累积动量，动量传递系数1=PI/mU。水平横线上是层裂碎片的质量和动量，层裂大大提高了抛射物的质量，如图10（a）所示。但层裂碎片速度较低，贡献动量较小如图10（b）所示。在不同速度下，使用AUTODYN模拟预测的陨石坑直径和深度分别与试验相差4.7%10.9%和3.7%8%。数值模拟得到的动量传递系数与试验的最大误差为10.2%。通过与不同速度的试验结果对比，验证了动量传递系数的统计方法、所选算法和模型参数的合理性。在此基础上，模拟了510 km/s这4种不同撞击速度撞击玄武岩靶，研究撞击速度对抛射物质量、速度分布和动量传递系数的影响，结果如表4所示。随着撞击速度提高，弹丸和靶板可能会发生融化和气化，Shock状态方程不能直接计算得到温度或焓，只能间接估算，因此本文未考虑融化和气化对抛射物产生和质量速度分布的影响。之前的数值模拟结果表明，在10 km/s以下，弹丸和靶板发生融化和气化的比例较低，Shock状态方程是适用的27-28。但在更高速度的情况下，需考虑熔化和气化问题，选择可以描述物质的凝聚态（固液）、膨胀态（气）和混合态的状态方程29，研究融化和气化产物对动量传递的影响。为防止高速下靶板被撞碎，将靶板尺寸扩大1倍。当陨石坑直径达到几米时，在陨石坑表面形成大面积的损伤，损伤区域的抛射物无法获得足够的速度离开小行星表面30。因此，将小型试验和模拟结果外推到小行星尺度需要修改试验和模拟结果，去除层裂机制产生的抛射物对动量传递的贡献。105104103102101100105104103102101100050100150200累计抛射物质量 M(v)/m抛射速度(v/U)（a）速度大于v/U无量纲抛射物质量M(v)/m的分布（b）1对无量纲抛射物速度v/U的依赖关系2.30 km/s 3.47 km/s 3.90 km/s层裂开挖0.51.01.51撞击速度(v/U)2.30 km/s 3.47 km/s 3.90 km/s剪切抛射形成抛射物贡献的动量层裂动量 图 10 不同速度下的抛射物的质量和速度分布Fig.10 Mass-velocity distribution of ejecta at different velocities 2.4 数值计算结果及分析 2.4.1 抛射物的质量与速度分布规律v速度大于 的累积抛射物质量M（v）分布曲线如图11所示，抛射物的速度逐渐变小，抛射物累计质量分布趋于平缓，接近抛射物质的总质量。不同撞击速度下累计抛射物质量M（v）分布规律相似，撞击速 表 4 数值模拟与试验结果的比较Table 4 Comparison between numerical simulation andexperimental results速度/（kms1）成坑直径/mm成坑深度/mm试验 模拟误差试验 模拟 误差试验 模拟 误差2.3042444.7%109.28.0%1.961.7614%3.4752487.7%1211.2 6.7%2.392.29 4.2%3.90554910.9%13.5143.7%2.512.41 3.9%424深空探测学报（中英文）2023年13度越大，抛射物的累计质量越大。撞击速度从5 km/s增加到10 km/s，从1.64增加到2.72（见表5）。抛射物的质速度分布函数可简化为幂律分布10，图11中黑色直线斜率等于，是耦合参数，黑色直线斜率等于1.65。103104103102101100102101100101102103累积抛射物质量 M(v)/m撞击速度 v/U5.0 km/s 6.5 km/s 8.0 km/s 10.0 km/s vM(v)图 11 速度大于 抛射物的质量分布规律vFig.11 Distribution of ejecta with velocity over 表 5 动量传递系数的相似律模型计算和模拟结果Table 5 Scaling law model and simulation results of momentumenhancement factor 撞击速度/（kms1）相对误差/%相似律模拟值52.822.646.86.53.163.131.083.473.372.9103.863.585.0 2.4.2 动量传递相似律相似律模型广泛用于将试验室尺度的试验和模拟结果外延至行星级别，强度区间的动量相似律10为1=KsKvs(n2(n1)UY/31()13(7)UYn=3KsKvs其中：、分别为弹丸密度、撞击速度、靶强度和密度；为耦合参数；、为材料相关的常数。Ks=MemUY/3()31(8)Kvs=v/Y/(9)mMevv其中：为弹丸质量；为速度大于 的抛射物质量；为抛射物的最小速度。vKs=0.29Kvs=0.47速度大于 的抛射物质量如图12所示，玄武岩动量相似关系常数，。将弹丸和靶板的参数带入到式（7），得到撞击密实小行星的动量传递U曲线（单位km/s）：=1+0.007 2U0.65(10)11105.0 km/s6.5 km/s8.0 km/s10.0 km/s1000.10.11051041030.010.01v/Y/mU()33v1()MeY/v图 12 速度大于 的抛射物质量vFig.12 Mass with ejection velocity 在10 km/s撞击速度范围内模拟的结果（图13）和相似律结果的最大误差在10%之间（见表5）。因此，在10 km/s以内的速度范围内，动量增强系数 随撞击速度的0.65次方的增大而增大，相律模型可以很好地预测动量增强系数。0246810012345678910撞击速度/(kms1)模拟值Y=146 MPa相似律值Y=146 MPa相似律值Y=4.2 MPa 图 13 密实小行星动量增强系数与撞击速度的关系Fig.13 Relationship between momentum enhancement factor andimpact velocity YxmaxY=Y0(L0/L)1/iiL0Y0L0根据威布尔理论，脆性岩石的强度随着其质量长度的增大而减小，其中 通常为2331，靶板质量长度为时靶板强度大小为。将实验室靶板质量长度外推至近地小天体大小，需要考虑岩石强度的变化。因此，考虑近地小天体尺寸的动量传递因子 为1=(01)(LL0)312i(11)0=0.55i=21其中：为试验室条件下动量传递系数的大小，由式(11)可知，对于，随近地小天体尺第 4 期刘文近，等：超高速动能撞击密实小行星的动量传递规律425U寸的0.16次方增大而增大。将试验室岩石的尺寸12 cm增加到140 m，近地小天体岩石的强度从146.0 MPa降低到4.2 MPa，则撞击140 m密实小行星的动量传递曲线（见图13）（单位ms1）为=1+0.015 6U0.65(12)3 结论本文利用二级轻气炮进行了6 mm铝弹丸以24 km/s超高速撞击玄武岩靶试验，获得了在不同速度下撞击玄武岩的动量传递系数。将模拟结果和试验结果进行了比较，验证了计算模型和参数的准确性。模拟和相似律模型结果均表明，密实小行星的动量增强系数明显依赖撞击速度。撞击岩石目标数值模拟结果与相似律模型计算结果吻合较好，相似律模型可以很好地预测动量增强结果，在10 km/s以内的速度范围内，动量增强系数 随撞击速度的0.65次方的增大而增大。参考文献参考文献 PERNA D，BARUCCI M A，FULCHIGNONI M.The near-Earthobjects and their potential threat to our planetJ.The Astronomy andAstrophysics Review，2013，21（1）：1-28.1 DEARBORN D S P，BRUCK SYAL M，BARBEE B W，et al.Optionsand uncertainties in planetary defense：impulse-dependent response andthe physical properties of asteroidsJ.Acta Astronautica，2020，166：290-305.2 DALY R T，ERNST C M，BARNOUIN O S，et al.Successful kineticimpact into an asteroid for planetary defenseJ.Nature，2023，616：443-447.3 PELTON J N，ALLAHDADI F.Handbook of cosmic hazards andplanetary defenseM.London：Springer，2015.4 SUGIMOTO Y，RADICE G，CERIOTTI M.Hazardous near Earthasteroid mitigation campaign planning based on uncertain informationon fundamental asteroid characteristicsJ.Acta Astronautica，2014（103）：333-357.5 GAULT D E，SHOEMAKER M E，MOORE H J.Spray ejected fromthe lunar surface by meteoroid impactR.Moffett Field，Calif：AmesResearch Center，1963.6 HOERTH T，SCHFER F，HUPFER J，et al.Momentum transfer inhypervelocity impact experiments on rock targetsJ.ProcediaEngineering，2015，103：197-204.7 CHOUREY S，KOSCHNY D，ROTT M，et al.Determining themomentum transfer in regolith-like targets using the TUM/LRT electro-thermal acceleratorJ.Planetary and Space Science，2020，194：105112.8 FLYNN G J，DURDA D D，MOLESKY M J，et al.Momentum transferin hypervelocity cratering of meteorites and meteorite analogs：implications for orbital evolution and kinetic impact deflection ofasteroidsJ.International Journal of Impact Engineering，2020，136：103437.9 HOLSAPPLE K A，HOUSEN K R.Momentum transfer in asteroidimpacts.I.theory and scalingJ.Icarus，2012，221（2）：875-887.10 RADUCAN S D，DAVISON T M，LUTHER R，et al.The role ofasteroid strength，porosity and internal friction in impact momentumtransferJ.Icarus，2019，329：282-295.11 RAINEY E S G，STICKLE A M，CHENG A F，et al.Impact modelingfor the Double Asteroid Redirection Test（DART）missionJ.International Journal of Impact Engineering，2020，142：103528.12 施尚春，陈攀森，黄跃.高速弹丸的磁感应测速方法J.高压物理学报，1991，5（3）：205-214.SHI S C，CHEN P S，HUANG Y.Velocity measurement of magnetinduced system for projectile J.Explosion and Shock Waves，1991，5（3）：205-214.13 唐恩凌，李志秋.二级轻气炮中组合式弹托分离器的应用J.沈阳理工大学学报，2012，31（1）：84-87.TANG E L，LI Z Q.The application of combined type dissociator ofprojectile from its seat in two-stage light gas gun J.Journal ofShenyang Ligong University，2012，31（1）：84-87.14 POELCHAU M H，KENKMANN T，HOERTH T，et al.Impactcratering experiments into quartzite，sandstone and tuff：the effects ofprojectile size and target properties on spallationJ.Icarus，2014，242：211-224.15 STICKLE A M，SYAL M B，CHENG A F，et al.Benchmarking impacthydrocodes in the strength regime：implications for modeling deflectionby a kinetic impactorJ.Icarus，2020（338）：113446.16 SUETSUGU R，TANAKA H，KOBAYASHI H，et al.Collisionaldisruption of planetesimals in the gravity regime with iSALE code：comparison with SPH code for purely hydrodynamic bodiesJ.Icarus，2018，314：121-132.17 PIERAZZO E，ARTEMIEVA N，ASPHAUG E，et al.Validation ofnumerical codes for impact and explosion cratering：impacts onstrengthless and metal targetsJ.Meteoritics&Planetary Science，2008，43（12）：1917-1938.18 REN S，ZHANG Q，GAO F，et al.Impact resistance mechanism ofreactive material bumper for spacecraft whipple shield：experimentsand numerical simulationsJ.Aerospace Science and Technology，2022，126：107646.19 LU Y Y，ZHANG Q M，XUE Y J，et al.Hypervelocity penetration ofconcrete targets with long-rod steel projectiles：experimental andtheoretical analysisJ.International Journal of Impact Engineering，2021，148：103742.20 SAITO T，KAIHO K，ABE A，et al.Hypervelocity impact ofasteroid/comet on the oceanic crust of the earthJ.International Journalof Impact Engineering，2008，35（12）：1770-1777.21 SAITO T，KAIHO K，ABE A，et al.Numerical simulations ofhypervelocity impact of asteroid/comet on the EarthJ.InternationalJournal of Impact Engineering，2006，33（1-12）：713-722.22 SYAL B，J O M，MILLER P L.Deflection by kinetic impact：sensitivity to asteroid propertiesJ.Icarus，2016，269：50-61.23 邸德宁，陈小伟.碎片云SPH方法数值模拟中的材料失效模型J.爆炸与冲击，2018，38（5）：948-956.DI D N，CHEN X W.Material failure models in SPH simulation ofdebris cloudJ.Explosion and Shock Waves，2018，38（5）：948-956.24 AHRENS T J，JOHNSON M L.Shock wave data for rocksM.Washington，DC.：American Geophysical Union，1995，35-44.25426深空探测学报（中英文）2023年 JOHNSON G R，COOK W H.A constitutive model and data for metalssubjected to large strains，high strain-rates and high temperaturesC/Proceedings of Seventh International Symposium on Ballistics.TheHague，The Netherlands：s.n.，1983.26 WALKER J D，CHOCRON S.Momentum enhancement inhypervelocity impactJ.International Journal of Impact Engineering，2011，38（6）：1-7.27 RADUCAN S D，JUTZI M，DAVISON T M，et al.Influence of theprojectile geometry on the momentum transfer from a kinetic impactorand implications for the DART missionJ.International Journal ofImpact Engineering，2022，162：104147.28 WNNEMANN K，COLLINS G S，OSINSKI G R.Numericalmodelling of impact melt production in porous rocksJ.Earth andPlanetary Science Letters，2008，269（3-4）：530-539.29 LUTHER R，ZHU M，COLLINS G，et al.Effect of target propertiesand impact velocity on ejection dynamics and ejecta depositionJ.30Planetary and Space Science，2018，53（8）：1705-1732.HOLSAPPLE K A.On