Поблочный текст программы баллистического проектирования ЛА класса Поверхность – Воздух управляемого по методу пропорциональной навигации
Московский авиационный институт
(Национальный исследовательский институт)
Аэрокосмический факультет
Кафедра: «Авиационно-ракетных комплексов»
Курсовая работа
По дисциплине: «Прикладная информатика»
На тему: «Разработка модульного ряда программ баллистического проектирования ЛА с различными способами управления»
Студент: Вавилов Михаил Андреевич
Группа: 60-205С
Преподаватель: доцент кафедры «Авиационные космические системы»: Петраш Виктор Яковлевич
Москва, 2015 год
Баллистический расчёт – это определение траекторных характеристик путём численного интегрирования уравнений движения совместно с уравнениями, приближенно описывающими систему управления ЛА.
Баллистическое проектирование – это баллистический расчёт ЛА с определением потребного запаса топлива, необходимого для реализации полученной траектории.
Укрупненная блок-схема баллистического расчета
Исходные данные программы баллистического проектирования ЛА класса Поверхность - Воздух по методу пропорциональной навигации
v0
ЗHAЧEHИE V ПPИ T = T0
50 [М/С
tt0
ЗHAЧEHИE TT ПPИ T = T0
60 [ГPAД]
vtr
УCKOPEHИE ЛA HA УЧACTKE PAЗГOHA
150 [M/C**2]
xk
KOOPДИHATA KOHЦA ПOЛETA
30000 [M]
vc
CKOPOCTЬ ЦEЛИ
500 [M/C]
tk
ПOЛHOE BPEMЯ ПOЛETA
35 [C]
hc
BЫCOTA ПOЛETA ЦEЛИ
15000 [M]
mut1
OTH.MACCA TOПЛИBA HA 1-M PEЖИME PAБOTЫ PДTT
0.15 [Б/P]
mut
OTHOCИTEЛЬHAЯ MACCA TOПЛИBA ЛA
0.4 [Б/P]
dt
ШAГ ЧИCЛEHHOГO ИHTEГPИPOBAHИЯ
0.2 [C]
ttc
УГOЛ HAKЛOHA TPAEKTOPИИ ЦEЛИ
3.14 [PAД]
k3
HABИГAЦИOHHAЯ ПOCTOЯHHAЯ
0.5 [Б/P]
taud
BPEMЯ PAБOTЫ MAPШEBOГO ДBИГATEЛЯ
20 [C]
p
УДEЛЬHAЯ HAГPУЗKA HA KPЫЛЬЯ
500 [KГ/M**2]
pud
УДEЛЬHЫЙ ИMПУЛЬC TЯГИ MAPШEBOГO ДBИГATEЛЯ
2500 [M/C]
H
BЫCOTA ПOЛETA
15000 [M]
JAER
ПPИЗHAK ПEPBOГO OБPAЩEHИЯ K ПOДПPOГPAMME AER
0 [Б/P]
JS
ПPИЗHAK AЭPOДИHAMИЧECKOЙ CXEMЫ ЛA
3 [Б/P]
M0
CTAPTOBAЯ MACCA ЛA (ИЛИ MACCA BTOPOЙ CTУПEHИ)
500 [KГ]
ROCP
CPEДHЯЯ ПЛOTHOCTЬ KOMПOHOBKИ ЛA
1200 [KГ/M**3]
LLF
УДЛИHEHИE KOPПУCA
15 [Б/P]
LLH
УДЛИHEHИE ЗAOCTPEHHOЙ HOCOBOЙ ЧACTИ KOPПУCA
3 [Б/P]
RHC
OTHOCИT.PAДИУC ЗATУПЛEHИЯ HOCOBOЙ ЧACTИ KOPПУCA
0 [Б/P]
IK
ЧИCЛO KOHCOЛEЙ KPЫЛЬEB
4 [Б/P]
CC
OTHOCИTEЛЬHAЯ TOЛЩИHA ПPOФИЛЯ KPЫЛA
0.03 [Б/P]
LLK
УДЛИHEHИE ДBУX KOHCOЛEЙ KPЫЛЬEB
1.2 [Б/P]
ETAK
CУЖEHИE KOHCOЛEЙ KPЫЛЬEB
4 [Б/P]
KPЫЛA
TX1 TAHГEHC УГЛA CTPEЛOBИДH.ПO ЗAДHEЙ KPOMKE
0 [Б/P]
PC
TЯГOBOOPУЖEHHOCTЬ ЛA
10 [Б/P]
V
CKOPOCTЬ ЛA
450 [M/C]
AA
УГOЛ ATAKИ
0 [ГPAД]
Поблочный текст программы баллистического проектирования ЛА класса Поверхность – Воздух управляемого по методу пропорциональной навигации
C+ v0,TT,VC,vtr,xk,tk,hc,mut1,ttc,TTO
TT0=TT0/57.3
C- t,v,tt,h,x,ff,r,mu
702 JAER=0
C/nachpv/Начальные условия для ЛА поверхность-воздух
c начальные условия
v=v0
t=v/vtr
h=0.5*v*t*sin(tt)
x=h/tan(tt)
ff=atan((hc-h)/(xk-x-vc*(tk-t)*cos(ttc)))
r=(hc-h)/sin(ff)
mu=1-mut1
g=9.81
print 10
write (17,10)
10 format (2x,' t ',2x,' v ',2x,' teta ',2x,' h ',2x,
*' x ',2x,' mu ',2x,' ny ',2x,'aa')
C+ K3
C- NY
C/UPRAVPV/управляющая перегрузка
c управляющая перегрузка
15 rt=vc*cos(ff-ttc)-v*cos(ff-tt)
fft=1/r*(v*sin(ff-tt)-v*sin(ff-ttc))
ny=k3*abs(rt)*fft
C+ TAUD,PUD,MUT
C-PC,MUC
C/UPRAVDV/УПРАВЛЕНИЕ ДВИГАТЕЛЕМ
c двигатель
if (t.lt.taud) then
muc=mut/taud
pc=muc*pud/g
else
muc=0
pc=0
end if
C- AA,VT,TTT,MUTT,NT,XT
C/DVIGEN/УРАВНЕНИЯ ДВИЖЕНИЯ
c уравнения движения
aa=g*mu*ny/((muc*pud/57.3)+(cya*roh*v*v/(2*p)))
vt=1/mu*(muc*pud-roh*v*v*cx/(2*p))-9.81*sin(tt)
ttt=g/v*(ny-cos(tt))
mutt=-muc
ht=v*sin(tt)
xt=v*cos(tt)
C+ DT
C- V,TT,MU,H,X,R,FF
C/INTEGRPV/ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ
c интегрирование
v=v+vt*dt
tt=tt+ttt*dt
mu=mu+mutt*dt
h=h+ht*dt
x=x+xt*dt
r=r+rt*dt
ff=ff+fft*dt
c вывод стоки результатов
print 3,t,v,tt,h,x,mu,ny,aa
write(17,3) t,v,tt,h,x,mu,ny,aa
3 format (2x,f5.1,2x,f9.2,2x,f9.5,2x,f9.2,2x,f9.2,2x,f9.5,2x,f9.4,