The Gitlab server is succesfully updated to version 14.4.0

Many users are receiving emails regarding excessive amounts of log-in attempts, which are brute-force attempts to crack your password. We are working on a 2FA solution, so no action is needed yet. If you receive emails regarding "sign-in from new location", please check the IP address and if unknown, notify frza@dtu.dk immediately.

Commit 0d5bcb50 authored by Frederik Zahle's avatar Frederik Zahle
Browse files

fixing PMtarget in old OpenMDAO model

parent 16d90cea
......@@ -16,7 +16,7 @@ class SEAMBladeStructure(Component):
lifetime_cycles = Float(1.e7, iotype='in', desc='Equivalent lifetime cycles')
wohler_exponent_blade_flap = Float(iotype='in', desc='Blade flap fatigue Wohler exponent')
PMtarget = Float(iotype='in', desc='')
PMtarget_blades = Float(1., iotype='in', desc='')
rotor_diameter = Float(iotype='in', units='m', desc='Rotor diameter')
MaxChordrR = Float(iotype='in', units='m', desc='Spanwise position of maximum chord')
......@@ -189,7 +189,7 @@ class SEAMBladeStructure(Component):
W = Ine/(thick/2.)
sfat = np.maximum(1.e-6, self.safety_factor_blade*1.e3*Mfat_flap/ W / 1.e6)
PM = self.lifetime_cycles/(pow(10, (self.stress_limit_fatigue_blade - self.wohler_exponent_blade_flap*np.log10(sfat))))
return abs(PM - self.PMtarget) / norm
return abs(PM - self.PMtarget_blades) / norm
# Solving for t in edge direction, fatigue
def solve_tfat_edge(self, t, C, thick, Mfat_edge, norm):
......@@ -197,7 +197,7 @@ class SEAMBladeStructure(Component):
W = Ine/(C/2.)
sfat = np.maximum(1.e-6, self.safety_factor_blade*1.e3*Mfat_edge/W/1.e6)
PM = self.lifetime_cycles/(pow(10, (self.stress_limit_fatigue_blade - self.wohler_exponent_blade_flap*np.log10(sfat))))
return abs(PM - self.PMtarget) / norm
return abs(PM - self.PMtarget_blades) / norm
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