=====Klasse Welt===== 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)