2022年小行星高速碰撞问题与建模收集.pdf

1 小行星高速碰撞问题与建模名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 5 页 - - - - - - - - - 2 摘 要： 宇宙中分布着大量的小行星，地球在运行过程中与太阳系中的小行星发生碰撞是可能的 这种碰撞是高速的，它对地球的生态环境和人类生存的影响有多大是我们非常关心的问题 我们通过建立一个小行星高速碰撞地球南极的模型，分析碰撞发生的全过程，将该过程分为进入大气过程与碰撞过程以及后期影响两部分，对能量的分配和吸收作出假设分析估计碰撞对南极大陆的直接影响关键字 : 小行星；高速碰撞；建模1 问题的提出我们都知道， 宇宙中分布着大量的小行星，在太阳系中， 小行星随时可能与运行过程中的地球发生碰撞，这种碰撞往往都是高速的，而这种高速的碰撞会对地球的生态环境和人类的生存带来相当的影响在宇宙中， 这种随时的碰撞是人们所不能阻止的，所以人们的注意力自然转移到如何避免或降低这种高速碰撞为人们所带来的影响这一问题上为了解决这一问题， 我们可以从数学的角度建立一个合理的模型来模拟这种碰撞，然后对其进行全方位深层次的分析， 从而估计出碰撞所带来的直接影响和潜在影响，为更有针对性地降低和解决这一问题提供一个精确且合理的平台。

2 问题分析我们将这个复杂过程分为三个部分：11 进入大气过程小行星从外太空高速进入大气层，收到空气阻力的作用会有一定的减速1.2碰撞过程高速运动的小行星体带有巨大的动能接触南极冰层后，与冰层发生剧烈碰撞，在地面上留下一个巨大的陨石坑，同时溅射出大量的物质巨大的动能迅速转化为内能，产生高温高压，瞬间将大量的冰层融化为水与蒸汽此外， 还有一部分能量转化为溅射物质的动能以及物体变形的动能地面吸收能量后引发地震，海啸等一系列次生灾害1.3 后期影响碰撞的后期影响主要集中灾害在生态环境的影响，包括融水对海平面高度的影响，陨石中物质流入海洋的影响，溅射物质的影响等等3 问题假设（1）入射的小行星形状近似为一个球体，其半径为1000m; （2） 小行星掉落的地点在地球的南极；（3）小行星进入大气层后沿着一条直线运动，不发生偏移；（4） 小行星进入大气层时不发生爆炸，与地面接触后也不发生反弹；（5）小行星主要由铁和镍构成，这基本是陨石中密度最大的物质，我们可以估计出最坏的情况；（6）南极冰层平均厚度约为2000m; 4 符号系统名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 5 页 - - - - - - - - - 3 star小行星的密度V0 小行星进入大气层时的速度dstar 小行星的直径小行星撞击地面的角度T0 冰层的初始温度Mstar 小行星的质量HH E E 5 模型建立5.1 进入大气过程我们首先计算一下小行星的质量：Mstar= star dstar 3/6 ,小行星的初始动能Estar= MstarV02/2。

在小行星进入大气层的过程里，不可避免的要受到空气阻力的影响，损失动能， 摩擦生热根据这个我们来计算一下空气阻力的影响，通过查找资料我们发现空气阻力的计算公式为R = CAV2式中 V代表这个物体的速度，A代表该物体与空气的接触面积，为空气的密度C是一个与空气相关的常数，一般为0.2到 0.4 kg/m3我们想计算出物体与空气磨擦时，动能转化为热能的上限我们假设大气层的厚度为100km 100km以内的大气随着高度的减小，其密度的对数 （以10为底） 是一个线性线性增加的过程由此我们假设大气的归一化密度为air=10ad+b( 即假设海平面处的空气密度为1) ，式中的 d为距海平面的高度，a和b为两个系数现在我们先来计算一下 a和b由于 100km 相对于地球的半径6400km来说很小，所以当小行星在大气层运动时可以把地面看作一个水平面小行星进入大气层时沿着一个角度，所以从 100km 处大气层到海平面经过的路程为S0=105/sin 在海平面处，air=10ad+b=1，可得 ad+b=0，而 d=0，故 b也等于 0在距海平面 S0处，空气密度大约是海平面的10-6，10-6=10a*S0+b, 所以 -6=a*105/sina,可得a=-6 sin/105。

