Determinação das características balístico-temporais do movimento do centro de massa de um paraquedista que aterra de um avião

1. Introdução


Para determinar as características balístico-temporais do movimento do centro de massa do paraquedista, é preciso escolher um modelo matemático simplificado, bastante acessível para a pesquisa analítica e, ao mesmo tempo, preservar as características mais características do objeto original.

Para construir modelos matemáticos simplificados do movimento do paraquedista, é realizada uma análise, determinação, sistematização de parâmetros constantes e de tempo.

Atualmente, não existem métodos regulares e suficientemente justificados para a construção de modelos matemáticos não lineares, para resolver problemas particulares, com a preparação correta dos sistemas iniciais de equações diferenciais não lineares, métodos numéricos para resolvê-los podem fornecer resultados bastante adequados.

O objetivo desta publicação é compilar e resolver numericamente sistemas de equações diferenciais que descrevem todos os estágios do movimento de um paraquedista pousando de um avião, levando em consideração o efeito das mudanças de altitude e temperatura da densidade de massa do ar.

Características balístico-temporais do movimento de um paraquedista


Parâmetros variáveis ​​constantes e limitados incluem:

N é a altura do paraquedista;
V0 - velocidade da aeronave;
k - peso, altura do paraquedista;
g é a aceleração da gravidade;
ρ é a densidade do ar;
T é a temperatura do ar.

Parâmetros temporários (variáveis) incluem:

tn - hora de pouso,
w é a velocidade do vento;
V é a velocidade do paraquedista;
u é a velocidade dos fluxos ascendentes (descendentes);
d - desvio (distância da projeção ao solo do ponto de ejeção até o ponto de aterrissagem);
C é o coeficiente de arrasto do objeto de pouso;
F - seção média do objeto de pouso.

Estágios de salto


O primeiro estágio é uma queda livre após a separação da aeronave:



O segundo estágio é uma diminuição de um pára-quedas estabilizador:



A principal propriedade do paraquedas estabilizador é a estabilização do paraquedista na posição mais conveniente para o paraquedas principal.

A terceira etapa - encher a cúpula do paraquedas principal:



O quarto estágio é uma diminuição do paraquedas aberto:



Compilação de um sistema de equações diferenciais para todas as etapas do salto de paraquedas


Escolhemos um sistema de coordenadas fixo OXY centrado no ponto de ejeção O. O eixo OX coincide com a direção do componente horizontal da velocidade da aeronave. O eixo OY é direcionado verticalmente para cima na direção oposta à velocidade vertical do paraquedista.

Vamos assumir que o movimento do pára-quedista é plano e ocorre no plano OXY. Esse modelo de salto pode ser considerado um modelo de salto em clima calmo, sem levar em consideração a influência do vento.

Acreditamos que, além do peso, o paraquedista é afetado pela força da resistência do ar proporcional ao quadrado da velocidade do paraquedista:

,
onde: , - densidade do ar, coeficiente de arrasto C, seção F do corpo.

Com o aumento da altura, a temperatura do ar muda:



A temperatura mínima já é atingida a uma altitude de 10 km. e é de -55 ° C. A densidade do ar também depende da pressão. Portanto, ao calcular a balística de um salto de para-quedas, é conveniente usar a seguinte fórmula para determinar a densidade do ar [1]:

,
onde K / m; - temperatura ao nível do mar; y é a altura em m; - densidade do ar em y = 0; .

Listagem
# -*- coding: utf8 -*- from numpy import* import matplotlib.pyplot as plt beta=1.25*10**-4 alfa=6.5*10**-3 T0=300 y=arange(0,10000,1) x=[exp(-beta*w*T0/(T0-alfa*w)) for w in y] plt.title("  \n    ") plt.plot(x,y) plt.ylabel('  .') plt.xlabel(' ') plt.grid(True) plt.show() 




Na prática dos cálculos, o quadrado do crescimento é tomado como a seção mediana; o valor de C é encontrado na tabela [2]:



Θ indica o ângulo de inclinação da trajetória. Sob as suposições feitas para os componentes , vetores de velocidade V, temos:



Dividindo por m os lados esquerdo e direito das equações do sistema resultante e denotando através de r, obtemos:

(1)

Escrevemos as equações de movimento de um pára-quedista na forma de um sistema de equações diferenciais para as funções V, θ, y (t), x (t).

Dado que:

,
,
e tempo de diferenciação: , levando em consideração o sistema de equações (1), obtemos:

,
.

Assim, nas condições iniciais:

nós temos o seguinte sistema de equações diferenciais:



Solução numérica de um sistema de equações diferenciais (2) usando Python


Para resolver (2), reescrevemos da seguinte forma, introduzindo forças de arrasto estabilizadoras controladas pelo tempo e pela densidade do ar



e principal



pára-quedas, respectivamente, multiplicados pelas funções de controle do tempo ks (t) e ko (t):



,

onde: –– tempo de queda livre do paraquedista; - o tempo de operação do pára-quedas de estabilização até a abertura do principal.

(3)

