Моделирование движения самолёта класса тяжелых маневренных истребителей палубного базирования типа Су-33 в продольном и боковом канале, страница 4

Задачей одной из работ было определение статических характеристик проектируемого самолёта, в частности: определение угла aбал, а также определение балансировочных положений органов управления. Для выполнения поставленной задачи была разработана следующая программа

Текст программы, написанной в среде САПР «Динамика» для моделирования продольного движения проектируемого самолёта.

<                 M o d e l

<                    o f

<               A I R C R A F T

<                 T-10KUB

<             LONGITUDINAL CHANNEL

<

 EN::V:VZAD:

 EN::H:HZAD:

 EN::A:ALFA:

 EN::WZ:WZZA:

<

 EN::RDR:RDR0:

 EN::RDL:RDL0:

 EN::FIST:STAB:

 EN::VPK:VPPK:

 EN::FLV:15.:

 EN::ELV:15.:

 EN::CMAN:CMA0:

<

 EN::DDXT:DXT0:

 EN::ONEI:-1.:  

 EN::NUL:0.:

 EN::FIVE:5.:

<                    A t m o s p h e r e

 ATHV:H,VMS:AM,Q,R,AZV,PH,TH::

 AVIN:Q:VIN::

 AQSG:AM,PH:QSG::

 WVPR:QSG:VPR::

 APTM:PH:PST::

 ZMZV:QSG,PST:AME::

 EN::CVME:^CVME:

 KK:QSG:QSDU:CQSD:

<          M A I N  C H A N N E L

<  L o n g i t u d i n a l  c h a n n e l

<            F l y - b y - W i r e

<                S D U - 10 M K - K U

 KK:V:VMS:0.27777778:

 KK:WZ:WZR:0.017452:

 KK:A:AR:0.017452:

<  E l e v a t o r  a c t u a t o r s

 DN1:FIST:FI:-20.,-20.,15.,15.:

 KK:FI:F:1.:

<<<   C a n a r d  c h a n n e l (PGO)

 DN1:FIST:FEN:-40.,-40.,30.,30.:

 KS:FEN,ONEI:DC1:2:

 DI:DC1:DC4:11.:

 KK:DC4:DC5:-1.5:

 DN2:A:DCAL:-6.5,5.,20.,-15.,35.,-50.:

 DN1:DCAL:DCA1:-25.,-25.,0.,0.:

 KP1:DCAL,DCA1,CMAN:DCA:0.5:

 KN1:DCA,AME:DC7:0.4,1.,1.,0.:

 KS:DC7,DC5:DC8:2:

 DN1:DC8:DCEN:-50.,-50.,0.,0.:

 KK:DCEN:DC:1.:

<<<      S i n c h r o  f l a p e r o n s'  c h a n n e l  

 DN3:A:FFLA:2.,0.,6.,5.,15.,5.,35.,-15.:

 KP1:FFLA,FLV,VPK:SFLA:0.5:

<  F l a p e r o n  a c t u a t o r s

 KK:SFLA:FE:1.:

<<<      S i n c h r o  e l e r o n' s  c h a n n e l 

 KP1:NUL,ELV,VPK:ELPR:0.5:

 KK:ELPR:EL:1.:

<<<   L e a d i n g  e a g e  c h a n n e l

 DN1:A:DNVP:8.,4.,15.,30.:

 DN1:A:DNPK:3.,-4.,10.,-30.:

 KP1:DNPK,DNVP,VPK:DNO:0.5:

<  N o z z e l s  a c t u a t o r s

 EN::FPL:0.:

 EN::FPR:0.:

<         E n g i n s  m o d e l ( v e r.2 )

 DN3:RDL:RSGL:0.,0.15,68.,0.65,75.,0.75,106.,1.:

 DN2:RSGL:PL:0.15,^PMG,0.65,^PM,1.,^PF:

 DN3:RDR:RSGR:0.,0.15,68.,0.65,75.,0.75,106.,1.:

 DN2:RSGR:PR:0.15,^PMG,0.65,^PM,1.,^PF:

<        R i g h t  p a r t s  o f  e q u a t i o n s

<           ( F o r s e s  a n d  M o m e n t s )

<

 ZZXT:VPK:ZZI,XTT,G:GT,ZAGR,dIz:

 KS:XTT,DDXT:XT:2:

 Z10L:A,F,DC,AM,Q,WZR,VMS,PL,PR,DNO,FE,EL,FPL,FPR,XT,VPK,G:

ZM,ZMAT,XN,YN,ZMS,PX,PY,CXB,CYB:62.04,4.64,0.,0.,1.,1.:

<       E X I T S  F O R  F L I G H T  S I M U L A T I O N

 EXIT:EXI:A,AR,WZ,WZR,ZMS,ZM,ZMAT,XN,YN,PX,PY,CXB,CYB,

H,F,FE,EL,DC,DNO,PL,PR,

AME,AM,VPR,VIN,V,QSG,QSDU,

XT,G,FIST,VPK:

К данному файлу в системе САПР «Динамика» подключался файл вычисления правых частей уравнений движения.

Текст файла с вычислением правых частей уравнений движения (сил и моментов)

Z10L

 Z10L:A,F,DC,AM,Q,WZR,V,PDL,PDR,DNO,FE,EL,FSLV,FSPR,XT,VPK,G:

ZM,ZMAT,XN,YN,ZMS,PX,PY,CXB,CYB:S,BA,ASOP,AIRB,FHS1,FHS2:

      IF(_V.LE.0.001)_V=1.

      _BAV=_BA/_V

      _QSG=_Q*_S/_G

      _AR=_A/57.3

      _QSB=_Q*_S*_BA

      _XP=6.2-(_XT-0.29)*_BA

      _PDV=_PDL+_PDR

      _FRL=_FSLV/57.3

      _FRR=_FSPR/57.3

      _CALC=COS(_ASOP/57.3)

C

      _CCYF=1.

      _CMZF=1.

      IF(_VPK.LT.0.5) GOTO _20

      _FF=_F

      IF(_F.LT.5.AND._F.GE.0.) _F=5.

      IF(_F.LT.0.AND._F.GT.-5.) _F=-5.

      _CCYF=^CCYF

      _CMZF=^CMZF

      _F=_FF