Diseño Conceptual de Aeronaves
Diseño Conceptual de Aeronaves
1
Documento I
Índice
3
ÍNDICE GENERAL
Índice general
I Índice 3
Índice general 5
Índice de figuras 12
Índice de tablas 18
II Memoria 21
1. Objeto 23
2. Alcance 25
3.1.2. RDS-win . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.1.4. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
4. Normas y referencias 31
4.3. Bibliografía . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
5. Definiciones y abreviaturas 37
5.1. Definiciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
5
ÍNDICE GENERAL
6. Requisitos de diseño 47
7. Análisis de soluciones 49
7.4.1. Dimensionamiento . . . . . . . . . . . . . . . . . . . . . . . . 55
7.5.1.1. Fuselaje . . . . . . . . . . . . . . . . . . . . . . . . . 57
7.5.1.2. Ala . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.5.1.3. Cola . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
6
ÍNDICE GENERAL
8. Resultados finales 93
9. Planificación 97
V Presupuesto 113
VI Anexos 119
7
ÍNDICE GENERAL
8
ÍNDICE GENERAL
9
ÍNDICE GENERAL
Ñ.12.Mobiliario . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
10
ÍNDICE GENERAL
11
ÍNDICE DE FIGURAS
Índice de figuras
15. Dibujo aclaratorio sobre el ángulo de flecha Λ para un ala en planta. [7] 74
17. Ilustración de la diferencia entre dos polares, una con perfiles con cur-
vatura y otra sin curvatura. (Uncambered=sin curvatura/Cambered=con
curvatura). [1] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
12
ÍNDICE DE FIGURAS
32. Representación del ala trapezoidal de referencia, con todos los pará-
metros que la definen. [6] . . . . . . . . . . . . . . . . . . . . . . . . . 144
13
ÍNDICE DE FIGURAS
39. Variación del parámetroKf0 en función del ángulo de deflexión δf [5]. 162
41. Obtención del valor Kf0 para el caso de estudio [5]. . . . . . . . . . . . 163
dα
42. Variación del parámetro dδ0fl en función del ángulo de deflexión del
flap, y de la relación de cuerdas cf /c, [5] . . . . . . . . . . . . . . . . 163
dα0l
43. Obtención del valor dδf
para el caso de estudio [5]. . . . . . . . . . . 164
47. Variación del ángulo de entrada en pérdida para los flaps tipo Fowler,
double y triple slotted (αstall )c0 /c en función del ángulo de deflexión, [5]167
14
ÍNDICE DE FIGURAS
67. Variación del parámetro ∆1 en función del ratio de cuerdas, para flaps
tipo slotted plain y Split [6]. . . . . . . . . . . . . . . . . . . . . . . . 203
15
ÍNDICE DE FIGURAS
16
ÍNDICE DE FIGURAS
17
ÍNDICE DE TABLAS
Índice de tablas
10. Tabla para introducir los datos de la misión de vuelo, Hoja2 del Excel
TabladeTodoslosdatos.xlsx del programa de diseño . . . . . . . . . . . 122
13. Tabla con los valores típicos de los coeficientes del cálculo de cola
horizontal y vertical [1]. . . . . . . . . . . . . . . . . . . . . . . . . . 148
16. Variación del parámetro ∆1 en función del espesor relativo y del ratio
de cuerdas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 201
17. Variación del parámetro ∆2 en función del espesor relativo y del ratio
de cuerdas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 202
18
ÍNDICE DE TABLAS
18. Valores típicos según aeronave, de los factores de carga límites, [1].
Tabla incorporada en el Excel TablaEstructuras.xlsx. . . . . . . . . . . 209
19
Documento II
Memoria
21
1. OBJETO
1. Objeto
Es por ello que si se quiere hacer un diseño conceptual sin tener un programa
informático a mano que pueda ayudar al diseñador, la tarea en cuestión puede
hacerse cuesta arriba.
23
2. ALCANCE
2. Alcance
Por otra parte, las herramientas que se han propuesto y usado a lo largo del
trabajo son habituales en los procedimientos y procesos de diseño y de cálculo en la
ingeniería: Excel y MATLAB.
25
3. ANTECEDENTES Y JUSTIFICACIÓN DEL PROYECTO
Si bien, ninguno de estos programas tienen como idea global la del diseño de una
aeronave, son de gran utilidad para el estudio de cada una de las partes mencionadas.
Si se quiere un programa informático que de verdad englobe todas las partes, o al
menos las imprescindibles en el diseño de una aeronave, se tiene en el mercado cierto
número de ellos. A continuación se mencionarán algunos de estos programas, a fin
de tener una idea de los programas existentes y de comparar los mismos al programa
que se pretende desarrollar en el presente proyecto.
Este es un programa basado en el libro del autor Jack D. Mattingly [4], siendo
su libro de diseño de aeronaves, un libro centrado principalmente en el estudio de
motores. Dicho programa es más bien de diseño de motorización de aeronaves.
3.1.2. RDS-win
Este programa informático esta basado en el libro del autor Raymer [1], ya
comentado anteriormente, la empresa que comercializa el software es American Ins-
titute of Aeronautics and Astronautics (AIAA). A diferencia del anterior caso, este
programa si sirve para el diseño global de una aeronave, ya que en un mismo pro-
27
3. ANTECEDENTES Y JUSTIFICACIÓN DEL PROYECTO
grama engloba aerodinámica, estructuras, una parte de diseño tipo CAD, estudios
de motores, de pesos, estabilidad etc...
Es por tanto un programa que se centra en la idea del presente proyecto, y por
tanto una herramienta a tener en cuenta en el estudio.
28
3. ANTECEDENTES Y JUSTIFICACIÓN DEL PROYECTO
29
3. ANTECEDENTES Y JUSTIFICACIÓN DEL PROYECTO
3.1.4. Conclusiones
Tanto el RDS-win como el AAA son programas que se deben tener en cuenta
como herramientas en el diseño de aeronaves, sin embargo el coste de los mismos es
muy elevado frente al uso que se les puede dar por parte de estudiantes de ingeniería.
Por otro lado los programas comentados, sobre todo el AAA program, tienen una
gran complejidad de uso. De este último se imparten clases a fin de poder utilizar
de forma correcta el programa, con el sobrecoste que ello conlleva para el usuario.
Se puede afirmar por tanto que la herramienta informática que se pretende desa-
rrollar, puede llegar a ser una alternativa para estudiantes o usuarios que quieran
un programa, aunque mucho mas simple, mas fácil de utilizar y sin coste alguno.
30
4. NORMAS Y REFERENCIAS
4. Normas y referencias
31
4. NORMAS Y REFERENCIAS
32
4.3. Bibliografía
[2] Dr. Chuan-Tau Edward Lan, Dr. Jan Roskam (2008) Aiplane Aerodyna-
mics and Performance. Lawrence, Kansas, DAR corporation, 4a Ed.
[3] E. Torenbeek (1982) Synthesis of Subsonic Airplane Design. Delft, The Net-
herlands, Delft Univesity Press.
[13] http://aircraftdesign.nuaa.edu.cn/pd-2004/DownLoad/AEDsys%20User%
20Guide.pdf, Visitado por ultima vez el día 01/08/2019.
[14] https://www.fzt.haw-hamburg.de/pers/Scholz/HOOU/AircraftDesign_7_
WingDesign.pdf, Visitado por ultima vez el día 03/08/2019.
[16] https://web.mit.edu/drela/Public/web/xfoil/
33
[17] J.Messeguer, A.Sanz (2011) Aerodinámica Básica.. Madrid, España, Garce-
ta, 2a Ed.
[19] Callaghan J.G. Aerodynamic Prediction Methods for Aircraft al Low speeds
with Mechanical High Lift Devices. Paper No. 2 in Prediction Methods for Air-
craft Aerodynamics Characteristics, AGARD-LS-67, (1974).
[20] Lowry, J.G. and Polhamus, E.C. A method for predicting Lift Increments
due to flap Deflection al Low Angles of Attack in Incompressible Flow. NACA
TN3911, (1957).
[23] Dr. Jan Roskam (2002) Aiplane design Part VI. Lawrence, Kansas, Design
Analysis and Research (DAR) Corporation.
[24] http://www.utias.utoronto.ca/wp-content/uploads/2017/09/
Potter-Strut_Braced_Wing.pdf University of Toronto Institute for Ae-
rospace Studies
[25] Dr. Jan Roskam (1984) Methods for Estimating Drag Polars of Subsonic
Airplanes. Lawrence, Kansas, Design Analysis and Research (DAR) Corpora-
tion, 4a Ed.
[26] Parkinson, John B., Roland E. Olson, and Rufus O. House (1939)
Hydrodynamic and Aerodynamic Tests of aFamily of Models of Seaplane Floats
with Varying Angles of Dead Rise.NACA TN-716.
[27] Young AD. (1947) The Aerodynamic Characteristics of Flaps.R.M. No. 2622,
British A. R. C.
[33] http://airfoiltools.com/search/index
34
[34] https://m-selig.ae.illinois.edu/ads/coord_database.html
[35] https://sourceforge.net/projects/xflr5/files/
35
5. DEFINICIONES Y ABREVIATURAS
5. Definiciones y abreviaturas
5.1. Definiciones
Carga alar: Se define la carga alar como un parámetro aerodinámico que re-
laciona el peso total de la aeronave con respecto a la superficie de referencia de la
aeronave o superficie de sustentación.
Crucero: Etapa de vuelo recto y nivelado, en el cual, las cuatro fuerzas prin-
cipales resistencia/empuje y sustentación/peso se encuentran en equilibrio. En esta
etapa se optimiza el combustible gastado para una distancia recorrida, o para una
velocidad dada.
Diedro del ala: Ángulo que mide la elevación sobre la horizontal de las alas.
37
5. DEFINICIONES Y ABREVIATURAS
motor de una aeronave, ya sea por desplazar el aire hacia atrás mediante hélices, o
por expulsar el aire mediante un motor a reacción. En general en aeronaves tipo jet
o aerorreactores se habla de empuje, mientras que en aeronaves tipo hélice se utiliza
la potencia como medida de lo que nos proporciona el motor.
Estrechamiento del ala: Relación entre la cuerda del perfil alar en la punta
del ala dividido por la cuerda del perfil alar del encastre. Es decir mide la variación
de la longitud de cuerda a lo largo del ala.
Factor de carga: Se define como la relación que existe entre la fuerza total en
dirección vertical que actúa sobre el avión y el peso del mismo.
Flecha del ala: Ángulo que forman el ala con respecto al eje transversal del
avión.
Mecánica del vuelo: Área de la ingeniería aeronáutica que estudia las actuacio-
nes de las aeronaves, así como la estabilidad y controlabilidad tanto estáticas como
dinámicas. En el presente proyecto al hablar de mecánica de vuelo se relacionará
solamente con las actuaciones.
Perfil: Forma del área transversal del ala, que al desplazarse a través del aire es
capaz de crear a su alrededor una distribución de presiones que genere sustentación.
Peso en vacío: El peso en vacío o vacío operativo es el peso del avión sin carga
útil y sin combustible. Por otro lado, incluye todos los líquidos del avión necesarios
38
5. DEFINICIONES Y ABREVIATURAS
para su operación.
Torsión del ala: Variación, a lo largo de la envergadura del ala, del ángulo
formado por una línea del perfil con respecto a un plano de referencia normal al de
simetría.
39
5. DEFINICIONES Y ABREVIATURAS
40
5. DEFINICIONES Y ABREVIATURAS
41
5. DEFINICIONES Y ABREVIATURAS
42
5. DEFINICIONES Y ABREVIATURAS
43
5. DEFINICIONES Y ABREVIATURAS
44
5. DEFINICIONES Y ABREVIATURAS
45
5. DEFINICIONES Y ABREVIATURAS
46
6. REQUISITOS DE DISEÑO
6. Requisitos de diseño
47
7. ANÁLISIS DE SOLUCIONES
7. Análisis de soluciones
Tal y como se comenta en el libro de referencia del autor Raymer [1], el diseño de
una aeronave se puede establecer de forma general en tres niveles: diseño conceptual,
diseño preliminar y diseño detallado. En el primer nivel se plantean aspectos básicos
y fundamentales de la aeronave como el peso y tamaño de la aeronave, se llega
a realizar un pequeño sketch de la aeronave, y se consiguen obtener ciertos datos
sobre las actuaciones de la aeronave. Con estos datos es posible analizar y simular
el comportamiento de la nave en distintas actuaciones.
El segundo, se trata del diseño preliminar. En este nivel del diseño, los cambios
que se realicen en la aeronave deben ser pequeños, y no alterar en gran medida al
diseño elegido en el primer nivel. Se realizan estudios más específicos de la aeronave,
tales como estudios estructurales, de estabilidad, de actuaciones (de una forma más
especifica que el caso anterior) y se diseñan nuevas partes de la aeronave que no se
habían realizado con anterioridad, véase el caso del tren de aterrizaje, por ejemplo.
Por último, el tercer nivel corresponde al diseño detallado, en este caso el dise-
ño de la aeronave ya no se debe tocar, se trata ahora de detallar las partes de la
aeronave, es decir, en el caso del ala de la aeronave, teniendo ya todas las caracterís-
ticas de diseño tamaño y aerodinámica, en esta fase se deberá seccionar y establecer
sus partes, larguerillos, costillas, número de tanques de combustible en dichas alas,
integración de motores (si los hubiera), conductos para el flujo de combustible, com-
ponentes eléctricos etc..., se trata por tanto, de la última fase de diseño antes de
comenzar la producción de la aeronave.
Se debe tener en cuenta que se busca realizar una herramienta de diseño siguiendo
la estructura de la asignatura de aeronaves ya comentada, de esta forma las divisiones
están ordenadas de la misma forma que se realiza en dicha asignatura y se trata
de una estructura similar a la que se sigue en los principales libros de diseño de
aeronaves, y que por tanto se entiende que es la mejor forma posible de llevar un
diseño.
49
7. ANÁLISIS DE SOLUCIONES
TablaGeometría
.xlsx
T/W T/W
W/S Geometría.m TablaAero W/S
W0 .xlsx W0
1ª aproximación n+1ª aproximación
Aero.m
TablaEstructuras
.xlsx
FuncionCalculos
Aerodinamicos.
m
Estructuras.m
TablaDiseñoAvanzado
.xlsx
FuncionCalculoPeso
EnVacio.m
MecVuelo.m T/W
W/S
NO
¿La aproximación n
Parte2.m coincide con la n+1?
SI
Por último, los bloques de color verde, corresponden a los archivos de Excel, que
son archivos de introducción de datos y visualización de resultados.
50
7. ANÁLISIS DE SOLUCIONES
Se repetirán tantas veces este proceso, hasta que los valores obtenidos en la
aproximación n+1 sea igual a la obtenida en la aproximación n, entendiendo n
como el número de veces que se ha repetido el programa.
Una vez explicado de forma general el esquema interno del programa (figura 4)
y su interpretación, se puede explicar el programa en sí, atendiendo a dos partes;
por un lado se encuentra la parte que se conoce como diseño básico y por otro el
diseño avanzado.
Una vez explicado de forma general el esquema del programa mediante un dia-
grama, es necesario realizar una descripción conceptual de dicho programa, a fin de
comprender la disposición de los bloques.
Téngase en cuenta que los nombres de diseño básico y diseño avanzado son una
forma de explicar el desarrollo seguido para el diseño de una aeronave, tal y como
se realiza en la asignatura de aeronaves, y tal como siguen algunos de los libros mas
importantes de diseño de aeronaves, véase el caso de [7], [4] o [6].
51
7. ANÁLISIS DE SOLUCIONES
Una vez realizada la parte primera (diseño básico), ya se tiene una idea general
de la aeronave a diseñar, por tanto en este momento el usuario, se encuentra en
disposición de especificar de una forma más detallada ciertas características de la
aeronave, que en la primera parte no era necesario especificar, puesto que la mayoría
de cálculos o datos tratados están basados en datos estadísticos o empíricos de
aeronaves ya existentes.
RFP
En esta parte, puesto que se utilizan ecuaciones de mecánica del vuelo, será
necesario tener una comprensión más amplia de la aeronave, como por ejemplo la
aerodinámica, el motor que utiliza la aeronave (tipo jet o hélice), y otros datos más
específicos de la misión de vuelo (distancia de despegue, de aterrizaje, altura de
vuelo, etc), que tal como se pueden observar en el anterior diagrama expuesto, ver
figura 4, corresponden a los bloques que se encuentran entre MecVuelo.m y Parte1.m.
52
7. ANÁLISIS DE SOLUCIONES
Segundo dimensionamiento
Finalmente una vez realizado el RFP, y teniendo una comprensión sobre cómo
será la aeronave a diseñar, se podrá realizar un segundo dimensionamiento, ahora sí
utilizando cálculos propios y no utilizando datos estadísticos que se habían utiliza-
do en la primera parte. A pesar de esto, el usuario debe tener claro que el primer
dimensionamiento es bastante válido para tener un primer diseño conceptual de la
aeronave. Así se puede observar en libros como Raymer [1], donde el dimensiona-
miento se queda solamente en la primera parte expuesta anteriormente.
En el siguiente apartado se explicará de una forma mas detallada cada uno de los
módulos que componen el programa, explicando los cálculos que se realizan en cada
uno de los desarrollos seguidos para su obtención; se debe tener en cuenta que a lo
largo de este apartado se intentará hacer referencia en todo momento al programa de
diseño, a fin de que se tenga una idea del posicionamiento de los módulos y tener una
idea mas precisa de como utilizar el programa; de todas formas, como ya se comentó
anteriormente, se deja en el Anexo Q, una guía de como utilizar el programa.
Por otra parte, definida la propulsión de la aeronave, también será necesario es-
tablecer de una forma no muy detallada ciertos valores característicos de la aeronave
a diseñar; puesto que el usuario en esta parte del diseño no tiene una idea fija de
la aeronave, en el programa informático diseñado se deja, a fin de tener una ayuda,
una tabla de aeronaves similares a las existentes en el mercado.
53
7. ANÁLISIS DE SOLUCIONES
54
7. ANÁLISIS DE SOLUCIONES
Por otra parte, se debe tener en cuenta que para esta primera parte del diseño se
aportan para todas las aeronaves valores de referencia, que se han obtenido de libros
de referencia [1] y [2], a fin de que el usuario pueda realizar una primera aproximación
sencilla de la aeronave a diseñar; por tanto a lo largo de todo el desarrollo que se
realizará a continuación, se debe tener en cuenta que el usuario tiene ya en primera
instancia todos los datos, aunque puede cambiarlos, y sencillamente solo tendrá que
elegir el tipo de aeronave que quiera y definir su misión de vuelo (se explicará en la
parte de dimensionamiento).
7.4.1. Dimensionamiento
Se establece en esta parte, como objetivo principal, el cálculo del peso máximo
al despegue de la aeronave y el peso del combustible necesario para la misión de
vuelo. Es en esta parte del diseño donde el usuario tendrá que comenzar a delimitar
su diseño.
Una vez definido esto, se plantea la misión de vuelo; puesto que la herramienta
realizada para el diseño de aeronaves abarca distintos tipos de aeronaves, es necesario
delimitar el tipo de misión en dos partes, por un lado una parte fija, que se supondrá
común para cualquier aeronave, y por otro una parte que variará en función del tipo
de misión que quiera plantear el usuario.
Las partes que componen la misión de vuelo en el diseño básico, al igual que
la forma en la que se aportan en la herramienta informática, queda definida en el
anexo A del presente proyecto.
En esta parte del módulo de MATLAB Parte1.m, donde se podrá obtener una
tabla que muestre tanto la carga alar como la relación empuje peso; dos variable de
gran importancia en el diseño de aviones.
55
7. ANÁLISIS DE SOLUCIONES
(ver apartado 7.2.2), a fin de afinar la obtención del empuje necesario, y de este
modo poder elegir el motor que necesitará el avión de una forma más eficiente.
En este caso se tienen dos vertientes, en función del tipo de propulsión. Es decir,
en el caso de que la aeronave tenga propulsión tipo jet se hablará de la relación
empuje-peso (T /W ) y en el caso de tener un motor propulsado por hélice el término
utilizado será Potencia-peso (P/W ); cualquiera de estos dos términos definen bási-
camente lo mismo, sin embargo, de este modo le será más fácil al usuario comprender
el tamaño del motor, puesto que para motores con hélices la mayoría de fabricantes
hablan de potencia y para motores a reacción se utiliza el empuje.
Una vez finalizado el diseño básico, se tiene en estos momentos una aproximación
de la aeronave, sin embargo, en ningún caso se ha especificado la geometría de la
aeronave ni se han estudiado ciertos aspectos necesarios en un diseño conceptual.
Se tiene por tanto en esta parte del diseño, los valores que se han incorporado,
ver tabla 1, y los valores obtenido en el módulo Parte1.m, es decir, una aproximación
del peso total de la aeronave, la relación empuje peso y la carga alar, lo que supone
tener el empuje y la superficie alar.
Una vez realizada la primera parte, donde se han obtenido datos comentados, es
posible ahora tener una estimación no excesivamente precisa, pero al menos válida,
del tamaño de la aeronave a diseñar.
56
7. ANÁLISIS DE SOLUCIONES
Se debe tener en cuenta que en este momento el usuario deberá antes de ejecutar
el módulo Geometría incorporar ciertos valores en el Excel TablaGeometria.xlsx.
Por otra parte, todos los datos y ecuaciones incorporadas a la hoja del Excel
del módulo TablaGeometría.xlsx, se encuentran explicadas detalladamente en el
anexo del presente proyecto, Anexo Q.
El diseño gráfico del fuselaje, las alas, colas y resto de componentes de la aeronave
no son objeto del trabajo, si bien se pueden mencionar ciertos programas que pueden
ser de utilidad al usuario a fin de completar y tener un sketch o boceto lo más
adecuado posible a la realidad; dichos programas pueden ser OpenVsp de código
libre, o el programa de diseño gráfico Catia, que si bien es de pago permite realizar
la mayoría de funciones mediante la versión de estudiante.
De este modo a partir de los datos que se obtengan de esta parte del programa se
podrá realizar el diseño gráfico de la aeronave; por otra parte, muchos de los valores
de los parámetros que se aportarán a continuación serán necesarios más adelante en
el estudio de de pesos y para el módulo de aerodinámica
7.5.1.1 Fuselaje
57
7. ANÁLISIS DE SOLUCIONES
que supone obtener la longitud del fuselaje de forma sencilla. Dicha ecuación se
plantea a continuación:
L = aW0 C (1)
Donde a y C son dos parámetros que dependen del tipo de aeronave elegida, y que
aparecen recogidos en la hoja 1 del Excel TablaGeometria.xlsx, L será la longitud del
fuselaje calculada y W0 es el peso total que se calculó en la Parte1.m. Los valores de a
como de C, se obtienen directamente de la tabla 2 donde se establecen sus valores en
función del tipo de aeronave, y que también se encuentra integrada dentro del Excel
del programa. Téngase en cuenta que la columna que aparece en dicha tabla vacía
corresponde al valor de la longitud de fuselaje calculado y que una vez ejecutado el
módulo de la parte1.m aparecerá el valor en la fila correspondiente a tipo de aeronave
elegida.
Tabla 2: Datos para el cálculo de la longitud del fuselaje según tipo de aeronave
elegida, Hoja1 del Excel TablaGeometria.xlsx del programa de diseño
58
7. ANÁLISIS DE SOLUCIONES
Otro factor importante será la propia estructura del fuselaje, como ya es sabido
la mayoría de fuselajes de aeronaves se encuentran presurizados, es por tanto que
la mayoría de los fuselajes presentan secciones circulares o cónicas, a fin de repartir
de forma más uniforme las fuerzas generadas por las diferencias de presión entre el
interior de la cabina y la parte exterior donde navega la aeronave; por otra parte las
fuerzas de sustentación y de resistencia también afectan al propio fuselaje, por lo
que también la estructura del mismo debe ser la adecuada para ser estructuralmente
resistente y evitar su rotura catastrófica.
Por último, en la geometría del fuselaje se puede definir lo que se conoce como
fineness ratio o relación de finura del fuselaje, la cual es la relación entre la longitud
del fuselaje y su diámetro máximo.
De este modo se puede determinar el diámetro máximo del fuselaje, puesto que
según [1], teóricamente, para un volumen interno fijo, la resistencia en régimen subsó-
nico se minimiza con una relación de finura de aproximadamente 3.0, encontrándose
la mayoría de los aviones convencionales en estos valores.
Se puede utilizar por tanto una relación de finura del fuselaje obtenida de datos
históricos de otras aeronaves, junto con la estimación de la longitud, para desarrollar
el diseño inicial del fuselaje, siempre teniendo en cuenta todo lo comentado en el
apartado.
Por otra parte, una vez que se tenga en mente el ancho máximo y longitud del
fuselaje, será necesario establecer el mismo de una manera específica, puesto que
para el módulo de aerodinámica serán necesario calcular la superficie mojada del
fuselaje (Swetted , parámetro que será explicado mas adelante) a fin de calcular la
resistencia generada por la estructura del fuselaje del avión.
Por ello, para el cálculo de la superficie, y también como ayuda a la creación del
fuselaje, se puede seccionar el fuselaje en cuatro partes, tal y como aparece recogido
en el libro de referencia [6], de este modo cada una de estas cuatro partes corresponde
a un cuerpo de revolución, de los cuales se calculará la superficie exterior (también
se calculará el volumen del fuselaje a mayores),y se sumarán obteniendo la superficie
total del fuselaje. El cálculo comentado aparece recogido en el anexo C del presente
proyecto.
7.5.1.2 Ala
Hasta este momento, lo único que se tenía en relación al ala era el Aspect Ratio
(A), que se trata de una relación entre la envergadura y la superficie alar, si bien
una vez que se realice la primera parte, se tiene el valor aproximado de carga alar,
por lo que ahora mismo es posible definir de una forma más adecuada dicha ala,
se pueden obtener los valores de la envergadura y se pueden empezar a dar otros
valores como la relación de cuerdas o el valor de la cuerda media.
59
7. ANÁLISIS DE SOLUCIONES
Se debe tener en cuenta que el valor de la superficie obtenida incluye el área que
se extiende hasta la línea central del fuselaje.
Los valores que tome el usuario sobre estos parámetros se basarán por tanto en
valores históricos de otras aeronaves (como es el caso del uso del Aspect Ratio),
o a partir de estudios estadísticos y utilizando las consideraciones teóricas que se
aportarán a continuación.
Por otra parte, al igual que se explicó en el caso del fuselaje, en el ala también son
necesarios tener en cuenta ciertos compromisos estructurales en el diseño, puesto que
en la mayoría de aeronaves el combustible se guarda en las alas (lo cual implica zonas
con peso y que requieren de ciertos espesores) y los motores de ciertas aeronaves
también se incorporan en las alas.
Teniendo en cuenta esto último comentado, las alas de la aeronave no solo será
un elemento que soporte la fuerza de sustentación, sino también una zona que debe
soportar de la incorporación de elementos con gran peso y que requieren que estas se
encuentren reforzadas de forma adecuada, integrando componentes como costillas
y larguerillos; si bien, tal y como se ha dicho en el estudio, estos detalles no se
están teniendo en cuenta en el programa de diseño, pero ciertamente a medida que
se avance en el diseño de una aeronave, el usuario se tendrá que enfrentar a estos
compromisos, y por tanto se deja en esta sección estos factores a fin de poder tenerlos
en cuenta en esta primera fase del diseño conceptual.
Las ecuaciones utilizadas para el cálculo de los parámetros comentados del ala
se encuentran recogidos en el anexo D del presente proyecto.
7.5.1.3 Cola
60
7. ANÁLISIS DE SOLUCIONES
Puesto que en Raymer [1] ya se dispone de una tabla con datos empíricos de
coeficientes de volumen de cola adimensionalizados, en función del tipo de aeronave,
la adimensionalización seguida será la misma que la que se expone en el libro de
referencia, que supone por tanto que, para el coeficiente de cola vertical adimensio-
nalizado se tenga que:
LV T SV T
CV T = (2)
bw S w
LH T SH T
CHT = (3)
C w Sw
De este modo una vez definidos estos dos coeficientes, se puede obtener las su-
perficies de la cola horizontal y vertical, a partir de:
CH T C w Sw
SHT = (4)
LHT
CV T bw Sw
SV T = (5)
LV T
Donde faltan por obtener los valores de los brazos LV T y LH T y el valor de los
coeficientes; la obtencion de los mismos aparece recogido en el anexo E del presente
61
7. ANÁLISIS DE SOLUCIONES
proyecto.
Tal y como ocurría en el anterior módulo, las fórmulas que se utilicen en los dis-
tintos apartados explicados, aparecerán recogidas en el anexo del presente proyecto
y serán incorporadas al programa de diseño, con lo que el usuario solo tendrá la
necesidad de aportar los valores de los parámetros pedidos.
Estudio de perfiles
62
7. ANÁLISIS DE SOLUCIONES
Una muestra de ello es la siguiente página web [33], donde se tienen una gran
cantidad de perfiles con sus correspondientes características. Dichos datos son es-
timaciones realizadas a partir del conocido programa Xfoil. También es de reseñar
[34], donde aparece una amplia lista con datos de la geometría de los perfiles, bási-
camente es una tabla con ejes X e Y donde se establecen una serie de puntos que
uniéndolos dan la forma de dicho perfil, lo cual es de gran ayuda si el usuario quiere
realizar el estudio aerodinámico por sí mismo (más adelante se dará información
para poder llevar a cabo esto de forma rápida).
63
7. ANÁLISIS DE SOLUCIONES
resistencia.
64
7. ANÁLISIS DE SOLUCIONES
Más adelante se desarrollará de una forma más extensa la polar que se establecerá
del avión y que será la que se utilizará a la hora de realizar el módulo de diseño
avanzado.
65
7. ANÁLISIS DE SOLUCIONES
dCl tmax
Cl α = = 1, 8π 1 + 0, 8 (6)
dα c
Por otra parte, siguiendo con los parámetros recogidos en la figura 8, α0l co-
rresponde al ángulo donde la sustentación es cero, en el caso de una placa plana se
demuestra que α0l coincide con el 0,0 de la gráfica, y al imponer curvatura al perfil
este término se desplaza hacia la izquierda en el eje x.
Una vez explicados los parámetros que son necesarios para seguir con el estudio
aerodinámico, se deja en el anexo F del presente proyecto una pequeña explicación
con los pasos que debe seguir el usuario a fin de escoger el perfil adecuado.
Por otra parte, el usuario debe tener en cuenta que muchos de las aeronaves hoy
en día utilizan los dispositivos conocidos como dispositivos hiper-sustentadores. Son
dispositivos que aumentan notablemente la sustentación que genera el ala, aumen-
tando la cuerda del perfil y su curvatura; la mayoría de ellos se utilizan básicamente
para las misiones de despegue y aterrizaje, a fin de acortar distancias de pista y
facilitar la maniobra al piloto.
66
7. ANÁLISIS DE SOLUCIONES
67
7. ANÁLISIS DE SOLUCIONES
Se tiene por tanto ahora mismo los siguientes parámetros, a partir del perfil que
ya se ha elegido:
Clα
α0l
Clmax
t/cmax
Antes de comenzar el estudio del ala es necesario calcular los datos aerodiná-
micos del perfil con la configuración sucia, para ello se añadirá a la sustentación
ya calculada la parte correspondiente al deflectar algún tipo de dispositivo hiper-
sustentador, por otra parte, se realizará a continuación una pequeña explicación
sobre los distintos tipos de dispositivos existentes y las variaciones que producen en
la sustentación.
1. Normal (Plain)
El flap simple aumenta su curvatura cuando se deflecta, aunque disminuye leve-
mente la superficie alar.
2. Flap de intradós (Split)
Este flap aumenta la curvatura, pero genera un aumento de resistencia muy gran-
de incluso para deflexiones pequeñas
3. Flap ranurado (Slotted)
Similar al flap normal, pero cuando se deflecta presenta una ranura que comunica
el intradós con el extradós que funciona como un dispositivo de control de la capa
límite. Además de aumentar la curvatura, el aire pasa de intradós a extradós,
aportando energía a la capa límite y retrasando su desprendimiento. Aumenta la
sustentación y la resistencia
4. Flap Fowler
Aumenta la curvatura, controla la capa límite y aumenta la superficie alar. Au-
menta la sustentación y no mucho la resistencia. Normalmente se usa del tipo
birranurado y trirranurado.
68
7. ANÁLISIS DE SOLUCIONES
De este modo, incorporando un flap simple o de intradós, ver figura 10, se puede
apreciar como la curva de sustentación en el caso del plain y split flap se desplaza
69
7. ANÁLISIS DE SOLUCIONES
70
7. ANÁLISIS DE SOLUCIONES
cada uno de los dispositivos, ver tabla 3. El usuario deberá elegir el dispositivo
hiper-sustentador teniendo en cuenta aeronaves similares al tipo de avión diseñado.
En el caso de los nombres mostrados en la tabla 3, los dispositivos de borde de
ataque Krueger y leading-edge flaps se denotarán en el programa como Slat1 y Slat2
respectivamente, ver tabla 15.
Tabla 3: valores del incremento de sustentación para distintos tipos de dispositivos
hipersustentadores, (leading edge devices= flaps de borde de ataque) [10]
Una vez explicado como es el cálculo de la sustentación del perfil para la con-
figuración sucia, la incorporación de los mismos en el programa y las ecuaciones
desarrolladas aparecen recogidas en el anexo G del presente proyecto.
Estudio de alas
71
7. ANÁLISIS DE SOLUCIONES
Dicha relación servirá más adelante para poder relacionar la sustentación del
perfil con el del ala.
72
7. ANÁLISIS DE SOLUCIONES
Parece ser que esta aproximación es adecuada para una gran cantidad de aerona-
ves, si bien la realidad es que para alas con cierta flecha y alargamiento elevado dicha
aproximación deja de ser efectiva, es por ello necesario buscar una aproximación más
adecuada a fin de poder manejar distintas aeronaves a diseñar. Siguiendo esta idea
se ha buscado en distintos libros de diseño y cálculo de aviones la mejor opción
para estimar este coeficiente y finalmente se han unido varios métodos empíricos de
obtención de este coeficiente que se mostrarán a continuación.
Por un lado lo primero que hay que realizar es separar dos métodos de obtención
de este coeficiente de sustentación máximo, según el alargamiento o Aspect Ratio que
tenga el ala; de este modo se dividirán los dos métodos en alas con gran alargamiento
y alas con bajo alargamiento, se establecerán los tipos de alargamiento en función
de las siguientes condiciones obtenidas de los libros de referencia [5] y [1]:
4
A> → Alargamiento elevado
(C1 + 1) cos ΛLE
4
A≥ → Alargamiento bajo
(C1 + 1) cos ΛLE
Donde el parámetro C1 depende del taper ratio λ del ala tal y como aparece
recogido en la gráfica de la figura 14. A fin de poder incorporar este dato en el
programa informático, se utilizará el Excel para obtener una regresión de dicha
curva dejando el parámetro C1 como una función del taper ratio, que el usuario ya
habría incorporado en el módulo de geometría.
C1 = 80, 959λ6 − 253, 71λ5 + 299, 83λ4 − 157, 24λ3 + 28, 336λ2 + 1, 946λ
Por otra parte se puede ver en la figura 15 cual es el ángulo de flecha Leading
Edge ΛLE .
Se supone que el usuario ya tiene un Sketch del ala del avión de diseño a partir
73
7. ANÁLISIS DE SOLUCIONES
Figura 14: Representación del parámetro C1 y C2 en función del taper ratio λ. [5]
Figura 15: Dibujo aclaratorio sobre el ángulo de flecha Λ para un ala en planta. [7]
De este modo se tiene el coeficiente de sustentación máximo del ala (CLmax ) según
tenga gran o bajo alargamiento y el ángulo de entrada en perdida del ala para la
configuración limpia.
Para el caso de calcular dichos valores para la configuración sucia (con flaps
desplegados), se seguirá la misma idea seguida en los perfiles. Se deja en la figura
16 representado como variará la curva de sustentación del ala en función del uso de
dispositivos hipersustentadores y que parámetros intervienen en esta nueva curva.
74
7. ANÁLISIS DE SOLUCIONES
Sf lap
∆(CLmax )F lap = 0,9∆(Clmax )F lap KΛ
S
Se debe tener en cuenta que en este caso el coeficiente (CLmax )Sucia es puramente
el del ala con flaps, por tanto queda por estudiar la sustentación que genera la
aeronave en completo.
75
7. ANÁLISIS DE SOLUCIONES
Dónde:
SCola SCanard
(CL0 )T otal = (CL0 )Ala + (CL0 )Cola + (CL0 )Canard (11)
S S
Dónde:
Además se tendrá que tener en cuenta los parámetros de (CL0 )T otal que aparecen
en la ecuación 11, dichos valores se obtendrán a partir de los datos de perfiles del ala,
de la cola y del canard que el usuario deberá establecer en el módulo TablaAero.xlsx.
76
7. ANÁLISIS DE SOLUCIONES
2πAef f Sexposed
CLα = r F (13)
A2 β 2
tan (Λmax,t )2
Sref
2 + 4 + η2 1 + β2
2πAef f Sexposed Scola ∂
(CLα )T otal = r F+ CLαCola 1−
2 2
tan (Λmax,t )2
Sref S ∂α
2 + 4 + Aη2β 1 + β2
SCanard ∂u
+ CLαCanard 1+ (14)
S ∂α
2πAef f Sexposed
(CLα )Sucia = r F
A2 β 2
tan (Λmax,t )2
Sref
2 + 4 + η2 1 + β2
0
CLα c cf 2
. 1+ Kb Kc 1 − 0 sin δf − 1 + ∆CLαCanard + ∆CLαCola (15)
Cl α c c
Una vez obtenido todos estos parámetros ya se tienen todos los parámetros de
sustentación de la aeronave que eran necesarios, llegados a este punto es necesario
obtener otros parámetros fundamentales para el diseño de la aeronave y a tener en
cuenta en el módulo avanzado; que son los parámetros de la resistencia que genera
el avión. Dicha resistencia se modelará a través de lo que se conoce como polar
del avión, que es la gráfica que representa el coeficiente de resistencia en función del
coeficiente de sustentación y que se nombró al comienzo del módulo de aerodinámica.
Tal y como se comentó existen dos formas de representar la polar, la forma más
simple y más conocida es la siguiente:
CD = CD0 + KCL 2
77
7. ANÁLISIS DE SOLUCIONES
A fin de poder tener en cuenta otros tipos de alas diferentes a las elípticas, se
añade un término a esta resistencia inducida que se denomina factor de Oswalld (e)
y que tiene una gran cantidad de formulaciones empíricas; básicamente reúne ciertas
correcciones a esta resistencia teniendo en cuenta estrechamiento y flecha del ala.
De este modo el parámetro K se puede definir como:
1
K= (16)
πAe
CD = CD 0 − K2 CL + K1 CL 2 (17)
A fin de ejemplificar esto se deja a continuación, ver figura 17, una ilustración
donde se muestran dos polares, una para perfiles sin curvatura y otra para un perfil
con curvatura.
78
7. ANÁLISIS DE SOLUCIONES
Figura 17: Ilustración de la diferencia entre dos polares, una con perfiles con curva-
tura y otra sin curvatura. (Uncambered=sin curvatura/Cambered=con curvatura).
[1]
Si tenemos en cuenta la polar sin curvatura, se puede observar que para un ángulo
de ataque cero, se obtiene la mínima resistencia en el punto de corte de la polar con
el eje de abscisas (CL = 0), de este modo se puede decir que CD0 = CDmin y que
por tanto el resto de la curva depende del coeficiente de sustentación al cuadrado al
tener una forma parabólica y ser simétrica respecto al eje de abscisas.
Se puede establecer que la diferencia entre los parámetros CDmin y CD0 será lo
que decline el uso de uno u otro método (cuanto menor sea la diferencia la polar no
compensada tendrá mayor validez) y es obvio que tiene que ver con la curvatura del
perfil.
79
7. ANÁLISIS DE SOLUCIONES
La polar compensada se puede escribir de esta forma a fin de calcular las cons-
tantes que la componen[7]:
2
CD = CDmin + KCL − CLmin−drag
El usuario debe tener en cuenta que aproximar la resistencia que genera la aero-
nave a una polar, es una aproximación a fin de simplificar los métodos de cálculo;
por tanto se puede destacar que existen métodos más complejos a fin de calcular
esta resistencia. Dichos métodos se basan en el estudio de la resistencia inducida a
partir de teoría potencial donde se tiene en cuenta las interferencias entre los distin-
tos elementos de la aeronave y modelos más exactos en el cálculo de la resistencia
parasitaria; dichos métodos de cálculo se pueden obtener a partir de los libros de
referencia [23] y [2].
80
7. ANÁLISIS DE SOLUCIONES
Es por ello que se hace necesario definir las componentes de coeficiente de re-
sistencia mínima CDmin , la constante de resistencia inducida K y el coeficiente de
sustentación de resistencia mínima CLmin−drag .
Por tanto mediante el uso del método de Component Buildup, se calculan los
coeficientes de resistencia mínima. De este modo, siguiendo el esquema de [1], [2]
y [8], se tiene que para régimen subsónico la resistencia mínima se puede escribir
como:
P
k Cfk IF k F Fk Swetk
(CDmin )T otal = + CDT ren aterrizaje + CDP rotuberancias + CDSoporte alas
S
+CDV entana cabina + CDf lotadores (18)
Por tanto a través de la polar del perfil (obtenida a partir de XFLR5 o mediante
datos obtenidos de páginas de referencia) el usuario podrá establecer el ángulo de
ataque en el que se produce la sustentación de resistencia mínima. Bastaría con
sustituir en la ecuación de la sustentación de la aeronave completa, ecuación (10),
el ángulo de ataque establecido y determinar el CLmin−drag .
Todo esto se ha evaluado para la configuración limpia (sin flaps). A fin de tener
la resistencia que se genera en el despegue y aterrizaje (momento donde se utilizarán
los flaps) se añadirá un término nuevo a la polar ya calculada, que será CDf laps .
81
7. ANÁLISIS DE SOLUCIONES
Donde:
bf
∆CDf laps = 0, 0023 ∆f lap
b
El propio Raymer [1] dice que esta aproximación es bastante aproximada, esto
se puede ver a simple vista. En primer lugar usa el mismo tipo de fórmula para
cualquier tipo de flap ya sea de borde de ataque o de borde de salida (ya se ha
visto que las variaciones que producen uno y otro tipo no son las mismas) y por
otro lado solo tiene en cuenta la resistencia parasitaria, cuando sería lógico pensar
que también debería depender del coeficiente de sustentación, puesto que los flaps
se alojan en las alas y la resistencia inducida es casi el 95 % dependiente de las alas
(o al menos así se ha aproximado).
La estimación depende tanto del tipo de flap como de la cuerda del mismo cf
ver figura 43, la deflexión δf y la envergadura del flap bf , ver figura 55.
Sf lap
∆CDf laps = ∆1 ∆2 (19)
Sref
Donde Sf lap es la superficie del flap (ver figura 55) y Sref la superficie calculada
en el diseño básico.
Finalmente se optará por utilizar esta última fórmula, la forma de obtener las
dos componentes ∆1 y ∆2 queda desarrollada en el anexo L del presente proyecto.
82
7. ANÁLISIS DE SOLUCIONES
Se debe tener en cuenta que al igual que ocurre a lo largo de las etapas del
proceso del programa, los datos calculados en el módulo de aerodinámica ya se
habrán pasado al Excel y al MATLAB de estructuras.
Es por ello que este módulo se puede dividir en tres partes que se explicarán a
continuación. Una primera parte llamada envolvente de vuelo, la segunda relacionada
con la envolvente de vuelo para ráfagas y la tercera denominada estimación de peso
en vacío.
Para poder realizar esta parte será necesario haber realizado el módulo de diseño
y el de aerodinámica, a fin de poder obtener los datos necesarios para cumplimentar
las ecuaciones que aparecerán en este módulo y que se explicarán a continuación.
Por tanto el usuario deberá aportar los factores de carga límite en función del
tipo de aeronave que se vaya a diseñar; de este modo en el Excel TablaEstructuras.xlsx
se aportará una tabla con los valores típicos de factores de carga n según tipo de
aeronave, ver tabla 18.
Una vez explicado lo que se quiere estudiar, a fin de ilustrar el resultado que
se quiere obtener, se deja a continuación, Figura 18, un ejemplo de un gráfico de
envolvente de vuelo.
83
7. ANÁLISIS DE SOLUCIONES
Figura 18: Ejemplo de envolvente de vuelo, con los limites estructurales indicados
mediante colores (parte roja, fallo estructural). [9]
Las cargas experimentadas cuando la aeronave encuentra una ráfaga fuerte pue-
den exceder las cargas de maniobra en algunos casos. Para una aeronave de trans-
porte que vuela cerca de tormentas eléctricas o que se encuentra con turbulencias de
aire a gran altura, no es extraño experimentar factores de carga debidos a ráfagas
que van desde un 1,5 negativo a un positivo de más de 3,5, [1].
Se puede explicar este fenómeno entendiendo que al verse expuesto a una ráfaga,
el avión varía su ángulo de ataque y por tanto también varía la sustentación; ver
figura 20.
2W/S
0,88
CgCL ρ
α
2W/S
5,3+
CgCL ρ
n3/4 = 1 ± α
(20)
498W/S
84
7. ANÁLISIS DE SOLUCIONES
Para crucero:
A menos de h=20.000 ft:
85
7. ANÁLISIS DE SOLUCIONES
V = Vcrucero
U = 50 f t/s
V = Vcrucero
U = 66, 67 − 0, 000833h f t/s
Para picado:
A menos de h=20.000 ft:
V = VD
U = 25 f t/s
V = VD
U = 33, 34 − 0, 000417h f t/s
La gráfica que se obtiene para la envolvente de vuelo para ráfagas se deja repre-
sentada en la figura 19.
Una vez realizado el módulo de diseño básico se tiene una idea preliminar del
tamaño de la aeronave; se tiene por ahora el valor de la superficie alar y del empuje
necesario de la aeronave.
Cabe recordar que para calcular estos valores se han utilizado aproximaciones
basadas en aeronaves similares a la de diseño.
Lo que se trata en este módulo es de entrar en detalle en una de las partes que
se establecieron en la parte de diseño básica; dicha parte es la correspondiente a la
de la estimación del peso en vacío de la aeronave We .
86
7. ANÁLISIS DE SOLUCIONES
Una vez calculadas cada una de las partes que componen el peso en vacío de la
aeronave, ya se estará en disposición de poder calcular el centro de masas de la ae-
ronave, punto de la aeronave que será imprescindible para determinar la estabilidad
de la aeronave. Se recalca que en esta parte solo se van a ver las componentes corres-
pondientes al peso en vacío, por lo que el usuario deberá tener en cuenta que otros
pesos como son los pasajeros de la aeronave, los pilotos, la carga que transportará
o el combustible que debe llevar no se contempla en los cálculos de este módulo;
todos estos pesos mencionados son valores que el usuario deberá integrar para el
cálculo de centro de masas por su propia cuenta, excepto el peso de combustible que
se calculará mas adelante.
De este modo a partir de la ecuación 40, se obtienen para cada caso dos formu-
laciones; seguidas por los autores Nicolai y Raymer.
Nótese que para cada ecuación se establecerán con los subíndices R y N en función
del autor Raymer y Nicolai respectivamente, la elección de las ecuaciones de uno
u otro autor se ha realizado comparando resultados con aeronaves similares, por lo
que en el código del programa se pueden encontrar ecuaciones de un autor o del otro
en función del tipo de aeronave elegida por el usuario.
Una vez planteadas cada una de las ecuaciones, se debe de tener en cuenta que
las estimaciones hechas son correspondientes a aeronaves sin tener en cuenta el uso
de ningún componente de material compuesto.
Puesto que estos materiales son de gran uso en la industria aeronáutica y cada
vez el porcentaje de uso de este material en las aeronaves es mayor, se establecerá
una forma de aproximar el peso de la aeronave teniendo en cuenta el uso de dichos
compuesto. Para ello se puede utilizar la tabla 4, donde se establece un porcentaje
de reducción frente al peso de cada estructura.
87
7. ANÁLISIS DE SOLUCIONES
Tabla 4: Porcentaje de reducción del peso en vacío según componente. [1] y [2]
Una vez llevado a cabo el estudio del diseño geométrico, aerodinámico y estruc-
tural de la aeronave; siguiendo el esquema del programa, ver figura 4, se llega a uno
de las partes del diseño avanzado comentadas en el apartado 7.2.2.
Este apartado corresponde a lo que se conoce como RFP y tal como ya se había
comentado se trata de una forma de limitar los posibles valores de empuje y superficie
alar de la aeronave, generando una gráfica donde se relaciona la carga alar con
respecto al empuje/peso; básicamente es un proceso similar al realizado en el diseño
básico, sin embargo, ahora ya no se obtiene un único valor limitante de carga alar
y de empuje/peso, si no que ambos parámetros están relacionados mediante una
ecuación.
Para llevar a cabo el RFP se utilizan distintas actuaciones de mecánica del vuelo,
véase el caso de despegue, crucero, ascenso o aterrizaje.
Para llevar a cabo este módulo el usuario deberá entrar en la Hoja1 del Excel
TablaDiseñoAvanzado.xlsx a fin de poder ejecutar el módulo de MATLAB MecVue-
lo.m.
Una vez ejecutado este módulo el usuario a partir de la gráfica del RFP podrá
escoger un valor de T /W y W/S que difiere del obtenido en el módulo Parte1.m. Este
nuevo valor de las variables comentadas deberán ser incorporadas en la Hoja2 del
Excel TablaDiseñoAvanzado.xlsx. A fin de completar el esquema de diseño de la figura
4 y poder ejecutar el módulo de MATLAB Parte2.m que se definirá mas adelante en
el apartado 7.5.5.
88
7. ANÁLISIS DE SOLUCIONES
Se llega finalmente a la última etapa del proceso de diseño del programa que
es la correspondiente a Parte2.m. Esta parte tal y como se comentó al comienzo de
este proyecto, se trata del segundo dimensionamiento o dimensionamiento del diseño
avanzado.
Por último otro aspecto que variará en este segundo dimensionamiento será la
parte correspondiente a los consumos de combustible.
89
7. ANÁLISIS DE SOLUCIONES
Tabla 5: Tabla para introducir los datos de la segunda misión de vuelo, Hoja2 del
Excel TablaDiseñoAvanzado.xlsx del programa de diseño
Por otra parte la incorporación de los datos no se hará en una tabla aparte, sino
que se realiza en la propia tabla 5. Los datos necesarios y como incorporarlos se
encuentra recogido en el manual del programa del Anexo Q del presente proyecto,
al igual que la forma de establecer la misión de vuelo.
Para el cálculo de las fracciones de peso son necesarios ciertos datos de cada
etapa de vuelo, dichos datos son los que el usuario ya ha incorporado en la misión
de vuelo, ver tabla 5. Los cálculos que se llevan a cabo para establecer las fracciones
de peso de combustible se explicarán a continuación.
Etapas no variables:
• Del inicio:
1-Encendido
2-Taxi
90
7. ANÁLISIS DE SOLUCIONES
3-Despegue
4-Ascenso con aceleración
• De final:
5-Descenso
6-Espera
7-Aterrizaje
Antes de explicar detalladamente las ecuaciones de mecánica del vuelo con las
que se han obtenido las fracciones de peso en cada etapa, se debe destacar que al
igual que en el caso del módulo de MecVuelo.m, el módulo Parte2.m también se
encuentra relacionado con el archivo funciónCálculoAerodinámicos.m y también tiene
en el código incorporado las tablas de atmósfera estándar.
Por otra parte, a fin de poder utilizar las ecuaciones obtenidas del libro de re-
ferencia [4], se ha tenido que recurrir a un modelo del consumo de combustible
especifico instalado o también denominado (TSFC), se trata de un parámetro simi-
lar al utilizado en el diseño básico en el caso de los motores tipo hélice, veanse las
ecuaciones 26 y 28, donde se utiliza este parámetro denominado en este caso Cbhp .
Las unidades de TSFC están en 1/hora.
91
7. ANÁLISIS DE SOLUCIONES
√
T SF C = (C1 + C2 M ) θ (21)
En el caso de que el usuario esté diseñando una aeronave tipo hélice, el consumo
de combustible especifico instalado sigue la siguiente regresión:
√
T SF C = (0, 18 + 0, 8M ) θ
√
T SF C = (0, 45 + 0, 54M ) θ
Otro parámetro que aparecerá a lo largo del desarrollo de las ecuaciones de las
fracciones de peso es δ que se trata de la relación de la presión con la presión
estándar.
92
8. RESULTADOS FINALES
8. Resultados finales
Por otra parte, los parámetros que ciertamente suponen un resultado final, son
los referentes al dimensionamiento y a la carga alar y relación empuje peso (RFP),
puesto que son estos los que definen en cierta manera el diseño conceptual de la aero-
nave y que se obtienen al realizar todo el proceso de diseño del programa informático
diseñado.
Es por ello que se deja a continuación un breve resumen de los resultado finales
que se obtienen de estas dos partes; del RFP y del dimensionamiento final del diseño
avanzado.
Por un lado, a fin de obtener los resultados del RFP el usuario deberá rellenar la
Hoja1 del Excel TablaDiseñoAvanzado.xlsx con los datos de actuación de la aeronave
que está diseñando, los datos y como debe establecerlos se explica detalladamente
en el manual del programa que se encuentra en el anexo Q del presente proyecto.
Siguiendo las curvas obtenidas en la figura 21, se tiene que la curva azul os-
curo representa la actuación de despegue, la curva roja crucero, la curva amarilla
velocidad de ascenso, la verde giro sostenido, la azul clara el aterrizaje, la de color
93
8. RESULTADOS FINALES
Figura 21: Gráfica obtenida del RFP en el módulo de MATLAB MecVuelo.m del
programa de diseño
negro supone una actuación establecida como ascenso tras despegue, a fin de ver
las limitaciones de la aeronave una vez terminado el despegue, y por último la recta
de color magenta representa la velocidad de entrada en pérdida donde la carga alar
2
de esta componente no es mas que W/S = 0, 5ρVStall CLmax , donde la velocidad de
entrada en pérdida es aquella impuesta al comienzo del diseño básico, ver tabla 1.
El usuario una vez obtenida esta gráfica, ver figura 21, deberá elegir un punto de
la gráfica por encima de las curvas obtenidas, a fin de seleccionar el punto de diseño
de la aeronave, ver figura 22.
Así por ejemplo en el caso de la figura 22, se tiene eligiendo ese punto un valor
de aproximadamente W/S ≈19 y de T /W ≈0,7.
94
8. RESULTADOS FINALES
Es por ello que una vez que se realiza un primer diseño es necesario volver a
repetir todo el proceso, a fin de afinar el diseño.
Los resultados obtenidos de carga alar, relación empuje peso y de peso total de
la aeronave aparecerán en la hoja resultados del Excel TablaDiseñoAvanzado.xlsx,
pero también serán enviados al Excel TablaGeometria.xlsx, mediante PasarDatosMo-
duloGeometria.m, y así el usuario podrá continuar con el proceso de diseño de la
aeronave, ver figura 4.
Una vez realizado el proceso repetidas veces y variando los parámetros que se
quiera, a fin de mejorar el diseño se llegará a un punto donde las variaciones que
se produzcan entre una iteración y otra sean pequeñas o casi inexistentes. En ese
punto el diseño conceptual de la aeronave se puede dar por cerrado.
95
9. PLANIFICACIÓN
9. Planificación
Para la realización del presente proyecto se han establecido una serie de tiempos
de trabajo, fraccionados según la parte correspondiente del mismo.
Índice: 2 horas.
Planos: 4 horas.
Presupuesto: 3 horas.
Por otra parte se deja a continuación un diagrama Gantt con el tiempo de cada
una de las tareas y su organización de realización, ver figura 23.
97
9. PLANIFICACIÓN
1. Presupuesto.
2. Pliego de condiciones
99
Documento III
Planos
101
No aplica
103
Documento IV
Pliego de condiciones
105
Pliego de condiciones
El presente Pliego de Condiciones tiene por objeto definir las fases del proyecto,
fijar las condiciones técnicas y económicas, así como las condiciones generales y
contractuales que han de cumplirse en la ejecución de las tareas previstas.
Condiciones Generales
107
Realizar una primera aproximación del diseño de la herramienta informáti-
ca, llegando a definir los aspectos básicos de dicha herramienta así como un
diagrama o esquema de la estructura interna.
Cumplir con la programación inicial definida y aprobada por las partes para
la ejecución del proyecto.
108
Cumplir con lo estipulado en el artículo 24 “Coordinación de actividades em-
presariales” de la Ley 31/95 de Prevención de Riesgos Laborales, en especial
con el requisito de información al autor del proyecto, sobre los riesgos y medi-
das de protección asociadas a los posibles trabajos a desarrollar dentro de las
instalaciones de la empresa peticionaria.
En este apartado se exponen las actividades y las condiciones bajo las que se
prevé la realización, desarrollo y ejecución del proyecto y la implantación del Sistema
de Gestión de Proveedores en la empresa.
4. Diseño de las partes del diseño básico: Realización de los documentos de Excel
y Matlab correspondientes.
9. Embellecer aquellos ficheros destinados al uso por parte del usuario del pro-
grama.
Entregable 1
109
Entregable 2
Entregable 3
Estudio del funcionamiento global del programa, realización de guía de uso del
programa y mejorar el diseño de uso para el usuario.
Penalizaciones
Si tras el vencimiento del plazo establecido para cada uno de los Entregables,
estos no han sido recepcionados en el plazo de 2 días, se penalizará con 1 500.00
euros cada 7 días de demora.
Primas
Mantenimiento
Cualquier mejora o modificación que tenga que ejecutarse por necesidad o man-
tenimiento correrá a cargo del contratista.
Contratista
La persona que realice el proyecto debe haber adquirido las competencias del
Grado en Ingeniería Aeroespacial y poder ejercer la profesión regulada de Ingeniero
Técnico Aeronáutico.
110
Contrato
Causas administrativas.
Condiciones particulares
Todas las fases del diseño relativas a este proyecto, así como la posterior pro-
ducción, estarán dirigidas por el mismo Ingeniero Técnico Aeronáutico que realizó
el presente proyecto de la herramienta informática de diseño.
111
Documento V
Presupuesto
113
Costes directos
Hardware
Precio Meses de uso Coste
Ordenador 1.200 e 6 150 e
Software
Precio Meses de uso Coste
Catia V5 700 e/año 0,1 6e
Matlab 800 e/año 6 400 e
Microsoft Excel 10 e/mes 6 60 e
Total 616 e
Costes indirectos
En la siguiente tabla 8 se muestran los costes que influyen de una manera indi-
recta en la realización de este proyecto.
115
Tabla 8: Costes indirectos
116
Documento VI
Anexos
119
A. MISIÓN DE VUELO EN EL DISEÑO BÁSICO
Una vez definido esto, se plantea la misión de vuelo; puesto que la herramienta
realizada para el diseño de aeronaves abarca distintos tipos de aeronaves, es necesario
delimitar el tipo de misión en dos partes, por un lado una parte fija, que se supondrá
común para cualquier aeronave, y por otro una parte que variará en función del tipo
de misión que quiera plantear el usuario.
Las etapas al final de cada misión: descenso, espera (por normativa, siempre se
establece en una misión de vuelo antes del aterrizaje una espera) y por último
aterrizaje, taxi y apagado.
Por otro lado, las etapas que podrán variar a lo largo de las distintas misiones
son las etapas que se denominan de vuelo, ver figura 24, y que corresponden al
crucero (dependerá de distancia recorrida), espera (dependerá el tiempo de espera
que realice la aeronave), descenso o ascenso (etapas que se pueden encontrar en
medio de la misiones de crucero de la aeronave) y de carga o descarga (para el caso
de la misión de carga de pago variable, en este caso se modelará como una etapa
donde no hay perdida de combustible, pero sin embargo, la masa de la aeronave
pierde o gana la cantidad de carga de pago variable).
El usuario, una vez definida la misión de vuelo con las correspondientes etapas,
deberá establecer esta información en el Excel TabladeTodoslosdatos.xlsx , la forma
121
A. MISIÓN DE VUELO EN EL DISEÑO BÁSICO
Tabla 10: Tabla para introducir los datos de la misión de vuelo, Hoja2 del Excel
TabladeTodoslosdatos.xlsx del programa de diseño
Una vez planteada la misión de vuelos, con las etapas definidas, es necesario
plantear las ecuaciones y datos estadísticos utilizados en cada etapa, a fin de que
el usuario tenga conocimiento de los datos que necesitará aportar al programa para
obtener el cálculo del peso total y de combustible.
We
W0 = Wcrew + Wf p + Wvp + Wf + W0 (22)
W0
Por otra parte, Wf se obtendrá a partir del sumatorio de los pesos de combustible
consumidos en cada etapa. Para ello:
X
Wf = 1, 06 Wf i W0 (23)
Donde el sumatorio de pesos queda multiplicado por un factor de 1,06 que supone
imponer un 6 % más de combustible a fin de tener en cuenta posibles emergencias
en vuelo.
122
A. MISIÓN DE VUELO EN EL DISEÑO BÁSICO
Wi
Wf i = 1 − Wi (24)
Wi−1
De este modo, a partir de los valores de las fracciones de peso de cada etapa WWi−1
i
se podrá obtener el valor del peso de combustible consumido en cada etapa, y con
ello el peso total de combustible.
Tabla 11: Datos para el cálculo del dimensionamiento de la aeronave, Hoja1 del
Excel TabladeTodoslosdatos.xlsx del programa de diseño
Fracciones de peso,según etapa de misión Valores de eficiencia Fracción de peso en vacío
INICIAL FINAL Crucero Espera W e/W o
Indicador Encendido Taxi Despegue Ascenso Descenso Aterrizaje CL/CD (CR) C Rendimiento Propulsivo CL/CD(ES) C Rendimiento Propulsivo A c
Homebuilts
1 0,998 0,998 0,998 0,995 0,995 0,995 9 0,7 0,7 11 0,6 0,65 0,86 -0,05
(Metal/wood)
Homebuilts
2 0,998 0,998 0,998 0,995 0,995 0,995 9 0,7 0,7 11 0,6 0,65 0,91 -0,05
(Composite)
Single
3 0,995 0,997 0,998 0,993 0,992 0,993 9 0,6 0,8 11 0,6 0,7 2,36 -0,18
engine
Twin
4 0,992 0,996 0,996 0,992 0,99 0,992 9 0,6 0,82 10 0,6 0,72 1,51 -0,1
engine
Agricultural 5 0,996 0,995 0,996 0,999 0,998 0,998 7 0,6 0,82 9,5 0,6 0,72 0,74 -0,03
Businessjet 6 0,99 0,995 0,995 0,99 0,98 0,992 12 0,7 - 13 0,5 - 1,51 -0,1
Regional
7 0,99 0,995 0,995 0,985 0,985 0,995 12,5 0,6 0,85 15 0,5 0,77 0,96 -0,05
TBP’s
Transport
8 0,99 0,99 0,995 0,99 0,98 0,992 15 0,7 - 16 0,5 - 1,02 -0,06
jets
Mil.
9 0,99 0,99 0,99 0,99 0,98 0,995 9 0,6 0,82 11 0,6 0,77 1,59 -0,1
Trainers
Transporte
10 0,99 0,99 0,995 0,99 0,98 0,992 14 0,6 0,82 15 0,6 0,77 0,93 -0,07
Militar (Hélice)
Transporte
11 0,99 0,99 0,995 0,99 0,98 0,992 14 0,7 - 15 0,5 - 0,93 -0,07
Militar (Jet)
Anfibio
12 0,992 0,99 0,99 0,99 0,985 0,99 11,5 0,6 0,82 14,5 0,6 0,77 1,09 -0,05
(Hélice)
Anfibio
13 0,992 0,99 0,99 0,99 0,985 0,992 11,5 0,7 - 14,5 0,5 - 1,09 -0,05
(Jet)
Ecuaciones crucero
Wi −RC
= exp (25)
Wi−1 V (L/D)
123
A. MISIÓN DE VUELO EN EL DISEÑO BÁSICO
Wi −RCbhp
= exp (26)
Wi−1 550ηp (L/D)
Ecuaciones de espera
Wi −RE
= exp (27)
Wi−1 L/D
Wi −EV Cbhp
= exp (28)
Wi−1 550ηp (L/D)
Tal y como se puede observar, las fracciones de peso tienen ecuaciones distintas
en estos dos casos, en función del tipo de propulsor de la aeronave (ténganse en
cuenta que al principio del apartado se hace hincapié en la elección del motor de la
aeronave a diseñar).
En este momento el usuario deberá, por tanto, a fin de completar los valores
de las fracciones de peso, aportar el valor de la duración de la espera, la velocidad
de crucero de la aeronave, además de valores estadísticos de rendimiento propulsivo
y de consumos específicos de los motores, estos dos valores recogidos en la propia
tabla del programa a fin de facilitar al propio usuario, a partir de la referencia [2].
Téngase en cuenta que dichos valores se pueden modificar si se considera oportuno.
Tal y como se puede ver en estas ecuaciones 25, 26, 27 y 28, aparece un término
conocido como eficiencia aerodinámica, L/D, el usuario se enfrenta en este momento
al primer dato aerodinámico de la aeronave que quiere diseñar. Al igual que ocurría
con otros de los datos, este parámetro también aparece incorporado en la propia
tabla, a partir del libro de referencia [1].
124
A. MISIÓN DE VUELO EN EL DISEÑO BÁSICO
aeronave, se puede llegar a obtener un valor de L/Dmax aproximado, por otro lado
si se busca obtener la eficiencia aerodinámica de una forma más precisa se puede
recurrir a un programa informático, denominado OpenVSP (de código libre), donde
se realiza un boceto de la aeronave que se busca diseñar, y con él se obtienen los
datos geométricos de la aeronave, para finalmente solamente tener que utilizar la
gráfica de referencia [1] que aparece en la figura 25.
Crucero Espera
Jet 0,866L/Dmax L/Dmax
Hélice L/Dmax 0,866L/Dmax
We
= AW0 c Kvs (29)
W0
125
A. MISIÓN DE VUELO EN EL DISEÑO BÁSICO
Establecidas las etapas de la misión de vuelo y todos los datos comentados an-
teriormente, bastará con ejecutar el programa de MATLAB correspondiente a la
Parte1.m, para obtener los valores finales de peso máximo despegue W0 , peso final
de combustible Wf t y el peso de combustible gastado en cada etapa Wf i . El proceso
de cálculo se basa en un proceso de iteración, donde al aportar un valor al W0 de
la ecuación 22 se genera un resultado con los términos a la derecha del igual, que
supone continuar con el cálculo hasta que se obtienen valores iguales o muy próximos
en ambos lados de la ecuación.
We
W0 (i) = Wcrew + Wf p + Wvp + Wf t (i − 1) + (i − 1) W0 (i − 1) (30)
W0
We
(i − 1) = A[W0 (i − 1)]c Kvs
W0
N
X +1
Wf t (i − 1) = 1, 06 (1 − C(k + 2)) W (k + 1)
k=1
W (k + 1) = W (k) + Wvp
126
B. RELACIÓN EMPUJE PESO Y CARGA ALAR
Estático
Máxima velocidad
T
= aMmax C (31)
W
P
= aVmax C (32)
W
Crucero/Despegue
T T WCrucero TDespegue
= (33)
W Despegue W Crucero WDespegue TCrucero
127
B. RELACIÓN EMPUJE PESO Y CARGA ALAR
P VCrucero T WCrucero PDespegue
= (34)
W Despegue ηp 550 W Crucero WDespegue PCrucero
Ascenso
T p
= G + 2 CD0 /Aeπ (35)
W
Donde G = VV ertical /VCrucero , CD0 resistencia parasitaria, e es el factor de Oswall y
A el Aspect Ratio.
P p VCrucero
= G + 2 CD0 /Aeπ (36)
W 550ηp
Estático
Entrada en pérdida
128
B. RELACIÓN EMPUJE PESO Y CARGA ALAR
Para este caso se utilizará la siguiente ecuación para determinar la carga alar,
que variará en función de los siguientes parámetros:
W
= 0, 5ρVStall 2 CLmax (37)
S
Despegue
Para este caso se utilizará la siguiente ecuación para determinar la carga alar,
que variará en función de si la aeronave es de tipo jet o tipo hélice, de este modo se
tiene que:
Para jet
W T
= (T OP )σCLT O (38)
S W
W P
= (T OP )σCLT O (39)
S W
129
B. RELACIÓN EMPUJE PESO Y CARGA ALAR
Aterrizaje
W SAterrizaje − Sa
= σCLmax (40)
S 80
Crucero
Para jet
p
W/S = q πAeCD0 /3 (41)
p
W/S = q πAeCD0 (42)
Ascenso
D T
= −G (43)
W W
qSCC 0 + qS CL 2 /πAe
D qCC 0 W
= = + (44)
W W W/S qSπAe
130
B. RELACIÓN EMPUJE PESO Y CARGA ALAR
q 2
T T
W W
−G± W
− G − (4CD0 /πAe)
= (45)
S 2/qπAe
Donde todos los parámetros que aparecen ya han sido definidos anteriormente.
131
C. CÁLCULO DE LA GEOMETRÍA DEL FUSELAJE
Se parte para definir el fuselaje de 4 volúmenes de revolución, ver figura 26, que
son un paraboloide de base circular, un cilindro y troncos de cono que forman parte
de la base de cola.
De esta forma, se puede calcular la superficie del fuselaje a partir de las superficies
que generan dichos cuerpos de revolución y los parámetros que los definen, ver figura
27; se debe tener en cuenta que el usuario podrá variar los parámetros de cada
volumen, por tanto si quiere que la base de cola acabe en punta, bastaría con hacer
cero el diámetro d2 , o si se quiere modelar el fuselaje teniendo en cuenta solamente
un tronco de cono se puede hacer cero las longitudes L4 o L3 .
Figura 27: Sección del fuselaje modelado con Catia, estableciendo parámetros que
lo definen.
Teniendo en cuenta la figura 27, se tiene por tanto los siguientes parámetros:
133
C. CÁLCULO DE LA GEOMETRÍA DEL FUSELAJE
Longitud L1
Diámetro D
Cilindro (2)
Longitud L2
Longitud L3
Longitud L3
LF uselaje = L1 + L2 + L3 + L4 (46)
s
2
D 2 − d1 2
πD 2 D π(D + d1 ) 2
SF uselaje = 4L 1 + + πDL 2 + L 3 +
12L 2 4 2 4
| 1
| {z }
{z } Cilindro | {z }
P arboloide T roncodecono
s
π(d2 + d1 ) d2 − d2 2
+ L4 2 + 1
2 4
| {z }
T roncodecono
134
C. CÁLCULO DE LA GEOMETRÍA DEL FUSELAJE
Con lo que reagrupando términos, se tiene que la superficie del fuselaje es:
" # s
2 1,5 3
D 2 − d1 2
πD 1 2 D D 2(D + d1 ) 2
SF uselaje = 4L1 + − − D + 4L2 + L3 +
4 3L1 2 4 8 D 4
s
2(d2 + d1 ) d2 − d2 2
+ L4 2 + 1 (47)
D 4
De la misma forma se puede calcular también el volumen del fuselaje, siendo este
a partir del libro de referencia [6]:
π 2 L1 L3 L3 2
L4 2 2
Vf uselaje = D ( + L2 + ) + Dd1 + d1 + d2 d1 + d1 + d2
4 2 3 3 3
(48)
135
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
Siguiendo los parámetros que se tienen de la figura 29, se define ctip es la cuerda
en la punta, croot es la cuerda en el encastre o línea media del fuselaje y b/2 es la
semi envergadura del ala.
137
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
a partir de los parámetros de figura comentada, sin más que aplicando trigonometría
básica para el caso de los dos trapecios que se forman (correspondientes al ala entera),
con lo que se tiene que:
b
S = (croot + ctip ) (49)
2
Por otra parte, es necesario definir una cuerda media adimensionalizada, o cuerda
media geométrica (M GC), que supone definir una cuerda, que en muchas ocasiones
coincide con la cuerda media aerodinámica (M AC) que es la cuerda donde se loca-
lizan las fuerzas de sustentación, resistencia y momentos que se generan desde una
perspectiva 2D.
Por otra parte, puesto que el caso de estudio se encuentra en régimen subsónico
(se sabe que por teoría el centro aerodinámica de un perfil en subsónico se encuentra
a 1/4 de la cuerda), el cuarto del M AC será el centro aerodinámico del ala, y
por ello su posición es de importancia en cuanto a consideraciones de estabilidad
longitudinal.
M GC = M AC (50)
Nótese que este último supuesto no siempre es cierto, el que no sean iguales
dependerá en mayor medida de la flecha, torsión del ala y del propio ángulo de
ataque.
138
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
Z b/2
2
M GC = C = C 2 dy (51)
S 0
2 croot ctip
C= croot + ctip + (52)
3 croot + ctip
λ2 + λ + 1
2
C = croot (53)
3 λ+1
b 1 + 2λ
Y = (54)
6 1+λ
En el caso de una cola vertical el valor de Y será el doble del calculado anterior-
139
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
2(λ − 1)
c(y) = croot 1 + y (56)
b
Se debe tener en cuenta que la longitud y posición de esta cuerda media aerodi-
námica solo es válida para un ala trapezoidal. Esta construcción no se puede aplicar
en aviones que presenten una discontinuidad en la cuerda a lo largo del ala, o bien
si se utiliza, debe ser considerada como una aproximación.
La forma más cómoda para obtener todo estos valores, sería en todo caso, y tal
como se ha explicado al comienzo de este apartado, mediante un programa de diseño
gráfico, como alguno de los mencionados.
D.2. Estrechamiento
2S 2b
croot = = (58)
b(1 + λ) (1 + λ)A
140
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
D.3. Alargamiento
b2 b
A= = (59)
S Cmedia
croot + ctip
Cmedia = (60)
2
Este parámetro A cuanto más grande es, supone una mayor eficiencia aerodiná-
mica, mayor momento flector, mayor peso del ala, mayor efecto suelo y una menor
maniobrabilidad en balance [11].
141
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
Es fácil ver que teniendo en este caso la superficie alar S (a partir del dimensiona-
miento) y establecido
√ un valor del Aspect Ratio es posible establecer la envergadura
del ala b = AS, que básicamente es la medida del ala de punta a punta.
El parámetro de envergadura también puede que tenga que ser limitado por
parte del usuario. Por ejemplo en el caso de ciertas limitaciones de aeropuertos
donde estará la aeronave de diseño, al exceder los límites de ancho de pista, o para
entrar la propia aeronave en ciertos hangares.
λ2 + λ + 1
4b
C= (61)
3A λ+1
Por otra parte, a mayor espesor, se tiene un menor Mach de divergencia de resis-
tencia, un menor peso del ala (aunque suene extraño, para espesores más pequeños,
se requieren mayores componentes estructurales, y por tanto supone un mayor peso
final del ala) y una mayor capacidad de combustible. [11]
En primer lugar, para entender lo que supone incorporar flecha y espesor al ala
de una aeronave, es necesario comprender los parámetros que se han tratado. En este
caso se ha definido un parámetro denominado Mach de divergencia, el cual puede
parecer que no afecte al diseño de una aeronave en régimen subsónico, pero esto no
es así; para aeronave en régimen subsónico alto (Mı́nf>0,6 ) se tienen zonas del ala
(se debe tener en cuenta que el flujo en las zonas del extradós es acelerado) donde
se alcanzan valores del Mach cercanos o que llegan al valor sónico (M = 1). En el
momento en el que esto ocurre, el Mach al que vuela la aeronave se denomina Mach
crítico (Mcr ), donde se tiene siempre que el Mach crítico es menor que 1, puesto que
si no se estaría hablando de régimen supersónico.
142
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
De este modo es fácil entender que para una aeronave que vuelen en régimen sub-
sónico alto, será necesario diseñar un ala con cierta flecha, a fin de poder aumentar
el valor del Mach de divergencia, a la que la resistencia aumenta tan bruscamente;
es por ello que en la mayoría de aeronaves convencionales de transporte de pasajeros
presentan cierta flecha, a fin de evitar este fenómeno en zonas en las que se vuele a
Mach subsónico elevados.
Por otra parte, el espesor del ala se elegirá en función de las necesidades que
tenga el usuario a la hora de diseñar la aeronave (guardar combustible, cálculo de
pesos), si bien, lo que limitará en gran medida la elección del mismo es en la parte
aerodinámica (perfil del ala), que se discutirá más adelante; se puede resumir su uso
como:
Por un lado se debe definir el ángulo de flecha del borde de ataque ΛLE , y a
partir de este se puede definir cualquier ángulo respecto a la cuerda, ver figura 34,
a partir de la siguiente fórmula:
4 1−λ
tan Λn = tan Λm − (n − m) (62)
A 1+λ
143
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
4 1 1−λ
tan Λc/4 = tan ΛLE − ( − 0)
A 4 1+λ
2S
Que teniendo en cuenta la definición de cuerda en el encastre croot = b(1+λ)
=
2b
(1+λ)A
se tiene que:
croot 1 − λ
Λc/4 = arctan tan ΛLE −
2b 1+λ
Figura 32: Representación del ala trapezoidal de referencia, con todos los parámetros
que la definen. [6]
D.5. Torsión
Esta torsión, giro relativo de las cuerdas, se suele hacer gradualmente desde el
encastre a la punta del ala. Si bien este parámetro no se tendrá en cuenta en el diseño
conceptual, se debe saber que es un parámetro utilizado en la aviación general para
variar la forma de distribución de sustentación de un ala asemejándola a la de un
ala elíptica, y se utiliza principalmente en alas con flecha para alejar el comienzo de
la entrada en pérdida de la punta del ala.
D.6. Diedro
Se puede definir otro parámetro del ala que es el diedro. Tiene influencia en
la estabilidad lateral del avión y por tanto, no se requiere su incorporación en la
144
D. CÁLCULO DE LA GEOMETRÍA DEL ALA
Figura 33: Representación del ángulo de disedro de una aeronave, ilustración obte-
nida de [15]
El ángulo de incidencia se escoge de modo que la resistencia del avión sea mínima
en alguna condición de vuelo, normalmente en crucero. Su valor suele estar entre 0o
y 2o . Este parámetro es requerido para el estudio de estabilidad de la aeronave, por
lo que, de momento, en esta fase del diseño no es necesario [15].
145
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
CH T C w Sw
SHT =
LHT
CV T bw Sw
SV T =
LV T
147
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
bien a fin de calcular la superficie de este elemento se debe seguir la misma idea que
se ha explicado con los anteriores casos. El usuario tendrá al igual que en todos los
restantes casos los parámetros en el Excel del módulo de geometría, donde podrá
incorporar este valor, y se dejará a continuación una definición de los parámetros
que históricamente se han utilizado.
Para los coeficientes se tiene la siguiente tabla 13. El usuario debe tener en cuenta
que dicha tabla ya aparece incorporada en el módulo de geometría, y su elección
está definida en función del tipo de aeronave que se haya elegido.
Tabla 13: Tabla con los valores típicos de los coeficientes del cálculo de cola horizontal
y vertical [1].
Queda por determinar el método para obtener en un primer caso los brazos o
distancias del estabilizado horizontal y vertical con respecto al ala de la aeronave.
LT = LV T = LH T
Con lo que las ecuaciones 4 y 5 quedan como:
bw Sw CV T
SV T =
LT
C w S w CH T
SH T =
LT
148
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
El método que se explicará se basa en tres pasos y sigue la idea del libro de
referencia [6]:
Paso 1
El usuario debe tener en cuenta que la forma como halla estructurado el fuselaje,
será la misma que utilice ahora y que por tanto es en este momento cuando podrá
establecer la posición del ala en el fuselaje. La siguiente figura 36 se muestra como
ejemplo de diseño, como se puede observar en este caso se tiene que el ala se encuentra
en el volumen del cilindro, por tanto el largo del brazo LT será la suma del largo
del tronco (3) más la suma del tronco (4) menos 3/4 de la cuerda del estabilizador
horizontal y más el largo que corresponde de 1/4 de la cuerda del ala hasta la unión
tronco de cono-cilindro, que será la distancia que podrá variar en función del largo
del brazo.
Figura 36: Ilustración del brazo de cola LT en función de los volúmenes de revolución
definidos, modelado en Catia.
149
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
denominaran como S las superficies mojadas que se van a calcular a fin de simplificar
terminología.
s
π(D + d1 ) D 2 − d1 2
SF uselajeCola = πDX
| {z } + L3 2 +
2 4
Cilindro | {z }
T roncodecono
s
π(d2 + d1 ) d2 2 − d1 2 3
+ L4 2 + − crootHT
2 4 4
| {z }
T roncodecono
Donde X es la distancia del cuarto de cuerda del ala hasta el tronco de cono (3),
se tiene por tanto que:
s s
2
D2 − d1 d2 2 − d1 2
LT = X + L3 2 + + L4 2 + (63)
4 4
s s s
2 2 2
D2 − d1 d2 − d1 π(D + d1 ) D 2 − d1 2
SF uselajeCola = πD LT − L3 2 − + L4 2 + + L3 2 +
4 4 2 4
| {z }
| {z } T roncodecono
Cilindro
s
π(d2 + d1 ) d2 2 − d1 2 3
+ L4 2 + − crootHT (64)
2 4 4
| {z }
T roncodecono
bw Sw CV T C w Sw CH T
SCola = 2SV T + 2SHT = 2 + (65)
LT LT
150
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
s s s
2 2 2
D2 − d1 d2 − d1 π(D + d1 ) D 2 − d1 2
SLT = πD LT − L3 2 − + L4 2 + + L3 2 +
4 4 2 4
| {z }
| {z } T roncodecono
Cilindro
s
d2 2 − d1 2 3
π(d2 + d1 ) 2 bw Sw CV T C w Sw CH T
+ L4 + − crootHT +2 + (66)
2 4 4 LT LT
| {z }
T roncodecono
dSLT
=0
dLT
s
2Sw (CH T C w + bw CV T )
LT min = (67)
πD
De este modo ya se tiene el valor del brazo que hace mínima la superficie mojada.
Para calcular cuánto valdrá la distancia del cuarto de cuerda del ala hasta el
tronco de cono (3), basta con sustituir en la ecuación 63 con lo que se obtiene que:
s s s
2 2 2
2Sw (CH T C w + bw CV T ) D2 − d1 d2 − d1
X= − L3 2 + + L4 2 + (68)
πD 4 4
Paso 2
bw Sw CV T
SV T = q
2Sw (CH T C w +bw CV T )
πD
C w S w CH T
SH T = q
2Sw (CH T C w +bw CV T )
πD
Paso 3
151
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
v
C w Sw CH T
p u
bHT = AHT SH T = tAHT q (69)
u
2Sw (CH T C w +bw CV T )
πD
r
C w Sw CH T
AHT q 2Sw (C C w +bw C
HT V T)
πD
CmediaHT = (70)
AHT
Donde:
crootHT + ctipHT
Cmedia =
2
r
C w Sw CH T
AHT q 2Sw (C C w +bw C V T)
2bHT HT
πD
crootHT = =2
(1 + λHT )AHT (1 + λHT )AHT
y la cuerda en la punta:
r
C w Sw CH T
AHT q 2Sw (C C w +bw C
HT V T)
πD
ctipHT = λHT crootHT = 2λHT
(1 + λHT )AHT
v
u bw Sw CV T
bV T = tAV T q (71)
u
2Sw (CH T C w +bw CV T )
πD
r
AV T q 2Sw (Cbw SwCCwV+bT w C
HT V T)
πD
CmediaV T = (72)
AV T
152
E. CÁLCULO DE LA GEOMETRÍA DE LA COLA
153
F. SELECCIÓN DEL PERFIL DEL ALA
L=W
L = 1/2ρVc2 SCL
2W/S
CLCrucero =
ρVc2
155
F. SELECCIÓN DEL PERFIL DEL ALA
CLCrucero
CLCruceroAla =
0, 95
CLCruceroAla CLCrucero
Cl i = =
0, 9 0, 855
156
F. SELECCIÓN DEL PERFIL DEL ALA
mismo, pero en este caso para el caso del coeficiente de sustentación máximo, este
coeficiente tal y como se vio anteriormente se obtiene en el momento que la aeronave
entra en pérdida.
Se dejará por tanto en la hoja del Excel del módulo TablaAero.xlsx el valor que
se obtiene de coeficiente de sustentación máxima en configuración limpia y del co-
eficiente de sustentación ideal. De este modo, a partir de la gráfica expuesta, ver
figura 37, el usuario podrá hacerse una idea de que tipo de perfil elegir. Si bien esto
último es una ayuda a fin de elegir el perfil; el usuario podrá seguir estos pasos o
elegir por su cuenta el perfil que prefiera.
157
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
En el caso que solamente haya flaps en el borde de salida bastará con imponer
a la sustentación de las flaps de borde de ataque como cero, la ecuación por tanto
quedaría como:
c0 ∆c cf
=1+ (74)
c cf c
Donde aparecen los términos ratio cuerda flap y cuerda perfil cf /c . Dicho ratio
relaciona la cuerda del flap sin deflectar con la cuerda del perfil también sin deflectar,
normalmente este valor suele estar entre 0,15 y 0,25 según [1], (se tomará como valor
el de 0,25 al ser el más común en la mayoría de aeronaves).
159
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
De este modo a partir de todo lo expuesto se puede establecer una tabla con
valores medios de dichas variables, a fin de que el usuario tenga una referencia de
160
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
que parámetros utilizar para cada tipo de dispositivo. Se tiene por tanto las siguientes
relaciones de cuerda en la tabla 15 que aparece incluida en el Excel del módulo de
aerodinámica.
Tabla 15: Valores del incremento de sustentación para distintos tipos de dispositivos
hiper-sustentadores, (BA= flaps de borde de ataque, BS= flaps de borde de salida)
[1]
El usuario debe tener por tanto para el cálculo del nuevo ángulo dos parámetros,
por un lado el ángulo de deflexión del flap δf que aparece recogido en la tabla 14 y
el ratio de cuerdas flap/perfil cf /c, que se ha definido con un valor de 0,25.
dCl 1
(∆α0l )P lain = − δf Kf0 (75)
dδf Clα
161
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
Figura 39: Variación del parámetroKf0 en función del ángulo de deflexión δf [5].
sin más que trazando una línea perpendicular al eje x en la coordenada de 60o
y hallando el corte con la curva correspondiente al valor cf /c=0,25 (véase figura
41).
Por otro lado para obtener el valor de dCl /dδf se utiliza la figura 40; se tiene el
valor de cf /c=0,25 con lo que se pueden obtener los datos de dCl /dδf en función
del espesor relativo del perfil t/c (se debe tener en cuenta que el valor de espesor
relativo se toma para el valor de espesor máximo, es decir, t/c = tmax /c, a partir
de una ecuación regresiva realizada en Excel, por lo que se tiene finalmente que:
dCl
= 47, 538(t/c)3 − 13, 834(t/c)2 + 3, 4006(t/c) + 3, 8208
dδf
Reagrupando todo se obtiene finalmente a partir de 75 que para flap tipo Plain:
162
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
Figura 41: Obtención del valor Kf0 para el caso de estudio [5].
dα0l
∆α0l = − δf (76)
dδf
Donde δf es el ángulo de deflexión que para el caso de flap tipo Fowler es de 30o ,
y para el caso de los tipos Single Double o triple Slotted se impondrá un valor de
40o de deflexión.
dα
Por otra parte, el valor de dδ0fl , se puede obtener a través de la gráfica expuesta
en la figura 42, teniendo en cuenta que el valor de cf /c es de 0,25.
dα0l
Figura 42: Variación del parámetro dδf
en función del ángulo de deflexión del flap,
y de la relación de cuerdas cf /c, [5]
163
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
donde el valor en color azul corresponde al flap tipo Fowler, y el color verde
corresponde a los flap tipo Slotted.
dα0l
Figura 43: Obtención del valor dδf
para el caso de estudio [5].
2π30 2π40
(∆α0l )F owler = −0, 45 = −0,24 (∆α0l )Slotted = −0, 45 = −0,24
360 360
1
(∆α0l )Split = − ∆lf Split
k (77)
Cl α
Donde la variación de la sustentación generada por el flap tipo split ∆Clf Split
y el parámetro k se pueden obtener de las gráficas de la figura 44.
El parámetro k por
tanto tomará un valor de 1,1 mientras que para obtener el
parámetro ∆Clf Split se volverá a realizar una tendencia de datos en Excel para
obtener la ecuación regresiva.
Con lo que finalmente queda como:
164
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
165
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
Con lo que se obtiene todos los datos necesarios para representar la curva de
sustentación del perfil con dispositivos hiper-sustentadores. Se debe tener en cuenta
que en esta aproximación se ha establecido la pendiente de la curva de sustenta-
ción del perfil como constante, hecho que no es del todo real, puesto que el uso de
166
G. CONTRIBUCIÓN DE LOS DISPOSITIVOS HIPERSUSTENTADORES A
LA SUSTENTACIÓN DEL PERFIL
Figura 47: Variación del ángulo de entrada en pérdida para los flaps tipo Fowler,
double y triple slotted (αstall )c0 /c en función del ángulo de deflexión, [5]
(Clmax )Limpia + ∆ (Clmax )F lapb.s + ∆ (Clmax )F lapb.a
(Clα )c0 /cSucia = (80)
(αstall )Limpia + (∆αstall )c0 /c − α0lLimpia + ∆α0l
167
H. CÁLCULO DEL COEFICIENTE DE SUSTENTACIÓN MAXIMO Y
ÁNGULO DE ENTRADA EN PÉRDIDA DEL ALA
4
A> → Alargamiento elevado
(C1 + 1) cos ΛLE
4
A≥ → Alargamiento bajo
(C1 + 1) cos ΛLE
De este modo se pasará a explicar en cada caso cada uno de los métodos y su
fundamento teórico.
CLmax
CLmax = Clmax + ∆CLmax
Clmax
169
H. CÁLCULO DEL COEFICIENTE DE SUSTENTACIÓN MAXIMO Y
ÁNGULO DE ENTRADA EN PÉRDIDA DEL ALA
Se tienen que hallar por tanto los términos comentados; en primer lugar Clmax
es un valor conocido, al ya haber realizado el estudio del coeficiente de sustentación
máximo del perfil con configuración limpia. Por otra parte para obtener el valor de
CLmax
Cl
se puede recurrir a la gráfica de la figura 49.
max
C
Figura 49: Variación del cociente CLl max en función de la flecha del ala medida res-
max
pecto al eje del borde de ataque ΛLE y al parámetro ∆y, (Mach=0,2) [5].
Por otro lado para hallar el valor de ∆CLmax se ha de utilizar varias gráficas en
función de la flecha ΛLE . Dichas gráficas se pueden observar en la figura 50.
170
H. CÁLCULO DEL COEFICIENTE DE SUSTENTACIÓN MAXIMO Y
ÁNGULO DE ENTRADA EN PÉRDIDA DEL ALA
Figura 50: Gráficas para distintos valores de ΛLE donde se observa la variación
∆CLmax en función del parámetro ∆y y el número de Mach [1].
CLmax
(αstall )Ala = α0L + + ∆αCLmax (81)
CLα
Donde el primer y segundo término son los términos que se tendrían en el cálculo
del ángulo de entrada en pérdida siguiendo la misma idea que la llevada a cabo en
perfiles (curva de sustentación lineal).
Por tanto para el cálculo de (αstall )Ala se tienen todos los parámetros, excepto
∆αCLmax , que se obtiene a partir de la gráfica de la figura 51.
171
H. CÁLCULO DEL COEFICIENTE DE SUSTENTACIÓN MAXIMO Y
ÁNGULO DE ENTRADA EN PÉRDIDA DEL ALA
Figura 51: Variación del parámetro ∆αCLmax en función de la flecha medida desde
el borde de ataque ΛLE y al parámetro ∆y [5].
172
H. CÁLCULO DEL COEFICIENTE DE SUSTENTACIÓN MAXIMO Y
ÁNGULO DE ENTRADA EN PÉRDIDA DEL ALA
Figura 52: Variación del parámetro ∆CLmax en función del parámetro (C2 +
1)A tan ΛLE y el número de Mach, [5]
Para el cálculo del ángulo de entrada en perdida para alas con bajo alargamien-
to se utiliza la siguiente fórmula, que sigue una forma similar a la expuesta del
coeficiente de sustentación máximo:
Figura 53: Variación del parámetro de ángulo de ataque para coeficiente de susten-
tación máximo “base” (αCLmax )base en función del parámetro (C1 + 1)A/β cos ΛLE
[5].
173
H. CÁLCULO DEL COEFICIENTE DE SUSTENTACIÓN MAXIMO Y
ÁNGULO DE ENTRADA EN PÉRDIDA DEL ALA
p
β= 1 − M∞ 2 → P ara regimen subsnico
Para el parámetro (αCLmax )base se puede, al igual que se hizo con los parámetros
C1 y C2 , obtener una ecuación recursiva en función del parámetro (C1 +1)A/ cos ΛLE
que queda de la siguiente forma:
(αCLmax )base = 0, 444[(C1 + 1)A/ cos ΛLE ]6 − 5, 0505[(C1 + 1)A/ cos ΛLE ]5
+21, 251[(C1 + 1)A/ cos ΛLE ]4 − 38, 993[(C1 + 1)A/ cos ΛLE ]3
+26, 373[(C1 + 1)A/ cos ΛLE ]2 − 5, 0622[(C1 + 1)A/ cos ΛLE ] + 35
Figura 54: Variación del parámetro ∆αCLmax en función del parámetro (C2 +
1)A tan λLE , el número de Mach y el parámetro A cos λLE (1 + 4λ2 ) [5].
174
I. VARIACIÓN DEL COEFICIENTE DE SUSTENTACIÓN MÁXIMO DEL
ALA DEBIDO A DISPOSITIVOS HIPERSUSTENTADORES
Sf lap
∆(CLmax )F lap = 0,9∆(Clmax )F lap KΛ
S
Por último falta por determinar el parámetro Sf lap /S, donde S obviamente es la
superficie alar, ya calculada anteriormente en el módulo de geometría del programa,
y Sf lap es la superficie del flap, entendiendo esta no como la superficie del propio
flap sino que es la superficie del ala al disponer de un flap. A fin de ejemplificar de
una forma más exacta la definición de Sf lap se deje en la figura 55 un ejemplo de
cómo obtener esa superficie para dispositivos de borde de ataque y borde de salida.
Sf lapb.s
∆(CLmax )F lap = 0, 9∆(Clmax )F lapb.s (KΛ) b.s
S
Sf lapb.a
+0, 9∆(Clmax )F lapb.a (KΛ) b.a
S
175
I. VARIACIÓN DEL COEFICIENTE DE SUSTENTACIÓN MÁXIMO DEL
ALA DEBIDO A DISPOSITIVOS HIPERSUSTENTADORES
176
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
dCL
CLα = (82)
dα
2πA
CLα = r (83)
A2 β 2 tan (Λmax,t )2
2+ 4+ η2
1+ β2
Cl α
η=
2π/β
Por último el ángulo Λmax,t es la flecha del ala en la cuerda donde el perfil tiene
el máximo espesor y A el alargamiento o Aspect Ratio del ala.
Por otra parte, a fin de tener la pendiente de sustentación no solo del ala, también
del resto de componentes de la aeronave a diseñar, es necesario añadir a la anterior
ecuación una serie de parámetros que mejorarán la aproximación para la aeronave,
donde se tendrá en cuenta los posibles winglets y el fuselaje.
177
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
2πAef f Sexposed
CLα = r F (84)
A2 β 2
tan (Λmax,t )2
Sref
2 + 4 + η2 1 + β2
Los nuevos términos que se han añadido están relacionados con el resto de elemen-
to del avión que no se habían tenido en cuenta. La superficie Sexposed es la superficie
de referencia del ala (recordar la definición del módulo de geometría) menos la parte
cubierta por el fuselaje.
Por otro lado se encuentra el parámetro Aef f , se trata del “nuevo” valor del
Aspect Ratio o alargamiento del ala debido a lo que se conocen como Winglets o
endplates, ver figura 56; estos dispositivos son elementos que se disponen al final
de las alas a fin de mejorar el comportamiento del ala de la aeronave y reducir la
resistencia inducida generada por las estela turbillonaria; básicamente reducen la
velocidad de perturbación y con ello la resistencia que genera la sustentación del
ala.
178
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
aeropuertos.
Por otra parte supone un beneficio para el consumo de combustible, puesto que
el alcance de la aeronave aumenta al verse reducida la resistencia. Por tanto es un
elemento que la mayoría de aerolíneas quiere en sus aviones, al tener una relación
de coste/beneficio bastante elevado. Es por ello que durante los últimos años la
mayoría de aeronaves del mercado llevan este tipo de dispositivos y multitud de
estudios aerodinámicos en el diseño de estos dispositivos han ido apareciendo.
Si se trata de Winglets:
Aef f = 1, 2A (86)
Si se trata de Endplates:
179
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
Scola ∂
∆CLαCola = CLαCola 1−
S ∂α
∂
El parámetro 1 − ∂α se trata de un parámetro que trata de definir el fenómeno
de variación del flujo de la corriente. Para ello se establece lo que se conoce como
una tasa de cambio del ángulo que ve la cola de la aeronave con respecto al ángulo
de ataque de la aeronave.
Dicha tasa se puede estimar a partir de una ecuación basada en datos empíricos
de un túnel de viento [22], [7], y se escribe como:
0,25
∂ 21 C 10 − 3λ zh
= 0,725 1−
∂α A lh 7 b
180
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
Para el caso del canard se tiene una ecuación similar a la de la cola, según [7]:
SCanard ∂u
∆CLαCanard = CLαCanard 1+
S ∂α
Donde en este caso u es el ángulo de la corriente aguas arriba del ala, y por
tanto la tasa de variación de este ángulo con respecto al ángulo de ataque es a partir
de una formulación empírica similar a la anterior [22] y [7]:
−(1,04+6A−1,7 )
∂u lc
= 0, 3A0,3 − 0, 33
∂α c
Donde todos los parámetros son conocidos, excepto lc que se obtiene de la figura
58 y CLαCanard que es la pendiente de la curva de sustentación de canard, que se
calculará de la misma forma que la calculada en el caso del ala.
2πAef f Sexposed Scola ∂
(CLα )T otal = r F+ CLαCola 1−
A2 β 2
tan (Λmax,t )2
Sref S ∂α
2 + 4 + η2 1 + β2
SCanard ∂u
+ CLαCanard 1+ (88)
S ∂α
En el caso de que se tengan flaps desplegados, configuración sucia, existe una co-
rrección que permite obtener la variación de la pendiente de la curva de sustentación
que se genera por desplegar los flaps. La fórmula se obtiene de [3]:
∆f CL0 c0
cf 2
(CLα )Sucia = (CLα )Limpia 1+ 1 − 0 sin δf − 1 (89)
∆f Cl0 c c
181
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
CLδ
∆f CL0 CLα CLα
= Kb C
∆f Cl0 Cl α lδ
Clα
C
Donde CLl α es un parámetro conocido, Kb es un factor que mide la efectividad
α
del flap en función de la envergadura; dicho parámetro se puede obtener a partir de
la figura 59.
Figura 59: Variación del parámetro Kb en función del ratio de envergaduras bf /(b/2)
y el taper ratio [20].
El parámetro que aparece en la figura 59, se puede apreciar en la figura 55, donde
el valor de bf es la diferencia entre bf0 y bfi , de este modo habrá que evaluar dos
puntos en la gráfica de Kb , lo que quedaría como:
Kb = Kb 1 − Kb 2
Donde Kb1 está evaluado en bf0 y Kb2 evaluado en bfi , con lo que la ecuación
recursiva de la gráfica obtenida en Excel queda como:
182
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
!4 !3 !2 !
bf0 bf0 bf0 bf0
Kb = −0, 7954 b
+ 1, 1101 b
− 0, 7685 b
+ 1, 4425 b
2 2 2 2
!4 !3 !2 !
bfi bf i bfi bf i
− −0, 7954 b
+ 1, 1101 b
− 0, 7685 b
+ 1, 4425 b
2 2 2 2
CL
δ
CL
Por último el parámetro α es nombrado el factor de cuerda del flap que
Cl
δ
Cl
α
también se puede escribir según [3] como Kc , este parámetro se puede calcular a
partir de la gráfica de la figura 60; sin embargo, hay que hallar primero el valor de
(aδ )Cl a fin de saber que curva escoger de la gráfica de la figura 60.
(aδ )Cl = 0, 6
183
J. CÁLCULO DE LA VARIACIÓN DE LA PENDIENTE DE LA CURVA DE
SUSTENTACIÓN
CLδ
CLα
= 1, 1924A−0,067
Clδ
Clα
2πAef f Sexposed
(CLα )Sucia = r F
A2 β 2
tan (Λmax,t )2
Sref
2 + 4 + η2 1 + β2
0
CLα c cf 2
. 1+ Kb Kc 1 − 0 sin δf − 1 + ∆CLαCanard + ∆CLαCola (90)
Cl α c c
184
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
P
k Cfk IF k F Fk Swetk
(CDmin )T otal = + CDT ren aterrizaje + CDP rotuberancias + CDSoporte alas
S
+CDV entana cabina + CDf lotadores
Dónde CDT ren aterrizaje corresponde al tren de aterrizaje, se debe tener en cuenta
que esta aportación solo aparecerá en la polar del avión en las actuaciones corres-
pondientes al aterrizaje y despegue.
La resistencia CDSoporte alas se añade en el caso de que las alas tengan un soporte
que las una al fuselaje; dicho soporte suele encontrarse en aeronaves con ala alta, y
que mejoran de esa forma la resistencia estructural de la propia ala, sin embargo en
detrimento de esto, habrá que añadir la resistencia que este elemento genera.
185
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
En el caso de CDV entana cabina , este término hace referencia a la parte correspon-
diente al parabrisas de la cabina de pilotos, ya que suele ocurrir que la forma del
parabrisas suele romper con la suavidad de la superficie del propio fuselaje, y por
tanto se genera una resistencia que es necesario cuantificar.
1, 328 0, 455
Cfk = %Laminark √ + %T urbulentok (91)
Re (log10 Re) (1 + 0, 144M 2 )0,65
2,58
186
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
Por tanto, una vez definidas cada uno de los parámetros de la componente, se
pasará a desarrollarlos para cada una de las componentes k de la aeronave, que para
el programa serán un total de siete y que aparecen como subíndices de la ecuación
18:
Donde:
1, 328 0, 455
CfAla = 0, 2 q + 0, 8 2,58
ρV C
µ
(log10 ρV C
µ
) (1 + 0, 144M 2 )0,65
IF Ala = 1, 1
" 4 #
0, 6 t t
1, 34M 0,18 (cos Λ)tmax 0,28
F FAla = 1+ + 100
(x/c)tmax c c
187
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
Donde:
1, 328 0, 455
CfHorizontal = 0, 2 q + 0, 8 2,58
ρV C h
µ
(log10 ρV C h
µ
) (1 + 0, 144M 2 )0,65
IF Horizontal = 1, 05
" 4#
0, 6 t t
1, 34M 0,18 (cos Λ)tmaxh 0,28
F FHorizontal = 1+ + 100
(x/c)tmaxh c h c h
SwetHorizontal = 2Sh
Donde:
1, 328 0, 455
CfV ertical = 0, 2 q + 0, 8 2,58
ρV C v
µ
(log10 ρV C v
µ
) (1 + 0, 144M 2 )0,65
IF V ertical = 1, 05
188
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
" 4#
0, 6 t t
1, 34M 0,18 (cos Λ)tmaxv 0,28
F FV ertical = 1+ + 100
(x/c)tmaxv c v c v
Donde:
1, 328 0, 455
CfCanard = 0, 3 q + 0, 7 2,58
ρV C c
µ
(log10 ρV C c
µ
) (1 + 0, 144M 2 )0,65
IF Canard = 1, 05
" 4#
0, 6 t t
1, 34M 0,18 (cos Λ)tmaxc 0,28
F FCanard = 1+ + 100
(x/c)tmaxc c c c c
SwetCanard = 2Sc
Donde:
1, 328 0, 455
CfF uselaje = 0, 15 q + 0, 85 2,58
ρV lf us
ρV lf us
µ
(log10 µ
) (1 + 0, 144M 2 )0,65
Siendo lf us la longitud del fuselaje sin tener en cuenta la parte del fuselaje de la
cola.
189
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
IF F uselaje = 1
60 ff us 60 ff us
F FF uselaje = 1 + + + 0, 4 1 + + %Rectangular
ff us 3 400 ff us 3 400
Se deja a modo de ejemplo la siguiente imagen, ver figura 63, donde se aprecia
la geometría de esta componente.
Donde:
1, 328 0, 455
CfN acelle = 0, 3 q + 0, 7 2,58
ρV lN
µ
(log10 ρVµlN ) (1 + 0, 144M 2 )0,65
IF N acelle = 1, 3
0, 35
F FN acelle = 1 +
fN acelle
190
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
Se deja a modo de ejemplo la siguiente imagen, ver figura 63, donde se aprecia
la geometría de esta componente.
Donde:
1, 328 0, 455
CfP ylon = 0, 3 q + 0, 7 2,58
ρV lp
ρV lp
µ
(log10 µ
) (1 + 0, 144M 2 )0,65
IF P ylon = 1, 5
0, 35
F FP ylon = 1+
fP ylon
191
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
CDT ren aterrizaje = CDT ren aterrizajeT rasero + CDT ren aterrizajeDelantero (92)
192
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
1
∆CDsT rasero = 0, 05328e5,615SA /(dw) → Con carenado
2
Dónde SA es la superficie de las barras soporte del tren de aterrizaje (ver figura
64)
• Sin flaps:
0, 003294W 0,785
CDT ren aterrizajeT rasero =
jet Sref
• Con flaps:
0, 001792W 0,785
CDT ren aterrizajeT rasero =
jet−f laps Sref
193
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
dw
CDT ren aterrizajeDelantero = ∆CDsDelantero
Sref
Donde para calcular el término ∆CDsDelantero habrá que utilizar una gráfica de
[25] donde aparecen una serie de curvas que establecen la variación de esta resistencia
en función de dos parámetros nuevos, ver figura 65, a distancia del morro del avión
hasta el tren delantero y e que es la altura del suelo a la parte de abajo del fuselaje.
e 2 e
∆CDsDelantero = 0, 0538 + 0, 0372 + 0, 1859
d d
e 4 e 3 e 2 e
∆CDsDelantero = 0, 043 − 0, 407 + 1, 4176 − 1, 8748 + 1, 3147
d d d d
Para a/d ≥ 2, 5:
e 2 e
∆CDsDelantero = 0, 1066 − 0, 0751 + 0, 3584
d d
194
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
" 2 #
t t Lc
CDSoporte alas = 2CfSoporte alas 1 + +
c c Sref
1, 328 0, 455
CfSoporte alas = 0, 4 q + 0, 6 2,58
ρV c
µ
(log10 ρVµ c ) (1 + 0, 144M 2 )0,65
Por su parte, el parámetro CDV entana cabina también se ha obtenido de [6] con lo
que queda como:
Amax
CDV entana cabina = ∆CDscabina
Sref
Finalmente, en el caso de diseñar una aeronave tipo anfibia que requiera del uso
de flotadores.
Se calcula la resistencia de los mismos según tres tipos de flotadores de los cuales
se obtuvieron las formulas empíricas, ver figura 66, de [26].
Dichas formulas depende del ángulo de ataque, sin embargo los términos que
multiplican al ángulo de ataque son muy pequeños, y a fin de tener un valor único
de resistencia se establecerá como un parámetro constante, de forma que, para cada
uno de los flotadores (A, B y C) se tiene que:
195
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
Amax
CDF lotadorA = N 0, 2142
Sref
Amax
CDF lotadorB = N 0, 1771
Sref
Amax
CDF lotadorC = N 0, 1981
Sref
Figura 66: Flotadores tipo A, B y C de los cuales se sacaron los datos para obtener
la resistencia [6].
196
K. CÁLCULO DEL COEFICIENTE DE RESISTENCIA MÍNIMA DE LA
POLAR DE LA AERONAVE
197
L. CÁLCULO DE LA CONSTANTE DE RESISTENCIA INDUCIDA
e = 1, 78 1 − 0, 045A0,68 − 0, 64
Por otra parte, a fin de tener en cuenta otros parámetros que afecta a la re-
sistencia inducida, se puede establecer a mayores el añadir Winglets o Endplates,
sustituyendo el Aspect Ratio, por el Aspect Ratio efectivo definido anteriormente,
ver ecuaciones 86 y 87.
33(h/b)1,5
Kef ecto suelo =K (93)
1 + 33(h/b)1,5
199
M. DETERMINACIÓN DE LAS COMPONENTES ∆1 Y ∆2 EN LA
CONTRIBUCIÓN DE LA RESISTENCIA
M.1. Cálculo de ∆1
Tabla 16: Variación del parámetro ∆1 en función del espesor relativo y del ratio de
cuerdas.
Se puede establecer que el valor de ∆1 haciendo una media para cada caso es:
201
M. DETERMINACIÓN DE LAS COMPONENTES ∆1 Y ∆2 EN LA
CONTRIBUCIÓN DE LA RESISTENCIA
M.2. Cálculo de ∆2
Para el caso de ∆2 es más complicado obtener un valor fijo, las gráficas en este
caso sí que difieren, ver figura 68.
Además cada una de las curvas se encuentra en función de la deflexión del flap,
el cual variará en función de la maniobra en la que se encuentra la aeronave, sin
embargo siguiendo lo establecido en la tabla 14, se elegirán los ángulos máximos de
deflexión, para luego establecer un porcentaje en el caso de que el flap no se deflecte
completamente.
Tabla 17: Variación del parámetro ∆2 en función del espesor relativo y del ratio de
cuerdas.
202
M. DETERMINACIÓN DE LAS COMPONENTES ∆1 Y ∆2 EN LA
CONTRIBUCIÓN DE LA RESISTENCIA
Figura 67: Variación del parámetro ∆1 en función del ratio de cuerdas, para flaps
tipo slotted plain y Split [6].
Se establece por tanto que para el parámetro ∆2 , realizando una media de cada
caso:
Slotted flap
Sf lap
∆CDf lap Slotted
= 0, 0542
Sref
Split flap
Sf lap
∆CDf lap Split
= 0, 156
Sref
Slotted flap
Sf lap
∆CDf lap plain
= 0, 16
Sref
Para el caso del flap tipo Fowler al no tener referencia con el anterior método,
se puede utilizar el método seguido por Corke [28], donde se tiene un valor fijo para
dos tipos de deflexiones.
203
M. DETERMINACIÓN DE LAS COMPONENTES ∆1 Y ∆2 EN LA
CONTRIBUCIÓN DE LA RESISTENCIA
Figura 68: Variación del parámetro ∆2 en función del ángulo de deflexión del flap y
del espesor relativo del perfil [6].
Por otro lado a fin de tener una idea de si la aproximación anterior de Young
se asemeja a la comentada y así poder validar el uso de la estimación de los flaps
tipo Fowler, se expondrán el resto de componentes estimadas por Corke para las
deflexiones calculadas, (se supondrá como valor de Sf lap /Sref = 0, 5 para poder
comparar los resultados):
1. Slotted flaps
∆CDf lap SlottedY oung
= 0, 0271
∆CDf lap SlottedCorke
= 0, 03
2. Split flaps
∆CDf lap SplitY oung
= 0, 078
∆CDf lap SplitCorke
= 0, 1
Los valores difieren en mayor medida para este caso. Se debe tener en cuenta
que en el caso de Corke el valor es fijo, sin embargo en el método aproximado
al de Young depende fundamentalmente de las superficies de los flaps, y esa
será la diferencia entre ambos métodos.
204
M. DETERMINACIÓN DE LAS COMPONENTES ∆1 Y ∆2 EN LA
CONTRIBUCIÓN DE LA RESISTENCIA
3. Plain flaps
∆CDf lap P lainY oung
= 0, 08
∆CDf lap P lainCorke
= 0, 1
4. Fowler flaps
Por último se dejarán los valores según Corke para los Fowler flaps y que serán
los que se utilizarán al ser los únicos que se tienen:
Para deflexión de 30o
∆CDf lap F owlerCorke
= 0, 032
∆CDf lap F owlerCorke
= 0, 083
205
N. CÁLCULO DE LA ENVOLVENTE DE VUELO
Los datos que debe establecer, al igual que la forma de establecerlos, aparece de
forma detallada en el manual del programa, en el Anexo Q del presente proyecto.
207
N. CÁLCULO DE LA ENVOLVENTE DE VUELO
Donde CLT otalf lapsmax es el coeficiente de sustentación máximo con los flaps
extendidos al máximo, lo que implica la fase de aterrizaje de la aeronave
• Velocidad de pérdida con flaps Vs1
s
2W0
VS1 = S
ρ0 CLf lapsmax
√
VA = VS n1
VD > Vcmin
VDF ≤ Vcmin
208
N. CÁLCULO DE LA ENVOLVENTE DE VUELO
Velocidad no sobrepasable VN E
Se debe cumplir para esta velocidad que:
(
≥ 0, 9VH
VN E p
≥ 4, 77 W0 /S
Tabla 18: Valores típicos según aeronave, de los factores de carga límites, [1]. Tabla
incorporada en el Excel TablaEstructuras.xlsx.
U
≈
∆α = arctan U/V |{z}
V
U/V 1
209
N. CÁLCULO DE LA ENVOLVENTE DE VUELO
∆L ρV U CLα
∆n = =
W 2W/S
Kg U V CLα
n3/4 = 1 ±
498W/S
Donde:
νg 0, 88
Kg =
5, 3 + νg
2W/S
νg =
CgCLα ρ
Agrupando todos los términos, se tiene que el factor de carga debido a ráfagas
es:
2W/S
0,88
CgCL ρ
α
2W/S
5,3+
CgCL ρ
n3/4 = 1 ± α
498W/S
210
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
Ñ.1. Ala
Para el caso del ala de la aeronave se pueden establecer las siguientes ecuaciones:
0,6 −0,3
AAla 100t/c
WALAR = 0, 036SAla 0,768
q 0,006 0,04
λ (nz W0 )0,49
cos Λc/4 2 cos Λc/4
WALAR = 0, 0051SAla 0,649 (t/c)−0,4 (1 + λ)0,1 cos Λc/4 −1 (nz W0 )0,557 Aala 0,5 Sf laps 0,1
Dónde:
nz es el factor limite último que según Raymer [1], se puede definir como el factor
límite máximo multiplicado por 1,5.
λ es el tapper ratio.
211
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
0,043 −0,12
AHT 100t/c
WHTR = 0, 016SHT 0,896
q 0,168
λHT −0,02
(nz W0 )0,414
cos Λc/4HT 2 cos Λc/4HT
212
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
−0,25
Fw
WHTR = 0, 00379Kuht 1+ W0 0,639 nz 0,1 SHT 0,75 lHT −1 (0, 3lHT )0,704
bHT
0,1
−1 0,166 Se
(cos ΛHT ) AHT 1+
SHT
Dónde:
213
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
0,357 −0,49
0,873 AV T 0,122 0,039 100t/c
WV TR = 0, 073SV T q λHT
cos Λc/4V T 2 cos Λc/4V T
(nz W0 )0,376 (1 + 0, 2Fcola )
−0,726
WV TN = 0, 19[(nz W0 )0,363 (1 + Fcola )0,5 SV T 1,089 Mo 0,601 lV T
0,217
Sr
1 + SV T AV T 0,337 (1 + λV T )0,363 cos ΛV T −0,484 ]1,014
Dónde:
214
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
Ñ.4. Fuselaje
−0,072
lF U S
WF U SR = 0, 052SF U S 1,086
lF −0,051 (nz W0 )0,177 q 0,241 +11, 9(Vp ∆P )0,271
dF U S
0,71
−2 0,283 −3 0,95 lF
WF U SN = 10, 43(q10 ) W0 10
wF
0,1
0,5 lF U S
WF U SR = 0, 3280kdoor KLg (W0 nz ) lF 0,25
SF U S 0,302
(1 + Kws )0,04
dF U S
Dónde:
215
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
0,409 0,845
0,768 lm 566 ln
WT renAterrizajeR = 0, 095(nl W0 ) + 0, 125(nl W0 )
12 12
0,601
0,684 lm
WT renAterrizajeN = 0, 054(nl W0 )
12
Para aviones comerciales y militares, en régimen subsónico:
0,84
WT renAterrizajeN = 62, 1 W0 10−3
216
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
WT renAterrizajeR = 0, 0106Kmp lm 0,4 Nmw 0,321 Nmss −0,5 W0 0,888 nl 0,25 VStall 0,1
+0, 032Knp ln 0,5 Nnw 0,45 W0 0,646 nl 0,2
Dónde:
Kmp parámetro de valor 1,126 para tren de aterrizaje tipo “Kneeling”, capacidad
de bajar la altura para facilitar la carga de la aeronave; si no tiene valor 1.
Ñ.6. Motorización
Donde:
WM otor es el peso del motor sin instalar; dicho peso se obtiene de otras ecuaciones,
217
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
que dependerán del tipo de motor escogido, que se obtienen de un listado de un gran
número de motores [6] y que se plantearán a continuación:
P otencia − 21, 55
WM otor =
0, 5515
Motor turbopropulsor (tipo hélice):
P otencia − 110, 7
WM otor =
2, 631
Motor turbofan (tipo jet):
Empuje − 153, 6
WM otor =
5, 646
0,8
1,536 0,371 nz W0
WCT RLR = 0, 053lF U S b
104
WHY DR = 0, 001W0
218
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
0,576
WHY DN = 56, 01 W0 q10−5
WAVN = Npil 15 + 0, 032W0 10−3 +NM otor 7, 4 + 0, 046W0 10−3 +0, 15W0 10−3
219
Ñ. ECUACIONES DEL CÁLCULO DE PESO EN VACÍO POR COMPONENTE
Dónde:
Ñ.12. Mobiliario
WM obR = 0, 0582W0 − 65
Dónde:
220
O. ECUACIONES DE MECÁNICA DEL VUELO
Eje horizontal X:
T −D =0 (94)
Eje horizontal Y :
L−W =0 (95)
Donde:
W es el peso de la nave.
Por otra parte, se conoce de aerodinámica que las fuerzas D y L se pueden definir
por:
1
L = CL ρSV 2
2
1
D = CD ρSV 2
2
221
O. ECUACIONES DE MECÁNICA DEL VUELO
Donde:
Para establecer relaciones más simples entre las variables y, sobre todo, para
determinar valores empíricos se puede relacionar la sustentación con la resistencia,
obteniendo lo que se conoce como la eficiencia aerodinámica:
L CL
=
D CD
1 1
T = D = CD ρSV 2 = CD0 − K2 CL + K1 CL 2 ρSV 2
2 2
222
O. ECUACIONES DE MECÁNICA DEL VUELO
T 1 1
CD0 − K2 CL + K1 CL 2 ρ V2
=
W 2 (W/S)
q = 1/2ρV 2
Por otra parte, teniendo en cuenta que la actuación seguida por la aeronave en
este apartado es la de crucero, teniendo en cuenta la ecuación 95, se tiene que:
L = W = qCL S
T 1 K1
= CD0 q − K2 + (W/S) (96)
W (W/S) q
223
O. ECUACIONES DE MECÁNICA DEL VUELO
αTSL 1 K1
= CD0 q + (β(WT O /S)) − K2
βWT O β(WT O /S) q
TSL β 1 K1
= CD0 q + (β(WT O /S)) − K2 (97)
WT O α β(WT O /S) q
O.2. Despegue
224
O. ECUACIONES DE MECÁNICA DEL VUELO
Cada una de las fases anteriores presenta distintas expresiones y algunas de ellas
distintas variantes para su análisis y estudio.
En primer lugar es necesario resaltar que las ecuaciones que se presentan durante
las fases del despegue son de tipo transitorio debido a la existencia de variaciones
de la velocidad y de la altura de la aeronave.
Por tanto, las expresiones estacionarias que se obtengan representan valores me-
dios o valores de carácter empírico. Las expresiones básicas durante las fases del
despegue se tienen por:
Eje horizontal X:
dVx
Tx − D − Rf = m (98)
dt
dx
Vx =
dt
Eje horizontal Y :
dVy
Ty + L − W = m (99)
dt
dy
Vy =
dt
225
O. ECUACIONES DE MECÁNICA DEL VUELO
Donde:
m es la masa de la aeronave
V es la velocidad VT O ya mencionada
O.2.1. Rodadura
Vy = 0, y = 0; Vx = V , x = s, Tx = T ; Ty = 0.
dV
T − D − Rf = m
dt
ds
V =
dt
Se realizará para esta actuación dos tipos de aproximaciones, a fin de ver qué
ecuación utilizar en el programa de MATLAB.
1a Aproximación
226
O. ECUACIONES DE MECÁNICA DEL VUELO
dV T
T =m →V = t
dt m
ds 1 T 2 1m 2
V = →s= t = V
dt 2m 2T
2 W 2 T 1 VT O 2
SG = VT O → = (100)
2g T W 2g SG
2
2 VStall 2 W 2
VT2O = KT O = KT O (101)
2 S ρCLmax
Y de sustituir en 100:
KT O 2
T W
= (102)
W S gSG CLmax
Y en consecuencia:
KT O 2
W
SG = (103)
S (T /W )gCLmax
Con lo que se puede ver que la expresión anterior establece que, si SG está defi-
nida, la relación entre (T /W ) respecto a (W/S) es una recta.
2a Aproximación
227
O. ECUACIONES DE MECÁNICA DEL VUELO
T − D − Rf
V = t
m
1 T − D − Rf 2 1 m 2
s= t = V
2 m 2T
" #
1 1−ϕ
SG = − (W/S) ln C
(104)
ϕρg (T /W − µT O ) kLmax
2
TO
Donde:
1
k=
ϕρg
CLmax
k0 =
KT O 2
Resulta:
228
O. ECUACIONES DE MECÁNICA DEL VUELO
T 1 − ϕ k(W/S)
SG
= e + µT O (105)
W k0
En este caso (T /W ) también se observa que crece con (W/S) aunque no lineal-
mente.
Vy = 0, y = 0; Ty = 0
SR = VT O tR
s
W 2
SR = tR KT O (106)
S ρCLmax
229
O. ECUACIONES DE MECÁNICA DEL VUELO
s
W 1 W 2
SG + SR = + tR KT O (107)
S (T /W )gCLmax S ρCLmax
Para el primer caso a), se considera que la aeronave realiza una trayectoria cir-
cular con un radio RC .
230
O. ECUACIONES DE MECÁNICA DEL VUELO
dV
T − D − R − Wv = m
dt
dδ
V =
dt
V2
(L − Wq ) = m
RC
dθ dθ
w= → V RC =
dt dt
δR = RC θT R = VT R tT R (108)
Siendo Wq = W cos θ
V2 W 2 1 W 2 1
(L − Wq ) = m → RC = V = V
RC g L − Wq g L − W cos θ
1 1
RC ≈ V 2 (109)
g n−1
231
O. ECUACIONES DE MECÁNICA DEL VUELO
V2
δR = RC θT R = θT R (110)
(n − 1)g
Se debe tener en cuenta que en la fase de despegue los flaps no estarán extendidos
al máximo, tal y como se explicó en el módulo de aerodinámica; el porcentaje de
deflexión y por tanto el valor del coeficiente de sustentación máximo dependerá del
tipo de flaps, en función de la tabla 14.
Eje horizontal X:
T − D = W sin γ (111)
Eje horizontal Y :
L − W cos γ = 0 (112)
232
O. ECUACIONES DE MECÁNICA DEL VUELO
Por otra parte, al igual que se ha descrito en el apartado del crucero, se tienen las
relaciones de la resistencia con respecto a la sustentación (polar del avión), ver ecua-
ción 17, y sabiendo que el ascenso es con velocidad de subida constante, se tiene que:
D
= CD0 − K2 CL + K1 CL 2
qs
dV
= 0 → γ = cte
dt
Por último, tal y como describe el libro de referencia [4], en la subida de la aero-
nave se puede suponer que γ 1, por lo que cos γ ≈ 1, y en consecuencia se tiene
que L ≈ W , misma expresión que se obtenía con el crucero; de este modo realizando
el mismo procedimiento seguido en el apartado de la actuación de crucero se tiene
la siguiente ecuación, que relaciona el empuje-peso con la carga alar de la aeronave:
T 1 1 2 1
= CD0 q − K2 q + sin γ + CL + K1 qCL (113)
W (W/S) (W/S) (W/S)
TSL β 1 K1
= CD0 q + (β(WT O /S)) − K2 + sin γ (114)
WT O α β(WT O /S) q
El valor de α que varía con el tiempo, se obtiene de una gráfica del libro de
referencia [1]; donde se establece una serie de curvas de como varía el empuje con la
altura en función del tipo de motor de la aeronave, y que se han incorporado en el
programa de MATLAB obteniendo una ecuación regresiva de dichas curvas.
233
O. ECUACIONES DE MECÁNICA DEL VUELO
Como ya es sabido, una aeronave para realizar un giro tendrá que utilizar los
alerones dispuestos en las alas, por lo que al realizar dicho giro, la aeronave gira con
respecto al eje del fuselaje, teniendo en ese caso las alas un ángulo con respecto a
un plano paralelo al suelo, denominado ángulo de alabeo (bank angle) φ tal y como
aparece representado en la figura 70.
Eje horizontal X:
Eje horizontal Y :
L − Wn = 0 (116)
s 2
ΩV
n= 1+ (117)
g
234
O. ECUACIONES DE MECÁNICA DEL VUELO
o bien:
s 2
V
n= 1+ (118)
gRg
Donde:
Rg es el radio de giro.
Ω velocidad angular.
g gravedad.
TSL β 1 K1 2
= CD0 q + n (β(WT O /S)) − K2 n (119)
WT O α β(WT O /S) q
O.5. Aterrizaje
235
O. ECUACIONES DE MECÁNICA DEL VUELO
s
2β
SF R = tF R KT D (120)
ρCLmax (WT O /S)
En este caso el valor del coeficiente de sustentación máximo será el calculado con
el máximo de deflexión.
Siguiendo con el desarrollo del despegue, ecuación 104, se tiene según [4] para la
parte correspondiente al frenado:
β 1+ϕ
SB = (WT O /S) ln (121)
ϕρg −αTSL
+ µ
CLmax
βWT O TO k 2
TD
" #
β 1+ϕ
SB = (WT O /S) ln C
(122)
ϕρg µT O kLmax
2
TD
236
O. ECUACIONES DE MECÁNICA DEL VUELO
Por último, según [4] se tiene que la fase de descenso se puede escribir como:
2β Kobs 2 − KT D 2
WT O CLmax 2hobs
SA = 2 2 + (123)
CD Kobs 2 + KT D 2
ρgCD Kobs + KT D S
237
P. ECUACIONES DE CONSUMO DE COMBUSTIBLE EN CADA ETAPA DE
VUELO (DISEÑO AVANZADO)
P.1. Taxi
√ α
Wf TSL
= 1 − C1 θ ∆t (124)
Wi T axi β WT O
C1 es una de las constantes del modelo del consumo del motor T SF C, ver ecua-
ción 21.
Wf
βT axi = =1
Wi Encendido
P.2. Despegue
-Aceleración despegue:
Wf
n o
C VT O
− T SF
=e g 1−u
(125)
Wi Aceleracion despegue
239
P. ECUACIONES DE CONSUMO DE COMBUSTIBLE EN CADA ETAPA DE
VUELO (DISEÑO AVANZADO)
β WT O q S
u= ϕ + µT O (126)
α TSL β WT O
Donde:
q es la presión dinámica.
Wf Wf
βAceleracion despegue =
Wi Encendido Wi T axi
-Rotación:
Wf α TSL
= 1 − T SF C ∆tR (127)
Wi Rotacion β WT O
Puesto que el cálculo de esta variable variará a partir de las etapas de vuelo, ya
que el usuario podrá modificar el número de etapas y la posición de las mismas, no
se volverá a realizar el cálculo de dicha variable en la explicación de los cálculos; si
bien el propio código de MATLAB se tiene en cuenta esta posibilidad y se aplica un
proceso de condiciones para que el cálculo del parámetro se haga de forma correcta
cualquiera que sea la la misión de vuelo.
En esta etapa, la aeronave acaba de despegar, por tanto los datos iniciales de la
etapa se tienen, de la etapa de despegue, como es el caso de la velocidad inicial en
este ascenso.
Se trata de una etapa de ascenso acelerado, por lo que los datos necesarios a
mayores serán la altura final, la velocidad final y el ángulo de ascenso. Con estos
tres valores se puede calcular la fracción de peso, a partir de la siguiente ecuación:
240
P. ECUACIONES DE CONSUMO DE COMBUSTIBLE EN CADA ETAPA DE
VUELO (DISEÑO AVANZADO)
∆(h+V 2 /2g)
T SF C
−
Wf
astd 1− CD α β WT O
( )
=e CL TSL
(128)
Wi Ascenso despegue
Donde las únicas variables no vistas hasta ahora son astd que es la velocidad del
sonido (atmosfera estandar) y ∆(h + V 2 /2g), que se trata de la variación de energía
potencial y cinética en la etapa de ascenso, con lo que dicha variación se puede
definir como:
Donde haeropuerto es la altura del aeropuerto de despegue, Vf inal yhf inal son la
velocidad y altura al finalizar la etapa respectivamente.
− T SF C ∆h
Wf
astd 1− CD β WT O
α( )
=e CL TSL
(129)
Wi Ascenso
∆V 2 /2g
T SF C
−
Wf
astd 1− CD αβ WT O
( )
=e CL TSL
(130)
Wi Crucero acelerado
P.6. Crucero
241
P. ECUACIONES DE CONSUMO DE COMBUSTIBLE EN CADA ETAPA DE
VUELO (DISEÑO AVANZADO)
( " # )
C1 0,5ρVCrucero 2 CD 2K1
WT O
0 ∆S
−M +C2 + +K2
Wf Crucero ( WST O ) ρVCrucero 2 S astd
=e (131)
Wi Crucero
P.7. Espera
√
√
Wf
= e{−(C1 +MEspera C2 ) θ( 4CD 0 K1 −K2 )tEspera }
(132)
Wi Espera
Donde en este caso, el Mach de espera no se establece como dato que aporte el
usuario, sino que la velocidad de espera se calculará a fin de obtener la óptima de
consumo mínimo.
v
u 2β WT O
MEspera =t
u
q
γPN ive lmar δ CD 0 S
K1
P.8. Descenso
Wf TSL
= 1 − T SF C tduracin descenso (133)
Wi Descenso WT O
Donde:
242
P. ECUACIONES DE CONSUMO DE COMBUSTIBLE EN CADA ETAPA DE
VUELO (DISEÑO AVANZADO)
√ C n2πN V
Wf −(C1 +Mgiro C2 ) θ D √ 2 giro
=e CL g n −1
(134)
Wi Giro
Donde:
Mgiro = Vgiro /a, siendo Vgiro la velocidad de giro y a la velocidad del sonido a la
altura donde se realiza el giro.
243
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Antes de comenzar el programa se deben tener todos los archivos en una única
carpeta. Los archivos que lo componen son los siguientes:
Archivos Excel
En lo referente a los archivos tipo Excel, los cuales se utilizarán para incorporar
datos al programa y también como una referencia donde ver los datos que se van
teniendo a lo largo del diseño de la aeronave, son los siguientes:
TablaTodosLosDatos.xlsxm
TablaGeometria.xlsx
TablaAero.xlsx
TablaEstructuras.xlsx
TablaDiseñoAvanzado.xlsx
Archivos MATLAB
TFGParte1.m
TFGModuloGeometria.m
TFGModuloAero.m
TFGModuloEstructuras.m
TFGModuloGeometria.m
245
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
TFGMecVuelo.m
TFGParte2.m
Y por otra parte aquellos archivos que se vinculan a los anteriores archivos de
MATLAB comentados, y que funcionan como funciones, bloques azules.
PasarDatosAlModuloDeGeometria.m
FuncionCalculosAerodinamicos.m
FunciónCalculoPesoEnVacio.m
Una vez explicado esto, existen una serie de características a tener en cuenta.
Es importante tener en cuenta que una vez se modifique algún archivo de Excel
es necesario guardarlo y cerrarlo antes de ejecutar el programa de MATLAB, si no
este no funcionará y dará algún tipo de error.
Q.2.1. En el Excel
Las celdas de color verde son aquellas que el usuario deberá rellenar de forma
obligatoria según los datos que requiera la aeronave de diseño.
Las celdas de color rojo son las que no se deben modificar, puesto que son datos
ya establecidos, necesarios para el buen funcionamiento del programa.
Por otra parte las celdas de color blanco se pueden modificar puesto que no se
encuentran vinculadas directamente con los cálculos del MATLAB, si bien pueden
ser datos que se han obtenido al ejecutar el programa, o datos que se aportan como
ayuda. Es por tanto elección del usuario borrarlas o modificarlas.
246
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Existe alguna celda adicional que será necesario explicar, si bien se realizará dicha
explicación a continuación, a medida que se expliquen cada uno de los archivos del
programa.
Q.2.2. En el MATLAB
Se puede dividir el programa en dos partes diferenciadas, por una parte lo que se
conoce como diseño básico, ver figura 73, y por otro el diseño avanzado, ver figura
79. La descripción conceptual de dichas partes aparece recogido en el apartado 7.2
del presente proyecto.
Se debe tener en cuenta que para llegar a la parte de diseño avanzado es necesario
llevar a cabo el diseño básico, a fin de tener un feedback de datos en el diseño.
247
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Por ello, la explicación del uso del programa se realizará a partir de las dos partes
definidas y siguiendo el orden de los archivos que el usuario deberá ir utilizando para
realizar el diseño de la aeronave.
Figura 73: Parte del diagrama de flujo de estructura del programa, correspondiente
al diseño básico
En la primera hoja de dicho Excel, ver figura 74, aparece recogida una tabla
donde aparecen en la primera columna los tipos de aeronaves que podrá elegir el
usuario según el tipo de aeronave que vaya a utilizar y una segunda columna con
el indicador o número asociado a cada aeronave. Véase de ejemplo, si lo que se
quiere es diseñar un aeronave de pasajeros el usuario podría elegir una aeronave
tipo Transport jet o bien si quiere que su aeronave sea de hélices podría utilizar
Twin engine.
Tal y como se puede observar las aeronaves tipo hélice aparecen en una celda
azul y las tipo jet sin colorear.
248
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
a cada tipo de aeronave. Los otros datos que aparecen en celdas sin colorear se
pueden modificar, siempre y cuando el usuario tenga un referencia más exacta de
dichos datos para su aeronave específica de diseño.
La Hoja4, ver figura 75, se trata de una serie de tablas donde se guardan una
gran cantidad de datos para cada tipo de aeronave, a fin de que el usuario tenga
una referencia de los valores que suelen aparecer para cada tipo de aeronave. Dichos
datos aparecen recogidos en libros de referencia [1] [2] y [3].
El recuadro de color gris sirve para trasladar los datos que el usuario haya in-
corporado en la Hoja3 y se pasen a la Hoja4, en la fila correspondiente al tipo de
aeronave elegida y de este modo tener estos datos guardados para futuros cambios
en algunas de esas especificaciones. Dicho recuadro gris no funcionará si no se tienen
activadas las macros en Excel.
La Hoja3 se trata de la primera tabla que deberá rellenar el usuario, ver figura
76.
Por un lado para el cálculo de pesos (peso total, peso en vacío, peso de com-
bustible) se necesitan datos como es la elección del tipo de aeronave, incorporando
el número de indicador de la aeronave, además de la velocidad de crucero de la
aeronave y de espera y algunos pesos de diseño.
Por otro lado, se encuentran los datos necesarios para el cálculo de la relación
empuje peso T /W y de la carga alar W/S.
Tal y como se puede ver muchos de estos datos son estimaciones, puesto que en
este momento el usuario se encuentra en la fase inicial de diseño, es por ello que se
ha implementado la Hoja4, para que el usuario tenga una referencia de qué datos
poner.
249
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Por último, la Hoja2 contiene una tabla para incorporar la misión de vuelo de
la aeronave, es decir las etapas de vuelo que seguirá la aeronave de diseño, véase la
figura 77.
Para determinar la misión de vuelo se establecen dos partes, por un lado las
casillas de color rojo corresponden a filas que no se deben modificar, puesto que
estas son fijas para cualquier tipo de misión de vuelo, el usuario solo deberá modificar
aquellas filas que tienen las casillas de color verde con el nombre VUELO.
250
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Cada una de esas cinco etapas tiene asociadas unos valores que se deben rellenar
en la columna correspondiente, tal y como se puede ver en la figura 77 y deben estar
colocadas entre las fases inicial y final que aparecen en color rojo.
De este modo se tiene ascenso (no es necesario incorporar ningún dato), descenso
(no es necesario incorporar ningún dato), crucero (distancia en millas náuticas),
espera (tiempo en horas), carga (valor de carga positivo en libras) y descarga (valor
de carga negativo en libras).
Véase un ejemplo, en el caso de la figura 77, se puede ver una misión tipo,
donde se utilizan los cinco tipos de etapas de Vuelo, si por ejemplo se quiere realizar
una misión más simple, donde la aeronave solo realice un crucero, una espera y un
descenso bastaría con eliminar el resto de etapas, con lo que el resultado que se
tendrá es el que se puede ver en la figura 78.
Tal y como se puede ver, para establecer la misión solo se deben tener nociones
básicas de Excel, y tener en cuenta que se deben modificar las celdas de color verde.
Por otra parte el número de etapas de la segunda columna también se deben de
variar según el orden en el que se vayan realizando.
251
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Se tienen ahora valores estimados de peso total W0 , peso en vacío We , peso del
combustible utilizado en la misión Wf , además de los valores de T /W y W/S. Con
todo ello es posible calcular la superficie de referencia de la aeronave S y el empuje
estimado que necesitará la aeronave; dichos valores no se tienen en primera estancia,
sino que mediante una función generada en MATLAB se pasan los datos obtenidos
al resto de ficheros Excel que corresponden al diseño avanzado y que se explicarán
a continuación.
S: Superficie de referencia
252
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Figura 79: Parte del diagrama de flujo de estructura del programa, correspondiente
al diseño avanzado
Si bien, tal y como se comentó anteriormente, algunos de los datos que se habían
calculado en el diseño básico ya están incorporadas en los ficheros de Excel que se
describirán a continuación.
Q.3.2.1 Geometría
La Hoja1 sirve para determinar medidas del fuselaje de la aeronave, las tablas
que se tienen en esta hoja se puede observar en la figura 80.
Siguiendo con la forma de proceder, se deben rellenar las celdas de color verde,
que aparecen recogida en la tabla que se puede ver a continuación, ver figura 81.
253
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Se rellanarán los valores de las celdas de color verde, en el caso de que la ae-
ronave de diseño no presente canard, el usuario deberá poner ceros en las celdas
correspondientes a esta parte de la aeronave.
Por último, la Hoja3, ver figura 83, corresponde al cálculo del tren de aterrizaje,
de góndola y pylon de los motores de la aeronave, del soporte de alas, en el caso
254
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
de que por ejemplo, la aeronave de diseño tenga alas altas y requiera de algún tipo
de soporte de las mismas y también de flotadores, en el caso de que la aeronave sea
tipo anfibio.
Q.3.2.2 Aerodinámica
Como segunda parte del diseño avanzado, el usuario deberá abrir el Excel Tabla-
Aero.xlsx, que consta de tres Hojas.
La Hoja3 del Excel de aerodinámica contiene los datos que se han obtenido de
la parte correspondiente a geometría descrita anteriormente, ver figura 84.
255
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
La Hoja1, ver figura 85, contiene datos sobre ciertas variables de los flaps que
tendrá la aeronave, en este caso al igual que ocurre con la Hoja3 no es necesario
cumplimentar ningún tipo de dato por parte del usuario.
Otro de los valores que se debe tener en consideración son los de las variables
Clideal y Clmax , con los que se podrá elegir el tipo de perfil que debe llevar el ala de
la aeronave, y que será necesario para cumplimentar los datos de la Hoja2 que se
explicarán a continuación.
Por último, también aparece recogida una tabla de parámetros de decisión, que
se rellenará automáticamente al ejecutar el archivo de MATLAB más adelante. De
esta forma el usuario tendrá en disposición los valores de cada una de las decisiones
que ha realizado, al igual que se explicó anteriormente, ver figura 85.
Tal y como se puede observar, solo será necesario rellenar la parte correspondiente
al perfil del ala, de la cola vertical, horizontal y del canard (se rellanará con ceros si
la aeronave no tiene). Para cumplimentar estos datos el usuario deberá elegir cada
uno de los perfiles, y obtener los datos que se piden.
256
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
257
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Q.3.2.3 Estructuras
De este modo se tiene en la primera hoja, ver figura 87, una serie de tablas
relacionadas con el cálculo de la gráfica de envolvente de vuelo, y que se basa en
la estimación de velocidades limites, establecidas según estándares de diseño de
aeronaves.
En el caso de la Hoja2, las tablas que incorpora esta hoja están relacionadas con
especificaciones de las componentes de la aeronave, con la finalidad del cálculo del
peso en vacío, para ello el usuario deberá incorporar ciertos valores a las variables
que aparecen en casillas de color verde, ver figura 88.
También aparecen ciertos datos necesarios para realizar dicha gráfica, como son
los coeficientes de sustentación máximos.
Por último, la Hoja4 presenta los valores de pesos calculados, ver figura 90 como
ejemplo de resultados.
258
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
259
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Es por ello que lo primero que debe realizar el usuario para continuar el proceso
de diseño en la herramienta informática será abrir el Excel TablaDiseñoAvanzado.xlsx,
Excel consta de tres hojas.
La primera Hoja de dicho archivo de Excel contiene una serie de tablas, que
son necesarias para realizar el cálculo de W/S y de T /W , es decir para realizar los
cálculos de las ecuaciones de mecánica del vuelo. Será en esta hoja, ver figura 91,
donde el usuario deberá rellenar las casillas de color verde, excepto la tabla de que
aparece recogida en la figura 92 y que se rellenará mas adelante.
Dicha tabla se debe completar una vez ejecutado el archivo de MATLAB TFG-
MecVuelo.m. Se debe tener en cuenta que este programa es de los más extensos en
260
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
cuanto a código y que por tanto tardará cierto tiempo en realizar todos los cálculos.
Siguiendo las curvas obtenidas en la figura 93 se tiene que la curva azul oscuro
representa la actuación de despegue, la curva roja crucero, la curva amarilla ve-
locidad de ascenso, la verde giro sostenido, la azul clara el aterrizaje, la de color
negro supone una actuación establecida como ascenso tras despegue, a fin de ver
las limitaciones de la aeronave una vez terminado el despegue, y por último la recta
de color magenta representa la velocidad de entrada en pérdida donde la carga alar
2
de esta componente no es mas que W/S = 0, 5ρVStall CLmax , siendo la velocidad de
entrada en pérdida aquella impuesta al comienzo del diseño básico, ver tabla 1.
261
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
El usuario una vez obtenida esta gráfica, ver figura 93, deberá elegir un punto de
la gráfica por encima de las curvas obtenidas, a fin de seleccionar el punto de diseño
de la aeronave, ver figura 94.
Así por ejemplo, en el caso de la figura 94, se tiene eligiendo ese punto un valor
262
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
Una vez hecho esto se está en disposición de realizar la parte del dimensiona-
miento. Para ello se tendrá que realizar, de una forma similar al diseño básico, la
misión de vuelo.
Para establecer la misión de vuelo, el usuario deberá entrar en la hoja2 del Excel
TablaDiseñoavanzado.xlsx, ver figura 95, y completar la misión de vuelo de la misma
forma ya explicada, moviendo y eligiendo las filas que relacionan el tipo de etapa
de la misión, dejando fijas aquellas etapas de FINAL e INICIAL, que aparecen de
color rojo.
En este caso, tal y como se puede ver en el figura 95, existe un mayor número de
etapas para poder elegir, y que requieren de un mayor número de datos para quedar
definidas.
Una vez completada la misión, y rellenado las casillas de los valores de W/S y
T /W , el usuario podrá ejecutar el archivo de MATLAB TFGParte2.m, con lo cual se
calculará finalmente el valor de peso total de la aeronave, y el peso de combustible
necesario para realizar la misión de vuelo definida.
263
Q. MANUAL DE LA HERRAMIENTA INFORMÁTICA DESARROLLADA
De este modo el usuario una vez hecho todo lo explicado, ya habrá realizado
todos los procesos a fin de utilizar la herramienta de diseño, si bien esta será la
primera aproximación; se deberá continuar con el mismo proceso reiterativamente,
del modo descrito en la figura 4, hasta que lleguen a converger los valores obtenidos
de W/S y T /W y W0 en dos iteraciones consecutivas.
Se debe tener en cuenta que una vez que se ejecute el último archivo de MATLAB
TFGParte2.m, los datos calculados se registrarán automáticamente en el Excel de geo-
metría, mediante la función PasarDatosModuloDeGeomatria.m, con lo que el usuario,
en el proceso reiterativo, solo tendrá que comenzar en el diseño avanzado, realizando
la parte de diseño básico solamente la primera vez.
264