Trabalho de Máquinas Parte

[pic 1]

RELATÓRIO TÉCNICO

DISCIPLINA: Máquinas Elétricas I

  2º semestre de 2015

Leonardo Henrique Silva     120950063

Professor: Marco Túlio Alves Êvo        

São João del-Rei

1.  de Setembro de 2015

1. Algoritmo feito no Matlab.

clc;

clear all;

close all;

%% Trabalho de Máquinas Parte 1

% Rotina para análise dos parâmetros de um MIT

% Dados do motor

Peixo = 1.5;            % Potência no eixo em Hp

f = 60;                 % Frequência de alimentação em Hz

n = 1700;               % Velocidade nominal em Rpm

In = 5;                 % Corrente nominal em A

v = 220;                % Tensão nominal em V

np = 4;                 % Número de polos

nf = 3;                 % Número de fases

Ip_In = 5.1;            % Relação corrente de partida e nominal

Io = 3 ;                % Corrente em A

Re = 2.5 ;              % Resistencia do estator em Ohm

Rr = 1.76 ;             % Resistencia do rotor em Ohm

Xm = 43;                % Reatância de magnetização em Ohm

Xr = 1.3;               % Reatância do rotor em Ohm

Xe = 1.3;               % Reatância do estator em Ohm

Prot = 100;             % Perdas rotacionais em W

% Cálculos

s = 0.001:0.0002:1;

Z2 = Rr./s + j*Xr;

Z2n = (Z2 * j*Xm)./(Z2 + j*Xm);

Zt = Re + j*Xe + Z2n;

V1 = v / sqrt(3);

Ie = V1./Zt;

ns = (120/np) * f;

ws = (ns/f)*2*pi;

n = (1-s).*ns;

% Equivalente de Thevenin do estator

Zth = ((Re+j*Xe)*(j*Xm))/(Re+j*(Xe+Xm));

Rth = real(Zth);

Xth = imag(Zth);

Vth = abs(V1*(j*Xm/sqrt((Re^2+j*(Xe+Xm)^2))));

% Calculo do Torque

Ir = abs(Vth./sqrt((Rth+ Rr./s).^2 + (Xr+Xth)^2));

Tmec = (nf*Ir.^2*Rr./s)./ws;

modIe = abs(Ie);

modIr = abs(Ir);

Imag = (V1 - Ie.*(Re + j*Xe))/(j*Xm);

modImag = abs(Imag);

Pg = 3*Ir.^2*Rr./s;

Wm = (1-s).* ws;

Peixo = (1-s).* Pg - Prot;

Teixo = Peixo./Wm;

Fp = cos(angle(Zt));

Pent = 3* V1 .* modIe .* Fp;

Rend = Peixo./Pent;

 

% Plotando os gráficos

figure;

plot(n,Tmec);

grid on;

xlabel('Velocidade [rpm]');

ylabel('Torque [N.m]');

title('Torque x Velocidade','FontSize',13);

 

figure;

plot(n,modIe);

grid on;

xlabel('Velocidade [rpm]');

ylabel('Ie [A]');

title('Corrente de estator x Velocidade','FontSize',13);

 

figure;

plot(n,Fp);

grid on;

xlabel('Velocidade [rpm]');

ylabel('FP');

title('Fator de potencia x Velocidade','FontSize',13);

 

figure;

plot(n,Rend);

xlabel('Velocidade [rpm]');

ylabel('Rendimento');

axis([0 1800 0 1.2]);

title('Rendimento x Velocidade','FontSize',13);

grid on;

 

figure;

plot(n,modIr);

grid on;

xlabel('Velocidade [rpm]');

ylabel('Ir [A]');

title('Corrente de rotor x Velocidade','FontSize',13);

 

figure;

plot(n,modImag);

grid on;

xlabel('Velocidade [rpm]');

ylabel('Imag [A]');

title('Corrente de magnetização x Velocidade','FontSize',13);

 

figure;

plot(n,Vth);

grid on;

xlabel('Velocidade [rpm]');

ylabel('Vth [V]');

title('Tensão de Thevenin x Velocidade','FontSize',13);