Análisis de vibraciones forzadas para una viga

Sebastián Castrillón & Cristian Johan Rodríguez Salazar

sebastian012396@hotmail.com & cjohanrs@gmail.com
Facultad Tecnológica- Universidad Distrital Francisco José de Caldas

Resumen - En el siguiente documento se presenta el análisis Donde:

de las vibraciones forzadas en una viga con perfil I sometida
a una carga puntual de 250 kN y con apoyos simples en los fn = Frecuencia natural de la viga [Hz]
extremos, se determina los comportamientos de la viga de
Cn = Número característico del modo de vibración
manera analítica, conceptual y con ayuda de ANSYS. Por
último se evalúan los resultados y se halla el margen de error E = Módulo de Young [Pa]
entre los diferentes métodos. I = Momento de inercia de la sección transversal [m^4]
L = Longitud de la viga, elemento inercial [m]
Palabras clave - Error relativo, análisis, software, sistema, m = Masa del elemento inercial [kg]
viga, frecuencia, vibración, variables, geometría, material.
B. Frecuencia natural en Radianes por segundo.
Abstract -The following document presents the analysis of
the forced vibrations in a beam with profile I subjected to a Recorrido en radianes que realiza el movimiento en un
point load of 250kN and with simple supports at the ends, the segundo [Rad/s]
behavior of the beam is determined analytically, conceptual
and with the help of ANSYS . Finally, the results are Se define como:
evaluated and the margin of error between the different Wn=(2 π *fn) (2)
methods is found.
Donde:
Keywords - Relative error, analysis, software, system,
fn = Frecuencia natural de la viga [Hz]
beam, frequency, vibration, variables, geometry, material.
Wn = Frecuencia natural de la viga [Rad/s]
1. INTRODUCCIÓN
C. Viga con extremos simplemente apoyados
El análisis de vibraciones forzadas es importante ya que
nos permite determinar causas, efectos y con ello las
condiciones vibratorias existentes en cualquier elemento
que está sometido a una carga P(t). Es indispensable tener
conocimiento del software ANSYS ya que cuenta
con herramientas muy efectivas e interesantes al calcular
frecuencias en una viga. El propósito del artículo es
determinar el margen de error entre los diferentes métodos Fig. 2.1. Modo de viga extremos simplemente apoyados con su
respectivo Cn
que se emplean en el cálculo de frecuencias forzadas en
una viga como también invitarlo a hacer uso de la
herramienta informática ANSYS .

2. MARCO TEÓRICO

A. Frecuencia natural (fn)

Frecuencias que pueden existir durante las vibraciones

libres, además son las frecuencias definidas a las que el
cuerpo vibra después de la perturbación inicial.
La frecuencia natural de una viga está dada por la siguiente
expresión.

𝐸∗𝐼
𝑓𝑛 = 𝐶𝑛√ (1)
𝑚∗𝐿3 Fig. 2.2. Cotas de la viga W250x80
 Por último se diferenciará cada uno de los datos obtenidos
por cada uno de los métodos en una tabla correspondiente
a los temas indicados, y seguidamente se realizará el
cálculo del error relativo entre estos, con el fin de
evidenciar la magnitud en que varía cada uno de los datos
Fig. 2.3. Dimensiones de la viga W250x80 y poder proponer o sugerir los posibles factores que
ocasionan dicha variación.
D. Vibración forzada por una fuerza constante.
5 . PROCEDIMIENTO
Se estudia la excitación de un sistema de grado de libertad
sometido a un tipo de carga escalón. La ecuación que De forma analítica y modelando la viga en ANSYS
representa la solución completa para el sistema es: WORKBENCH, determine el valor de la
primera frecuencia naturales en [rad/s] y en [Hz] para la
𝜁𝑃𝑜 viga W250x80.
𝑋(𝑡) = 𝑒 −𝜁𝑊𝑛𝑡 [− 𝑆𝑒𝑛 𝑊𝑎𝑡 −
𝑚𝑊𝑛𝑊𝑎
𝑃𝑜 𝑃𝑜
𝐶𝑜𝑠 𝑊𝑎𝑡] + (3) Evalúe solo desplazamientos en el plano XY.
𝑚𝑊𝑛2 𝑚𝑊𝑛2