Lista completa do programa ajustado para alterações na densidade do ar
 # -*- coding: utf8 -*- from numpy import* from scipy.integrate import odeint import matplotlib.pyplot as plt g=9.81#   /2 rou0=1.29#  /3 Cp=0.3#    Cpar=0.75#     mp=120#    100+20=120 . s=0.8#   2 ss=1.5#   2 Cpars=0.35#     S=70#   2 tsp=40#    c. tssp=3 #      c. tpr=210#     . h=4000.0#    . V0=100.0#    / beta=1.25*10**-4 alfa=6.5*10**-3 T0=300 def ks(t): if t<= tsp: z=0 else: z=1 return z def ko(t): if t<= tsp+tssp: z=0 else: z=1 return z def f(y,t): y1,y2,y3,y4=y rou=rou0*exp(-beta*y1*T0/(T0-alfa*y1)) rp=0.5*Cp*rou*s/mp rpars=0.5*Cpars*rou*ss/mp rpar=0.5*Cpar*rou*S/mp return [-y2*sin(y3),g*sin(y3)-rp*(y2**2)-ko(t)*rpar*(y2**2)-ks(t)*rpars*(y2**2),g*cos(y3)/y2,y2*cos(y3)] t=arange(0,tpr,1) y0=[h,100,0,0] [y1,y2,y3,y4]=odeint(f, y0,t, full_output=False).T plt.figure() plt.title('     : %s .. \n   : %s /. \n   : %s /'%(h,V0,round(y2[tpr-tsp-tssp],0))) plt.plot(y1,y2*sin(y3),label='  .: %s c. \n     : %s . \n     : %s . '%(tsp,tssp,tpr-tsp-tssp)) plt.xlabel('  .') plt.ylabel('   /') plt.legend(loc='best') plt.grid(True) plt.figure() plt.title("-         \n   ") plt.plot(y1,y4) plt.xlabel('  .') plt.ylabel('     .') plt.grid(True) plt.figure() plt.title("     ") plt.plot(t,y2*sin(y3)) plt.xlabel('  .') plt.ylabel('   /.') plt.grid(True) plt.show() 


Temos:







Levar em consideração o ar rarefeito levou a um aumento na taxa de queda livre e mudou a natureza da trajetória nesta seção.

Este problema pode ser resolvido com a ajuda de um sistema de duas equações diferenciais, que são dadas abaixo (excluindo pára-quedas e alterações na densidade do ar):



A mudança nas forças de resistência e na densidade do ar é mostrada na lista abaixo do spoiler, levando em consideração o exposto acima e sem explicações
# - * - codificação: utf8 - * -
de importação numpy *
de scipy.integrate import odeint
importar matplotlib.pyplot como plt
m = 100
r0 = 1,3
c1 = 0,3
c2 = 0,6
c3 = 0,5
c4 = 0,75
S = 70
s = 0,8
ss = 1,5
g = 9,8
colher de chá = 6
tsbp = 10
tp = 90,0
h = 1000.0
beta = 1,25 * 10 ** - 4
alfa = 6,5 * 10 ** - 3
T0 = ​​300
def ks (t):
se t <tsp:
z = 0
mais:
z = 1
retornar z
def ko (t):
se t <tsp + tsbp:
z = 0
mais:
z = 1
retornar z
# dy1 / dt = y2
# dy2 / dt = g- (k1 * y2 ** 2) / m- (k2 * y2) / m- (ks (t) * k3 * y2 ** 2) / m- (ko (t) * k4 * y2 ** 2) / m
def f (y, t):
y1, y2 = y
r = r0 * exp (-beta * y1 * T0 / (T0-alfa * y1))
k1 = 0,5 * r * c1 * s
k2 = 0,5 * r * c2 * s
k3 = 0,5 * r * c3 * ss
k4 = 0,5 * r * c4 * S
return [-y2, g- (k1 * y2 ** 2) / m- (k2 * y2) / m- (ks (t) * k3 * y2 ** 2) / m- (ko (t) * k4 * y2 ** 2) / m]
t = arange (0,0, tp)
y0 = [h, 0,0]
[y1, y2] = nó (f, y0, t, saída_computador = Falso) .T
plt.title ("Paraquedismo de 1000 e 800 metros")
plt.plot (y1, y2, label = 'Altitude 1000 m')
h = 800,0
colher de chá = 6
tsbp = 2
tp = 80,0
def ks (t):
se t <tsp:
z = 0
mais:
z = 1
retornar z
def ko (t):
se t <tsp + tsbp:
z = 0
mais:
z = 1
retornar z
def f (y, t):
y1, y2 = y
r = r0 * exp (-beta * y1 * T0 / (T0-alfa * y1))
k1 = 0,5 * r * c1 * s
k2 = 0,5 * r * c2 * s
k3 = 0,5 * r * c3 * ss
k4 = 0,5 * r * c4 * S
return [-y2, g- (k1 * y2 ** 2) / m- (k2 * y2) / m- (ks (t) * k3 * y2 ** 2) / m- (ko (t) * k4 * y2 ** 2) / m]
t = arange (0,0, tp)
y0 = [h, 0,0]
[y1, y2] = nó (f, y0, t, saída_computador = Falso) .T
plt.plot (y1, y2, label = 'Altura 800 m')
plt.xlabel ('Altura em metros')
plt.ylabel ('Velocidade de gravação em m / s.')
plt.legend (loc = 'melhor')
plt.grid (True)
plt.show ()

Temos:



Conclusão


São determinadas as características balístico-temporais do movimento do centro de massa de um paraquedista que aterra de um avião.

Referências


  1. Pressão atmosférica.
  2. Gerasimenko I.A. Treinamento aéreo: livro didático. M: Military Publishing, 1986. Parte 1, p. 32.

Source: https://habr.com/ru/post/pt413807/


All Articles