Clculo de la Curva Polar de un SAAB JAS 39 GRIPPEN

V Ao Ingeniera Aeronutica, Aerodinmica

Resumen
En el presente trabajo se obtendr la curva polar de la aeronave SAAB JAS 39 GRIPEN. Material areo polivalente creado especialmente para suplir la necesidad de cazas de la fuerza area sueca. Presentado originalmente por la empresa SAAB en la dcada de los 80, el JAS es una avin con caractersticas similares a las del F-16. Para poder obtener la curva polar de esta aeronave, se utilizarn aproximaciones geomtricas, con el fin de obtener valores cercanos a la real superficie mojada. Con estos valores y ecuaciones presentados durante el ramo de aerodinmica I, se obtendrn datos y curvas necesarios para llegar a obtener el objetivo de este trabajo.

1. Introduccin 2. Marco Terico

En el estudio del ramo de aerodinmica, se hace imperante el anlisis de efectos fenomenolgicos que ocurren en el vuelo de una aeronave. Dentro del diseo, es de especial importancia obtener ciertos datos que permitan mostrar los coeficientes de sustentacin y arrastre de tal aeronave, es por esto que se obtiene una curva en donde se enfrentan ambos resultados. En este Paper, se demostrara como obtener dicha curva para el avin SAAB JAS 39 GRIPPEN. Para la obtencin de la Curva Polar del avin JAS 39 se utilizarn los conocimientos adquiridos en el ramo de Aerodinmica. A continuacin se detalla la nomenclatura a usar en el desarrollo:
Tabla 1 Nomenclatura utilizada en el desarrollo del trabajo. Nomenclatura CD CL : Coeficiente de resistencia. : Coeficiente de sustentacin. : Densidad. : Velocidad. : Viscosidad. : ngulo de flecha en el espesor mximo. : ngulo de flecha para articulaciones de flaps. : Razn de aspecto. : Envergadura. : Superficie. : Factor de eficiencia de envergadura. : ngulo de ataque. : Distancia 1/4 de cuerda de superficie de ala y canard

tmax HL
Fig. 1 JAS 39 GRIPPEN AR b S e

El Saab 39 Gripen, designado JAS 39 Gripen por la Fuerza Area Sueca, es un caza polivalente monomotor desarrollado en los aos 1980 por el consorcio de empresas suecas IG-JAS, liderado por la compaa aeroespacial Saab. El fin era equipar a la Fuerza Area de Suecia con un caza, un avin de ataque y un avin de reconocimiento areo, en sustitucin de los aviones Saab 35 Draken y Saab 37 Viggen, y eventualmente al entrenador SK 60 cuando su sustituto fue cancelado. Con un nivel similar al del F16 estadounidense.

Swet

: Razn de ahusamiento. : Superficie mojada.

Mecnica de Vuelo

Sexposed t Cfe e0 h c w l

: Superficie en el plano. : Espesor mximo. : Coeficiente de resistencia de friccin equivalente. : Coeficiente de Oswald. : Altura. : Cuerda. : Ancho. : Longitud.

Despus de haber analizado con el Profesor (CDE (I) Ricardo Sez) distintas posibilidades, se tom en comn acuerdo la utilizacin de dos perfiles distintos para canard y ala. Obtenidas las curvas del libro THEORY OF WING SECTIONS, IRA [Link], respectivamente se utilizaron los perfiles NACA 64-206 [fig. 2] y NACA 64-210 (que presenta resultados en la pendiente similares al F-16)[fig. 3]. Para calcular la razn de aspecto, tanto para el ala y canard es: (1)

Adems, es necesario determinar el rea mojada de la aeronave descrita anteriormente para lo cual, usaremos los siguientes datos mostrados en la Tabla 1(Fig. 2):
Tabla 2 Especificaciones del avin JAS 39 GRIPEN

Para comenzar el clculo de sustentacin de la aeronave, lo primero que debe considerarse es el factor de eficiencia, esto es: (2)

