1      SUBROUTINE SUPHB
2C
3C***  CALCULATES SUPERSONIC HORIZONTAL TAIL-BODY CL,CD,CLA AND CMA
4C
5      REAL MACH,NXX,KWB,KBW,IVBW  ,KKBW,KKWB
6      DIMENSION ROUTID(2),Q1211A(3),Q1211B(3),Q31210(3),Q1212A(3),
7     1          Q1212B(3),Q2118A(3),Q2118B(3),Q2137A(3),Q2137B(3)
8      DIMENSION T4312B(8),D4312B(8)
9      DIMENSION CDW(20),CDB(20),WTYPE(4),IVBW(20),GAMMA(20)
10      DIMENSION T4337A(11),D4337A(24),T4337B(12),D4337B(20)
11      DIMENSION T4312A(11),D4312A(11)
12      DIMENSION TFIG10(11),DKWB10(11),DKBW10(11)
13      DIMENSION T4218A(18),DL218A(54),DR218A(54),T4218B(14),DL218B(42),
14     1          DR218B(42)
15      DIMENSION LGH(4),VAR(4),CD(20),CN(20),CA(20),CL(20),CLB(20),
16     1          CLW(20),ALPHAB(20),CDL(20)
17      DIMENSION CM(20)
18C
19      COMMON /CONSNT/ PI,DEG,UNUSED,RAD
20      COMMON /FLGTCD/ FLC(160)
21      COMMON /SYNTSS/ SYNA(19)
22      COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF
23      COMMON /HTI/    W(2),SPANS,SPAN,W1,CR,W2(8),TYPE
24      COMMON /BODYI/  NXX,XCOOR(20),S(20),P(20),R(20),ZU(20),ZL(20),
25     1                BNOSE,BTAIL,RLN
26      COMMON /BDATA/  BD(762)
27      COMMON /HTDATA/ AHT(195)
28      COMMON /OVERLY/ NLOG,NMACH,I,NALPHA
29      COMMON /SUPBOD/ SBD(229)
30      COMMON /SUPWBB/  SWB(61),SHB(61)
31      COMMON /SUPWH/  SLG(141), STG(141)
32      COMMON /IBODY/  PBODY, BODY(400)
33      COMMON /IHT/    PHT, HT(380)
34      COMMON /IBH/    PBH, BH(380)
35C
36      EQUIVALENCE (CA(1),BH(81)),(CDW(1),HT(1)),(CMA,BH(121)),
37     1(DXCG,AHT(173)),(CLAB,SBD(18)),(KKWB,SHB(2)),
38     2(XACN,SHB(3)),(CDOWB,SHB(4)),(DELXW,BD(89)),(ARSTAR,AHT(7)),
39     3(CD(1),BH(1)),(DD,SHB(5)),(XW,SYNA(6)),(CLB(1),BODY(21)),
40     4(ALIW,SYNA(8)),(BETA,SHB(6)),(CLAW,HT(101)),
41     5(XACW,STG(134)),(CLABW,SHB(7)) ,(XACBW,SHB(8))
42      EQUIVALENCE (CLA,BH(101)), (CDL(1),STG(53))
43      EQUIVALENCE (CL(1),BH(21)),(FA,SHB(9)),(CLI,SHB(10)),
44     1(KBW,SHB(11)),(CDOB,SBD(124)),(IVBW(1),SHB(12)),(RKBW,SHB(32)),
45     2(CLAWB,SHB(33)),(TANL ,AHT(62)),(CRSTAR,AHT(10)),
46     3(TAPEXP,AHT(27)),(CN(1),BH(61)),(FN, SHB(34)),(CDB(1),BODY(1)),
47     4(CLW(1),HT(21)),(KWB,SHB(35)),(XAC,SHB(36)), (CDOW,STG(80)),
48     5(DN,SBD(4)),(D1,SBD(5)),(KKBW,SHB(37)),(RLAP,SHB(38))
49      EQUIVALENCE (XACA,SHB(39)),(GAMMA(1),SHB(40)),
50     1(TRINO,SHB(60)),(XCPLN,SHB(61)),(ALPHAB(1),BD(255))
51     2,(CM(1),BH(41))
52C
53      DATA ROUTID/4HSUPW,4HB   /,
54     2 Q31210/4H4.3.,4H1.2-,4H10  /,Q1212A/4H4.3.,4H1.2-,4H12A /,
55     3 Q2118A/4H4.2.,4H2.1-,4H23A /,Q1212B/4H4.3.,4H1.2-,4H12B /,
56     4 Q2118B/4H4.2.,4H2.1-,4H23B /,Q2137A/4H4.3.,4H2.2-,4H37A /,
57     5                              Q2137B/4H4.3.,4H2.2-,4H37B /
58      DATA WTYPE/4HSTRA,4HDOUB,4HCRAN,4HCURV/
59C
60C                   FIGURE 4.3.1,2-10 KWB
61C
62      DATA TFIG10/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/
63      DATA DKWB10/1.0,1.08,1.16,1.26,1.36,1.46,1.56,1.67,1.78,1.89,2.0/
64C
65C                   FIGURE 4.3.1.2-10 KBW
66C
67      DATA DKBW10/0.0,.13,.29,.45,.62,.80,1.0,1.22,1.45,1.70,2.0/
68C
69C               FIGURE 4.2.2.1-23A (LEFT SIDE)
70C
71      DATA T4218A
72     1  / 0., 0.2, 0.4, 0.6, 0.8, 1.0 , 3*0.0,
73     2    0.0, 0.4, 0.8, 1.2, 1.6, 2.0, 3.0, 4.0, 5.0 /
74      DATA DL218A
75     1 /  0.543,  0.542,  0.541,  0.540,  0.534,  0.526,
76     2    0.400,  0.409,  0.418,  0.430,  0.441,  0.450,
77     3    0.305,  0.328,  0.350,  0.369,  0.387,  0.400,
78     4    0.238,  0.265,  0.295,  0.318,  0.339,  0.356,
79     5    0.198,  0.221,  0.246,  0.274,  0.298,  0.320,
80     6    0.160,  0.185,  0.210,  0.239,  0.262,  0.288,
81     7    0.065,  0.095,  0.122,  0.150,  0.177,  0.210,
82     8    0.000,  0.005,  0.035,  0.062,  0.089,  0.130,
83     9    0.000,  0.000,  0.000,  0.000,  0.002,  0.050  /
84C
85C               FIGURE 4.2.2.1-23A (RIGHT SIDE)
86C
87      DATA DR218A
88     1  / 0.445,  0.464,  0.485,  0.500,  0.518,  0.526,
89     2    0.448,  0.455,  0.460,  0.460,  0.459,  0.450,
90     3    0.460,  0.449,  0.438,  0.424,  0.412,  0.400,
91     4    0.450,  0.430,  0.412,  0.394,  0.375,  0.356,
92     5    0.432,  0.410,  0.388,  0.365,  0.343,  0.320,
93     6    0.420,  0.394,  0.369,  0.340,  0.314,  0.288 ,
94     7    0.388,  0.354,  0.322,  0.278,  0.244,  0.210,
95     8    0.357,  0.314,  0.273,  0.216,  0.171,  0.130,
96     9    0.325,  0.274,  0.225,  0.154,  0.100,  0.050  /
97C
98C              FIGURE 4.2.2.1-23B (LEFT SIDE)
99C
100      DATA T4218B
101     1  / 0., 0.2, 0.4, 0.6, 0.8, 1.0  , 0.,
102     2    0.0, 0.5, 1.0, 2.0, 3.0, 4.0, 5.0 /
103      DATA DL218B
104     1 /  .665,  .665,  .665,  .665,  .665,  .665,
105     2    .425,  .492,  .539,  .550,  .550,  .550,
106     3    .330,  .370,  .405,  .438,  .459,  .470,
107     4    .184,  .215,  .250,  .284,  .318,  .350,
108     5    .060,  .097,  .133,  .170,  .206,  .240,
109     6    .000,  .000,  .044,  .083,  .127,  .170,
110     7    .000,  .000,  .000,  .020,  .063,  .105  /
111C
112C              FIGURE 4.2.2.1-23B (RIGHT SIDE)
113C
114      DATA DR218B
115     1 /  .665,  .665,  .665,  .665,  .665,  .665,
116     2    .480,  .500,  .519,  .536,  .546,  .550,
117     3    .338,  .388,  .430,  .458,  .471,  .470,
118     4    .338,  .372,  .394,  .394,  .375,  .350,
119     5    .410,  .375,  .341,  .308,  .272,  .240,
120     6    .377,  .338,  .294,  .251,  .211,  .170,
121     7    .342,  .300,  .246,  .194,  .146,  .100   /
122C
123C                   FIGURE 4.3.1.2-12A (KBW)
124C
125      DATA T4312A/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/
126      DATA D4312A/0.0,.11,.21,.31,.41,.51,.6,.7,.8,.9,1.0/
127C
128C                   FIGURE 4.3.1.2-12B (KWB)
129C
130      DATA T4312B/.015,.1,.2,.3,.4,.6,.8,.975/
131      DATA D4312B/1.,.975,.956,.947,.941,.950,.978,1.0/
132C
133C                   FIGURE 4.3.2.2-37A
134C
135      DATA T4337A/0.0,.4,.8,1.2,1.6,2.,2.4,2.8,
136     1           .1,1.0,999999./
137      DATA D4337A/
138     1.5,.72,.900,1.08,1.24,1.39,1.53,1.68,
139     2.5,.72,.910,1.09,1.25,1.41,1.57,1.72,
140     3.5,.73,.920,1.11,1.27,1.43,1.59,1.74/
141C
142C                   FIGURE 4.3.2.2-37B
143C
144      DATA T4337B/0.0,.1,.2,.3,.4,.5,.6,.8,1.0,2.8,
145     1            0.2,999999./
146      DATA D4337B/
147     10.5,.56,.595,.62,.64,.65,.66,.669,.669,.671,
148     20.5,.54,.578,.60,.62,.638,.649,.66,.669,.671/
149C
150      NX=NXX+.5
151      DCYL=(DN+D1)/2.
152      MACH=FLC(I+2)
153      BETA=SQRT(MACH**2-1.)
154      RLB=XCOOR(NX)
155      DD=2.0*(SPAN-SPANS)
156      TANLE=TANL
157      IF(TANLE.EQ.0.0)TANLE=.00001
158C
159C   ***SUPERSONIC H.T.-BODY LIFT CURVE SLOPE,BODY IN PRESENCE OF H.T.***
160C                         NON-TRIANGULAR H.T.S
161C
162      IF(TAPEXP.EQ.0.0)GO TO 1050
163      ARG1=BETA*ARSTAR*(1.0+TAPEXP)
164      ARG2=1.+TANLE/BETA
165      TRINO=ARG1*ARG2
166      IF(TRINO.LE.4.)GO TO 1030
167 1000 OLE=AHT(34)
168      DX=RLB-XW-DD*AHT(38)-CRSTAR
169      IF(DX.LE.(-CRSTAR))RKBW=0.
170      IF(DX.LE.(-CRSTAR))GO TO 1020
171      CALL INTKBW(MACH,OLE,CRSTAR,DD,DX,RKBW,RXAC)
172 1020 KBW=RKBW/(RAD*BETA*(SREF/AHT(3))*CLAW*(TAPEXP+1.)*(2.*SPAN/DD-1.))
173      GO TO 1040
174 1030 CONTINUE
175      LGH(1)=11
176      VAR(1)=DD/(2.*SPAN)
177C
178C                   FIGURE 4.3.1.2-10 KBW
179C
180      CALL INTERX(1,TFIG10,VAR,LGH,DKBW10,KBW,11,11,
181     1            0,0,0,0,0,0,0,0,Q31210,3,ROUTID)
182 1040 CONTINUE
183      GO TO 1060
184C
185C  ***SUPERSONIC H.T.-BODY LIFT CURVE SLOPE,BODY IN PRESENCE OF H.T.***
186C                         TRIANGULAR H.T.
187C
188 1050 CONTINUE
189      ARG=BETA*ARSTAR
190      IF(ARG.GT.1.)GO TO 1000
191      GO TO 1030
192 1060 CONTINUE
193C
194C  ***SUPERSONIC H.T.-BODY LIFT CURVE SLOPE
195C
196      ALBO=BD(81)
197      IF(BD(81).EQ.UNUSED)ALBO=0.0
198      DO 1070 J=1,NALPHA
199 1070 ALPHAB(J)=FLC(J+22)+ALBO
200      VAR(1)=(SPAN-SPANS)/SPAN
201      LGH(1)=11
202C
203C                   FIGURE 4.3.1.2-10 KWB
204C
205      CALL INTERX(1,TFIG10,VAR,LGH,DKWB10,KWB,11,11,
206     1            0,0,0,0,0,0,0,0,Q31210,3,ROUTID)
207      CLAWB=CLAW*KWB
208      CLABW=CLAW*KBW
209      CLA=CLABW+CLAWB+CLAB
210C
211C  ***SUPERSONIC H.T.-BODY LIFT AT ANGLE OF ATTACK***
212C
213      IF(TYPE.NE.WTYPE(1))GO TO 1100
214      IF(ALIW.EQ.0.0.OR.(ALIW.EQ.UNUSED))GO TO 1080
215C
216C                   FIGURE 4.3.1.2-12A (KBW) INCIDENCE
217C
218      CALL INTERX(1,T4312A,VAR,LGH,D4312A,KKBW,11,11,
219     1            0,0,0,0,0,0,0,0,Q1212A,3,ROUTID)
220C
221C                   FIGURE 4.3.1.2-12B (KWB) INCIDENCE
222C
223      LGH(1)=8
224      CALL INTERX(1,T4312B,VAR,LGH,D4312B,KKWB,8,8,
225     1            0,0,0,0,0,0,0,0,Q1212B,3,ROUTID)
226 1080 CLI=CLAW*ALIW
227      DETCL=0.0
228      XCBO4=XW+AHT(161)
229      BMBS=SPAN-SPANS
230      CALL BODOWG(BD(255),XCBO4,BMBS,SPAN,AHT(27),
231     1            IVBW,GAMMA,NALPHA)
232      ARG1=KWB+KBW
233      ARG2=(KKWB+KKBW)*CLI
234      ARG3=CLAW*(DD/(2.*SPAN))
235      DO 1090 J=1,NALPHA
236         DETCL=ARG3*ALPHAB(J)*IVBW(J)*GAMMA(J)
237 1090 CL(J)=CLB(J)+ARG1*(CLW(J)-CLI)+ARG2+DETCL
238C
239C  ***SUPERSONIC CENTER OF PRESSURE FOR BODY NOSE AND FORBODY***
240C
241 1100 CONTINUE
242      DELXW=(SPAN-SPANS)*TANLE*COS(ALIW/RAD)
243      RLAP=XW+DELXW-RLN
244      ARG1=RLN
245      IF(RLAP.LT.0.)ARG1=RLN+RLAP
246      IF(RLAP.LT.0.)RLAP=0.0
247      FA=RLAP/DCYL
248      FN=ARG1/DCYL
249      VAR(1)=BETA/FN
250      VAR(2)=FA/FN
251      LGH(1)=6
252      LGH(2)=9
253      IF(BNOSE.EQ.1.)GO TO 1120
254C
255C                   FIGURE 4.2.2.1-23A LEFT SIDE(XCP) OGIVE
256C
257      IF(VAR(1).GT.1.)GO TO 1110
258      CALL INTERX(2,T4218A,VAR,LGH,DL218A,XCPLN,9,54,
259     1   0,0,0,0,0,0,0,0,Q2118A,3,ROUTID)
260      GO TO 1140
261 1110 VAR(1)=1./VAR(1)
262C
263C                   FIGURE 4.2.2.1-23A RIGHT SIDE(XCP) OGIVE
264C
265      CALL INTERX(2,T4218A,VAR,LGH,DR218A,XCPLN,9,54,
266     1   0,0,0,0,0,0,0,0,Q2118A,3,ROUTID)
267      GO TO 1140
268 1120 CONTINUE
269C
270C                   FIGURE 4.2.2.1-23B LEFT SIDE(XCP) CONE
271C
272      LGH(2)=7
273      IF(VAR(1).GT.1.)GO TO 1130
274      CALL INTERX(2,T4218B,VAR,LGH,DL218B,XCPLN,7,42,
275     1            0,0,0,0,0,0,0,0,Q2118B,3,ROUTID)
276      GO TO 1140
277 1130 VAR(1)=1./VAR(1)
278C
279C                   FIGURE 4.2.2.1-23B RIGHT SIDE(XCP) CONE
280C
281      CALL INTERX(2,T4218B,VAR,LGH,DR218B,XCPLN,7,42,
282     1            0,0,0,0,0,0,0,0,Q2118B,3,ROUTID)
283 1140 CONTINUE
284      XACN=(XCPLN-1.)*(ARG1+RLAP)/CBARR
285C
286C  ***SUPERSONIC H.T.-LIFT CARRYOVER ON BODY***
287C
288      VAR(1)=BETA*DD/CRSTAR
289      VAR(2)=BETA/TANLE
290      ARG=(XW+CR)/RLB
291      IF(ARG.GT.1.)GO TO 1150
292C
293C                   FIGURE 4.3.2.2-37A(XAC)B(W)
294C
295      LGH(1)=8
296      LGH(2)=3
297      CALL INTERX(2,T4337A,VAR,LGH,D4337A,XACA,8,24,
298     1            0,0,0,0,1,0,0,0,Q2137A,3,ROUTID)
299      GO TO 1160
300C
301C                   FIGURE 4.3.2.2-37B(XAC)B(W)
302C
303 1150 LGH(1)=10
304      LGH(2)=2
305      CALL INTERX(2,T4337B,VAR,LGH,D4337B,XACA,10,20,
306     1            0,0,0,0,0,0,0,0,Q2137B,3,ROUTID)
307 1160 XACBW=XACA*CRSTAR/CBARR
308      XACW=XACW*CRSTAR/CBARR
309C
310C  ***SUPERSONIC H.T.-BODY PITCHING MOMENT SLOPE***
311C
312      DNUM=XACN*CLAB+XACW*CLAWB+XACBW*CLABW
313      DNOM=CLAB+CLAWB+CLABW
314      XAC  =DNUM/DNOM
315      CMA=(DXCG/CBARR-XAC)*CLA
316C
317C  ***SUPERSONIC H.T.-BODY DRAG***
318C
319      CDOWB=CDOW+CDOB
320      DO 1170 J=1,NALPHA
321         COSA=COS(FLC(J+22)/RAD)
322         SINA=SIN(FLC(J+22)/RAD)
323         CD(J)=CDOW+CDL(J)+CDB(J)
324         CN(J)=CL(J)*COSA+CD(J)*SINA
325C
326C   COMPUTE DELTA XCP WING-BODY
327C
328         DXCPWB=DXCG/CBARR-SLG(134)
329C
330C   COMPUTE DELTA XCP BODY-WING
331C
332         DXCPBW=DXCG/CBARR-SWB(8)
333C
334C   COMPUTE DELTA CN VORTEX
335C
336         DCNV=SWB(11+J)*SWB(39+J)*SWB(5)/(2*SPAN)*BD(254+J)*HT(101)
337C
338C   COMPUTE WING-BODY PITCHING MOMENT
339C
340         BH(40+J)=BODY(40+J)+HT(60+J)*SWB(35)*DXCPWB
341     1            +HT(101)*SYNA(4)*SWB(2)*DXCPWB/RAD
342     2            +HT(60+J)*SWB(11)*DXCPBW
343     3            +HT(101)*SYNA(4)*SWB(37)*DXCPBW/RAD
344     4            +DCNV*DXCPBW
345     5            +HT(80+J)*(SYNA(3)-SYNA(5))/CBARR
346 1170 CA(J)=CD(J)*COSA-CL(J)*SINA
347      RETURN
348      END
349