Scribd
Cargar
29 vistas7 páginas

Inicio de Program1

El documento es un código en MATLAB que simula el movimiento de un mecanismo de manivela y biela, calculando ángulos, velocidades, aceleraciones y fuerzas a lo largo de 360 grados. Utiliza entradas como la longitud de la manivela, la longitud de la biela, y las masas de ambos componentes para realizar cálculos y generar gráficos que representan las relaciones entre las variables. Los resultados se visualizan en múltiples figuras que muestran las variaciones de ángulos, velocidades, aceleraciones y fuerzas en función del ángulo de la manivela.

Cargado por

Rodrigo Diaz Meza
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
29 vistas7 páginas

Inicio de Program1

El documento es un código en MATLAB que simula el movimiento de un mecanismo de manivela y biela, calculando ángulos, velocidades, aceleraciones y fuerzas a lo largo de 360 grados. Utiliza entradas como la longitud de la manivela, la longitud de la biela, y las masas de ambos componentes para realizar cálculos y generar gráficos que representan las relaciones entre las variables. Los resultados se visualizan en múltiples figuras que muestran las variaciones de ángulos, velocidades, aceleraciones y fuerzas en función del ángulo de la manivela.

Cargado por

Rodrigo Diaz Meza
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

% Inicio de programa

clear all
clc
lm=input('lm:');
lb=input('lb:');
n=360;
V1=input('Velocidad de la manibela:');
mm=input('Masa de la manibela:');
mb=input('Masa de la biela:');
mc=mb/3;
aa=360/n;
i=1;
g=9.81;
uk=0.1;
j=1;
na=0;
while j<=n
om=((2*pi)/360)*na ;
r=((lm/lb)*sin(om));
if (r<=1&&r>=-1)
ob=(360/(2*pi))*(asin((lm/lb)*sin(om)));
Om(i,j)=na;
Ob(i,j)=ob;
S(i,j)=(lm*cos(om)+(lb*cos(ob*((2*pi)/360))));
else
Ob(i,j)=0;
S(i,j)=0;
end
%calculo de velocidades
Vm(i,j)=V1;
o1=om;
o2=r;
M1=[lb*cos(o2) (0);
lb*sin(o2) (1)];
kv=[(lm*cos(o1)*V1);
(-lm*sin(o1)*V1)];
V=inv(M1)*kv;
Vb(i,j)=V(1);
Vs(i,j)=V(2);
%calculo de aceleracione
Am(i,j)=0;
M2=[(lb*cos(o2)) (0);
(lb*sin(o2)) (1)];
ka=[(-lm*sin(o1)*(V1^2))+(lm*cos(o1)*0)+(lb*sin(o2)*((V(1))^2));
(-lm*cos(o1)*(V1^2))-(lm*sin(o1)*0)-(lb*cos(o2)*((V(1))^2))];
A=inv(M2)*ka;
Ab(i,j)=A(1);
As(i,j)=A(2) ;
%grafico de la fuerza
if (om>=(pi/2) && om<=pi)
F(i,j)=((5000/(pi/2))*(om-(pi/2)));
else
if (om>pi && om<=((pi)*(3/2)))
F(i,j)=(5000/(-(pi/2))*(om-((3*pi)/2)));
else
F(i,j)=0;
end
end
%análisis inercial
amx=-0.5*lm*(V1^2)*cos(o1);
amy=-0.5*lm*(V1^2)*sin(o1);
abx=0;
aby=0;
Im=((1/12)*mm*(lm^2));
Ib=((1/12)*mb*(lb^2));
X=((Im*0)-((lm/2)*sin(o1)*mm*amx)+((lm/2)*cos(o1)*mm*amy)+
((lm/2)*mm*g*cos(o1)));
Y=((-Ib*A(1))-((lb/2)*sin(o2)*mb*abx)-((lb/2)*cos(o2)*mb*aby)-
((lb/2)*mb*g*cos(o2)));
M3=[1 0 1 0 0 0 0 0;
0 1 0 1 0 0 0 0;
0 0 (-lm*sin(o1)) (lm*cos(o1)) 0 0 0 -1;
0 0 -1 0 1 0 0 0;
0 0 0 -1 0 1 0 0;
0 0 (lb*sin(o2)) (lb*cos(o2)) 0 0 0 0;
0 0 0 0 -1 0 uk 0;
0 0 0 0 0 -1 1 0];
k3=[(mm*amx);
((mm*amy)+(mm*g));
(X);
(mb*abx);
((mb*aby)+(mb*g));
(Y);
((mc*A(2))+F(j));
(mc*g)];
B=inv(M3)*k3;
Ax(i,j)=B(1);
Ay(i,j)=B(2);
Bx(i,j)=B(3);
By(i,j)=B(4);
Cx(i,j)=B(5);
Cy(i,j)=B(6);
Nc(i,j)=B(7);
Ma(i,j)=B(8);
j=j+1;
na=na+aa;
end
Ob
S
figure(1)
subplot(121)
plot(Om,Ob,'b')
title('Om vs Ob')
xlabel('Angulo de la manivela')
ylabel('Angulo de la biela')
subplot(122)
plot(Om,S,'g')
title('Om vs S')
xlabel('Angulo de la manivela')
ylabel('Distancia S')
figure(2)
subplot(121)
plot(Om,Vb,'b')
title('Om vs Ob')
xlabel('Angulo de la manivela')
ylabel('Velocidad de la biela')
subplot(122)
plot(Om,Vs,'g')
title('Om vs S')
xlabel('Angulo de la manivela')
ylabel('Velocidad S')
figure(3)
subplot(121)
plot(Om,Ab,'b')
title('Om vs Ob')
xlabel('Angulo de la manivela')
ylabel('aceleracion de la biela')
subplot(122)
plot(Om,As,'g')
title('Om vs S')
xlabel('Angulo de la manivela')
ylabel('aceleracion S')
figure(4)
subplot(121)
plot(Om,F,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza F')
subplot(122)
plot(Om,Ma,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza F')
figure(5)
subplot(121)
plot(Om,Ax,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Ax')
subplot(122)
plot(Om,Ay,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Ay')
figure(6)
subplot(121)
plot(Om,Bx,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Bx')
subplot(122)
plot(Om,By,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza By')
figure(7)
subplot(121)
plot(Om,Cx,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Cx')
subplot(122)
plot(Om,Cy,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Cy')
figure(8)
plot(Om,Nc,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Nc')
% Inicio de programa
clear all
clc
lm=input('lm:');
lb=input('lb:');
n=360;
V1=input('Velocidad de la manibela:');
mm=input('Masa de la manibela:');
mb=input('Masa de la biela:');
mc=mb/3;
aa=360/n;
i=1;
g=9.81;
uk=0.1;
j=1;
na=0;
while j<=n
om=((2*pi)/360)*na ;
r=((lm/lb)*sin(om));
if (r<=1&&r>=-1)
ob=(360/(2*pi))*(asin((lm/lb)*sin(om)));
Om(i,j)=na;
Ob(i,j)=ob;
S(i,j)=(lm*cos(om)+(lb*cos(ob*((2*pi)/360))));
else
Ob(i,j)=0;
S(i,j)=0;
end
%calculo de velocidades
Vm(i,j)=V1;
o1=om;
o2=r;
M1=[lb*cos(o2) (0);
lb*sin(o2) (1)];
kv=[(lm*cos(o1)*V1);
(-lm*sin(o1)*V1)];
V=inv(M1)*kv;
Vb(i,j)=V(1);
Vs(i,j)=V(2);
%calculo de aceleracione
Am(i,j)=0;
M2=[(lb*cos(o2)) (0);
(lb*sin(o2)) (1)];
ka=[(-lm*sin(o1)*(V1^2))+(lm*cos(o1)*0)+(lb*sin(o2)*((V(1))^2));
(-lm*cos(o1)*(V1^2))-(lm*sin(o1)*0)-(lb*cos(o2)*((V(1))^2))];
A=inv(M2)*ka;
Ab(i,j)=A(1);
As(i,j)=A(2) ;
%grafico de la fuerza
if (om>=(pi/2) && om<=pi)
F(i,j)=((5000/(pi/2))*(om-(pi/2)));
else
if (om>pi && om<=((pi)*(3/2)))
F(i,j)=(5000/(-(pi/2))*(om-((3*pi)/2)));
else
F(i,j)=0;
end
end
%análisis inercial
amx=-0.5*lm*(V1^2)*cos(o1);
amy=-0.5*lm*(V1^2)*sin(o1);
abx=0;
aby=0;
Im=((1/12)*mm*(lm^2));
Ib=((1/12)*mb*(lb^2));
X=((Im*0)-((lm/2)*sin(o1)*mm*amx)+((lm/2)*cos(o1)*mm*amy)+
((lm/2)*mm*g*cos(o1)));
Y=((-Ib*A(1))-((lb/2)*sin(o2)*mb*abx)-((lb/2)*cos(o2)*mb*aby)-
((lb/2)*mb*g*cos(o2)));
M3=[1 0 1 0 0 0 0 0;
0 1 0 1 0 0 0 0;
0 0 (-lm*sin(o1)) (lm*cos(o1)) 0 0 0 -1;
0 0 -1 0 1 0 0 0;
0 0 0 -1 0 1 0 0;
0 0 (lb*sin(o2)) (lb*cos(o2)) 0 0 0 0;
0 0 0 0 -1 0 uk 0;
0 0 0 0 0 -1 1 0];
k3=[(mm*amx);
((mm*amy)+(mm*g));
(X);
(mb*abx);
((mb*aby)+(mb*g));
(Y);
((mc*A(2))+F(j));
(mc*g)];
B=inv(M3)*k3;
Ax(i,j)=B(1);
Ay(i,j)=B(2);
Bx(i,j)=B(3);
By(i,j)=B(4);
Cx(i,j)=B(5);
Cy(i,j)=B(6);
Nc(i,j)=B(7);
Ma(i,j)=B(8);
j=j+1;
na=na+aa;
end
Ob
S
figure(1)
subplot(121)
plot(Om,Ob,'b')
title('Om vs Ob')
xlabel('Angulo de la manivela')
ylabel('Angulo de la biela')
subplot(122)
plot(Om,S,'g')
title('Om vs S')
xlabel('Angulo de la manivela')
ylabel('Distancia S')
figure(2)
subplot(121)
plot(Om,Vb,'b')
title('Om vs Ob')
xlabel('Angulo de la manivela')
ylabel('Velocidad de la biela')
subplot(122)
plot(Om,Vs,'g')
title('Om vs S')
xlabel('Angulo de la manivela')
ylabel('Velocidad S')
figure(3)
subplot(121)
plot(Om,Ab,'b')
title('Om vs Ob')
xlabel('Angulo de la manivela')
ylabel('aceleracion de la biela')
subplot(122)
plot(Om,As,'g')
title('Om vs S')
xlabel('Angulo de la manivela')
ylabel('aceleracion S')
figure(4)
subplot(121)
plot(Om,F,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza F')
subplot(122)
plot(Om,Ma,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza F')
figure(5)
subplot(121)
plot(Om,Ax,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Ax')
subplot(122)
plot(Om,Ay,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Ay')
figure(6)
subplot(121)
plot(Om,Bx,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Bx')
subplot(122)
plot(Om,By,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza By')
figure(7)
subplot(121)
plot(Om,Cx,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Cx')
subplot(122)
plot(Om,Cy,'g')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Cy')
figure(8)
plot(Om,Nc,'b')
title('Om vs F')
xlabel('Angulo de la manivela')
ylabel('Fuerza Nc')

También podría gustarte