EPMN304
Electrical Machines 2
Computer Assignment 3
Submitted by :-
Ahmed Ayman - 1200782
1|Page
MATLAB CODE
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Givens %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
P=1000*10^3; % watt
Vline=4160; % volt
F=60; % Hz
n=865; % rpm
poles=8; % 2p=8
Prot=5500; % watt
R1=0.521; R2=1.32; X1=4.98; X2=5.32; Xm=136; %ohm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vphase=Vline/sqrt(3);
ns=(60*F)/(0.5*poles);
sfl=(ns-n)/ns;
ws=(2*pi*ns)/60;
w=(2*pi*n)/60;
V1eq=abs((Vphase*(Xm*1i))/(Xm*1i+X1*1i+R1));
Z1eq=(Xm*1i)*(R1+X1*1i)/(Xm*1i+X1*1i+R1);
R1eq=real(Z1eq);
X1eq=imag(Z1eq);
sm=R2/sqrt((R1eq)^2 + (X1eq+X2)^2);
sst=1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Full-Load Torque %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Tfl=(3*(V1eq)^2*R2)/(ws*sfl*((R1eq+R2/sfl)^2+(X1eq+X2)^2));
fprintf('Full Load Torque = %4.2f',Tfl); fprintf(' N-m \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Full-Load Current %%%%%%%%%%%%%%%%%%%%%%%%%%%%
Ifl=abs(V1eq/((R1eq+R2/sfl)+1i*(X1eq+X2)));
fprintf('Full Load Current = %4.2f',Ifl); fprintf(' A \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Break-Down Torque %%%%%%%%%%%%%%%%%%%%%%%%%%%%
Tmax=(3*(V1eq)^2*R2)/(ws*sm*((R1eq+R2/sm)^2+(X1eq+X2)^2));
fprintf('Break Down Torque = %4.2f',Tmax); fprintf(' N-m \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%% Locked-Rotor Torque %%%%%%%%%%%%%%%%%%%%%%%%%%%%
Tst=(3*(V1eq)^2*R2)/(ws*sst*((R1eq+R2/sst)^2+(X1eq+X2)^2));
fprintf('Locked-Rotor Torque = %4.2f',Tst); fprintf(' N-m \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%% Locked-Rotor Current %%%%%%%%%%%%%%%%%%%%%%%%%%%
Ist=abs(V1eq/((R1eq+R2/sst)+1i*(X1eq+X2)));
fprintf('Locked-Rotor Current = %4.2f',Ist); fprintf(' A \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Plotting Variables %%%%%%%%%%%%%%%%%%%%%%%%%%%
s=sfl:-0.01:0;
n_=1:length(s);
z=1:length(s);
I=1:length(s);
pf=1:length(s);
Pin=1:length(s);
I2=1:length(s);
Pcu2=1:length(s);
Pd=1:length(s);
2|Page
Po=1:length(s);
Eff=1:length(s);
Imagnitude=1:length(s);
I2magnitude=1:length(s);
Td=1:length(s);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:length(s)
n_(i)=(1-s(i))*ns;
z(i)=R1+(1i*X1)+((((R2/s(i))+(1i*X2))*(1i*Xm))/((R2/s(i))+(1i*X2)+(1i*Xm)));
I(i)=Vphase/z(i);
pf(i)=cos(0-atan(imag(I(i))/real(I(i))));
Imagnitude(i)=(sqrt(((real(I(i)))^2)+((imag(I(i)))^2)));
Pin(i)=3*Vphase*Imagnitude(i)*pf(i);
I2(i)=I(i)*((1i*Xm)/((1i*Xm)+(R2/s(i))+(1i*X2)));
I2magnitude(i)=(sqrt((real(I2(i))^2+(imag(I2(i)))^2)));
Pcu2(i)=3*((I2magnitude(i))^2)*R2;
Pd(i)=Pcu2(i)*((1-s(i))/s(i));
Po(i)=Pd(i)-Prot;
Eff(i)=((Po(i))/(Pin(i)))*100;
Td(i)=Pd(i)/(2*pi*n_(i)/60);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Plotting %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%% Output Power & Supply Current %%%%%%%%%%%%%%%%%%%%%%%
figure
plot(Po,Imagnitude, 'b-', 'LineWidth', 1);
xlabel('Output Power in Watt'),ylabel('Supply Current');
grid on;
%%%%%%%%%%%%%%%%%%%% Output Power & Input Power %%%%%%%%%%%%%%%%%%%%%%%%%%
figure
plot(Po,Pin, 'b-', 'LineWidth', 1);
xlabel('Output Power in Watt'),ylabel('Input Power in Watt');
grid on;
%%%%%%%%%%%%%%%%%% Output Power & Developed Torque %%%%%%%%%%%%%%%%%%%%%%%
figure
plot(Po,Td, 'b-', 'LineWidth', 1);
xlabel('Output Power in Watt'),ylabel('Developed Torque');
grid on;
%%%%%%%%%%%%%%%%%%%% Output Power & Power Losses %%%%%%%%%%%%%%%%%%%%%%%%%
figure
plot(Po,Pcu2+Prot, 'b-', 'LineWidth', 1);
xlabel('Output Power in Watt'),ylabel('Power Losses in Watt');
grid on;
%%%%%%%%%%%%%%%%%% Output Power & Input Power Factor %%%%%%%%%%%%%%%%%%%%%
figure
plot(Po,pf, 'b-', 'LineWidth', 1);
xlabel('Output Power in Watt'),ylabel('Power Factor');
grid on;
%%%%%%%%%%%%%%%%%%%%% Output Power & Effeciency %%%%%%%%%%%%%%%%%%%%%%%%%%
figure
plot(Po,Eff,'b-', 'LineWidth', 1);
xlabel('Output Power in Watt'),ylabel('Efficiency %');
grid on;
3|Page
Matlab Output
These rules were used to calculate the torque and current.
The full load torque and current were calculated at S full load.
The break-down torque, which is the maximum torque at STm, was calculated.
The locked rotor torque, which is the starting torque at S=1, was also calculated
4|Page
Graphs Output
Supply Current
As the supply current increases, so does the output power, as they are directly proportional.
Input Power
As the input power increases, so does the output power, as they are directly proportional.
4|Page
Power Factor
The input power factor rises until it reaches the maximum output power under rated conditions.
Efficiency
Efficiency rises with output power until it reaches the maximum power opposite to Smax, beyond
which it decreases.
5|Page
