1 SUBROUTINE LIFTFP 2C 3C*** COMPUTES INCREMENTAL WING LIFT DUE TO FLAPS 4C 5 COMMON /OVERLY/ NLOG,NMACH,M,NALPHA 6 COMMON /FLGTCD/ FLC(93) 7 COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 8 1 HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN, 9 2 SUPERS,SUBSON,TRANSN,HYPERS, 10 3 SYMFP,ASYFP,TRIMC,TRIM 11 COMMON /CONSNT/ PI,DR,UNUSED,RAD 12 COMMON /FLAPIN/ F(116) 13 COMMON /POWR/ PW(104),FLP(189) 14 COMMON /OPTION/ SREF,CBARR,RUFF,BLREF 15 COMMON /WINGI/ WINGIN(77) 16 COMMON /SUPWH/ FCM(287) 17 COMMON /HTI/ HTIN(131) 18 COMMON /WINGD/ A(195) 19 COMMON /HTDATA/ AHT(195) 20 COMMON /IWING/ PWING, WING(400) 21 COMMON /IHT/ PHT, HT(380) 22 COMMON /SBETA/ STB(135), TRA(108), TRAH(108) 23 REAL KFPRM,MACH 24 REAL KSWEEP 25 LOGICAL TRANSL 26 LOGICAL EXPDCL 27 DIMENSION ROUTID(2),DSCLMX(10),ALDAG(40) 28 DIMENSION CPI(10),CPO(10),DALPDE(5) 29 DIMENSION SDCL(10),ADCAD(10) 30 DIMENSION RK2(10),DELCL2(10),DCLMAX(10) 31 DIMENSION VAR(4),LGH(4),RKB(5),CF(5),ETA(5),CHRD(5),DKB(4), 32 1 CFOC(5),DELTA(10),DELCL(5),CP(5),ALDAVG(4),ADCADS(5), 33 2 DCL(10),DCLK(5) 34 DIMENSION SWF(4),DELCLA(40),CPOCF(4),CLAD(10),CFACTR(4),CFACT(10) 35 DIMENSION CLDTHY(5),CLDOCT(5),ALPHAD(5) 36 DIMENSION ITRANS(5) 37 DIMENSION X2128A(11),X1128A(3),Y1128A(33),F1128A(3) 38 DIMENSION X1125A(8),X2125A(8),Y1125A(64),F1125A(3) 39 DIMENSION X1125B(16),X2125B( 6),Y1125B(96),F1125B(3) 40 DIMENSION X11126(7),X21126(14),Y11126(98),F11126(3) 41 DIMENSION X11419(3),X21419(12),Y61419(36),F61419(3) 42 DIMENSION X11127(5),X21127(13),Y11127(65),F11127(3) 43 DIMENSION X11150(7),Y11150(7),F11150(3) 44 DIMENSION X11418(10),X21418(10),Y61418(100),F61418(3) 45 DIMENSION X137AD(17),Y137AD(17),F137AD(3) 46 DIMENSION X137B2(2),Y137B2(2),F137B2(3) 47 DIMENSION X137B1(13),Y137B1(13),F137B1(3) 48 DIMENSION X137AB(12),Y137AB(12),F137AB(3) 49 DIMENSION X137AC(12),Y137AC(12),F137AC(3) 50 DIMENSION X1138A(12),Y138A1(12),Y138A2(12),Y138A3(12),Y138A4(12), 51 1 F138A1(3),F138A2(3),F138A3(3),F138A4(3) 52 DIMENSION X1138B(9),Y138B1(9),Y138B2(9),F138B1(3),F138B2(3) 53 DIMENSION X1418A(12),Y1418A(12),F1418A(3) 54 DIMENSION CF2(5),X61142(11),Y61142(11),F61142(3), 55 1 X6143A(14),Y6143A(21),F6143A(3),CF2OC(5),DF2(10), 56 2 CAPI(10),CAPO(10),X6143B(14),Y6143B(21),F6143B(3) 57 DIMENSION X11147(20),Y11147(40),F11147(3) 58 LOGICAL FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 59 1 HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN, 60 2 SUPERS,SUBSON,TRANSN,HYPERS, 61 3 SYMFP,ASYFP,TRIMC,TRIM 62 EQUIVALENCE (CFOC(1),FLP(62)),(CF2I,F(115)),(CF2O,F(116)), 63 1 (CAPI(1),F(85)),(CAPO(1),F(95)),(DF2(1),F(105)) 64 EQUIVALENCE (FTYPE,F(17)),(BOF,F(15)),(BIF,F(14)),(CFI,F(12)) 65 1,(CFO,F(13)),(RDELTA,F(16)),(SDCL(1),F(19)),(CPI(1),F(39)) 66 2,(CPO,F(49)),(TANPHE,F(11)),(DELTA(1),F(1)) 67 EQUIVALENCE (FLP(1),ETA(1)),(CHRD(1),FLP(6)),(CF(1),FLP(11)) 68 1,(DKB(1),FLP(20)),(SWF(1),FLP(24)),(CP(1),FLP(28)) 69 2,(CLOCLT,FLP(33)),(CLDOCT(1),FLP(34)),(CLDTHY(1),FLP(39)), 70 3(ALDAVG(1),FLP(16)),(DALPDE(1),FLP(54)),(DELCLA(1),FLP(110)) 71 4,(CFACT(1),FLP(71)),(ADCADS(1),FLP(67)),(DCL(1),WING(201)) 72 5,(DCLMAB,FLP(102)),(RK2(1),FLP(91)),(DELCL2(1),FLP(44)) 73 6,(CLAD(1),WING(241)) 74 EQUIVALENCE (DELN4,FLP(60)),(DELCL(1),FCM(283)) 75 EQUIVALENCE (ALDAG(1),FLP(150)),(CFOCA,FLP(61)) 76 1 ,(RK1,FLP(101)),(RK3,FLP(103)),(DSCLMX(1),FLP(81)), 77 2 (KSWEEP,FLP(104)),(DCLMAX(1),WING(221)),(ALPHAD(1),FLP(105)) 78 DATA ROUTID/4HLIFT,4HFP / 79C 80C FIGURE 4.1.1.2-8A 81C 82 DATA X2128A 83 1 /0., 0.02, 0.04, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 84 2 0.18, 0.20 / 85 DATA X1128A 86 1 /6., 7., 8./ 87 DATA Y1128A 88 1 / 0.900, 0.878, 0.858, 0.836, 0.815, 0.794, 0.772, 0.750, 89 2 0.728, 0.708, 0.685, 90 3 0.950, 0.938, 0.924, 0.907, 0.894, 0.878, 0.860, 0.842, 91 4 0.822, 0.802, 0.780, 92 5 0.966, 0.957, 0.947, 0.936, 0.924, 0.910, 0.896, 0.880, 93 6 0.862, 0.842, 0.822 / 94 DATA F1128A/4H4.1.,4H1.2-,4H8A / 95C 96C FIGURE 6.1.1.1-39A 97C 98 DATA X2125A 99 1 / 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.40, 0.50 / 100 DATA X1125A 101 1 / 0.00, 0.02, 0.04, 0.06, 0.08, 0.10, 0.12, 0.15 / 102 DATA Y1125A 103 1 / 1.770, 2.500, 3.000, 3.460, 3.820, 4.160, 4.690, 5.140, 104 2 1.770, 2.515, 3.030, 3.500, 3.873, 4.220, 4.780, 5.240, 105 3 1.770, 2.530, 3.060, 3.540, 3.926, 4.290, 4.870, 5.350, 106 4 1.770, 2.545, 3.090, 3.580, 3.979, 4.350, 4.950, 5.460, 107 5 1.770, 2.560, 3.120, 3.620, 4.032, 4.400, 5.040, 5.560, 108 6 1.770, 2.575, 3.150, 3.660, 4.085, 4.480, 5.120, 5.690, 109 7 1.770, 2.590, 3.180, 3.700, 4.138, 4.550, 5.210, 5.790, 110 8 1.770, 2.600, 3.220, 3.740, 4.190, 4.620, 5.330, 5.960/ 111 DATA F1125A/4H6.1.,4H1.1-,4H39A / 112C 113C FIGURE 6.1.1.1-39B 114C 115 DATA X2125B 116 1 / 0.05, 0.10, 0.15, 0.20, 0.25, 0.50 / 117 DATA X1125B 118 1 / 0.70, 0.72, 0.74, 0.76, 0.78, 0.80, 0.82, 0.84, 0.86, 119 2 0.88, 0.90, 0.92, 0.94, 0.96, 0.98, 1.00 / 120 DATA Y1125B 121 1 / .356, .382, .409, .431, .452, .548, 122 2 .399, .426, .452, .477, .498, .583, 123 3 .442, .471, .499, .523, .543, .619, 124 4 .485, .521, .548, .569, .589, .659, 125 5 .530, .569, .594, .613, .630, .693, 126 6 .578, .614, .639, .657, .671, .729, 127 7 .619, .655, .678, .692, .709, .761, 128 8 .659, .696, .713, .733, .746, .793, 129 9 .700, .734, .750, .765, .778, .819, 130 A .742, .771, .789, .800, .810, .850, 131 B .784, .809, .824, .838, .843, .875, 132 C .826, .843, .860, .865, .873, .900, 133 D .865, .885, .895, .900, .903, .921, 134 E .910, .921, .928, .931, .933, .938, 135 F .951, .962, .964, .966, .967, .968, 136 G 1.000, 1.000, 1.000, 1.000, 1.000, 1.000/ 137 DATA F1125B/4H6.1.,4H1.1-,4H39B / 138C 139C FIGURE 6.1.1.1-40 140C CTAB4(1,I) = FLAP DEFLECTION ANGLE (DEG) 141C CTAB4(2,I) = FLAP TO A1RFOIL CHORD RATIO 142C KPRIME(I) = EMPIRICAL CORRECTION FACTOR 143C 144 DATA X21126 145 1 / 0., 10., 12., 14., 16., 18., 20., 23., 27., 30., 35., 40., 146 2 50., 60. / 147 DATA X11126 148 1 / 0.10, 0.15, 0.20, 0.25, 0.30, 0.40, 0.50 / 149 DATA Y11126 150 1 / 1.000, 1.000, 0.994, 0.989, 0.970, 0.938, 0.900, 0.829, 151 2 0.755, 0.722, 0.672, 0.641, 0.596, 0.562, 152 3 1.000, 1.000, 0.994, 0.989, 0.970, 0.937, 0.890, 0.809, 153 4 0.737, 0.698, 0.650, 0.618, 0.569, 0.531, 154 5 1.000, 1.000, 0.994, 0.989, 0.968, 0.936, 0.870, 0.783, 155 6 0.710, 0.673, 0.630, 0.595, 0.542, 0.500, 156 7 1.000, 1.000, 0.994, 0.989, 0.965, 0.935, 0.850, 0.740, 157 8 0.677, 0.644, 0.600, 0.569, 0.518, 0.480, 158 9 1.000, 1.000, 0.994, 0.989, 0.963, 0.905, 0.800, 0.700, 159 A 0.643, 0.610, 0.570, 0.541, 0.496, 0.461, 160 B 1.000, 1.000, 0.993, 0.969, 0.924, 0.860, 0.750, 0.656, 161 C 0.606, 0.579, 0.540, 0.513, 0.471, 0.440, 162 D 1.000, 1.000, 0.981, 0.943, 0.880, 0.790, 0.695, 0.625, 163 E 0.571, 0.542, 0.512, 0.490, 0.450, 0.423/ 164 DATA F11126/4H6.1.,4H1.1-,4H40 / 165C 166C FIGURE 6.1.4.1-15 167C TIKB(1,I) ETA 168C TIKB(2,I) LAMBDA 169C TDKB(I) K(B) 170C 171 DATA X21419/ 172 1 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 173 2 0.8, 0.85, 0.9, 1.0/ 174 DATA X11419/ 175 1 0.0, 0.5, 1.0/ 176 DATA Y61419/ 177 1 0.0, 0.160, 0.305, 0.440, 0.560, 0.670, 0.772, 0.860, 178 2 .930, .96, .981, 1.0, 179 3 0.0, 0.140, 0.270, 0.400, 0.515, 0.630, 0.735, 0.830, 180 4 0.912, 0.947, 0.972, 1.0, 181 5 0., .125, .255, .37, .49, .6, .705, .8, 182 6 0.885, 0.921, 0.955, 1.0/ 183 DATA F61419/4H6.1.,4H4.1-,4H15 / 184C 185C ----FIGURE 6.1.1.1-41 186C 187 DATA Y11127/-.373,-.373,-.368,-.345,-.32,-.296,-.268,-.23,-.19, 188 1-.165,-.151,-.142,-.140 , -.45,-.446,-.43,-.412,-.39,-.364,-.328, 189 2-.29,-.24,-.207,-.188,-.172,-.167 , -.512,-.51,-.491,-.48,-.458, 190 3-.43,-.398,-.356,-.30,-.247,-.229,-.20,-.190,-.55,-.55,-.53 ,-.52, 191 4-.5,-.472,-.44,-.396,-.34,-.28,-.25,-.224,-.21, -.60,-.60,-.585, 192 5-.57,-.55,-.521,-.495,-.455,-.4,-.33,-.29,-.25,-.23 / 193 DATA X21127/0.,10.,20.,25.,30.,35.,40.,45.,50.,55.,60.,70.,80./ 194 DATA X11127/.15,.20,.25,.30,.40/,F11127/4H6.1.,4H1.1-,4H41 / 195C 196C FIGURE 6.1.4.1-14 197C TIADER(1,I) A(W) 198C TIADER(2,I) (ALFA(DELTA)) CL 199C TDADER(I) (ALFA(DELTA)) CL RATIO 200C 201 DATA X21418/ 202 1 0.0, 0.25, 0.50, 1.0, 2.0, 3.0, 4.0, 5.0, 203 2 6.0, 8.0/ 204 DATA X11418/ 205 1 -1.0, -.9, -.8, -.7, -.6, -.5, -.4, -.3, 206 2 -.2, -.1/ 207 DATA Y61418/10*1., 208 1 1.11, 1.075, 1.06, 1.049, 1.03, 1.019, 1.011, 1.01, 209 2 1.009, 1.007, 210 3 1.2, 1.15, 1.125, 1.095, 1.06, 1.04, 1.03, 1.02, 211 4 1.019, 1.015, 212 5 1.365, 1.245, 1.20, 1.145, 1.095, 1.065, 1.05, 1.04, 213 6 1.032, 1.025, 214 7 1.5, 1.35, 1.285, 1.205, 1.132, 1.095, 1.072, 1.06, 215 8 1.05, 1.039, 216 9 1.7, 1.5, 1.4, 1.29, 1.185, 1.135, 1.10, 1.085, 217 A 1.071, 1.055, 218 B 2.0, 1.63, 1.53, 1.39, 1.25, 1.185, 1.145, 1.118, 219 C 1.1, 1.079, 220 D 2.4, 1.90, 1.70, 1.52, 1.35, 1.26, 1.20, 1.165, 221 E 1.142, 1.11, 222 F 2.73, 2.43, 2.10, 1.73, 1.49, 1.375, 1.30, 1.24, 223 G 1.2, 1.16, 224 H 3.08, 2.85, 2.55, 2.18, 1.79, 1.58, 1.46, 1.385, 225 I 1.325, 1.245/ 226 DATA F61418/4H6.1.,4H4.1-,4H14 / 227C 228C ----FIGURE 6.1.1.3-12AD SPLIT AND PLAIN 229C 230 DATA X137AD/0.,2.,4.,5.,6.,8.,9.,10.,11.,12.,14.,15.,16.,17.,18., 231 1 19.,20./,F137AD/4H6.1.,4H1.3-,4H12AD/, 232 2Y137AD/2*1.,.979,.95,.92,.82,.80,.82,.85,.91,1.09,1.19,1.31,1.43, 233 3 1.51,1.57,1.60/ 234C 235C ----FIGURE 6.1.1.3-12B2 2-SLOT AND FOWLER 236C 237 DATA X137B2/0.,30./,Y137B2/0.,1.2/,F137B2/4H6.1.,4H1.3-,4H12B2/ 238C 239C ----FIGURE 6.1.1.3-12B1 SPLIT,PLAIN AND 1-SLOT 240C 241 DATA X137B1/0.,2.,4.,6.,8.,10.,12.,14.,16.,20.,24.,28.,30./, 242 1Y137B1/0.,.2,.34,.47,.57,.65,.72,.78,.83,.92,.99,1.04,1.06/ , 243 2F137B1/4H6.1.,4H1.3-,4H12B1/ 244C 245C ----FIGURE 6.1.1.3-12AB AVERAGE 2-SLOT,FOWLER 246C 247 DATA X137AB/0.,2.,5.,7.,9.,11.,13.,15.,16.,17.,18.,19./, 248 1Y137AB/1.0,1.0,1.04,1.09,1.17,1.29,1.45,1.64,1.73,1.77,1.80,1.82/, 249 2F137AB/4H6.1.,4H1.3-,4H12AB/ 250C 251C ----FIGURE 6.1.1.3-12AC NACA 2-SLOT,NACA 1-SLOT 252C 253 DATA X137AC/0.,5.,6.,8.,10.,12.,14.,16.,17.,18.,19.,20./, 254 1Y137AC/1.,1.,1.02,1.08,1.17,1.30,1.47,1.67,1.71,1.73,1.715,1.68/, 255 2F137AC/4H16.1,4H.1.3,4H-12A/ 256C 257C ----FIGURE 6.1.1.3-13A 258C ----A1=FOWLER , A2=1-SLOT , A3=2-SLOT , A4=SPLIT AND PLAIN 259C 260 DATA X1138A /0.,5.,10.,15.,20.,25.,30.,35.,40.,45.,50.,60./, 261 1 Y138A1/.4,.5,.61,.71,.79,.87,.94,.98,4*1.0/, 262 2 Y138A2/.18,.33,.47,.59,.70,.79,.87,.93,.97,3*1.0/, 263 3 Y138A3/.18,.32,.44,.56,.66,.76,.84,.90,.95,.99,2*1.0/, 264 4 Y138A4/0.,.17,.33,.46,.57,.67,.76,.83,.87,.92,.95,1.0/, 265 5 F138A1/4H6.1.,4H1.3-,4H13A1 /, 266 6 F138A2/4H6.1.,4H1.3-,4H13A2 /, 267 7 F138A3/4H6.1.,4H1.3-,4H13A3 /, 268 8 F138A4/4H6.1.,4H1.3-,4H13A4 / 269C 270C ----FIGURE 6.1.1.3-13B 271C ----B1=1-SLOT , 4-BAR , FOWLER , CIRCULAR ARC 272C ----B2= 2-SLOT, 4-BAR MOTION 273C 274 DATA X1138B/0.,.2,.3,.4,.45,.55,.6,.8,1.0/,F138B1/4H6.1.,4H1.3-, 275 1 4H13B1 /,F138B2/4H6.1.,4H1.3-,4H13B2 /, 276 2 Y138B1/0.,.26,.39,.50,.57,.66,.70,.87,1.0/, 277 3 Y138B2/0.,.11,.23,.40,.52,.66,.70,.87,1.0/ 278C 279C ----FIGURE 6.1.4.1-14A 280C 281 DATA X1418A/0.,.05,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/, 282 1 Y1418A/0.,-.26,-.40,-.55,-.66,-.75,-.82,-.88,-.93,-.97,-.99,-1./, 283 2 F1418A /4H6.1.,4H4.1-,4H13A / 284 DATA X61142/0.,.05,.1,.15,.2,.25,.3,.35,.4,.45,.5/ 285 DATA Y61142/0.,.03,.042,.052,.06,.068,.072,.078,.082,.086,.09/ 286 DATA X6143A/15.,20.,30.,40.,50.,60.,70.,.1,.25,.4,4*0./ 287 DATA Y6143A/2*.66,.645,.505,.5,.4,.35,2*.73,.7,.645,.55,.43, 288 1 .38,2*.78,.745,.68,.585,.475,.4/ 289 DATA X6143B/0.,15.,20.,30.,40.,60.,90.,0.,10.,20.,4*0./ 290 DATA Y6143B/9*1.,.945,.89,.885,.875,.865,2*1.,.9,.79,.75, 291 1 .725,.72/ 292 DATA F61142/4H6.1.,4H1.1-,4H42 / 293 DATA F6143A/4H6.1.,4H1.1-,4H43A / 294 DATA F6143B/4H6.1.,4H1.1-,4H43B / 295 DATA X11147/5.,10.,15.,20.,30.,40.,50.,60.,70.,75., 296 1 .1,.2,.3,.4,6*0./ 297 DATA Y11147/-.36,-.295,-.24,-.22,-.195,-.18,-.17,-.155,-.14, 298 1 -.13,-.435,-.37,-.32,-.29,-.26,-.245,-.23,-.215,-.195,-.185, 299 2 -.57,-.48,-.415,-.375,-.33,-.305,-.28,-.27,-.255,-.25,-.6, 300 3 -.545,-.48,-.44,-.378,-.35,-.325,-.3,-.275,-.26/ 301 DATA F11147/4H6.1.,4H1.1-,4H47 / 302 DATA X11150/0.,.05,.1,.2,.3,.4,.5/ 303 DATA Y11150/0.,-.0005,-.0016,-.0044,-.0082,-.0108,-.02/ 304 DATA F11150/4H6.1.,4H1.1-,4H50 / 305 DATA ITRANS/2,3,4,7,8/ 306C 307 IF(HTPL)GO TO 1000 308C 309C -----HERE FOR H.L. OR CONTROL DEVICE ON WING 310C 311 CBAREX=A(16) 312 BO2=WINGIN(3) 313 TANTE=A(80) 314 TANLE=A(62) 315 CR=WINGIN(6) 316 TAPRI=A(25) 317 BSTRO2=WINGIN(2) 318 TANTEO=A(104) 319 TANLEO=A(86) 320 CB=WINGIN(5) 321 TAPRO=A(28) 322 SR=SREF 323 AW=A(120) 324 BTHEO=WINGIN(4) 325 TOVC=WINGIN(16) 326 TOVCO=WINGIN(66) 327 COSC4=A(67) 328 CLASEC=WINGIN(M+20) 329 CLAW=WING(101) 330 IF(TRANSN)CLAW=TRA(70) 331 IF(TRANSN)CLASEC=WINGIN(69)/.8 332 GO TO 1010 333 1000 CBAREX=AHT(16) 334 BO2=HTIN(3) 335 TANTE=AHT(80) 336 TANLE=AHT(62) 337 CR=HTIN(6) 338 TAPRI=AHT(25) 339 BSTRO2=HTIN(2) 340 TANTEO=AHT(104) 341 TANLEO=AHT(86) 342 CB=HTIN(5) 343 TAPRO=AHT(28) 344 SR=SREF 345 AW=AHT(120) 346 BTHEO=HTIN(4) 347 TOVC=HTIN(16) 348 TOVCO=HTIN(66) 349 COSC4=AHT(67) 350 CLASEC=HTIN(M+20) 351 CLAW=HT(101) 352 IF(TRANSN)CLAW=TRAH(70) 353 IF(TRANSN)CLASEC=HTIN(69)/.8 354 1010 CONTINUE 355C 356C*****CHECK FOR TRANSALATING DEVICE,IF SO CALCULATE SWF 357C 358 TRANSL=.FALSE. 359 IFTYPE=FTYPE+0.5 360 NDELTA=RDELTA+0.5 361 RF=FLC(M+42) 362 MACH=FLC(M+2) 363 IF(TRANSN)MACH=.6 364C 365C-----IF DELTA EQUALS ZERO,SET TO 0.01 DEG 366C 367 DO 1020 J=1,NDELTA 368 1020 IF(DELTA(J).EQ.0.0)DELTA(J)=0.01 369 DO 1030 K=1,5 370 1030 IF(IFTYPE.EQ.ITRANS(K))TRANSL=.TRUE. 371C 372C*****GEOMETRIC PARAMETERS PER FLAP SEGMENT 373C 374 DELN4=0.25*(BOF-BIF)/BTHEO 375 ETA(1)=BIF/BTHEO 376 VAR(1)=ETA(1) 377 CF(1)=CFI 378 ARG1=(CFI-CFO)/(4.*DELN4) 379 ARG2=(TANTE-TANLE)*BTHEO 380 ARG3=(BOF-BIF)*CR /4.0 381 ARG4=CR 382 ARG5=TAPRI 383 TC=TOVC 384 IF(BIF.LT.(BTHEO-BSTRO2))GO TO 1040 385 ARG2=(TANTEO-TANLEO)*BO2 386 ARG3=(BOF-BIF)*CB /4.0 387 ARG4=CB 388 ARG5=TAPRO 389 TC=TOVCO 390 1040 CONTINUE 391C 392C FIGURE 6.1.4.1-15 (KB) 393C 394 CALL TLINEX(X11419,X21419,Y61419,3,12,ARG5,ETA(1),RKB(1), 395 1 0,0,0,0,F61419,3,ROUTID) 396 CHRD(1)=ARG4+ETA(1)*ARG2 397C 398C FIGURE 6.1.4.1-14 (ALPHA)DELTA 399C 400 CFOC(1)=CF(1)/CHRD(1) 401 VAR(1)=CFOC(1) 402 LGH(1)=12 403 CALL INTERX(1,X1418A,VAR,LGH,Y1418A,ALPHAD(1),12,12, 404 1 0,0,0,0,0,0,0,0,F1418A,3,ROUTID) 405 DO 1050 K=2,5 406 N=K-1 407 ETA(K)=ETA(N)+DELN4 408 CF(K)=CFI-ARG1*(ETA(K)-ETA(1)) 409C 410C FIGURE 6.1.4.1-15 (KB) 411C 412 CALL TLINEX(X11419,X21419,Y61419,3,12,ARG5,ETA(K),RKB(K), 413 1 0,0,0,0,F61419,3,ROUTID) 414 DKB(N)=RKB(K)-RKB(N) 415 CHRD(K)=ARG4+ETA(K)*ARG2 416C 417C FIGURE 6.1.4.1-14 (ALPHA)DELTA 418C 419 CFOC(K)=CF(K)/CHRD(K) 420 VAR(1)=CFOC(K) 421 LGH(1)=12 422 CALL INTERX(1,X1418A,VAR,LGH,Y1418A,ALPHAD(K),12,12, 423 1 0,0,0,0,0,0,0,0,F1418A,3,ROUTID) 424 ALDAVG(N)=0.50*(ALPHAD(K)+ALPHAD(N)) 425 SWF(N)=ARG3*(2.-(1.-ARG5)*(ETA(N)+ETA(K))) 426 1050 CONTINUE 427C 428C*****CALCULATE SECTION LIFT INCREMENT FOR EACH SPANWISE SEGMENT 429C 430 EXPDCL=.FALSE. 431 NN=0 432C 433C FIGURE 4.1.1.2-8A 434C 435 ARG1=ALOG10(RF*CBAREX) 436 IF(IFTYPE.EQ.1)CALL TLINEX(X1128A,X2128A,Y1128A,3,11,ARG1,TANPHE, 437 1 CLOCLT,1,0,0,1,F1128A,3,ROUTID) 438 DO 1200 I=1,NDELTA 439 ARG8=DELTA(I)*CLASEC 440 IF(SDCL(I).NE.UNUSED)EXPDCL=.TRUE. 441 IF(EXPDCL.AND.(.NOT.TRANSL))GO TO 1190 442 CP(1)=CPI(I) 443 DO 1170 K=1,5 444 IF(K.EQ.1)GO TO 1060 445 IF(.NOT.TRANSL)GO TO 1060 446 ARG1= (CPI(I)-CPO(I))/(4.*DELN4) 447 CP(K)=CPI(I)-ARG1*(ETA(K)-ETA(1)) 448 1060 CONTINUE 449 IF(EXPDCL)GO TO 1150 450 ARGZ=ABS(DELTA(I)) 451 GO TO (1070,1090,1090,1150,1100,1110,1120,1130),IFTYPE 452 1070 CONTINUE 453C 454C*****PLAIN TRAILING-EDGE FLAPS 455C 456 IF(I.GT.1)GO TO 1080 457C 458C FIGURE 6.1.1.1-39B 459C 460 CALL TLINEX(X1125B,X2125B,Y1125B,16,6,CLOCLT,CFOC(K), 461 1 CLDOCT(K),0,0,0,0,F1125B,3,ROUTID) 462C 463C FIGURE 6.1.1.1-39A 464C 465 CALL TLINEX(X1125A,X2125A,Y1125A,8,8,TC,CFOC(K),CLDTHY(K), 466 1 0,0,0,0,F1125A,3,ROUTID) 467C 468C FIGURE 6.1.1.1-40 469C 470 1080 CALL TLINEX(X11126,X21126,Y11126,7,14,CFOC(K),ARGZ,KFPRM, 471 1 0,0,0,0,F11126,3,ROUTID) 472 DELCL(K)=DELTA(I)*CLDOCT(K)*CLDTHY(K)*KFPRM/RAD 473 GO TO 1140 474 1090 CONTINUE 475C 476C*****SINGLE-SLOTTED AND FOWLER TRAILING-EDGE FLAPS 477C FIGURE 6.1.1.1.-41 478C 479 CALL TLINEX(X11127,X21127,Y11127,5,13,CFOC(K),ARGZ, 480 1 ALPHAD(K),0,0,0,0,F11127,3,ROUTID) 481 DELCL(K)=-CLASEC*ALPHAD(K)*DELTA(I) 482 GO TO 1140 483 1100 CONTINUE 484C 485C*****SPLIT TRAILING-EDGE FLAPS 486C 487 VAR(1)=DELTA(I) 488 VAR(2)=CFOC(K) 489 LGH(1)=10 490 LGH(2)=4 491 CALL INTERX(2,X11147,VAR,LGH,Y11147,ALFAD,10,40, 492 1 1,1,0,0,1,1,0,0,F11147,3,ROUTID) 493 DELCL(K)=-CLASEC*ALFAD*DELTA(I) 494 GO TO 1140 495 1110 CONTINUE 496C 497C*****LEADING-EDGE FLAPS 498C 499 VAR(1)=CFOC(K) 500 LGH(1)=7 501 CALL INTERX(1,X11150,VAR,LGH,Y11150,CLDK,7,7, 502 1 0,0,0,0,1,0,0,0,F11150,3,ROUTID) 503 DELCL(K)=CLDK*DELTA(I) 504 GO TO 1140 505 1120 CONTINUE 506C 507C*****LEADING-EDGE SLATS 508C 509 VAR(1)=CFOC(K) 510 LGH(1)=7 511 CALL INTERX(1,X11150,VAR,LGH,Y11150,CLDK,7,7, 512 1 0,0,0,0,1,0,0,0,F11150,3,ROUTID) 513 DELCL(K)=CLDK*DELTA(I)*CP(K)/CHRD(K) 514 GO TO 1140 515 1130 CONTINUE 516C 517C*****KRUEGER L.E. FLAPS 518C 519 VAR(1)=CFOC(K) 520 LGH(1)=7 521 CALL INTERX(1,X11150,VAR,LGH,Y11150,CLDK,7,7, 522 1 0,0,0,0,1,0,0,0,F11150,3,ROUTID) 523 DELCL(K)=CLDK*DELTA(I)*CP(K)/CHRD(K) 524 1140 CONTINUE 525C 526C*****COMPUTE AVERAGE SECTION DELTA CL OVER SPANWISE SEGMENT 527C 528 N=K-1 529 IF(N.EQ.0) GO TO 1170 530 NN=NN+1 531 DELCLA(NN)=(DELCL(K)+DELCL(N))/2. 532 ALDAG(NN)=-DELCLA(NN)/ARG8 533C 534C*****COMPUTE AVERAGE CP/CK OVER SPANWISE SEGMENT 535C 536 IF(IFTYPE .EQ. 4) GO TO 1160 537 IF(.NOT.TRANSL)GO TO 1160 538 1150 CONTINUE 539 CPOCF(N)=(CP(K)/CHRD(K)+CP(N)/CHRD(N))/2. 540 CFACTR(N)=(CPOCF(N)-1.)*SWF(N)/SR 541 IF(EXPDCL)GO TO 1170 542 IF(IFTYPE .NE. 4) GO TO 1160 543C 544C*****HERE FOR SECTION DELTA-CL FOR DOUBLE SLOTTED FLAPS 545C 546 AARG1=(CF2I-CF2O)/(4.*DELN4) 547 AARG2=(CAPI(I)-CAPO(I))/(4.*DELN4) 548 CAPR=CAPI(I)-AARG2*(ETA(K)-ETA(1)) 549 CF2(K)=CF2I-AARG1*(ETA(K)-ETA(1)) 550 CF2OC(K)=CF2(K)/CHRD(K) 551 PHI1=DELTA(I)+ATAN(TANPHE)*RAD 552 PHI2=PHI1+DF2(I) 553 VAR(1)=PHI1 554 IF(CFOC(K) .LE. .10) VAR(2)=.10 555 IF(CFOC(K) .GE. .40) VAR(2)=.40 556 IF(CFOC(K) .GT. .10 .AND. CFOC(K) .LT. .40) VAR(2)=CFOC(K) 557 LGH(1)=7 558 LGH(2)=3 559 CALL INTERX(2,X6143A,VAR,LGH,Y6143A,ATEA1,7,21, 560 1 1,1,0,0,1,1,0,0,F6143A,3,ROUTID) 561 LGH(1)=11 562 VAR(1)=CFOC(K) 563 CALL INTERX(1,X61142,VAR,LGH,Y61142,CLDF1,11,11, 564 1 1,0,0,0,1,0,0,0,F61142,3,ROUTID) 565 VAR(1)=PHI2 566 LGH(1)=7 567 IF(CF2OC(K) .LE. .10) VAR(2)=.10 568 IF(CF2OC(K) .GE. .40) VAR(2)=.40 569 IF(CF2OC(K) .GT. .10 .AND. CF2OC(K) .LT. .4)VAR(2)=CF2OC(K) 570 CALL INTERX(2,X6143A,VAR,LGH,Y6143A,ATEA2,7,21, 571 1 1,1,0,0,1,1,0,0,F6143A,3,ROUTID) 572 LGH(1)=11 573 VAR(1)=CFOC(K) 574 CALL INTERX(1,X61142,VAR,LGH,Y61142,CLDF2,11,11, 575 1 1,0,0,0,1,0,0,0,F61142,3,ROUTID) 576 IF((CFOC(K)/CF2OC(K)) .LE. 0.60) DELCL(N)=ATEA1*CLDF1*DELTA 577 1 (I)*(1.+CFOC(K))+ATEA2*CLDF2*(DELTA(I)+DF2(I))*CP(K)/CHRD(K) 578 IF((CFOC(K)/CF2OC(K)) .LE. 0.60) GO TO 1140 579 VAR(1)=DF2(I) 580 VAR(2)=DELTA(I) 581 LGH(1)=7 582 LGH(2)=3 583 CALL INTERX(2,X6143B,VAR,LGH,Y6143B,ATEAT,7,21, 584 1 1,1,0,0,1,1,0,0,F6143B,3,ROUTID) 585 DELCL(N)=ATEA1*CLDF1*DELTA(I)*CAPR/CHRD(K)+ATEA2*ATEAT 586 1 *CLDF2*DF2(I)*(1.+(CP(K)-CAPR)/CHRD(K)) 587 GO TO 1140 588 1160 CONTINUE 589C 590C*****SPANWISE SEGMENT DELTA CL 591C FIGURE 6.1.4.1-8 592C 593 CALL TLINEX(X11418,X21418,Y61418,10,10,ALDAVG(N),AW, 594 1 ADCADS(N),0,0,0,1,F61418,3,ROUTID) 595 DCLK(N)=DELCLA(NN)*ADCADS(N)*DKB(N)*CLAW/CLASEC 596 1170 CONTINUE 597 IF(EXPDCL)GO TO 1180 598 DCL(I)=DCLK(1)+DCLK(2)+DCLK(3)+DCLK(4) 599 IF(.NOT.TRANSL)GO TO 1200 600 1180 CONTINUE 601 CFACT(I)=(CFACTR(1)+CFACTR(2)+CFACTR(3)+CFACTR(4))/4.0 602 CLAD(I)=CFACT(I)*CLAW+CLAW 603 IF(.NOT.EXPDCL)GO TO 1200 604 1190 CONTINUE 605C 606C*****HERE IF SECTION DELTA CL WAS INPUT 607C 608 VAR(1)=SDCL(I)/(CLASEC*DELTA(I)) 609 VAR(2)=AW 610C 611C FIGURE 6.1.4.1-8 612C 613 CALL TLINEX(X11418,X21418,Y61418,10,10,VAR(1),VAR(2),ADCAD(I), 614 1 0,0,0,1,F61418,3,ROUTID) 615 DCL(I)=SDCL(I)*ADCAD(I)*(RKB(5)-RKB(1))*CLAW/CLASEC 616 1200 CONTINUE 617 IF(IFTYPE.GE.6) GO TO 1380 618 DO 1370 I=1,NDELTA 619 IF(I.GT.1)GO TO 1270 620C 621C*****CALCULATE SECTION DELTA (CL)MAX 622C 623 TROFS=CHRD(5)/CHRD(1) 624 CBARFS=2./3.*CHRD(1)*(1.+TROFS*(1.+TROFS))/(1.+TROFS) 625 ETAFS=(1.+2.*TROFS)/(3.*(1.+TROFS)) 626 CFOCA=(CFI-ETAFS*(CFI-CFO))/CBARFS 627 VAR(1)=TC*100. 628 GO TO(1220,1210,1230,1230,1220,1380,1380,1380),IFTYPE 629 1210 CONTINUE 630C 631C FIGURE 6.1.1.3-12AC (SINGLE SLOT) (DCL)BASE 632C 633 LGH(1)=12 634 CALL INTERX(1,X137AC,VAR,LGH,Y137AC,DCLMAB ,12,12, 635 1 0,0,0,0,0,0,0,0,F137AC,3,ROUTID) 636 GO TO 1240 637 1220 CONTINUE 638C 639C FIGURE 6.1.1.3-12AD (SPLIT AND PLAIN) (DCL)BASE 640C 641 LGH(1)=17 642 CALL INTERX(1,X137AD,VAR,LGH,Y137AD,DCLMAB ,17,17, 643 1 0,0,0,0,0,0,0,0,F137AD,3,ROUTID) 644 GO TO 1240 645 1230 CONTINUE 646C 647C FIGURE 6.1.1.3-12AB (2-SLOT AND FOWLER) (DCL)BASE 648C 649 LGH(1)=12 650 CALL INTERX(1,X137AB,VAR,LGH,Y137AB,DCLMAB ,12,12, 651 1 0,0,0,0,0,0,0,0,F137AB,3,ROUTID) 652 1240 CONTINUE 653 VAR(1)=CFOCA *100. 654 GO TO(1250,1250,1260,1260,1250,1380,1380,1380),IFTYPE 655 1250 CONTINUE 656C 657C FIGURE 6.1.1.3-12B1 (SPLIT,PLAIN,AND 1-SLOT) K1 658C 659 LGH(1)=13 660 CALL INTERX(1,X137B1,VAR,LGH,Y137B1,RK1,13,13, 661 1 0,0,0,0,1,0,0,0,F137B1,3,ROUTID) 662 GO TO 1270 663 1260 CONTINUE 664C 665C FIGURE 6.1.1.3-12B2 (2-SLOT AND FOWLER) K1 666C 667 LGH(1)=2 668 CALL INTERX(1,X137B2,VAR,LGH,Y137B2,RK1,2,2, 669 1 0,0,0,0,1,0,0,0,F137B2,3,ROUTID) 670 1270 CONTINUE 671C 672C*****SET REFERENCE FLAP DEFLECTIONS FROM FIG.6.1.1.3-13A 673C 674 DREFF=40. 675 DREF1S=45. 676 DREF2S=50. 677 DREFSP=60. 678 VAR(1)=ABS(DELTA(I)) 679 LGH(1)=12 680 GO TO(1280,1290,1300,1310,1280,1380,1320,1380),IFTYPE 681 1280 CONTINUE 682C 683C FIGURE 6.1.1.3-13A (SPLIT AND PLAIN) K2 684C 685 CALL INTERX(1,X1138A,VAR,LGH,Y138A4,RK2(I),12,12, 686 1 0,0,0,0,0,0,0,0,F138A4,3,ROUTID) 687 GO TO 1330 688 1290 CONTINUE 689C 690C FIGURE 6.1.1.3-13A (1-SLOT) K2 691C 692 CALL INTERX(1,X1138A,VAR,LGH,Y138A2,RK2(I),12,12, 693 1 0,0,0,0,0,0,0,0,F138A2,3,ROUTID) 694 GO TO 1330 695 1300 CONTINUE 696C 697C FIGURE 6.1.1.3-13A (FOWLER) K2 698C 699 CALL INTERX(1,X1138A,VAR,LGH,Y138A1,RK2(I),12,12, 700 1 0,0,0,0,0,0,0,0,F138A1,3,ROUTID) 701 GO TO 1330 702 1310 CONTINUE 703C 704C FIGURE 6.1.1.3-13A (2-SLOT) K2 705C 706 CALL INTERX(1,X1138A,VAR,LGH,Y138A3,RK2(I),12,12, 707 1 0,0,0,0,0,0,0,0,F138A3,3,ROUTID) 708 GO TO 1330 709 1320 CONTINUE 710C 711C*****DELTA CLMAX FOR L.E. SLATS 712C 713 DCLMAX(I)=1.28*CFOC(K)/.18*COSC4**2*((BOF-BIF)/BTHEO)**2 714 GO TO 1370 715 1330 CONTINUE 716 IF(IFTYPE.EQ.1.OR.IFTYPE.EQ.5) GO TO 1350 717 IF(IFTYPE.EQ.4) GO TO 1340 718C 719C*****HERE FOR 1-SLOT OR FOWLER 720C 721 VAR(1)=ABS(DELTA(I))/DREFF 722 IF(IFTYPE.EQ.2)VAR(1)=ABS(DELTA(I))/DREF1S 723C 724C FIGURE 6.1.1.3-13B1 (1-SLOT AND FOWLER ) K3 725C 726 LGH(1)=9 727 CALL INTERX(1,X1138B,VAR,LGH,Y138B1,RK3,9,9, 728 1 0,0,0,0,0,0,0,0,F138B1,3,ROUTID) 729 GO TO 1360 730 1340 CONTINUE 731C 732C*****HERE FOR 2-SLOT 733C 734 VAR(1)=ABS(DELTA(I))/DREF2S 735C 736C FIGURE 6.1.1.3-13B2 (2-SLOT) K3 737C 738 LGH(1)=9 739 CALL INTERX(1,X1138B,VAR,LGH,Y138B2,RK3,9,9, 740 1 0,0,0,0,0,0,0,0,F138B2,3,ROUTID) 741 GO TO 1360 742 1350 RK3=1. 743 1360 CONTINUE 744 DSCLMX(I)=RK1*RK2(I)*RK3*DCLMAB 745C 746C*****INCREMENT IN MAXIMUM LIFT COEFFICIENT DUE TO TE FLAP DEFLECTION 747C 748 SWFT=SWF(1)+SWF(2)+SWF(3)+SWF(4) 749 KSWEEP=(1.-0.08*COSC4**2)*COSC4**0.75 750 DCLMAX(I)=DSCLMX(I)*SWFT*KSWEEP/SR 751 1370 CONTINUE 752 1380 CONTINUE 753 RETURN 754 END 755