Dept.

de Aeronaves y Vehículos Espaciales - Vibraciones y Aeroelasticidad 6-Jul-2023

UNIVERSIDAD POLITÉCNICA DE MADRID
Número de expediente
Apellidos
Nombre
EJERCICIO 1 (60 minutos)
Un avión de ala con echa nula realiza un viraje manteniendo la altura de vuelo, para lo cual deecta de forma antisimétrica
los alerones (ala izquierda, alerón giro positivo δc hacia abajo; ala derecha, alerón negativo −δc hacia arriba) adquiriendo un
ángulo de balance ϕ de la gura. El ángulo de balance produce una desviación de la sustentación respecto a la vertical y, para
compensar el peso de la aeronave y mantener la altura, se debe aumentar el ángulo de ataque lo que conlleva un aumento del
factor de carga. El incremento de ángulo de ataque se representa con αw y el incremento del factor de carga se representa con
∆Nz , ambos proporcionales al ángulo de balance ϕ con constantes k1 y k2 , respectivamente. El comportamiento aeroelástico
del ala izquierda y derecha se representan con sus respectivas secciones típicas al 75 % de la semienvergadura; la sección típica
del ala izquierda se representa en la gura inferior. En todo momento, se considera que el viraje del avión es muy lento y, por
tanto, se puede considerar un proceso estacionario.
Utilizando como referencia el vuelo inicial equilibrado a una presión dinámica de vuelo q∞ , se pide analizar el viraje posterior
contestando a los siguientes apartados:
1. (1.5 puntos) Formular la ecuación de equilibrio estático en torsión (respecto al eje elástico) que permite determinar la
torsión elástica αeL de la sección típica del ala izquierda, dejándola en función de las variables αw y ∆Nz que aparecen en
el viraje del avión. Recuérdese que, en el ala izquierda, la deexión de alerón es positiva δc (hacia abajo).
2. (1.5 puntos) Repetir el apartado anterior para determinar la torsión elástica del ala derecha αeR . Recuérdese que, en el
ala derecha, la deexión del alerón se debe expresar como −δc ya que es negativa (hacia arriba).
Como se ha mencionado anteriormente, el ángulo de ataque incremental del ala durante el viraje es αw = k1 ϕ y el incremento
del factor de carga viene dado por la expresión ∆Nz = k2 ϕ. Utilizando estas expresiones, responder a los siguientes apartados:
1. (1 punto) Adimensionalizar la ecuación de equilibrio en torsión del ala izquierda y obtener la torsión elástica del ala
izquierda αeL y el ángulo αeL + k1 ϕ. Deben aparecer, entre otros, la presión dinámica adimensionalizada q̂∞ = q∞ /qD ,
mgd
donde qD es la presión dinámica de divergencia del perl de la sección típica, y el parámetro µ = .
qD SeCLα
2. (1 punto) Repetir el apartado anterior para el ala derecha.
∆LL
3. (1 punto) Calcular la efectividad de mando en la sección típica del ala izquierda como ηL = .
q∞ SCLδ δc
4. (1 punto) Repetir el apartado anterior para el ala derecha, calculando la efectividad ηR .
5. (1 punto) Calcular la presión dinámica de inversión de mando en la sección típica del ala izquierda q̂R
L
, adimensionalizada
con qD .
6. (1 punto) Repetir al apartado anterior para el ala derecha, calculando la presión dinámica de inversión de mando q̂R
R
.
7. (1 punto) ¾Qué ala, izquierda o derecha, presenta la menor presión dinámica de inversión de mando? ¾Por qué?.
𝑳
𝜙

𝑐/2 x
h 𝐾ℎ
𝐾𝛼
AC
𝛼𝑒 EA

𝛼𝑤
𝐾𝛿
𝛿𝑐
𝑒 𝑑

U

ℎ, 𝛼, 𝛿𝑐 flexión, torsión elástica y rotación de la superficie de control, respectivamente