故在大气层的范围内，空气阻力R = C AV2=C*10aAV2空气阻力做功为W = -0sina/105Rdh= )101(6ln10sinCAV10625sin10ln61025CAV为了估计一下最坏的情况，我们取C=0.4, A为r2,此处 r取500m(因为小行星的面直径约名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 5 页 - - - - - - - - - 4 为1000m) 我们计算一下在不同角度下空气阻力做功占总动能的比例P，结果如下：角度 10 20 30 40 50 60 70 80 90 比例 P 0.163 0.083 0.057 0.044 0.037 0.033 0.030 0.029 0.028 从这个结果可以看到空气阻力的影响非常小，作用最强时也仅消耗总动能的0.163%，所以我们在下面的分析里忽略空气阻力的影响，即动能在进入大气层时不损失5.2 碰撞过程小行星以超高速撞击地球的南极点，根据碰撞遵循质量守恒原理、动量守恒原理和高速碰撞理论，强烈的冲击波将引起强大的冲击压力。

首先让我们初步估计一下撞击的能量：撞击产生的能量来自于小行星的动能根据假设， 令小行星的速度V0=30km/s, 小行星的密度 =5g/cm3,小行星的直径D=1000m,所以小行星的体积Mstar= star dstar 3/6， 质量 M=2.621012kg，因此，小行星的动能为E=11.78 1020J根据对大量碰撞所产生的能量的估计，直径为 1000m的小行星撞击地球所产生的能量相当于100M TNT 炸药所产生的能量（其中M为小行星质量） 1kg TNT炸药爆炸，会放出4560kJ 热量，记 C=4560kJ/kg ，则有 E=100M C11.94 1020J误差是由于忽略了小行星穿过大气层所损失的能量之后我们来计算一下冲击波压力Pc（其中 1为冰块密度， 1为水的密度， star为小行星平均密度， V0为小行星速度， C1为冲击波在冰块中的传播速度，C2为冲击波在行星内的传播速度， Vc为冲击波的反射速度，P0为冰块的初始压力，Pc为冲击波的反射压力，为与材料有关的弹塑性常量对于冰块，根据质量守恒和动量守恒，得：1 C1 =1（ C1 - Vc ） (1) （1 C1） C1 + P0 = 1（ C1 - Vc ） （C1 - Vc ）+ Pc （2）由（ 1） 、 （2）式得：Pc - P0 =1 C1 Vc (3）对于小行星Pc - P0 =2 C2（V0- Vc）（4）由（ 3） 、 （4）式得：Pc - P0 =（1 C1star C2） V0/ （1 C1 +star C2）（5）根据高速碰撞理论：C = （/ ）0.5（6）对于冰块 1 = 0.9 2；初始压力p00，由（1）式到（ 5）式，得 Pc约为 0.2 1016Pa 。

在这个地区， 温度在短时间内快速升高，大量冰块融化向各个方向扩散，同时大量水蒸气进入大气根据功能原理，冲击压力Pc对地球所做的功W为：W = 0.5Pcv (7) 其中 v是压缩形变因为是超高速碰撞，v和 Pc不是线性关系，两者服从下面的非线性关系：v = aPc2 + bPc (8）这里 a 和 b 是两个常数（ a0） 我们把由冲击波产生的能量分配为被南极大陆冰块吸收的能量E1, 深层未融化冰块吸收的能量 E2，升华成水蒸气的冰块吸收的能量，也就是说，E1、E2之和远大于E3、Ek之和，若总能量为 E，我们假定E1 + E2 = 80%E，E3 = 10%E，Ek = 10%E 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 5 页 - - - - - - - - - 5 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 5 页 - - - - - - - - - 。