Mathesis Wiki

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)