UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
Nombre y Apellido: Asdrubal Bravo C.I.:14.518.790
ASIGNACIÓ'N 2
Para realizar esta asignación se dispone de dos videos que se anexan y los
cuales presentan la información básica y el desempeño del sistema vibrante.
La asignación requiere:
La elaboración de un código (o rutina) en el lenguaje de su preferencia, para obtener
las respuestas a las preguntas abajo indicadas.
La entrega de los resultados debe llevar las gráficas solicitadas, el código generado
y que éste se encuentre completamente documentado. Una sección de análisis de
los datos y conclusiones también son requeridos.
La fecha de esta asignación es para el 02-03-21 y deberá ser entregada vía internet
el 08-03-21 antes de las 05:00 pm
Parte I: Sistema libre
1a. Observe ambos videos y determine la constante elástica de la barra
Del equilibrio estático las fuerzas que actúan sobre la masa deben ser cero de
modo, en este caso partidos del hecho de.
𝐤∆𝐱 − 𝐦𝐠 = 𝟎
A partir de los datos obtenidos del video 1 se tiene
𝐦𝟏 𝐠 𝐦𝟏 𝐠 𝟎, 𝟓𝟏𝟏 ∙ 𝟗, 𝟖 𝐍 𝐍
𝒌𝟏 = = = −𝟐
= 𝟏𝟒𝟑, 𝟎𝟖
𝚫𝐱 𝟏 𝐋𝟏 − 𝐋𝟎 (𝟏𝟑, 𝟓 − 𝟏𝟎) ∙ 𝟏𝟎 𝐦 𝐦
Para el video 2
𝐦𝟐 𝐠 𝐦𝟐 𝐠 𝟎, 𝟖𝟑𝟔 ∙ 𝟗, 𝟖 𝐍 𝐍
𝒌𝟐 = = = −𝟐
= 𝟏𝟒𝟑, 𝟕𝟑
𝚫𝐱 𝟐 𝐋𝟐 − 𝐋𝟎 (𝟏𝟓, 𝟕 − 𝟏𝟎) ∙ 𝟏𝟎 𝐦 𝐦
El promedio de ambas nos daría la constante del problema experimentalmente
𝐍
𝐤 𝟏 + 𝐤 𝟐 (𝟏𝟒𝟑, 𝟎𝟖 + 𝟏𝟒𝟑, 𝟑𝟔) 𝐦 𝐍
𝐤̅ = = = 𝟏𝟒𝟑, 𝟒𝟏
𝟐 𝟐 𝐦
Para el cálculo de este valor se hizo empleando el programa Octave y para ello se
generó la siguiente rutina.
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
clc; %borra la pantalla
clear all; %borra la memoria
%Ingresamos los datos experimentales o variables
m=input('Ingrese el valor de la masa: '); %ingresamos el valor de la masa para el caso considerado
fprintf('\n');
l0=input('Ingrese el valor de la longitud inicial: '); %ingresamos el valor de la longitud inicial (resorte
sin estirar)
fprintf('\n');
lf=input('Ingrese el valor de la longitud final: '); %ingresamos el valor de la longitud final (resorte
estirado con la masa)
fprintf('\n');
g=9.81;
fprintf('\n');
%Calculamos las relaciones
fprintf('a) El valor de la constante del resorte es:')
d=lf-l0;
k=(m*g)/d %constante del resorte del resorte
1b. ¿Puede aplicar una estrategia analítica para comprobar numéricamente el valor
obtenido en (1a) con los datos experimentales? En caso afirmativo, muéstrelo
gráficamente y explique el procedimiento empleado.
Un procedimiento oportuno para determinar la validez del resultado anterior es
comparar el valor obtenido de la frecuencia de oscilación para los casos (m1=0,511 g) y
(m2=0,836g) con el correspondiente con los datos mostrados en el video 2. En este
sentido, empleamos el valor promedio obtenido de k y a partir de la masa m2 se calcula la
frecuencia natural teórica y después lo comparamos con la frecuencia vista en el video
donde se muestran los datos de t = 4,988 seg necesarios para realizar los 10 ciclos. A
partir de este desarrollo sencillo, se puede comprobar si son iguales o no ambos
resultados.
El valor teórico de la frecuencia natural es el siguiente:
𝐍
𝛚𝟏 𝟏 𝐤̅ 𝟏 √𝟏𝟒𝟑, 𝟒𝟏 𝐦 𝟏𝟑, 𝟎𝟗𝟕𝒔−𝟏
𝐟𝟏 = = √ = = = 𝟐, 𝟎𝟖𝟒𝟓 𝒔−𝟏
𝟐𝛑 𝟐𝛑 𝐦𝟐 𝟐𝛑 𝟎, 𝟖𝟑𝟔𝐤𝐠 𝟐𝝅
Luego la frecuencia natural experimental se obtiene de la forma siguiente
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
𝐍 𝐜𝐢𝐜𝐥𝐨𝐬 𝟏𝟎
𝐟𝟏 = = = 𝟐, 𝟎𝟎𝟓 𝒔−𝟏
𝐓 𝟒, 𝟗𝟖𝟖𝐬
Se puede comprobar que ambos valores difieren entre si en un 4%, para fines
prácticos, se puede entonces comprobar que los valores experimentales y teóricos son
bastante buenos.
Para la comprobación de esta parte se desarrolló el siguiente código programa en
Octave
clc;
clear all;
% Ingresamos los datos
k=input('Ingrese el valor de la constante del resorte: '); %Ingresamos el valor de la constante del resorte
fprintf('\n');
m=input('Ingrese el valor de la masa: '); %Ingresamos el valor de la masa
fprintf('\n');
n=input('Ingrese el numero de Ciclos: '); %Ingresamos el número de ciclos u oscilaciones
fprintf('\n');
t=input('Ingrese el tiempo medido de las oscilaciones: ');
fprintf('\n');
g=9.81;
%Calculamos las relaciones
fprintf('a) El valor de la frecuencia natural teorica es:')
w0=sqrt(k/m);
f=w0/(2*pi)
T=1/f;
fprintf('\n');
fprintf('b) El valor de la frecuencia natural experimental es:')
f=n/t
Parte II: Sistema amortiguado libre
El sistema dado en el video 1 se somete a un proceso de fricción viscosa y el
tiempo registrado en segundos cuando la masa pasa por su punto de equilibrio se
muestra en la siguiente tabla:
0.39 6.47 12.53 18.87 25.13
1.28 7.58 13.59 19.92 26.16
2.26 8.55 14.60 21.00 27.26
3.17 9.52 15.63 21.95 28.28
4.21 10.53 16.78 23.17 29.26
5.36 11.59 17.77 24.02 30.24
Note que los datos van secuencialmente en orden creciente
en cada columna para un total de 30 registros
En tal sentido se pide:
2a. Grafique la amplitud amortiguada
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
La posición inicial de la masa 2 en el video 2 se ubica en 15,7cm suponiendo que
se empuja hasta una posición en cm y se suelta la masa, podemos verificar que el
sistema como tal alcanza valores comprendidos entre (17,00-17,2)cm. De aquí podemos
obtener las condiciones iniciales para esta parte:
x(0)=(17,2-15,7) cm= 1,5 cm=0,015 m
v̇(0)=0 m/s
Pasamos a encontrar el valor de a partir de los datos suministrados por la tabla
Con los tiempos obtenidos en la medición y que fueron registrados en la tabla, se
usaran para obtener el cuasi periodo 𝛕, tiempo que la masa tarda en alcanzar dos picos
𝛕
consecutivos (los cuales pueden ser un máximo o mínimo). Sea el promedio de las
𝟐
diferencias consecutivas de los tiempos que tarda la masa en pasar por su punto de
equilibrio; es decir:
𝟐𝟖
𝟏
𝛕= ∑(𝒕𝒊+𝟏 − 𝒕𝒊 ) = 𝟐, 𝟎𝟓𝟖𝟔𝟐
𝟐𝟖
𝒊=𝟏
𝛕
En consecuencia 𝟒 = 𝟎, 𝟓𝟏𝟒𝟔𝟓𝟓 será el tiempo que tarda la masa en ir del punto de
equilibrio a un extremo. Ahora bien el cuasi-periodo 𝛕 para el caso sub-amortiguado se
calcula empleando la siguiente relación:
2𝜋
𝛕=
𝜔 0 √ 1 − 𝜀2
β
Donde ∈= ω . Por tanto calculamos ese valor experimental realizando un despeje
o
2 2
2𝜋 2𝜋
∈= √1 − ( ) = √1 − ( ) = 0,973
𝜔0 𝛕 𝟏𝟑, 𝟐𝟗𝟖𝒔−𝟏 ∙ 𝟐, 𝟎𝟓𝟖𝟔𝟐𝟏𝒔
La solución de la EDO con amortiguación libre (ver Apéndice del problema 2) viene
dada por
𝐱(𝐭) = 𝐞−∈ωo 𝐭 (𝐀𝐜𝐨𝐬(𝛚𝐭) + 𝐁𝐬𝐞𝐧(𝛚𝐭))
Y su derivada es
𝐱̇ (𝐭) = −∈ ωo (𝐀𝐜𝐨𝐬(𝛚𝐭) + 𝐁𝐬𝐞𝐧(𝛚𝐭))𝐞−∈ωo 𝐭 + 𝛚(−𝐀𝐬𝐞𝐧(𝛚𝐭) + 𝐁𝐜𝐨𝐬(𝛚𝐭))𝐞−∈ωo 𝐭
Con
𝛚 = 𝛚𝟎 √𝟏 − 𝛆𝟐 = 𝟏𝟑, 𝟎𝟗𝟕 ∙ 𝟎, 𝟏𝟔𝟑𝟒 = 𝟐, 𝟏𝟒𝟎𝟏
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
Empleando las condiciones iniciales del problema x(0)=0,015m y v(0)=0 m/s
𝐀𝐜𝐨𝐬(𝟎) + 𝐁𝐬𝐞𝐧(𝟎) = 𝐱 𝟎
{
−∈ ωo (𝐀𝐜𝐨𝐬(𝟎) + 𝐁𝐬𝐞𝐧(𝟎)) + 𝛚(−𝐀𝐬𝐞𝐧(𝟎) + 𝐁𝐜𝐨𝐬(𝟎)) = 𝟎
Luego de resolver el sistema de ecuaciones anteriores, se obtuvo
∈ωo
𝐱(𝐭) = 𝐞−∈ωo 𝐭 𝐱 𝟎 (𝐜𝐨𝐬(𝛚𝐭) + 𝐬𝐞𝐧(𝛚𝐭)),
𝛚
∈
= 𝐞−∈ωo 𝐭 𝐱 𝟎 (𝐜𝐨𝐬(𝛚𝐭) + 𝐬𝐞𝐧(𝛚𝐭))
√𝟏 − 𝛆𝟐
Se comprueba que
𝐱(𝐭) = 𝟎, 𝟎𝟏𝟓𝐞−12,75𝐭 (𝐜𝐨𝐬(𝟐, 𝟏𝟒𝟎𝟏𝐭) + 𝟓, 𝟗𝟔𝐬𝐞𝐧(𝟐, 𝟏𝟒𝟎𝟏𝐭))𝒎
Para verificar las relaciones anteriores se empleó el siguiente programa en Octave
clc;
clear all;
%Ingresamos los datos del problema
k=143.41; %valor de la constante del resorte
m2=0.836; %valor de la masa 2
x0 = 0.016; % posicion inicial de la masa en t=0 en m
v0 = 0.000; % velocidad inicial de la masa en t =0 en m/s
t2= [0.39 1.28 2.26 3.17 4.21 5.36 6.47 7.58 8.55 9.52 10.53 11.59 12.53 13.59 14.6 15.63 16.78 17.77
18.87 19.92 21 21.95 23.17 24.02 25.13 26.16 27.26 28.28 29.26 30.24]; % tiempos en que la masa 2
pasa por el punto de equilibrio
% Calculo del cuasi-periodo
dif=zeros(1,size(t2,2)-2);
for i=1:size(t2,2)-2;
dif(i)= t2(i+2)-t2(i);
end
Td= sum(dif)/(size(t2,2)-2) % Cuasi-periodo Td
t=(0:Td/10:3*Td);
% calculo de frecuencia angular para la masa 2 a partir de la masa dada en
% el video 1 y la constante elastica calculada en el item 1a
w0 = sqrt(k/m2) % Frecuencia natural masa 2
% Calculo del parametro epsilon
e = sqrt(1-((2*pi)/(w0*Td))^2)
%Calculo de la pseudofrecuencia
w= w0*sqrt(1-e^2)
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
B= e*w0
%Ecuacion de la solucion del movimiento de la masa
y=x0.*exp(-B.*t).*(cosd(-w.*t)+B/w.*sind(-w.*t));
y1=x0.*exp(-B.*t);
%y2=-x0.*exp(-B.*t);
1) Grafica usando el programa Octave
Instrucciones empleadas en Octave para el desarrollo de la grafica
hold on
p1=plot(t,y);
set(p1,'Color','blue','LineWidth',1)
xlabel('t(segundos)');
ylabel('Oscilacion(metros)');
title('Sistema amortiguado libre');
grid on
legend('Especto de frecuencia');
axis([0,3,-0.001,0.017])
figure
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
2b. Considerando que el cuasi-período es el tiempo entre dos máximos (o mínimos)
de la solución de la ecuación de movimiento amortiguada. Determine los dos
coeficientes de dicha solución
A=0,015m
B= 0,0894m
𝛃 = 𝟏𝟐, 𝟕𝟑𝟗 𝐬−𝟏
2c. Obtenga el espectro de amplitud amortiguada
En este caso empleamos la ecuación
𝐲 = 𝐀𝐞−∈𝛚𝟎 𝐭 = 𝟎, 𝟎𝟏𝟓𝐞−𝟏𝟐,𝟕𝟑𝟗 𝐭
Instrucciones en Octave para el desarrollo de la grafica
p2=plot(t,y1);
set(p2,'Color','red','LineWidth',1)
grid on;
xlabel('t(segundos)');
ylabel('Amplitud(metros)');
title('Sistema amortiguado libre');
legend('Amplitud Amortiguada');
axis([0,3,-0.001,0.017])
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
hold off
2d. En base a la respuesta en (2c) determine el valor del coeficiente de
amortiguamiento.
Vemos que
𝐤𝐠
𝐛 = 𝟐𝐦𝛃 = 𝟐(𝟎, 𝟖𝟏𝟏)(𝟏𝟐, 𝟕𝟑𝟗)
𝐬
Se verifica
𝐤𝐠
𝐛 = 𝟐𝟏, 𝟐𝟗𝟗
𝐬
Parte III: Sistema amortiguado Forzado
Considere ahora que obliga al sistema a vibrar armónicamente a una frecuencia
̅𝑜 con una intensidad de
excitadora de 2%, 4% y 10% la frecuencia natural promedio 𝜔
20 (N), y se pide:
3a. Graficar el espectro de amplitud del sistema forzado amortiguado
La frecuencia promedio en este caso vendría dado por la suma de las dos de las
situaciones anteriores
𝑟𝑎𝑑
𝜔1 + 𝜔2 (16,73 + 13,09) 𝑠 𝑟𝑎𝑑
𝜔
̅= = = 14,91
2 2 𝑠
Por tanto los valores de la frecuencia excitadora serán los siguientes
𝐫𝐚𝐝
𝛀𝟏 = 𝟎, 𝟐𝟗𝟖
𝐬
𝐫𝐚𝐝
𝛀𝟐 = 𝟎, 𝟓𝟗𝟔
𝐬
𝐫𝐚𝐝
𝛀𝟑 = 𝟏, 𝟒𝟗𝟏
𝐬
𝐫𝐚𝐝
𝛀𝟒 = 𝟏𝟒, 𝟗𝟏
𝐬
La fuerza forzante se escribirá como sigue
𝐟(𝐭) = 𝐟𝐨 𝐜𝐨𝐬(𝜴𝒊 𝐭)
Con i=1, 2, 3 y 4 que corresponden a los casos de fuerza excitadora. Para la
amplitud del sistema masa resorte con amortiguamiento forzado, se usara la siguiente
ecuación.
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
𝑭𝒐
𝐱(𝐭) = (𝐀𝒄𝒐𝒔(𝝎𝐭) + 𝑩𝒔𝒆𝒏(𝝎𝐭)) + ( ) (𝝎𝟐 𝒄𝒐𝒔(𝛀𝐭) + 𝟐𝜷𝛀𝒄𝒐𝒔(𝛀𝐭))
𝒎((𝝎 ) + (𝟐𝜷𝛀)𝟐 )
𝟐 𝟐
Por otro lado, definimos
𝝎𝟐 = 𝝎𝒐 𝟐 − 𝛀𝟐
El valor de las constantes A y B se realizó de la misma forma como en el caso
del problema anterior.
Presentamos el conjunto de instrucciones en Octave para la generación de la
gráfica y cálculos necesarios
clear;
clear all;
%Datos del problema
w1=16.73;
w2=13.09;
m=0.836;
F0=20;
b=21.299;
x0=0.015;
v0=0.000;
%Parametros a usar
wn=(w1+w2)/2;
f=wn/(2*pi);
T=1/f;
w=wn*0.02;
z1=wn*0.04;
z2=wn*0.1;
z3=wn;
f=F0/m;
B=b/(2*m);
p=2*B*w; %frecuencia forzante al 2% de la frecuencia promedio
p1=2*B*z1; %frecuencia forzante al 4% de la frecuencia promedio
p2=2*B*z2; %frecuencia forzante al 10% de la frecuencia promedio
p3=2*B*z3; %frecuencia forzante igual a la frecuencia promedio
e=B/wn;
wd=wn*sqrt(1-(B/wn)^2); %frecuencia natural del sistema
s=sqrt(wn^2-w^2);
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
s1=sqrt(wn^2-z1^2);
s2=sqrt(wn^2-z2^2);
s3=sqrt(wn^2-z3^2);
%Ecuaciones a usar
t= 0:T/30:T;
M= f*(s^2)/(s^4+p^2);
M1= f*(s1^2)/(s1^4+p1^2);
M2= f*(s2^2)/(s2^4+p2^2);
M3 = f*(s3^2)/(s3^4+p3^2);
%Calculo de los coeficientes
A=x0-M;
C=(1/wd)*(v0+B*A-f*(p*w)/(s^4+p^2));
L=f/(s^4+p^2);
A1=x0-M1;
C1=(1/wd)*(v0+B*A1-f*(p1*s1)/(s1^4+p1^2));
L1=f/(s1^4+p1^2);
A2=x0-M2;
C2=(1/wd)*(v0+B*A2-f*(p2*s2)/(s2^4+p2^2));
L2=f/(s2^4+p2^2);
A3=x0-M3;
C3=(1/wd)*(v0+B*A3-f*(p3*s3)/(s3^4+p3^2));
L3=f/(s3^4+p3^2);
%Posicion
x=(A*cos(wd.*t)+C*sin(wd.*t)).*exp(-B.*t)+L*(s^2*cos(w.*t)+p*sin(w.*t));
x1=(A1*cos(wd.*t)+C1*sin(wd.*t)).*exp(-B.*t)+L1*(s1^2*cos(z1.*t)+p1*sin(z1.*t));
x2=(A2*cos(wd.*t)+C2*sin(wd.*t)).*exp(-B.*t)+L2*(s2^2*cos(z2.*t)+p2*sin(z2.*t));
x3=(A3*cos(wd.*t)+C3*sin(wd.*t)).*exp(-B.*t)+L3*(s3^2*cos(z3.*t)+p3*sin(z3.*t));
%Graficas
hold on
p=plot(t,x);
set(p,'Color','blue','LineWidth',1)
legend('Frecuencia forzante 2%')
p1=plot(t,x1);
set(p1,'Color','red','LineWidth',1)
p2=plot(t,x2);
set(p2,'Color','green','LineWidth',1)
p3=plot(t,x3);
set(p3,'Color','black','LineWidth',1)
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
xlabel('Tiempo (s)')
ylabel('Amplitud(m)')
title('Amplitud vs Tiempo')
grid on
axis([0,0.30,-0.01,0.15])
hold off
No se pudo verificar directamente la forma de la campana de Gauss en las tres
frecuencias
3b. Teóricamente en cuanto se afecta la amplitud en la energía respecto a la
amplitud del desplazamiento. En base a ello, puede determinar el efecto de ese
porcentaje en el desplazamiento e indicar su valor.
Por la presencia del término atenuador (el proporcional a la velocidad) entonces
esto afectara el circuito con la perdida de energía a través del rozamiento.
La relación de la energía con la amplitud viene dado por la expresión:
𝐅𝐨 𝟐𝛑𝛆𝛀
𝐖= ( )
𝐦 (√(𝛚𝐨 𝟐 − 𝛀𝟐 )𝟐 + (𝟐𝛃𝛀)𝟐 )
Para cada frecuencia se tiene
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
𝐅𝐨 𝟐𝛑𝛆𝛀
𝑾𝟏 = ( ) = 𝟐, 𝟔𝟏 𝐉
𝐦 (√(𝛚𝐨 𝟐 − 𝛀𝟐 )𝟐 + (𝟐𝛃𝛀)𝟐 )
𝐅𝐨 𝟐𝛑𝛆𝛀𝟏
𝐖𝟐 = = 𝟐, 𝟓𝟗 𝐉
𝟐 𝟐
𝐦 (√(𝛚𝐨 𝟐 − 𝛀𝟏 ) + (𝟐𝛃𝛀𝟏 )𝟐 )
( )
𝐅𝐨 𝟐𝛑𝛆𝛀𝟐
𝐖𝟑 = = 𝟐, 𝟓𝟒 𝐉
𝟐 𝟐
𝐦 (√(𝛚𝐨 𝟐 − 𝛀𝟐 ) + (𝟐𝛃𝛀𝟐 )𝟐 )
( )
La instrucción en Octave para el cálculo de la energía viene dada por
%Parametros a usar
U=(F0/m)*((2*pi*e*p)/(sqrt((wn^2-p^2)^2+(2*B*p)^2))) %Energia al 2%
U1=(F0/m)*((2*pi*e*p1)/(sqrt((wn^2-p1^2)^2+(2*B*p1)^2))) %Energia al 4%
U2=(F0/m)*((2*pi*e*p2)/(sqrt((wn^2-p2^2)^2+(2*B*p2)^2))) %Energia al 10%
3c. ¿Puede determinar el tiempo en que el sistema alcanza el 6% de
amortiguamiento? En caso afirmativo indique el valor.
La forma de determinar ese valor es comprando el punto donde la amplitud ha
decaído en un 6%
3d. Determine el valor de la amplitud cuando son iguales frecuencia natural promedio y
frecuencia excitadora.
𝑭𝒐
𝐱(𝐭) = ( )
𝒎 (√(𝝎𝒐 𝟐 − 𝜴𝟐 )𝟐 + (𝟐𝜷𝜴)𝟐 )
En este caso se produce la resonancia y ella ocurre cuando un sistema vibrante
es forzado con una frecuencia que produce la mayor amplitud posible. Generalmente
esto ocurre cuando la frecuencia de forzado es igual a la frecuencia natural. En otras
palabras
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
(𝝎𝒐 𝟐 − 𝛀𝟐 )𝟐 + (𝟐𝜷𝛀)𝟐 = 𝟎
Y usando la relación
𝝎𝟐 = 𝝎𝒐 𝟐 − 𝛀𝟐
Por tanto el fenómeno de resonancia y la situación con la mayor amplitud posible
se alcanza cuando
𝝎𝟐 = −(𝟐𝜷𝛀)𝟐
Discusión de resultados
Problema 1
El desarrollo del problema tuvo como finalidad hacer una descripción y
verificación de los datos obtenidos luego de observar los videos 1 y 2. En el caso del
video 1, se comprobó el sistema de oscilación libre. En este sentido se determinó el
valor de la constante del resorte aplicando el siguiente procedimiento. En primer lugar
se colocó la barra en reposo en una posición de equilibrio. Luego se colocó la masa 1
para determinar cuánto se bajaba la barra para encontrar su nueva posición de
equilibrio. Aplicando la segunda ley de Newton se determinó la constante del resorte
teniendo presente la ley de Hooke, además de considerar el límite elástico y que el
cuerpo en esa nueva posición permanece en reposo. Se realizó ese mismo
procedimiento para la segunda masa y se pudo comprobar que ambos valores para la
constante elástica son parecidos.
Se verifico y propuso un diseño experimental sencillo para encontrar el valor de
la frecuencia natural del sistema. Todo lo anterior se demostró numéricamente
empleando el programa Octave, por medio del cual se hizo una programación sencilla
para verificar todo lo anterior.
Problema 2
El desarrollo de este problema coincidió con la solución de un sistema masa
resorte con amortiguación libre. En este caso, partiendo de los tiempos obtenidos de la
medición de las diferentes oscilaciones (datos mostrados en la tabla), se procedió a
determinar el cuasi-período (el tiempo entre dos máximos o mínimos) con el fin de
poder encontrar un factor llamado épsilon el cual relacionaba la frecuencia natural del
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
sistema con el factor de atenuación.
Seguidamente, se determinaron los valores de b (constante multiplicativa de la
rapidez), y las constantes de integración A y B, que acompañaban a las funciones
trigonométricas, partiendo de las condiciones iniciales del problema para tener por tanto
la solución completa del problema. Se pudo verificar que este sistema se asemejaba un
poco al sistema subamortiguado que los valores de la frecuencia natural y la frecuencia
proporcionada por la amortiguación del sistema no diferían en mucho, por lo que la
oscilación se atenuaba muy rápido.
Problema 3
El desarrollo de este problema estableció la posibilidad de trabajar el sistema
masa-resorte con la presencia de un término forzante. De aquí se pudo verificar el caso
de resonancia y la disipación de energía producto a la fuerza de roce. En la fase de
resonancia la velocidad se encuentra en fase con la fuerza oscilante aplicada y de allí
se desprende que la energía de la fuerza que se aplica sobre el oscilador forzado tiene
un valor máximo.
La no observación de la campana en la amplitud en el problema puede ser
debido a los siguientes factores, error en la sustitución de los valores dentro de las
ecuaciones o uso incorrecto de algunas de ellas al emplear Octave.
Conclusiones
Se realizó el análisis teórico y revisión de los datos prácticos de tres problemas
relacionados con sistemas oscilatorios, específicamente el de masa-resorte
considerando los casos siguientes: oscilación libre, oscilación libre amortiguada y de
oscilación amortiguada forzada. En los dos últimos casos se pudo comprobar que se
trabajó con la situación del subamortiguado con idénticas condiciones iniciales; liberado
desde el reposo y recorriendo una distancia de 1,5 cm por debajo del punto de equilibrio
x(0) = 0,015m y del reposo x0(0) = 0m/s. Luego de realizar los calculo
correspondientes, el sistema pudo evolucionar mas a uno con el caso de
amortiguamiento critico que con subamortiguado, tal y como se pudo verificar en las
gráficas generadas en Octave.
Para el caso de las dos últimas situaciones que la amplitud máxima fue la inicial
y (excepto en el caso de resonancia en la oscilación amortiguada forzada) y en donde
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
para el caso de la oscilación libre amortiguada, pues dicha amplitud decrecía mientras
el tiempo pasaba.
El desarrollo de esta asignación permitió verificar lo obtenido teóricamente con lo
observado en los videos 1 y 2 con el fin de desarrollar competencias en el momento de
extraer información de fuentes verificables como los videos o las tablas. Finalmente, el
aprendizaje significativo se verifica luego de confrontar los datos experimentales con
elementos concretos para un buen análisis de lo observado.
Apéndice
Problema 1
La ecuación diferencial a resolver queda escrita como sigue:
𝐦𝒙̈ + 𝐤𝐱 = 𝟎
La cual escribimos
𝐤
𝒙̈ + 𝐱=𝟎
𝒎
Y hacemos las sustituciones
𝐤
𝝎𝒐 𝟐 =
𝒎
Para obtener la ecuación diferencial
𝒙̈ + 𝝎𝒐 𝟐 𝐱 = 𝟎
En este caso hacemos el cambio
𝐱 = 𝐀𝐞𝐤𝐭
Lo que permite transformarla en una ecuación de la forma
𝐤 𝟐 + 𝛚𝐨 𝟐 = 𝟎
La cual conduce a la solucion
𝐱(𝐭) = (𝐀𝐜𝐨𝐬(𝛚𝐭) + 𝐁𝐬𝐞𝐧(𝛚𝐭))
El cálculo de la constante K en ambos videos se procedió teóricamente de la
siguiente forma
𝐦𝐠 − 𝐤𝚫𝐱 = 𝟎
Problema 2
La ecuación diferencial a resolver queda escrita como sigue:
𝐦𝒙̈ + 𝐛𝒙̇ + 𝐤𝐱 = 𝟎
La cual escribimos
UNIVERSIDAD CENTRAL DE VENEZUELA * FACULTAD DE INGENIERÍA * DEPARTAMENTO DE FÍSICA APLICADA
POSTGRADO EN MECÁNICA TEÓRICA Y APLICADA
TEORÍA LINEAL DE VIBRACIONES (80803010)
SEMESTRE LECTIVO 3-2020
𝒃 𝐤
𝒙̈ + 𝒙̇ + 𝐱 = 𝟎
𝒎 𝒎
Y hacemos las sustituciones
𝐛 𝐤
𝛃 = 𝟐𝐦, 𝝎𝒐 𝟐 = 𝒎
Para obtener la ecuación diferencial
𝒙̈ + 𝟐𝜷𝒙̇ + 𝝎𝒐 𝟐 𝐱 = 𝟎
En este caso hacemos el cambio
𝐱 = 𝐀𝐞𝐤𝐭
Lo que permite transformarla en una ecuación de la forma
𝐤 𝟐 + 𝟐𝛃𝐤 + 𝛚𝐨 𝟐 = 𝟎
La idea es encontrar las raíces de esa ecuación
−𝟐𝛃 ± √𝟒𝛃𝟐 − 𝟒𝛚𝐨 𝟐
𝐤=
𝟐
La cual se reduce a
𝐤 = −𝛃 ± √𝛃𝟐 − 𝛚𝐨 𝟐
Problema 3
La fuerza forzante tiene la forma
𝐟(𝐭) = 𝐟𝐨 𝐜𝐨𝐬(𝛀𝐭)
La ecuación diferencial a resolver queda escrita como sigue:
𝐦𝒙̈ + 𝐛𝒙̇ + 𝐤𝐱 = 𝐟(𝐭)
La solución de la ecuación anterior
𝐱(𝐭) = 𝐱 𝐩 + 𝐱 𝐡
La solución homogénea (xh) viene dada por
𝐱(𝐭) = 𝐞−𝛃𝐭 (𝐀𝐜𝐨𝐬(𝛚𝐭) + 𝐁𝐬𝐞𝐧(𝛚𝐭))
La solución particular (xp) teniendo presente el término forzante es:
𝐱(𝐭) = 𝐂𝐞𝐢𝛀𝐭 ,
Donde C es un coeficiente a determinar, el cual es imaginario