Scribd
Upload
91 views

Q.No.2 Using MATLAB, Calculate and Plot The Torque-Speed Characteristics of The Above Motor Both With The Original Rotor Resistance and With The Rotor Resistance Doubled

This document contains MATLAB code to calculate and plot the torque-speed characteristics of an induction motor under different operating conditions: 1) Original and doubled rotor resistance. 2) Supply frequencies of 60 Hz and 40 Hz. 3) Line voltages of 460V and 360V. In each case, the code calculates the torque-speed curve and plots the results for each condition on the same graph for comparison. Key motor parameters and operating points like maximum torque are also calculated.

Uploaded by

Khalid Akram
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
91 views

Q.No.2 Using MATLAB, Calculate and Plot The Torque-Speed Characteristics of The Above Motor Both With The Original Rotor Resistance and With The Rotor Resistance Doubled

This document contains MATLAB code to calculate and plot the torque-speed characteristics of an induction motor under different operating conditions: 1) Original and doubled rotor resistance. 2) Supply frequencies of 60 Hz and 40 Hz. 3) Line voltages of 460V and 360V. In each case, the code calculates the torque-speed curve and plots the results for each condition on the same graph for comparison. Key motor parameters and operating points like maximum torque are also calculated.

Uploaded by

Khalid Akram
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 6

Q.No.

2
(a) Using MATLAB, calculate and plot the torque-speed characteristics of the
above motor both with the original rotor resistance and with the rotor
resistance doubled.
clear all;
clc;
close all;

r1 = 0.641; % Stato r res i s tance


x1 = 1.106; % Stato r r eactance
r2 = 0.332; % Rotor r es i s tance
x2 = 0.464;
xm = 26.3;
v_phase = 460 / sqrt (3);
n_sync = 1800;
w_sync = n_sync*2*pi/60;

v_th = v_phase * (xm/(sqrt(r1^2 + (x1+xm)^2)));


v_th2 = 1.5*v_phase * (xm/(sqrt(r1^2 + (x1+xm)^2)));
z_th = ((1i*xm)*(r1+1i*x1))/(r1 + 1i*(x1 + xm));
r_th = real(z_th);
x_th = imag(z_th);
s = (0:1:50) / 50;
s(1) = 0.001;
nm = (1 - s) * n_sync;

for ii = 1:51
t_indl(ii) = (3 * v_th^2 * r2 / s (ii)) / (w_sync * (( r_th + r2 / s(ii))^2 + (x_th + x2)^2));
end

for ii = 1:51
t_ind2(ii) = (3 * v_th^2 * (2*r2) / s (ii)) / (w_sync * (( r_th + 2*r2 / s(ii))^2 + (x_th + x2)^2));
end

plot (nm , t_indl, 'Color' , 'k' , 'LineWidth',2.0);


hold on;
plot (nm , t_ind2, 'Color', 'k' , 'LineWidth' ,2.0, 'LineStyle','-.');
xlabel('\itn_m','Fontweight','Bold ');
ylabel('\tau_(ind)','Fontweight','Bold') ;
title('Induction motor torque-speed characteristic' , ...
'Fontweight','Bold');
legend('R_r = R_2 \Omega','R_r = 2 \times R_2 \Omega');
grid on;
hold off ;
T_max = 3*v_th^2/(2*w_sync*(r_th+sqrt(r_th^2+(x_th+x2)^2)))
T_start = (3 * v_th^2 * (1*r2)) / (w_sync * (((r_th + 1*r2)^2 + (x_th + x2)^2)))

(b) Using MATLAB, calculate and plot the torque-speed characteristics of the
above motor both with a supply frequency of 60 Hz and 40 Hz.
clear all;
clc;
close all;

r1 = 0.641; % Stato r res i s tance


x1 = 1.106; % Stato r r eactance
r2 = 0.332; % Rotor r es i s tance
x2 = 0.464;
xm = 26.3;
P = 4;
f1 = 60;
f2 = 40;
v_phase = 460 / sqrt (3);n_sync1 = 120*f1/P;
n_sync2 = 120*f2/4;
w_sync1 = n_sync1*2*pi/60;
w_sync2 = n_sync2*2*pi/60;
n_sync = 120;
v_th = v_phase * (xm/(sqrt(r1^2 + (x1+xm)^2)));
v_th2 = 1.5*v_phase * (xm/(sqrt(r1^2 + (x1+xm)^2)));
z_th = ((1i*xm)*(r1+1i*x1))/(r1 + 1i*(x1 + xm));
r_th = real(z_th);
x_th = imag(z_th);
s = (0:1:50) / 50;
s(1) = 0.001;
nm1 = (1 - s) * n_sync1;
nm2 = (1 - s) * n_sync2;
for ii = 1:51
t_indl(ii) = (3 * v_th^2 * r2 / s (ii)) / (w_sync1 * (( r_th + r2 / s(ii))^2 + (x_th + x2)^2));
end

