Section 5.3 - Mars SmallSat Constellation - Deflection Maneuvers

[1]:
import numpy as np
from AMAT.maneuver import ProbeProbeDeflection
from AMAT.maneuver import ProbeOrbiterDeflection
[2]:
deflection1 = ProbeProbeDeflection(arrivalPlanet="MARS",
                                    v_inf_vec_icrf_kms=np.array([ 2.23930484,  1.20086474, -0.73683366]),
                                    rp_probe1=(3389.5+52)*1e3,  psi_probe1=1.5*np.pi, h_EI_probe1=120e3,
                                    rp_probe2=(3389.5+52)*1e3,  psi_probe2=1.375*np.pi, h_EI_probe2=120e3,
                                    r_dv_rp=481)

deflection2 = ProbeProbeDeflection(arrivalPlanet="MARS",
                                    v_inf_vec_icrf_kms=np.array([ 2.23930484,  1.20086474, -0.73683366]),
                                    rp_probe1=(3389.5+52)*1e3,  psi_probe1=1.375*np.pi, h_EI_probe1=120e3,
                                    rp_probe2=(3389.5+52)*1e3,  psi_probe2=1.25*np.pi, h_EI_probe2=120e3,
                                    r_dv_rp=414)

deflection3 = ProbeProbeDeflection(arrivalPlanet="MARS",
                                    v_inf_vec_icrf_kms=np.array([ 2.23930484,  1.20086474, -0.73683366]),
                                    rp_probe1=(3389.5+52)*1e3,  psi_probe1=1.25*np.pi, h_EI_probe1=120e3,
                                    rp_probe2=(3389.5+52)*1e3,  psi_probe2=1.125*np.pi, h_EI_probe2=120e3,
                                    r_dv_rp=346)

deflection4 = ProbeProbeDeflection(arrivalPlanet="MARS",
                                    v_inf_vec_icrf_kms=np.array([ 2.23930484,  1.20086474, -0.73683366]),
                                    rp_probe1=(3389.5+52)*1e3,  psi_probe1=1.125*np.pi, h_EI_probe1=120e3,
                                    rp_probe2=(3389.5+52)*1e3,  psi_probe2=np.pi, h_EI_probe2=120e3,
                                    r_dv_rp=278)

deflection5 = ProbeOrbiterDeflection(arrivalPlanet="MARS",
                                    v_inf_vec_icrf_kms=np.array([ 2.23930484,  1.20086474, -0.73683366]),
                                    rp_probe=(3389.5+52)*1e3,  psi_probe=np.pi, h_EI_probe=120e3,
                                    rp_space=(3389.5+250)*1e3, psi_space=np.pi,
                                    r_dv_rp=210)
[3]:
print('Deflection Maneuver #1')
print("Divert manuever DV: "+str(deflection1.dv_maneuver_vec)+ " m/s")
print("Divert manuever DV magnitude: "+str(round(deflection1.dv_maneuver_mag,3))+ " m/s")
print('---------------------------------')
print('Deflection Maneuver #2')
print("Divert manuever DV: "+str(deflection2.dv_maneuver_vec)+ " m/s")
print("Divert manuever DV magnitude: "+str(round(deflection2.dv_maneuver_mag,3))+ " m/s")
print('---------------------------------')
print('Deflection Maneuver #3')
print("Divert manuever DV: "+str(deflection3.dv_maneuver_vec)+ " m/s")
print("Divert manuever DV magnitude: "+str(round(deflection3.dv_maneuver_mag,3))+ " m/s")
print('---------------------------------')
print('Deflection Maneuver #4')
print("Divert manuever DV: "+str(deflection4.dv_maneuver_vec)+ " m/s")
print("Divert manuever DV magnitude: "+str(round(deflection4.dv_maneuver_mag,3))+ " m/s")
print('---------------------------------')
print('Deflection Maneuver #5')
print("Divert manuever DV: "+str(deflection5.dv_maneuver_vec)+ " m/s")
print("Divert manuever DV magnitude: "+str(round(deflection5.dv_maneuver_mag,3))+ " m/s")
print('---------------------------------')
Deflection Maneuver #1
Divert manuever DV: [0.50172394 0.76333674 4.54282315] m/s
Divert manuever DV magnitude: 4.634 m/s
---------------------------------
Deflection Maneuver #2
Divert manuever DV: [1.38318623 2.65307886 4.47180516] m/s
Divert manuever DV magnitude: 5.38 m/s
---------------------------------
Deflection Maneuver #3
Divert manuever DV: [2.35903472 4.80111454 3.57218211] m/s
Divert manuever DV magnitude: 6.432 m/s
---------------------------------
Deflection Maneuver #4
Divert manuever DV: [3.36281421 7.08473581 1.55925957] m/s
Divert manuever DV magnitude: 7.996 m/s
---------------------------------
Deflection Maneuver #5
Divert manuever DV: [ 0.02474457 -0.01157163  0.95082094] m/s
Divert manuever DV magnitude: 0.951 m/s
---------------------------------
[4]:
print("TOF from probe release to atm. entry interface: "+str(round(deflection1.TOF_probe1,6))+ " days")
print("TOF from probe release to atm. entry interface: "+str(round(deflection2.TOF_probe1,6))+ " days")
print("TOF from probe release to atm. entry interface: "+str(round(deflection3.TOF_probe1,6))+ " days")
print("TOF from probe release to atm. entry interface: "+str(round(deflection4.TOF_probe1,6))+ " days")
print("TOF from probe release to atm. entry interface: "+str(round(deflection5.TOF_probe,6))+ " days")
TOF from probe release to atm. entry interface: 7.002802 days
TOF from probe release to atm. entry interface: 6.013292 days
TOF from probe release to atm. entry interface: 5.009726 days
TOF from probe release to atm. entry interface: 4.00719 days
TOF from probe release to atm. entry interface: 3.006257 days