Dónde: Considere:
ζ= Factor de amortiguamiento - Densidad = 7855 [kg/m^3]
Po=Fuerza externa [N] - Módulo de Young (E) = 200 x 10^9 [Pa]
Wn= Frecuencia [Rad/seg] - Poisson Ratio = 0.33
Wa= Frecuencia amortiguada [Rad/seg] - Longitud = 2270 [mm] = 2.27 [m]
m= Masa de la viga [kg]
5.1. Método analítico.
E. Factor de amortiguamiento (ζ)
A continuación se presenta el desarrollo analítico para el
Para muchas aplicaciones, excepto cuando el cálculo de la primera frecuencia natural.
amortiguamiento es introducido de forma deliberada en
el sistema, el factor de amortiguamiento es menor de 5.1.1. Momento de inercia y área total.
0.1.
Y
3. HERRAMIENTAS

 Software MATLAB R2015b.

 Software ANSYS WORKBENCH
 Guía, vibración forzada- Ing. Wilmer
Cruz Guayacundo.

4. METODOLOGÍA

Se realizará un método analítico conceptual con el uso de

Fig. 5.1. División de la viga con origen en cero (0)
las ecuaciones algebraicas consultadas en la bibliografía
suministrada por el docente para la obtención de la primera El centroide de la viga teniendo en cuenta que es
frecuencia natural en unidades de Hertz [Hz] y radianes totalmente simétrica es 128 mm
por segundo [rad/s]. Posteriormente se analizará la viga
con la ayuda del software ANSYS, según las instrucciones Y(mm) A(mm2) A*Y(mm3) d(mm) Ixx (mm4)
observadas en la sesión de clase, con el fin de obtener 1 7.8 3978 31028.4 120.2 80673.84
2 128 2113.12 270479.36 0 8898883.644
nuevamente la frecuencia natural, antes mencionada.
3 248.2 3978 987339.6 120.2 80673.84
Tabla 5.1. Calculo de los momentos de inercia de cada división.
Después de hallar la primera frecuencia a partir de los dos
(2) métodos descritos anteriormente, con ayuda de Veamos que:
ANSYS y MATLAB se determinará las respuestas de la
viga para las razones de amortiguamiento ζ= 0,01 ζ= 0,05 A=A1+A2+A3= 10069.12mm2 =0.010m2
y ζ= 0,20 (Amplitud) a partir de gráficas. Ixx=Ixx’1+Ixx’2+Ixx’3= 124008837.6mm4=1.24x10-4m4
Nota: Tenga en cuenta que se debe tener conocimientos
previos al teorema de los ejes paralelos, con el fin de
calcular el momento de inercia total de la viga.

5.1.2. Volumen y masa de la viga.

V=A*L=0.010m2*2.27m= 0.0227m3

m= D*V= 7855 kg/m3*0.0227m3 = 178.31 kg

5.1.3. Frecuencia natural en (Hz)

Fig. 5.3. Puntos para determinar número de secciones y longitud
total de la viga.
Según la ecuación (1) la frecuencia natural se determina
de la siguiente manera:
b) Asignar el perfil con las medidas indicadas
anteriormente.

200𝑥109 ∗ 1.24𝑥10−4
𝑓𝑛 = 1.57√ = 171.20 𝐻𝑧
178.31 ∗ 2.273

Recuerde que Cn= 1.57 según figura 1.

5.1.3. Frecuencia natural en (rad/s)

Según la ecuación (2)

Wn =(2π*171.20 Hz) = 1075.67 rad/s

5.2. Modelamiento a partir de ANSYS.

Fig. 5.4. Sección transversal de la viga.

A continuación, se describe un paso a paso para poder

calcular las frecuencias naturales de la viga con ayuda de
ANSYS WORKBENCH:

MODAL

5.2.1. Enginnering Data:

a) Determinar propiedades físicas de la viga, Fig. 5.5. Tabla de dimensiones de la sección transversal de la
viga en milímetros.
densidad, módulo de Young y poisson ratio. (ver
figura 5.2)
c) Una vez generado el perfil se determina y
establece de forma predeterminada y uniforme
por la extensión total de la viga.

