atmosphere.py

Функции расчёта параметров атмосферы
mail@pastecode.io avatar
unknown
plain_text
a year ago
3.3 kB
1
Indexable
Never
#Импорт данных из data.py
from mainProgect import data
import math

r = data.r
g0 = data.g0
R = data.R
i = data.i
h = data.h

# Функции атмосферы
def geopotencial_v():
    H = (r * h) / (r + h)  # Геопотенциальная высота
    return H
def yskorenie_s_p():
    g = g0 * (r / (r + h)) ** 2  # Ускорение свободного падения
    return g
def parametry_t():
    H = geopotencial_v()
    # Параметры, для определения температуры и давления
    if (H >= -2001) and (H < 0):     # 3-4 стр. ГОСТа
        Hx = -2000
        batm = -0.0065
        Tx = 301.15
        px = 127774
    elif (H >= 0) and (H <= 11000):  # 5-15 стр. ГОСТа
        Hx = 0
        batm = -0.0065
        Tx = 288.15
        px = 101325
    elif (H > 11000) and (H < 20000):  # 15-24 стр. ГОСТа
        Hx = 11000
        Tx = 216.65
        batm = 0
        px = 22632
    elif (H >= 20000) and (H < 32000):  # 25-36 стр. ГОСТа
        Hx = 20000
        batm = 0.001
        Tx = 216.65
        px = 5474.87
    elif (H >= 32000) and (H < 47000):  # 37-43 стр. ГОСТа
        Hx = 32000
        batm = 0.0028
        Tx = 228.65
        px = 868.014
    elif (H >= 47000) and (H < 51000):  # 43-46 стр. ГОСТа
        Hx = 47000
        Tx = 270.65
        batm = 0
        px = 110.906
    elif (H >= 51000) and (H < 71000):  # 45-52 стр. ГОСТа
        Hx = 51000
        batm = -0.0028
        Tx = 270.650
        px = 66.9384
    elif (H >= 71000) and (H < 85000):  # 51-53 и 159 стр. ГОСТа
        Hx = 71000
        batm = -0.002
        Tx = 214.65
        px = 3.95639
    return dict(Hx = Hx, batm = batm, Tx = Tx, px = px)
def temperatura():
    p = parametry_t()
    H = geopotencial_v()
    T = p["Tx"] + p["batm"] * (H - p["Hx"])       # Температура
    return T
def davlenie():
    p = parametry_t()
    H = geopotencial_v()
    T = temperatura()
    if p["batm"] != 0:  # Давление
        p = 10 ** (math.log10(p["px"]) - (g0 / (p["batm"] * R)) * math.log10(T / p["Tx"]))
    else:
        p = 10 ** (math.log10(p["px"]) - ((0.434294 * g0 * (H - p["Hx"])) / (T * R)))
    return p
def plotnost():
    p = davlenie()
    T = temperatura()
    ro = p / (R * T)  # Плотность
    return ro
def skorost_z():
    T = temperatura()
    a = 20.046796 * math.sqrt(T)  # Скорость звука
    return a
def vivod():
    H = geopotencial_v()
    g = yskorenie_s_p()
    T = temperatura()
    p = davlenie()
    ro = plotnost()
    a = skorost_z()
    print("Геометрическая высота: " + "{:.0f}".format(h) + ' м')
    print("Геопотенциальная высота: " + "{:.0f}".format(H) + ' м')
    print("Ускорение свободного падения: "+str(round(g,3))+' м/с^2')
    print("Температура: "+str(round(T/100,3))+' * 10^2 K')
    print("Давление: "+str(round(p/10000,3))+' * 10^4'+' Па')
    print('Плотность: '+str(round(ro*10,3))+' * 10^-1'+' кг/м^3')
    print("Скорость звука: "+str(round(a/100,3))+' * 10^2 м/с')
    print(' ')