Ind Motors Assignment

Download as pdf or txt
Download as pdf or txt
You are on page 1of 7

Assignment – Single Phase

Induction Motor

Erandi M.A.U
200161V
With Starting Winding
%% With Starting Winding

% Motor parameters
R1 = 1.8;
X1 = 2.4;
R2 = 2.5;
X2 = 2.4;
R1_Aux = 10;
X1_Aux = 0.495;
Xm = 60;
V_m = 120;
V_a = 120;
N_s = 1800;

%% Turns ratio
a = sqrt(X1_Aux / X1);

%%Synchronous speed in rpm


W_s = (N_s * 2 * pi) / 60;

%% Calulating Impedance
Z1_aux = R1_Aux + (1i * X1_Aux);
Z1_motor = R1 + (1i * X1);

%% Variables for speed and torque


Speed = [];
Torque = [];

%% Motor speed with starting winding


for Slip_1 = 1:-0.002:0
N = N_s - (Slip_1 * N_s);
W = (N * 2 * pi) / 60;
Z_bwd = (((R2 / (2 - Slip_1)) + (1i * X2)) * (1i * Xm))
/ ((R2 / (2 - Slip_1)) + (1i * X2) + (1i * Xm));
Z_fwd = (((R2 / Slip_1) + (1i * X2)) * (1i * Xm)) /
((R2 / Slip_1) + (1i * X2) + (1i * Xm));

I_aux = (V_a*(Z1_motor + Z_fwd + Z_bwd) -


((V_m*a*(Z_fwd-Z_bwd))*1i))/((Z1_motor + Z_fwd+
Z_bwd)*(Z1_aux + (a^2)*(Z_fwd + Z_bwd)) - ((a^2)*(Z_fwd -
Z_bwd)^2));
I_motor = (V_m + (I_aux*a*(Z_fwd -
Z_bwd))*1i)/(Z1_motor + Z_fwd + Z_bwd);
%Power
AG_Power = ((abs(I_motor)^2 + abs(I_aux * a)^2) *
(real(Z_fwd) - real(Z_bwd))) + 2 * a * abs(I_aux) *
abs(I_motor) * real(Z_fwd + Z_bwd) * sin(angle(I_aux) -
angle(I_motor));

%Torque
T = AG_Power / W_s;

Speed(end+1) = N;
Torque(end+1) = T;

if N == 0.7 * N_s
break
end
end

% Motor speed after disconnecting the starting winding


for S2 = Slip_1:-0.001:0
N = N_s - (S2 * N_s);
W = (N * 2 * pi) / 60;

% Calculate impedances
Z_fwd = ((R2 / S2) + (1i * X2)) * (1i * Xm) / ((R2 /
S2) + (1i * X2) + (1i * Xm));
Z_bwd = ((R2 / (2 - S2)) + (1i * X2)) * (1i * Xm) /
((R2 / (2 - S2)) + (1i * X2) + (1i * Xm));

% Calculate current and power


I_1 = V_m / (R1 + (1i * X1) + 0.5 * Z_fwd + 0.5 *
Z_bwd);
AG_PowerF = abs(I_1)^2 * (0.5 * real(Z_fwd));
AG_PowerB = abs(I_1)^2 * (0.5 * real(Z_bwd));
AG_Power = AG_PowerF - AG_PowerB;

% Calculate torque
T = AG_Power / W_s;

Speed(end+1) = N;
Torque(end+1) = T;
end
%% Torque-speed characteristics
figure
plot(Speed, Torque, 'LineWidth', 2)
title('Torque-Speed Characteristic');
xlabel('Rotor speed (rpm)')
xlim([0 2000])
ylabel('Torque (Nm)')
ylim([0 5])
grid on
Without Starting Winding
%% Without starting winding %%%%%%%%%%%%%%%%%%%%%

%Calculating slip values for 0 to 1800


N_m = linspace(0, N_s-1, N_s);
Slip_2 =(N_s - N_m)/ N_s;

%% forword and backword Impedances


Z_fwd2 = ((R2./Slip_2 + 1i*X2)*(1i*Xm))./(R2./Slip_2 +
1i*X2 + 1i*Xm);
Z_bwd2 = ((R2./(2-Slip_2) + 1i*X2)*(1i*Xm))./(R2./(2-
Slip_2) + 1i*X2 + 1i*Xm);

%% Input current
I_1 = V_m./(R1 + 1i*X1 + 0.5*Z_fwd2 + 0.5*Z_bwd2);

%% calculating the air-gap power


AG_PowerF = (abs(I_1).^2)*0.5.*real(Z_fwd2);
AG_PowerB = (abs(I_1).^2)*0.5.*real(Z_bwd2);
AG_Power = AG_PowerF - AG_PowerB;

%% Calculate Torque
T_2 = AG_Power./(N_s*((2*pi)/60));

%% Torque speed characteristic


figure;
plot(N_m, T_2,'LineWidth', 2);
title('Torque-Speed Characteristic');
xlabel('Speed (rpm)');
ylabel('Induced Torque (N-m)');
grid on;

You might also like