Fig. 5.2 Propiedades del material

5.2.2. Geometry:

a) Definir el número de secciones, en que se dividirá

Fig. 5.6. Extensión total de la viga
la viga y mediante puntos, extender líneas hasta
completar la longitud total de la viga establecida.
 5.2.3. Model: 5.2.4. Resultados ANSYS WORKBENCH:

a) Se define el material

Fig. 5.7. Extensión de la malla.

Fig. 5.11. Frecuencia naturales en [Hz] a partir de ANSYS.
b) Se crea una malla y se dimensiona para que se
extienda a lo largo de la viga. Nota: Se selecciona la frecuencia que más se aproxima
al análisis Modal.

5.2.5. Frecuencia natural en [Rad/s]:

Utilizando la ecuación (2) se determina la frecuencia

en [Rad/s]:

Wn=(2 π *179.71Hz)=1129,15 Rad/s

Fig.5.8. Extensión de la malla. TRANSIENT STRUCTURAL

c) Se apoyan dos extremos de la viga 5.2.6. Transient Structural

Primero se realiza un paso a paso similar al MODAL al

definir propiedades y la geometría de la viga.

a) Se ajusta el tiempo de análisis en la viga

Fig. 5.9. Fijación de los extremos de la viga.

.
d) Se permite y concede el desplazamiento
únicamente la dirección vertical (y), según el eje
coordenado seleccionado.

Fig. 5.12. Tiempos de análisis en la viga

b) Se aplica una carga de 250 kN en la mitad de la

viga y dos apoyos simples a cada extremo de la
viga.

Fig.
5.10. Desplazamiento en dirección vertical

e) Por último, se resuelve y soluciona.

Fig. 5.13. Carga de 250kN en una viga con apoyos simples en los
extremos.
 c) Se crea un desplazamiento dejando únicamente
la componente ‘Y’ libre.

Fig. 5.18. Gráfica para ζ=0,20

Fig. 5.14. Viga con desplazamiento libre en el eje Y

d) Se genera la deformación direccional orientado 5.3. Graficas a partir de MATLAB.

en el componente Y, y se resuelve el análisis
transitorio. A continuación se presenta el desarrollo analítico para el
modelamiento de tres (3) gráficas que describen el
comportamiento de tres (3) factores de amortiguamiento,
ζ=0,01, ζ=0,05 y ζ=0,2.

t=0:0.001:0.16 %Intervalo de tiempo (s)

