Hier werden die Unterschiede zwischen zwei Versionen gezeigt.
Beide Seiten der vorigen Revision Vorhergehende Überarbeitung Nächste Überarbeitung | Vorhergehende Überarbeitung | ||
ss2021:project5:flugzeugsimulation:protokoll [2021/07/08 16:24] emilyeberhardt |
ss2021:project5:flugzeugsimulation:protokoll [2021/07/08 16:29] (aktuell) p_hilsberg |
||
---|---|---|---|
Zeile 31: | Zeile 31: | ||
* Wie spielen die Kräfte und Anstellwinkel mit dem Auftrieb zusammen? | * Wie spielen die Kräfte und Anstellwinkel mit dem Auftrieb zusammen? | ||
* Besprechung der Darstellung | * Besprechung der Darstellung | ||
+ | * [[ss2021:project5:Flugzeugsimulation:Klasse Welt]] | ||
{{:ss2021:project5:flugzeugsimulation:flugphasen_01.svg.png?200|}} | {{:ss2021:project5:flugzeugsimulation:flugphasen_01.svg.png?200|}} | ||
- | Quelle: https://de.wikipedia.org/wiki/Flugphase | + | Quelle: https://de.wikipedia.org/wiki/Flugphase |
- | + | ||
- | **Klasse Welt** | + | |
- | + | ||
- | <code python> | + | |
- | import math as m | + | |
- | + | ||
- | class Welt(): | + | |
- | + | ||
- | def __init__(self): | + | |
- | + | ||
- | self.position = 90 # in m über NN | + | |
- | self.temperatur_nn = 273,15 + 15 # in K (15 Grad) auf NN | + | |
- | + | ||
- | self.welt_masse = 5.97223*10**24 # Erdmasse in kg | + | |
- | self.welt_radius = 6378137 # Erdradius in m | + | |
- | + | ||
- | self.welt_Gkonst = 6.6739*10**(-11) # Gravitationskonstante in m^3/(kg*s^2) | + | |
- | self.welt_Rkonst = 287.058 # in J/kg*K (spez. Gaskonstante trockene Luft) | + | |
- | self.welt_mol = 0.02896 # in kg/mol, mittlere molare Masse | + | |
- | + | ||
- | def temperatur(self, hoehe=12500): | + | |
- | if hoehe == 0: | + | |
- | return 273.15+15 # 15 Grad Celsius | + | |
- | elif hoehe < 10001: | + | |
- | return self.temperatur(hoehe = 0)-(hoehe/140) | + | |
- | else: # Temperatur ab 10.000m sinkt nicht mehr | + | |
- | return self.temperatur(hoehe = 10000) | + | |
- | + | ||
- | def ortsfaktor(self, hoehe=12500): | + | |
- | return self.welt_Gkonst*self.welt_masse/(self.welt_radius+hoehe)**2 | + | |
- | + | ||
- | def schwerkraft(self, hoehe=12500, masse=200000): | + | |
- | return masse*self.ortsfaktor(hoehe) | + | |
- | # eigentlich von flugzeugmasse abhängig (vernachlässigbar (?)) | + | |
- | + | ||
- | def luftdruck(self, hoehe=12500): | + | |
- | if hoehe == 0: | + | |
- | return 101325 # in Pa bei 15 Grad auf NN | + | |
- | elif hoehe < 10001: | + | |
- | # return self.luftdruck(0)*m.exp(-hoehe*(self.welt_mol*self.ortsfaktor(hoehe)/(self.welt_Rkonst*self.temperatur(hoehe)))) | + | |
- | return self.luftdruck(0)*m.exp(hoehe/-8000) | + | |
- | else: | + | |
- | return self.luftdruck(10000) | + | |
- | + | ||
- | def luftdichte(self, hoehe=12500): | + | |
- | if hoehe < 10001: | + | |
- | return self.luftdruck(hoehe)*self.welt_mol/(self.ortsfaktor(hoehe)*self.temperatur(hoehe)) | + | |
- | else: | + | |
- | return self.luftdichte(10000) | + | |
- | + | ||
- | def widerstand(self, hoehe=12500, flaeche=361.6, v=250, widerstandsbeiwert = 0.08): | + | |
- | return (1/2)*self.luftdichte(hoehe)*widerstandsbeiwert*flaeche*v**2 | + | |
- | # Luftwiderstandsbeiwert stromlinienförmiger Objekte (z.b. Flugzeug) | + | |
- | # am besten die einzelnen parameter des flg ersetzen durch def widerstand(self, objekt=A330) | + | |
- | + | ||
- | </code> | + | |
Zeile 107: | Zeile 51: | ||
* Überarbeitung, Anpassung der Klassen | * Überarbeitung, Anpassung der Klassen | ||
* Code vorläufig sortiert | * Code vorläufig sortiert | ||
- | * Klasse Flugzeug | + | * [[ss2021:project5:Flugzeugsimulation:Klasse Flugzeug]] |