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 4631f16a authored by Frederik Zahle's avatar Frederik Zahle
Browse files

use np.trapz to integrate volume of blade

parent c6192b08
......@@ -18,32 +18,29 @@ class SEAMBladeStructure(Component):
WohlerExpFlap = Float(iotype='in', desc='')
PMtarget = Float(iotype='in', desc='')
rotor_diameter = Float(iotype='in', units='m', desc='') #[m]
# RootChord = Float(iotype='in', units='m', desc='') #[m]
# MaxChord = Float(iotype='in', units='m', desc='') #[m]
MaxChordrR = Float(iotype='in', units='m', desc='') #[m]
rotor_diameter = Float(iotype='in', units='m', desc='Rotor diameter')
# RootChord = Float(iotype='in', units='m', desc='')
# MaxChord = Float(iotype='in', units='m', desc='')
MaxChordrR = Float(iotype='in', units='m', desc='Spanwise position of maximum chord')
overallMaxFlap = Float(iotype='in', desc='')
overallMaxEdge = Float(iotype='in', desc='')
TIF_FLext = Float(iotype='in', desc='') # Tech Impr Factor _ flap extreme
TIF_EDext = Float(iotype='in', desc='')
overallMaxFlap = Float(iotype='in', desc='Max blade root flap moment')
overallMaxEdge = Float(iotype='in', desc='Max blade root edge moment')
TIF_FLext = Float(iotype='in', desc='Technology improvement factor flap extreme')
TIF_EDext = Float(iotype='in', desc='Technology improvement factor edge extreme')
FlapLEQ = Float(iotype='in', desc='')
EdgeLEQ = Float(iotype='in', desc='')
TIF_FLfat = Float(iotype='in', desc='')
FlapLEQ = Float(iotype='in', desc='Blade root flap lifetime eq. moment')
EdgeLEQ = Float(iotype='in', desc='Blade root edge lifetime eq. moment')
TIF_FLfat = Float(iotype='in', desc='Technology improvement factor FlapLEQ')
sc_frac_flap = Float(iotype='in', desc='') # sparcap fraction of chord flap
sc_frac_edge = Float(iotype='in', desc='') # sparcap fraction of thickness edge
sc_frac_flap = Float(iotype='in', desc='spar cap fraction of chord')
sc_frac_edge = Float(iotype='in', desc='spar cap fraction of thickness')
SF_blade = Float(iotype='in', desc='') #[factor]
Slim_ext_blade = Float(iotype='in', units='MPa', desc='')
Slim_fat_blade = Float(iotype='in', units='MPa', desc='')
AddWeightFactorBlade = Float(iotype='in', desc='') # Additional weight factor for blade shell
BladeDens = Float(iotype='in', units='kg/m**3', desc='density of blades') # [kg / m^3]
BladeCostPerMass = Float(iotype='in', desc='') #[e/kg]
HubCostPerMass = Float(iotype='in', desc='') #[e/kg]
SpinnerCostPerMass = Float(iotype='in', desc='') #[e/kg]
AddWeightFactorBlade = Float(iotype='in', desc='Additional weight factor for blade shell')
BladeDens = Float(iotype='in', units='kg/m**3', desc='density of blades')
BladeWeight = Float(iotype = 'out', units = 'kg', desc = 'BladeMass' )
RootChord = Float(iotype = 'out', units = 'm', desc = 'blade root chord') # 07/09/2015 added for HubSE model
......@@ -162,16 +159,12 @@ class SEAMBladeStructure(Component):
tfinal_edge = np.maximum(text_edge, tfat_edge)
self.tfinal_flap = tfinal_flap
self.tfinal_edge = tfinal_edge
# Calculating different costs and mass
for i in range(0,int(self.Nsections)):
if i>0:
volumen = volumen+self.AddWeightFactorBlade*(r[i]-r[i-1])*((2*self.sc_frac_flap*C[i]*tfinal_flap[i]
+2*self.sc_frac_edge*thick[i]*tfinal_edge[i])+(2*self.sc_frac_flap*C[i-1]*tfinal_flap[i-1]
+2*self.sc_frac_edge*thick[i-1]*tfinal_edge[i-1]))/2.
self.BladeWeight = self.BladeDens*volumen
self.volume = np.trapz(2 * (self.sc_frac_flap * C * tfinal_flap + \
self.sc_frac_edge * thick * tfinal_edge), r)
self.volume *= self.AddWeightFactorBlade
self.BladeWeight = self.BladeDens * self.volume
print 'volumen', volumen
print 'BladeWeigth', self.BladeWeight
#
# rotor = self.BladeWeight*self.BladeCostPerMass/1e6 # Meuro
......
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