def __init__(self, *args, **kwargs): #Constants used Cts = AL_BF.ConstantsBook() #Constants for planets # Define bodies involved self.Spacecraft = AL_2BP.Spacecraft() self.earthephem = pk.planet.jpl_lp('earth') self.marsephem = pk.planet.jpl_lp('mars') self.jupiterephem = pk.planet.jpl_lp('jupiter') self.venusephem = pk.planet.jpl_lp('venus') self.sun = AL_2BP.Body('Sun', 'yellow', mu=Cts.mu_S_m) self.earth = AL_2BP.Body('Earth', 'blue', mu=Cts.mu_E_m) self.mars = AL_2BP.Body('Mars', 'red', mu=Cts.mu_M_m) self.jupiter = AL_2BP.Body('Jupiter', 'green', mu=Cts.mu_J_m) self.venus = AL_2BP.Body('Venus', 'brown', mu=Cts.mu_V_m) departure = kwargs.get('departure', "Earth") arrival = kwargs.get('arrival', "Mars") if departure == 'Earth': self.departure = self.earth self.departureephem = self.earthephem elif departure == 'Mars': self.departure = self.mars self.departureephem = self.marsephem if arrival == 'Mars': self.arrival = self.mars self.arrivalephem = self.marsephem elif arrival == 'Jupiter': self.arrival = self.jupiter self.arrivalephem = self.jupiterephem elif arrival == 'Venus': self.arrival = self.venus self.arrivalephem = self.venusephem elif arrival == 'Earth': self.arrival = self.earth self.arrivalephem = self.earthephem # Settings self.color = ['blue', 'green', 'black', 'red', 'yellow', 'orange'] # Choose odd number to compare easily in middle point self.Nimp = kwargs.get('Nimp', 11)
from FitnessFunction import Fitness import AstroLibraries.AstroLib_Basic as AL_BF from AstroLibraries import AstroLib_2BP as AL_2BP from AstroLibraries import AstroLib_Ephem as AL_Eph from AstroLibraries import AstroLib_Plots as AL_Plot from AstroLibraries import AstroLib_Trajectories as AL_TR from AstroLibraries import AstroLib_ShapingMethod as AL_Sh import LoadConfigFiles as CONFIG # # In[] # ############################################## # ############# VALIDATION ##################### # ############################################## Cts = AL_BF.ConstantsBook() def pykepephem(): # earth = pk.planet.jpl_lp('earth') # marsephem = pk.planet.jpl_lp('mars') earth = pk.planet.jpl_lp('earth') K_0 = earth.eph(0) # print(K_0) earth.mu_central_body print(earth.mu_self) elem = earth.osculating_elements(pk.epoch(2345.3, 'mjd2000')) print(elem) print(np.array(elem))