def tier2Averaged():
baseCost = 290000000
# Method Definitons
IBS = removal.orbitalRemovalMethod("IBS", (3000, 5000), (300, 1000), 1.71, baseCost + 20000000)
robotArmPropellant = removal.orbitalRemovalMethod("Robot Arm -- Propellant", (65, 3800), (750, 900), 5, baseCost + 10000000)
robotArmEDT = removal.orbitalRemovalMethod("Robot Arm -- EDT", (500, 3400), (800, 1400), 1, baseCost)
EDDE = removal.orbitalRemovalMethod("EDDE", (2, 50), (800, 1000), 357, baseCost - 30000000)
# Data Processing
orbitalObjects = dataProcessing.junkList("spaceobjects.csv")
# Base ROI Lists
ROIlistIBS = []
ROIlistPropellant = []
ROIlistEDT = []
ROIlistEDDE = []
MissionRisklistIBS = []
MissionRisklistPropellant = []
MissionRisklistEDT = []
MissionRisklistEDDE = []
# Sort List
orbitalObjects.sort(key=operator.attrgetter('weightedRisk'))
orbitalObjects = risk.normalizeWeightedRisk(orbitalObjects)
orbitalObjects.reverse()
for i in range(0, 1000):
years = 5
IBSresults = removal.debrisRemoval(IBS, orbitalObjects, years)
Propellantresults = removal.debrisRemoval(robotArmPropellant, orbitalObjects, years)
EDTresults = removal.debrisRemoval(robotArmEDT, orbitalObjects, years)
EDDEresults = removal.debrisRemoval(EDDE, orbitalObjects, years)
ROIlistIBS.append(cost.ROI(IBSresults, IBS))
ROIlistPropellant.append(cost.ROI(Propellantresults, robotArmPropellant))
ROIlistEDT.append(cost.ROI(EDTresults, robotArmEDT))
ROIlistEDDE.append(cost.ROI(EDDEresults, EDDE))
MissionRisklistIBS.append(risk.missionRisk(IBSresults))
MissionRisklistPropellant.append(risk.missionRisk(Propellantresults))
MissionRisklistEDT.append(risk.missionRisk(EDTresults))
MissionRisklistEDDE.append(risk.missionRisk(EDDEresults))
import numpy as np
ROIaverageIBS = np.mean(ROIlistIBS)
ROIaveragePropellant = np.mean(ROIlistPropellant)
ROIaverageEDT = np.mean(ROIlistEDT)
ROIaverageEDDE = np.mean(ROIlistEDDE)
MissionRiskIBS = np.mean(MissionRisklistIBS)
MissionRiskPropellant = np.mean(MissionRisklistPropellant)
MissionRiskEDT = np.mean(MissionRisklistEDT)
MissionRiskEDDE = np.mean(MissionRisklistEDDE)
print "ROI Average IBS: " + str(ROIaverageIBS) + " Mission Risk: " + str(MissionRiskIBS)
print "ROI Average Propellant: " + str(ROIaveragePropellant) + " Mission Risk: " + str(MissionRiskPropellant)
print "ROI Average EDT: " + str(ROIaverageEDT) + " Mission Risk: " + str(MissionRiskEDT)
print "ROI Average EDDE: " + str(ROIaverageEDDE) + " Mission Risk: " + str(MissionRiskEDDE)
def tier3():
baseCost = 290000000
IBS = removal.orbitalRemovalMethod("IBS", (3000, 5000), (300, 1000), 1.71, baseCost + 20000000)
robotArmPropellant = removal.orbitalRemovalMethod("Robot Arm -- Propellant", (65, 3800), (750, 900), 5, baseCost + 10000000)
robotArmEDT = removal.orbitalRemovalMethod("Robot Arm -- EDT", (500, 3400), (800, 1400), 1, baseCost)
EDDE = removal.orbitalRemovalMethod("EDDE", (2, 50), (800, 1000), 357, baseCost - 30000000)
orbitalObjects = dataProcessing.junkList("spaceobjects.csv")
# Base ROI Lists
ROIlistIBS = []
ROIlistPropellant = []
ROIlistEDT = []
ROIlistEDDE = []
# Combination ROI Lists
ROIlistIBS_Propellant = []
ROIlistIBS_EDT = []
ROIlistIBS_EDDE = []
ROIlistPropellant_EDT = []
ROIlistPropellant_EDDE = []
ROIlistEDT_EDDE = []
MissionRisklistIBS = []
MissionRisklistPropellant = []
MissionRisklistEDT = []
MissionRisklistEDDE = []
MissionRisklistIBS_Propellant = []
MissionRisklistIBS_EDT = []
MissionRisklistIBS_EDDE = []
MissionRisklistPropellant_EDT = []
MissionRisklistPropellant_EDDE = []
MissionRisklistEDT_EDDE = []
for i in range(0, 1):
# Sort List
orbitalObjects.sort(key=operator.attrgetter('weightedRisk'))
orbitalObjects = risk.normalizeWeightedRisk(orbitalObjects)
orbitalObjects.reverse()
# Risk Analysis
for years in range(0, 200):
IBSresults = removal.debrisRemoval(IBS, orbitalObjects, years)
Propellantresults = removal.debrisRemoval(robotArmPropellant, orbitalObjects, years)
EDTresults = removal.debrisRemoval(robotArmEDT, orbitalObjects, years)
EDDEresults = removal.debrisRemoval(EDDE, orbitalObjects, years)
IBS_Propellanresults = removal.debrisRemoval(robotArmPropellant, IBSresults, years)
IBS_EDTresults = removal.debrisRemoval(robotArmEDT, IBSresults, years)
IBS_EDDEresults = removal.debrisRemoval(EDDE, IBSresults, years)
Propellant_EDTresults = removal.debrisRemoval(robotArmEDT, Propellantresults, years)
Propellant_EDDEresults = removal.debrisRemoval(EDDE, Propellantresults, years)
EDT_EDDEresults = removal.debrisRemoval(EDDE, EDTresults, years)
ROIlistIBS.append(cost.ROI(IBSresults, IBS))
ROIlistPropellant.append(cost.ROI(Propellantresults, robotArmPropellant))
ROIlistEDT.append(cost.ROI(EDTresults, robotArmEDT))
ROIlistEDDE.append(cost.ROI(EDDEresults, EDDE))
ROIlistIBS_Propellant.append(cost.ROI(IBS_Propellanresults, IBS, robotArmPropellant))
ROIlistIBS_EDT.append(cost.ROI(IBS_EDTresults, IBS, robotArmEDT))
ROIlistIBS_EDDE.append(cost.ROI(IBS_EDDEresults, IBS, EDDE))
ROIlistPropellant_EDT.append(cost.ROI(Propellant_EDTresults, robotArmPropellant, robotArmEDT))
ROIlistPropellant_EDDE.append(cost.ROI(Propellant_EDDEresults, robotArmPropellant, EDDE))
ROIlistEDT_EDDE.append(cost.ROI(EDT_EDDEresults, robotArmEDT, EDDE))
MissionRisklistIBS.append(risk.missionRisk(IBSresults))
MissionRisklistPropellant.append(risk.missionRisk(Propellantresults))
MissionRisklistEDT.append(risk.missionRisk(EDTresults))
MissionRisklistEDDE.append(risk.missionRisk(EDDEresults))
MissionRisklistIBS_Propellant.append(risk.missionRisk(IBS_Propellanresults))
MissionRisklistIBS_EDT.append(risk.missionRisk(IBS_EDTresults))
MissionRisklistIBS_EDDE.append(risk.missionRisk(IBS_EDDEresults))
MissionRisklistPropellant_EDT.append(risk.missionRisk(Propellant_EDTresults))
MissionRisklistPropellant_EDDE.append(risk.missionRisk(Propellant_EDDEresults))
MissionRisklistEDT_EDDE.append(risk.missionRisk(EDT_EDDEresults))
for debris in orbitalObjects:
debris.scramble()
print ROIlistIBS
print ROIlistPropellant
print ROIlistEDT
print ROIlistEDDE
print ROIlistIBS_Propellant
print ROIlistIBS_EDT
print ROIlistIBS_EDDE
print ROIlistPropellant_EDT
print ROIlistPropellant_EDDE
print ROIlistEDT_EDDE
print MissionRisklistIBS
print MissionRisklistPropellant
print MissionRisklistEDT
print MissionRisklistEDDE
print MissionRisklistIBS_Propellant
print MissionRisklistIBS_EDT
print MissionRisklistIBS_EDDE
print MissionRisklistPropellant_EDT
print MissionRisklistPropellant_EDDE
print MissionRisklistEDT_EDDE
def tier3():
baseCost = 290000000
IBS = removal.orbitalRemovalMethod("IBS", (3000, 5000), (300, 1000), 1.71,
baseCost + 20000000)
robotArmPropellant = removal.orbitalRemovalMethod(
"Robot Arm -- Propellant", (65, 3800), (750, 900), 5,
baseCost + 10000000)
robotArmEDT = removal.orbitalRemovalMethod("Robot Arm -- EDT", (500, 3400),
(800, 1400), 1, baseCost)
EDDE = removal.orbitalRemovalMethod("EDDE", (2, 50), (800, 1000), 357,
baseCost - 30000000)
orbitalObjects = dataProcessing.junkList("spaceobjects.csv")
# Base ROI Lists
ROIlistIBS = []
ROIlistPropellant = []
ROIlistEDT = []
ROIlistEDDE = []
# Combination ROI Lists
ROIlistIBS_Propellant = []
ROIlistIBS_EDT = []
ROIlistIBS_EDDE = []
ROIlistPropellant_EDT = []
ROIlistPropellant_EDDE = []
ROIlistEDT_EDDE = []
MissionRisklistIBS = []
MissionRisklistPropellant = []
MissionRisklistEDT = []
MissionRisklistEDDE = []
MissionRisklistIBS_Propellant = []
MissionRisklistIBS_EDT = []
MissionRisklistIBS_EDDE = []
MissionRisklistPropellant_EDT = []
MissionRisklistPropellant_EDDE = []
MissionRisklistEDT_EDDE = []
for i in range(0, 1):
# Sort List
orbitalObjects.sort(key=operator.attrgetter('weightedRisk'))
orbitalObjects = risk.normalizeWeightedRisk(orbitalObjects)
orbitalObjects.reverse()
# Risk Analysis
for years in range(0, 200):
IBSresults = removal.debrisRemoval(IBS, orbitalObjects, years)
Propellantresults = removal.debrisRemoval(robotArmPropellant,
orbitalObjects, years)
EDTresults = removal.debrisRemoval(robotArmEDT, orbitalObjects,
years)
EDDEresults = removal.debrisRemoval(EDDE, orbitalObjects, years)
IBS_Propellanresults = removal.debrisRemoval(
robotArmPropellant, IBSresults, years)
IBS_EDTresults = removal.debrisRemoval(robotArmEDT, IBSresults,
years)
IBS_EDDEresults = removal.debrisRemoval(EDDE, IBSresults, years)
Propellant_EDTresults = removal.debrisRemoval(
robotArmEDT, Propellantresults, years)
Propellant_EDDEresults = removal.debrisRemoval(
EDDE, Propellantresults, years)
EDT_EDDEresults = removal.debrisRemoval(EDDE, EDTresults, years)
ROIlistIBS.append(cost.ROI(IBSresults, IBS))
ROIlistPropellant.append(
cost.ROI(Propellantresults, robotArmPropellant))
ROIlistEDT.append(cost.ROI(EDTresults, robotArmEDT))
ROIlistEDDE.append(cost.ROI(EDDEresults, EDDE))
ROIlistIBS_Propellant.append(
cost.ROI(IBS_Propellanresults, IBS, robotArmPropellant))
ROIlistIBS_EDT.append(cost.ROI(IBS_EDTresults, IBS, robotArmEDT))
ROIlistIBS_EDDE.append(cost.ROI(IBS_EDDEresults, IBS, EDDE))
ROIlistPropellant_EDT.append(
cost.ROI(Propellant_EDTresults, robotArmPropellant,
robotArmEDT))
ROIlistPropellant_EDDE.append(
cost.ROI(Propellant_EDDEresults, robotArmPropellant, EDDE))
ROIlistEDT_EDDE.append(cost.ROI(EDT_EDDEresults, robotArmEDT,
EDDE))
MissionRisklistIBS.append(risk.missionRisk(IBSresults))
MissionRisklistPropellant.append(
risk.missionRisk(Propellantresults))
MissionRisklistEDT.append(risk.missionRisk(EDTresults))
MissionRisklistEDDE.append(risk.missionRisk(EDDEresults))
MissionRisklistIBS_Propellant.append(
risk.missionRisk(IBS_Propellanresults))
MissionRisklistIBS_EDT.append(risk.missionRisk(IBS_EDTresults))
MissionRisklistIBS_EDDE.append(risk.missionRisk(IBS_EDDEresults))
MissionRisklistPropellant_EDT.append(
risk.missionRisk(Propellant_EDTresults))
MissionRisklistPropellant_EDDE.append(
risk.missionRisk(Propellant_EDDEresults))
MissionRisklistEDT_EDDE.append(risk.missionRisk(EDT_EDDEresults))
for debris in orbitalObjects:
debris.scramble()
print ROIlistIBS
print ROIlistPropellant
print ROIlistEDT
print ROIlistEDDE
print ROIlistIBS_Propellant
print ROIlistIBS_EDT
print ROIlistIBS_EDDE
print ROIlistPropellant_EDT
print ROIlistPropellant_EDDE
print ROIlistEDT_EDDE
print MissionRisklistIBS
print MissionRisklistPropellant
print MissionRisklistEDT
print MissionRisklistEDDE
print MissionRisklistIBS_Propellant
print MissionRisklistIBS_EDT
print MissionRisklistIBS_EDDE
print MissionRisklistPropellant_EDT
print MissionRisklistPropellant_EDDE
print MissionRisklistEDT_EDDE
def tier2Averaged():
baseCost = 290000000
# Method Definitons
IBS = removal.orbitalRemovalMethod("IBS", (3000, 5000), (300, 1000), 1.71,
baseCost + 20000000)
robotArmPropellant = removal.orbitalRemovalMethod(
"Robot Arm -- Propellant", (65, 3800), (750, 900), 5,
baseCost + 10000000)
robotArmEDT = removal.orbitalRemovalMethod("Robot Arm -- EDT", (500, 3400),
(800, 1400), 1, baseCost)
EDDE = removal.orbitalRemovalMethod("EDDE", (2, 50), (800, 1000), 357,
baseCost - 30000000)
# Data Processing
orbitalObjects = dataProcessing.junkList("spaceobjects.csv")
# Base ROI Lists
ROIlistIBS = []
ROIlistPropellant = []
ROIlistEDT = []
ROIlistEDDE = []
MissionRisklistIBS = []
MissionRisklistPropellant = []
MissionRisklistEDT = []
MissionRisklistEDDE = []
# Sort List
orbitalObjects.sort(key=operator.attrgetter('weightedRisk'))
orbitalObjects = risk.normalizeWeightedRisk(orbitalObjects)
orbitalObjects.reverse()
for i in range(0, 1000):
years = 5
IBSresults = removal.debrisRemoval(IBS, orbitalObjects, years)
Propellantresults = removal.debrisRemoval(robotArmPropellant,
orbitalObjects, years)
EDTresults = removal.debrisRemoval(robotArmEDT, orbitalObjects, years)
EDDEresults = removal.debrisRemoval(EDDE, orbitalObjects, years)
ROIlistIBS.append(cost.ROI(IBSresults, IBS))
ROIlistPropellant.append(
cost.ROI(Propellantresults, robotArmPropellant))
ROIlistEDT.append(cost.ROI(EDTresults, robotArmEDT))
ROIlistEDDE.append(cost.ROI(EDDEresults, EDDE))
MissionRisklistIBS.append(risk.missionRisk(IBSresults))
MissionRisklistPropellant.append(risk.missionRisk(Propellantresults))
MissionRisklistEDT.append(risk.missionRisk(EDTresults))
MissionRisklistEDDE.append(risk.missionRisk(EDDEresults))
import numpy as np
ROIaverageIBS = np.mean(ROIlistIBS)
ROIaveragePropellant = np.mean(ROIlistPropellant)
ROIaverageEDT = np.mean(ROIlistEDT)
ROIaverageEDDE = np.mean(ROIlistEDDE)
MissionRiskIBS = np.mean(MissionRisklistIBS)
MissionRiskPropellant = np.mean(MissionRisklistPropellant)
MissionRiskEDT = np.mean(MissionRisklistEDT)
MissionRiskEDDE = np.mean(MissionRisklistEDDE)
print "ROI Average IBS: " + str(
ROIaverageIBS) + " Mission Risk: " + str(MissionRiskIBS)
print "ROI Average Propellant: " + str(
ROIaveragePropellant) + " Mission Risk: " + str(MissionRiskPropellant)
print "ROI Average EDT: " + str(
ROIaverageEDT) + " Mission Risk: " + str(MissionRiskEDT)
print "ROI Average EDDE: " + str(
ROIaverageEDDE) + " Mission Risk: " + str(MissionRiskEDDE)