for ii = 1:51
t_ind2(ii) = (3 * v_th^2 * r2 / s (ii)) / (w_sync2 * (( r_th + r2 / s(ii))^2 + (x_th + x2)^2));
end

plot (nm1 , t_indl, 'Color' , 'k' , 'LineWidth',2.0);


hold on;
plot (nm2 , t_ind2, 'Color', 'k' , 'LineWidth' ,2.0, 'LineStyle','-.');
xlabel('\itn_m','Fontweight','Bold ');
ylabel('\tau_(ind)','Fontweight','Bold') ;
title('Induction motor torque-speed characteristic' , ...
'Fontweight','Bold');
legend('f = 60 Hz','f = 40 Hz');
grid on;
hold off ;

T_max1 = 3*v_th^2/(2*w_sync1*(r_th+sqrt(r_th^2+(x_th+x2)^2)))
T_max2 = 3*v_th^2/(2*w_sync2*(r_th+sqrt(r_th^2+(x_th+x2)^2)))
T_start1 = (3 * v_th^2 * (1*r2)) / (w_sync1 * (((r_th + 1*r2)^2 + (x_th + x2)^2)))
T_start2 = (3 * v_th^2 * (1*r2)) / (w_sync2 * (((r_th + 1*r2)^2 + (x_th + x2)^2)))

(c) Using MATLAB, calculate and plot the torque-speed characteristics of the
above motor both with a line voltage of 460 V and 360 V.
clear all;
clc;
close all;

r1 = 0.641; % Stato r res i s tance


x1 = 1.106; % Stato r r eactance
r2 = 0.332; % Rotor r es i s tance
x2 = 0.464;
xm = 26.3;
v_phase1 = 460 / sqrt (3);
v_phase2 = 360 / sqrt (3);
n_sync = 1800;
n_sync2 = 120*40/4;
w_sync = n_sync*2*pi/60;

v_th1 = v_phase1 * (xm/(sqrt(r1^2 + (x1+xm)^2)));


v_th2 = v_phase2 * (xm/(sqrt(r1^2 + (x1+xm)^2)));
z_th = ((1i*xm)*(r1+1i*x1))/(r1 + 1i*(x1 + xm));
r_th = real(z_th);
x_th = imag(z_th);
s = (0:1:50) / 50;
s(1) = 0.001;
nm = (1 - s) * n_sync;

for ii = 1:51
t_indl(ii) = (3 * v_th1^2 * r2 / s (ii)) / (w_sync * (( r_th + r2 / s(ii))^2 + (x_th + x2)^2));
end

for ii = 1:51
t_ind2(ii) = (3 * v_th2^2 * r2 / s (ii)) / (w_sync * (( r_th + r2 / s(ii))^2 + (x_th + x2)^2));
end

plot (nm , t_indl, 'Color' , 'k' , 'LineWidth',2.0);


hold on;
plot (nm , t_ind2, 'Color', 'k' , 'LineWidth' ,2.0, 'LineStyle','-.');
xlabel('\itn_m','Fontweight','Bold ');
ylabel('\tau_(ind)','Fontweight','Bold') ;
title('Induction motor torque-speed characteristic' , ...
'Fontweight','Bold');
legend('V_{line} = 460 V','V_{line} = 360 V');
grid on;
hold off ;

T_max1 = 3*v_th1^2/(2*w_sync*(r_th+sqrt(r_th^2+(x_th+x2)^2)))
T_max2 = 3*v_th2^2/(2*w_sync*(r_th+sqrt(r_th^2+(x_th+x2)^2)))
T_start1 = (3 * v_th1^2 * (1*r2)) / (w_sync * (((r_th + 1*r2)^2 + (x_th + x2)^2)))
T_start2 = (3 * v_th2^2 * (1*r2)) / (w_sync * (((r_th + 1*r2)^2 + (x_th + x2)^2)))

You might also like