Commit 681de231 authored by Frederik Zahle's avatar Frederik Zahle
Browse files

fixing PMtarget in old OpenMDAO model

parent 5cb42725
......@@ -23,7 +23,7 @@ class SEAMTower(Component):
stress_limit_extreme_tower = Float(iotype='in', units='MPa', desc='Tower ultimate strength')
stress_limit_fatigue_tower = Float(iotype='in', units='MPa', desc='Tower fatigue strength')
safety_factor_tower = Float(iotype='in', desc='Tower loads safety factor')
PMtarget = Float(1., iotype='in', desc='')
PMtarget_tower = Float(1., iotype='in', desc='')
wohler_exponent_tower = Float(4., iotype='in', desc='Tower fatigue Wohler exponent')
tower_z = Array(iotype='out', desc='Tower discretization')
......@@ -73,7 +73,7 @@ class SEAMTower(Component):
W = (np.pi / 32.) * (D**4 - (D - 2. * t)**4) / D
sfat = self.safety_factor_tower * 1.e3 * Meq/ W / 1.e6
PM = self.lifetime_cycles / (pow(10, (self.stress_limit_fatigue_tower - self.wohler_exponent_tower*np.log10(sfat))))
return abs(PM - self.PMtarget)
return abs(PM - self.PMtarget_tower)
def plot(self, fig):
"""
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment