1 SUBROUTINE SUBPAH 2C 3C THIS ROUTINE YIELDS THE FOLLOWING WING DERIVATIVES- 4C SUBSONIC CLQ 5C CMQ 6C CLADOT 7C CMADOT 8C 9 COMMON /OVERLY/ NLOG,NMACH,IM,NALPHA 10 COMMON /CONSNT/ PI,DR,UNUSED,RAD 11 COMMON /OPTION/ SR,CBARR,RUFF,BLREF 12 COMMON /HTDATA/ A(195) 13 COMMON /HTI/ WINGIN(154) 14 COMMON /FLGTCD/ FLC(93) 15 COMMON /SYNTSS/ SYNA(19) 16 COMMON /BDATA/ XBD(300),DYN(213) 17 COMMON /IHT/ PW,WING(380) 18 COMMON /SBETA/ XST(243),TRA(108),STB(135) 19 COMMON /WHAERO/ XCD(106),C(51) 20 COMMON /EXPER/ KEXP(105),KWING 21 COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 22 1 HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN, 23 2 SUPERS,SUBSON,TRANSN,HYPERS, 24 3 SYMFP,ASYFP,TRIMC,TRIM 25 LOGICAL FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 26 1 HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN, 27 2 SUPERS,SUBSON,TRANSN,HYPERS, 28 3 SYMFP,ASYFP,TRIMC,TRIM 29 LOGICAL KWING 30 REAL MACH,LAMDA,MACOE 31 DIMENSION ROUTID(2) 32 DIMENSION WTYPE(4) 33 DIMENSION X21416(8),Y71416(8),Q71416(3) 34 DIMENSION X21428(9),Y71428(9) 35 DIMENSION I21416(4),I21428(4) 36 DIMENSION Q71428(3) 37 EQUIVALENCE(WING(101),CLA),(WING(121),CMA),(DYN(21),DCMCLO) 38 EQUIVALENCE(COSLC4,A(67)),(TANLC4,A(68)),(ASTRW,A(7)) 39 1 ,(SW,A(3)),(MACOE,A(16)),(CRSTR,A(10)) 40 2,(LAMDA,A(27)),(DXCG,A(173)) 41 EQUIVALENCE(CLQWT,WING(201)),(CLADWT,WING(241)) 42 1,(CMQMFB,DYN(1)) 43 EQUIVALENCE(CLQWSB,WING(201)),(CMQ2,DYN(2)),(CMQWSB,WING(221)) 44 1,(CLG,DYN(4)),(CLADWL,WING(241)),(CMOG,DYN(6)),(CMADPP,DYN(7)) 45 2,(CMADWL,WING(261)) 46 DATA WTYPE /4HSTRA,4HDOUB,4HCRAN,4HCURV/ 47 DATA ROUTID /4HSUBP,4HAW / 48 DATA Q71416 /4H7.1.,4H4.1-,4H6 / 49 DATA X21416 /0.,.25,.5,1.5,2.25,2.75,3.25,4./ 50 DATA Y71416 /0.,.005,.014,.065,.095,.112,.125,.14/ 51 DATA X11416 /1./ 52 DATA I21416 /8,3*0/ 53 DATA I21428 /9,3*0/ 54 DATA Q71428 /4H7.1.,4H4.2-,4H8 / 55 DATA X21428 /0.,.25,.5,1.25,1.75,2.25,2.75,3.25,4./ 56 DATA Y71428 /0.,.003,.008,.033,.047,.057,.066,.072,.077/ 57C 58 TYPE=WINGIN(15) 59 IF(TRANSN.AND.(TYPE.NE.WTYPE(1)))GO TO 1020 60 MACH=FLC(IM+2) 61 XACCR=C(6) 62 DCMDCL=-DCMCLO*CBARR/MACOE 63 DCMCLQ=-C(7)*CBARR/MACOE 64 IF(TRANSN)DCMCLQ=(TRA(105)*CRSTR-(SYNA(1)-SYNA(2)))/MACOE 65 IF(KWING)DCMCLQ=-(CMA/CLA)*CBARR/MACOE 66 IF(TRANSN)DCMDCL=(TRA(95)*CBARR-DXCG)/MACOE 67 CCLAL=WING(101) 68 CLASM=WINGIN(21)*SQRT(ABS(1.-FLC(3)**2)) 69 IF(TRANSN)CLASM=WINGIN(69) 70C 71C ---- WING PITCHING DERIVATIVES,CLQ,CMQ ---- 72C 73 CLQWSB=(.5+2.*DCMCLQ)*CCLAL*MACOE/CBARR 74 IF(TRANSN)CLQWT=CLQWSB 75 IF(TRANSN)MACH=TRA(6) 76 CMQ2=-0.7*CLASM*COSLC4*(ASTRW*(.5*DCMDCL+2.*DCMDCL**2)/(ASTRW+2.* 77 1 COSLC4)+ASTRW**3*TANLC4**2/((ASTRW+6.*COSLC4)*24.)+.125)*SW/SR 78 2 *(MACOE/CBARR)**2 79 IF(MACH.GT.0.2)GO TO 1000 80 CMQWSB=CMQ2 81 GO TO 1010 82 1000 BEE=SQRT(1.-MACH**2*COSLC4**2) 83 SAVE=ASTRW**3*TANLC4**2 84 CMQWSB=((SAVE/(ASTRW*BEE+6.*COSLC4)+3./BEE)/(SAVE/(ASTRW+6.*COSLC4 85 1 )+3.))*CMQ2 86 IF(TRANSN)CMQMFB=CMQWSB 87 IF(TRANSN)MACH=FLC(IM+2) 88 IF(MACH.GT.1.)GO TO 1020 89C 90C ---- WING ACCELERATION DERIVATIVES,CLAD,CMAD ---- 91C 92 1010 IF(LAMDA.NE.0.)GO TO 1020 93 BETA=SQRT(1.-MACH**2) 94 TEST=BETA*ASTRW 95 IF(TEST.GT.4.)GO TO 1020 96 IF(TRANSN)XACCR=TRA(105) 97C 98C -- FIGURE 7.1.4.1-6 -- 99C 100 CALL INTERX(1,X21416,TEST,I21416,Y71416,BPIAO2,8,8, 101 1 0,0,0,0,0,0,0,0,Q71416,3,ROUTID) 102 CLG=BPIAO2*PI*ASTRW/(-2.*BETA**2) 103 EXPID=SW*MACOE/(SR*CBARR) 104 CLADWL=(1.5*XACCR*CCLAL*SR/SW+CLG/19.1 )*EXPID 105 IF(TRANSN)CLADWT=CLADWL 106C 107C -- FIGURE 7.1.4.2-8 -- 108C 109 CALL INTERX(1,X21428,TEST,I21428,Y71428,BCMPAW,9,9, 110 1 0,0,0,0,0,0,0,0,Q71428,3,ROUTID) 111 CMOG=BCMPAW*ASTRW*0.5*PI/BETA**2 112 CMADPP=(-2.53125*XACCR**2*CCLAL*SR/SW+.078532*CMOG)*EXPID 113 1 *MACOE/CBARR 114 CMADWL=CMADPP+(DXCG/CBARR)*CLADWL 115 1020 CONTINUE 116 RETURN 117 END 118