1 SUBROUTINE VTDRAG(I) 2C 3C *** CALCULATES VERTICAL DRAG COEFFICIENT 4C 5 DIMENSION ROUTID(2) 6 DIMENSION Q15127(3),Q5128B(3),Q128BD(5) 7 COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF 8 COMMON /FLGTCD/ FLC(93) 9 COMMON /WHAERO/ C2(51),D2(55),CHT2(51),DHT2(55),DVT(55) 10 COMMON /VTDATA/ AVT(195) 11 COMMON /VTI/ VT(100) 12 DIMENSION D(55),A(195),B(93) 13 EQUIVALENCE (A(1),AVT(1)),(D(1),DVT(1)),(B(1),FLC(1)) 14 DIMENSION WTYPE(4) 15 DIMENSION X228BD(9),X128BD(4),Y28BD(36) 16 DIMENSION X228B(11), X128B(4), Y28B(44) 17 DIMENSION X27M(4),X27I(4),C27(6) 18 DATA ROUTID/4HVTDR,4HAG / 19 DATA Q15127/4H4.1.,4H5.1-,4H27 /,Q5128B/4H4.1.,4H5.1-,4H28-B/, 20 1 Q128BD/4H4.1.,4H5.1-,4H28-B,4H(DAS,4HHED)/ 21 DATA X27M /0.0,1.0,2.0,3.0/ 22 DATA X27I /1.57780,1.67221,1.98509,2.28874/ 23 DATA I27 /0/ 24C 25C--------- FIGURE 4.1.5.1-28B ---- 26C---- X228B=COS(SWEEP ANGLE(T/C)MAX) X128B= MACH NUMBER 27C---- Y28B= (R)LS 28C 29 DATA WTYPE / 4HSTRA ,4HDOUB ,4HCRAN ,4HCURV / 30 DATA X228B /0.5,.55,.60,.65,.70,.75,.80,.85,.90,.95,1.0/ 31 DATA X128B /0.9,.80,.60,.25/ 32 DATA Y28B/ 1.1,1.13,1.17,1.20,1.24,1.27,1.3,1.33,1.34,1.35,1.36 , 33 11.0,1.04,1.08,1.11,1.15,1.18,1.21,1.23,1.25,1.25,1.26 , 0.88, .92, 34 2.96,1.0,1.04,1.08,1.11,1.13,1.14,1.14,1.15 , 0.81,.85,.89,.925, 35 3.96,1.0,1.03,1.05,1.06,1.06,1.07/ 36C 37C--------- FIGURE 4.1.5.1-28B-DASHED --------- 38C---- X228BD= COS(T/C)MAX X128BD= MACH NUMBER 39C---- Y28BD= (R)LS 40C 41 DATA X228BD /0.45,.65,.70,.75,.80,.85,.90,.95,1.0/ 42 DATA X128BD /0.9,.80,.60,.25/ 43 DATA Y28BD /1.20,1.20,1.24,1.27,1.30,1.33,1.34,1.35,1.36 , 1.11, 44 11.11,1.15,1.18,1.21,1.23,1.25,1.25,1.26 , 1.0,1.0,1.04,1.08,1.11, 45 21.13,1.14,1.14,1.15 , 0.925,0.925,.96,1.0,1.03,1.05,2*1.06,1.07/ 46C 47 D(2)=12.0*A(16)/ROUGFC 48C 49C--------- DETERMINE WING CONFIGURATION ---------------------- 50C IF NOT STRAIGHT-TAPERED, GO TO PAGE 3-2 51C 52 IF(VT(15).NE.WTYPE(1)) GO TO 1010 53C 54C--------- PAGE 3-1 CALCULATIONS ---------------------- 55C********* FOR STRAIGHT-TAPER WING CONFIGURATION ************ 56C--------- COMPUTE REYNOLDS NUMBER, (RN) 57C 58 1000 D(17)=A(16)*A(129) 59C 60C--------- COMPUTE RL FROM EQUATION FOR FIG.4.1.5.1-27 61C********* FIGURE 4.1.5.1-27 EQUATION 62C 63 CALL TBFUNX(B(I+2),CEPT,DYDX,4,X27M,X27I,C27,I27,MI,NG, 64 1 0,0,Q15127,3,ROUTID) 65 D(14)=D(2)**1.0482*10.0**CEPT 66C 67C********* FIGURE 4.1.5.1-26 EQUATIONS (SUBROUTINE FIG26 ) 68C 69 IF( D(14).LT.D(17)) D(17)=D(14) 70 CALL FIG26 (D(17),B(I+2),D(10) ) 71C 72C********* FIGURE 4.1.5.1-28-B 73C 74 CALL TLINEX(X128B,X228B,Y28B,4,11,B(I+2),A(178),D(13), 75 1 0,2,0,2,Q5128B,3,ROUTID) 76 CAPL=2.0 77 IF(VT(18).GE.0.30) CAPL=1.20 78 D(20)=D(10)*(1.+CAPL*VT(16)+100.*VT(16)**4)*D(13)*(2.*A(3)/SREF) 79 GO TO 1030 80C 81C ***** INBOARD CALCUALTIONS ***** 82C 83 1010 D(18)=A(15)*A(129) 84C 85C********* FIGURE 4.1.5.1-27 EQUATION 86C 87 CALL TBFUNX(B(I+2),CEPT,DYDX,4,X27M,X27I,C27,I27,MI,NG, 88 1 0,0,Q15127,3,ROUTID) 89 D(15)=D(2)**1.0482*10.0**CEPT 90C 91C********* IF RLI.LT.RNI SET RNI=RLI 92C 93 IF( D(15).LT.D(18) ) D(18)= D(15) 94C 95C********* FIGURE 4.1.5.1-26 SUBROUTINE 96C 97 CALL FIG26 ( D(18),B(I+2),D(11) ) 98C 99C--------- FIGURE 4.1.5.1-28-B (DASHED) 100C 101 CALL TLINEX(X128BD,X228BD,Y28BD,4,9,B(I+2),A(190),D(23), 102 1 0,0,0,0,Q128BD,5,ROUTID) 103C 104C********* IF XTI GREATER TAN 0.30, THEN CAPL= 2.0 * 105C 106 1020 CAPL= 1.20 107 IF(VT(18).LT.0.30) CAPL=2.00 108 D(21)=D(11)*(1.+CAPL*VT(16)+100.*VT(16)**4)*D(23)*(2.*A(1)/SREF) 109C 110C ***** OUTBOARD CALCULATIONS ***** 111C 112 D(19)=A(17)*A(129) 113C 114C********* FIGURE 4.1.5.1-27 EQUATION 115C 116 D(16)=D(15) 117C 118C********* IF RLO .LT. RNO, SET RNO= RLO 119C 120 IF(D(16).LT.D(19)) D(19)=D(16) 121C 122C********* FIGURE 4.1.5.1-26 SUBROUTINE 123C 124 CALL FIG26( D(19),B(I+2),D(12) ) 125C 126C********* FIGURE 4.1.5.1-28-B 127C 128 CALL TLINEX(X128B,X228B,Y28B,4,11,B(I+2),A(184),D(24), 129 1 0,2,0,2,Q5128B,3,ROUTID) 130C 131C********* IF XTO GREATER THAN 0.30, THEN CAPL= 2.0 132C 133 CAPL=1.20 134 IF(VT(66).LT.0.30) CAPL=2.00 135 D(22)=D(12)*(1.+CAPL*VT(65)+100.*VT(65)**4)*D(24)*(2.*A(2)/SREF) 136C 137C********* COMPUTE CDO FROM PAGE 3-2 138C 139 D(20)=D(21)+D(22) 140 1030 CONTINUE 141 RETURN 142 END 143