Envergadura Alar Superficie Alar Superficie Estabilizador Horizontal ngulo de flecha

27,88 ft 183.9 ft2 34.93 ft2 3

Luego con los datos obtenidos en la curva de sustentacin de las cartas NACA, especficamente la pendiente (

El autor, para poder obtener la geometra de esta aeronave se vio apoyado por un plano esquemtico encontrado en internet, el cual utiliz debido a que presenta dimensiones similares a las reales en escala.

(3)

Se reemplaza este valor en la ecuacin (4) (4) De la geometra de la aeronave se debe obtener el valor de , as tambin y cavg, esto es: (5) (6) Posteriormente se calcula la variacin de downwash respecto del ngulo de ataque. Para ello se tiene lo siguiente:
Fig. 2 Plano esquemtico

Esta aeronave, al ser de un consorcio europeo, present algunas dificultades para poder obtener informacin de los respectivos perfiles alares. El autor envi diversos correos a las empresas GRIPEN y ARA (Aircraft Research Association), sin respuesta alguna.

( (7)

Mecnica de Vuelo

El coeficiente de arrastre se calcular de la siguiente manera: (12) Continuando con el clculo de la resistencia debida a la sustentacin, se calcula el coeficiente de Oswald, en nuestro caso, para LE< 30: (13) Entonces con el valor conseguido se logra: (14) (8) Obtenido el resultado, se calcular el valor mximo para el coeficiente de sustentacin de la aeronave ya sea en despegue como en aterrizaje, esto es: (9) (10) Siendo h.l. el ngulo de flecha del lugar geomtrico de las articulaciones del flap. Una vez conseguido el a para el despegue y el aterrizaje y con la ecuacin siguiente se puede calcular el coeficiente de sustentacin mximo utilizable: (11) Luego, para la obtencin de la resistencia parsita del avin se utilizar el mtodo de aproximacin mediante figuras geomtricas sencillas, como se muestra en las siguientes figuras: Y con esto se logra obtener el coeficiente de resistencia de la aeronave en funcin de CL: (18) Debido a que este coeficiente de resistencia es vlido hasta velocidades inferiores a , se procede a obtener el y y supersnico: (19) (20)
( )

Fig. 3 Distancia lc

Luego se calcula la pendiente de la curva CL v/s , lo cual se har de la siguiente manera: ms la contribucin del canard.

Posteriormente, se calcular un nmero de Reynolds para un M dado y una condicin asl: (15) Con esto se genera una tabla de datos para generar la curva CL v/s CD. Con los datos entregados por la tabla se calcula el CD0 del avin, para ello se utiliza la siguiente ecuacin: (16) Luego se calcula el valor de k2: (17)

(21) (22)

Fig. 4 Esquema geomtrico

Ya obtenido el valor de la superficie mojada de la aeronave se pueden realizar clculos que servirn para la elaboracin de la curva polar.

Con tal funcin se obtendrn los datos para generar la curva de coeficiente de resistencia de Onda VS M.

Mecnica de Vuelo

Se procede a obtener

supersnico.
( )

(23)

Para calcular la variacin del downwash respecto al ngulo de ataque, es necesaria la razn de ahusamiento, la cuerda media, la distancia vertical entre los planos de ala y estabilizador horizontal y la distancia entre cuartos de cuerda entre la superficie alar y el estabilizador horizontal respectivamente:

Ya obtenidos todos los datos anteriores, se procede a calcular la resistencia de perfil del avin, siendo este la suma del coeficiente de resistencia parasita ms el coeficiente de resistencia de onda: (24) Esta ecuacin da origen a una familia de curvas polares para diferentes nmeros de Mach para el SAAB JAS 39 GRIPEN.

De (7):

Calculando la pendiente de la curva CL vs para el avin completo, se reemplaza en (8):

3. Resultados
Reemplazando los datos expuestos anteriormente, se puede obtener una cantidad de datos que se utilizarn para el clculo de la curva polar de la aeronave seleccionada: Reemplazando en (1) obtenemos la razn de aspecto, tanto para el ala como para el estabilizador: AR=2,802; ARc=2,39 Con estos valores se obtiene el factor de eficiencia para las alas, considerando que el avin posee distintos perfiles en la raz y en la punta; para ello se promedia el espesor mximo. Para el NACA 2418, el espesor mximo corresponde al 18% de la cuerda y para el NACA 2412, al 12% de la misma; el promedio resulta 15% de la cuerda ubicado en el 30% de sta, como en todos los perfiles NACA de 4 dgitos. Reemplazando en (2): Una vez conseguido el a para el despegue y el aterrizaje y con la ecuacin siguiente se puede calcular el coeficiente de sustentacin mximo utilizable, reemplazando en (11): Con estos valores obtenidos, se calcula el valor mximo para el coeficiente de sustentacin del avin, tanto en configuracin de despegue como aterrizaje, sustituyendo en (9) y en (10):

Ahora, con los modelos aproximados de la aeronave se calcular la resistencia parsita, para ello se utilizar la tabla 1: Tabla 1 Datos presentados por similitud geomtrica.
Superficie s1 y 2 s3 y 4 s5 Cilindros c1 c2 Conos o1 Largo Alto1 6,56 4,27 3,28 5,90 6,56 Largo Alto1 27,88 27,88 3,28 3,28 3,28 3,28 3,28 4,58 4,58 Largo 34,44 3,28 Envergadura 10,49 4,27 7,21 c root 16,40 7,21 11,15 Alto 4,92 3,92 c tip 3,28 1,64 2,29 Ancho 3,28 2,62 t max/c S wet S- INTERSE S WET NETO 0 0 0 370,94 206,87 103,10 0,04 370,94 0,04 206,87 0,06 103,10

De las cartas NACA se obtiene la pendiente de la curva CL vs , dicha pendiente se calcular a partir de los perfiles alares y se ampliar para el estabilizador horizontal.

S wet S- INTERSE S WET NETO 443,61 16,85 0 0 443,61 16,85

Ancho 1 Alto2 3,28 6,56 3,60 2,62 2,62 Ancho 1 1,64 1,64 0,00 3,28 3,28 0,00 0,00

Ancho 2 S wet S- INTERSE S WET NETO 0,00 5,58 2,62 0,00 0,00 33,80 62,64 32,92 16,68 18,55 0 0 0 0 0 33,80 62,64 32,92 16,68 18,55

Sustituyendo valores en (4):

o2 o3 1/2o1 1/2o2 Paralelo p1 p2 Total

137,17 137,17

0 0

137,17 137,17 1580,30

Mecnica de Vuelo

A continuacin se calcular CDmin con la ecuacin (12), considerando Cfe= 0.0035, quedando entonces:

Luego continuamos calculando el factor de eficiencia de Oswald, sustituyendo en (13):

Entonces, reemplazando los nuevos datos obtenidos se reemplaza en (14):

Fig. 5 Cd en funcin de Cl

Teniendo los datos, calculamos el nmero de Reynolds para una condicin asl y nmero de Mach M=0.2:

Calculando la derivada de la ecuacin encontramos el valor mnimo, es decir: CLminD=0.0176 Luego, reemplazamos en (16):

Debemos recordar que tenemos para hacer un clculo de coeficiente de arrastre tenemos dos perfiles alares, sin embargo observando los valores en ambos grficos de momento la diferencia es despreciable, por lo que se considerar el perfil NACA 2418 que tiene un coeficiente de momento ligeramente mayor: De la carta del perfil NACA 2418 (Fig.) y ocupando un valor de Reynolds aproximado Re=9000000, se genera la tabla siguiente: Tabla 2 Obtencin de datos necesarios para obtener la curva polar.

Y calculamos k2 (17):

Con lo cual el coeficiente de resistencia del avin es (en funcin del CL):

Para variados CL:

Cl -0,2 -0,1 0 0,1 0,2 0,3

cd 0,0063 0,006 0,0055 0,0045 0,0045 0,0043

k1 Cl^2 0,0054 0,0013 0 0,0013 0,0054 0,0121

Cd 0,0117 0,0073 0,0055 0,0058 0,0099 0,0164

Fig. 6 curva polar JAS 39

Debido a que este coeficiente de resistencia es vlido hasta velocidades inferiores a , se procede a obtener el y y (19, 20, 21, 22):

Mecnica de Vuelo

Para se procede a formar la curva a variados M, que nos permitirn observar la resistencia de Onda.

Fig. 9 Resistencia VS velocidad en Mach

Fig. 7 resistencia de onda a variados nmeros de Mach

Posterior a esto, el coeficiente de resistencia total del avin, en toda su envolvente de vuelo, viene definido por la siguiente relacin de CD. Esto da origen a una familia de curvas Polares con diferentes nmeros de MACH. (Fig.10)

Obtenemos de Mach.

supersnico para variados nmeros

Fig. 8 variacin de k1 frente a la velocidad

Finalizado esto se obtiene coeficiente de resistencia de perfil. El cual se obtiene con la siguiente ecuacin y es la suma del coeficiente de resistencia parsita ms el coeficiente de resistencia de onda:

Fig. 10 Familia de curvas polares

4. Discusin
Durante este trabajo se ha obtenido la curva polar del avin JAS 39 GRIPEN, segn las ecuaciones y forma planteadas por el ramo de MEC. DE VUELO. Pero los datos obtenidos correspondieron a aproximaciones, debido a ciertos factores: Las formas geomtricas son similares al original, por lo cual se supone un rea mojada similar a la real. El problema de no poder obtener la informacin de los perfiles que cubren la superficie alar y canard, produjo que el autor, buscara perfiles similares, esto causa que los coeficientes no sean los reales.

Todos estos problemas y la reduccin en los clculos junto con sus aproximaciones, alejan los resultados, en un porcentaje importante, de la realidad.

5. Conclusiones
En el presente trabajo se desarroll el clculo de la curva polar del avin JAS 39 GRIPPEN utilizando

como referencia el documento entregado por el profesor de la asignatura, en donde se presenta un ejemplo de la obtencin de la curva polar del F-16. Para poder obtener los datos de las superficies se utilizaron planos esquemticos encontrados en internet. Por lo que los valores obtenidos de superficie mojada del avin mediante el clculo de figuras geomtricas se suponen cercanos a los reales, al no obtener informacin propia del avin. Diferentemente a los resultados reales, se puede apreciar que se obtienen curvas similares a las del F-16, este avin (JAS 39), fue creado con caractersticas similares, por lo que se puede apreciar resultados no tan escapados a los del F-16. Al final del trabajo se obtienen los resultados de curvas polares para diferentes MACH, aqu se puede apreciar el salto que ocurre al aumentar el MACH y todo el arrastre que puede traer con l, suponiendo con esto, un aumento en el consumo de combustible, que pueda afectar algunas performances como su autonoma. Si se pasa de un rgimen subsnico (0.2M)a transnico (0.8M), el coeficiente de arrastre producido aumenta en una magnitud de 0,2 con un coeficiente de sustentacin mayor de 1.5, no variando tanto como se puede observar para un rgimen supersnico donde se

Mecnica de Vuelo

genera una aumenta de su resistencia del orden de 0,4 (desde transnico). En el grfico de familia de curvas (FIG.10), se puede apreciar este efecto.

Referencias
Janes All the Wolrds Aircraft 2008-2009. Abbott, H., Theory of Wing Section. Saez, R., Apuntes Curso Aerodinmica. 2010. .