Образовательный проект «SnakeProject» Михаила Козлова

⇒ Python ⇐

CISCO

Voice(Asterisk\Cisco)

Microsoft

Powershell

SQL\T-SQL

FreeBSD and Nix

1С

Общая

WEB Разработка

ORACLE SQL \ JAVA

Мото

Стрельба, пневматика, оружие

Саморазвитие и психология

Баллистический калькулятор с учетом ветра и притяжения на python

Баллистический калькулятор с учетом ветра и притяжения на python

Сделан расчет влияния сопротивления воздуха на движение снаряда

В коде имеются 2 разных примера на сверхзвуке и дозвуке

Сверхзвук рассчитан проверен до меня умными людьми, дозвук проверял сам

Выводятся 2 графика, с учетом и без учета сопротивления и притяжения

from pylab import * import math #Пример, гаубица Д-20 #152мм снаряд, начальная скорость 655м/с, масса 43.5кг #Коэффициент сопротивления 0.219 найден в интернете #Калибр мм (диаметр пули или снаряда) v1_cal = 152 #Масса v1_m = 43.5 #Начальная скорость м/с v1_v = 655 #Угол запуска в градусах v1_a_gr = 40 #Коэффициент сопротивления v1_Cd = 0.219 #Пример, пневматика 4.5мм, 3 джоуля, 100м/с, пуля 0.5гр v2_cal = 4.5 v2_m = 0.0005 v2_v = 100 v2_a_gr = 0.41 #Реальный коэффициент на дозвуке домножения ~33 v2_Cd = 0.023 * 33 #Плотность воздуха v_rho = 1.225 def rasch(cal, m, v, a_gr, Cd, rho, weapon): #Притяжение м/с g = 9.8 #Площадь поперечного сечения снаряда A = math.pi * pow(cal / 1000, 2.0) / 4 #Сила сопротивления alpha = rho * Cd * A / 2.0 #Гравитация beta = alpha / m #Заданные значения точки вылета на графике X0 = 0.0 Y0 = 0.0 #Формула горизонтальной и вертикальной скорости #Vx=V*cos(a), Vy=V*sin(a) #Где Vx - горизонтальная скорость (м/с) #Vy - вертикальная скорость (м/с) #V - величина скорости (м/с) #a - угол изменения скорости (градусы, python принимает в радианах) a_rad = math.pi * a_gr / 180 Vx0 = v * math.cos(a_rad) Vy0 = v * math.sin(a_rad) #Шаг времени steps = 1000 t_HIT = 2.0*Vy0/g dt = t_HIT / steps #График без drag X_ND = list() Y_ND = list() for i in range(steps+1): X_ND.append(X0 + Vx0 * dt * i) Y_ND.append(Y0 + Vy0 * dt * i - 0.5 * g * pow(dt * i,2.0)) #График с drag X_WD = list() Y_WD = list() Vx_WD = list() Vy_WD = list() #Добавили начальные значения на график X_WD.append(X0) Y_WD.append(Y0) Vx_WD.append(Vx0) Vy_WD.append(Vy0) stop = 0 print(weapon) print("Начальная скорость по x: " + str(Vx0) + "Начальная скорость по y: " + str(Vy0)) for i in range(1,steps+1): if stop != 1: #Скорости в динамике speed = pow(pow(Vx_WD[i-1],2.0)+pow(Vy_WD[i-1],2.0),0.5) print( "Время: " + str(dt * i) +"Скорость: " + str(speed) ) #С сопротивлением и притяжением Vx_WD.append(Vx_WD[i-1] * (1.0 - beta * speed * dt)) Vy_WD.append(Vy_WD[i-1] + ( - g - beta * Vy_WD[i-1] * speed) * dt) #Без сопротивления X_WD.append(X_WD[i-1] + Vx_WD[i-1] * dt) Y_WD.append(Y_WD[i-1] + Vy_WD[i-1] * dt) # Прекратить расчеты по достижении 0 на графике if Y_WD[i] <= 0.0: #print( "Нс: " + str(pow(Vx_WD[i-1],2.0)) + "Пс: " + str(pow(Vy_WD[i-1],2.0)) ) stop = 1 #График без drag plot(X_ND, Y_ND) #График с drag plot(X_WD, Y_WD) show() if __name__ == '__main__': rasch(v1_cal, v1_m, v1_v, v1_a_gr, v1_Cd, v_rho, "Д-20") rasch(v2_cal, v2_m, v2_v, v2_a_gr, v2_Cd, v_rho, "3-Дж")

Комментарии пользователей