Polar
Polar
Polar
NACIONAL
ESIME TICOMAN
Aerodinámica
Grupo: 5AM1
Figura 2.- Dibujo con las tres vistas del avión Douglas C-47 Sytrain.
Inicialmente se obtuvieron los siguientes datos del manual de mantenimiento T.O. 1C-47-2 del avión
tanto para el ala como para el estabilizador horizontal y se calcularon otros datos como la superficie,
alargamiento y conicidad que se utilizaran para el desarrollo del método Anderson.
Para el Ala:
Envergadura: 28.956m
𝐶𝑅 = 4.318m
𝐶𝑃 = 1.4224m
Perfil de raíz: NACA 2215
Perfil de punta: NACA 2206
(𝐶𝑅 + 𝐶𝑃 )𝑏 (4.318 + 1.4224)28.956
𝑆𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 = = = 83.109𝑚 2
2 2
𝑏2 28.9562
𝐴𝑙𝑎𝑟𝑔𝑎𝑚𝑖𝑒𝑛𝑡𝑜 (𝐴𝑅) = = = 10.2114
𝑆 82.109
𝐶𝑃 1.4224
𝐶𝑜𝑛𝑖𝑐𝑖𝑑𝑎𝑑 (𝜆) = = = 0.3294
𝐶𝑅 4.318
Con estos datos podemos llevar a cabo el método Anderson, para obtener las graficas
𝐶𝐿 𝑣𝑠 𝛼 𝑦 𝐶𝐷 𝑣𝑠 𝛼, tanto para el ala como para el estabilizador horizontal.
Método Anderson
Para el Ala.
1.8
CL vs Estación
1.6
1.4
0.8 Cl = 0.8
0.6 Cl = 1
0.4 Cl = 1.1
0.2 Cl = 1.19
0
0 0.2 0.4 0.6 0.8 1 1.2
Estación
Figura 3.- Método Anderson: gráfica CL vs estación para el ala.
0.09393 1
𝑎𝑎𝑙𝑎 = (0.996) = 0.080053
57.3(. 09393) °
1+( )
𝜋(10.2114)
Conociendo la pendiente de levantamiento del ala podemos utilizar la siguiente fórmula obtenida
también del reporte NACA 572 para calcular el valor de CL para cada ángulo de ataque.
Para conocer el valor del ángulo de cero levantamiento del ala utilizamos la siguiente ecuación:
𝛼𝐿=0 𝐴𝑙𝑎 = 𝛼𝐿=0 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 𝑝𝑒𝑟𝑓𝑖𝑙𝑒𝑠 + 𝐽𝜀
Se realiza una aproximación parabólica para conocer cómo se comporta la gráfica cuando deja de
ser lineal. De igual forma se obtienen valores promedios para 𝛼𝑆0 y 𝛥𝛼0 con los datos de ambos
perfiles.
𝑎𝑒 0.09393
𝛼𝛥0 𝑎𝑙𝑎 = ∗ 𝛥𝛼0 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 = ∗ 6.5 = 7.6733°
𝑎𝑎𝑙𝑎 0.08005
Para conocer el valor del CL del ala para este ángulo utilizamos la fórmula que ya usamos
anteriormente:
𝐶𝐿 𝑎𝑙𝑎 = 𝑎𝑎𝑙𝑎 (𝛼 + |𝛼𝐿=0 𝑎𝑙𝑎 |)
Sustituimos:
𝐶𝐿 𝑎𝑙𝑎 = 0.08005(6.2499 + |−1.9892|) = 0.6595
Conociendo estos valores podemos obtener la ecuación general de una parábola considerando que:
𝑉𝑒𝑟𝑡𝑖𝑐𝑒 (13.923, 1.19)
58.8764 = 4𝑃 (−0.5305)
58.8764 = −2.122𝑃
𝑃 = −27.7457
Una vez que conocemos el valor del foco podemos determinar la ecuación de la parábola que nos
describirá el comportamiento no lineal de la gráfica 𝐶𝐿 𝑣𝑠 𝛼.
(𝑥 − ℎ)2 = 4𝑃 (𝑦 − 𝑘)
𝑥2 27.846 61.7804
𝑦=− − 𝑥+
110.98 110.98 110.98
𝑦 = −0.009𝑥 2 + 0.2509𝑥 − 0.5566
Con esta ecuación obtendremos los valores de CL para cada valor de ángulo de ataque en la zona no
lineal de la gráfica es decir a partir de un ángulo de ataque de 6°, obteniendo la tabla que se muestra
a continuación.
𝐶𝐿2𝑎𝑙𝑎
𝐶𝐷𝑖 =
𝜋𝐴𝑅 𝑢
De las gráficas en el reporte NACA 572 leemos el valor de 𝑢 = 0.98 y obtenemos la siguiente tabla:
CL vs α CD vs α
1.60 0.30
1.20
0.80 0.20
CD
0.40
CL
0.00 0.10
-15 -5 5 15 25
-0.40
0.00
-0.80
α -15 -10 -5 0 α 5 10 15 20 25
Figura 4 y 5 .- Gráficas CL vs α y CD vs α del ala.
Para el estabilizador horizontal.
1.6
CL vs Estación
1.4
1.2
CL max Perfil
1
CL
Cl = 0.8
0.8
Cl = 1
0.6
Cl = 1.2
0.4
Cl = 1.3
0.2 Cl = 1.41
0
0 0.2 0.4 0.6 0.8 1 1.2
Estación
Figura 6.- Método Anderson: gráfica CL vs estación para el estabilizador horizontal.
De la gráfica podemos observar que el CLMax para el estabilizador horizontal será de 1.41.
Del reporte NACA 572 leemos los siguientes valores de gráficas para obtener la pendiente de
levantamiento del estabilizador horizontal:
𝐸 = 1.1
𝑓 = 0.998
1 𝑎𝑒
𝑎𝑝𝑒𝑟𝑓𝑖𝑙 = 0.104774 ° 𝑎𝑎𝑙𝑎 = 𝑓 57.3𝑎𝑒
1+( )
𝜋𝐴𝑅
𝑎𝑝𝑒𝑟𝑓𝑖𝑙 0.104774 1
𝑎𝑒 = = = 0.095249
𝐸 1.1 °
0.095249 1
𝑎𝐸𝐻 = (0.998) = 0.073570
57.3(. 095249) °
1+( )
𝜋(5.9481)
Conociendo la pendiente de levantamiento del estabilizador horizontal podemos utilizar la
siguiente fórmula obtenida también del reporte NACA 572 para calcular el valor de CL para cada
ángulo de ataque.
𝐶𝐿 𝐸𝐻 = 𝑎𝐸𝐻 (𝛼 + |𝛼𝐿=0 𝐸𝐻 |)
Para conocer el valor del ángulo de cero levantamiento del ala utilizamos la siguiente ecuación:
𝛼𝐿=0 𝐸𝐻 = 𝛼𝐿=0 𝑝𝑒𝑟𝑓𝑖𝑙 + 𝐽𝜀
Sabemos que el EH no tiene ningún tipo de torcimiento, por lo que el termino 𝐽𝜀 se vuelve cero,
entonces el valor del ángulo de cero levantamiento del estabilizador horizontal será igual al valor
del ángulo de cero levantamiento del perfil:
Se realiza una aproximación parabólica para conocer cómo se comporta la gráfica cuando deja de
ser lineal.
𝑎𝑒 0.095249
𝛼𝛥0 𝐸𝐻 = ∗ 𝛥𝛼0 = ∗ 10 = 12.9467°
𝑎𝐸𝐻 0.073570
Para conocer el valor del CL del estabilizador horizontal para este ángulo utilizamos la fórmula que
ya usamos anteriormente:
𝐶𝐿 𝐸𝐻 = 𝑎𝐸𝐻 (𝛼 + |𝛼𝐿=0 𝐸𝐻 |)
Sustituimos:
𝐶𝐿 𝐸𝐻 = 0.07357(6.4731 + |−2.262|) = 0.6426
Conociendo estos valores podemos obtener la ecuación general de una parábola considerando que:
𝑉𝑒𝑟𝑡𝑖𝑐𝑒 (19.4198, 1.41)
167.6170 = 4𝑃 (−0.7674)
167.6170 = −3.0696𝑃
𝑃 = −54.6054
Una vez que conocemos el valor del foco podemos determinar la ecuación de la parábola que nos
describirá el comportamiento no lineal de la gráfica 𝐶𝐿 𝑣𝑠 𝛼.
(𝑥 − ℎ)2 = 4𝑃 (𝑦 − 𝑘)
𝑥2 38.8396 69.1541
𝑦=− − 𝑥+
218.4216 218.4216 218.4216
𝑦 = −0.00457𝑥 2 + 0.1778𝑥 − 0.3166
Con esta ecuación obtendremos los valores de CL para cada valor de ángulo de ataque en la zona no
lineal de la gráfica, es decir a partir de un ángulo de ataque de 6°, obteniendo la tabla que se muestra
a continuación.
α (°) CL EH α (°) CL EH
-10 -0.5693 6 0.6078
-9 -0.4957 7 0.7068
-8 -0.4221 8 0.8171
-7 -0.3486 9 0.9184
-6 -0.2750 10 1.0107
-5 -0.2014 11 1.0940
-4 -0.1279 12 1.1683
-3 -0.0543 13 1.2336
-2 0.0193 14 1.2899
-1 0.0928 15 1.3372
0 0.1664 16 1.3755
1 0.2400 17 1.4048
2 0.3136 18 1.4090
3 0.3871 19 1.4163
4 0.4607 20 1.4120
5 0.5343 21 1.4084
Tabla 3.- Método Anderson: Valores de CL para cada ángulo de ataque del estabilizador horizontal.
Una vez que conocemos los valores de CL podemos calcular los valores de CD del estabilizador
horizontal de la siguiente manera:
𝐶𝐷 = 𝐶𝐷𝑂 + 𝐶𝐷𝑖
Donde:
𝐶𝐷𝑂 = 𝐶𝐷 𝑙𝑒í𝑑𝑜 𝑑𝑒 𝑙𝑎 𝑔𝑟á𝑓𝑖𝑐𝑎 𝑑𝑒𝑙 𝑝𝑒𝑟𝑓𝑖𝑙
𝐶𝐿2𝐸𝐻
𝐶𝐷𝑖 =
𝜋𝐴𝑅 𝑢
De las gráficas en el reporte NACA 572 leemos el valor de 𝑢 = 0.999 y obtenemos la siguiente tabla:
1.75 CL vs α 0.30 CD vs α
1.40
1.05 0.20
0.70
CD
CL
0.35
0.10
0.00
-15 -5 5 15 25
-0.35
0.00
-0.70
α -20 -10 0 α 10 20 30
Figura 7 y 8 .- Gráficas CL vs α y CD vs α del estabilizador horizontal.
Con estos valores y calculando el valor de la estela producida por el ala en la columna 7 podemos
obtener las primeras 15 columnas de datos que se muestran a continuación:
𝑆𝐵𝐴𝑅 19.3624
𝐶𝐷 ′𝐵𝐴𝑅 = 0.04 ( ) = 0.04 ( ) = 0.009
𝑆𝑎𝑙𝑎𝑟 83.1095
Figura 9 .- Dibujo con las dimensiones aprox. de la barquilla.
Figura 10 .- Dibujo del tren de aterrizaje principal señalando cada una de sus partes.
# de C'D elemento
Elemento Observaciones Superficie m2 CD elemento
elementos referido
1 2 Cilindro 0.00524 1.18 0.00015
2 4 Cilindro 0.04489 1.30 0.00281
3 4 Cilindro 0.03445 1.23 0.00203
4 4 Cilindro 0.01122 1.17 0.00063
5 2 Llanta 0.23587 1.16 0.00658
6 8 Sec. Rectangular 0.01123 1.25 0.00135
7 2 Cilindro 0.02806 1.27 0.00086
8 4 Cilindro 0.01348 1.15 0.00075
9 2 Cilindro 0.07489 1.28 0.00231
10 2 Cilindro 0.08987 1.34 0.00290
C'D TANariz = 0.02036
Tabla 5.- CD de cada uno de los elementos del tren de aterrizaje principal.
Figura 11 .- Dibujo del tren de aterrizaje de cola señalando cada una de sus partes.
# de C'D elemento
Elemento Observaciones Superficie m2 CD elemento
elementos referido
1 1 Llanta 0.01960 1.19 0.00028
2 1 Sec. Rectangular 0.01430 1.28 0.00022
Fuselaje
Para obtener el CD parasita del fuselaje se utilizó la gráfica de la página 512 del libro Theory
of Flight”. R.Von Mises, para un fuselaje de sección transversal circular.
CD vs α
0.10
y = 0.0003x2 + 0.0004x + 0.0418
0.08
0.05
CD
0.03
0.00
-10 -5 0 α 5 10 15
Con los valores obtenidos anteriormente podemos completar las columnas 16-21 del método.
Para calcular el incremento en el valor del coeficiente de levantamiento referido a la superficie alar,
utilizamos la fórmula para aletas simples y “Split”:
Se obtuvieron las dimensiones aproximadas en base al siguiente dibujo a escala, estimando como
una cuerda máxima de 3.72m y una cuerda mínima de 1.43m para la superficie con aletas, con una
semienvergadura de 7.3914m.
Sfi Sfd
(3.72 + 1.43)(7.3914)
𝑆𝑓𝑖 = 𝑆𝑓𝑑 = = 19.032𝑚2
2
𝑆𝑓 = 𝑆𝑓𝑖 + 𝑆𝑓𝑑
𝑆𝑓 38.065
𝛥𝐶𝐿′ = 0.9 ( ) = 0.9 ( ) = 0.4122
𝑆𝐴 83.109
Sabemos que la cuerda de la aleta es de 0.785m, la cuerda aerodinámica media del ala es: 3.1136m
y el ángulo de deflexión de la aleta (δ) es de 60° como se indicó en el archivo, por lo tanto:
Con esto, podemos completar las ultimas columnas de las tablas y construir así las polares de
nuestro avión.
La tabla completa con los datos obtenidos es la siguiente:
αF (°) αA (°) CLA CDoA CDiA CDA δ (°) αEH (°) CLEH CL'EH CDoEH CDiEH CDEH CD'EH
-4 -2 -0.0009 0.0075 0.0000 0.0075 -0.0023 -4.00 -0.1277 -0.0128 0.0150 0.0009 0.0159 0.0016
-3 -1 0.0792 0.0065 0.0002 0.0067 0.2065 -3.21 -0.0695 -0.0070 0.0130 0.0003 0.0133 0.0013
-2 0 0.1592 0.0055 0.0008 0.0063 0.4153 -2.42 -0.0113 -0.0011 0.0100 0.0000 0.0100 0.0010
-1 1 0.2393 0.0060 0.0018 0.0078 0.6241 -1.62 0.0469 0.0047 0.0100 0.0001 0.0101 0.0010
0 2 0.3193 0.0065 0.0033 0.0098 0.8329 -0.83 0.1051 0.0105 0.0090 0.0006 0.0096 0.0010
1 3 0.3994 0.0070 0.0051 0.0121 1.0416 -0.04 0.1634 0.0164 0.0080 0.0014 0.0094 0.0009
2 4 0.4795 0.0075 0.0073 0.0148 1.2504 0.75 0.2216 0.0222 0.0080 0.0026 0.0106 0.0011
3 5 0.5595 0.0085 0.0100 0.0185 1.4592 1.54 0.2798 0.0280 0.0070 0.0042 0.0112 0.0011
4 6 0.6396 0.0105 0.0131 0.0236 1.6680 2.33 0.3380 0.0339 0.0070 0.0061 0.0131 0.0013
5 7 0.7587 0.0120 0.0184 0.0304 1.9787 3.02 0.3887 0.0390 0.0060 0.0081 0.0141 0.0014
6 8 0.8746 0.0140 0.0245 0.0385 2.2809 3.72 0.4400 0.0441 0.0060 0.0104 0.0164 0.0016
7 9 0.9725 0.0185 0.0302 0.0487 2.5363 4.46 0.4948 0.0496 0.0050 0.0131 0.0181 0.0018
8 10 1.0524 0.0330 0.0354 0.0684 2.7446 5.26 0.5531 0.0554 0.0060 0.0164 0.0224 0.0022
9 11 1.1143 0.0545 0.0397 0.0942 2.9061 6.09 0.6148 0.0616 0.0060 0.0202 0.0262 0.0026
10 12 1.1582 0.0695 0.0429 0.1124 3.0206 6.98 0.7044 0.0706 0.0070 0.0266 0.0336 0.0034
11 13 1.1841 0.0865 0.0448 0.1313 3.0881 7.91 0.8077 0.0810 0.0080 0.0350 0.0430 0.0043
12 14 1.1920 0.0985 0.0454 0.1439 3.1087 8.89 0.9078 0.0910 0.0090 0.0441 0.0531 0.0053
13 15 1.1819 0.1105 0.0447 0.1552 3.0824 9.92 1.0034 0.1006 0.0110 0.0539 0.0649 0.0065
14 16 1.1538 0.1235 0.0426 0.1661 3.0091 10.99 1.0933 0.1096 0.0120 0.0640 0.0760 0.0076
15 17 1.1077 0.1370 0.0392 0.1762 2.8889 12.11 1.1760 0.1179 0.0140 0.0741 0.0881 0.0088
16 18 1.0436 0.1525 0.0348 0.1873 2.7217 13.28 1.2502 0.1253 0.0160 0.0837 0.0997 0.0100
17 19 0.9615 0.1710 0.0296 0.2006 2.5076 14.49 1.3143 0.1317 0.0180 0.0925 0.1105 0.0111
18 20 0.8614 0.1890 0.0237 0.2127 2.2465 15.75 1.3669 0.1370 0.0210 0.1001 0.1211 0.0121
19 21 0.7433 0.2050 0.0177 0.2227 1.9385 17.06 1.4063 0.1410 0.0230 0.1059 0.1289 0.0129
20 22 0.6072 0.2200 0.0118 0.2318 1.5836 18.42 1.4309 0.1434 0.0270 0.1097 0.1367 0.0137
Tabla 8a.- Valores obtenidos por el método Neo POLAR Signed IV.
CDTS/A y CDTC/A y
CLTS/A CD'EV CD'BAR CD'TA CD'FUS 1.05CD'PAR ΔCL' ΔCD' CLTC/A TA TA CDTC/TA CDTC/A
-0.0137 0.0005 0.0093 0.0212 0.0450 0.0798 0.4122 0.0872 0.3985 0.0677 0.1761 0.0889 0.1549
0.0722 0.0005 0.0093 0.0212 0.0433 0.0780 0.4122 0.0872 0.4844 0.0648 0.1732 0.0860 0.1520
0.1581 0.0005 0.0093 0.0212 0.0422 0.0768 0.4122 0.0872 0.5703 0.0630 0.1714 0.0842 0.1502
0.2440 0.0005 0.0093 0.0212 0.0417 0.0763 0.4122 0.0872 0.6562 0.0640 0.1724 0.0852 0.1512
0.3299 0.0005 0.0093 0.0212 0.0418 0.0764 0.4122 0.0872 0.7421 0.0659 0.1744 0.0871 0.1532
0.4158 0.0005 0.0093 0.0212 0.0425 0.0772 0.4122 0.0872 0.8280 0.0690 0.1774 0.0902 0.1562
0.5017 0.0005 0.0093 0.0212 0.0438 0.0785 0.4122 0.0872 0.9139 0.0732 0.1816 0.0944 0.1604
0.5875 0.0005 0.0093 0.0212 0.0457 0.0805 0.4122 0.0872 0.9997 0.0789 0.1874 0.1001 0.1662
0.6734 0.0005 0.0093 0.0212 0.0482 0.0831 0.4122 0.0872 1.0856 0.0868 0.1952 0.1080 0.1740
0.7977 0.0005 0.0093 0.0212 0.0513 0.0864 0.4122 0.0872 1.2099 0.0970 0.2054 0.1182 0.1842
0.9187 0.0005 0.0093 0.0212 0.0550 0.0903 0.4122 0.0872 1.3309 0.1092 0.2176 0.1304 0.1964
1.0221 0.0005 0.0093 0.0212 0.0593 0.0948 0.4122 0.0872 1.4343 0.1241 0.2326 0.1453 0.2114
1.1078 0.0005 0.0093 0.0212 0.0642 0.0999 0.4122 0.0872 1.5200 0.1494 0.2578 0.1706 0.2366
1.1759 0.0005 0.0093 0.0212 0.0697 0.1057 0.4122 0.0872 1.5881 0.1813 0.2898 0.2025 0.2686
1.2288 0.0005 0.0093 0.0212 0.0758 0.1121 0.4122 0.0872 1.6410 0.2067 0.3151 0.2279 0.2939
1.2651 0.0005 0.0093 0.0212 0.0825 0.1192 0.4122 0.0872 1.6773 0.2336 0.3420 0.2548 0.3208
1.2830 0.0005 0.0093 0.0212 0.0898 0.1268 0.4122 0.0872 1.6952 0.2549 0.3633 0.2761 0.3421
1.2825 0.0005 0.0093 0.0212 0.0977 0.1351 0.4122 0.0872 1.6947 0.2756 0.3840 0.2968 0.3628
1.2634 0.0005 0.0093 0.0212 0.1062 0.1440 0.4122 0.0872 1.6756 0.2965 0.4049 0.3177 0.3837
1.2256 0.0005 0.0093 0.0212 0.1153 0.1536 0.4122 0.0872 1.6378 0.3175 0.4259 0.3387 0.4047
1.1689 0.0005 0.0093 0.0212 0.1250 0.1638 0.4122 0.0872 1.5811 0.3399 0.4483 0.3611 0.4271
1.0932 0.0005 0.0093 0.0212 0.1353 0.1746 0.4122 0.0872 1.5054 0.3650 0.4734 0.3862 0.4522
0.9984 0.0005 0.0093 0.0212 0.1462 0.1860 0.4122 0.0872 1.4106 0.3897 0.4981 0.4109 0.4769
0.8843 0.0005 0.0093 0.0212 0.1577 0.1981 0.4122 0.0872 1.2965 0.4125 0.5209 0.4337 0.4997
0.7506 0.0005 0.0093 0.0212 0.1698 0.2108 0.4122 0.0872 1.1628 0.4351 0.5435 0.4563 0.5223
Tabla 8b.- Valores obtenidos por el método Neo POLAR Signed IV.
Polares obtenidas:
Polar del avión “limpio”. (aletas y tren de aterrizaje retraídos).
CL vs CD
Avion limpio
1.40
1.20
1.00
0.80
CL
0.60
0.40
0.20
0.00
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40
CD
-0.20
CL vs CD
Avión sucio
1.80
1.60
1.40
1.20
1.00
CL
0.80
0.60
0.40
0.20
0.00
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50
CD
Polar del avión sólo con tren de aterrizaje extendido.
CL vs CD
Solo Tren extendido
1.40
1.20
1.00
0.80
CL
0.60
0.40
0.20
0.00
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40
-0.20 CD
CL vs CD
Solo aletas extendidas
1.80
1.60
1.40
1.20
1.00
CL
0.80
0.60
0.40
0.20
0.00
0.00 0.05 0.10 0.15 0.20 C 0.25 0.30 0.35 0.40 0.45 0.50
D