QJ 1029-1986 弹道式导弹水下弹道计算.pdf

1、Q.J 中华人民共和噩航天工业雷部标准QJ 1 029-86 弹道式导弹水下弹道计算1986 -08-1 a发布1987-05-01实施中华人民共和国戴天工业部批准中华人民共和国航天工业部部标准(1Jl029-86 弹道式导弹水下弹道计算本标准适用于水下潜艇发射或水于限定发射的弹道式导弹的水F弹i直讨算c其中，簿肉分段只造什于.单筒采用适配器作为配合结构的情况I 术语、符号潜能纵摇角(Js:潜艇级验。王军在发射坐标系垂直匮XOY(射面)上的投影与(-OX)寄自之夹角.OX，投影在(-OX!之上先jf-潜慈偏航角在潜艇纵被OX，与XOY二fjljj之夹角，OX，在射画之左为正。潜能横摇角r，艇体

2、唔主标酒X，Oy，与通过OX，轴的铅垂丽之夹角，OZ3在X，OZ平在言之下为正=导弹街仰角j-r为正c相对向你角11 S单体纵幸自OX，在艇体坐标面X，OY卡的投影与OY，辅之夹角.OX，投影领i句糙足对为正Q相对锅章主角手，弹体级辙。豆豆与X，OY军平面之夹窍.ox，在X，OY，之左为主。相对滚动角:弹体坐标铀(-OZ 1 )与x，oz，平面之夹角(-OZl)在X，OZ，之后为正3也射击方位都以jf北方向j.零度，顿时引转为正).度A l,A , 浪向(以正离方向为零度.I额防针转为主忡，度Ar一一滚动摩擦阻尼系数航天工业部1986-08-18发布1987-05-01实施1 2 QJ1029

3、-86 AI 流向4以正北方向为零度，顺时针转为正)，度Aw一风向以正离方向为零度，顺时针转为正)，度a一导弹质心至F适配器中点之距离，米B一适量E器抗压刚度，牛/米B 滚动压差阻尼系数b-发射街口至中适配器F端之距离，米CXl , CYJ CZ1 液体(或空气)铀向毛法i甸和销商力系数C旱，Ci液体附加力系数，1/度3一发射篱内径.米. e-Xs，在艇的纵对称面内，指向上为正苦。，Z，轴按右手直角坐标系确定。3 QJ 1029-86 2. 1.3 .量篝系。1X 1 Y I ZI 固连于弹体，原点01取在导弹质心上O. Xl袖沿弹体纵辙，指向头部z01Y，袖在弹体纵对称窗内，垂直于01X.辙

4、，指向上(p第1Il象跟线70 , z，辘安右手直角坐标系确定也2. 1.4 速度坐标系。.X, Y市ZE毒连于弹体，原点。1取在导弹质心上3()l X，辙沿导弹速度方向，O. Y.辙在导弹纵对称面肉，垂直于O，X穹轴，当速度矢量为水平时，指向上号。.Z.轴按右手直角坐标系确定。2.2 坐标变换2. 2. 1 发射坐标系与弹体坐标系的关系用矩阵表示(见图1) (XI)(飞Y. =A, - I Y Z ,/ Z) f司qQL O.哼lEzflE勾0.a.j I (1.1Q.jJ J r 盹3剧?=卜坷回归回f冲向Y f5:Y肯世勾在中酬步主伊f凶俨睛f刷Yf5马1占/f.马f4斟ff:;tf忖务

5、r剧f-$.时俨1 r-马Y1.-盔. (1) 制fCldI x,/! 因1发射在企标系与弹体坐标系的关系4 QJ 1 029-86 2.2.2 发徽坐标系与震体坐标泵的关系用矩萍表示见图2), ()=叫)(， b， b巧lAz_a, I 1., t,. .l-l I - l lI jl. J r -硝冉制于;=叫捆在-叫5仰视( - S.:.g,S.:.r. -(j( S:呢硝4几f主E 如唔晦结剧e，硝草+5，斗号军/;占4-约加S:.f忱f.如军飞4盼圳，一主酬皂剧r.I x 图2发射坐标系与童医体坐标系的关系2.2. 3 爱苦坐标系与速度垒挥系的关系用矩薄表示(见图3), 5 QJ 1

6、 029-86 G)=Zv ( ) (C 乌C号飞Az-!.v I乌C，Cj I 飞马1CJL 乙1$I f伊，cs(J 一1-时间碍于;.f; ，.斗8f1.tcef6五(TCI啄41-(6) 5乌8乡章制剧在ff，8时间仍es二.r.千元占19J锦4图3发射坐标系与速度坐标系的关系2.2.4 锺像坐标系与弹体坐标系的关系用矩阵表示(见图4) ()A， (iJ -S-r 1 . Yv I. (3 ) 耐副主/飞v1QJ 1 029-86 IJI 比J3飞A1.-芦，=14.1 ll I b adii J拎JJ ( -:贯副代=ft8多1，叫五句均句L叫件均内创仲、III-Jf章，仇，OV丁，

7、4HMrMl ry晤的44 ，缸，咱，哨你咣沪、如r叫叫响和加倒JkiehH叫J,-手t守ft主?( 主理4艇体坐标系与弹体坐标系的关系2. 2. 5 速度坐标系与弹体坐标系前关系用矩阵表示(见图5) G.)=叫:)7 /己刽制#=1-均由。仿古 ;乌卢f;二.0(.臼当琪。QJ 1 029-86 -但坦，:p，5hdF占81. - .(l) 田f ,/, Xl 回5速度坐标系与弹体坐标系的关系3运动程经根据导弹受力情况，水下弹道分为z篱中第一段事筒中第二段g水下毅和空中段。3 _ 1 篝Z争第一段运动方程经本段从导弹起动jlJ第二排适配器脱落为止之认为导弹与发射箭之间没离相对焦运动和横码及侧

8、向运动，严密地被约束在发射筒中c其运动方程为28 r产ofscosYs-女主约t，-A-hj乌产95- h, 9$ lwr，-是E司令炕叫Hjs十叭!乡=#F鸣叫争鸣用。!专=均产凡Mt-均巧rr吟g仇护+吨pfj-. -. (7) (K) (机吗J=al鸣叫冯f号，旨，主=鸟，峰1+司呼呼产也吃吗=.Yt?t!仲乳铲=8.阳川mj扫飞=吨%嗣必f唔号叶咆叫.) j 骂肚=t川唔号号+呜L!JP s c . S 毛=-穹可r，-书号f骂=吨.t电tlv土(除非吃产以-位夺走#=mGJ在QJ 1 029-串串9= r.白白f持)，O llT-QfCS乌恃)，主吼叫时讲荒草可=脱ELf-争)ffm

9、z;唯忐(11，川鸟飞J1. .(f) - (10) (11) . - . - -. (12) . (3) 9 QJ 1 029-86 用以下各式计算W，宝，y筐，Wz!，(/，礼Y，Vx!.Vy1 .Vz盒，的部值zi!x， k cos Y, - b.j, lwy:=其-1,95 1比，f，五叫+A,J / , i号=rcf;(:磊i中=-I) ?手写书号问，wtM-u, t约在扒qJ阶KOr-阳1白QJ 1 029-86 其=飞-aG生产飞马-(吻专吃fc)J，宁坞，亏电民-际，达-(，吗夺了咆r.).tj乌=守吗?以-(吟s今气与)JIt(，-吨EF际)d乌)4夕，王)J4=吨鸣，Ar鸣

10、，也-电J.j1-. -. - -(.2.Z) t=气r，均4广吆七吨JYr吟户均，d-4.J-吗J，- - - - -(.5) j町乌在没苦11Y，-纳=-军SYr-认.-.-.4J frHA-465 Wt，t冯唱非-抱A呜、e引铲LU弘臂品堪码:沟、组品吗4iRyg=马15._句号z叶叫苦与RZ1 -cJ.S，咄-ci，iJ呈S叫手吐鲁5号. - - . . - (3J 向1M.飞r.间号-只呼叫G巧(;.-A.)PEL 5412 字号叫均知约卢.lL!:1E二号已， (削)-豆子州-ALI-X.l) 县乡=主-主墨与p，.-q (1-毛只=哇t+号号耳饵阶叩-z认仲叫宁生;主=号令贵阳制

11、.)-:(.4;叫i忖翩叫f彷传叭占圳J冲句妒归叭川川护问叫呻-才飞叩tlr/d俨吨(11.龟f副晖蚓哽钊f优儿协为与)-a，品乒仇，均岛4A吗马?吨(，峭11M.川-xhJ15ihr守!我l乓7月闸V -有i内rR严严)令m?认斗是44P(小叫=5-，o哈哈;二飞irA民心时F咕叫呐耐(川CF原与(1.忻忑矿-忑)科叫吗片tMRJ，-F1F也jCJ;叮古的斗呐t， t /1ri, t F JF 1,. J 气喘;(A-AJ16 .、纠J. .?) . .IJ 年牛牛.(品。、，Jl J f _ . . (2) -(63) QJ 1029-86 度相剧痛姆UUA明指运或段HhQ Jl成=均1的

12、矶结ya引aAHaAMd吨，、oe事主仲件仲HU算陆凉，Lm计LIN品h也hs、民mw终斗川厅川(算!ltdew-aE二UH寸斗fe?咱/)段叫叫吨以1飞-V飞8中WMMTaVJA民手市奇左且tlui+-14?hv一哨AMV乌龟干冉afz:-AfE忡组初-ivLFMnzttLUL程算附件峙州内户户hAF俨VZ上本. . (64-) 6J 4 i十篝说胡4.1 计算初筐根据需要可按两种情况进行计算z一种从导弹起动开始计算g一科从弹尾出街开始计算。4.1.1 从弹起动开始计算防.其初始条件需给出以下参数数值:Xo , Y.2、A2耳、hl(l、Yw、Awo4.2方程组中常数的说明V /IiS. V

13、AS、fs、P由燃气发生器内弹道数据绘出播值衰。a苟岳、z。、Gq、ll.G、1.1，、1z.1，.L楠、s.x、XzB、K.f、D、D】圄哀、JxJ，、Jz、fxs、fz.f.1HS、，.玄，、H、r，、H，等也常数表给出。4.3 it体、空气动力来擎的说明17 QJ 1 029-86 军;:-, _=.十Cx盆、Cy、CZ1、mx，、my，、mz，. mY1、E221、C，二、c豆、4、X等由流体动力或空气动力系数表给出。Cy，、m:ll(CZi ，mYt)荡和1(1，61)查取，取a(，6)的符号cCX1 . C . Xy (C;t )按时叫1)查取，取正号。mz (m卖按悔和料i(Vf

14、A蚓#掷a1刽1查取，取a的符号。mx.由a，事事rtl查理呈，当rt在I.11象器重树，mXI 0 空中段的血泪按MJ.t查取，取a她符号3导弹自发射筒入水对，叮、A、B、F、XFCCY1cz宝、c、ttf盒、ml;l、mY1 mz宝、cmt、cmZ等在水系数表绘出，均按EzJZ-查取。导弹出水时，j-.A、B、F.XF. Cl电o智、杂可. C(Yc.宝、m又主、古町、m.m町、m., .F .Mn ，M峙，M，辜的子程序段cA. 2 .2. 4 SUBROUTINE DWV (MMS.DY) 计算卒y，.z.Vx，审y.T z. V1爪l，VZ1;盒，饨，叫空，wz室，Vz兔的子程序段c

15、A. 2 .2. 5 S t:BROUTINE DA YIN 打印各量约中间贮存子程序段。A. 2. 2. 6 SUBROUTINE AT (. A,K; 计算坐标变换矩萍A，D的子程序段C当K=o 1ft计算J生，K=1时计算B.DoA.2.2.7日t:BROUTINEATD (., AD, Sl,S 2 ) 20 QJ102曹86计算坐标变换短萍导数A，骂，D，的子程序段。81= 1. 82=IR去计算A81=-1，82=-1对计算B81= 1. 82=-1时计算DoA. 2. 2. 8 8UBROUTINE OPG (T ) 根据时间T计算比兔弘件，.，J,.r s.f s，及?啤号子程

16、序段。A, 2. 2. 9 l1l帘的描值.求导及积分子程序LAGR12 一元等距捶筐.阶导数及二阶导数。LAGR13 一无不等I!插筐。LAQ 元或组插值。LAGR22 二元拯筐。8P3 样条插值及求导sRKI Runge-Kutta法积分。A.2.3 打印符号与任务书符号对黑打印符号任务书符号x X 王Y Z z X. X y. Y z z G:I辈Cv Q q W v y FI. 1各P8. 非GM I VXI VXl VYI VYI V ZI 飞lZlVXI. 飞可= 0 ，插值表从古.20开始造p可。反举CYl为例，其他各量烈1比办理。A.2.4.3 关于A.B.D矩阵的导戴矩肆.B

17、.b若假设22 QJ 1 029-86 (句)A= I al a，_生:(.1 司Jal) =知b.l-b) 年弘p=! J. J/,- JH 经计算得(-，如+刊俨51户问弘T斗硝护A= 1-士也十，/，(/orf卡Ylln机I作乌f4F问于几乎;t问气-钊J叭产Sf-YI夕由11叶aJfdP弘一户句cPJy-y马f认十元占45MSg41-钊3弘&古切需主=1马b.t-古LOr+Y，4，-41，1-*二4+，t古儿身材-r1 飞时才古Jc.riOf-均I-54才4/:号-弘t句叫一与L( 1,dL +抖均生II-t，J1硝贯古柄吭飞如= -1, t-;/:11咐!t74d加州生内J;li+j

18、) -11，刑!I句-队Ij，d1J-刘严?以Lt斗的可-l/A.2.4.4 iP.击.r的计算公式1= w,.:y t乌叫什(zY-咳，)弘铲23 A. 3 15IJI BYYQ 24 QJ 1 029-86 +=气叫一可以-I句尹=乌，+JS今中争中响?计算程序FORTRAN(4) BIAN YI CHENG XU DY MASTER WATR 苏FORTYIl YU，瞅盼税制俯REAL LAD，Ll，L2，LO，LP，J豆，JY，LPM，MDIMENSJON Y (25), DY(25),TV AYS(2份，CVA YS(28,4),CCXl(l7), . CCYI(I7), 也争CCY

19、WZ(17)，CCZWY(17), CMZl(1)， CMYl(l7), C羽YWY(l7)，CDP. (17), .C范X(24，6)，CllXG:1(24), CM XAF(6)，CLll(2日，CL22(2日，CI，33(25)，备CIA剑21)，CL.55(25).CL66位份，CL35(2剖，CL26恒的，CAP(25)，CBP(21)，也.CF(24 ),CXFH(2l),CXF(21), CLl1l (3在)，CL222俗的，CL333(3的，善CL444(30)，CL555(30)，CL6饰。份，CL355(30)，CL揭6(30)，CF(3份，普CX(3的，CAPP(3的，C

20、BPP(30).CKCX(21).CKC引26，CKMZ(2i)，.CKMX(2窍，CKMZWZ(27)，CKCYWZ(28).CKCYWH(2剖，CKXYZH刮35)，磐CKXM(35，份，CMXlI(衍，OIXA(日，CMXl衍，日.CAF(27l.CCXP(27)，. CCYP(27),CCY1P(2), CX Y P(27), ALWL(15，衍，CAL5(22)，CWL5(2窍，.TAU(匀，H(纱，LAD(窍，PHO(2)，AI(2)，STA(26)，A3，3)，衍，纱，. CFL(15).CWL(l5),SA(2),CA(2), PL(幻，TL(2)，:-.r(2日，C:-.r(

21、25)，蕃ARl(24毡，7)，AR2(240，7), AR3(240, 7), AR4(240, 7), AR5(240, 7), AR6(240,7), 善AR7(240.5)，AR8(24告，5)，P t:HLIC TV A YS,CV A YS，CCXl，CCYl，CCYWZ.CCZ轨，Y，CMZl，CHYl，*CMYWY, *CDP，C刻豆，CMX(;M，CMXAF，CLl1，CL22，CI，33，CL44，CI，55，CL66，CL35.* CL26,CAP,CBP, CF ,CX FH, C X F,CL1 1 1, CL222,CL333,CL444,CL555, 铸C1,66

22、6, C L355, C 1，却也CFH，CXB，CAPP，CPP，CKC豆，CKCY，CKlrZ，*CKlIX. 铃CKllZWZ，CKCYWZ，CKCY毛，VH，CKXYZH，CKX刻，CMXM，CMXA，*ClIXl,CAF, * CCXP，CCYP，CCYIP，CXYP，AL币L，CAL5，C矶iL5，CFL，C币rLSA，CA，*PL,TL, *Y，DY，TAU，狂，LAD.PHO，AI，A，F.VP,Q, BTP,AFP.XL, QJ 1 029-86 *面区，N，M，A毛，VW，OSM，PSM，GSM，OSO，PSG，GSO，VSO，可SD，在白，OS，PS*GS, OSD,PS

23、D,GSD, HS,AR1,AR2,AR3,AR4,AR5,AR6, AR7,AR8, . AA,BB,LP ,LO.L2, RXS.RZS,X G,XZ,S l1AX.GO,G,JX,JY, RHS, ROO, *.DD,FF, *LP剂，DMAX.CON，PY2，PY1，PY，VWX, V韵Z，CON1，CON2，CO交3DATA AA.BB，LP，L2，Lt，LO，RXSRZS，XG，XZ，SMA豆，G-白，G，JX，JY，. RHS，ROO，DD，FF，LP室，DMAX，CON，PY2，PY，PY1/等3.309，5.87,4.149, 10. 22, 10.62,10.407 .11

24、.6, 1.45，6.121，6.25兽，替1.5垂，9.80665，16459.853，427.674，9750.017，104.6293，0.125，1.49，善。.3，.028，1.4，57.29578，6.2831853，3.1415927，1.57白79631IJK=l XL=LO-XG M=G/GO READ(2告，905)COX1，CON2，CON3制记FORMAT(3FO.o) READ(20,I) TVAYS,CVAYS,CCXl, CCYl ，ccn可Z，CCZWY，CMZl，* CMY1, CMYWY, . CDP，Cl1X，CMXGM，CMXAF，CLll，CL22，C

25、L33，CL44，CL55，CL6毡，CL35，每CL26，CAP，CBP，CF，CXFH，CXF，CL111，CL222，CL333.CL44哇，CL555. CL666, CL355,CL266,CFB,CXB, CAP P, CBPP, CK CX, CKCY ，CK盖fZ，*CKM豆，.CKMZWZ，CKCY韵TZ，CKCYW置，C主运XYZH，CKXM，CMXI岖，CMXA，CMXl. CAF, CCXP,CCYP,CCYlP ,CXYP,ALWL, CAL5,CWL5 1 FORl1AT(242FO.OI纷8FO.0/243FO.O/246FO.0/180FO.O/215FO./昏

26、251FO.0!135FO.0/164F币.0)PAFSE 2 READ(20,3)VSTA,STA FORl斗室AT(27F市.告:IISTA=VSTA+O.5 WRITE(40,55) (STA(I),I=I,16) 55 FOR羽AT(! 1 / i 11/1/ ! 11总3豆，13H1NITIALV ALUEII 普35X，3HXO=，F军.435X，3HYO=，F哇，4/35X.3HZO=,F9.41 普33X，5HTAOl =，F9.4/35X，3HHl=，F9.4/33豆，5HLMD1=，F9，垂/*33豆，5HFIOl= ，F9.4!35豆，3日Al=，F9.4/33X，5廷T

27、A02=，F9.41善35豆，3HH2=，F9.4f33X，5日LMD2=, F9.4/3.lX，5日FI02=，F9.4/普35X，3HA2=，F9.4/33X.5HWLAL=，F百.4!35X.3H飞W=，F9.41磐35豆，3日A轨i=，F古.4)25 QJ102窜86IF(NSTA.NE.25) GOTO 57 WRITE(40,56) (ST A(窍，1= 17,25) 56 FOR班主T(34X.4HOSM= ，F9.4f34豆，4日PSM=，F9.4* 34X.4HGSM = ，F9.4/34豆，4HOSO=，F9.4!34豆，4HPSO=,F9.4! 铸34X，4HGSO=，F

28、9.4!34X，4HVS如=，F9.4!34X，4HVS.= ,F9.4 *35X,3HHO= ,F9,41 / 1111 11111) GOTO 59 57 WRITE(40，5剖(STA(I)，I=17,26) 58 FORMAT(33X，5HVXIO=;F9.4!33豆.5HVYlO=.F9.419 26 普33五，5HVZIO=，F9.垂134X，4HFIO=，F9.在i34X，4HPSO=，F9.4j * 34X，4HGMO=，F9.4/33X.5HOXIO=，F9.4/33豆，5HOYIO=.F9.4/* 33X,5HOZIO= .F9.435豆.3HHO=，F9.411111!1

29、1/ /) y =0. AR1=0. AR2=怠AR3=往AR4=非，AR=O. AR岳=0.AR7=0. AR8=在.Y(ll)=STA(l) Y(l2), Y(l4)= ST生(2)Y (l3)=STA(3) TAU(I)=STA(4) TAU(2)=STA(9) H(1 )=STA(5) H(2)=STA(l0) LAD(l)=STA(6) LAD(2)=STA(1日PHO(1)=STA(7)/CON PHO(2)=STA(12)/CON AI(l)=STA位)/CONAI(2) =STA(13)/CON 实K=STA(14)+.5DO 81 =1,2 PL(I)=PY2/LAD(I)

30、TL(I)=PY2 /TAU(I) CA(I) = COS(AI(I) 8 SA(I)=SIN(AI(I) IF(NK.EQ.O)GOTO 10 IF(技K.GT.4)GOTO 10 DO 91=1,15 QJ 1 029-66 CFL()=ALWL(I,2. NK- l) 9 CWL(I)=ALWL役，2.NK)10飞W=STA(15)AW =STA(16)/CON vwx=-V飞可.COS(PY-A矶)VWZ=VW养SI去(PY-AW)IF(NST A.:-I E. 25) GOTO 4 08M =STA(l7)/CON P8M =STA(l8)/CON G8M =STA(l9)/CON

31、OSO=STA(20) PSc =8TA(21) GS =STA(22) VSo =STA(23) VSD=STA(24) HO =STA(25) CALL OPG(Y(1) CALL AT(OS.PS.GS,CTS,STS,CPS,SPS,CGS,SGS,B,I) Y(5)=PSD善CGS- B(2,3) . OSD Y(6)=GSD-B仪，3)祷。SDY(7)=PSD.SGS+B(3，刻.OSDY(9)=ASI去(-B位，3)CP=COS(Y(9) Y(B)=ASIN(B(2，到/CP)Y(lO)=ASI交(B(l，3)/CP)Y(2)=-B(2，1)铸VSO-Y(6)铸RZS-Y(7).

32、 RXS Y(3)=-B(1，。每VSO斗争Y(5).RZS-Y仔善.21Y(4)=B(3,1)*VSO + Y(6)*.21+ Y(5)* RXS lI = 1 AFP,BTP,Q,VP,F,HS=O. GOTO 5 垂Y(2) = STA(17) 可(3)= ST A(18) 27 Y(4)=STA(l9) Y(日=STA(23)jCONY(6)=STA(24)jCON Y(7)=STA(25)/CON Y彷)=STA(20)/CONY份)=STA(21)/CONY(l创=STA(22)jCON交=3QJ 1029-86 5 CALL AT(Y(B), Y(9), Y(10),CFF,SF

33、,CP,SP,CG,SG,A,O) HH=，创CALL RKl(25，HH，O，宜，DY，BN，CN)II YB=Y(12)-BB赞A(I，2)YD=Y(12)-L1*A役，2)YE=HS YF=DY(l2) CALL RK1(25.在H，l，Y，DY，BN，CN)COE=100 IF(Y(l2).GT.50.) COE=20. IF(ABS(COE* Y(1)-IFIX(COE祷Y(l)+O.OO5).GT.0.005)GOTO 50 CALL DAYI买50 IF(N.NE.l) GOTO 16 YA=Y(12)-BB椿A(1，2)IF(ABS(Y A).LT.Ol) GOTO 14 IF

34、(YA祷YB)12,12,11 12 HH=-HH/10. GOTO 11 14 CALL DA YIN HH =.01 * (IFIX(Y(l)铸l告0.)+1) - Y(ll N=2 CALL RKl(25,HH.l,Y,DY, BN,CN) CALL DAYIN HH=.Ol GOTO 11 16 TFN.NE.2) GOTO 20 YC= Y(l2)-Ll骨A(1，2)28 IF(ABS(YC).LT.Ol) GOTO 18 IF(YC* YD) 17,17,11 QJ 1029-86 17廷H=-HH/I0.GOTO 11 18 Y(到=Y(2)+2.27CALL DAYI只?面=3

35、HH =.01 * (IFIX(￥(1) * 100.)+ 1)-Y(I) CALL RK1(25，HH，O，宫，DY，BN，CN)CALL RK1(25,HH,I,Y,DY,BN,CN ) CALL DAYI到HH=.OI GOTO 11 20 IF(N.NE.3) GOTO 24 IF(ABS(HS).LT.Ol)GOTO 22 IF(HS备YE)21,21,11 21 HH=-HH/I0. GOTO 11 22 CALL DA YDr N=4 F,HS=O. HH =.oI * (IFIX(Y(1)善100.)+1)-Y(I) CALL RK1(25,HH,l ,Y.DY,BN,CN)

36、CALL DAYIN HH=.OI GO011 24 IF(ABS(DY(12).LT.Ol) COTO 2埠IF(DY(12) * YF) 25,25,11 25 HH=-HH/1O. GOTO 11 26 CALL DA YIN DO 130 1=1,3 12=40*(Z*I-l) 11=12-39 WRITE(4告，31)(ARl(K,L),L= 1,7),(AR1(K +40,L),L=1,),K = 11,12) WRITE(40,32) (AR2(K,L),L= 1，7)，(AR2(区+40，L)，L=I，)，K= 11,12) 矶lRITE(4告，33)(AR3(K,L),L=

37、1, 7),(AR3(K +40,L), L= 1,7), K = 11, (2) 巩IRITE(4串，34)(AR4(K，L)，L=1，7)，(AR垂(K+40,L),L= 1. 7).K =11.12) 哥市RITE但告，35)(AR5(K，L)，L=I，7)，(AR到K+40，L)，L=I，7)，豆=Il，I2)29 QJ 1029-66 WRITE(40,36) (AR6(区，L)，L = 1, 7),(AR6(K + 40,L),L = 1, 7),K =口，12)130飞，vRITE(40，38)(AR8(茧，L)，L=I，5)，(AR民K+ 40, L),L= 1,5),K= 1

38、1,12) 31 FORMAT(! I j I j jllX，1HT，9豆，lHX，9X，1HY,10X, 1 HZ，8X，2豆豆.RX，*2日Y刊* 8X，2HZ.，15X，IHT，9豆，IH豆，9豆，IHY，10X，lHZ，8X，2豆豆.，RX，2HY.，8X，铃2日Z，* !40(lFI5,4, FI0.4,Fll.4,4FI0.4,FI6.4,FI0.4,FIL4,4FI0.4)j j I11 /!) 32 FORMAT(! i /11 IllX，IHT，8X，3HG:vIC，8X，IHQ，9X，2HWY，7豆，4HFI.，* 6X,4H PS.，6豆，4HGM.，14X，IHT，8X

39、，3HGMC，8X，lHQ，9豆，2HWY，*7豆，4 HFI.，6X，4HPS.，6豆，4HGM.!.40(jF15.4，FI0.4，Fl1.3，Fl0.4，3Fl告.3，FI6.岳，FI0A，FIL3，Fl告.4，普3FlO.3)ill!J!ll 33 FORMAT(!fllf/12X，lHT，8X，3HVX1，7X，3HVY1，7豆，3HVZ1，6X，4HVX1.，9号6X，4HVYl.，6豆，4H飞Zl.，14豆，IHT，8X，3日VX1，7X，3HVY1，7X，3H可Zl，.6X，4HVX1.，6X，4H1.，6豆，4H飞iZl.1. 40(!FI6.4，3Fl白.也，1FI0.3，

40、F16.4，3FI0，4，3FI0.3)/11 1 I / f) 34 FOR览AT(!/ 1 / I /12豆，IHT，8豆，2HFI，8X，2HPS，8豆，2HGM，7豆，3HFl.，4舍7X，3HPS.7豆，3HGM.，14豆，1廷T，8X，2豆FI，8X，2HPS，8豆，2HGM，7X，3HFI.，* 7X,3HPS.,7X,3HGM./ 善40(lF16.4，3FI0.4，3FlO.3，F16.4,3FI0.4,3FIO.3)111 / / 11) 35 FORMAT(/ 1111 /12豆，lHT，7X，4HOMX1，6X，4HOMY1，毡X，4HOMZ1，普豆，5日OMX1.，5

41、X，5H 0 M Yl.，5豆，5HO双ZI.，13X，lHT，7豆，4HOMXl，普岳豆，4HOllYl，6X，4亘OMZ1，6豆，5HOMXL，5X，5HOlIYl.，5X，5HOMZL/祷40(1FI6.4,6FI0.3,FI6.4,6FIO.3)111/1 /) 36 FORMAT(j I i 11112豆，1HT ,8X,2HAF,8X,3HAFl,7X,2HBT,8X,3HBTl, .8豆，lHO，8X，3HSEG，14X，lHT，8X，2HAF，8X，3HAF1，X，2HBT，8豆，铸3HBTl，8豆，1日0，8X，3HSEG/* 4O(iF16.4,6FI0.4,FI6.4,6

42、F10.4)111! I! /) 38 FORMAT(/ / 111117豆，lHT，l1X，IHV，12豆，2HVl，l1X，lHF，12X，2HHS，*21X，1日T，ll豆，IHV，12X，2HV1，l1X，lHF，12X，2HHSI善40(lF2L4，2FI3.4，FI3.2，F13.4，F22.4，2F13.4，F13.2，F13.4lllllll/l IF(IJK.EG.13) GOTO 913 IJK=IJK+l GOTO 2 913 STOP END 30 QJ 1029- 86 SUROUTINE DERY(MM,Y,DY) REAL LO ,LP,LPM,JX,JY, KK

43、,KCX, MACH, MRX1, MRY1, MRZ1, *MX，l岖，L2，祷L11，L22，L33，王A4，L55，I.66，MLX1，MLYl，MLZI，LFDIMENSION YUIM)，Dycn日，A(3，3)，AD(3，刻，B(3，剖，BD衍，坊，D口，剑，替TD(3，3)，TV A YS(2柿，CVA YS(28,4),PV(4),CCX1(l7),CCXP(27),CCYP *(27), * CCYlP(27)，CAF(27)，CMXM传)，CMXA(时，CMXl(6，5)DIMENSTO买CXYP(27)PUBLIC L11, L22, L33,L44, L55, L66,

44、 CXYP, A,AD PUBLIC CCX1,TV AYS,CV AYS,CCXP, CCYP,CCY1P, CAF,CMXM, * CMXA, CMXl, *MK，窍，M，OS，PS，GS，OSD，PSD，GSD，OSDD，PSDD，GSDD，昏FLX1，FLY1,FLZ1, R 11, RZl, MLXl, ML Y1, MLZl, MRXl,l1RYI, *MRZ1,LF, * KCX, VSO, VSD, VP , F,Q, CP,SP,CG ,SG, *AA，LO，L2，LP，RXS，RZS，XZ，SMA豆，GO，G，JX，JY，ROO营养DD，FF，LPM，PY1，PY,CON

45、,AAFP,QSM, V飞鸟7豆.，.可iVZ，AFP，BTPAMG(VV,SS)= VV +8S/9.8秘65DY=O. DY(l)=l. CALL AT(Y(8), Y(9), Y(lO)，CF，SF，CP，SP，CG，SG，A，的DY(8)=(Y(6)* 8G+ Y(7)骨CG)/CPDY(9)=Y(6)祷CG-1(7)*8G DY(l创=Y(5)+8P，导DY(8)C ALL A TD( CF,8F, CP ,8P, CG,SG, DY (8), D1 (9), D1 (10),A, AD .1., L) DO 41 1=1,3 DO垂1J=I,3 41 DY(I+IO)=DY(l+l

46、+ A(J,U *Y(J + 1) DY(4)=Y(2) D1(15)= 1(3) D 1(16)= 1(4) IF(N.GT.3) GOTO 30 VS=V80+VSD骨Y(1)CALL LAQ(28，也Y(1 ),TVAY8,CV AYS,PV) CALL OPG(Y(l) CALL AT(OS,P8, G8 ，CS，ST8，CPS，SP鼠CGS，SGS，B,!) 31 QJ 1029-86 CALL ATD(CTS,STS,CPS,SPS,CGS,SGS,OSD,PSD,GSD,B,BD, -1., -1.) OXSD=GSDD-B(I，3)赞OSDD- BD(1,3) * OSD OY

47、SD= PSDD * CGS-PSD普GSD*SGS-B(2，窍*OSDD-BD(2，3)铸OS号。ZSD=- PSDD * SGS-PSD * GSD善CGS-B(33)*OSDD-BD(3，3)棒*OSDIF(N.GT.l) GOTO 10 DY(5)=OYSD DY(6)=OXSD DY(7)= -OZSD l=Y(5)祷Y(5)D2=Y(6)* Y(创D3= Y(7)* Y(7) DY(2)=Y(5)祷(Y(6)*RXS-Y(7)*RZS)-PV(3)普(D21卡D3)-棒B(2，I)* VSD-BD(2,1)* VS-DY(6) * RZS-DY(7)*RXS+PV(l) DY(3)=2.铸Y(7)*PV(2)+ Y(6)警(Y(5)铸PV(3)-Y(7)* RZS)-RXS善(Dl善+D3-BI1.1)* VSD-Bnll.I)骨VS+DY(5)