Performance del Beechcraft B60 Duke
Performance del Beechcraft B60 Duke
PRIMER ENTREGA.
Determinación de la Performance.
SACKS, Ulises. - [Link]@[Link]
TATO, Manuela. - [Link]@[Link]
TORRES, Anibal. - [Link]@[Link]
Resumen
En el presente informe, se determinan las caracterı́sticas de performance del avión Beechcraft B60
Duke. Estas caracterı́sticas corresponden a vuelo estacionario simétrico, rectos y nivelado, vuelo en
trepada y giro nivelado coordinado. Se obtienen, además, el alcance y autonomı́a para diferentes condi-
ciones de carga, como ası́ también los parámetros correspondientes a las etapas más crı́ticas de vuelo
que son el despegue y aterrizaje
Palabras claves: Beechcraft B60 Duke, Estacionario, Alcance, Autonomı́a, Despegue, Aterrizaje,
Performance.
Abstract
In this report, the performance characteristics of the Beechcraft B60 Duke aircraft are determined. These
characteristics correspond to symmetrical, straight and level hover, climb flight and coordinated level turn.
In addition, the range and autonomy are obtained for different load conditions, as well as the parameters
corresponding to the most critical flight stages, which are take-off and landing.
Keyswords: Beechcraft B60 Duke, Stationary, Range, Autonomy, Take-off, Landing, Performance.
20 de noviembre de 2020
Índice
1. Introducción. 5
3. Análisis de Performance. 5
3.1. Vuelo estacionario, recto, simétrico y nivelado. . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3.1.1. Potencia. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.1.2. Rendimiento de la hélice. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.1.3. Expresión de la envolvente de vuelo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1.4. Velocidad de pérdida. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1.5. Techo absoluto. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.2. Vuelo estacionario, recto y simétrico (vuelo en trepada). . . . . . . . . . . . . . . . . . . . . . 12
3.2.1. Tasa de ascenso (R.C.). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.2.2. Gradiente de ascenso. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.2.3. Tasa de ascenso máxima. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.2.4. Techo de servicio. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.3. Alcance y autonomı́a. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.3.1. Combustible consumido en la etapa de trepada. . . . . . . . . . . . . . . . . . . . . . 15
3.3.2. Alcance en vuelo crucero a velocidad y altura constantes. . . . . . . . . . . . . . . . . 17
3.3.3. Diagrama de carga paga. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.3.4. Autonomı́a en vuelo crucero a velocidad y altura constantes. . . . . . . . . . . . . . . 19
3.4. Giro coordinado, nivelado y estacionario. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.4.1. Factor de carga. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.4.2. Radio y velocidad de giro. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.4.3. Obtención de los parámetros correspondientes a la aeronave en estudio. . . . . . . . 21
3.5. Despegue y aterrizaje. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.5.1. Despegue. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.5.2. Aterrizaje. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
4. Conclusión. 29
Abreviaturas
∆ Alargamiento.
b Envergadura [m]
S Superficie alar [m2 ].
CD Coeficiente de resistencia.
CD0 Coeficiente de resistencia a sustentación nula.
CL Coeficiente de sustentación.
CLM AX Coeficiente de sustentación máxima.
D Resistencia aerodinámica [N].
L Fuerza de sustentación [N].
K Coeficiente de resistencia inducido.
e Factor de Ostwald [N].
n Factor de carga.
Pe Potencia en el eje del motor [W].
Pd Potencia disponible [W].
Preq Potencia requerida [W].
T Empuje del motor [N].
W Peso de la aeronave [N].
Wv Peso operativo vacio de la aeronave [kg].
WP L Peso de carga útil [kg].
WP L Peso de carga útil [kg].
Wf Peso de combustible [kg].
Wf,climb Peso de combustible en trepada [kg].
cp Consumo especı́fico de combustible [N/kW hr].
Xb , Zb Ejes X y Z del sistema de ejes cuerpo.
Xs , Zs Ejes X y Z del sistema de ejes estabilidad.
α Ángulo de incidencia de la aeronave [◦].
θ Ángulo de elevación [◦].
γ Ángulo de trayectoria [◦].
φ Ángulo de inclinación [◦].
φT Ángulo de [◦].
V Velocidad de la aeronave relativa al aire [m/s].
Vmin Velocidad mı́nima de vuelo [m/s].
Vmax Velocidad mı́nima de vuelo [m/s].
Vs Velocidad de pérdida [m/s].
VRC Velocidad de reación de ascenso [m/s].
VRh−V,max Velocidad de alcance máximo para vuelo crucero a velocidad y altura constantes [m/s].
VEh−V,max Velocidad de autonomı́a máxima para vuelo crucero a velocidad y altura constantes [m/s].
ρ Densidad del aire [kg/m3 ].
η Rendimiento de la hélice.
ζ Fracción de combustible.
Em Eficiencia aerodimámica.
habs Techo absoluto [m].
hs Techo de servicio [m].
1. Introducción.
El Beechcraft 60 Duke es un avión de ala fija bimotor fabricado por Beechcraft. El avión tiene tren de
aterrizaje retráctil y cabina presurizada. Los dos motores de pistón están sobrealimentados, y los sobre-
alimentadores también presurizan la cabina con aire sangrado.
3. Análisis de Performance.
Se considera a la aeronave como un sólido rı́gido bajo la acción de fuerzas externas de origen aero-
dinámico (sustentación y resistencia), gravitatorio (peso) y propulsivo (empuje).
Se asume que las fuerzas antes mencionadas actúan en el centro de gravedad de la aeronave.
Para la descripción del movimiento de la aeronave, se considera a la mı́sma como una masa puntual
afectada por las fuerzas externas.
De acuerdo a las hipótesis establecidas, las ecuaciones de movimiento para vuelo estacionario, recto
y simétrico con potencia (ec. 1 y 2) se obtienen a partir de diagrama de cuerpo libre de la Figura 1.
Figura 1: Defnición de fuerzas y sistemas de referencia en vuelo estacionario, simétrico, recto y nivelado.
Como se trata de un avion convencional, es válido asumir que (α + φT ) es un valor pequeño, por lo que
sin(α + φT ) ≈ 0 y cos(α + φT ) ≈ 1. Esto simplifica las ecuaciones de movimiento quedando:
T − D − W sin γ = 0 (3)
L − W cos γ = 0 (4)
La condición de vuelo nivelado establece que el ángulo de trayectoria sea nulo, es decir γ = 0. Dicha
condición simplifica aún más las ecuaciones de movimiento que quedan expresadas en las ec.5 y 6. La
Figura 2 muestra el diagrama de cuerpo libre correspondiente a esta condión de vuelo.
Figura 2: Defnición de fuerzas y sistemas de referencia en vuelo estacionario, simétrico, recto y nivelado.
T = D = Pd /V (5)
L=W (6)
3.1.1. Potencia.
Como la planta propulsora de la aeronave en estudio corresponde a dos hélices impulsadas por mo-
tores alternativos, la potencia disponible (Pd ) se asume independiente de la velocidad (V ) y se escribe en
función de la potencia en el eje (Pe ) y el rendimiento de la hélice (ηp ). Reescribiendo la ec.5 de acuerdo a
lo anteriormente expuesto y explicitando la resistencia aerodinámica (D):
1 W2
Pd = ηp · Pe = ρV 3 SCD0 + 2K (7)
2 ρV S
Potencia disponible.
Como ya se mencionó, la potencia disponible se escribe en térmisos del rendimiento de la hélice y
de la potencia en el eje del motor. A su vez, la potencia en el eje varı́a con altura de vuelo (Figura 3)
y el rendimiento de la hélice es función de la velocidad y altura de vuelo.
Como el rendimiento es función tanto de la velocidad de vuelo como de la altura, este se obtiene a
través de un método iterativo (Numerical Method for Estimation of Propeller Effciences) que se detalla a
continuación:
η · Pe
T = (9)
V
Donde:
5. Se repite el procedimiento desde el punto 2 con el nuevo valor del rendimiento ηi hasta la conver-
gencia.
La siguiente tabla muestra los resultados obtenidos del rendimiento de la hélices para las alturas y
velocidades de vuelo consideradas.
Al graficar los valores de la tabla anterior (Figura 5) se observa que la variación del rendimiento de la
hélice respecto a la altura de vuelo es mı́nima, por lo que se considera válido asumir dicha dependencia
como nula y considerar al rendimiento como función solo de la velocidad de vuelo.
A fin de obtener una expresión útil del rendimiendo en función de la velocidad, se realiza una interpo-
lación para obtener el siguiente polinomio:
ηp (V ) = aV 4 + bV 3 + cV 2 + dV + e (12)
Donde:
a = −1,97584520192429 e−9
b = 1,23988715377169 e−6
c = −0,000271045382728803
d = 0,026153527439733
e = −0,000316985463625294
Afectando por un factor de 2 a la ec.8 y reemplazandola, junto con la ec.12, en la ec.7; luego de
reordenar se obtiene la expresión general de la envolvente de vuelo:
Donde los coeficientes A, B, C, F , G e I son funciones de la altura de vuelo y están dados por las
siguientes expresiones:
A = 2 · a · Pe (h)
1
B = 2 · b · Pe (h) − ρSCD0
2
C = 2 · c · Pe (h)
F = 2 · d · Pe (h)
G = 2 · e · Pe (h)
2KW 2
I=−
ρS
La ec.13 tiene como solución 5 raı́ces, de las cuales se consideran solo las reales positivas ya que las
restantes carecen de sentido fı́sico. Dichas raices corresponden a a las velocidades Vmı́n y Vmáx .
La velocidad mı́nima (Vmı́n ) obtenida de la ec.13 no puede ser inferior a la velocidad de pérdida de la
aeronave, caso contrario se considera esta última como velocidad mı́nima en la envolvente de vuelo. La
velocidad de pérdida es función de la altura de vuelo a través de la densidad del aire (ρ) y está dada por
la siguiente expresión:
s
2W
Vs = (14)
ρSCLmáx
El techo absoluto puede definirse como la altura para la cual las velocidades mı́nima y máxima de
vuelo coinciden, constituyendo dicho punto la altura máxima de vuelo. Dicho de otra manera, es la altura
máxima para velocidad de trepada nula. La obtención de este valor se realiza por medio de un software
de cálculo, el cual entrega el siguiente valor:
Al igual que la velocidad de pérdida, este valor lı́mite se incluye en la gráfica de la envolvente de vuelo
(Figura 8).
La tasa de ascenso es el incremento, en el tiempo, de la altura de vuelo y está dada en función del
ángulo de trayectoria (γ) por la siguiente expresión:
dh
R.C. = = V sin γ (15)
dt
Introduciendo la ec.15 en la ec.3:
(T − D)
R.C. = ·V (16)
W
La ec.16 se puede escribir en función de la potencia disponible y la potencia requerida (Preq,c ) como:
(Pd − Preq,c )
R.C. = (17)
W
Recordando que la potencia disponible Pd se escribe en función del rendimiento de la hélice ηp (ec.12)
y la potencia en el eje del motor P e (ec.8), y que además la potencia requerida para el ascenso (Preq,c )
está dada por resistencia aerodinámica D para dicha condición de vuelo, la ec.17 se puede escribir como:
2W (1 − sin2 γ)1/2
CL,c = (19)
ρSV 2
Reemplazando la ec.19 en la ec.18, se llega a la siguiente expresión cuadrática cuya incógnita es sin γ:
Donde:
2KW 2 1 2KW 2
A= , B = ηp (V )Pe (h) − ρ(h)V 3 SCD0 − (21)
ρ(h)V S 2 ρ(h)V S
Al afectar los valores obtenidos de la ec.20 por la velocidad de vuelo (V ) se tiene la tasa de ascenso
R.C.(h, V ).
Este valor corresponde a aquel para el cual la expresión de la tasa de ascenso (R.C.(h, V )) presenta
un extremo absoluto. Dicho valor se obtiene igualando a cero la primera derivada de la tasa de ascenso
respecto a la velocidad de vuelo V .
d
[R.C.(h, V )] = 0 (23)
dV
De lo obtenido en la ec.23 se despeja la velocidad, que es la correspondiente a la tasa de ascenso
máxima ([Link] ), para luego reemplazarla en la ec.20 y llegar ası́ a una expresión de la tasa de ascenso
que es función de la altura de vuelo h.
Los valores de [Link] y [Link] para la aeronave en estudio se obtienen aplicando el desarrollo
antes descripto en forma iterativa con la ayuda de un software de cálculo. Dichos valores se incluyen en la
gráfica de la envolvente de vuelo (Figura 8).
Este valor se define como la altura a la cual se garantiza una tasa se ascenso máxima, determinada
por normativa, en función del tipo de planta propulsora y las condiciones operativas de la misma. Para una
planta propulsora compuesta por dos motores alternativos, ambos en condición operativa, se establece
una tasa de ascenso máxima de 100 f pm (Bibliografı́a [1]). Mientras que para la condición de un motor
inoperativo, se establece una tasa de ascenso máxima de 50 f pm (Bibliografı́a [7]).
El consumo de combustible Wf,climb durante la etapa de vuelo en trepada puede obtenerse como:
Z t2 Z h2
Ẇf,climb
Wf,climb = Ẇf,climb dt = dh (24)
t1 h1 R.C.
h1 + h2
h1−2 =
2
0
3. Se propone un valor inicial de consumo de combustible Wf,climb .
5. Se obtiene la tasa de ascenso R.C.i1−2 evaluada en el peso W1−2 y la altura h1−2 . La velocidad a
la cual se evalúa esta relación de ascenso es la correspondiente a la máxima velocidad de trepada
para esa condición de peso y altura.
i+1
7. Se obtiene el nuevo valor de combustible consumido Wf,climb como:
!
i
Ẇf,climb
i+1
Wf,climb = (h2 − h1 )
R.C.i
Una vez obtenido el peso del combustible consumido durante la trepada, se puede calcular el tiem-
po requerido para la misma mediante la relación de ascenso correspondiente, a través de la siguiente
expresión:
h2 − h1
t1−2 = (25)
R.C.1−2
Los valores obtenidos mediante los desarrollos descriptos se incluyen en la Tabla 3.3.3.
El análisis del alcance se realiza considerando un vuelo cuasi-estacionario nivelado con viento en
reposo, por lo que son válidas las ecuaciones 5 y 6 de movimiento.
Se considera además, que la variación en el peso de la aeronave (dW ) se debe únicamente a la
variación de consumo de combustible (dWf ) por lo tanto:
dWf = −dW (27)
Se define el consumo especı́fico de combustible como:
" #
flujo de combustible N
cp = (28)
Pe kW · hr
Este valor se obtiene a partir de la siguiente tabla:
El cálculo del alcance de la aeronave en estudio, para altura y velocidad de vuelo constantes, se realiza
para las tres condiciones de carga que se detallan a continuación:
Carga paga máxima.
Carga máxima de combustible con carga paga.
Carga de combustible máxima sin carga paga.
Los valores obtenidos se muestran en la siguiente tabla, con la que se genera el diagrama de carga
paga de la Figura 9.
La autonomı́a se define como el tiempo transcurrido en el cual una aeronave permanece en vuelo para
un consumo de combustible y una condición de vuelo dada, mientras que la autonomı́a especı́fca (SE) es
definida por la duración del vuelo por unidad de peso de combustible consumida.
Z Wf,2
dt 1000
SE = = 0 −→ E= SE dWf (34)
dWf cp Pe Wf,1
Reemplazando la ec. 7 en la ec. 34, se obtiene la expresión general (todas las condiciones de vuelo)
de la autonomı́a para la aeronave en estudio.
Z W2
1000 CL dW
E [hr] = − (35)
Wf cp V CD W
Bajo la condicion de vuelo crucero a velocidad y altura constantes, la ec.35 toma la siguiente forma:
2ηp Em −1 E1 ζ
Eh−V [hr] = 1000 tan (36)
cp V 2Em [1 − KE1 CL,1 ζ]
La sigiente tabla presenta los valores de alcance obtenidos para las condiciones de carga establecidas
en el inciso anterior.
T −D =0 (38)
WV 2
W − L cos φ = 0 con: L sin φ = (39)
G Rt
1. Se elige una velocidad de vuelo dentro del rango de velocidades de vuelo mı́nima y máxima del vuelo
recto y nivelado en función de la altura h.
2. Se calcula el coeficiente de sustentación CL,L correspondiente a la condición de vuelo recto y nive-
lado para dicha velocidad.
3. Se obtiene el factor de carga máximo limitado por la aerodinámica, como:
CL,max
nmax,CL ,max = (41)
CL,L
6. El factor de carga máximo se obtiene como el mı́nimo entre los factores de carga permitidos por la
aerodinámica, la potencia disponible y la carga estructural.
nmax = min(nmax,CL ,max , nmax,Pav , nstr ) (44)
El radio de giro Rt y la velocidad de giro ψ̇, se obtienen a partir de las ecuaciones de movimiento y
relaciones cinemáticas, como:
V 2W V2 V2 1
Rt = = = √ (45)
gL sin φ g tan φ g n2 − 1
√
V g n2 − 1
ψ̇ = = (46)
Rt V
Se calcularon los factores de carga para cuatro alturas (0m, 3000m, 6000m y 9000m). Para el rango de
velocidades, se tomaron 100 puntos desde la velocidad de pérdida (Vs ) hasta la velocidad máxima (Vmax )
dadas por la envolvente de vuelo (Figura 8). Con los valores del factor de carga, se obtuvieron los radios
y la velocidades de giro correspondientes a cada altura (Figuras 11 y12).
En ambos casos se observa que hay un quiebre en la pendiente de la curva, que corresponde a un
cambio en la limitación del factor de carga. Hacia la izquierda, la limitación es por la aerodinámica mientras
que a la derecha se debe a la potencia disponible.
Se destaca que para el radio de giro, las curvas tienden a infinito para la velocidad de pérdida y la
máxima. De un modo similar, la velocidad de giro tiende a cero en los valores lı́mite.
También se graficaron el radio de giro mı́nimo y la velocidad de giro máxima para diferentes alturas
(Figuras 13 y 14).
3.5.1. Despegue.
Distancia en tierra.
La distancia en tierra y el tiempo en el cual se raliza esta fase se determina mediante la siguientes
expresiones:
Z VR Z VR
V ± VW V
SN GR = dV = g dV (47)
∓VW ag 0 W [T − Dg − µg (W − Lg ) − Wφ ]
Z VR Z VR
dV dV
tN GR = = g (48)
∓VW a g 0 W [T − D g − µ g (W − Lg ) − Wφ ]
Donde:
1 µg
Lg = ρSCL,g V 2 con: CL,g = = cte. (49)
2 2K
1
Dg = ρSCD,g V 2 con: CD,g = CD0 + KCL,g = cte. (50)
2
η p Pe
T = (51)
V
Distancia de Transición.
El coeciente de sustentacion en la transicion (CLT R ) puede ser escrito como:
2W
CLT R = 2 + ∆CL (56)
ρSVLO
Con el coeficiente de sustentación CLT R se puede calcular el coeficiente de resistencia CDT R a través
de la curva polar, y con éste último a su vez la resistencia de transicion DT R :
Distancia de Trepada.
La etapa final del despegue se calcula a través de las siguientes expresiones:
hscreen − hT R
SCL = = 50,24 [m] (64)
tan θCL
2SCL
tCL = = 1,31 [s] (65)
VLO + V2
Donde:
ST R ST R θCL
hT R = RT R − RT R cos θCL = (1 − cos θCL ) ≈ (66)
sin θCL 2
3.5.2. Aterrizaje.
Distancia de transición.
Distancia de rotación.
Distancia en tierra.
En la Figura 16, se puede observar la maniobra de aterrizaje. Durante el descenso desde la altura del
obstáculo hasta la etapa de flare, la velocidad varı́a desde un valor igual a la velocidad de aproximación
VA hasta un valor igual a la velocidad a la cual es realizado el flare VF L . Usualmente se asume que
VF L ≈ 0,95VA .
VF2L
hf lare = Rf lare (1 − cos γ̄A ) con: Rf lare = (73)
g(nF L − 1)
Distancia de Transición.
Durante el flare, la distancia recorrida SLT R y el tiempo empleado tLT R se determinan como:
Done:
VLR ∓VW
V ± VW
Z
SLR = dV =
VT D ∓VW ag
Z VLR
V
= g dV (79)
VT D W [T − Dg − µgbrake (W − Lg ) − Wφ ]
Z ∓VW
V ± VW
SLN GR = dV =
VLR ∓VW ag
Z 0
V
= g dV (80)
VLR W [T − Dg − µgbrake (W − Lg ) + Nn (µgbrake − µg ) − Wφ ]
Para calcular esta distancia es necesario determinar la velocidad de contacto en pista VT D , para ello
puede utilizarse la siguiente expresión basada en el balance de energı́a:
v" #
u
u W W
VT D = t VF2L + hf lare W − (D − T )T D SLT R = 40,45 [m/s] (81)
2g 2g
Donde:
C
1 D
(D − T ) = W γ̄A + (82)
2 CL T D
Al igual que para el cálculo de SN GR , las ecuaciones 79 y 80 se resuelven por integracióon numéerica
utilizando el método del trapecio, obteniendo los siguientes resultados:
4. Conclusión.
Comparando los resultados obtenidos con los datos técnicos de la aeronave (Tabla 2), se observa una
diferencia en los valores, la cual se debe principalmente a la aproximación de la curva polar (CD0 + KCL2 ),
además de las simplificaciones realizadas en las hipótesis. No obstante, dichas simplificaciones permiten
realizar un rápido análisis de las prestaciones de una aeronove que resultan en una buena aproximación,
dando lugar, en el caso de que se trate de un anteproyecto, a realizar modificaciones en el diseño.
Bibliografı́a
[1] J. Roskam and C. E. Lan. Airplane Aerodynamics and Performance, DARcorporation, 1997
[5] [Link]
[6] [Link]