𝑑, 𝑒 distancia del eje elástico al centro de gravedad y al centro aerodinámico, respectivamente
𝑚𝑔 peso por unidad de envergadura
𝐾ℎ , 𝐾𝛼 , 𝐾𝛿 rigideces elásticas de flexión, torsión y rotación de superficie de control, respectivamente
𝐶𝐿𝛼 , 𝐶𝑀𝐴𝐶 pendiente del coeficiente de sustentación y momento del perfil
𝐶𝐿𝛿 , 𝐶𝑀𝐴𝐶𝛿 pendiente del coeficiente de sustentación y momento del perfil asociados a la rotación de superficie de control
SOLUCIÓN:

Ala izquierda: 𝑞 𝑆𝐶 𝛿

𝑞 𝑆𝐶 𝛼 𝛼 𝑞 𝑆𝑐𝐶 𝛿 𝐾

𝛼 EA
AC

𝛼 𝑘 𝜙
𝛿
𝑒 𝑑
k W𝜙
U

Ala derecha:
𝑞 𝑆𝐶 𝛼 𝛼
𝑞 𝑆𝑐𝐶 𝛿 𝐾
AC
𝛼 EA
𝑞 𝑆𝐶 𝛿 𝛿

𝛼 𝑘 𝜙

𝑒 𝑑
k 𝑊𝜙
U

Apartado 1-1 (1.5 puntos)

Kα αeL = q∞ SeCLα αw + αeL + q∞ SeCLδ δc + q∞ ScCM ACδ δc + mgd · ∆Nz

Apartado 1-2 (1.5 puntos)

Kα αeR = q∞ SeCLα αw + αeR − q∞ SeCLδ δc − q∞ ScCM ACδ δc + mgd · ∆Nz

Apartado 2-1 (1 punto)

Kα αeL = q∞ SeCLα k1 ϕ + αeL + q∞ SeCLδ δc + q∞ ScCM ACδ δc + mgd · k2 ϕ

 
c CM ACδ
(Kα − q∞ SeCLα ) αeL
= q∞ SeCLδ 1 + δc + q∞ SeCLα k1 ϕ + mgd · k2 ϕ
e CLδ
 
CLδ c CM ACδ mgd
(qD − q∞ ) αeL = q∞ 1+ δc + q∞ k1 ϕ + k2 ϕ
CLα e CLδ SeCLα
 
CLδ c CM ACδ
q̂∞ 1+ δc + q̂∞ k1 ϕ + µk2 ϕ
L CLα e CLδ mgd
αe = , donde µ =
1 − q̂∞ qD SeCLα

   
CLδ c CM ACδ CLδ c CM ACδ
q̂∞ 1+ δc + q̂∞ k1 ϕ + µk2 ϕ q̂∞ 1+ δc + (k1 + µk2 ) ϕ
CLα e CLδ CLα e CLδ
k1 ϕ + αeL = k1 ϕ + =
1 − q̂∞ 1 − q̂∞

Apartado 2-2 (1 punto)

 
CLδ c CM ACδ mgd
(qD − q∞ ) αeR
= −q∞ 1+ δc + q∞ k1 ϕ + k2 ϕ
CLα e CLδ SeCLα
 
CLδ c CM ACδ
−q̂∞ 1+ δc + q̂∞ k1 ϕ + µk2 ϕ
CLα e CLδ
αeR =
1 − q̂∞
   
CLδ c CM ACδ CLδ c CM ACδ
−q̂∞ 1+ δc + q̂∞ k1 ϕ + µk2 ϕ −q̂∞ 1+ δc + (k1 + µk2 ) ϕ
CLα e CLδ CLα e CLδ
k1 ϕ + αeR = k1 ϕ + =
1 − q̂∞ 1 − q̂∞
Apartado 2-3 (1 punto)

∆LL = q∞ SCLα k1 ϕ + αeL + q∞ SCLδ δc

 
c CM ACδ
q̂∞ 1+
∆LL CLα k1 ϕ + αeL e CLδ CLα k1 + µk2 ϕ
ηL = =1+ =1+ + =
q∞ SCLδ δc CLδ δc 1 − q̂∞ CLδ 1 − q̂∞ δc
c CM ACδ
1 + q̂∞
e CLδ CLα k1 + µk2 ϕ
= +
1 − q̂∞ CLδ 1 − q̂∞ δc

