{ "cells": [ { "cell_type": "markdown", "id": "7c826511", "metadata": {}, "source": [ "# Section 6.3 - Venus SmallSat Constellation - Deflection Manuevers" ] }, { "cell_type": "code", "execution_count": 1, "id": "b569c934", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from AMAT.maneuver import ProbeProbeDeflection\n", "from AMAT.maneuver import ProbeOrbiterDeflection" ] }, { "cell_type": "code", "execution_count": 2, "id": "c0baf7d4", "metadata": {}, "outputs": [], "source": [ "deflection1 = ProbeProbeDeflection(arrivalPlanet=\"VENUS\",\n", " v_inf_vec_icrf_kms=np.array([-3.26907094, 0.67649494, -1.0697747 ]),\n", " rp_probe1=(6051.8+103.85)*1e3, psi_probe1=1.05988*np.pi, h_EI_probe1=150e3,\n", " rp_probe2=(6051.8+103.85)*1e3, psi_probe2=1.03000*np.pi, h_EI_probe2=150e3,\n", " r_dv_rp=215.5)\n", "\n", "deflection2 = ProbeProbeDeflection(arrivalPlanet=\"VENUS\",\n", " v_inf_vec_icrf_kms=np.array([-3.26907094, 0.67649494, -1.0697747 ]),\n", " rp_probe1=(6051.8+103.85)*1e3, psi_probe1=1.03000*np.pi, h_EI_probe1=150e3,\n", " rp_probe2=(6051.8+103.85)*1e3, psi_probe2=1.00000*np.pi, h_EI_probe2=150e3,\n", " r_dv_rp=189.6)\n", "\n", "deflection3 = ProbeProbeDeflection(arrivalPlanet=\"VENUS\",\n", " v_inf_vec_icrf_kms=np.array([-3.26907094, 0.67649494, -1.0697747 ]),\n", " rp_probe1=(6051.8+103.85)*1e3, psi_probe1=1.00000*np.pi, h_EI_probe1=150e3,\n", " rp_probe2=(6051.8+103.85)*1e3, psi_probe2=0.97000*np.pi, h_EI_probe2=150e3,\n", " r_dv_rp=164)\n", "\n", "deflection4 = ProbeProbeDeflection(arrivalPlanet=\"VENUS\",\n", " v_inf_vec_icrf_kms=np.array([-3.26907094, 0.67649494, -1.0697747 ]),\n", " rp_probe1=(6051.8+103.85)*1e3, psi_probe1=0.97000*np.pi, h_EI_probe1=150e3,\n", " rp_probe2=(6051.8+103.85)*1e3, psi_probe2=0.94000*np.pi, h_EI_probe2=150e3,\n", " r_dv_rp=138.5)\n", "\n", "deflection5 = ProbeOrbiterDeflection(arrivalPlanet=\"VENUS\",\n", " v_inf_vec_icrf_kms=np.array([-3.26907094, 0.67649494, -1.0697747 ]),\n", " rp_probe=(6051.8+103.85)*1e3, psi_probe=0.94000*np.pi, h_EI_probe=150e3,\n", " rp_space=(6051.8+400)*1e3, psi_space=1.0*np.pi,\n", " r_dv_rp=112.5)" ] }, { "cell_type": "code", "execution_count": 3, "id": "a0e24b36", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Deflection Maneuver #1\n", "Divert manuever DV: [-0.75224316 -4.61251793 0.6152991 ] m/s\n", "Divert manuever DV magnitude: 4.714 m/s\n", "---------------------------------\n", "Deflection Maneuver #2\n", "Divert manuever DV: [-0.67485826 -5.31732722 0.23444224] m/s\n", "Divert manuever DV magnitude: 5.365 m/s\n", "---------------------------------\n", "Deflection Maneuver #3\n", "Divert manuever DV: [-0.56133357 -6.15017052 -0.27012275] m/s\n", "Divert manuever DV magnitude: 6.182 m/s\n", "---------------------------------\n", "Deflection Maneuver #4\n", "Divert manuever DV: [-0.40090375 -7.21278622 -0.9523719 ] m/s\n", "Divert manuever DV magnitude: 7.286 m/s\n", "---------------------------------\n", "Deflection Maneuver #5\n", "Divert manuever DV: [ 0.35802126 17.78715771 4.02682792] m/s\n", "Divert manuever DV magnitude: 18.241 m/s\n", "---------------------------------\n" ] } ], "source": [ "print('Deflection Maneuver #1')\n", "print(\"Divert manuever DV: \"+str(deflection1.dv_maneuver_vec)+ \" m/s\")\n", "print(\"Divert manuever DV magnitude: \"+str(round(deflection1.dv_maneuver_mag,3))+ \" m/s\")\n", "print('---------------------------------')\n", "print('Deflection Maneuver #2')\n", "print(\"Divert manuever DV: \"+str(deflection2.dv_maneuver_vec)+ \" m/s\")\n", "print(\"Divert manuever DV magnitude: \"+str(round(deflection2.dv_maneuver_mag,3))+ \" m/s\")\n", "print('---------------------------------')\n", "print('Deflection Maneuver #3')\n", "print(\"Divert manuever DV: \"+str(deflection3.dv_maneuver_vec)+ \" m/s\")\n", "print(\"Divert manuever DV magnitude: \"+str(round(deflection3.dv_maneuver_mag,3))+ \" m/s\")\n", "print('---------------------------------')\n", "print('Deflection Maneuver #4')\n", "print(\"Divert manuever DV: \"+str(deflection4.dv_maneuver_vec)+ \" m/s\")\n", "print(\"Divert manuever DV magnitude: \"+str(round(deflection4.dv_maneuver_mag,3))+ \" m/s\")\n", "print('---------------------------------')\n", "print('Deflection Maneuver #5')\n", "print(\"Divert manuever DV: \"+str(deflection5.dv_maneuver_vec)+ \" m/s\")\n", "print(\"Divert manuever DV magnitude: \"+str(round(deflection5.dv_maneuver_mag,3))+ \" m/s\")\n", "print('---------------------------------')" ] }, { "cell_type": "code", "execution_count": 4, "id": "010582f4", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "TOF from probe release to atm. entry interface: 4.007605 days\n", "TOF from probe release to atm. entry interface: 3.500867 days\n", "TOF from probe release to atm. entry interface: 3.001488 days\n", "TOF from probe release to atm. entry interface: 2.506011 days\n", "TOF from probe release to atm. entry interface: 2.003604 days\n" ] } ], "source": [ "print(\"TOF from probe release to atm. entry interface: \"+str(round(deflection1.TOF_probe1,6))+ \" days\")\n", "print(\"TOF from probe release to atm. entry interface: \"+str(round(deflection2.TOF_probe1,6))+ \" days\")\n", "print(\"TOF from probe release to atm. entry interface: \"+str(round(deflection3.TOF_probe1,6))+ \" days\")\n", "print(\"TOF from probe release to atm. entry interface: \"+str(round(deflection4.TOF_probe1,6))+ \" days\")\n", "print(\"TOF from probe release to atm. entry interface: \"+str(round(deflection5.TOF_probe,6))+ \" days\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.16" } }, "nbformat": 4, "nbformat_minor": 5 }