1 SUBROUTINE JETPWE 2C 3C ----SUBROUTINE COMPUTES JET POWER EFFECTS 4C 5 COMMON /POWER/ XX(12),AIETLJ,NENGSJ,THSTCJ,JIALOC,JEVLOC,JEALOC, 6 1 JINLTA,JEANGL,JEVELO,AMBTMP,JESTMP,JELLOC,JETOTP, 7 2 AMBSTP,JERAD 8 COMMON /POWR/ ATP,DCLT(20),XBARIN,XINOCR,DEUDA,EPSLON,ATJ, 9 1 DCLNJ(20),XEP,ZJP,XJP,XHP,AIN,VIN,TINOTJ,VJPOVI, 10 2 ZJPORJ,DE,ZJPOBH,YTOB2H,DEBODE,ZJPXHP,SRTPCO, 11 3 ZJDEXH,COMP1,PTEOPI,RJPORJ,RJP,DXPORJ,DXP,XEPC, 12 4 XHPC,ZTP,ZJPRJP,DCLHE(20),ZBART,DCMT(20),XL, 13 5 DCMNJ(20),DLH,DCME(20) 14 COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 15 1 HEAD,PRPOWR,JETPOW,LOASRT 16 COMMON /CONSNT/ PI,DEG,UNUSED,RAD 17 COMMON /SYNTSS/ XCG,XW,ZW,ALIW,ZCG,XH,ZH,ALIH,XV,VERTUP,HINAX, 18 1 XVF,SCALE,ZV,ZVF,YV,YF,PHIV,PHIF 19 COMMON /OVERLY/ NLOG,NMACH,M,NALPHA 20 COMMON /FLGTCD/ FLC(93) 21 COMMON /OPTION/ SR,CBARR 22 COMMON /WINGD/ A(195),B(49) 23 COMMON /HTDATA/ AHT(195) 24 COMMON /WHAERO/ C(51),D(55) 25 COMMON /WINGI/ WINGIN(77) 26 COMMON /HTI/ HTIN(131) 27 COMMON /IWING/ PWING, WING(400) 28 COMMON /IHT/ PHT, HT(380) 29 COMMON /IDWASH/ PDWASH,DWASH(60) 30 COMMON /IPOWER/ PPOWER, DCDPOW(20),DCLPOW(20),DCMPOW(20) 31 LOGICAL VERTUP 32 LOGICAL FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 33 1 HEAD,PRPOWR,JETPOW,LOASRT 34 LOGICAL LOGCOM(11) 35 REAL MACH(20) 36 REAL NENGSJ,JIALOC,JEVLOC,JEALOC,JINLTA,JEANGL,JEVELO,JESTMP, 37 1 JELLOC,JETOTP,JERAD 38 EQUIVALENCE (CDL(1),D(36)),(CLW(1),WING(21)) 39 EQUIVALENCE (ALSCHD(1),FLC(23)),(MACH(1),FLC(3)) 40 EQUIVALENCE (A(120),AR),(ALPHA0,A(134)) 41 EQUIVALENCE (LOGCOM(1),FLTC) 42 EQUIVALENCE (BO2H,HTIN(4)),(CLALPH,HT(101)),(QHOQI(1),DWASH(1)), 43 1 (XBARRH,AHT(161)) 44 DIMENSION CLW(20),CDL(20) 45 DIMENSION ALSCHD(20) 46 DIMENSION ROUTID(2),Q44161(2),Q46128(2),Q46129(2),Q46131(2), 47 1 Q6132A(3),Q6132B(3),Q46114(2) 48 DIMENSION QHOQI(20) 49 DIMENSION X16114(6),X26114(8),Y46114(48) 50 DIMENSION X16117(13),X26117(4),Y46117(52) 51 DIMENSION X6118A(8),Y6118A(8) , C6118A(6) 52 DIMENSION X6118B(8),Y6118B(8),C6118B(6) 53 DIMENSION X6119A(9),Y6119A(9),C6119A(6) 54C 55C ----TABLES FOR UPWASH GRADIENT AT PLANE OF SYMMETRY FOR UNSWEPT 56C ----WINGS. FIGURE 4.4.1-73 57C 58 DIMENSION X14161(4),X24161(10),Y44161(10,4) 59 DATA Q44161/4H4.4.,4H1-73/,Q46128/4H4.6.,4H1-28/,ROUTID/4HJETP, 60 1 4HWE /, Q46129/4H4.6.,4H1-29/,Q46131/4H4.6.,4H1-31/, 61 2 Q6132A/4H4.6.,4H1-32,4HA /, Q46114/4H4.6.,4H1-28/, 62 3 Q6132B/4H4.6.,4H1-32,4HB / 63 DATA X14161/4.,6.,9.,12./ 64 DATA X24161/.25,.4,.5,.6,.7,.8,1.0,1.2,1.6,2.0/ 65 DATA Y44161/1.08,.545,.40,.31,.24,.20,.13,.10,.06,.04 , 1.18,.680, 66 1.52,.40,.32,.27,.19,.15,.10,.08 , 1.30,.81,.62,.49,.40,.34,.25,.20 67 2,.13,.12 , 1.4,.88,.70,.56,.445,.39,.30,.24,.165,.14 / 68C 69C -----FIGURE 4.6.1-28 70C -----EFFECTIVE MEAN DOWNWASH ACTING ON HORIZONTAL TAIL 71C 72 DATA X16114/0.0,0.3,0.5,0.7,0.9,1.1/ 73 DATA X26114/0.0,0.2,0.4,0.6,0.8,1.0,1.2,1.4/ 74 DATA Y46114/0.0,0.54,0.765,0.860,0.915,0.950,0.969,0.988, 75 10.0,0.51,0.721,0.825,0.890,0.923,0.950,0.970, 76 20.0,0.450,0.675,0.790,0.860,0.900,0.922,0.945, 77 30.0,0.395,0.595,0.725,0.810,0.861,0.892,0.914, 78 40.0,0.300,0.522,0.660,0.748,0.810,0.859,0.890, 79 50.0,0.200,0.415,0.570,0.682,0.758,0.810,0.850/ 80C 81C -----FIGURE 4.6.1-29 82C -----VARIATION OF EQUIVALENT VELOCITY RATIO,WITH ACTUAL VELOCITY RA 83C 84 DATA X16117/0.46,0.42,0.38,0.34,0.32,0.30,0.28,0.26,0.24,0.22, 85 1.20,0.16,0.12/ 86 DATA X26117/1.0,6.0,14.0,24.0/ 87 DATA Y46117/1.0,4.25,9.66,16.6, 1.0,4.10,9.21,15.7, 1.0,3.90,8.80, 88 114.9, 1.0,3.85,8.35,14.10, 1.,3.6,8.1,13.7, 1.,3.55,7.89,13.3, 89 21.0,3.40,7.60,12.90, 1.0,3.30,7.40,12.40, 1.0,3.20,7.10,12.0, 90 31.0,3.10,6.80,11.40, 1.0,3.00,6.50,11.00, 1.0,2.72,5.90,9.85, 91 41.0,2.42,5.15,8.55/ 92C 93C -----FIGURE 4.6.1-31 94C -----DOWNWASH INCREMENT DUE TO A SUBSONIC JET IN A SUBSONIC STREAM 95C 96 DATA I6118A/0/ 97 DATA X6118A/0.0,0.4,0.8,1.2,1.6,2.0,2.4,2.6/ 98 DATA Y6118A/0.0,0.34,0.63,0.88,1.1,1.31,1.51,1.61/ 99C 100C -----FIGURE 4.6.1-32A 101C -----EQUIVALENT JET ORIFICE RADIUS FOR SUPERSONIC JET IN A SUBSONIC 102C -----STREAM 103C 104 DATA I6118B/0/ 105 DATA X6118B/3.0,4.0,6.0,8.0,10.0,12.0,14.0,16.0/ 106 DATA Y6118B/1.08,1.30,1.70,2.02,2.23,2.41,2.55,2.70/ 107C 108C -----FIGURE 4.6.132B 109C -----DOWNSTREAM DISPLACEMENT OF EQUIVALENT JET ORIFICE LOCATION 110C -----FOR A SUPERSONIC JET IN A SUBSONIC STREAM 111C 112 DATA I6119A/0/ 113 DATA X6119A/1.0,1.2,1.4,1.6,1.8,2.0,2.2,2.4,2.6/ 114 DATA Y6119A/0.0,1.0,1.6,2.07,2.40,2.65,2.86,3.04,3.16/ 115C 116 TMPINK(T)=273.+.556*(T-492.) 117 TPCNE=NENGSJ*THSTCJ 118 COSAIW=COS(DEG*ALIW) 119 IF(NENGSJ.EQ.1.)GO TO 1010 120 IF(WINGIN(2).EQ.UNUSED) GO TO 1010 121 IF(JELLOC .LE.WINGIN(4)-WINGIN(2)) GO TO 1000 122 CRJ=WINGIN(5)-(WINGIN(5)-WINGIN(1))*((JELLOC -(WINGIN(4) 123 1 -WINGIN(2)))/WINGIN(2)) 124 XBARJ=(XW+A(38)*(WINGIN(4)-WINGIN(2))+A(86)*(JELLOC -(WINGIN(4) 125 1 -WINGIN(2)))-JIALOC)*COS(DEG*ALIW)/CRJ 126 GO TO 1020 127 1000 WINGIN(1)=SAVE 128 WINGIN(1)=WINGIN(5) 129 1010 CRJ=WINGIN(6)-(WINGIN(6)-WINGIN(1))*JELLOC /WINGIN(4) 130 XBARJ=(XW+A(38)*JELLOC +CRJ/4.0-JIALOC)*COS(DEG*ALIW)/CRJ 131 1020 CONTINUE 132 SAVE=WINGIN(1) 133C 134C ----FIG 4.4.1-73 135C 136 CALL TLINEX(X14161,X24161,Y44161,4,10,AR,XBARJ,DEUDA, 137 1 2,2,2,2,Q44161,2,ROUTID) 138 AINNE2=2.0*JINLTA*NENGSJ/SR 139 SINAIT=SIN(DEG*AIETLJ) 140 COSAIT=COS(DEG*AIETLJ) 141C 142C ----GEOMETRIC COMPUTATIONS IF TAIL PRESENT. 143C 144 IF(.NOT.HTPL)GO TO 1090 145 COSAIH=COS(DEG*ALIH) 146 XEP=(XH+XBARRH*COSAIH-JEALOC) 147 TANTHJ=TAN(DEG*JEANGL) 148 ZT=JEVLOC+JERAD*SINAIT/TANTHJ 149 ZJP=XEP*SINAIT+(ZH-ZT)*COSAIT 150 XEP=XEP/COSAIT 151 XJP=4.6*JERAD 152 XHP=XJP+XEP 153 AIN=49.0*SQRT(AMBTMP) 154 VIN=MACH(M)*AIN 155 VJOVIN=JEVELO/VIN 156 XEPORJ=XEP/JERAD 157 ZJPOBH=ZJP/(2.0*BO2H) 158 YTOB2H=JELLOC/BO2H 159C 160C ----FIGURE 4.6.1-28 EFF. MEAN DOWNWASH ACTING ON HORT. TAIL 161C 162 CALL TLINEX(X16114,X26114,Y46114,6,8,YTOB2H,ZJPOBH,DEBODE, 163 1 0,0,2,2,Q46114,2,ROUTID) 164 KASE=3 165 IF(JEVELO.GE.AIN)GO TO 1030 166 KASE=1 167 IF(XEPORJ.GT.16.0)KASE=2 168 1030 GO TO (1040,1060,1070),KASE 169C 170C ----COMPUTE PART OF CASE 1 WHICH IS NOT A FUNCTION OF ALPHA 171C ----SUBSONIC FREE STREAM , SUBSONIC JET . HORIZONTAL TAIL .LT. 8 172C ----EXIT DIAMETERS DOWNSTREAM OF JET EXIT PLANE. 173C 174 1040 TINOTJ=TMPINK(AMBTMP)/TMPINK(JESTMP) 175C 176C ----FIGURE 4.6.1-29 VAR. OF EQUIV. VEL. RATIO TO ACTUAL VEL. RATIO 177C 178 CALL TLINEX(X16117,X26117,Y46117,13,4,TINOTJ,VJOVIN,VJPOVI, 179 1 2,0,2,1,Q46129,2,ROUTID) 180 IF(KASE.EQ.3) GO TO 1050 181 ZJPORJ=ZJP/JERAD 182C 183C ----FIGURE 4.6.1-31 A-C 184C ----FIGURE 4.6.1-31 DOWNWASH INCREMENT DUE TO A SUBSONIC JET IN A 185C ----SUBSONIC STREAM. 186C 187 1050 CALL FG6115(DE,XEPORJ,ZJPORJ,VJPOVI) 188 GO TO 1080 189C 190C ----COMPUTE PART OF CASE II WHICH IS NOT A FUNCTION OF ALPHA 191C ----SUBSONIC FREE STREAM SUBSONIC JET,HORIZONTAL TAIL LOCATED 192C ---MORE THAH 8 EXIT DIAMETERS DOWNSTREAM OF JET EXIT PLANE. 193C 194 1060 ZJPXHP=ZJP/XHP 195 SRTPCO=SR*THSTCJ/(XHP**2) 196C 197C ----FIGURE 4.6.1-31 DOWNWASH INCREM. SUBSON JET IN SUB. STREAM 198C 199 CALL TBFUNX(SRTPCO,ZJDEXH,Q,8,X6118A,Y6118A,C6118A,I6118A,MI,NG, 200 1 0,1,Q46131,2,ROUTID) 201 DE=ZJDEXH/ZJPXHP 202 GO TO 1080 203C 204C ----CASE III SUBSONIC FREE STREAM,SUPERSONIC JET,HORTZ TAIL IS 205C DOWNSTREAM OF THE FULLY EXPANDED FLOW. 206C 207 1070 PTEOPI=JETOTP/AMBSTP 208C 209C ----FIGURE 4.6.1-32A 210C 211 CALL TBFUNX(PTEOPI,RJPORJ,Q,8,X6118B,Y6118B,C6118B,I6118B,MI,NG, 212 1 2,2,Q6132A,3,ROUTID) 213 RJP=JERAD*RJPORJ 214C 215C ----FIGURE 4.6.1-32B 216C 217 CALL TBFUNX(RJPORJ,DXPORJ,Q,9,X6119A,Y6119A,C6119A,I6119A,MI,NG, 218 1 0,2,Q6132B,3,ROUTID) 219 DXP=DXPORJ*JERAD 220 XEPC=XEP+DXP 221 XHPC=XEPC+RJP/TANTHJ 222 ZTP=JEVLOC-DXP*SINAIT 223 ZJP=XHPC*SINAIT+COSAIT*(ZH-ZTP) 224 ZJPOBH=ZJP/(2.0*BO2H) 225 CALL TLINEX(X16114,X26114,Y46114,6,8,YTOB2H,ZJPOBH,DEBODE, 226 1 0,0,2,2,Q46128,2,ROUTID) 227 ZJPORJ=ZJP/RJP 228 XEPORJ=XEPC/RJP 229 GO TO 1040 230 1080 COMP1=CLALPH*DE*DEBODE 231C 232C ----COMPUTE JET POWER EFFECTS ON PITCHING MOMENT DUE TO OFFSET OF 233C ----THRUST AXIS 234C 235 1090 ZBART=COSAIT*(JEVLOC-ZCG)+SINAIT*(JEALOC-XCG) 236 DCMT(1)=-THSTCJ*ZBART*NENGSJ/CBARR 237 1100 CONTINUE 238 DO 1110 J=1,NALPHA 239C 240C ----COMPUTE LIFT COMPONENT DUE TO THRUST VECTOR (DCL)T 241C 242 ATP=ALSCHD(J)+AIETLJ 243 SINATP=SIN(DEG*ATP) 244 DCLT(J)=TPCNE*SINATP 245C 246C ----COMPUTE LIFT DUE TO INCLINATION OF THRUST LINE 247C 248 EPSLON=DEUDA*(ALIW+ALSCHD(J)-ALPHA0) 249 ATJ=ALSCHD(J)+AIETLJ+EPSLON 250 SINATJ=SIN(DEG*ATJ) 251 DCLNJ(J)=AINNE2*SINATJ 252C 253C ----COMPUTE LIFT DUE TO HORIZONTAL TAIL IF IT IS PRESENT 254C 255 DCLHE(J)=0.0 256 IF(HTPL)DCLHE(J)=COMP1*QHOQI(J) 257C 258C ----COMPUTE TOTAL LIFT DUE TO POWER 259C 260 DCLPOW(J)=DCLNJ(J)+DCLT(J)+DCLHE(J) 261C 262C ----COMPUTE JET POWER EFFECTS ON PITCHING MOMENT DUE TO NORMAL 263C ----FORCE ON THE ENGINE INLET. 264C 265 XL=JIALOC-XCG 266 DCMNJ(J)=-DCLNJ(J)*XL/CBARR 267C 268C ----COMPUTE JET POWER EFFECTS ON PITCHING MOMENT DUE TO THE CHANGE 269C ----IN ANGLE OF ATTACK OF THE HORIZONTAL TAIL. 270C 271 DLH=XH+XBARRH*COSAIH-XCG 272 DCME(J)=-DCLHE(J)*DLH/CBARR 273C 274C ----TOTAL JET POWER EFFECTS ON PITCHING MOMENT 275C 276 DCMPOW(J)=DCMT(1)+DCMNJ(J)+DCME(J) 277C 278C ------A SINGLE ENGINE JET CHANGES THE BASE PRESSURE AND THUS 279C ALTERS CDO HOWEVER NO METHODS ARE DEFINED FOR THIS CHANGE IN DRAG 280C -------JET EFFECTS DO NOT CAUSE A CHANGE IN INDUCED DRAG 281C BECAUSE THEY DO NOT MODIFY THE AERODYNAMIC LIFT OF THE WING 282C 283 DCDPOW(J) = 0.0 284 1110 CONTINUE 285 RETURN 286 END 287