• 基于离散元法的陨石撞击月球熔岩管数值模拟分析

    Numerical Simulation Analysis of Meteorite Impact on Lunar Lava Tubes Based on the Discrete Element Method

    • 为了探究月球熔岩管受陨石撞击时的结构稳定性,基于离散元法(Discrete Element Method,DEM)建立二维数值模型,模拟不同撞击速度及角度的陨石对不同截面形状的月球熔岩管的动态撞击情景。模型采用平行黏结接触模型表征月岩力学特性,通过撞击坑形态、应力时程曲线及接触力链演化分析表明:截面形状影响微弱,三种截面熔岩管的撞击坑形态与应力分布基本一致,表明应力波传播受几何构型影响较小;撞击速度主导破坏,速度增至30 km/s时,撞击坑深度与宽度较撞击速度20 km/s时分别增加136.11%和172.47%,应力极值有所提升,熔岩管顶胶结键破坏形成贯通剪切带,结构失稳风险显著提升;撞击角度调制损伤路径分别是45°倾斜撞击坑深度较垂直撞击减小26.7%,能量水平耗散削弱深部破坏,但高速下结构稳定性仍存在安全隐患。数值模拟验证了陨石撞击时,陨石速度>角度>截面形状对月球熔岩管管顶岩体稳定性的影响优先级,揭示了应力波-颗粒介质的损伤传播机制,为月球熔岩管基地选址提供参考。

       

      Abstract: To investigate the structural stability of lunar lava tubes subjected to meteoroid impacts, a two-dimensional numerical model was developed using the Discrete Element Method (DEM). Dynamic impact scenarios were simulated, where meteoroids with varying impact velocities and incident angles strike lava tubes with different cross-sectional geometries. The parallel bond contact model was employed to characterize the mechanical behavior of lunar basaltic rock at the particle scale. Analyses of crater morphology, stress-time histories, and contact force chain evolution indicate that the influence of cross-sectional shape is relatively minor. The three tube geometries studied exhibit comparable crater profiles and stress distributions, suggesting that the propagation of impact-induced stress waves is only weakly affected by geometric configuration. Impact velocity is identified as the primary factor governing structural damage. When the impact velocity increases to 30 km/s, the crater depth and width increase by 136.11% and 172.47%, respectively, compared to those at an impact velocity of 20 km/s. The associated peak stresses also rise significantly, and extensive breakage of cemented bonds occurs near the tube crown, forming a continuous shear band that substantially increases the risk of structural instability. The incident angle mainly influences the damage pattern: for an oblique impact at 45°, the crater depth decreases by 26.7% relative to that under vertical incidence, and the reduced energy transmission mitigates deep-seated damage. Nevertheless, at high velocities, significant stability hazards persist even for oblique impacts. The numerical results demonstrate that the relative importance of impact parameters to the stability of the lava-tube roof rock mass follows the order: impact velocity > impact angle > cross-sectional geometry. The simulations further reveal the coupled mechanisms of stress-wave propagation and particle-scale damage evolution in granular rock media. This study provides scientific insights into the vulnerability of lunar lava tubes under meteoroid bombardment and offers a theoretical reference for site selection and safety evaluation of future lunar lava-tube habitats.

       

    /

    返回文章
    返回