1 SUBROUTINE TRSONJ(NZ) 2C 3C*** COMPUTES TRANSONIC H.T. LIFT SLOPE,CLMAX,ALPHA CLMAX,AND CD0 4C*** BODY LIFT AND MOMENT SLOPES,DRAG AT ANGLE OF ATTACK 5C 6 DIMENSION ROUTID(2) 7 DIMENSION Q3243A(3),Q3243B(3),Q3244B(3),Q3244A(3),Q3244C(3),Q34171 8 1 (4),Q34172(4),Q3418A(3),Q34181(3),Q34182(3),Q3419B(3),Q3419C(3), 9 2 Q3419A(3),Q15127(4),Q5129R(3),Q5129L(3),Q23124(3),Q23126(3) 10 REAL K,MACH,MFB,MT(5) 11 REAL MFB0 12 DIMENSION T429L(15),D429L(56),T429R(14),D429R(35),T424(16), 13 1 D424(39) 14 DIMENSION WTYPE(4),C17(6) ,TR(11),DR(11) 15 DIMENSION T426(16),D426(55) 16 DIMENSION X(11),Y(11),C(6) 17 DIMENSION T418A(10),D418A(10), T18B1(3),D18B1(20),C18B1(60), 18 1 T18B2(10), D18B2(10) ,C18B2(100) 19 DIMENSION T419A(22),D419A(77) ,T419B(13),D419B(13),T419C(26), 20 1 D419C(130) 21 DIMENSION T43A(15),D43A(56),T43B(13),D43B(40),T44A(13),D44A(42), 22 1 T44B(13),D44B(42),T44C(6),D44C(6) 23 DIMENSION XMTD(15),CDW1(15),CDW2(15),VAR(4),LGH(4) 24 DIMENSION CLAMT(5),ALPHA(20) ,CDJB(20) 25 DIMENSION X27M(4),X27I(4),C27(6) 26 COMMON /CONSNT/ PI,DEG,UNUSED,RAD 27 COMMON /OVERLY/ NLOG,NMACH,M,NALPHA,IG,NF 28 COMMON /OPTION/ SR,CBARR,RUFF,BLREF 29 COMMON /FLGTCD/ FLC(93) 30 COMMON /BDATA/ BD(762) 31 COMMON /HTDATA/ A(195),B(49) 32 COMMON /HTI/ WINGIN(77) 33 COMMON /SUPBOD/ SBD(227) 34 COMMON /SBETA/ STB(243),TRA(108) 35 COMMON /IBODY/ PB, BODY(400) 36 COMMON /IHT/ PW, WING(380) 37 COMMON /IBH/ PBW, BW(380) 38 COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 39 1 HEAD,PRPOWR,JEQPOW,LOASRT,TVTPAN, 40 2 SUPERS,SUBSON,TRANSN,HYPERS 41 LOGICAL FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 42 1 HEAD,PRPOWR,JEQPOW,LOASRT,TVTPAN, 43 2 SUPERS,SUBSON,TRANSN,HYPERS 44 EQUIVALENCE (TOC,WINGIN(16)),(ARSTAR,A(7)),(SWEPE2,A(71)) 45 1,(COSLE,A(61)),(TAPR,A(27)),(CBAR,A(16)),(RNFS,A(129)), 46 2(XOC,WINGIN(18)),(SRSTAR,A(3)),(COSC4,A(67)),(SS,BD(93)) 47 3,(DB,SBD(6)),(DMAX,SBD(120)),(XLB,SBD(2)),(SB,BD(57)) 48 4,(TANC2,A(74)) 49 EQUIVALENCE (CLA14,TRA(1)),(CLAMO,WINGIN(69)),(K,TRA(3)) 50 1,(MFB0,TRA(5)),(MFB,TRA(6)),(AOC,TRA(7)),(CFBCT,TRA(8)),(BETAFB, 51 2TRA(9)),(CLAFBT,TRA(10)),(MACH,TRA(4)),(CLAFB,TRA(12)),(CLAA,TRA 52 3(13)),(BOC,TRA(14)),(CLAB,TRA(15)),(MT(1),TRA(16)),(CLAMT(1),TRA 53 4(21)),(CLA,WING(101)),(C1,TRA(27)),(ARATIO,TRA(28)),(BU4,TRA(29)), 54 5(CLMAX6,TRA(30)),(ACLBAS,TRA(31)),(DACMA6,TRA(32)),(C3,TRA(33)) 55 EQUIVALENCE (CLAB14,SBD(18)),(CMAB14,SBD(110)),(CD014,SBD(124)), 56 1(CMAB6,BODY(121)),(CFB6,BD(92)),(CD0B6,BD(61)),(CDBB6,BD(60)) 57 EQUIVALENCE (DALCM,TRA(34)),(DCLMAX,TRA(35)),(ALCLM6,TRA(36)), 58 1(ALCLMT,TRA(37)),(CLMAXT,TRA(38)),(RLCOFF,TRA(39)),(RNN,TRA(40)), 59 2(RL,TRA(41)),(CF,TRA(42)),(CDW2(1),TRA(43)), 60 3(CDW,TRA(67)),(CDF,TRA(68)),(CD0W,WING(1)),(CLAW6,TRA(70)) 61 4,(CLABM6,BODY(101)),(CLABD,BODY(101)),(CMAB,BODY(121)) 62 5,(CDPB,TRA(80)),(CDBFIG,TRA(81)),(CDBB,TRA(76)),(CDWB,TRA(77)) 63 6,(CD0B,TRA(78)),(CDJB(1),BODY(1)),(CD0WB,BW(1)) 64 7,(CDFB,TRA(79)),(DCLA14,TRA(82)) 65 EQUIVALENCE (TANLE,A(62)) 66 LOGICAL LOWAR 67 DATA ROUTID /4HTRSO,4HNI / 68 DATA Q3243A/4H4.1.,4H3.2-,4H53A /,Q3243B/4H4.1.,4H3.2-,4H53B /, 69 1Q3244B/4H4.1.,4H3.2-,4H54B /,Q3244A/4H4.1.,4H3.2-,4H54A /,Q3244C/ 70 24H4.1.,4H3.2-,4H54C /,QCLA/4HCLA /,Q34171/4H4.1.,4H3.4-,4H24B , 4H 71 3C1 /,Q34172/4H4.1.,4H3.4-,4H24B ,4HC2 /,Q3418A/4H4.1.,4H3.4-, 4H 72 425A /,Q34181/4H4.1.,4H3.4-,4H25B1/,Q34182/4H4.1.,4H3.4-,4H25B2/, 73 5Q3419B/4H4.1.,4H3.4-,4H26B /,Q3419C/4H4.1.,4H3.4-,4H26C /,Q3419A 74 6/4H4.1.,4H3.4-,4H26A /,Q15127/4H4.1.,4H5.1-,4H27 ,4HCEPT/,Q5129R 75 7/4H4.1.,4H5.1-,4H29R /,Q5129L/4H4.1.,4H5.1-,4H29L /,QCDW/4HCDW /, 76 8Q23124/4H4.2.,4H3.1-,4H24 /,Q23126/4H4.2.,4H3.1-,4H26 / 77C 78C FIGURE 4.1.5.1-27 79C 80 DATA X27I /1.57780,1.67221,1.98509,2.28874/ 81 DATA X27M /0.0,1.0,2.0,3.0/ 82 DATA I27 /0/ 83 DATA WTYPE /4HSTRA,4HDOUB,4HCRAN,4HCURV/ 84 DATA I17,IN/0,0/ 85 DATA NP/11/ 86C 87C FIGURE 4.1.3.4-24B 88C 89 DATA TR /0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1./ 90 DATA DR /0.,.21,.5,.9,1.08,1.05,1.,.94,.9,.86,.85/ 91 DATA Y /0.,.225,.47,.496,.43,.32,.21,.125,.075,.0475,.0/ 92 DATA X /0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1./ 93C 94C ----4.1.3.4-25(A) 95C 96 DATA T418A /0.,.4,.8,1.2,1.6,2.0,2.4,2.8,3.2,3.6/ 97 DATA D418A /3*35.,32.,28.,25.,23.2,22.,21.5,21./ 98C 99C ----4.1.3.4-25(B) 100C 101 DATA T18B1/.2,.4,.6/ 102 DATA D18B1/4.5,5.,5.5,6.,6.5,7.,7.5,8.,8.5,9.,9.5,10.,10.5,11., 103 111.5,12.,12.5,13.,13.5,14./ 104 DATA C18B1/0.,.5,.9,1.4,1.9,2.5,3.3,4.,4.6,5.6,6.4,7.3,8.2,9.2,10. 105 1,11.,12.,13.,14.,15. , 0.,.2,.4,.7,1.2,1.7,2.4,3.,3.7,4.6,5.2,6., 106 26.9,7.8,8.6,9.5,10.4,11.4,12.3,13.5 , 0.,0.,.1,.2,.5,.7,1.,1.3,1.6 107 3,2.,2.5,3.,3.6,4.3,4.9,5.5,6.2,7.,7.6,8.5/ 108 DATA T18B2/0.,2.,3.,4.,5.,6.,7.,8.,9.,30./ 109 DATA C18B2 / 10.,8.5,6.9,5.5,4. 110 1,2.6,1.5,.7,.1,0. , 8.7,7.3,5.3,4.2,2.6,1.4,0.5,-.2,-.5,0. , 7.5, 111 25.9,4.2,2.5,1.2,0.,-.7,-1.1,-.8,0. , 5.5,3.4,1.6,0.,-1.3,-2.1,-2.5 112 3,-2.,-.8,0. , 3.,.7,-1.4,-3.3,-4.3,-4.3,-3.1,-2.,-.8,0. , .3,-2.5, 113 4-4.7,-5.8,-5.3,-4.3,-3.1,-2.,-.8,0. ,-2.2,-5.,-6.7,-6.3,-5.3,-4.3, 114 5-3.1,-2.,-0.8,0. , -3.3,-6.6,-7.2,-6.3,-5.3,-4.3,-3.1,-2.,-.8,0., 115 6 -4.2,-7.,-7.2,-6.3,-5.3,-4.3,-3.1,-2.,-.8,0. , -8.5,-7.9,-7.2, 116 7 -6.3,-5.3,-4.3,-3.1,-2.,-.8,0. / 117 DATA D18B2 /.0,.5,1.,1.5,2.,2.5,3.,3.5,4.,4.5/ 118C 119C FIGURE 4.1.3.4-26B 120C 121 DATA T419B 122 1 / 0., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13., 14., 16. / 123 DATA D419B 124 1 / 1., 1., .93, .8, .55, .4, .3, .24, .2, .16, .14, .12, .12 / 125C 126C FIGURE 4.1.3.2-53A 127C 128 DATA T43A/ 129 1 0., 2., 4., 6., 8., 10., 12., 14. , 130 2 0., 1., 2., 3., 4., 6., 8./ 131 DATA D43A 132 1/8*1.0, 133 2 6*1.0,.98,.9, 134 3 5*1.0,.983,.89,.855, 135 4 4*1.0,.939,.888,.857,.835, 136 5 2*1.0,.958,.904,.869,.842,.82,.802, 137 6 1.0, .952, .908, .872, .842, .816, .80, .787, 138 7 1.0, .952, .908, .872, .842, .816, .80, .787 / 139C 140C FIGURE 4.1.3.2-53B 141C 142 DATA T43B 143 1 / 0.0, 0.17453, 0.34907 , 0.52360, 0.69813, 144 2 1.04720, 1.22173 , 1.57080 , 145 3 .799,.85,.9,.95,1.000/ 146 DATA D43B/ 147 1 .8, .807, .82, .841, .867, .918, .945, 1.0, 148 2 .85, .855, .865, .88, .9, .938, .958, 1.0, 149 3 .90, .903, .91, .92, .932, .958, .974, 1.0, 150 4 .95, .952, .955, .96, .967, .98, .987, 1.0, 151 5 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,1.0, 1.0 / 152C 153C FIGURE 4.1.3.2-54A 154C 155 DATA T44A/ 156 1 0., 2., 4., 6., 8., 12., 14. , 157 2 1., 2., 3., 4., 6., 8. / 158 DATA D44A/ 159 1 1.09, 1.088, 1.056, 1.008, .962, .872, .828, 160 2 .98, 1.088, 1.1, 1.05, .992, .882, .83, 161 3 .85, 1., 1.075, 1.06, 1.005, .89, .835, 162 4 .74, .9, 1.015, 1.07, 1.025, .895, .837, 163 54*1.15,1.07,.915,.915, 164 65*1.12,.94,.94/ 165C 166C FIGURE 4.1.3.2-54B 167C 168 DATA T44B/ 169 1 0., 4., 6., 8., 10., 12., 14. , 170 2 1., 2., 3., 4., 6., 8. / 171 DATA D44B/ 172 1 -.14, -.14, -.14, -.14, -.14, 0., .44, 173 2 -.07, -.07, -.07, -.07, 0., .31, .64, 174 3 .05, .05, .05, .06, .23, .49, .75, 175 4 .06, .06, .06, .15, .35, .58, .79, 176 5 .08, .08, .14, .29, .48, .66, .84, 177 6 .09, .15, .24, .39, .555, .72, .89 / 178C 179C FIGURE 4.1.3.2-54C 180C 181 DATA T44C 182 1/0., 2., 4., 6., 8., 16. / 183 DATA D44C 184 1 / -.04, .01, .075, .13, .15, .15 / 185C 186C FIGURE 4.1.3.4-26A 187C 188 DATA T419A 189 1 / 0.8, 1.0, 1.2, 1.6, 2.0, 2.4, 3.0, 4*0., 190 2 .6, .7, .8, .9, .95, 1.0, 1.05, 1.1, 1.15, 1.2, 1.4 / 191 DATA D419A 192 1 / 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 193 2 -.050,-.050,-.050,-.030,-.015,-.010, 0.0 , 194 3 -.115,-.100,-.083,-.047,-.018, 0.0 , .025, 195 4 -.145,-.110,-.080,-.020, .022, .050, .097, 196 5 -.120,-.080,-.030, .050, .100, .140, .204, 197 6 -.060, 0.0 , .060, .150, .210, .250, .325, 198 7 .090, .130, .180, .257, .310, .330, .360, 199 8 .113, .165, .210, .280, .320, .345, .360, 200 9 .108, .160, .205, .284, .320, .344, .360, 201 A .096, .145, .190, .273, .320, .343, .355, 202 B .020, .050, .100, .200, .300, .330, .340/ 203C 204C FIGURE 4.1.3.4-26C 205C 206 DATA T419C 207 1 / 0.0, 0.4, 0.8, 1.0, 1.2, 1.6, 2.0, 2.4, 3.0, 5.0 ,3*0., 208 2 .6,.7,.75,.8,.85,.9,.95,1.0,1.05,1.1,1.15,1.2,1.4 / 209 DATA D419C 210 1 / 0.0, 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0, 0.0 , 0.0 , 211 2 -1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00, 212 3 -1.25,-1.22,-1.20,-1.20,-1.20,-1.10,-1.10,-1.00,-1.00,-0.85, 213 4 -2.05,-1.75,-1.50,-1.40,-1.30,-1.20,-1.10,-1.00,-1.00,-0.95, 214 5 -2.20,-1.90,-1.60,-1.50,-1.40,-1.20,-1.10,-1.00,-0.90,-0.63, 215 6 -1.85,-1.60,-1.40,-1.30,-1.20,-1.00,-0.80,-0.70,-0.60,-0.40, 216 7 -1.80,-1.38,-1.00,-0.90,-0.70,-0.50,-0.30,-0.10, 0.0 , 0.10, 217 8 -0.60,-0.20, 0.20, 0.30, 0.50, 0.80, 1.00, 1.20, 1.40, 2.00, 218 9 1.90, 2.40, 3.00, 3.20, 3.50, 3.90, 4.30, 4.60, 5.00, 6.20, 219 A 2.40, 4.05, 5.50, 6.20, 7.00, 8.00, 8.90, 10.0, 10.9, 13.9, 220 B 1.60, 4.04, 6.90, 8.00, 9.10, 10.6, 11.9, 13.0, 14.3, 18.2, 221 C -1.20, 3.90, 7.60, 8.90, 10.2, 12.1, 13.5, 14.7, 15.9, 18.8, 222 D -11.0,-0.60, 8.10, 10.5, 12.0, 14.4, 16.2, 17.5, 18.8, 21.3 / 223C 224C FIGURE 4.1.5.1-29 LEFT SIDE (SUBSONIC) 225C 226 DATA T429L/ 227 10.0,0.4,0.6,0.8,1.0,1.2,1.4,1.475, 228 20.,.5,1.0,1.5,2.0,3.0,4.0/ 229 DATA D429L/ 230 A8*0.0, 231 10.90,0.71,0.55,0.39,0.21,0.10,0.00,0.0, 232 22.34,2.03,1.70,1.05,0.51,0.21,0.01,0.0, 233 33.01,2.92,2.68,1.99,1.00,0.32,0.01,0.0, 234 43.30,3.27,3.12,2.50,1.26,0.42,0.03,0.0, 235 53.49,3.43,3.38,2.80,1.52,0.51,0.06,0.0, 236 63.61,3.58,3.52,3.14,1.83,0.52,0.09,0.0/ 237C 238C FIGURE 4.1.5.1-29 RIGHT SIDE (SUPERSONIC) 239C 240 DATA T429R/ 241 10.0,.4,.8,1.2,1.4,0.0,0.0, 242 20.,.5,1.0,1.5,2.0,3.0,4.0/ 243 DATA D429R/ 244 A5*0.0, 245 10.90,1.08,1.18,1.28,1.31, 246 22.34,2.53,2.68,2.75,2.78, 247 33.01,3.08,3.10,3.12,3.12, 248 43.30,3.32,3.32,3.30,3.28, 249 53.49,3.49,3.48,3.42,3.39, 250 63.61,3.61,3.60,3.56,3.51/ 251C 252C FIGURE 4.2.3.1-24 253C 254 DATA T424/ 255 1.75,.85,.90,.95,1.0,1.015,1.030,1.050,1.075,1.125,1.15,1.2,1.3, 256 2.011,.067,.122/ 257 DATA D424/ 258 1.011,.014,.018,.031,.069,.040,.020,6*0.0, 259 2.067,.071,.078,.092,.140,.150,.120,.086, .067, .050,.047,.045,.05, 260 3.122,.130,.139,.155,.204,.220,.236,.210, .168, .110, .101,.100,.1/ 261 DATA T426/ 262 16.0,7.0,8.0,10.,12.,14.,16.,18.,20.,24.,26., 263 21.0,1.025,1.05,1.1,1.2/ 264 DATA D426/ 265 1.118,.084,.065,.041,.029,.020,.016,.012,.010,.007,.006, 266 2.147,.110,.088,.061,.043,.033,.026,.022,.019,.015,.013, 267 3.166,.120,.097,.068,.050,.039,.031,.026,.022,.015,.013, 268 4.186,.138,.110,.076,.057,.044,.034,.028,.022,.015,.013, 269 5.200,.152,.120,.080,.059,.046,.036,.028,.022,.015,.013/ 270C 271C *** WING ALONE *** 272C *** WING LIFT CURVE SLOPE *** 273C 274 MACH=FLC(NZ+2) 275 IF(WINGIN(15).NE.WTYPE(1))GO TO 1100 276 CALL TRNHT(CLA14, DCLA14) 277 K=CLAMO*RAD/(2.0*PI) 278 BETA6=0.80 279C 280C FIGURE 4.1.3.2-53A FORCE BREAK MACH NUMBER=MFB ( 0 SWEEP) 281C 282 VAR(1)=TOC*100. 283 VAR(2)=ARSTAR 284 LGH(1)=8 285 LGH(2)=7 286 CALL INTERX(2,T43A,VAR,LGH,D43A,MFB0,8,56, 287 1 2,2,0,0,2,-1,0,0,Q3243A,3,ROUTID) 288C 289C FIGURE 4.1.3.2-53B MFB0 CORRECTED FOR SWEEP 290C 291 VAR(1)=SWEPE2 292 VAR(2)=MFB0 293 LGH(1)=8 294 LGH(2)=5 295 CALL INTERX(2,T43B,VAR,LGH,D43B,MFB,8,40, 296 1 0,2,0,0,2,0,0,0,Q3243B,3,ROUTID) 297 IF(NF .LT. 0) RETURN 298C 299C FIGURE 4.1.3.2-54B (A/C) 300C 301 VAR(1)=TOC*100. 302 VAR(2)=ARSTAR 303 LGH(1)=7 304 LGH(2)=6 305 CALL INTERX(2,T44B,VAR,LGH,D44B,AOC,7,42, 306 1 0,2,0,0,1,2,0,0,Q3244B,3,ROUTID) 307C 308C FIGURE 4.1.3.2-54A CORRECTION TO LIFT CURVE SLOPE 309C 310 CALL INTERX(2,T44A,VAR,LGH,D44A,CFBCT,7,42, 311 1 2,1,0,0,2,2,0,0,Q3244A,3,ROUTID) 312C 313C 314C 315 BETAFB = SQRT(1.0-MFB**2) 316 IF(MFB .GT. 0.98) BETAFB = 0.0 317 BB = BETAFB 318 ARG1 = 2.0*PI*ARSTAR/RAD 319 ARG2=(ARSTAR/K)**2 320 ARG3=BB**2+TANC2**2 321 ARG4 = 2.0+SQRT(ARG2*ARG3+4.0) 322 CLAFBT = ARG1/ARG4 323C 324 BB = 0.8 325 ARG2 = (BB*ARSTAR/K)**2 326 ARG3 = 0.0 327 IF(ARG2 .GT. 0.0) ARG3 = 1.0+(TANC2/BB)**2 328 ARG4 = 2.0+SQRT(ARG2*ARG3+4.0) 329 CLAW6 = ARG1*SRSTAR/(ARG4*SR) 330C 331 XM = 0.75 332 IF(XM .GT. (MFB-0.1)) XM = MFB-0.1 333 IF(XM .GT. MACH) XM=MACH 334 BB = SQRT(1.0-XM**2) 335 ARG2=(ARSTAR/K)**2 336 ARG3=BB**2+TANC2**2 337 ARG4 = 2.0+SQRT(ARG2*ARG3+4.0) 338 CLAW7 = ARG1*SRSTAR/(ARG4*SR) 339 DCLA7 = XM*ARG1*ARG2*SRSTAR/SR 340 DCLA7 = DCLA7/(ARG4**2*SQRT(ARG2*ARG3+4.0)) 341 CLAFB=CLAFBT*CFBCT *SRSTAR/SR 342 CLAA=(1.-AOC)*CLAFB 343C 344C FIGURE 4.1.3.2-54C (B/C) 345C 346 VAR(1)=TOC*100. 347 LGH(1)=6 348 CALL INTERX(1,T44C,VAR,LGH,D44C,BOC,6,6, 349 1 2,0,0,0,1,0,0,0,Q3244C,3,ROUTID) 350 CLAB=(1.-BOC)*CLAFB 351 MT(1)=XM 352 MT(2)=MFB 353 MT(3)=MFB+.07 354 MT(4)=MFB+.14 355 MT(5)=1.40 356 CLAMT(1)=CLAW7 357 CLAMT(2)=CLAFB 358 CLAMT(3)=CLAA 359 CLAMT(4)=CLAB 360 CLAMT(5)=CLA14 361 CALL TRANF(5, MT, CLAMT, DCLA7, DCLA14, MACH, CLA) 362C 363C *** CLASSIFICATION OF ASPECT RATIO *** 364C 365 CALL TBFUNX(A(27), C1 ,DYDX,NP,X,Y,C,IN,MI,NG, 366 1 0,0,Q34171,4,ROUTID) 367 ARG1=(C1+1.)*COSLE 368 ARATIO=A(128)/ARG1 369C 370C FIGURE 4.1.3.4-24B (C1) AND (C2) 371C 372 VAR(1)=TAPR 373 LGH(1)=11 374 CALL INTERX(1,TR,VAR,LGH,DR,C2,11,11, 375 1 0,0,0,0,0,0,0,0,Q34172,4,ROUTID) 376 LOWAR=.FALSE. 377 IF(ARSTAR.LT.ARATIO)LOWAR=.TRUE. 378C 379C *** IF LOW ASPECT RATIO CALCULATE CLMAX AND (ALPHA)CLMAX *** 380C 381 IF(.NOT.LOWAR)GO TO 1070 382 BU4=(C1+1.)*ARSTAR*COSLE/BETA6 383 ARG=(C2+1.)*ARSTAR*TANLE 384 A(160)=ARG 385 CALL CLMXB1(BU4,CLMAX6,A,WINGIN) 386C 387C *** CALCULATE (ALPHA)CLMAX AT MACH=0.60 *** 388C FIGURE 4.1.3.4-25A 389C 390 CALL TBFUNX(BU4,ACLBAS,DYDX,10,T418A,D418A,C17,I17,MI,NG, 391 1 0,0,Q3418A,3,ROUTID) 392 IF(ARG.LE.4.5)GO TO 1050 393C 394C FIGURE 4.1.3.4-25B 395C 396 CALL TLINEX(T18B1,D18B1,C18B1,3,20,MT(1),ARG,DACMA6, 397 1 0,0,2,1,Q34181,3,ROUTID) 398 GO TO 1060 399 1050 ARG1=ARSTAR*COSLE*(1.+4.*TAPR**2) 400 CALL TLINEX(T18B2,D18B2,C18B2,10,10,ARG1,ARG,DACMA6, 401 1 0,2,2,0,Q34182,3,ROUTID) 402 1060 ALCLM6=ACLBAS+DACMA6 403C 404C *** ADJUST PARAMETERS TO TRANSONIC MACH NUMBER *** 405C FIGURE 4.1.3.4-26B 406C 407 VAR(1)=ARG 408 LGH(1)=13 409 CALL INTERX(1,T419B,VAR,LGH,D419B,C3,13,13, 410 1 0,0,0,0,0,0,0,0,Q3419B,3,ROUTID) 411C 412C FIGURE 4.1.3.4-26C DELTA (ALPHA)CLMAX 413C 414 VAR(1)=BU4*BETA6 415 VAR(2)=MACH 416 LGH(1)=10 417 LGH(2)=13 418 CALL INTERX(2,T419C,VAR,LGH,D419C,DALCM ,13,130, 419 1 2,0,0,0,2,2,0,0,Q3419C,3,ROUTID) 420C 421C FIGURE 4.1.3.4-26A DELTA CLMAX 422C 423 VAR(1)=C3*VAR(1) 424 LGH(1)=7 425 LGH(2)=11 426 CALL INTERX(2,T419A,VAR,LGH,D419A,DCLMAX,11,77, 427 1 2,0,0,0,2,2,0,0,Q3419A,3,ROUTID) 428 ALCLMT=ALCLM6+DALCM 429 CLMAXT=CLMAX6+DCLMAX*(SRSTAR/SR) 430 1070 CONTINUE 431C 432C *** ZERO ANGLE OF ATTACK DRAG,CD0 *** 433C FIGURE 4.1.5.1-27 (EQUATION FOR RLCOFF) 434C 435 ARG=12.*CBAR/RUFF 436 CALL TBFUNX(0.6,CEPT,DYDX,4,X27M,X27I,C27,I27,MI,NG, 437 1 0,0,Q15127,4,ROUTID) 438 RLCOFF=ARG**1.0482*10.0**CEPT 439 RNN=CBAR*RNFS 440 IF(RLCOFF.LT.RNN)RNN=RLCOFF 441 CALL FIG26(RNN,0.60,CF) 442 RL=1.2 443 IF(XOC.LT.0.30)RL=2.0 444 CDF=CF*(1.+RL*TOC)*2.*SRSTAR/SR 445 ARG=SQRT(COSC4) 446 XMTD(1)=0.60/ARG 447 XMTD(2)=0.65/ARG 448 XMTD(3)=0.70/ARG 449 XMTD(4)=0.75/ARG 450 XMTD(5)=0.80/ARG 451 XMTD(6)=0.85/ARG 452 XMTD(7)=0.90/ARG 453 XMTD(8)=0.925/ARG 454 XMTD(9)=0.950/ARG 455 XMTD(10)=0.975/ARG 456 XMTD(11)=1.000/ARG 457 XMTD(12)=1.050/ARG 458 XMTD(13)=1.100/ARG 459 XMTD(14)=1.200/ARG 460 XMTD(15)=1.400/ARG 461 ARG2=TOC**.3333 462 VAR(2)=ARSTAR*ARG2 463 ARG3=COSC4**2.5 464 DO 1090 I=1,15 465 ARG4=ABS((XMTD(I)*ARG)**2-1.) 466 VAR(1)=SQRT(ARG4)/ARG2 467C 468C FIGURE 4.1.5.1-29 469C 470 IF(XMTD(I)*ARG.LE.1.)GO TO 1080 471 LGH(1)=5 472 LGH(2)=7 473 CALL INTERX(2,T429R,VAR,LGH,D429R,CDW1(I),7,35, 474 1 0,0,0,0,0,2,0,0,Q5129R,3,ROUTID) 475 GO TO 1090 476 1080 LGH(1)=8 477 LGH(2)=7 478 CALL INTERX(2,T429L,VAR,LGH,D429L,CDW1(I),8,56, 479 1 0,0,0,0,-1,2,0,0,Q5129L,3,ROUTID) 480 1090 CDW2(I)=CDW1(I)*TOC**1.666*ARG3*SRSTAR/SR 481 CALL TRANF(15, XMTD, CDW2, 0.0, 0.0, MACH, CDW) 482 CD0W=CDW+CDF 483C 484C *** BODY ALONE *** 485C 486 1100 IF(.NOT.BO)RETURN 487C 488C ***BODY LIFT AND MOMENT SLOPES *** 489C 490 ARG1=(MACH-.60)/.80 491 CLABD=CLABM6+(CLAB14-CLABM6)*ARG1 492 CMAB=CMAB6+(CMAB14-CMAB6)*ARG1 493C 494C *** BODY DRAG *** 495C 496 CDFB=CFB6*SS/SR 497 CDPB=CD0B6-CDBB6-CDFB 498 IF(MACH.GE.1.0.AND.MACH.LE.1.2)CDPB=CDPB*(1.-(MACH-1.0)/.2) 499 IF(MACH.GT.1.2)CDPB=0.0 500 IF(DB.LT.0.3*DMAX)DB=0.3*DMAX 501 ARG=(DB/DMAX)**2 502C 503C FIGURE 4.2.3.1-24 (CDBT) 504C 505 VAR(1)=MACH 506 VAR(2)=CDBB6/ARG*4.0*SR/(PI*DMAX**2) 507 LGH(1)=13 508 LGH(2)=3 509 CALL INTERX(2,T424,VAR,LGH,D424,CDBFIG,13,39, 510 1 0,2,0,0,0,2,0,0,Q23124,3,ROUTID) 511 CDBB=CDBFIG*ARG*PI*DMAX**2/(4.0*SR) 512 IF(MACH.LT.1.0)GO TO 1110 513 VAR(1)=XLB/DMAX 514 VAR(2)=MACH 515 LGH(1)= 11 516 LGH(2)=5 517 CALL INTERX(2,T426,VAR,LGH,D426,CDWB,11,55, 518 1 0,0,0,0,2,2,0,0,Q23126,3,ROUTID) 519 CDWB=CDWB*PI*DMAX**2/(4.*SR) 520 GO TO 1120 521 1110 CDWB=0.0 522 1120 CD0B=CDFB+CDPB+CDBB+CDWB 523 IF(MACH.GT.1.2)GO TO 1130 524 GO TO 1140 525 1130 CD0B=CD0B+(CD014-CD0B)*(MACH-1.2)/0.2 526C 527C *** BODY DRAG AT ANGLE OF ATTACK *** 528C 529 1140 DO 1150 J=1,NALPHA 530 1150 CDJB(J)=CD0B+(FLC(J+22)/RAD)**2*SB/SR 531 CD0WB=CD0B+CD0W 532 TRA(73)=CD0WB 533 RETURN 534 END 535