class Test_HubWPACT(unittest.TestCase):
def setUp(self):
self.hub = HubWPACT()
self.hub.blade_mass = 17740.0 # kg
self.hub.rotor_diameter = 126.0 # m
self.hub.blade_number = 3
self.hub.blade_root_diameter = 3.542
AirDensity= 1.225 # kg/(m^3)
Solidity = 0.0517
RatedWindSpeed = 11.05 # m/s
self.hub.rotor_bending_moment = (3.06 * pi / 8) * AirDensity * (RatedWindSpeed ** 2) * (Solidity * (self.hub.rotor_diameter ** 3)) / self.hub.blade_number
def test_functionality(self):
self.hub.run()
self.assertEqual(round(self.hub.hub_system_mass,1), 44709.1)
class Test_HubWPACT(unittest.TestCase):
def setUp(self):
self.hub = HubWPACT()
self.hub.blade_mass = 17740.0 # kg
self.hub.rotor_diameter = 126.0 # m
self.hub.blade_number = 3
self.hub.blade_root_diameter = 3.542
AirDensity = 1.225 # kg/(m^3)
Solidity = 0.0517
RatedWindSpeed = 11.05 # m/s
self.hub.rotor_bending_moment = (3.06 * pi / 8) * AirDensity * (
RatedWindSpeed**
2) * (Solidity *
(self.hub.rotor_diameter**3)) / self.hub.blade_number
def test_functionality(self):
self.hub.run()
self.assertEqual(round(self.hub.hub_system_mass, 1), 44709.1)
hub = HubWPACT()
hub.blade_mass = 17740.0 # kg
hub.rotor_diameter = 126.0 # m
hub.blade_number = 3
hub.blade_root_diameter = 3.542
AirDensity = 1.225 # kg/(m^3)
Solidity = 0.0517
RatedWindSpeed = 11.05 # m/s
hub.rotor_bending_moment = (3.06 * pi / 8) * AirDensity * (
RatedWindSpeed**2) * (Solidity *
(hub.rotor_diameter**3)) / hub.blade_number
# 2 ----------
# 3 ----------
hub.run()
# 3 ----------
# 4 ----------
print "Estimate of Hub Component Sizes for the NREL 5 MW Reference Turbine"
print "Hub Components"
print ' Hub: {0:8.1f} kg'.format(hub.hub.mass) # 31644.47
print ' Pitch system: {0:8.1f} kg'.format(hub.pitchSystem.mass) # 17003.98
print ' Nose cone: {0:8.1f} kg'.format(hub.spinner.mass) # 1810.50
print 'Hub system total: {0:8.1f} kg'.format(hub.hub_system_mass) # 50458.95
print ' cm {0:6.2f} {1:6.2f} {2:6.2f} [m, m, m]'.format(
hub.hub_system_cm[0], hub.hub_system_cm[1], hub.hub_system_cm[2])
print ' I {0:6.1f} {1:6.1f} {2:6.1f} [kg*m^2, kg*m^2, kg*m^2]'.format(
hub.hub_system_I[0], hub.hub_system_I[1], hub.hub_system_I[2])
print
# NREL 5 MW turbine
hub = HubWPACT()
hub.blade_mass = 17740.0 # kg
hub.rotor_diameter = 126.0 # m
hub.blade_number = 3
hub.blade_root_diameter = 3.542
AirDensity= 1.225 # kg/(m^3)
Solidity = 0.0517
RatedWindSpeed = 11.05 # m/s
hub.rotor_bending_moment = (3.06 * pi / 8) * AirDensity * (RatedWindSpeed ** 2) * (Solidity * (hub.rotor_diameter ** 3)) / hub.blade_number
# 2 ----------
# 3 ----------
hub.run()
# 3 ----------
# 4 ----------
print "Estimate of Hub Component Sizes for the NREL 5 MW Reference Turbine"
print "Hub Components"
print ' Hub: {0:8.1f} kg'.format(hub.hub.mass) # 31644.47
print ' Pitch system: {0:8.1f} kg'.format(hub.pitchSystem.mass) # 17003.98
print ' Nose cone: {0:8.1f} kg'.format(hub.spinner.mass) # 1810.50
print 'Hub system total: {0:8.1f} kg'.format(hub.hub_system_mass) # 50458.95
print ' cm {0:6.2f} {1:6.2f} {2:6.2f} [m, m, m]'.format(hub.hub_system_cm[0], hub.hub_system_cm[1], hub.hub_system_cm[2])
print ' I {0:6.1f} {1:6.1f} {2:6.1f} [kg*m^2, kg*m^2, kg*m^2]'.format(hub.hub_system_I[0], hub.hub_system_I[1], hub.hub_system_I[2])
print
# 4 ---------