Apartado 2-4 (1 punto)

∆LR = q∞ SCLα k1 ϕ + αeR − q∞ SCLδ δc

 
c CM ACδ
−q̂∞ 1+
∆LR CLα e CLδ CLα k1 + µk2 ϕ
k1 ϕ + αeR = −1 +


ηR = = −1 + + =
q∞ SCLδ δc CLδ 1 − q̂∞ CLδ 1 − q̂∞ δc
c CM ACδ
−1 − q̂∞
e CLδ CLα k1 + µk2 ϕ
= +
1 − q̂∞ CLδ 1 − q̂∞ δc

Apartado 2-5 (1 punto)

L c CM ACδ
1 + q̂R
e CLδ CLα k1 + µk2 ϕ
L
+ L δ
=0
1 − q̂R CLδ 1 − q̂R c
CLα ϕ
−1 − (k1 + µk2 )
LcCM ACδ CLα ϕ L CLδ δc
1+ q̂R + (k1 + µk2 ) = 0 ⇒ q̂R =
e CLδ CLδ δc c CM ACδ
e CLδ

Apartado 2-6 (1 punto)

c CM ACδ
−1 − q̂∞
e CLδ CLα k1 + µk2 ϕ
+ =0
1 − q̂∞ CLδ 1 − q̂∞ δc
CLα ϕ
−1 + (k1 + µk2 )
c CM ACδ CLα ϕ R CLδ δc
−1 − q̂∞ + (k1 + µk2 ) = 0 ⇒ q̂R =
e CLδ CLδ δc c CM ACδ
e CLδ

Apartado 2-7 (1 punto)

La inversión de mando ocurre antes en el ala derecha debido a que el incremento de ángulo de ataque αw y el incremento
de factor de carga ∆Nz asociados al viraje empeoran su efectividad de mando.
 Dept. de Aeronaves y Vehı́culos Espaciales - Vibraciones y Aeroelasticidad 6-Jul-2023
UNIVERSIDAD POLITÉCNICA DE MADRID

Número de expediente
Apellidos
Nombre

EJERCICIO 2 (60 minutos)

Enunciado: Se desea analizar un diseño mitigador del flameo que consiste en la instalación de una
masa puntual en el borde de ataque. Para ello se estudia la sección tı́pica del ala caracterizada por
una cuerda c y con un espesor t y densidad volumétrica ρ. En el borde de ataque se añade una
masa puntual de valor βρct. El eje elástico se encuentra en la posición mostrada en la figura, las
rigideces a flexión y torsión de la sección tı́pica son Kh y Kα respectivamente y se asume el mismo
amortiguamiento estructural en ambos grados de libertad g = gh = gα .
Se pide:
1. (2 puntos) Determinar la expresión de la energı́a cinética del sistema en función de los parámetros
del sistema (ρ, c, t, β, Λ1 ) y definir a partir de ella la masa del sistema M , el momento de inercia
Iα y el momento estático Sα .
2. (1 punto) Escribir la expresión de la función de disipación y la expersión de la energı́a potencial
en función de las frecuencias naturales en flexión, ωh , y torsión ωα , M e Iα .
3. (1 punto) Las ecuaciones de la dinámica del sistema indicando las fuerzas generalizadas sobre
los grados de libertad h(t) y α(t), Qh y Qα respectivamente.
4. (1.5 puntos) Asumir movimiento armónico de frecuencia ω y adimensionalizar las ecuaciones
anteriores para que aparezcan los parámetros habituales (µ, xα , rα , k).
5. (2 puntos) Formular las fuerzas generalizadas para aerodinámica cuasi-estacionarias desarrolla-
das por Fung y expresarlas en el dominio de la frecuencia en forma adimensional acorde al punto
anterior para que aparezcan los parámetros habituales (µ,k).
Para el caso β = 1, Λ1 = −1/2:
6. (0.5 puntos) Expresar las ecuaciones en formato matricial de la forma
h ii  
h̃/(c/2)
hh i
h i
ωα 2 1
M̂ij − (1 + ig) ω K̂ij + 2πµk2 Q̂ij = ...
α̃

7. (2 puntos) Analizar la existencia de flameo para este caso.

SOLUCIÓN:

Apartado 1
Z c/2
1 h  c i2 1 h  c c i2 1 1
T = (ρ1tdx) ḣ + α̇ x − Λ1 + βρct ḣ + α̇ − − Λ1 = M ḣ2 + Iα α̇2 + Sα ḣα̇
−c/2 2 2 2 2 2 2 2
Z c/2
M= ρ1tdx + βρct = (1 + β) ρct = (1 + β) m
−c/2
c/2
Λ21 β
Z  
c 2  c c 2 2 1 2
Iα = ρ1tdx x − Λ1 + βm − − Λ1 = mc + + (1 + Λ1 )
−c/2 2 2 2 12 4 4
Z c/2  c  c c c
Sα = ρ1tdx x − Λ1 + βm − − Λ1 = −m [Λ1 + β (1 + Λ1 )]
−c/2 2 2 2 2
Apartado 2
1 gKh 2 1 gKα 2 1 gM ωh2 2 1 gIα ωα2 2
D= ḣ + α̇ = ḣ + α̇
2 ω 2 ω 2 ω 2 ω
1 1 1 1
U = Kh h2 + Kα α2 = M ωh2 h2 + Iα ωα2 α2
2 2 2 2
Apartado 3
          
M Sα ḧ g M ωh2 0 ḣ M ωh2 0 h Qh
+ 2 + 2 =
Sα Iα α̈ ω 0 Iα ωα α̇ 0 Iα ωα α Qα
Apartado 4
 ( h̃
)
 ω 2  ( h̃
) (
Qh
)
1 xα α ωh2 /ωα2 0 M (c/2)ω 2
− c/2 + (1 + ig) c/2 =
xα rα2 α̃ ω 0 rα2 α̃
Qα
M (c/2)2 ω 2

Apartado 5
( ) ( )
Qh   1

 h̃
M (c/2)ω 2 1 −4πik
−4π
1 + 2
− Λ1
ik  c/2
Qα =
M (c/2)2 ω 2 2πµk 2 4π 21 + Λ1 ik 4π 1
2
+ Λ1 1 + 1
2
− Λ1 ik − iπk α̃
Apartado 6
Para este caso xα = 0, rα2 = 5/12
"  " 2 #  # ( )  
ω h̃
1 0 2 h
0 −4πik −4π (1 + ik) 0
− (1 + ig) ωωα 1


ωα + c/2 =
0 rα2 0 rα2
2πµk2 0 −iπk α̃ 0

Apartado 7
Para que exista flameo la matriz del sistema de ecuaciones del apartado anterior debe ser singular,
tener determinante nulo:
" #
 ω 2  ω 2 2i  ω 2 i

α h 2 2 α
1 − (1 + ig) − rα − rα (1 + ig) − =A×B =0
ω ωα µk ω 2µk
Para que sea nulo han de serlo A o B, tanto su parte real como su parte imaginaria.
 ω 2  ω 2
α h
Re(A) = 1 − = 0 → ω 2 = ωh2
ω ωα
 ω 2  ω  2 2 2
α h
Im(A) = −g − = −g − = 0 → No tiene solución si g > 0, k > 0
ω ωα µk µk
 ω 2
α
Re(B) = rα2 − rα2 = 0 → ω 2 = ωα2
ω
 ω 2 1 1
2 α
Im(B) = −rα g − = −rα2 g − = 0 → No tiene solución si g > 0, k > 0
ω 2µk 2µk
Por lo tanto no puede existir flameo para el sistema analizado (β = 1, Λ1 = −1/2).