Page 224 - Solutions Manual to accompany Electric Machinery Fundamentals

n_f = 2700; % Number of turns on field
f_ar0 = 1000; % Armature reaction @ 110 A (A-t/m)

% Calculate the armature current for each load.
i_a = i_l - v_t / (r_f + r_adj);

% Now calculate the internal generated voltage for
% each armature current.
e_a = v_t - i_a * r_a;

% Calculate the armature reaction MMF for each armature
% current.
f_ar = (i_a / 55) * f_ar0;

% Calculate the effective field current with and without
% armature reaction. Ther term i_f_ar is the field current
% with armature reaction, and the term i_f_noar is the
% field current without armature reaction.
i_f_ar = v_t / (r_f + r_adj) - f_ar / n_f;
i_f_noar = v_t / (r_f + r_adj);

% Calculate the resulting internal generated voltage at
% 1200 r/min by interpolating the motor's magnetization
% curve.
e_a0_ar = interp1(if_values,ea_values,i_f_ar);
e_a0_noar = interp1(if_values,ea_values,i_f_noar);

% Calculate the resulting speed from Equation (9-13).
n_ar = ( e_a ./ e_a0_ar ) * n_0;
n_noar = ( e_a ./ e_a0_noar ) * n_0;

% Calculate the induced torque corresponding to each
% speed from Equations (8-55) and (8-56).
t_ind_ar = e_a .* i_a ./ (n_ar * 2 * pi / 60);
t_ind_noar = e_a .* i_a ./ (n_noar * 2 * pi / 60);

% Plot the torque-speed curves
figure(1);
plot(t_ind_noar,n_noar,'b-','LineWidth',2.0);
hold on;
plot(t_ind_ar,n_ar,'k--','LineWidth',2.0);
xlabel('\bf\tau_{ind} (N-m)');
ylabel('\bf\itn_{m} \rm\bf(r/min)');
title ('\bfShunt DC Motor Torque-Speed Characteristic');
legend('No armature reaction','With armature reaction');
axis([ 0 125 800 1600]);
grid on;
hold off;

218