Commit 06c8cb88 authored by Frederik Zahle's avatar Frederik Zahle
Browse files

changed another varname

parent cbad7f60
......@@ -36,8 +36,8 @@ class SEAMBladeStructure(Component):
sc_frac_edge = Float(iotype='in', desc='spar cap fraction of thickness')
SF_blade = Float(iotype='in', desc='')
Slim_ext_blade = Float(iotype='in', units='MPa', desc='Ultimate strength')
Slim_fat_blade = Float(iotype='in', units='MPa', desc='Fatigue strength')
stress_limit_extreme_blade = Float(iotype='in', units='MPa', desc='Ultimate strength')
stress_limit_fatigue_blade = Float(iotype='in', units='MPa', desc='Fatigue strength')
AddWeightFactorBlade = Float(iotype='in', desc='Additional weight factor for blade shell')
blade_material_density = Float(iotype='in', units='kg/m**3', desc='Average density of blade materials')
......@@ -175,22 +175,22 @@ class SEAMBladeStructure(Component):
Ine = (2./3.)*self.sc_frac_flap*C*t**3-self.sc_frac_flap*C*thick*t**2+self.sc_frac_flap*C*thick**2/2.*t
W = Ine/(thick/2.)
sext = self.SF_blade*1e3*Mext_flap/W/1e6
# print 'text_flap', abs(sext - self.Slim_ext_blade)
return abs(sext - self.Slim_ext_blade) / norm
# print 'text_flap', abs(sext - self.stress_limit_extreme_blade)
return abs(sext - self.stress_limit_extreme_blade) / norm
#Solving for t in edge direction, extremes
def solve_text_edge(self, t, C, thick, Mext_edge, norm):
Ine = (2/3.)*self.sc_frac_edge*thick*t**3-self.sc_frac_edge*thick*C*t**2+self.sc_frac_edge*thick*C**2/2.*t
W = Ine/(C/2.);
sext = self.SF_blade*1.e3*Mext_edge/W/1.e6
return abs(sext - self.Slim_ext_blade) / norm
return abs(sext - self.stress_limit_extreme_blade) / norm
# Solving for t in flap direction, fatigue
def solve_tfat_flap(self, t, C, thick, Mfat_flap, norm):
Ine = (2/3.)*self.sc_frac_flap*C*t**3-self.sc_frac_flap*C*thick*t**2+self.sc_frac_flap*C*thick**2/2.*t
W = Ine/(thick/2.)
sfat = np.maximum(1.e-6, self.SF_blade*1.e3*Mfat_flap/ W / 1.e6)
PM = self.lifetime_cycles/(pow(10, (self.Slim_fat_blade - self.wohler_exponent_blade_flap*np.log10(sfat))))
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
# Solving for t in edge direction, fatigue
......@@ -198,7 +198,7 @@ class SEAMBladeStructure(Component):
Ine = (2/3.)*self.sc_frac_edge*thick*t**3-self.sc_frac_edge*thick*C*t**2+self.sc_frac_edge*thick*C**2/2.*t
W = Ine/(C/2.)
sfat = np.maximum(1.e-6, self.SF_blade*1.e3*Mfat_edge/W/1.e6)
PM = self.lifetime_cycles/(pow(10, (self.Slim_fat_blade - self.wohler_exponent_blade_flap*np.log10(sfat))))
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
def plot(self, fig):
......
......@@ -25,8 +25,8 @@ def configure():
top.sc_frac_edge = 0.8 # sparcap fraction of thickness edge
top.SF_blade = 1.1 #[factor]
top.Slim_ext_blade = 200.0
top.Slim_fat_blade = 27.
top.stress_limit_extreme_blade = 200.0
top.stress_limit_fatigue_blade = 27.
top.AddWeightFactorBlade = 1.2 # Additional weight factor for blade shell
top.blade_material_density = 2100. # [kg / m^3]
......
Supports Markdown
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