SISTEMA MASA-RESORTE-FRICCIÓN CON DOS GRADOS DE LIBERTAD
Enviado por Rimma • 15 de Noviembre de 2017 • 1.432 Palabras (6 Páginas) • 601 Visitas
...
L = 6.0; %ámbito de elevación permisible en m
dP = 600e-3; %diámetro del tambor de cable en m
mG = 3.35; %relación de transmisión (r2)/(r1)
g=9.81;%gravedad terreste
m1=me+mc;
A=pi*((dc/2)^2);
%-----------------------------------------------------------------------
l=L; %l cte
k=A*E/l; %Resorte lineal(con l cte)
%Constantes totales
JT=J1+(J2/(mG^2));
BOT=bo1+(bo2/(mG^2));
KOT=k*((dP/2)^2)/(mG^2);
C=k*dP/(2*mG);
%Con T=640Nm
T=640;
dzdt(1)=z(2);
dzdt(2)=(1/JT)*(-BOT*z(2)-KOT*z(1)+T+C*z(3));
dzdt(3)=z(4);
dzdt(4)=(1/(m1+m2))*(-bL*z(4)-k*z(3)-(m1-m2)*g+C*z(1));
dzdt=dzdt'; %trasposición de filas a columnas
end
Archivo lab4b.m
function dzdt=lab4b(t,z)
%----------------------------------------------------------------------
%Parte genérica de todos los items
%Constantes dadas, ENTRADAS del programa
me = 350; %masa del elevador en kg
mc = 850; %máxima carga que puede llevar el elevador en kg
m2 = 750; %masa del contrapeso en kg
bL = 860; %fricción en rieles del ascensor en N/(m/s)
J1 = 1.12; %inercia rotacional del motor + engrane 1 en kg*m^2
J2 = 4.34; %inercia rotacional del tambor + engrane 2 en kg*m^2
bo1 = 5.58e-2; %fricción en cojinetes del motor en N*m/(rad/s)
bo2 = 17.4e-2; %fricción en cojinetes del tambor en N*m(rad/s)
dc = 12e-3; %diámetro de los cables de elevador en m
E = 200e9; %módulo de elasticidad del acero en Pa
L = 6.0; %ámbito de elevación permisible en m
dP = 600e-3; %diámetro del tambor de cable en m
mG = 3.35; %relación de transmisión (r2)/(r1)
g=9.81;%gravedad terreste
m1=me+mc;
rp=dP/2;
A=pi*((dc/2)^2);
%-----------------------------------------------------------------------
l=L; %l cte
k=A*E/l; %Resorte lineal(con l cte)
%Constantes totales
JT=J1+(J2/(mG^2));
BOT=bo1+(bo2/(mG^2));
KOT=k*((dP/2)^2)/(mG^2);
C=k*dP/(2*mG);
%Con Tmotor(t) para v1=0,8m/s, recibida de "lab4ejecuta"
% global Tmotor;
% Corte de la cuerda con m2 en t=5s
if t>5
m2=0;
end
Archivo lab4d.m
function dzdt=lab4d(t,z)
%----------------------------------------------------------------------
%Parte genérica de todos los items
%Constantes dadas, ENTRADAS del programa
me = 350; %masa del elevador en kg
mc = 850; %máxima carga que puede llevar el elevador en kg
m2 = 750; %masa del contrapeso en kg
bL = 860; %fricción en rieles del ascensor en N/(m/s)
J1 = 1.12; %inercia rotacional del motor + engrane 1 en kg*m^2
J2 = 4.34; %inercia rotacional del tambor + engrane 2 en kg*m^2
bo1 = 5.58e-2; %fricción en cojinetes del motor en N*m/(rad/s)
bo2 = 17.4e-2; %fricción en cojinetes del tambor en N*m(rad/s)
dc = 12e-3; %diámetro de los cables de elevador en m
E = 200e9; %módulo de elasticidad del acero en Pa
L = 6.0; %ámbito de elevación permisible en m
dP = 600e-3; %diámetro del tambor de cable en m
mG = 3.35; %relación de transmisión (r2)/(r1)
g=9.81;%gravedad terreste
m1=me+mc;
A=pi*((dc/2)^2);
%-----------------------------------------------------------------------
k=A*E/L-z(3); %Resorte progresivo(con k cte)
%Constantes totales
JT=J1+(J2/(mG^2));
BOT=bo1+(bo2/(mG^2));
KOT=k*((dP/2)^2)/(mG^2);
C=k*dP/(2*mG);
%Con T=520Nm
T=520;
...