1 SUBROUTINE SUBLAT(STB,WING,WINGIN,BW,A,B,AHT,AVT,VT,VTIN, 2 1 VTPL,TVTPAN,HTPL,WGPL,ITY) 3C 4C THIS SUBROUTINE COMPUTES THE SUBSONIC LATERAL STABILITY DERIVATIVES 5C 6 COMMON /OVERLY/ NLOG,NMACH,I,NALPHA,IG 7 COMMON /BDATA/ BD(762) 8 COMMON /IBODY/ PBODY,BODY(400) 9 COMMON /SYNTSS/ SYNA(19) 10 COMMON /VTI/ XXXX(154), TVTIN(8) 11 COMMON /HTI/ HTIN(154) 12 COMMON /FLOLOG/ FLTC,OPTI,BO,XX(14),TRANSN 13 COMMON /FLGTCD/ FLC(160) 14 COMMON /CONSNT/ PI,DEG,UNUSED,RAD 15 COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF 16C 17 LOGICAL FLTC,OPTI,BO,WGPL,TVTPAN,TRANSN,HTPL,VTPL 18C 19 DIMENSION STB(135),WING(380),WINGIN(100),BW(380),A(195) 20 DIMENSION B(49), AHT(195), AVT(195), VT(380), VTIN(100) 21 DIMENSION ROUTID(2),Q12127(3),Q2128A(3),Q2128B(3),Q2130B(3), 22 1 Q2130A(3),Q12129(3),Q12131(3),Q52318(3),Q22126(3), 23 2 Q1122A(3),Q1122B(3),Q1122C(3),Q1122D(3),Q1124A(3), 24 3 Q1124B(3),Q1124C(3) 25 DIMENSION WTYPE(4) 26 DIMENSION C4(6),C5(6),C6(6),C7(6),C8(6),C9(6),C1(6),C2(6) 27 DIMENSION X327(3),X227(8),X127(5),Y27(120) 28 DIMENSION X128A(7),X228A(10),Y28A(70) 29 DIMENSION X128B(3),X228B(9),Y28B(27) 30 DIMENSION X129(3),X229(9),X329(3),Y29(81) 31 DIMENSION X130A(5),X230A(9),Y30A(45) 32 DIMENSION X130B(4),X230B(9),Y30B(36) 33 DIMENSION X1526(7),X2526(9),Y526(63) 34 DIMENSION X131(4),X231(11),X331(3),X431(3),Y31(396),T31A(132), 35 1 T31B(132),T31C(132) 36 DIMENSION X158A(10),X258A(2),Y58A(20),X158B(5),X258B(3),Y58B(15), 37 1 X158C(5),X258C(2),Y58C(10) 38 DIMENSION X5322D(13),Y5322D(13),X122A(2),X222A(19),Y22A(38),X122B( 39 14),X222B(13),Y22B(52),X5322C(11),Y5322C(11),X5324A(11),Y5324A(11), 40 2X124B(2),X224B(10),Y24B(20),X124C(5),X224C(11),Y24C(55) 41C 42 EQUIVALENCE(Y31(1),T31A(1)),(Y31(133),T31B(1)),(Y31(265),T31C(1)) 43C 44 DATA ROUTID/4HSUBL,4HAT /,Q12127/4H5.1.,4H2.1-,4H27 /, 45 1 Q2128A/4H5.1.,4H2.1-,4H28-A/,Q2130B/4H5.1.,4H2.1-,4H30-B/, 46 2 Q2128B/4H5.1.,4H2.1-,4H28-B/,Q2130A/4H5.1.,4H2.1-,4H30-A/, 47 3 Q12129/4H5.1.,4H2.1-,4H29 /,Q12131/4H5.1.,4H2.1-,4H31 /, 48 4 Q52318/4H5.2.,4H3.1-,4H8 /,Q22126/4H5.2.,4H2.1-,4H26 /, 49 5 Q1122A/4H5.3.,4H1.1-,4H22A /,Q1124A/4H5.3.,4H1.1-,4H24A /, 50 6 Q1122B/4H5.3.,4H1.1-,4H22B /,Q1124B/4H5.3.,4H1.1-,4H24B /, 51 7 Q1122C/4H5.3.,4H1.1-,4H22C /,Q1124C/4H5.3.,4H1.1-,4H24C /, 52 8 Q1122D/4H5.3.,4H1.1-,4H22D / 53 DATA WTYPE/4HSTRA,4HDOUB,4HCRAN,4HCURV/ 54C 55C--------- FIGURE 5.1.2.1-27 56C 57 DATA X327/0.0,0.5,1.0/ 58 DATA X127/1.0,2.0,4.0,6.0,8.0/ 59 DATA X227/-20.,0.0,20.,30.,40.,50.,55.,60./ 60 DATA Y27/.0014,0.0,-.00125,-.002,-.0027,-.0036,-.004,-.0044 , 61 1.0015,0.0,-.00145,-.0022,-.003,-.0041,-.005,-.00595 , .0016,0.0, 62 2-.0016,-.0024,-.0033,-.0047,-.0057,-.0071 , .0016,0.0,-.0016, 63 3-.0024,-.0035,-.0049,-.006,-.0074 , .0016,0.0,-.0016,-.0027,-.0035 64 4,-.0049,-.006,-.0074 , .0012,0.0,-.0012,-.0019,-.0026,-.0034, 65 5-.0039,-.0044 , .0013,0.0,-.0013,-.0021,-.003,-.0043,-.00515, 66 6-.0064 , .0015,0.0,-.0014,-.0024,-.0036,-.005,-.00605,-.0075 , 67 7.00165,0.0,-.0016,-.0025,-.0038,-.0054,-.0066,-.0082 , .0018,0.0, 68 8-.00175,-.0027,-.004,-.0058,-.007,-.0089 , .00105,0.0,-.001,-.0016 69 9,-.0023,-.003,-.0035,-.0038 , .0012,0.0,-.0013,-.0021,-.0031, 70 A-.00435,-.00505,-.0062 , .0014,0.0,-.00165,-.00245,-.0036,-.0052, 71 B-.0061,-.0078 , .00167,0.0,-.0017,-.0028,-.004,-.00595,-.00715, 72 C-.009 , .0018,0.0,-.0018,-.00295,-.0042,-.0062,-.0078,-.010/ 73C 74C--------- FIGURE 5.1.2.1-28-A ------------ 75C 76 DATA X128A/2.,3.,4.,5.,6.,8.,10./ 77 DATA X228A/0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,0.95/ 78 DATA Y28A/8*1.0,0.995,0.990 , 6*1.0,1.01,2*1.03,1.02 , 2*1.0,1.01, 79 11.015,1.025,1.05,1.08,1.09,1.10,1.10 , 1.0,1.01,1.015,1.02,1.05, 80 21.09,1.115,1.16,1.20,1.21 , 1.0,1.01,1.02,1.04,1.07,1.12,1.17,1.24 81 3,1.32,1.36 , 1.0,1.01,1.05,1.07,1.12,1.18,1.27,1.40,1.58,1.70 , 82 41.0,1.02,1.05,1.10,1.15,1.23,1.37,1.54,1.84,2.08/ 83C 84C--------- FIGURE 5.1.2.1-28-B ------------ 85C 86 DATA X128B/0.0,0.50,1.0/ 87 DATA X228B/1.0,1.5,2.0,2.5,3.0,4.0,5.0,6.0,8.0/ 88 DATA Y28B/-.0058,-.00345,-.00235,-.00145,-.001,-.00045,-.00025, 89 1.00005,.0004 , -.008,-.00555,-.004,-.003,-.00235,-.0014,-.001, 90 2-.00065,-.0002 , -.0113,-.008,-.00595,-.00465,-.0037,-.00255, 91 3-.00182,-.00147,-.00097/ 92C 93C--------- FIGURE 5.1.2.1-29 ------------ 94C--------- NOTE% INPUT ABSOLUTE VALUE OF THE ------------ 95C--------- SWEEP ANGLE FOR FIG. 5.1.2.1-29 ------------ 96C--------- ( THE X129 ARGUMENT ) ------------ 97C 98 DATA X329/0.0,0.5,1.0/ 99 DATA X129/0.0,40.0,60.0/ 100 DATA X229/0.0,1.0,2.0,3.0,4.0,5.0,6.0,8.0,10./ 101 DATA Y29/0.0,-.000052,-.000088,-.00011,-.000134,-.000153,-.000168, 102 1-.00019,-.0002 , 0.0,-.000048,-.000085,-.000108,-.000128,-.000141, 103 2-.000153,-.000173,-.000178 , 0.0,-.00004,-.000073,-.000095, 104 3-.000108,-.000119,-.000127,-.000135,-.000138 , 0.0,-.000052, 105 4-.000098,-.000132,-.000162,-.000186,-.000208,-.00024,-.00026 , 106 50.0,-.00005,-.000096,-.000124,-.00015,-.00017,-.000188,-.000217, 107 6-.00023 , 0.0,-.00005,-.000087,-.000111,-.000129,-.000142, 108 7-.000153,-.000166,-.000170 , 0.0,-.00005,-.000096,-.000133, 109 8-.000167,-.000193,-.000216,-.000252,-.00028 , 0.0,-.00005,-.000095 110 9,-.000129,-.000155,-.000178,-.000197,-.000225,-.000245 , 0.0, 111 A-.00005,-.000088,-.000113,-.000132,-.000147,-0.000159,-.000172, 112 B-.00018/ 113C 114C--------- FIGURE 5.1.2.1-30-A 115C 116 DATA X130A/2.,4.,6.,8.,10./ 117 DATA X230A/0.0,.2,.4,.5,.6,.7,.8,.9,.95/ 118 DATA Y30A/1.0,1.01,1.018,1.02,1.023,1.03,1.04,1.05,1.057 , 1.0, 119 11.012,1.03,1.045,1.06,1.085,1.118,1.16,1.19 , 1.0,1.015,1.045, 120 21.07,1.1,1.14,1.197,1.27,1.33 , 1.0,1.018,1.05,1.085,1.125,1.19, 121 31.26,1.39,1.485 , 1.0,1.02,1.058,1.097,1.148,1.215,1.325,1.495, 122 41.635/ 123C 124C--------- FIGURE 5.1.2.1-30-B ------------ 125C 126 DATA X130B/0.0,0.4,0.6,1.0/ 127 DATA X230B/3.,4.,5.,6.,7.,8.,9.,10.,11./ 128 DATA Y30B/-.0000192,-.0000222,-.0000238,-.0000231,-.0000230, 129 1-.0000241,-.000026,-.0000284,-.0000328 , -.000022,-.0000287, 130 2-.0000323,-.0000335,-.0000339,-.0000342,-.000035,-.000037,-.000042 131 3 , -.0000233,-.00003,-.0000335,-.000035,-.0000366,-.000037, 132 4-.0000375,-.00004,-.000047 , -.0000233,-.00003,-.0000335,-.000035, 133 5-.0000366,-.000037,-.0000375,-.00004,-.000047/ 134C 135C--------- FIGURE 5.1.2.1-31 ------------ 136C 137 DATA X431 /0.0,0.5,1.0/ 138 DATA X331 /2.,4.,8./ 139 DATA X131 /-40.,0.0,40.,60./ 140 DATA X231 /0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/ 141 DATA T31A/0.,-.0001,-.00026,-.00056,-.001,-.0015,-.00217,-.0028, 142 1-.00336,-.0039,-.00435 , 0.,-.0001,-.00032,-.0006,-.0011,-.00176, 143 2-.0025,-.0032,-.00385,-.00438,-.0048 , 0.,-.0001,-.0003,-.00063, 144 3-.00117,-.0018,-.00256,-.0033,-.00395,-.00448,-.00496 , 0.0, 145 4-.0001,-.0003,-.0006,-.00115,-.0017,-.00236,-.0031,-.0037,-.00417, 146 5-.00450 , .0,-.0001,-.00025,-.00078,-.00153,-.0024,-.00335,-.0043, 147 6-.0051,-.00586,-.0066 , 0.,-.0001,-.00038,-.001,-.00183,-.0029, 148 7-.00405,-.0052,-.00633,-.00726,-.0079 , 0.,-.0001,-.0004,-.00102, 149 8-.00186,-.00303,-.00425,-.0054,-.00655,-.00736,-.00775 , 0.,-.0001 150 9,-.00040,-.001,-.00177,-.00275,-.0038,-.00476,-.00565,-.00635, 151 A-.00653 , 0.,-.00008,-.00057,-.00145,-.00255,-.00374,-.00484, 152 B-.00585,-.00682,-.0077,-.00845 , 0.,-.0001,-.00063,-.00165,-.003, 153 C-.00445,-.00625,-.00786,-.0093,-.01040,-.01110 , 0.,-.00013,-.0009 154 D,-.00195,-.00325,-.00483,-.00655,-.00815,-.0095,-.01035,-.0106 , 155 E0.,-.00013,-.0009,-.00195,-.00325,-.0046,-.00585,-.00697,-.00787, 156 F-.0084,-.00842/ 157 DATA T31B/0.,-.00015,-.0003,-.00065,-.0011,-.0017,-.0025,-.00325, 158 1-.004,-.00453,-.005 , 0.0,-.00015,-.00037,-.00076, 159 2-.00128,-.0019,-.00276,-.00355,-.00426,-.00493,-.00536 , 0., 160 3-.00015,-.00037,-.00076,-.00128,-.0019,-.00276,-.00355,-.00426, 161 4-.0048,-.0052 , 0.0,-.00014,-.00036,-.0007,-.00125,-.0019,-.0027, 162 5-.0034,-.004,-.00452,-.00508 , 0.,-.0001,-.0004,-.00096,-.00163, 163 6-.00255,-.0038,-.0051,-.00637,-.00743,-.00825 , 0.,-.0002,-.00065, 164 7-.00125,-.00215,-.0032,-.0045,-.00595,-.0072,-.00845,-.00935 , 0., 165 8-.0002,-.00065,-.00125,-.00215,-.0032,-.0045,-.00595,-.0072, 166 9-.00823,-.009 , 0.,-.0001,-.0005,-.00113,-.00185,-.00285,-.004, 167 A-.0051,-.00605,-.0069,-.00745 , 0.,-.00015,-.0006,-.00125,-.00215, 168 B-.00327,-.0047,-.00645,-.0084,-.0100,-.01140 , 0.,-.0003,-.00094, 169 C-.00192,-.00315,-.0046,-.0068,-.0089,-.0109,-.01255,-.0141 , 0., 170 D-.0003,-.0009,-.00195,-.00325,-.00482,-.00675,-.00855,-.01015, 171 E-.01155,-.01265 , 0.,-.00015,-.0009,-.00176,-.00285,-.00405,-.0053 172 F,-.00665,-.0078,-.00885,-.0097/ 173 DATA T31C/0.,-.00015,-.00038,-.0007,-.00126,-.0019,-.00268,-.00344 174 1,-.00418,-.00484,-.0053 , 0.,-.00016,-.0004,-.00076,-.0013,-.002, 175 2-.00288,-.0037,-.00445,-.0051,-.0056 , 0.,-.00015,-.00035,-.0007, 176 3-.00125,-.0019,-.00278,-.00356,-.0043,-.005,-.0054 , 0.,-.00012, 177 4-.0003,-.00067,-.0012,-.0018,-.00257,-.0033,-.00396,-.00458, 178 5-.00505 , 0.,-.00015,-.0005,-.00115,-.0018,-.00273,-.00385,-.00515 179 6,-.0066,-.0079,-.0088 , 0.,-.00015,-.00052,-.0012,-.0021,-.0033, 180 7-.00486,-.00627,-.00756,-.0087,-.00972 , 0.,-.00015,-.00052,-.0012 181 8,-.0021,-.0033,-.00465,-.006,-.00728,-.00836,-.00925 , 0.,-.00015, 182 9-.0005,-.00116,-.0019,-.00285,-.004,-.0051,-.0061,-.00705,-.0079 , 183 A-0.,-.00018,-.0006,-.00125,-.0021,-.00335,-.00505,-.0071,-.00915, 184 B-.01100,-.01255 , 0.,-.00025,-.0009,-.0018,-.0031,-.0048,-.007, 185 C-.0094,-.01145,-.0133,-.01475 , 0.,-.0003,-.0009,-.00195,-.00327, 186 D-.0048,-.00697,-.0089,-.01057,-.01205,-.0132 , 0.,-.00012,-.00073, 187 E-.00155,-.00265,-.004,-.00536,-.0067,-.00795,-.009,-.010/ 188C 189C FIGURE 5.2.2.1-26 190C 191 DATA X1526/4.0,4.5,5.0,5.5,6.0,7.0,8.0/ 192 DATA X2526/0.0,.2,.4,.6,.8,1.0,1.2,1.4,1.6/ 193 DATA Y526/1.0,1.0,1.0,1.0,1.0,1.0,1.0,.99,.97, 194 11.0,1.0,1.0,1.0,1.0,1.0,.98,.948,.911, 195 21.0,1.0,1.0,1.0,.997,.971,.933,.883,.827, 196 31.0,1.0,1.0,.991,.963,.922,.870,.811,.746, 197 41.0,1.0,.995,.970,.932,.884,.829,.764,.695, 198 51.0,1.0,.977,.944,.899,.845,.780,.715,.641, 199 61.0,.985,.960,.921,.870,.812,.745,.670,.592/ 200C 201C FIGURE 5.2.3.1-8A 202C 203 DATA X158A/20.,14.,10.,8.,7.,6.,5.,4.,3.,2.5/ 204 DATA X258A/.2,.8/ 205 DATA Y58A/.1,1.88, 206 1.40,2.21, 207 2.74,2.60, 208 3.98,2.80, 209 41.30,3.13, 210 51.61,3.50, 211 62.00,3.88, 212 72.50,4.40, 213 82.99,5.00, 214 93.45,5.40/ 215C 216C FIGURE 5.2.3.1-8B 217C 218 DATA X158B/.8,1.0,1.2,1.4,1.6/ 219 DATA X258B/0.0,3.0,6.0/ 220 DATA Y58B/0.0,2.35,4.68, 221 10.0,3.00,6.00, 222 20.0,3.60,7.25, 223 30.0,4.18,8.50, 224 40.0,4.79,9.50/ 225C 226C FIGURE 5.2.3.1-8C 227C 228 DATA X158C/.5,.6,.8,1.,2./ 229 DATA X258C/0.0,6.0/ 230 DATA Y58C/-.00048,.00251, 231 1-.00048,.0035, 232 2-.00048,.00477, 233 3-.00048,.00559, 234 4-.00048,.00641/ 235C 236C FIGURE 5.3.1.1-22A 237C 238 DATA X122A/1.0,.6/ 239 DATA X222A/0.,.125,.25,.5,.75,1.,1.25,1.5,1.75,2.,2.25,2.5, 240 13.00,3.25,3.50,3.75,4.00,5.00,7.00/ 241 DATA Y22A/0.,.40,.720,.990,1.19,1.32,1.40,1.46,1.50,1.51,1.48,1.42 242 1,1.27,1.21,1.17,1.13,1.10,1.04,1.02, 243 20.,.70,.940,1.18,1.35,1.46,1.54,1.60,1.63,1.64,1.60,1.53,1.36, 244 31.28,1.21,1.16,1.13,1.06,1.02/ 245C 246C FIGURE 5.3.1.1-22B 247C 248 DATA X122B/.5,.6,.7,.8/ 249 DATA X222B/0.0,-.2,-.3,-.4,-.5,-.6,-.65,-.7,-.75,-.8,-.85,-.9,-1./ 250 DATA Y22B/1.05,.94,.90,.87,.86,.87,.90,.93,.98,1.06,1.16,1.29,1.7, 251 11.15,1.00,.95,.90,.89,.90,.92,.96,1.01,1.08,1.18,1.31,1.70, 252 21.22,1.05,.99,.94,.92,.92,.95,.98,1.03,1.10,1.20,1.33,1.70, 253 31.29,1.09,1.02,.97,.94,.94,.96,1.0,1.06,1.12,1.22,1.36,1.70/ 254C 255C FIGURE 5.3.1.1-22C 256C 257 DATA X5322C/0.0,.2,.4,.6,.7,.8,.9,1.2,1.4,1.6,2.0/ 258 DATA Y5322C/0.0,.29,.52,.70,.77,.83,.87,.98,1.04,1.07,1.13/ 259C 260C FIGURE 5.3.1.1-22D 261C 262 DATA X5322D/0.0,.50,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0/ 263 DATA Y5322D/5*.75,.835,.92,6*1.0/ 264C 265C FIGURE 5.3.1.1-24A 266C 267 DATA X5324A/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/ 268 DATA Y5324A/1.5,1.33,1.197,1.098,1.03,1.0,1.03,1.098,1.197,1.33, 269 11.5/ 270C 271C FIGURE 5.3.1.1-24B 272C 273 DATA X124B/20.0,0.0/ 274 DATA X224B/.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0/ 275 DATA Y24B/.90,1.64,2.19,2.63,2.95,3.22,3.44,3.61,3.75,3.86, 276 11.0,1.72,2.30,2.75,3.11,3.39,3.62,3.82,3.98,4.11/ 277C 278C FIGURE 5.3.1.1-24C 279C 280 DATA X124C/.2,.4,.6,.8,1.0/ 281 DATA X224C/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/ 282 DATA Y24C/1.0,.90,.818,.741,.68,.631,.593,.563,.542,.527,.510, 283 11.0,.93,.87,.825,.784,.752,.730,.710,.696,.688,.680, 284 21.0,.95,.911,.882,.859,.844,.830,.820,.812,.805,.800, 285 31.0,.968,.941,.923,.908,.898,.890,.883,.881,.880,.879, 286 41.0,.979,.962,.954,.949,.943,.938,.935,.935,.931,.931/ 287 DATA I4,I5,I6/3*0/ 288C 289 IF(.NOT. WGPL) GO TO 1255 290 YA27 =UNUSED 291 YA28A=UNUSED 292 YA28B=UNUSED 293 YA29 =UNUSED 294 YA30A=UNUSED 295 YA30B=UNUSED 296 YA31I=UNUSED 297 YA31O=UNUSED 298C 299C LATERAL STABILITY DERIVATIVES 300C 301 IF(TRANSN)GO TO 1190 302C 303C ***** WING ALONE CALCULATIONS ***** 304C 305C WING ALONE SIDEFORCE DERIVATIVE 306C 307 BDAT=2.*WINGIN(4) 308 ANUM1=6.0*A(44)*A(42) 309 DEN1=PI*A(120)*(A(120)+4.*A(43)) 310 BB=SQRT(1.-(FLC(I+2)*A(43))**2) 311 ANUM2=A(120)+4.0*A(43) 312 DEN2=A(120)*BB+4.0*A(43) 313 DO 1050 J=1,NALPHA 314 STB(J+15)=(1./RAD)*(B(J+2)**2)*(ANUM1/DEN1)-.0001*ABS(STB(122)) 315 IF(B(J+2).EQ.0.0)GO TO 1000 316 STB(J+35)=(ANUM2/DEN2)*(STB(J+15)/B(J+2)) 317 WING(J+140)=STB(J+35)*WING(J+20) 318 GO TO 1010 319 1000 STB(J+35)=0.0 320 WING(J+140)=0.0 321 1010 CONTINUE 322C 323C WING ALONE YAWING MOMENT DERIVATIVE 324C 325 IF(J.GT.1)GO TO 1020 326 APART=1.0/(4.0*PI*A(120)) 327 BPART=ABS(A(44))/(PI*A(120)*(A(120)+4.0*A(43))) 328 CPART=A(43)-0.5*A(120)-A(120)**2/(8.*A(43)) 329 EPART=(A(120)+4.*A(43))/(A(120)*BB+4.*A(43)) 330 ANUM3=(A(120)*BB)**2+4.*A(120)*BB*A(43)-8.0*A(43)*A(43) 331 DEN3=A(120)*A(120)+4.0*A(120)* A(43)-8.0*A(43)*A(43) 332 1020 CONTINUE 333 IF(WING(J+60).EQ.0.0)GO TO 1030 334 DPART=6.0*(WING(J+40)/WING(J+60)/CBARR)*(ABS(A(42))/A(120)) 335 GO TO 1040 336 1030 DPART=0.0 337 1040 CONTINUE 338 STB(J+75)=(APART-BPART*(CPART+DPART))/RAD 339 WING(J+160)=EPART*(ANUM3/DEN3)*STB(J+75)*WING(J+20)**2 340 WING(J+160)=WING(J+160)*BDAT*SREF/(BLREF*A(3)) 341 1050 CONTINUE 342C 343C--------- TEST FOR CONFIGURATION AND UNIFORM DIHEDRAL ANGLE ------- 344C 345C WING ALONE ROLLING MOMENT DERIVATIVE 346C (FOR STRAIGHT WING CONFIGURATION WITH UNIFORM DIHEDRAL) 347C 348 IF(WINGIN(15).NE.WTYPE(1)) GO TO 1120 349 IF(A(120).GE.1.0) GO TO 1070 350 DO 1060 J=1,NALPHA 351 WING(J+180)=-2./(RAD*3.*A(120))*WING(J+20)-STB(122)*A(120)/ 352 1 (6.*RAD**2)*(A(4)/SREF) 353 1060 WING(J+180)=WING(J+180)*BDAT/ BLREF 354 GO TO 1180 355 1070 CONTINUE 356 XA2= FLC(I+2)*A(73) 357 X1ARG=ABS(A(70)) 358 XA1=A(120)/A(73) 359C 360C --------- FIGURE 5.1.2.1-27 ---------- 361C 362 CALL TLIN3X(X127,X227,X327,Y27,5,8,3,A(120),X1ARG,A(118),YA27, 363 1 0,1,0,2,2,0,Q12127,3,ROUTID) 364 IF(A(70).LT.0.0) YA27=ABS(YA27) 365C 366C --------- FIGURE 5.1.2.1-28-A ---------- 367C 368 CALL TLINEX(X128A,X228A,Y28A,7,10,XA1,XA2,YA28A, 369 1 0,2,1,2,Q2128A,3,ROUTID) 370C 371C --------- FIGURE 5.1.2.1-28-B ---------- 372C 373 CALL TLINEX(X128B,X228B,Y28B,3,9,A(118),A(120),YA28B, 374 1 0,0,0,0,Q2128B,3,ROUTID) 375C 376C --------- FIGURE 5.1.2.1-30-B ---------- 377C 378 CALL TLINEX(X130B,X230B,Y30B,4,9,A(118),A(120),YA30B, 379 1 0,2,0,2,Q2130B,3,ROUTID) 380 DO 1080 J=1,NALPHA 381 1080 STB(J+95)=(YA27*YA28A+YA28B)*WING(J+20)*SREF/A(4)+WINGIN(11)* 382 1 (A(68))*YA30B 383 IF(WINGIN(12).NE.UNUSED.AND.WINGIN(12).NE.0.0)GO TO 1100 384C 385C --------- FIGURE 5.1.2.1-30-A ---------- 386C 387 CALL TLINEX(X130A,X230A,Y30A,5,9,XA1,XA2,YA30A, 388 1 2,2,2,2,Q2130A,3,ROUTID) 389C 390C --------- FIGURE 5.1.2.1-29 ---------- 391C 392 CALL TLIN3X(X129,X229,X329,Y29,3,9,3,X1ARG,A(120),A(118),YA29, 393 1 0,0,0,2,2,0,Q12129,3,ROUTID) 394 DO 1090 J=1,NALPHA 395 WING(J+180)=(STB(J+95)+WINGIN(13)*YA29*YA30A)*A(4)/SREF 396 WING(J+180)=WING(J+180)*BDAT/BLREF 397 1090 CONTINUE 398 GO TO 1180 399C 400C WING ALONE ROLLING MOMENT DERIVATIVE 401C (FOR STRAIGHT WING CONFIGURATION WITH NON-UNIFORM DIHEDRAL) 402C 403 1100 CONTINUE 404 STB(7)=A(131)/(2.0*PI)*RAD 405 X1ARG=ATAN(A(68)/B(2))*RAD 406 X3ARG=B(2)*A(120)/STB(7) 407 AKOVB=STB(7)/B(2) 408C 409C --------- FIGURE 5.1.2.1-31 ---------- (INBOARD) 410C 411 CALL TLIN4X(X131,X231,X331,X431,Y31,4,11,3,3,X1ARG,STB(2),X3ARG, 412 1 A(118),YA31I,0,0,1,0,2,2,1,0,Q12131,3,ROUTID) 413C 414C --------- FIGURE 5.1.2.1-31 ---------- (OUTBOARD) 415C 416 CALL TLIN4X(X131,X231,X331,X431,Y31,4,11,3,3,X1ARG,STB(3),X3ARG, 417 1 A(118),YA31O,0,0,1,0,2,2,1,0,Q12131,3,ROUTID) 418 DO 1110 J=1,NALPHA 419 1110 WING(J+180)=(STB(J+95)+(YA31O-YA31I)*AKOVB*WINGIN(14)/RAD+YA31I/ 420 1 RAD*AKOVB*WINGIN(13))*A(4)/SREF*BDAT/BLREF 421 GO TO 1180 422C 423C WING ALONE ROLLING MOMENT DERIVATIVE 424C DOUBLE DELTA OR CRANKED WING CONFIGURATION 425C (INBOARD PANEL ROLLING MOMENT) 426C 427 1120 CONTINUE 428 X1ARG=ABS(A(70)) 429 STB(135)=-2.0/(3.0*RAD*A(163)) 430C 431C --------- FIGURE 5.1.2.1-27 ----------- 432C 433 CALL TLIN3X(X127,X227,X327,Y27,5,8,3,A(163),X1ARG,A(25),STB(132), 434 1 0,1,0,2,2,0,Q12127,3,ROUTID) 435 IF(A(70).LT.0.0) YA27=ABS(YA27) 436C 437C --------- FIGURE 5.1.2.1-28-A ---------- 438C 439 XA1I=A(163)/A(73) 440 XA2I=A(73)*FLC(I+2) 441 CALL TLINEX(X128A,X228A,Y28A,7,10,XA1I,XA2I,STB(134), 442 1 0,2,1,2,Q2128A,3,ROUTID) 443C 444C --------- FIGURE 5.1.2.1-28-B ---------- 445C 446 CALL TLINEX(X128B,X228B,Y28B,3,9,A(25),A(163),STB(133), 447 1 0,0,0,0,Q2128B,3,ROUTID) 448C 449C (OUTBOARD PANEL ROLLING MOMENT) 450C 451 STB(131)=-2.0/(3.*RAD*A(168)) 452 X1ARG=ABS(A(94)) 453C 454C --------- FIGURE 5.1.2.1-27 ---------- 455C 456 CALL TLIN3X(X127,X227,X327,Y27,5,8,3,A(168),X1ARG,A(169),STB(128), 457 1 0,1,0,2,2,0,Q12127,3,ROUTID) 458 IF(A(94).LT.0.0) YA27=ABS(YA27) 459C 460C --------- FIGURE 5.1.2.1-28-A ---------- 461C 462 XA1O=A(168)/A(97) 463 XA2O=A(97)*FLC(I+2) 464 CALL TLINEX(X128A,X228A,Y28A,7,10,XA1O,XA2O,STB(130), 465 1 0,2,1,2,Q2128A,3,ROUTID) 466C 467C --------- FIGURE 5.1.2.1-28-B ---------- 468C 469 CALL TLINEX(X128B,X228B,Y28B,3,9,A(169),A(168),STB(129), 470 1 0,0,0,0,Q2128B,3,ROUTID) 471 STB(58)=A(171)*(A(1)/SREF)+A(172)*(A(167)/SREF) 472 IF(A(163).GE.1.0) GO TO 1130 473 STB(127)=A(171)*(A(1)/SREF)*STB(135)*(WINGIN(3)-WINGIN(2))/ 474 1 WINGIN(4) 475 GO TO 1140 476 1130 STB(127)=A(171)*(A(1)/SREF)*(STB(132)*STB(134)+STB(133))* 477 1(WINGIN(3)-WINGIN(2))/WINGIN(4) 478 1140 IF(A(168).GE.1.0)GO TO 1150 479 STB(126)=A(172)*(A(167)/SREF)*STB(131) 480 GO TO 1160 481 1150 STB(126)=A(172)*(A(167)/SREF)*(STB(128)*STB(130)+STB(129)) 482 1160 CONTINUE 483 DO 1170 J=1,NALPHA 484 1170 WING(J+180)=(WING(J+20)/STB(58))*(STB(127)+STB(126))*BDAT/BLREF 485 1180 CONTINUE 486 STB(68)=YA27 487 STB(71)=YA28A 488 STB(70)=YA28B 489 STB(67)=YA29 490 STB(66)=YA30A 491 STB(69)=YA30B 492 STB(64)=YA31I 493 STB(65)=YA31O 494 IF(BO)GO TO 1190 495 GO TO 1330 496C 497C ***** WING BODY CALCULATIONS ***** 498C 499C WING-BODY SIDEFORCE DERIVATIVE(ALL PLANFORMS) 500C 501 1190 CONTINUE 502 ARG=2.0*STB(1)/STB(72) 503 IF(ARG.LT.0.0)GO TO 1200 504C 505C FIGURE 5.2.1.1-7 506C 507 STB(57)=1.0+.49*ARG 508 GO TO 1210 509 1200 STB(57)=1.0-.85*ARG 510 1210 BW(141)=-STB(57)*BODY(101)-0.0001*ABS(STB(122)) 511C 512C WING-BODY YAWING MOMENT DERIVATIVE(ALL PLANFORMS) 513C 514 ARG1=FLC(I+42)*BD(1) 515 STB(56)=1.+ALOG(1.E-6*ARG1)/4.86 516 ARG11=STB(59)/STB(123) 517 ARG12=SQRT(STB(61)/STB(60)) 518 ARG13=BD(1)**2/STB(62) 519 ARG14=SYNA(1)/BD(1) 520C FIGURE 5.2.3.1-8 521C 522 CALL TLINEX(X158A,X258A,Y58A,10,2,ARG13,ARG14,YDUMY, 523 1 2,1,2,1,Q52318,3,ROUTID) 524 CALL TLINEX(X158B,X258B,Y58B,5,3,ARG12,YDUMY,YDUMY2, 525 1 2,0,2,1,Q52318,3,ROUTID) 526 CALL TLINEX(X158C,X258C,Y58C,5,2,ARG11,YDUMY2,STB(15), 527 1 2,0,2,1,Q52318,3,ROUTID) 528 BW(161)=-STB(15)*STB(56)*(STB(62)*BD(1))/(SREF*BLREF) 529 IF(.NOT.TRANSN)GO TO 1220 530 RETURN 531 1220 CONTINUE 532C 533C WING-BODY ROLLING MOMENT DERIVATIVE(ALL PLANFORMS,UNIFORM AND NON- 534C UNIFORM DIHEDRAL) 535C 536 ANUM=1.2*STB(4)*STB(72)*SQRT(A(120)) 537 ADEN=RAD*2.0* WINGIN(4)**2 538 STB(14) = ANUM/ADEN 539 STB(13)=-.0005*SQRT(A(120))*STB(122)*(STB(72)/(2.0*WINGIN(4)))**2 540C 541C FIGURE 5.2.2.1-26 542C 543 ARG1=A(120)/A(49) 544 ARG2=STB(63)/(2.0* WINGIN(4)) 545 CALL TLINEX(X1526,X2526,Y526,7,9,ARG1,ARG2,STB(8), 546 1 0,0,0,1,Q22126,3,ROUTID) 547C 548 IF(WINGIN(15).EQ.WTYPE(1))GO TO 1230 549 APART= STB(8)*(STB(132)*STB(134)+STB(128)*STB(130))+STB(133) 550 1 +STB(129) 551 FPART=0.0 552 EPART=STB(13) 553 GO TO 1240 554 1230 CONTINUE 555 APART=STB(68)*STB(71)*STB(8)+STB(70) 556 FPART=WINGIN(11)*A(68)*STB(69) 557 EPART=WING(181)*BLREF/BDAT-(STB(96)-STB(13))*A(4)/SREF 558 1240 CONTINUE 559 DO 1250 J=1,NALPHA 560 1250 BW(J+180)=(BW(J+20)*APART+EPART+(STB(14)+FPART)*A(4)/SREF)* 561 1 BDAT/BLREF 562 1255 CONTINUE 563C 564C *****WING-BODY-TAIL CALCULATIONS***** 565C 566 IF(VTPL.OR.TVTPAN)GO TO 1260 567 GO TO 1330 568 1260 CONTINUE 569 VT(141) = 0.0 570 VT(161) = 0.0 571 DO 1270 J=1,NALPHA 572 1270 VT(J+180) = 0.0 573 IF(.NOT.VTPL) GO TO 1310 574C 575C SINGLE VERTICAL PANEL ON PLANE OF SYMMETRY 576C FIGURE 5.3.1.1-22A 577C 578 ARG2=VTIN(4)/(2.0*(VTIN(4)-VTIN(3))) 579 CALL TLINEX(X122A,X222A,Y22A,2,19,AVT(118),ARG2,STB(120), 580 1 0,0,0,0,Q1122A,3,ROUTID) 581C 582C FIGURE 5.3.1.1-22D 583C 584 CALL TBFUNX(ARG2,STB(118),DYDX,13,X5322D,Y5322D,C6,I6,MI,NG, 585 1 0,0,Q1122D,3,ROUTID) 586 ZH = SYNA(7)-((HTIN(4)-HTIN(3))*AHT(62)+AHT(30)-AHT(16)/4.0)* 587 1 SIN(SYNA(8)/RAD) 588C 589 IF(.NOT.HTPL.OR.ZH.LT.0.0) GO TO 1280 590C 591C FIGURE 5.3.1.1-22B 592C 593 ARG1=STB(9)/STB(10) 594 ARG2=-ZH/ VTIN(4) 595 CALL TLINEX(X122B,X222B,Y22B,4,13,ARG1,ARG2,STB(121), 596 1 0,0,2,0,Q1122B,3,ROUTID) 597C 598C FIGURE 5.3.1.1-22C 599C 600 ARG=(AHT(2)+AHT(119))/AVT(4) 601 CALL TBFUNX(ARG,STB(119),DYDX,11,X5322C,Y5322C,C5,I5,MI,NG, 602 1 0,0,Q1122C,3,ROUTID) 603C 604 GO TO 1290 605 1280 STB(121)=0.0 606 STB(119)=0.0 607 1290 STB(116)=STB(120)*AVT(120)*(1.+STB(119)*(STB(121)-1.)) 608 BPART=3.06*((AVT(4)/A(4))/(1.+A(43))) 609 CPART=.4* STB(4)/STB(59)+ .009*A(120) 610 STB(117)=0.724+BPART+CPART 611 ANUM=2.*PI*STB(116) 612 APART=STB(116)**2/(VTIN(I+20)*RAD/(2.*PI))**2 613 BPART= 1.+AVT(50)**2/B(2)**2 614 STB(5)=ANUM/(2.+SQRT(APART*BPART+4.)) 615 VT(141)=-STB(118)*STB(5)*STB(117)*AVT(4)/(RAD*SREF) 616 VT(161)=-VT(141)*STB(11)/BLREF 617 DO 1300 J=1,NALPHA 618 1300 VT(J+180)=(VT(141)*(STB(12)*COS(FLC(J+22)/RAD)-STB(11)* 619 1 SIN(FLC(J+22)/RAD)))/BLREF 620 IF(SYNA(18).EQ.UNUSED.AND.SYNA(19).EQ.UNUSED) GO TO 1310 621C 622C TWIN VERTICAL TAIL CORRECTION 623C 624 IF(ITY.EQ.0) TVCORR=2*(COS(SYNA(18)/RAD))**2 625 IF(ITY.EQ.1) TVCORR=2*(COS(SYNA(19)/RAD))**2 626 VT(141)=TVCORR*VT(141) 627 VT(161)=TVCORR*VT(161) 628 DO 1305 J=1,NALPHA 629 VT(J+180)=TVCORR*VT(J+180) 630 1305 CONTINUE 631 1310 CONTINUE 632 IF(.NOT.TVTPAN) GO TO 1330 633C 634C FIGURE 5.3.1.1-24A 635C 636 ARG =TVTIN(1)/TVTIN(2) 637 CALL TBFUNX(ARG ,STB(75),DYDX,11,X5324A,Y5324A,C4,I4,MI,NG, 638 1 0,0,Q1124A,3,ROUTID) 639 STB(6)=TVTIN(2)**2/TVTIN(5) 640C 641C FIGURE 5.3.1.1-24B 642C 643 ARG2=STB(75)*STB(6) 644 STB(116) = ARG2 645 ARG1=TVTIN(6) 646 CALL TLINEX(X124B,X224B,Y24B,2,10,ARG1,ARG2,STB(73), 647 1 0,2,0,2,Q1124B,3,ROUTID) 648C 649C FIGURE 5.3.1.1-24C 650C 651 ARG1=TVTIN(4)/BD(1) 652 ARG2=TVTIN(3)/TVTIN(2) 653 CALL TLINEX(X124C,X224C,Y24C,5,11,ARG1,ARG2,STB(74), 654 1 2,0,0,0,Q1124C,3,ROUTID) 655C 656 TVT=-STB(74)*STB(73)*2.*TVTIN(5)/(SREF*RAD) 657 VT(141)=VT(141)+TVT 658 VT(161)=VT(161)-TVT*TVTIN(7)/BLREF 659C 660 DO 1320 J=1,NALPHA 661 1320 VT(J+180)=VT(J+180)+TVT*(TVTIN(8)*COS(FLC(J+22)/RAD)-TVTIN(7)* 662 1 SIN(FLC(J+22)/RAD))/BLREF 663 1330 CONTINUE 664 RETURN 665 END 666