%**********************************************************
%
% Synchronous Generator
%
% ECE 520 Torque Angles for a Round Rotor S/G
%
% April 29, 2008
%
%**********************************************************

clc;
clear;

%**********************************************************
% I. Machine Parameters

H = 2.5;
w_eB = 2*60*pi;
T_max = 2.1;
T_m = 0.7;

%---------------------------------------------------------------------
% Calculate Torque Angles

delta_0 = asin(T_m/T_max);

delta_cr = acos((T_m/T_max)*(pi-2*delta_0) - cos(delta_0)); 

t_cr = sqrt(4*H*(delta_cr-delta_0)/(w_eB*T_m));

%---------------------------------------------------------------------
% Output Results

OPF = fopen('RoundRotorAngles.dat','w');

fprintf(OPF,'\n\n delta_0 = %5.4f radians \n', delta_0);
fprintf(OPF,' delta_0 = %5.4f degrees \n\n', delta_0*180/pi);

fprintf(OPF,'\n\n delta_cr = %5.4f radians \n', delta_cr);
fprintf(OPF,' delta_cr = %5.4f degrees \n\n', delta_cr*180/pi);

fprintf(OPF,'\n\n t_cr = %5.4f radians \n', t_cr);


fclose(OPF);