rotorconstantesolve

Prévia do material em texto

09/05/20 20:50 C:\Users\ABNER\D...\ROTORCONSTANTESOLVE.m 1 of 2
 clear all
 clc
 
 tspan=[0 2.0];
 y0=[0;0;0;0];
 [t,x]=ode45('ROTORCONSTANTE',tspan,y0);
 
% Dados do equipamento
D = 12 * 10^-3; % diâmetro do eixo [m]
d = 100* 10^-3; % diâmetro do disco [m]
Le = 850 * 10^-3; % comprimento do eixo [m]
Ld = 15 * 10^-3; % largura do disco [m]
r1 = D/2; % raio do eixo 
r2 = d/2; % raio do disco
% Dimenões Físicas do conjunto
me = 1.5 * 10^-3; % momento de excentricidade
ds = 7850; % Kg/m^3 densidade
v1 = pi*r1^2*Le ; %m^3 volume do eixo
v2 = pi*r2^2*Ld ; %m^3 volume do disco
m1 = ds*v1 ; % Kg massa do eixo
m2 = ds*v2; % Kg massa do disco
m =m2; % Massa 
% Constantes do Conjunto
E =200 * 10^9 ; % 
I = pi*(D^4 ) / 64; %Inércia do Eixo
k = 48*E*I /( Le^3) ; % Rigidez do Eixo
c = 10^-3 * k; % Amortecimento do Eixo
g = 9.804;
Wn = sqrt(k/m);
e = me/m;
 
% Coeficientes de multiplicação
n = 1.0;
% Colocar a mesma ressonância da Function aqui
W = n*Wn ;
 
 
 
 
figure(1)
plot3( x(:,1),-t,x(:,3), x(:,1) + e*cos(W*t),-t,x(:,3) + e*sin(W*t) )
grid
xlabel('tempo')
ylabel(' W e G Point')
title(' Fenômenos de Centragem ')
legend('W','G')
hold on; 
 
figure(2)
plot( t, x(:,1),t, x(:,1) + e*cos(W*t) )
grid
xlabel('tempo')
ylabel(' W e G Point ')
title(' Fenômenos de Centragem Uy e Wy Lateral ')
legend('W','G')
09/05/20 20:50 C:\Users\ABNER\D...\ROTORCONSTANTESOLVE.m 2 of 2
hold on; 
 
figure(3)
plot( t, x(:,3),t, x(:,3) + e*sin(W*t) )
grid
xlabel('tempo')
ylabel(' W e G Point')
title(' Fenômenos de Centragem Uz e Wz Lateral')
legend('W','G')
hold on; 
 
figure(4)
plot( x(:,1), x(:,3) , x(:,1) + e*cos(W*t) , x(:,3) + e*sin(W*t) )
grid
xlabel('tempo')
ylabel(' W e G Point')
title(' Fenômenos de Centragem Frontal ')
legend('W','G')
hold on;