1 SUBROUTINE HBTRAN(I) 2C 3C*** COMPUTES (XAC/CBARR)B(H) FOR MACH=1.40,TRANSONIC H.T.-BODY CLALPHA 4C 5 DIMENSION ROUTID(2),Q1211A(3),Q1211B(3),Q31210(3),Q2137A(3), 6 1 Q2137B(3) 7 REAL MACH,NXX,KWB,KBW,IVBW ,KKBW,KKWB 8 LOGICAL FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 9 1 HEAD,PRPOWR,JEQPOW,LOASRT,TVTPAN, 10 2 SUPERS,SUBSON,TRANSN,HYPERS 11 DIMENSION CDW(20),CDB(20),WTYPE(4),IVBW(20),GAMMA(20) 12 DIMENSION T4337A(11),D4337A(24),T4337B(12),D4337B(20) 13 DIMENSION T4311A(38),D4311A(285),DUMYA(150),DUMYB(135) 14 DIMENSION TFIG10(11),DKWB10(11),DKBW10(11) 15 DIMENSION T4311B(24),D4311B(135) 16 DIMENSION LGH(4),VAR(4),CD(20),CN(20),CA(20),CL(20),CLB(20), 17 1 CLW(20),ALPHAB(20) 18 COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 19 1 HEAD,PRPOWR,JEQPOW,LOASRT,TVTPAN, 20 2 SUPERS,SUBSON,TRANSN,HYPERS 21 COMMON /CONSNT/ PI,DEG,UNUSED,RAD 22 COMMON /OVERLY/ NLOG,NMACH,M,NALPHA 23 COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF 24 COMMON /FLGTCD/ FLC(93) 25 COMMON /SYNTSS/ XCG,AA(4),XW 26 COMMON /BODYI/ NXX,XCOOR(20) 27 COMMON /HTI/ WINGIN(154) 28 COMMON /HTDATA/ A(195) 29 COMMON /BDATA/ BD(762) 30 COMMON /SUPBOD/ SBD(227) 31 COMMON /SUPWBB/ SSS(61), SWB(61) 32 COMMON /SBETA/ STB(243), TRA(108) 33 COMMON /IBODY/ PB, BODY(400) 34 COMMON /IHT/ PW, WING(380) 35 COMMON /IBH/ PBW, BW(380) 36 EQUIVALENCE (D4311A(1),DUMYA(1)),(D4311A(151),DUMYB(1)) 37 EQUIVALENCE (CA(1),BW(81)),(CR,WINGIN(6)),(CDW(1),WING(1)), 38 1 (KKWB,SWB(2)),(SPANS,WINGIN(3)),(ARSTAR,A(7)), 39 2 (CD(1),BW(1)),(DD,SWB(5)),(CLB(1),BODY(21)), 40 3 (SPAN,WINGIN(4)),(XACBW,SWB(8)) 41 EQUIVALENCE (CL(1),BW(21)),(KBW,SWB(11)),(TANL,A(62)), 42 1 (IVBW(1),SWB(12)),(RKBW,SWB(32)),(CRSTAR,A(10)), 43 2 (TAPEXP,A(27)),(CN(1),BW(61)),(CDB(1),BODY(1)), 44 3 (CLW(1),WING(21)),(KWB,SWB(35)), 45 4 (DN,SBD(4)),(D1,SBD(5)),(KKBW,SWB(37)) 46 EQUIVALENCE (XACA,SWB(39)),(GAMMA(1),SWB(40)), 47 1 (TRINO,SWB(60)),(ALPHAB(1),BD(255)) 48 EQUIVALENCE (CLAB,BODY(101)),(CLAW,WING(101)),(CLA,BW(101)), 49 1 (CLABW,TRA(72)),(CLAWB,TRA(71)) 50 DATA Q1211A/4H4.3.,4H1.2-,4H11A /,Q2137A/4H4.3.,4H2.2-,4H37A /, 51 1 Q1211B/4H4.3.,4H1.2-,4H11B /,Q2137B/4H4.3.,4H2.2-,4H37B /, 52 2 Q31210/4H4.3.,4H1.2-,4H10 /,ROUTID/4HWBTR,4HAN / 53 DATA WTYPE/4HSTRA,4HDOUB,4HCRAN,4HCURV/ 54C 55C FIGURE 4.3.1,2-10 KWB 56C 57 DATA TFIG10/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/ 58 DATA DKWB10/1.0,1.08,1.16,1.26,1.36,1.46,1.56,1.67,1.78,1.89,2.0/ 59C 60C FIGURE 4.3.1.2-10 KBW 61C 62 DATA DKBW10/0.0,.13,.29,.45,.62,.80,1.0,1.22,1.45,1.70,2.0/ 63C 64C FIGURE 4.3.1.2-11A KBW(WITH AFTERBODY) 65C 66 DATA T4311A/ 67 1.0,.2,.4,.6,.8,1.,1.2,1.4,1.6,1.8,2.,2.4,2.8,3.2,4., 4*0., 68 20.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.,1.2,1.6,2.,3.,4.,8.,10.,999999./ 69 DATA DUMYA/15*0.,.75,.58,.47,.38,.33,.32,.3,.3,.29, 70 1 .28, .27, .25, .23, .21, .19, 71 2 1.45, 1.2, 1., .87, .78, .75, .7, .68, .64, .6, .58, .55, 72 3 .5,.48,.4, 73 4 1.95, 1.7, 1.48, 1.33, 1.2, 1.12, 1.07, 1., .97, .92, 74 5 .88,.82,.75,.72,.64, 75 6 2.5, 2.2, 1.93, 1.74, 1.58, 1.5, 1.42, 1.38, 1.3, 1.25, 76 7 1.2,1.1,1.,.93,.8, 77 8 2.95, 2.6, 2.25, 2.05, 1.89, 1.8, 1.7, 1.65, 1.57, 1.5, 78 9 1.48,1.36,1.27,1.2,1.03, 79 A 3.3, 2.92, 2.6, 2.35, 2.18, 2.04, 1.95, 1.86, 1.8, 80 B 1.73,1.7,1.59,1.49,1.4,1.21, 81 C 3.7, 3.32, 2.99, 2.7, 2.5, 2.33, 2.23, 2.12, 2.03, 82 D 1.97,1.9,1.8,1.7,1.6,1.42, 83 E 4.1, 3.64, 3.25, 2.97, 2.77, 2.6, 2.45, 2.32, 2.23, 84 F 2.15,2.1,1.95,1.85,1.76,1.62, 85 G 4.3, 3.83, 3.49, 3.2, 2.99, 2.8, 2.65, 2.5, 2.42, 2.33, 86 H 2.29,2.14,2.,1.89,1.74/ 87 DATA DUMYB/4.6,4.08,3.69,3.35,3.14,2.95,2.8, 88 1 2.69, 2.6, 2.49, 2.4, 2.25, 2.16, 2.02, 1.85, 89 2 5., 4.37, 3.98, 3.64, 3.39, 3.19, 3.03, 2.89, 2.8, 90 3 2.7, 2.61, 2.43, 2.33, 2.21, 2., 91 4 5.45, 4.75, 4.29, 3.95, 3.7, 3.5, 3.34, 3.2, 3.1, 3., 92 5 2.91,2.75,2.6,2.48,2.28, 93 6 5.9, 5.1, 4.63, 4.25, 4., 3.8, 3.65, 3.5, 3.33, 3.23, 94 7 3.14,2.95,2.8,2.67,2.4, 95 8 6.2, 5.5, 5., 4.65, 4.35, 4.15, 3.98, 3.8, 3.67, 96 9 3.53,3.4,3.24,3.1,2.92,2.68, 97 A 6.6, 5.75, 5.25, 4.9, 4.6, 4.4, 4.22, 4.05, 3.91, 98 B 3.78,3.68,3.45,3.26,3.1,2.83, 99 C 7.05, 6.15, 5.6, 5.2, 4.95, 4.7, 4.49, 4.3, 4.13, 100 D 4.,3.85,3.68,3.5,3.34,3.1, 101 E 7.4, 6.35, 5.8, 5.4, 5.12, 4.85, 4.67, 4.49, 4.33, 102 F 4.2,4.07,3.85,3.65,3.5,3.24, 103 G 8., 6.73, 6.2, 5.8, 5.5, 5.25, 5., 4.82, 4.67, 104 H 4.5,4.35,4.12,3.9,3.8,3.59/ 105C 106C FIGURE 4.3.1.2-11B KBW(NO AFTERBODY) 107C 108 DATA T4311B/0.,.2,.4,.6,.8,1.,1.2,1.4,1.6,1.8,2.,2.4,2.8,3.2,4., 109 1 0.,.2,.4,.6,.8,1.,2.,4.,999999./ 110 DATA D4311B/15*0.,1.2,1.,.75,.58,.45,.37,.3,.28,.24,.2,.18, 111 1 .15,.13,.11,.1, 112 2 2.4, 1.85, 1.42, 1.1, .92, .77, .66, .6, .51, .48, .41, 113 3 .31,.26,.22,.22, 114 4 3.5, 2.6, 2., 1.58, 1.28, 1.07, .91, .78, .7, .6, .52, 115 5 .42,.39,.33,.3, 116 6 4.3, 3.1, 2.35, 1.8, 1.45, 1.2, 1.04, .92, .83, .75, 117 7 .68,.55,.48,.4,.38, 118 8 5., 3.65, 2.78, 2.22, 1.82, 1.5, 1.3, 1.15, 1., .92, 119 9 .82,.69,.58,.5,.42, 120 A 5.75, 4.55, 3.66, 2.9, 2.3, 1.92, 1.62, 1.42, 1.25, 1.1, 121 B .97,.82,.7,.63,.52, 122 C 6.7, 5.25, 4.18, 3.32, 2.63, 2.2, 1.88, 1.62, 1.45, 123 D 1.28,1.15,1.,.85,.75,.62, 124 E 7.6, 6.2, 4.91, 3.95, 3.18, 2.6, 2.22, 1.91, 1.7, 1.5, 125 F 1.35,1.12,.97,.85,.7/ 126C 127C FIGURE 4.3.2.2-37A 128C 129 DATA T4337A/0.0,.4,.8,1.2,1.6,2.,2.4,2.8, 130 1 .1,1.0,999999./ 131 DATA D4337A/ 132 1.5,.72,.900,1.08,1.24,1.39,1.53,1.68, 133 2.5,.72,.910,1.09,1.25,1.41,1.57,1.72, 134 3.5,.73,.920,1.11,1.27,1.43,1.59,1.74/ 135C 136C FIGURE 4.3.2.2-37B 137C 138 DATA T4337B/0.0,.1,.2,.3,.4,.5,.6,.8,1.0,2.8, 139 1 0.2,999999./ 140 DATA D4337B/ 141 10.5,.56,.595,.62,.64,.65,.66,.669,.669,.671, 142 20.5,.54,.578,.60,.62,.638,.649,.66,.669,.671/ 143C 144 NX=NXX+.5 145 DCYL=(DN+D1)/2. 146 MACH=FLC(I+2) 147 RLB=XCOOR(NX) 148 DD=2.0*(SPAN-SPANS) 149 TANLE=TANL 150 IF(TANLE.EQ.0.0)TANLE=.00001 151 IF(MACH.EQ.1.)BETA=.0000001 152 IF(MACH .EQ. 1.) GO TO 1030 153 IF(MACH.GT.1.)BETA=SQRT(MACH**2-1.) 154 IF(MACH.LT.1.)BETA=SQRT(1.-MACH**2) 155C 156C ***SUPERSONIC WING-BODY LIFT CURVE SLOPE,BODY IN PRESENCE OF WING*** 157C NON-TRIANGULAR WINGS 158C 159 IF(TAPEXP.EQ.0.0)GO TO 1050 160 ARG1=BETA*ARSTAR*(1.0+TAPEXP) 161 ARG2=1.+TANLE/BETA 162 TRINO=ARG1*ARG2 163 IF(TRINO.LE.4.)GO TO 1030 164 1000 LGH(1)=15 165 ARG=(XW+CR)/RLB 166 VAR(1)=BETA* DD/CRSTAR 167 VAR(2)=BETA/TANLE 168 IF(ARG.GT.1.)GO TO 1010 169 LGH(2)=19 170C 171C FIGURE 4.3.1.2-11A 172C 173 CALL INTERX(2,T4311A,VAR,LGH,D4311A,RKBW,19,285, 174 1 2,0,0,0,1,0,0,0,Q1211A,3,ROUTID) 175 GO TO 1020 176 1010 LGH(2)=9 177C 178C FIGURE 4.3.1.2-11B 179C 180 CALL INTERX(2,T4311B,VAR,LGH,D4311B,RKBW,15,135, 181 1 2,0,0,0,1,0,0,0,Q1211B,3,ROUTID) 182 1020 KBW=RKBW/(RAD*BETA*(SREF/A(3))*CLAW*(TAPEXP+1.)*(2.*SPAN/DD-1.)) 183 GO TO 1040 184 1030 CONTINUE 185 LGH(1)=11 186 VAR(1)=DD/(2.*SPAN) 187C 188C FIGURE 4.3.1.2-10 KBW 189C 190 CALL INTERX(1,TFIG10,VAR,LGH,DKBW10,KBW,11,11, 191 1 0,0,0,0,0,0,0,0,Q31210,3,ROUTID) 192 1040 CONTINUE 193 GO TO 1060 194C 195C ***SUPERSONIC WING-BODY LIFT CURVE SLOPE,BODY IN PRESENCE OF WING*** 196C TRIANGULAR WING 197 1050 CONTINUE 198 ARG=BETA*ARSTAR 199 IF(ARG.GT.1.)GO TO 1000 200 GO TO 1030 201 1060 CONTINUE 202C 203C ***SUPERSONIC WING-BODY LIFT CURVE SLOPE 204C 205 ALBO=BD(81) 206 IF(BD(81).EQ.UNUSED)ALBO=0.0 207 DO 1070 J=1,NALPHA 208 1070 ALPHAB(J)=FLC(J+22)+ALBO 209 VAR(1)=(SPAN-SPANS)/SPAN 210 LGH(1)=11 211C 212C FIGURE 4.3.1.2-10 KWB 213C 214 CALL INTERX(1,TFIG10,VAR,LGH,DKWB10,KWB,11,11, 215 1 0,0,0,0,0,0,0,0,Q31210,3,ROUTID) 216 CLAWB=CLAW*KWB 217 CLABW=CLAW*KBW 218 CLA=CLABW+CLAWB+CLAB 219C 220C ***SUPERSONIC WING-LIFT CARRYOVER ON BODY*** 221C 222 BETA=.98 223 VAR(1)=BETA*DD/CRSTAR 224 VAR(2)=BETA/TANLE 225 ARG=(XW+CR)/RLB 226 IF(ARG.GT.1.)GO TO 1080 227C 228C FIGURE 4.3.2.2-37A(XAC)B(W) 229C 230 LGH(1)=8 231 LGH(2)=3 232 CALL INTERX(2,T4337A,VAR,LGH,D4337A,XACA,8,24, 233 1 0,0,0,0,1,0,0,0,Q2137A,3,ROUTID) 234 GO TO 1090 235C 236C FIGURE 4.3.2.2-37B(XAC)B(W) 237C 238 1080 LGH(1)=10 239 LGH(2)=2 240 CALL INTERX(2,T4337B,VAR,LGH,D4337B,XACA,10,20, 241 1 0,0,0,0,0,0,0,0,Q2137B,3,ROUTID) 242 1090 XACBW=XACA*CRSTAR/CBARR 243 RETURN 244 END 245