F=250000; %Fuerza en N
E=200e9; %Modulo de elasticidad en Pa
A=0.010; %Área en m^2
L=2.27; %Longitud en m
Fig. 5.15. Deformación direccional orientado en el eje Y
V=A*L; %Volumen de la viga en m^3
D=7855; %Densidad en kg/m^3
5.2.6. Análisis y resultados en Transient Structral
m=D*V; %Masa en kg
Ix=1.24e-4; %Momento de inercia in^4
Cn=1.57; %Cte de modo 1
Cn2=6,21; %Cte de modo 2
Cn3=14.14; %Cte de modo 3
fn=Cn*sqrt((E*Ix)/(m*L^3)); %primera frecuencia
natural en hertz
fn2=Cn2*sqrt((E*Ix)/(m*L^3)); %segunda frecuencia
natural en hertz
fn3=Cn3*sqrt((E*Ix)/(m*L^3)); %tercera frecuencia
natural en hertz
wn=fn*2*pi(); %primera frecuencia natural en rad/s
wn2=fn2*2*pi(); %segunda frecuencia natural en rad/s
Fig. 5.16. Gráfica para ζ=0,01 wn3=fn3*2*pi(); %tercera frecuencia natural en rad/s
c1=0.01; %Cte amortiguamiento 1
c2=0.05; %Cte amortiguamiento 2
c3=0.20; %Cte amortiguamiento 3
wa1=wn*sqrt(1-c1^2);%frecuencia amortiguada en rad/s
wa2=wn2*sqrt(1-c2^2);%frecuencia amortiguada en rad/s
wa3=wn3*sqrt(1-c3^2);%frecuencia amortiguada en rad/s
Ymax=(F*(L^3))/(48*E*Ix); %deflexión máxima en m
Po=Ymax*m*fn^2;
Po2=Ymax*m*fn2^2;
Po3=Ymax*m*fn3^2;
hold on;
A=(-(c1*Po)/(m*wa1*wn));
B=(-Po)/(m*(wn^2));
Fig. 5.17. Gráfica para ζ=0,05
x=exp(-
c1*wn*t).*((A*sin(wa1*t))+(B*(cos(wa1*t))))+((Po/(m* El error relativo en porcentaje, está definido por la
(wn^2)))); %ecuación del sistema 1 siguiente expresión:
plot(t,x,'r')
hold on; 𝐵−𝐴
A=(-(c2*Po2)/(m*wa2*wn2)); %𝐸𝑅𝑅𝑂𝑅 = | | ∗ 100 (4)
𝐴
B=(-Po2)/(m*(wn2^2));
x=exp(- Donde:
c2*wn2*t).*((A*sin(wa2*t))+(B*(cos(wa2*t))))+((Po2/(
m*(wn2^2)))); %ecuación del sistema 2 %ERROR= Error relativo en porcentaje [ % ]
plot(t,x,'g') B= Frecuencia natural con modelamiento en ansys [Hz] ó
hold on; [Rad/s]
A=(-(c3*Po3)/(m*wa3*wn)); A=Frecuencia natural de manera analítica en [Hz] ó
B=(-Po3)/(m*(wn3^2)); [Rad/s]
x=exp(-
c3*wn3*t).*((A*sin(wa3*t))+(B*(cos(wa3*t))))+((Po3/( Utilizando todas las frecuencias en [Hz] tenemos:
m*(wn3^2)))); %ecuación del sistema 3
plot(t,x,'b')
grid on; Modo %ERROR
title('COMPORTAMIENTO VIBRATORIO DE LA
VIGA') 1 4.17
xlabel('tiempo (s)')
Tabla 6.2.. Tabla del error relativo .
ylabel('Desplazamiento (m)')
legend('comportamiento con coeficiente de
0.01','comportamiento con coeficiente de
7. CONCLUSIONES.
0.05','comportamiento con coeficiente de 0.20')
Para el análisis modal, el porcentaje de error relativo entre
frecuencias naturales da 4,17%, este valor asegura un buen
diseño estructural cuando la viga se ve sometida a una
carga de 250kN ya que el margen de error no excede el
8%. Para la primera frecuencia natural, el valor hallado
analíticamente corresponde o se aproxima al modo seis (6)
de la frecuencia natural hallada en ANSYS como se puede
observar en la figura 5.11, que es el valor que se tomó para
realizar el análisis.

Las gráficas generadas a partir de ANSYS WORKBECH

y MATLAB dan como resultado un desplazamiento
Fig. 5.18. Gráfica para los tres factores de amortiguamiento. máximo aproximado de 6,15 mm y 1,2 mm
respectivamente para ζ=0,01, la diferencia no es mucha si
6. ANÁLISIS DE RESULTADOS tenemos en cuenta que la viga tiene dimensiones grandes
en comparación con los resultados (lo mismo pasa con
Una vez se hayan calculado las frecuencias naturales de la ζ=0,05 y ζ=0,20). Para un análisis más detallado y preciso
viga de manera analítica y con ayuda de ANSYS donde un error micro pueda afectar nuestro diseño, es
WORKBECH, se organizan los datos en la siguiente tabla: indispensable contar con una mayor precisión.

Analítico (A) ANSYS (B) 8. BIBLIOGRAFIA

Modo Cn Fn Wn Fn [1] Cruz. W. (s.f) Vibraciones Libres [Documento].

Wn [rad/seg]
[Hz] [rad/seg] [Hz] Bogotá D.C., Colombia: Universidad Distrital
Francisco José de Caldas.
1 1.57 171.68 1075.67 179.15 1129.15
Tabla 6.1. Frecuencias naturales de manera analítica y con ayuda de [2] Cruz. W. (s.f) Vibraciones forzadas [Documento].
ANSYS.
Bogotá D.C., Colombia: Universidad Distrital
Francisco José de Caldas.
Teniendo los valores de las frecuencias naturales se halla
el error relativo entre los resultados.