Program dM6
C—————6自由度弹丸外弹道飞行计算程序—————
C= = = 积分变量注释
C Y(1)————————————————弹丸飞行时间 T
C Y(2)————————————————弹丸飞行速度 V
C Y(3)————————————————弹道倾角 Q
C Y(4)————————————————质心的X轴坐标 X
C Y(5)————————————————质心的Y轴坐标 Y
C Y(6)————————————————质心的Z轴坐标 Z
C Y(7)————————————————弹丸自转角 r
C Y(8)————————————————弹丸自转角速度 dr/dt
C Y(9)————————————————偏角 Pasi1
C Y(10)————————————————偏角 Pasi2
C Y(11)————————————————摆动角 fai1
C Y(12)————————————————摆动角 fai2
C Y(13)————————————————摆动角速度 d(fai1)/dt
C Y(14)————————————————摆动角速度 d(fai2)/dt
C = = = = = = = = = = = = = =
C Pasi=Pasi1+i*Pasi2
C Fai=fai1+i*fai2
C = = = = =* * * * * * * = = = = = = * * * * * * * = = = = = =
REAL M0,Mzz0,My0,Kz0,Ia,Ic
Real Mzzb,Myb,Mxzb
DIMENSION Y(14),DY(14),YA(14),YB(14),Y1(14),Y2(14)
DIMENSION Y3(14)
COMMON/SHU/MMM
COMMON/WBB/I1,I3,I0
COMMON/WBC/N1,SXm(12),Cx0b(12),Cnb(12),Xpb(12)
COMMON/WBD/N2,DXm(12),Czb(12),Mzzb(12),Myb(12),Mxzb(12)
COMMON/WBE/M0,Ic,Ia,Cd0k,S,D,NH,Xcg
COMMON/WBF/Sg,Sd,Cx,Cy0,Mzz0,My0,Kz0
COMMON/WBP/FIP0,PSI0
C = = = = = = = = = = = = = = FORMAT-01 BEGAIN = = = = = = = = = = = =
1000 FORMAT(2(/),13X,8( ‘*’),2X, ‘TYPE OF PROGECTILE IS:’ ,
外弹道计算程序 来自淘豆网www.taodocs.com转载请标明出处.