Facultad de Ciencias de la Ingeniería y Tecnología FCITEC

Análisis modal para frecuencias naturales en un

perfil NACA 2415
Cuevas Aceves V.M.

Universidad Autónoma de Baja California

Facultad de Ciencias de la Ingeniería y Tecnología (FCITEC), UABC,
Tijuana, Baja California, México

valeria.cuevas@uabc.edu.mx

Resumen. - En este documento se analizará las frecuencias naturales de vibración de un perfil alar a través de
un análisis modal de manera analítica y numérica. Se realizo el modelo del perfil alar en SolidWorks para
posteriormente ser analizado en Ansysworkbench. Se compararon los resultados obtenidos de manera teórica y
simulada y se obtuvo que son aproximadamente muy cercanos con un error menor al 5%.

Palabras clave: Perfil; modal; frecuencia; vibración.

Abstract. - In this paper we will analyze the natural frequencies of vibration of an airfoil through an analytical
and numerical modal analysis. The airfoil was modeled in SolidWorks and then analyzed in Ansysworkbench. The
results obtained in a theoretical and simulated way were compared and it was obtained that they are approximately
very close with an error of less than 5%.

Key words: Profile; modal; frequency; vibration.

Para poder llevar a cabo este análisis modal,
1. Introducción se seleccionará un perfil alar NACA y de modelara a
través de un software CAD como lo es SolidWorks
Dentro de los distintos análisis que se les para posteriormente ser simulado en ANSYS
realizan a las estructuras, uno de ellos es el análisis Workbench y así poder comparar los resultados
modal el cual se utiliza para conocer la respuesta del obtenidos teóricamente con los obtenidos
sistema ante las vibraciones producidas en él, nos numéricamente.
permite evaluar las propiedades como la frecuencia
natural [1]que es lo que se llevara a cabo en este 2.1 Modelado del perfil alar en SolidWorks
reporte. Se realizará el análisis modal a un perfil alar
NACA 2415. El análisis de las vibraciones libres del El modelado del perfil alar se realizó
ala de un avión es uno de los problemas más difíciles utilizando la función de curvas por puntos, los datos
de abordar en el proceso de diseño de una aeronave. se adquirieron de una base de datos [3] obteniendo
Un aspecto fundamental en el examen de las así la sección transversal del perfil y se procedió a
respuestas a las ráfagas y los fenómenos darle una extrusión de 5 m, lo que nos deja un perfil
aeroelásticos es investigar las frecuencias y modos alar con una cuerda de 1 m y un largo de 5 m.
naturales. Las alas de una aeronave suelen estar
formadas por estructuras de paredes finas, no
uniformes y de forma arbitraria que suelen tener
ángulos de barrido y diedro y son cónicas. [2]

2. Metodología

1
Revista de Ciencias de la Ingeniería y Tecnología (RECIT). Vol. 1. No. 1. 2017

ℎ = 𝑚𝑎𝑥[(𝑍𝑢 + 𝑍𝑖 )⁄2] (9)

Donde:

K = coeficiente de proporcionalidad
c = cuerda del perfil alar
𝜏 = Relación de espesor
𝜀 = Relación de inclinación
t = máximo espesor
h = máxima inclinación
Figura 1. Perfil alar modelado en SolidWorks.
Para obtener los coeficientes de
2.2 Cálculos teóricos proporcionalidad tenemos:
Para poder calcular la frecuencia natural a 1 𝑐
través del análisis modal del perfil alar es necesario 𝐾𝐴 → ∫ [𝑍 − 𝑍𝑖 ] 𝑑𝑥 (10)
primero obtener el área y el momento de inercia del 𝑐2𝜏 0 𝑢
𝑐
perfil: 1 1
𝐾𝐼 → 4 2 2
∫ [(𝑍𝑢 − 𝑧̅)3
𝑐 𝜏(𝜏 + 𝜀 ) 𝑜 3
− (𝑍𝑖 − 𝑧̅)3 ] 𝑑𝑥 (11)

Al evaluar las expresiones anteriores para los

perfiles más comunes obtenemos que los
coeficientes se aproximan a un valor:

𝐾𝐴 ≅ 0.6 (12)
Figura 2. Variables para obtener el momento de inercia y 𝐾𝐼 ≅ 0.036 (13)
el área de la sección transversal de un perfil alar. [4]

