火星霍曼转移轨道卫星发射速度范围及近距偏差问题问询

我计划将一颗1062kg的卫星从地球发射至火星，目标是尽可能接近火星（理论可至地表600km处）但不进入环绕轨道，需确定霍曼转移轨道的发射速度范围。当前模拟情况如下：

• 采用步长31577s的Beeman算法进行太阳系模拟（模拟已验证正常）
• 通过hohmann_transfer_parameters函数计算霍曼转移ΔV
• 设置日地-日火夹角44-46度作为发射窗口
• 在launch_satellite函数中叠加0-5000m/s的附加速度进行多轮模拟

模拟后得到的发射速度与最小火星距离关系图显示，最小距离约为1×10^6km，远大于目标的约4385km（火星半径3785km+地表600km），恳请协助排查参数计算或发射时机设置中的问题。

相关代码片段 #

def main():
    # 获取附加速度范围
    velocity_range = np.linspace(0, 5000, 20)  # 附加速度示例范围
    min_distances = []
    velocities = []
    # 遍历附加速度进行模拟
    for additional_velocity in velocity_range:
        velocity_at_launch = sim_beeman.run(additional_velocity)
        velocities.append(velocity_at_launch)
        min_distance = sim_beeman.calculate_min_distance()
        min_distances.append(min_distance)
    # 绘图代码...
def hohmann_transfer_parameters(G, mass_sun, radius_earth, radius_mars):
    a_transfer = (radius_earth + radius_mars) / 2  
    V_earth = np.sqrt(G * mass_sun / radius_earth)
    V_transfer_perihelion = np.sqrt(G * mass_sun * (2 / radius_earth - 1 / a_transfer))
    Delta_V1 = V_transfer_perihelion - V_earth
    return Delta_V1
# 其余类函数代码及发射参数示例保留原文

可能的排查方向 #

• 霍曼转移参数计算：确认radius_earth和radius_mars是否为日心轨道平均半径（地球_1AU，火星1.524AU），单位是否统一（如均转换为米）；检查G和太阳质量的取值精度是否足够。
• 发射窗口角度：霍曼转移的最优相位角并非固定44-46度，需根据转移时间计算实时相位角——卫星从地球到火星的霍曼转移时间约为259天，此时火星需要在地球前方约44度（(2π/火星轨道周期)*转移时间），但需确认角度参考系（日心黄道面内的相对角度）是否正确。
• 附加速度的叠加方向：霍曼转移的ΔV需要沿地球公转方向切线方向叠加，若速度叠加方向错误（如径向或其他方向），会导致轨道偏离霍曼转移轨道，无法接近火星。
• 模拟步长影响：31577s约为8.77小时，步长过大可能导致轨道积分误差累积，尤其是在接近火星时的轨道计算精度不足，可尝试缩小步长验证。
• 最小距离计算逻辑：确认calculate_min_distance函数是否正确计算卫星与火星的实时距离，并取整个模拟周期内的最小值，而非某一时刻的距离。

内容的提问来源于stack exchange，提问作者Anay Chadha