Поблочный текст программы баллистического проектирования ЛА класса Поверхность – Воздух управляемого по методу пропорциональной навигации

Московский авиационный институт

(Национальный исследовательский институт)

Аэрокосмический факультет

Кафедра: «Авиационно-ракетных комплексов»

Курсовая работа

По дисциплине: «Прикладная информатика»

На тему: «Разработка модульного ряда программ баллистического проектирования ЛА с различными способами управления»

Студент: Вавилов Михаил Андреевич

Группа: 60-205С

Преподаватель: доцент кафедры «Авиационные космические системы»: Петраш Виктор Яковлевич

Москва, 2015 год

Баллистический расчёт – это определение траекторных характеристик путём численного интегрирования уравнений движения совместно с уравнениями, приближенно описывающими систему управления ЛА.

Баллистическое проектирование – это баллистический расчёт ЛА с определением потребного запаса топлива, необходимого для реализации полученной траектории.

Укрупненная блок-схема баллистического расчета

Исходные данные программы баллистического проектирования ЛА класса Поверхность - Воздух по методу пропорциональной навигации

Поблочный текст программы баллистического проектирования ЛА класса Поверхность – Воздух управляемого по методу пропорциональной навигации

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,

*2x,f9.1)

if (r.gt.0) then

t=t+dt

goto 15

end if

C+ Q

C- MUT

C/KORMUTPV/КОРРЕКЦИЯ СЕКУНДНОГО РАСХОДА

c коррекция от масс топлива

if(abs(t-tk).gt.DT) then

MUT=MUT*(T+Q)/(tk+Q)

print 34,MUT

write (17,34)MUT

34 FORMAT(30x,'коррекция MUT',3x,'MUT=',f7.4)

goto 702

endif

print 35,MUT

WRITE (17,35) MUT

35 FORMAT (30x,'MUT=',f7.4)

Поиск по сайту: