diff --git a/Model/Propulsion/__pycache__/power_comp.cpython-36.pyc b/Model/Propulsion/__pycache__/power_comp.cpython-36.pyc index 98ccae6..663225d 100644 Binary files a/Model/Propulsion/__pycache__/power_comp.cpython-36.pyc and b/Model/Propulsion/__pycache__/power_comp.cpython-36.pyc differ diff --git a/Model/Propulsion/power_comp.py b/Model/Propulsion/power_comp.py index 4b96418..0d83a25 100644 --- a/Model/Propulsion/power_comp.py +++ b/Model/Propulsion/power_comp.py @@ -1,9 +1,9 @@ from openmdao.api import ExplicitComponent +n = 1 # variable accounting for failed motor and e extra power class PowerComp(ExplicitComponent): - - + def setup(self): @@ -30,7 +30,7 @@ def compute(self, inputs, outputs): S = 3.14*(r**2) - T = W/8 # thrust required per motor + T = n*W/8 # thrust required per motor ct = T/( rho *S* (tip_speed**2) ) cp = ( 1.15*( (ct**(3/2)) /(2**.5) ) + (1/8)*0.1*cd0 ) PI = T * ((T/(2*S*rho))**.5) diff --git a/solver.py b/solver.py index 7d764ff..b9c98f7 100644 --- a/solver.py +++ b/solver.py @@ -42,8 +42,8 @@ prob.driver.options['disp'] = True model.add_design_var('r',lower=0.5,upper=1.5) -model.add_design_var('We/W0',lower=0.50,upper=0.80) -model.add_constraint('weight.W0',lower = 2000,upper=2500) +model.add_design_var('We/W0',lower=0.30,upper=0.70) # 30% - 70%, from lecture +model.add_constraint('weight.W0',lower = 1500,upper=2500) model.add_constraint('FOM.FM',equals=0.80) model.add_objective('FOM.PH')