1 SUBROUTINE WTLIFT(A,B,AIN,AOUT) 2C 3C *** CALCULATES WING OR TAIL LIFT CHARACTERISTICS 4C 5 COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF 6 COMMON /OVERLY/ NLOG,NMACH,IM,NALPHA,IG 7 COMMON /CONSNT/ PI, DEG, UNUSED, RAD 8 DIMENSION ROUTID(2),Q3417B(3),Q13242(3),Q3418A(3),Q3418B(3) 9 DIMENSION Q13415(3),Q3414A(3),Q3414B(3) 10 DIMENSION AMACH(5) , DYA(6) , SALE4(4), DCLTB(5,6,4) , DCAR(120) 11 DIMENSION A(195),B(49),AOUT(101),AIN(101) 12 DIMENSION DACLMX(13,4),DACLL(52) 13 DIMENSION SALE(13) , DELTAY(7) ,CLOCL(13,7),CLL(91) , DY(4) 14 DIMENSION TR(11),C2(11),C(6) 15 DIMENSION C1ABCS(10),ACLMX(10),C17(6) 16 DIMENSION DMN(3), C1TAB(20), DACL(60) 17 DIMENSION ACLE(10),DACLO(100) , C1TABO(10) 18 DIMENSION BA(12),CLOVCL(12) , CBA(6) 19 DIMENSION WTYPE(4) 20C 21 EQUIVALENCE (DCLTB(1,1,1),DCAR(1)) 22 EQUIVALENCE (CLOCL(1,1),CLL(1)) ,(DACLMX(1,1),DACLL(1)) 23C 24 LOGICAL ITYPTT 25 LOGICAL LOARAT 26C 27 DATA ROUTID /4HWTLI,4HFT / 28 DATA Q3417B /4H4.1.,4H3.4-,4H24B /, 29 1 Q13242 /4H4.1.,4H3.2-,4H52 /, 30 2 Q3418A /4H4.1.,4H3.4-,4H25A /, 31 3 Q3418B /4H4.1.,4H3.4-,4H25B /, 32 4 Q13415 /4H4.1.,4H3.4-,4H22 /, 33 5 Q3414A /4H4.1.,4H3.4-,4H21A /, 34 6 Q3414B /4H4.1.,4H3.4-,4H21B / 35C 36C ----4.1.3.4-22 37C 38 DATA AMACH/.2,.3,.4,.5,.6/ 39 DATA DYA /2.,2.25,2.5,3.,4.,4.5/ 40 DATA SALE4/0.,20.,40.,60./ 41 DATA DCAR/0.,2*-.02,2*0., 0.,-.13,-.19,-.2,-.19 , 0.,-.185,-.32, 42 1-.4,-.445 , 0.,-.21,-.36,-.45,-.5 , 0.,-.24,-.42,-.545,-.64 , 43 2-0.,-.24,-.455,-.605,-.72 , 0.,-.04,-.045,-.02,0. , 0.,-.095,-.165 44 3,-.22,-.25 , 0.,-.105,-.19,-.26,-.29 , 0.,-.125,-.225,-.3,-.345 , 45 4 0.,-.15,-.265,-.37,-.45 , 0.,-.15,-.29,-.41,-.51 , 0.,-.04,-.07, 46 5-.095,-.1 , 0.,-.043,-.080,-.108,-.123 , 0.,-.045,-.085,-.12,-.145 47 6 , 0.,-.05,-.095,-.14,-.19 , 0.,-.07,-.125,-.19,-.26 , 0.,-.07, 48 7-.14,-.215,-.3,4*0.,-.02,4*0.,-.022,4*0.,-.03,3*0.,-.02,-.07,2*0., 49 8-.02,-.055,-.085,0.,-.04,-.085,-.14,-.2/ 50C 51C ----4.1.3.4-21 A 52C 53 DATA SALE/0.,5.,10.,15.,20.,25.,30.,35.,40.,45.,50.,55.,60./ 54 DATA DELTAY/1.4,1.6,1.8,2.0,2.2,2.4,2.5/ 55 DATA CLL /.9,.915,.93,.95,.975,1.,1.03,1.07,1.1,1.15,1.19,1.25,1.3 56 1 , .9,.91,.92,.94,.96,.98,1.,1.02,1.05,1.08,1.11,1.15,1.19 , .9, 57 2 .91,.92,.93,.94,.95,.96,.965,.975,.99,1.,1.01,1.03 , 6*.9, 2*.89, 58 3.885,2*.88,.875,.87 , .9,.895,.89,.88,.87,.86,.85,.835,.82,.8,.78, 59 4.755,.73 , .9,.89,.88,.87,.85,.83,.81,.79,.76,.725,.695,.65,.59 , 60 5 .9,.89,.87,.86,.84,.82,.795,.77,.73,.7,.65,.59,.52 / 61C 62C ----4.1.3.4-21 B 63C 64 DATA DY/1.2,2.,3.,4./ 65 DATA DACLL/1.75,1.9,2.2,2.7,3.4,4.15,5.1,6.1,7.3,8.7,10.15,11.75, 66 113.3 , .1,.5,1.05,1.65,2.3,3.1,3.9,4.7,5.7,6.7,7.7,8.75,9.8, 67 21.2,1.4,1.7,2.,2.4,2.85,3.35,3.7,4.25,4.7,5.3,5.9,6.65, 68 32.2,2.1,2*2.,2.1,2.15,2.3,2.4,2.55,2.7,2.9,3.05,3.3/ 69C 70C ----4.1.3.4-24(B) 71C 72 DATA IN/0/ 73 DATA TR/0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1./ 74 DATA C2/0.,.21,.5,.9,1.08,1.05,1.,.94,.9,.86,.85/ 75C 76C ----4.1.3.4-25(A) 77C 78 DATA C1ABCS/0.,.4,.8,1.2,1.6,2.0,2.4,2.8,3.2,3.6/ 79 DATA ACLMX/3*35.,32.,28.,25.,23.2,22.,21.5,21./ 80C 81C ----4.1.3.4-25(B) 82C 83 DATA I17/0/ 84 DATA DMN/ .2,.4,.6/ 85 DATA C1TAB/4.5,5.,5.5,6.,6.5,7.,7.5,8.,8.5,9.,9.5,10.,10.5,11., 86 111.5,12.,12.5,13.,13.5,14./ 87 DATA DACL /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. 88 1,11.,12.,13.,14.,15. , 0.,.2,.4,.7,1.2,1.7,2.4,3.,3.7,4.6,5.2,6., 89 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 90 3,2.,2.5,3.,3.6,4.3,4.9,5.5,6.2,7.,7.6,8.5/ 91 DATA ACLE /0.,2.,3.,4.,5.,6.,7.,8.,9.,30./ 92 DATA DACLO / 10.,8.5,6.9,5.5,4. 93 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, 94 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 95 3,-2.,-.8,0. , 3.,.7,-1.4,-3.3,-4.3,-4.3,-3.1,-2.,-.8,0. , .3,-2.5, 96 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, 97 5-3.1,-2.,-0.8,0. , -3.3,-6.6,-7.2,-6.3,-5.3,-4.3,-3.1,-2.,-.8,0., 98 6 -4.2,-7.,-7.2,-6.3,-5.3,-4.3,-3.1,-2.,-.8,0. , -8.5,-7.9,-7.2, 99 7 -6.3,-5.3,-4.3,-3.1,-2.,-.8,0. / 100 DATA C1TABO /.0,.5,1.,1.5,2.,2.5,3.,3.5,4.,4.5/ 101C 102C ----4.1.3.2-52 103C 104 DATA BA/1.15,1.4,2.,2.2,3.,3.6,4.,5.,6.,7.,8.,9./ 105 DATA CLOVCL /1.25,1.2,1.12,1.1,1.04,1.,.985,.96,.95,3*.94/ 106 DATA IN4132 /0/ 107 DATA WTYPE /4HSTRA ,4HDOUB ,4HCRAN ,4HCURV / 108C 109C ----INITIALIZATION ENTRY 110C 111 ITYPTT = .TRUE. 112 LOARAT = .FALSE. 113 IF(A(7).LT.A(125)) LOARAT = .TRUE. 114 TYPE = AIN(15) 115C 116C ----COMPUTE CLA0 + CLAI FOR OTHER THAN STRAIGHT TAPERED 117C 118 IF(TYPE.EQ.WTYPE(1)) GO TO 1000 119C 120C ----CLA0 121C 122 TEMP1 = 2.*PI*A(168)*DEG 123 TEMP2 =(A(168)*DEG*2.*PI/A(131))**2 124 TEMP3 = A(98)**2 125 TEMP4=B(2)**2 126 A(172)=TEMP1/(2.+SQRT(TEMP2*(1.+TEMP3/TEMP4)+4.)) 127C 128C ----CLAI 129C 130 TEMP1 = 2.*PI*A(5)*DEG 131 TEMP2 = (A(5)*DEG*2.*PI/A(131))**2 132 TEMP3 = A(74)**2 133 A(171)=TEMP1/(2.+SQRT(TEMP2*(1.+TEMP3/TEMP4)+4.)) 134 1000 CONTINUE 135C 136C ----CURVED WING TEST 137C 138 IF(TYPE.EQ.WTYPE(4)) GO TO 1070 139 ITYPTT = .FALSE. 140C 141C ----4.1.3.4-24(B) 142C 143 CALL TBFUNX(A(27),A(159),DCDT,11,TR,C2,C,IN,MI,NG, 144 1 0,0,Q3417B,3,ROUTID) 145 A(160) = (A(159)+1.)*A(38)*A(7) 146C 147C ----COMPUTE BETA*S 148C 149 TEMP1 = 2.*PI*A(7)*DEG 150 TEMP2=A(7)*DEG*2.*PI/A(131) 151 TEMP2 = TEMP2**2 152 TEMP3 = A(50)**2 153 CALL TLINEX(DELTAY, SALE,CLL,7,13,AIN(17),A(34),A(145), 154 1 -1,0,-1,2,Q3414A,3,ROUTID) 155 CALL TLINEX(DY,SALE,DACLL,4,13,AIN(17),A(34),A(144), 156 1 -1,0,-1,2,Q3414B,3,ROUTID) 157 1010 IF(ITYPTT) GO TO 1070 158 TEMP4 = B(2)**2 159 AOUT(101)=(TEMP1/(2.+SQRT(TEMP2*(1.+TEMP3/TEMP4)+4.)))*A(3)/ 160 1 SREF 161 IF(TYPE.NE.WTYPE(3)) GO TO 1020 162C 163C ----HERE FOR CRANKED HIGH ASPECT RATIO WING-FINISH CLA COMPUTATION 164C ----FIGURE 4.1.3.2-52 165C 166 BAARG = A(7)* B(2) 167 CALL TBFUNX(BAARG,CLRAT,DYDX,12,BA,CLOVCL,CBA ,IN4132,M4132,N4132, 168 1 2,0,Q13242,3,ROUTID) 169 AOUT(101) = CLRAT * AOUT(101) 170 1020 CONTINUE 171 IF(LOARAT) GO TO 1030 172 CALL TLIN3X(DYA,AMACH,SALE4,DCAR,6,5,4,AIN(17),B(1),A( 34),TEMP5 173 1 ,0,-1,0,2,2,2,Q13415,3,ROUTID) 174 A(146)=A(145)*A(132) 175 B(44)=(A(145)*A(132)+TEMP5)*A(3)/SREF 176 B(43) = B(44)/AOUT(101)+B(49)+A(144) 177 GO TO 1060 178C 179C ----HERE FOR LOW ASPECT RATIO 180C 181 1030 C1P1AC = A(7)*A(124)/B(2) 182 CALL CLMXBS(C1P1AC,B(44),A,B,AIN) 183C 184C ----4.1.3.4-25(A) 185C 186 CALL TBFUNX(C1P1AC,ACLBAS,DYDX,10,C1ABCS,ACLMX,C17,I17,MI,NG, 187 1 0,0,Q3418A,3,ROUTID) 188 IF(A(160).LE.4.5) GO TO 1040 189 CALL TLINEX(DMN,C1TAB,DACL,3,20,B(1),A(160),DACMAX, 190 1 -1,0,2,0,Q3418B,3,ROUTID) 191 GO TO 1050 192 1040 TMP=A(7)*A(37)*(1.+4.*A(27)**2) 193 CALL TLINEX(ACLE,C1TABO,DACLO,10,10,TMP,A(160),DACMAX, 194 1 0,0,2,0,Q3418B,3,ROUTID) 195 1050 B(43)=ACLBAS+DACMAX 196 1060 CONTINUE 197 RETURN 198 1070 WRITE(6,1080) TYPE,LOARAT 199 1080 FORMAT(33H NO CLALPHA COMPUTATION- WING IS ,A4,18H LOW ASPECT RAT 200 1IO= ,L1) 201 RETURN 202 END 203