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

< KS:FL11,FLFF:FLEE:2:

<<

< KK:ELEE:ELEN:1.:     

< KK:FLEE:FLEN:1.:

< KK:DFEE:DFEN:1.:

< KK:RNEE:RNEN:1.:

<< 

< 

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

< 

< MP1:SLEV:FPL:0.,VFC,DFC,-15.,15.:

< MP1:SPR:FPR:0.,VFC,DFC,-15.,15.:

 EN::FPL:FPL0:

 EN::FPR:FPR0:

< 

<        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 )

< 

 ZIXT:VPK:XXI,YYI,ZZI,XYI,XTT,G:GT,ZAGR,dIx,dIy,dIz:

 KS:XTT,DDXT:XT:2:

< 

 ZKUB:A,F,FE,EL,FD,DE,ELR,DN,DNO,DC,FPR,FPL,B,WXR,WYR,GAM,VMS,Q,H,

AM,PL,PR,G,XT:ZN,XM,YM,CZ,DP,XMBT,YMBT,CZB,XMB,YMB:

62.04,14.7,32.,1.2,1.,0.,0.:

< 

<     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:EXT:A,AR,B,BR,GAM,GAMR,WX,WXR,WY,WYR,

V,VIN,VPR,AM,AME,Q,QSG,YMAO,YMIO,VPK,H,

PL,PR,DC,DNO,DN,EL,ELR,FE,DE,FI,FD,

XXI,YYI,ZZI,XYI,XT,G,ZN,XM,YM,CZ,DP,CZB,XMB,YMB:

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

 ZKUB:A,F,FE,EL,FD,DE,ELR,DN,DNO,DC,FPR,FPL,B,WXR,WYR,GAM,V,Q,H,

AM,PL,PR,G,XT:ZN,XM,YM,CZ,DP,XMBT,YMBT,CZB,XMB,YMB:

S,BL,ALC,ZDV,CMXW,DMXI,DMYI:

      _AR=_A/57.3

      _BL2V=_BL/(2.*_V)

      _QSG=_Q*_S/_G

      _QSL=_Q*_S*_BL

      _RPH=1.0-6.E-5*_Q

      _RQST=1.0-2.E-5*_Q

      _XP=6.2-(_XT-0.28)*_BL

      _CALC=COS(_ALC/57.3)

      _SALC=SIN(_ALC/57.3)

      _FRR=_FPR/57.3

      _FRL=_FPL/57.3

      _SFRR=SIN(_FRR)

      _CFRR=COS(_FRR)

      _SFRL=SIN(_FRL)

      _CFRL=COS(_FRL)

      _XMBT=0.

      _YMBT=0.

C******************** T-10MK **************************

C      Интерполирование                                                      *

C******************************************************

      _CZB=^CZBM*^CZBA

      _XMB=^XMB*^CMXB+_DMXI

      _YMB=^YMB*^CMYB+_DMYI

C******************************************************

C      Вычисление M_X и M_Y от элеронов                                           *

C******************************************************

      _CMEL=^CMEL

      _XMEL=(^XMEL/(-10.))*_CMEL*_ELR

      _YMEL=(^YMEL/(-10.))*_CMEL*_ELR

C******************** T-10MK **************************

C       Расчёт M_X и M_Y от дифференциального стабилизатора,       *

C                  флаперонов и руля направления                                  *

C******************************************************

_XMFD=((^XMFD/(-5.))*^CMXD*_FD*_RQST)+(^XMFL/(-10.))*_DE*^CMXF

      _YMFD=((^YMFD/5.)*^CMYD*_FD*_RQST)+(^YMFL/(-10.))*_DE*^CMYF

      _CRN=1.

      IF(_DN.GE.0.)_CRN=-1.

      _DN=_DN*_CRN

      _DMYN=1.1*^DMYN*_CRN*^CMYR

      _DMXN=1.1*^DMXN*_CRN*^CMXR

      _DCZN=1.1*^CZDN*_CRN*^CCZR*_DN

      _DN=_DN*_CRN

C******************************************************

C        Расчёт суммарных управляющих моментов                        *

C******************************************************

 _XM07=(_XMB*_B)+_XMFD+_DMXN*_RPH+_XMEL

      _YM07=(_YMB*_B)+_YMFD+_YMEL*_RPH+_DMYN

      _CZ=((_CZB)*_B)+_DCZN

C******************************************************

C        Расчёт моментов от вектора тяги                                          *

C******************************************************

  _DP=(_PL*_SFRL-_PR*_SFRR)

      _DPZ=_DP*_SALC

      _DPY=_DP*_CALC

      _YMSP=(_PR*_CFRR-_PL*_CFRL+_DPZ*_XP*1.1)/_QSL

      _XMSP=_DPY*1.3*_CALC/_QSL

      _ZNSP=_DPZ/_G

C******************************************************

C           Расчёт суммарных моментов                                              *

C******************************************************

 _XM11=_XM07+(((_CMXW*^XMWX*^XMXM*_WXR)+^XMWY*_WYR)*_BL2V)

      _YM12=_YM07+(^YMWY*^YMY*_WYR+^YMWX*_WXR)*_BL2V

      _YM14=((_PR-_PL)*_ZDV)

      _YM15=_YM14/_QSL

 _YM11=_YM12+_YM15

      _WXE=(_WXR*COS(_AR)-_WYR*SIN(_AR))*57.3

      _WYE=(_WYR*COS(_AR)+_WXR*SIN(_AR))*57.3

      _GAMA=_GAM*COS(_AR)

      _XM=_XM11*_QSL+_XMSP

      _YM=_YM11*_QSL+YMSP

      _ZN=_CZ*_QSG-_ZNSP