1 SUBROUTINE SUPHB 2C 3C*** CALCULATES SUPERSONIC HORIZONTAL TAIL-BODY CL,CD,CLA AND CMA 4C 5 REAL MACH,NXX,KWB,KBW,IVBW ,KKBW,KKWB 6 DIMENSION ROUTID(2),Q1211A(3),Q1211B(3),Q31210(3),Q1212A(3), 7 1 Q1212B(3),Q2118A(3),Q2118B(3),Q2137A(3),Q2137B(3) 8 DIMENSION T4312B(8),D4312B(8) 9 DIMENSION CDW(20),CDB(20),WTYPE(4),IVBW(20),GAMMA(20) 10 DIMENSION T4337A(11),D4337A(24),T4337B(12),D4337B(20) 11 DIMENSION T4312A(11),D4312A(11) 12 DIMENSION TFIG10(11),DKWB10(11),DKBW10(11) 13 DIMENSION T4218A(18),DL218A(54),DR218A(54),T4218B(14),DL218B(42), 14 1 DR218B(42) 15 DIMENSION LGH(4),VAR(4),CD(20),CN(20),CA(20),CL(20),CLB(20), 16 1 CLW(20),ALPHAB(20),CDL(20) 17 DIMENSION CM(20) 18C 19 COMMON /CONSNT/ PI,DEG,UNUSED,RAD 20 COMMON /FLGTCD/ FLC(160) 21 COMMON /SYNTSS/ SYNA(19) 22 COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF 23 COMMON /HTI/ W(2),SPANS,SPAN,W1,CR,W2(8),TYPE 24 COMMON /BODYI/ NXX,XCOOR(20),S(20),P(20),R(20),ZU(20),ZL(20), 25 1 BNOSE,BTAIL,RLN 26 COMMON /BDATA/ BD(762) 27 COMMON /HTDATA/ AHT(195) 28 COMMON /OVERLY/ NLOG,NMACH,I,NALPHA 29 COMMON /SUPBOD/ SBD(229) 30 COMMON /SUPWBB/ SWB(61),SHB(61) 31 COMMON /SUPWH/ SLG(141), STG(141) 32 COMMON /IBODY/ PBODY, BODY(400) 33 COMMON /IHT/ PHT, HT(380) 34 COMMON /IBH/ PBH, BH(380) 35C 36 EQUIVALENCE (CA(1),BH(81)),(CDW(1),HT(1)),(CMA,BH(121)), 37 1(DXCG,AHT(173)),(CLAB,SBD(18)),(KKWB,SHB(2)), 38 2(XACN,SHB(3)),(CDOWB,SHB(4)),(DELXW,BD(89)),(ARSTAR,AHT(7)), 39 3(CD(1),BH(1)),(DD,SHB(5)),(XW,SYNA(6)),(CLB(1),BODY(21)), 40 4(ALIW,SYNA(8)),(BETA,SHB(6)),(CLAW,HT(101)), 41 5(XACW,STG(134)),(CLABW,SHB(7)) ,(XACBW,SHB(8)) 42 EQUIVALENCE (CLA,BH(101)), (CDL(1),STG(53)) 43 EQUIVALENCE (CL(1),BH(21)),(FA,SHB(9)),(CLI,SHB(10)), 44 1(KBW,SHB(11)),(CDOB,SBD(124)),(IVBW(1),SHB(12)),(RKBW,SHB(32)), 45 2(CLAWB,SHB(33)),(TANL ,AHT(62)),(CRSTAR,AHT(10)), 46 3(TAPEXP,AHT(27)),(CN(1),BH(61)),(FN, SHB(34)),(CDB(1),BODY(1)), 47 4(CLW(1),HT(21)),(KWB,SHB(35)),(XAC,SHB(36)), (CDOW,STG(80)), 48 5(DN,SBD(4)),(D1,SBD(5)),(KKBW,SHB(37)),(RLAP,SHB(38)) 49 EQUIVALENCE (XACA,SHB(39)),(GAMMA(1),SHB(40)), 50 1(TRINO,SHB(60)),(XCPLN,SHB(61)),(ALPHAB(1),BD(255)) 51 2,(CM(1),BH(41)) 52C 53 DATA ROUTID/4HSUPW,4HB /, 54 2 Q31210/4H4.3.,4H1.2-,4H10 /,Q1212A/4H4.3.,4H1.2-,4H12A /, 55 3 Q2118A/4H4.2.,4H2.1-,4H23A /,Q1212B/4H4.3.,4H1.2-,4H12B /, 56 4 Q2118B/4H4.2.,4H2.1-,4H23B /,Q2137A/4H4.3.,4H2.2-,4H37A /, 57 5 Q2137B/4H4.3.,4H2.2-,4H37B / 58 DATA WTYPE/4HSTRA,4HDOUB,4HCRAN,4HCURV/ 59C 60C FIGURE 4.3.1,2-10 KWB 61C 62 DATA TFIG10/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/ 63 DATA DKWB10/1.0,1.08,1.16,1.26,1.36,1.46,1.56,1.67,1.78,1.89,2.0/ 64C 65C FIGURE 4.3.1.2-10 KBW 66C 67 DATA DKBW10/0.0,.13,.29,.45,.62,.80,1.0,1.22,1.45,1.70,2.0/ 68C 69C FIGURE 4.2.2.1-23A (LEFT SIDE) 70C 71 DATA T4218A 72 1 / 0., 0.2, 0.4, 0.6, 0.8, 1.0 , 3*0.0, 73 2 0.0, 0.4, 0.8, 1.2, 1.6, 2.0, 3.0, 4.0, 5.0 / 74 DATA DL218A 75 1 / 0.543, 0.542, 0.541, 0.540, 0.534, 0.526, 76 2 0.400, 0.409, 0.418, 0.430, 0.441, 0.450, 77 3 0.305, 0.328, 0.350, 0.369, 0.387, 0.400, 78 4 0.238, 0.265, 0.295, 0.318, 0.339, 0.356, 79 5 0.198, 0.221, 0.246, 0.274, 0.298, 0.320, 80 6 0.160, 0.185, 0.210, 0.239, 0.262, 0.288, 81 7 0.065, 0.095, 0.122, 0.150, 0.177, 0.210, 82 8 0.000, 0.005, 0.035, 0.062, 0.089, 0.130, 83 9 0.000, 0.000, 0.000, 0.000, 0.002, 0.050 / 84C 85C FIGURE 4.2.2.1-23A (RIGHT SIDE) 86C 87 DATA DR218A 88 1 / 0.445, 0.464, 0.485, 0.500, 0.518, 0.526, 89 2 0.448, 0.455, 0.460, 0.460, 0.459, 0.450, 90 3 0.460, 0.449, 0.438, 0.424, 0.412, 0.400, 91 4 0.450, 0.430, 0.412, 0.394, 0.375, 0.356, 92 5 0.432, 0.410, 0.388, 0.365, 0.343, 0.320, 93 6 0.420, 0.394, 0.369, 0.340, 0.314, 0.288 , 94 7 0.388, 0.354, 0.322, 0.278, 0.244, 0.210, 95 8 0.357, 0.314, 0.273, 0.216, 0.171, 0.130, 96 9 0.325, 0.274, 0.225, 0.154, 0.100, 0.050 / 97C 98C FIGURE 4.2.2.1-23B (LEFT SIDE) 99C 100 DATA T4218B 101 1 / 0., 0.2, 0.4, 0.6, 0.8, 1.0 , 0., 102 2 0.0, 0.5, 1.0, 2.0, 3.0, 4.0, 5.0 / 103 DATA DL218B 104 1 / .665, .665, .665, .665, .665, .665, 105 2 .425, .492, .539, .550, .550, .550, 106 3 .330, .370, .405, .438, .459, .470, 107 4 .184, .215, .250, .284, .318, .350, 108 5 .060, .097, .133, .170, .206, .240, 109 6 .000, .000, .044, .083, .127, .170, 110 7 .000, .000, .000, .020, .063, .105 / 111C 112C FIGURE 4.2.2.1-23B (RIGHT SIDE) 113C 114 DATA DR218B 115 1 / .665, .665, .665, .665, .665, .665, 116 2 .480, .500, .519, .536, .546, .550, 117 3 .338, .388, .430, .458, .471, .470, 118 4 .338, .372, .394, .394, .375, .350, 119 5 .410, .375, .341, .308, .272, .240, 120 6 .377, .338, .294, .251, .211, .170, 121 7 .342, .300, .246, .194, .146, .100 / 122C 123C FIGURE 4.3.1.2-12A (KBW) 124C 125 DATA T4312A/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/ 126 DATA D4312A/0.0,.11,.21,.31,.41,.51,.6,.7,.8,.9,1.0/ 127C 128C FIGURE 4.3.1.2-12B (KWB) 129C 130 DATA T4312B/.015,.1,.2,.3,.4,.6,.8,.975/ 131 DATA D4312B/1.,.975,.956,.947,.941,.950,.978,1.0/ 132C 133C FIGURE 4.3.2.2-37A 134C 135 DATA T4337A/0.0,.4,.8,1.2,1.6,2.,2.4,2.8, 136 1 .1,1.0,999999./ 137 DATA D4337A/ 138 1.5,.72,.900,1.08,1.24,1.39,1.53,1.68, 139 2.5,.72,.910,1.09,1.25,1.41,1.57,1.72, 140 3.5,.73,.920,1.11,1.27,1.43,1.59,1.74/ 141C 142C FIGURE 4.3.2.2-37B 143C 144 DATA T4337B/0.0,.1,.2,.3,.4,.5,.6,.8,1.0,2.8, 145 1 0.2,999999./ 146 DATA D4337B/ 147 10.5,.56,.595,.62,.64,.65,.66,.669,.669,.671, 148 20.5,.54,.578,.60,.62,.638,.649,.66,.669,.671/ 149C 150 NX=NXX+.5 151 DCYL=(DN+D1)/2. 152 MACH=FLC(I+2) 153 BETA=SQRT(MACH**2-1.) 154 RLB=XCOOR(NX) 155 DD=2.0*(SPAN-SPANS) 156 TANLE=TANL 157 IF(TANLE.EQ.0.0)TANLE=.00001 158C 159C ***SUPERSONIC H.T.-BODY LIFT CURVE SLOPE,BODY IN PRESENCE OF H.T.*** 160C NON-TRIANGULAR H.T.S 161C 162 IF(TAPEXP.EQ.0.0)GO TO 1050 163 ARG1=BETA*ARSTAR*(1.0+TAPEXP) 164 ARG2=1.+TANLE/BETA 165 TRINO=ARG1*ARG2 166 IF(TRINO.LE.4.)GO TO 1030 167 1000 OLE=AHT(34) 168 DX=RLB-XW-DD*AHT(38)-CRSTAR 169 IF(DX.LE.(-CRSTAR))RKBW=0. 170 IF(DX.LE.(-CRSTAR))GO TO 1020 171 CALL INTKBW(MACH,OLE,CRSTAR,DD,DX,RKBW,RXAC) 172 1020 KBW=RKBW/(RAD*BETA*(SREF/AHT(3))*CLAW*(TAPEXP+1.)*(2.*SPAN/DD-1.)) 173 GO TO 1040 174 1030 CONTINUE 175 LGH(1)=11 176 VAR(1)=DD/(2.*SPAN) 177C 178C FIGURE 4.3.1.2-10 KBW 179C 180 CALL INTERX(1,TFIG10,VAR,LGH,DKBW10,KBW,11,11, 181 1 0,0,0,0,0,0,0,0,Q31210,3,ROUTID) 182 1040 CONTINUE 183 GO TO 1060 184C 185C ***SUPERSONIC H.T.-BODY LIFT CURVE SLOPE,BODY IN PRESENCE OF H.T.*** 186C TRIANGULAR H.T. 187C 188 1050 CONTINUE 189 ARG=BETA*ARSTAR 190 IF(ARG.GT.1.)GO TO 1000 191 GO TO 1030 192 1060 CONTINUE 193C 194C ***SUPERSONIC H.T.-BODY LIFT CURVE SLOPE 195C 196 ALBO=BD(81) 197 IF(BD(81).EQ.UNUSED)ALBO=0.0 198 DO 1070 J=1,NALPHA 199 1070 ALPHAB(J)=FLC(J+22)+ALBO 200 VAR(1)=(SPAN-SPANS)/SPAN 201 LGH(1)=11 202C 203C FIGURE 4.3.1.2-10 KWB 204C 205 CALL INTERX(1,TFIG10,VAR,LGH,DKWB10,KWB,11,11, 206 1 0,0,0,0,0,0,0,0,Q31210,3,ROUTID) 207 CLAWB=CLAW*KWB 208 CLABW=CLAW*KBW 209 CLA=CLABW+CLAWB+CLAB 210C 211C ***SUPERSONIC H.T.-BODY LIFT AT ANGLE OF ATTACK*** 212C 213 IF(TYPE.NE.WTYPE(1))GO TO 1100 214 IF(ALIW.EQ.0.0.OR.(ALIW.EQ.UNUSED))GO TO 1080 215C 216C FIGURE 4.3.1.2-12A (KBW) INCIDENCE 217C 218 CALL INTERX(1,T4312A,VAR,LGH,D4312A,KKBW,11,11, 219 1 0,0,0,0,0,0,0,0,Q1212A,3,ROUTID) 220C 221C FIGURE 4.3.1.2-12B (KWB) INCIDENCE 222C 223 LGH(1)=8 224 CALL INTERX(1,T4312B,VAR,LGH,D4312B,KKWB,8,8, 225 1 0,0,0,0,0,0,0,0,Q1212B,3,ROUTID) 226 1080 CLI=CLAW*ALIW 227 DETCL=0.0 228 XCBO4=XW+AHT(161) 229 BMBS=SPAN-SPANS 230 CALL BODOWG(BD(255),XCBO4,BMBS,SPAN,AHT(27), 231 1 IVBW,GAMMA,NALPHA) 232 ARG1=KWB+KBW 233 ARG2=(KKWB+KKBW)*CLI 234 ARG3=CLAW*(DD/(2.*SPAN)) 235 DO 1090 J=1,NALPHA 236 DETCL=ARG3*ALPHAB(J)*IVBW(J)*GAMMA(J) 237 1090 CL(J)=CLB(J)+ARG1*(CLW(J)-CLI)+ARG2+DETCL 238C 239C ***SUPERSONIC CENTER OF PRESSURE FOR BODY NOSE AND FORBODY*** 240C 241 1100 CONTINUE 242 DELXW=(SPAN-SPANS)*TANLE*COS(ALIW/RAD) 243 RLAP=XW+DELXW-RLN 244 ARG1=RLN 245 IF(RLAP.LT.0.)ARG1=RLN+RLAP 246 IF(RLAP.LT.0.)RLAP=0.0 247 FA=RLAP/DCYL 248 FN=ARG1/DCYL 249 VAR(1)=BETA/FN 250 VAR(2)=FA/FN 251 LGH(1)=6 252 LGH(2)=9 253 IF(BNOSE.EQ.1.)GO TO 1120 254C 255C FIGURE 4.2.2.1-23A LEFT SIDE(XCP) OGIVE 256C 257 IF(VAR(1).GT.1.)GO TO 1110 258 CALL INTERX(2,T4218A,VAR,LGH,DL218A,XCPLN,9,54, 259 1 0,0,0,0,0,0,0,0,Q2118A,3,ROUTID) 260 GO TO 1140 261 1110 VAR(1)=1./VAR(1) 262C 263C FIGURE 4.2.2.1-23A RIGHT SIDE(XCP) OGIVE 264C 265 CALL INTERX(2,T4218A,VAR,LGH,DR218A,XCPLN,9,54, 266 1 0,0,0,0,0,0,0,0,Q2118A,3,ROUTID) 267 GO TO 1140 268 1120 CONTINUE 269C 270C FIGURE 4.2.2.1-23B LEFT SIDE(XCP) CONE 271C 272 LGH(2)=7 273 IF(VAR(1).GT.1.)GO TO 1130 274 CALL INTERX(2,T4218B,VAR,LGH,DL218B,XCPLN,7,42, 275 1 0,0,0,0,0,0,0,0,Q2118B,3,ROUTID) 276 GO TO 1140 277 1130 VAR(1)=1./VAR(1) 278C 279C FIGURE 4.2.2.1-23B RIGHT SIDE(XCP) CONE 280C 281 CALL INTERX(2,T4218B,VAR,LGH,DR218B,XCPLN,7,42, 282 1 0,0,0,0,0,0,0,0,Q2118B,3,ROUTID) 283 1140 CONTINUE 284 XACN=(XCPLN-1.)*(ARG1+RLAP)/CBARR 285C 286C ***SUPERSONIC H.T.-LIFT CARRYOVER ON BODY*** 287C 288 VAR(1)=BETA*DD/CRSTAR 289 VAR(2)=BETA/TANLE 290 ARG=(XW+CR)/RLB 291 IF(ARG.GT.1.)GO TO 1150 292C 293C FIGURE 4.3.2.2-37A(XAC)B(W) 294C 295 LGH(1)=8 296 LGH(2)=3 297 CALL INTERX(2,T4337A,VAR,LGH,D4337A,XACA,8,24, 298 1 0,0,0,0,1,0,0,0,Q2137A,3,ROUTID) 299 GO TO 1160 300C 301C FIGURE 4.3.2.2-37B(XAC)B(W) 302C 303 1150 LGH(1)=10 304 LGH(2)=2 305 CALL INTERX(2,T4337B,VAR,LGH,D4337B,XACA,10,20, 306 1 0,0,0,0,0,0,0,0,Q2137B,3,ROUTID) 307 1160 XACBW=XACA*CRSTAR/CBARR 308 XACW=XACW*CRSTAR/CBARR 309C 310C ***SUPERSONIC H.T.-BODY PITCHING MOMENT SLOPE*** 311C 312 DNUM=XACN*CLAB+XACW*CLAWB+XACBW*CLABW 313 DNOM=CLAB+CLAWB+CLABW 314 XAC =DNUM/DNOM 315 CMA=(DXCG/CBARR-XAC)*CLA 316C 317C ***SUPERSONIC H.T.-BODY DRAG*** 318C 319 CDOWB=CDOW+CDOB 320 DO 1170 J=1,NALPHA 321 COSA=COS(FLC(J+22)/RAD) 322 SINA=SIN(FLC(J+22)/RAD) 323 CD(J)=CDOW+CDL(J)+CDB(J) 324 CN(J)=CL(J)*COSA+CD(J)*SINA 325C 326C COMPUTE DELTA XCP WING-BODY 327C 328 DXCPWB=DXCG/CBARR-SLG(134) 329C 330C COMPUTE DELTA XCP BODY-WING 331C 332 DXCPBW=DXCG/CBARR-SWB(8) 333C 334C COMPUTE DELTA CN VORTEX 335C 336 DCNV=SWB(11+J)*SWB(39+J)*SWB(5)/(2*SPAN)*BD(254+J)*HT(101) 337C 338C COMPUTE WING-BODY PITCHING MOMENT 339C 340 BH(40+J)=BODY(40+J)+HT(60+J)*SWB(35)*DXCPWB 341 1 +HT(101)*SYNA(4)*SWB(2)*DXCPWB/RAD 342 2 +HT(60+J)*SWB(11)*DXCPBW 343 3 +HT(101)*SYNA(4)*SWB(37)*DXCPBW/RAD 344 4 +DCNV*DXCPBW 345 5 +HT(80+J)*(SYNA(3)-SYNA(5))/CBARR 346 1170 CA(J)=CD(J)*COSA-CL(J)*SINA 347 RETURN 348 END 349