Section 6.3 - Venus SmallSat Constellation - Deflection Manuevers

[1]:
import numpy as np
from AMAT.maneuver import ProbeProbeDeflection
from AMAT.maneuver import ProbeOrbiterDeflection
[2]:
deflection1 = ProbeProbeDeflection(arrivalPlanet="VENUS",
                                    v_inf_vec_icrf_kms=np.array([-3.26907094,  0.67649494, -1.0697747 ]),
                                    rp_probe1=(6051.8+103.85)*1e3,  psi_probe1=1.05988*np.pi, h_EI_probe1=150e3,
                                    rp_probe2=(6051.8+103.85)*1e3,  psi_probe2=1.03000*np.pi, h_EI_probe2=150e3,
                                    r_dv_rp=215.5)

deflection2 = ProbeProbeDeflection(arrivalPlanet="VENUS",
                                    v_inf_vec_icrf_kms=np.array([-3.26907094,  0.67649494, -1.0697747 ]),
                                    rp_probe1=(6051.8+103.85)*1e3,  psi_probe1=1.03000*np.pi, h_EI_probe1=150e3,
                                    rp_probe2=(6051.8+103.85)*1e3,  psi_probe2=1.00000*np.pi, h_EI_probe2=150e3,
                                    r_dv_rp=189.6)

deflection3 = ProbeProbeDeflection(arrivalPlanet="VENUS",
                                    v_inf_vec_icrf_kms=np.array([-3.26907094,  0.67649494, -1.0697747 ]),
                                    rp_probe1=(6051.8+103.85)*1e3,  psi_probe1=1.00000*np.pi, h_EI_probe1=150e3,
                                    rp_probe2=(6051.8+103.85)*1e3,  psi_probe2=0.97000*np.pi, h_EI_probe2=150e3,
                                    r_dv_rp=164)

deflection4 = ProbeProbeDeflection(arrivalPlanet="VENUS",
                                    v_inf_vec_icrf_kms=np.array([-3.26907094,  0.67649494, -1.0697747 ]),
                                    rp_probe1=(6051.8+103.85)*1e3,  psi_probe1=0.97000*np.pi, h_EI_probe1=150e3,
                                    rp_probe2=(6051.8+103.85)*1e3,  psi_probe2=0.94000*np.pi, h_EI_probe2=150e3,
                                    r_dv_rp=138.5)

deflection5 = ProbeOrbiterDeflection(arrivalPlanet="VENUS",
                                    v_inf_vec_icrf_kms=np.array([-3.26907094,  0.67649494, -1.0697747 ]),
                                    rp_probe=(6051.8+103.85)*1e3,  psi_probe=0.94000*np.pi, h_EI_probe=150e3,
                                    rp_space=(6051.8+400)*1e3, psi_space=1.0*np.pi,
                                    r_dv_rp=112.5)
[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.75224316 -4.61251793  0.6152991 ] m/s
Divert manuever DV magnitude: 4.714 m/s
---------------------------------
Deflection Maneuver #2
Divert manuever DV: [-0.67485826 -5.31732722  0.23444224] m/s
Divert manuever DV magnitude: 5.365 m/s
---------------------------------
Deflection Maneuver #3
Divert manuever DV: [-0.56133357 -6.15017052 -0.27012275] m/s
Divert manuever DV magnitude: 6.182 m/s
---------------------------------
Deflection Maneuver #4
Divert manuever DV: [-0.40090375 -7.21278622 -0.9523719 ] m/s
Divert manuever DV magnitude: 7.286 m/s
---------------------------------
Deflection Maneuver #5
Divert manuever DV: [ 0.35802126 17.78715771  4.02682792] m/s
Divert manuever DV magnitude: 18.241 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: 4.007605 days
TOF from probe release to atm. entry interface: 3.500867 days
TOF from probe release to atm. entry interface: 3.001488 days
TOF from probe release to atm. entry interface: 2.506011 days
TOF from probe release to atm. entry interface: 2.003604 days