Análisis Dinámico De Un Eslabón

 Análisis dinámico de un eslabón MATLAB
Un eslabón en rotación pura tiene una longitud de 15 cm (del punto O2 al punto P) y su centro de masa está ubicado justo a la mitad de la distancia entre elpunto O2 y el punto P. El eslabón tiene una masa de m2 = 5 kg y un momento de inercia de IG2 = 9.375 x 10-3 kg m2.
Considere que el eslabón en el tiempo t = 0 s tienen un ángulo Ɵ2 = 0°, unavelocidad angular de 𝜔2 = 376 rad/s y una desaceleración constante de α2 = -75.2 rad/s2. Durante el movimiento, el eslabón está sujeto a la acción de una fuerza 𝑃⃗ ={230[cos(𝜔𝑡)−1]}𝑖 +0𝑗
Obtenga la gráficade las componentes rectangulares de la fuerza ejercida por el perno en O2 y la gráfica del par requerido para mantener este movimiento durante el primer giro completo del eslabón.
Nota: Primero debecalcularse la aceleración lineal del punto G2.


% Análisis Dinámico de un Eslabón
% 29 Septiembre 2015
% Dinámica de Máquinas y Herramientas
% Luis Manuel Palacios
%% Datos iniciales
L=0.075;% [m] Distancia al centro de masa
m=5; % [kg] Masa del eslabón
Ig=9.375e-3; % [kg-m^2] Inercia Rotacional del Eslabón
omega_0=376; % [rad/s] Velocidad Angular Inicial
alfa=-75.2;% [rad/s^2] Aceleración Angular
%% Calculo de la aceleración lineal del Punto G
teta=(0:1:359);
teta_rad=pi()*teta/180;
a_n=omega.^2*L;
omega=sqrt(omega_0^2+2*alfa*teta_rad);
a_n=omega.^2*L;plot(teta,omega);
a_t=L*alfa;
a_x=a_n.*cos(teta_rad+pi)+abs(a_t).*cos(teta_rad+3*pi/2);
a_y=a_n.*sin(teta_rad+pi)+abs(a_t).*sin(teta_rad+3*pi/2);
%% Vectores de Posición
R12_x=L.*cos(teta_rad+pi);R12_y=L.*sin(teta_rad+pi);
Rp_x=L.*cos(teta_rad);
Rp_y=L.*sin(teta_rad);
%% Obtención de la fuerza sobre el Eslabón
Fy=0;
Fx=230*(cos(teta_rad)-1);
% Cálculo de las Fuerzas en el Perno O_2
F12_x=m.*a_x-Fx;F12_y=m.*a_y-Fy;
%% Obtención del Par Requerido Para Obtener el Movimiento
T=Ig*alfa-(Rp_x.*Fy-Rp_y.*Fx)+F12_x.*R12_y-F12_y.*R12_x;
%% Grafica de Resultados
figure (1)
plot(teta,F12_x,teta,F12_y)...