𝑐
Utilizando los datos adquiridos sobre el perfil
𝐴 = ∫ [𝑍𝑢 − 𝑍𝑖 ] 𝑑𝑥 (1) 2415 y las ecuaciones previamente vistas ya
0 podemos obtener el área de la sección transversal y
𝑐
1 1 el momento de inercia:
𝑧̅ = ∫ [𝑍𝑢2 − 𝑍𝑖2 ] 𝑑𝑥 (2)
𝐴 0 2
𝑐 Tabla 1. Datos obtenidos del perfil.
1
𝐼=∫ [(𝑍𝑢 − 𝑧̅)3 − (𝑍𝑖 − 𝑧̅)3 ] 𝑑𝑥 (3)
𝑜 3 𝒁𝒖 + 𝒁𝒊
c 𝒁𝒖 𝒁𝒊 𝒁𝒖 − 𝒁𝒊
𝟐
Donde: 1 0.00157 -0.00157 0.00314 0
0.9974 0.00223 -0.00187 0.0041 0.00018
𝑍𝑢 = Superficie superior
0.9893 0.00419 -0.00275 0.00694 0.00072
𝑍𝑖 = Superficie inferior
0.79536 0.04498 -0.02222 0.0672 0.01138
𝑧̅ = Superficie neutral
0.75154 0.05267 -0.02628 0.07895 0.013195
A = Área de sección transversal
I = Momento de inercia 0.7049 0.06025 -0.03048 0.09073 0.014885
0.65595 0.06753 -0.03472 0.10225 0.016405
Para poder aproximarnos de una manera más 0.60524 0.07431 -0.03893 0.11324 0.01769
correcta al área y momento de inercia tenemos las 0.55331 0.0804 -0.04298 0.12338 0.01871
siguientes relaciones: 0.39597 0.09272 -0.05272 0.14544 0.02
0.34448 0.09409 -0.05483 0.14892 0.01963
𝐴 = 𝐾𝐴 𝑐 2 𝜏 (4) 0.29469 0.09366 -0.05633 0.14999 0.018665
𝐼 = 𝐾𝐼 𝑐 4 𝜏(𝜏 2 + 𝜀 2 ) (5) 0.02145 0.0346 -0.02985 0.06445 0.002375
𝜏 = 𝑡⁄𝑐 (6) 0.00878 0.02319 -0.02103 0.04422 0.00108
𝜀 = ℎ⁄𝑐 (7) 0.00161 0.01161 -0.01106 0.02267 0.000275
𝑡 = 𝑚𝑎𝑥(𝑍𝑢 − 𝑍𝑖 ) (8)
2
Revista de Ciencias de la Ingeniería y Tecnología (RECIT). Vol. 1. No. 1. 2017

De la tabla obtenemos los valores máximos de

t y h usando las ecuaciones (8) y (9): 𝐸 = 69𝑥109 𝑁⁄𝑚2
𝜌 = 2700 𝑘𝑔⁄𝑚3
𝑡 = 0.14999 𝑚
ℎ = 0.02 𝑚 2.3 Simulación en ANSYS Workbench

Utilizando (6) y (7) tenemos que: Para realizar la simulación se utilizó el

módulo de Análisis modal que nos proporciona
0.14999 𝑚 Ansys, y se utilizaron las siguientes condiciones:
𝜏 = 𝑡 ⁄𝑐 = = 0.14999
1𝑚
0.02 𝑚 • Mallado de 50 mm.
𝜀 = ℎ ⁄𝑐 = = 0.02
1𝑚 • Soporte fijo en una de las caras del perfil.
Por lo tanto, ya podemos obtener el área y el • Desplazamiento a lo largo del perfil alar
momento de inercia utilizando las ecuaciones (4), o Componente en X: 0 (restringido)
(5), (12) y (13): o Componente en Y: libre
o Componente en Z: libre
𝐴 = 𝐾𝐴 𝑐 2 𝜏 = (0.6)(1 𝑚)2 (0.14999) • Desplazamiento remoto en la otra cara del
𝐴 = 0.089994𝑚2 perfil
𝐼 = 𝐾𝐼 𝑐 4 𝜏(𝜏 2 + 𝜀 2 ) o Componente en X: 0 (restringido)
𝐼 = (0.036)(1𝑚)4 (0.14999)[(0.14999)2 o Componente en Y: libre
+ (0.02)2 ] o Componente en Z: 0 (restringido)
𝐼 = 1.236355𝑥10−4 𝑚4 o Rotación en X: libre
o Rotación en Y: 0°
Ahora que ya tenemos estas variables o Rotación en Z: 0°
podemos proceder a calcular la frecuencia natural del • Se buscarán 10 modos para resolver.
perfil alar tomándolo como una viga que esta • La solución que se le aplicara al modelo es
empotrada en un extremo y libre en el otro: una deformación total.

𝐸𝐼
𝜔𝑛 = (𝛽𝑛 𝑙)2 √
𝜌𝐴𝐿4

Donde:

𝜔𝑛 = Frecuencia natural
𝛽𝑛 𝑙 = Angulo natural
E = Modulo elástico del material
𝜌 = Densidad del material
L = Longitud del perfil Figura 3. Perfil alar con un mallado de 50 mm y un total
de 4600 elementos.
Tabla 2. Valores de 𝛽𝑛 𝑙 para una viga voladiza [5].

n 𝜷𝒏 𝒍
1 1.8751
2 4.69409
3 7.8539
4 10.99557
5 14.1372
6 17.279 Figura 4. Soporte fijo en un extremo del perfil alar.

Como material se utilizará Aluminio 6061

para el cual:

3
Revista de Ciencias de la Ingeniería y Tecnología (RECIT). Vol. 1. No. 1. 2017

Tabla 5. Comparación de resultados.

Frecuencias naturales (Hz)

𝑴𝒐𝒅𝒐𝒔 Analítico Numérico % de error
1 4.1940 4.3847 4.3492
2 26.2841 27.368 3.9604
3 73.5797 76.14 3.3626
Figura 5. Desplazamiento a lo largo del perfil alar. 4 144.2191 147.84 2.4492
5 238.4045 241.57 1.3103
6 356.1433 355.88 0.0739

Figura 6. Desplazamiento remoto en el otro extremo del

perfil alar.

3. Resultados
3.1 Tablas Figura 7. Grafica de resultados.

Tabla 3. Frecuencias de vibración para perfil alar NACA 3.2 Modos en la simulación
2415 calculadas analíticamente.

n 𝜷𝒏 𝒍 𝝎𝒏 (𝒓𝒂𝒅⁄𝒔) 𝒇𝒏 (𝑯𝒛)
1 1.8751 26.3521 4.1940
2 4.69409 165.1484 26.2841
3 7.8539 462.3152 73.5797
4 10.99557 906.1559 144.2191
5 14.1372 1497.9399 238.4045
6 17.279 2237.7144 356.1433

Tabla 4. Frecuencias de vibración para perfil alar NACA Figura 8. Modo 1.

2415 calculadas numéricamente.

n 𝜷𝒏 𝒍 𝒇𝒏 (𝑯𝒛)
1 1.8751 4.3847
2 4.69409 27.368
3 7.8539 76.14
4 10.99557 147.84
5 14.1372 241.57
6 17.279 355.88

Figura 9. Modo 2.

4
Revista de Ciencias de la Ingeniería y Tecnología (RECIT). Vol. 1. No. 1. 2017

4. Conclusiones
Los valores obtenidos de manera analítica y
numéricamente se aproximan correctamente bajo un
porcentaje de error no mayor al 5%, lo que nos quiere
decir que el análisis teórico realizado es correcto. El
análisis modal nos permite conocer como responderá
el perfil alar bajo las condiciones a las que el entorno
donde se encuentren le apliquen. Para poder llevar a
Figura 10. Modo 3. cabo este análisis de un perfil alar fue necesario
analizarlo como una viga en voladizo ya que si
observamos como se encuentran los perfiles en un
avión, estos están empotrados en un extremo y libres
en el otro.

Referencias
[1] Grupo CTT, “¿Qué es y para qué sirve un
análisis modal?,” May 28, 2021.
Figura 11. Modo 4.
[2] M. BAYRAKTAR and A. DEMİRTAŞ,
“FREE VIBRATION ANALYSIS OF AN
AIRCRAFT WING BY CONSIDERING
AS A CANTILEVER BEAM,” Selcuk
University Journal of Engineering ,Science
and Technology, vol. 7, no. 1, pp. 12–21,
Mar. 2019, doi: 10.15317/Scitech.2019.178.
[3] Airfoil, “NACA 2415 - NACA 2415
airfoil.”
[4] “Area and Bending Inertia of Airfoil
Figura 12. Modo 5. Sections.”

Figura 13. Modo 6.