SISTEMA MASA-RESORTE-FRICCIÓN CON DOS GRADOS DE LIBERTAD

...

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;