Universidad Nacional de San Antonio Abad del

Cusco
Gráfico de Funciones
Máquinas Eléctricas - Grupo A

Edison Abado Ancco 145012

Gabriela Morales Soto 145022
Christian Rimachi Maccarcco 161771
Brayahan Espinoza Salas 160897
Watner Ochoa Núñez 171174
Javier Pauccar Queshuallpa 164255
Winer Andree Rios Rado 150427

8 de enero de 2020

1. Lugares Geométricos de la Máquina Síncrona

Las posibles condiciones de operación de la máquina síncrona circunscriben las diferentes
variables que definen su comportamiento, en ciertos rangos y figuras geométricas cuando
se representan gráficamente. Un determinado punto de operación queda definido mediante
un diagrama fasorial, pero la variación del factor de potencia, genera un lugar geométrico
para la fuerza electromotriz producida por el campo. El análisis de estos diagramas permite
evaluar las diferentes posibilidades y limitaciones en la operación de la máquina síncrona.
Algunos diagramas o lugares geométricos son suministrados por el fabricante de la máquina
debido a su utilidad e importancia. En otras ocasiones es necesario construirlos a partir de
la información disponible.
La figura siguiente presenta el lugar geométrico de la fuerza electromotriz de una máqui-
na síncrona con polos salientes (con sus datos respectivos), el círculo correspondiente a la
potencia aparente nominal y el círculo correspondiente a la fuerza electromotriz máxima.
El gráfico de la figura anterior (con sus respectivos datos) se obtiene mediante el programa
MATLAB que se lista a continuación:

1 %LUGARES GEOMETRICOS DE LA MAQUINA SINCRONA DE POLOS SALIENTES

2 Van = 1. ; %Tension nominal
3 Ian = 1. ; %Corriente nominal
4 San = 1. ; %Potencia aparente nominal
5 CosFin = 0.8 ; %Factor de potencia nominal
6 Xd = 1.2; %Reactancia del eje directo
7 Xq = 0.8; %Reactancia del eje de cuadratura
8 Re = 0.0; %Resistencia de las bobinas de armadura
9 Fi = 0: -0.01: -2* pi; %Vector de posibles fdp de operacion
10 Ia = Ian*( cos(Fi) + j*sin(Fi)); %Vector de corriente de armadura
11 D = Van+ (Re + j*Xq).*Ia; %Vector de fasores D para c/corriente
12 delta = atan2(imag(D), real(D)); %Vector de angulos de carga para cada corriente ←-
de armadura
13 Id = Ian*sin(delta -Fi).*( sin(delta) - j*cos(delta)); %Corriente del eje directo
14 Ef = D + j * (Xd - Xq).* Id; %Calculo de la fuerza electomotriz
15 Fin = acos(CosFin); %Angulo de factor de potencia nominal
16 Iax = Ian*( sin(cos(Fin))-j*sin(Fin)); %Fasor de la corriente nominal
17 Dx = Van + (Re + j*Xq).*Iax; %Fasor D para la condicion nominal

1
 Ingenieria Electrónica UNSAAC 2019-II

Figura 1: Lugares geométricos de una máquina síncrona de polos salientes

18 deltax = atan(imag(Dx)/real(Dx)); %Angulo de carga nominal

19 Idx = Ian*sin(deltax - Fin).*( sin(deltax) - j*cos(deltax)); %Corriente direcc. nominal
20 Efx = Dx + j*(Xd - Xq).*Idx; %Fasor de la fem maxima
21 Efmax = abs(Efx) * (cos(Fi) + j*sin(Fi)); %Magnitud de la fem maxima
22
23 axis('square '); %Def. de ejes cuadrados (grafica)
24 axis ([ -2.5 2.5 -2.5 2.5]); %Escalas del grafico en p.u.
25 xlabel('Eje Real ');
26 ylabel('Eje Imaginario ');

El siguiente código muestra la gráfica progresiva de cada punto para una variación de Fi,
y el gráfico se muestra en la figura (2)

1
2 %% % % % % % % % % % % % % % % % %--GRAFICA
% PROGRESIVA
- %-% % % % % % % % % % % % % % % % % % % % % % % % % % %
3
4 for k=1: length(Id)
5 Ef=D(k) + j* (Xd -Xq).*Id(k);
6 Ia=Ian*( cos(Fi(k)) + j*sin(Fi(k)));
7 hold on;
8 grid on;
9 Ejexx = [-2.5 2.5];
10 Ejexy = [0 0];
11 Ejeyx = [0 0];
12 Ejeyy = [2.5 -2.5];
13 plot (Ef , 'b.')
14 plot (Ia , 'r.');
15 plot (Efmax , 'g.');
16 plot (Ejexx , Ejexy , 'k', Ejeyx , Ejeyy , 'k');
17 pause (0.01); %insertar un delay
18 end
19
20 pause (1);
21 cla

Y para teminar se tiene el siguiente código que demuestra a la gráfica de la figura (3),
una gráfica estática con todos los datos necesarios.

1 hold on;
2 grid on;
3 Ef = D + j*(Xd -Xq).*Id; %Calculo de la fem
4 plot(Ef ,'b'); %Grafica del lugar geometrico de la fem a ←-
potencia aparente constante

1 LUGARES GEOMÉTRICOS DE LA MÁQUINA SÍNCRONA 2

 Ingenieria Electrónica UNSAAC 2019-II

Figura 2: Gráfica progresiva de las variables. Azul: FEM a potencia aparente constante, Rojo:
Potencia aparente onstante, Verde: FEM Máxima.

5 Ia = Ian*( cos(Fi) + j*sin(Fi)); %Vector de corrientes de armadura

6 plot(Ia , 'r'); %Grafica del lugar geometrico de la Sn
7 Efmax = abs(Efx)*( cos(Fi)+j*sin(Fi)); %Magnitud de la fem maxima
8 plot(Efmax , 'g'); %Grafica del lugar geometrico de la fem ←-
maxima
9 xlabel('Eje Real ');
10 ylabel('Eje Imaginario ');
11 Ejexx = [-2.5 2.5];
12 Ejexy = [0 0];
13 Ejeyx = [0 0];
14 Ejeyy = [2.5 -2.5];
15 plot(Ejexx , Ejexy , 'k', Ejeyx , Ejeyy , 'k'); %GRafica de los ejes
16 legend('FEM a Sn cte', 'Potencia aparente nominal ', 'FEM max.');
17 legend('Location ', 'southwest ')

Figura 3: Gráfica final de las variables rotuladas con su respectiva leyenda.

2. Curvas en V de la Máquina Síncrona

Estas curvas están formadas por una familia de gráficos que representan la relación entre
la magnitud de la corriente del estator Ia, y la corriente de campo o corriente de excitación If,

2 CURVAS EN V DE LA MÁQUINA SÍNCRONA 3

 Ingenieria Electrónica UNSAAC 2019-II

utilizando como parámetro diferentes valores de la potencia eléctrica. Estos gráficos se realizan
manteniendo la tensión de armadura en un valor constante, generalmente en su valor nominal.
Mediante el diagrama fasorial de la máquina sincrónica se pueden obtener directamente las
curvas en V de la máquina sincrónica. En la figura siguiente se han representado una familia
de curvas en “V” para una máquina sincrónica de polos salientes. Es interesante destacar la
linealidad de la curva cuando la potencia eléctrica es cero. En este caso todas las caídas de
tensión y fuerzas electromotrices coinciden exactamente con el eje cuadratura, obteniéndose
una relación entre la fuerza electromotriz y la corriente de armadura es lineal. La zona a la
derecha del gráfico corresponde a inyección de reactivos desde la máquina hacia la red y a
la izquierda de la característica de factor de potencia unitario, se consume potencia reactiva
inductiva desde la red eléctrica.

Figura 4: Curvas en V para una máquina sincrónica de polos salientes.

Las curvas en V fueron utilizadas en el pasado con la finalidad de evitar los laboriosos
cálculos fasoriales. Estas curvas permiten una rápida visualización de los límites operativos
de la máquina. Sin embargo, actualmente no es necesario el uso de esta herramienta debido
a las facilidades de cálculo disponibles. A continuación se lista un programa en el entorno
MATLAB que reproduce el gráfico de la figura (4). El algoritmo utiliza la rutina fsolve
incorporada en el entorno de programación, para determinar el valor del ángulo de carga en
cada punto de operación. El listado de este programa reproduce algunos aspectos prácticos
del cálculo de máquinas sincrónicas en régimen permanente:

1 %CURVAS EN 'V' DE LA MAQUINA SINCRONA DE POLOS SALIENTES

2 global c1 Ef Van xd xq Pe;
3 Van = 1. ; %Tension nominal
4 Ian = 1. ; %Corriente nominal
5 San = 1. ; %Potencia aparente nominal
6 CosFin = 0.8; %Factor de potencia nominal
7 xd = 0.9; %Reactancia del eje directo
8 xq = 0.6; %Reactancia del eje de cuadratura
9 Re = 0.0; %Resis. de las bobinas de armadura
10
11 Fin = acos(CosFin); %Angulo de fdp nominal
12 Iax = Ian*( cos(Fin)-j*sin(Fin)); %Fasor de la corriente nominal
13 Dx = Van + (Re + j*xq)*Iax; %Fasor D para la condicion nominal
14 deltax = atan(imag(Dx)/real(Dx)); %Angulo de carga nominal
15 Idx = Ian*sin(deltax -Fin).*( sin(deltax)-j*cos(deltax));
16 %Corriente directa nominal
17
18 Efx = Dx + j*(xd -xq).*Idx; %Fasor de la FEM maxima
19 Efmax = abs(Efx)*( cos(Fin) + j*sin(Fin)); %Magnitud de la FEM maxima
20 Ef = 0.1 : 0.01 : 2.0; %Rango de fuerzas electromotrices
21 kk = 0;
22 for Pe = 0: 0.2: 1.2 %Rango de potencias electricas
23 kk = kk + 1 ;

2 CURVAS EN V DE LA MÁQUINA SÍNCRONA 4

 Ingenieria Electrónica UNSAAC 2019-II

24 for l = 1 : length(Ef)
25 delta(l) = 0; %Inicializacion del Angulo de Carga
26 end
27 for c1 = 1 : 1 : length(Ef)
28 aaa = Pe*xd/(Ef(c1)*Van);
29 if abs(aaa) < 1.0
30 deltac1 = asin(aaa); %Estimacion inicial del Ang. de carga
31 deltac1 = fsolve(@(delta) Ef(c1)*Van*sin(delta)/xd+Van ^2*(xd -xq)*sin (2* delta)←-
/(2* xd*xq)-Pe , deltac1);
32 %Solucion de la ecuacion no lineal
33 Q = Ef(c1)*Van.* cos(deltac1)/xd -Van ^2*( xq*cos(deltac1).^2+ xd*sin(deltac1)←-
.^2/( xd*xq));
34 %Calculo de la potencia reactiva
35 S = sqrt(Pe*Pe+Q.*Q); %Calculo de la pontecia aparente
36 Ia(c1 , kk) = S/Van; %Calculo de la corriente de armadura
37 else
38 Ia(c1 , kk) = 0; %Corriente de armadura en casos no factibles
39 end
40 end
41 end
42 %Grafico de las curvas en 'V'
43 axis on;
44 plot(Ef ,Ia(:,1),'.',Ef ,Ia(:,2),'.',Ef , Ia(:,3),'.',Ef ,Ia(:,4),'.',Ef ,Ia(:,5),'.',Ef ,Ia←-
(:,6),'.',Ef ,Ia(:,7),'.');
45 hold on;
46 plot(abs(Efmax), 0.0, 'x');
47 hold off;
48 axis ([0 2 0 1.5]); %Definicion delos ejes del grafico
49 grid on;
50 legend('P=0','P=0.2 ','P=0.4 ','P=0.6 ','P=0.8 ','P=1.0 ','P=1.2 ');
51 legend('Location ', 'southeast ');
52 xlabel('Ef , If');
53 ylabel('|Ia|');

Las gráficas generadas por el código de MATLAB, se muestran en la figura (5).

Figura 5: Curvas en V de la máquina síncrona.

3. Análisis transitorio de la Máquina Síncrona

A continuación se reproduce el listado y los resultados de un algoritmo en el entorno
de programación MATLAB que permite evaluar las corrientes en coordenadas primitivas y
transformadas durante un cortocircuito brusco de la máquina sincrónica, considerando las
atenuaciones debidas a las resistencias de las bobinas. Se tiene como gráficas las figuras (6)
y (7).

1 w = 1.0; %Velocidad sincronica

2 re = 0.01; %Resistencia del estator

3 ANÁLISIS TRANSITORIO DE LA MÁQUINA SÍNCRONA 5

 Ingenieria Electrónica UNSAAC 2019-II

3 rf = 0.02; %Resistencia del campo

4 R = diag ([re re rf], 0); %matriz de resistencias
5 Ld = 1.0; %inductancia del eje d
6 Lq = 0.6; %inductancia del eje q
7 Lf = 1.0; %inductancia del campo
8 Ldf = 0.8; %inductancia mutua df
9 L = [Ld 0 Ldf; 0 Lq 0; Ldf 0 Lf]; %matriz de inductancias
10 G = w*[0 -Lq 0; Ld 0 Ldf; 0 0 0]; %matriz de generaci 'on
11 Lin = inv(L); %inversa de L
12 A = -Lin*(R + G); % matriz de transici 'on de estado
13 Ain = inv(A); %Inversa de la matriz A
14 [V, g] = eig(A); %Autovalores y autovectores
15 io = [0; 0; 1]; %Condiciones iniciales de prefalla
16 t = 0: 0.01:160; %tiempo de la soluci 'on
17 e1 = exp(g(1,1)*t); %Decaimiento del 1er autovalor
18 e2 = exp(g(2,2)*t); %Decaimiento del 2do autovalor
19 e3 = exp(g(3,3)*t); %Decaimiento del 3er autovalor
20 ip = -Ain*Lin *[0; 0; rf]; %soluci 'on particular p[i] = 0
21 ih = io -ip; %corrientes homog 'eneas en t=0+
22 k = inv(V)*ih; %coeficientes homog 'eneos
23 ke = diag(k, 0)*[e1; e2; e3];
24 i = V*ke; %soluci 'on homog 'enea
25 n = length(t); %n'umero de puntos temporales
26 for m = 1:n %soluci 'on total
27 i(1,m) = i(1,m) + ip(1);
28 i(2,m) = i(2,m) + ip(2);
29 i(3,m) = i(3,m) + ip(3);
30 end
31 %==================================== %
32 %Gr 'afico de las corriente Id , Iq e If
33 %==================================== %
34 plot(t, i(1,:), t, i(2,:), t, i(3,:));
35 legend('Id', 'Iq', 'If');
36 xlabel('Tiempo ');
37 ylabel('Corriente ');
38 title('Corrientes de corto circuito en coordenadas dqf');
39 grid minor;
40 pause (20);
41 cla
42 %==================================== %
43 %Gr 'afico de las coordenadas If e Ia
44 %==================================== %
45 ia = (i(1,:).* cos(w*t) - i(2,:).* sin(w*t));
46 plot(t, ia , t, i(3,:));
47 legend('Ia', 'If');
48 xlabel('Tiempo ');
49 ylabel('Corriente ');
50 title('Corrientes de corto circuito brusco en coordenadas primitivas ia e if');
51 grid minor;

Figura 6: Gráfica de corrientes de corto circuito en coordenadas dqf (Id, Iq, If ) para una
máquina síncrona.

3 ANÁLISIS TRANSITORIO DE LA MÁQUINA SÍNCRONA 6

 Ingenieria Electrónica UNSAAC 2019-II

Figura 7: Gráfica de corrientes de corto circuito brusco en coordenadas primitivas Ia e If .

4. Análisis transitorio de Máquina Síncrona con devanados

amortiguadores
Se tiene el siguiente código en MatLab y sus gráficas en las figuras (8), (9) y (10).

1 %An 'alisis transitorio de m'aquinas s'incronas con devanados amortiguadores

2 w = 1;
3 rd = 0.01; rq = 0.01; rad = 0.03;
4 raq = 0.03; rf = 0.01;
5 R = diag ([rd rq rad raq rf], 0);
6
7 Ld = 1.0 ; Lq = 0.6 ; Lad = 0.9 ; Laq = 0.5;
8 Lf = 1; Ldad = 0.85 ; Lqaq = 0.4; Ldaq = 0.4;
9 Lqad = 0.8; Lfad = 0.7; Ldf = 0.8;
10 L = [Ld 0 Ldad 0 Ldf; 0 Lq 0 Lqaq 0; Ldad 0 Lad 0 Lfad; 0 Lqaq 0 Laq 0;
11 Ldf 0 Lfad 0 Lf];
12 G = w*[0 -Lq 0 -Ldaq 0; Ld 0 Lqad 0 Ldf; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0];
13 Lin = inv(L);
14 A = -Lin*(R+G); %Matriz de transici 'on de estado
15 Ain = inv(A);
16 [V, g] = eig(A); %autovalores y autovectores
17 io = [0; 0; 0; 0; 1]; %condiciones iniciales
18 t = 0: 0.25 : 160;
19 e1 = exp(g(1,1)*t); e2 = exp(g(2,2)*t);
20 e3 = exp(g(3,3)*t); e4 = exp(g(4,4)*t);
21 e5 = exp(g(5,5)*t);
22 ip = -Ain*Lin *[0; 0; 0; 0; rf *1]; %Soluci 'on particular p[i]=0
23 ih = io - ip;
24 k = inv(V)*ih;
25 ke = diag(k, 0)*[e1; e2; e3; e4; e5];
26 i = V*ke;
27 n = length(t);
28 for m = 1: n
29 i(1,m) = i(1,m) + ip(1);
30 i(2,m) = i(2,m) + ip(2);
31 i(3,m) = i(3,m) + ip(3);
32 i(4,m) = i(4,m) + ip(4);
33 i(5,m) = i(5,m) + ip(5);
34 end
35 figure (1); %Corrientes dqf durante el corto circuito brusco en devanados amortiguadores
36 plot(t, i(1,:), t, i(2,:), t, i(5,:));
37 legend('Id', 'Iq', 'If');
38 ylabel('Corriente ');
39 xlabel('Tiempo ');
40 title('Corrientes dqf durante el corto circuito brusco con devanados amortiguados ');
41 grid minor;
42
43

4 ANÁLISIS TRANSITORIO DE MÁQUINA SÍNCRONA CON DEVANADOS 7

AMORTIGUADORES
 Ingenieria Electrónica UNSAAC 2019-II

44 figure (2); %Corrientes de los devanados amortiguadores

45 plot(t, i(3,:), t, i(4,:));
46 legend('Iad', 'Iaq');
47 ylabel('Corriente ');
48 xlabel('Tiempo ');
49 title('Corrientes de los devanados amortiguadores ');
50 grid minor;
51
52
53 %Corriente en coordenadas primitivas
54 ia = (i(1,:).* cos(w*t) - i(2,:).* sin(w*t));
55 figure (3);
56 plot(t, ia ,t, i(5,:));
57 legend('Ia', 'If');
58 ylabel('Corriente ');
59 xlabel('Tiempo ');
60 title('Corriente en coordenadas primitivas ');
61 grid minor;

Figura 8: Gráfica de corrientes dqf durante el corto circuito brusco de la máquina síncrona
con devanados amortiguadores.

Figura 9: Gráfica de corrientes de los devanados amortiguadores durante el corto circuito

brusco de la máquina síncrona.

4 ANÁLISIS TRANSITORIO DE MÁQUINA SÍNCRONA CON DEVANADOS 8

AMORTIGUADORES
Ingenieria Electrónica UNSAAC 2019-II

Figura 10: Gráfica de Corrientes en coordenadas primitivas durante el corto circuito brusco
de la máquina síncrona (Ia, If ).

5. Problemas: Opcionales
A continuación se presentan dos problemas de Máquina Síncrona.

5.1. Problema A
Un motor síncrono de 208V , 45 hp, un factor de potencia de 0,8 en adelanto, conectado
en D, a 60Hz, alimenta a una carga de 15 hp con un factor de potencia inicial de 0,85 en
retraso. La corriente de campo IF en estas condiciones es de 4,0A.

a) Dibuje el diagrama fasorial inicial del motor y encuentre los valores de IA y EA .

b) Dibuje el nuevo diagrama fasorial del motor si se incrementara el flujo del motor en
25 %. ¿Cuál es el valor de IA y EA y el factor de potencia del motor en este momento?

c) Suponga que el flujo en el motor varía linealmente con la corriente de campo IF .

Haga una gráfica de IA contra IF del motor síncrono con una carga de 15 hp.

5.1.1. Solución a)
Del ejemplo anterior se sabe que la potencia eléctrica de entrada, incluyendo todas las
pérdidas, es de Pent = 13,69kW . Ya que el factor de potencia del motor es de 0,85 en retraso,
el flujo de corriente en el inducido resultante es

Pent
IA = (1)
3Vφ cosθ
13,69kW
= 25,8 (2)
3(208V )(0,85)
El ángulo θ es cos−1 (0,85) = 31,8, por lo que la corriete del fasor IA es igual a

208 0o V − (j2,5Ω)(25,8 −31,8o A) = 208 0o V − 64,5 58,2o V = 182 −17,5V (3)

El diagrama fasorial resultante se muestra en la figura

5 PROBLEMAS: OPCIONALES 9
Ingenieria Electrónica UNSAAC 2019-II

Figura 11: Diagrama fasorial del problema 1

5.1.2. Solución b)
Si el motor se encuentra trabajando en la parte no saturada de la máuina y se incrementa
25 % el flujo φ, entonces EA 5Kϕω también se incrementará 25 %

EA2 = 1,25EA1 = 1,25(128V ) = 227,5V (4)

Sin embaro, la pottencia suministrada a la carga debe permanecer constante. Debido a
que la distancia debe ser constante en el diagrama fasorial del nivel de flujo originalñ al nuevo
de flujo, por lo tanto:

EA1 senδ1 = EA1 senc (5)

 EA

δ2 = sen−1 senδ1 |EA1 = sen−1 sen(−17,5)|182V o
(6)


2
227,5V = −13,9

Ahora se puede calcular la corriente en el inducido a partir de la ley de voltaje de Kirchhoff:

Vϕ − EA2
IA2 = (7)
jXS
208 0o V − 227,5 0o − 13,9o V
IA = = 22,5 13,2o A (8)
j2,5Ω
Finalmete, el factor de potencia del motor es:

F P = cos(13,2) = 0,974 en adelanto (9)

El diagrama fasorial resultante se muestra en la figura (11)

5.1.3. Solución c)
Debido a que se supone que el flujo varía linealmente con la corriente de campo, EA
también variará linealmente con ella. Se sabe que EA es 182 V para una corriente de campo
de 4,0A, por lo que se puede calcular EA para cualquier corriente de campo con la relación

EA2 IF2
= → EA2 = 45,5IF2 (10)
182V 4,0A
El ángulo de par δ de cualquier corriente de campo dada se puede calcular a partir del
hecho de que la potencia suministrada a la carga debe permanecer constante:0

EA1 sen(dδ1 ) = EA2 sen(δ2 ) (11)

Por lo que
 EA

δ2 = sen−1 senδ1 |EA1 (12)
2

5 PROBLEMAS: OPCIONALES 10
 Ingenieria Electrónica UNSAAC 2019-II

Esos dos datos nos dan el voltaje fasorial EA . una vez que se conoce EA , se puede calcular
la nueva corriente en el inducido por medio de la ley de voltaje de kirchhoff

Vϕ − EA2
IA2 = (13)
jXS
A continuación se muestra un archivo .m de MATLAB para calcular y hacer la gráfica de
IA e IF por medio de las ecuaciones anteriores; y su gráfico se muestra en la figura (12)

1 % Primero , inicializar los valores de la corriente de campo (21

2 % valores en el rango de 3.8 -5.8 A)
3 i_f = (3.8:0.1:5.8);
4 % Ahora inicializar todos los dem 'as valores
5 i_a = zeros (1 ,21); % Crear espacio para la matriz i_a
6 x_s = 2.5; % Reactancia s'incrona
7 v_phase = 208; % Voltaje de fase a 0 grados
8 delta1 = -17.5 * pi /180; % delta 1 en radianes
9 e_a1 = 182*( cos(delta1)+1i*sin(delta1));
10 % Calcular la corriente en el inducido para cada valor
11 for iia = 1:21
12 % Calcular la magnitud de e_a2
13 e_a2 = 45.5 * i_f(iia);
14 % Calcular delta2
15 delta2 = asin(abs(e_a1)/abs(e_a2)*sin(delta1));
16 % Calcular el fasor e_a2
17 e_a2 = e_a2 *( cos(delta2)+1i*sin(delta2));
18 % Calcular i_a
19 i_a(iia) = (v_phase -e_a2)/(1i*x_s);
20 end
21 % Hacer la gr 'afica de la curva en V
22 plot(i_f ,abs (i_a),'Color ', 'k','Linewidth ' ,2.0)
23 xlabel('Corriente de campo (A)','Fontweight ','Bold ');
24 ylabel('Corriente en el inducido (A)','Fontweight ','Bold ');
25 title ('Curva en V del motor sincrono ','Fontweight ','Bold ');
26 grid on

Figura 12: Diagrama fasorial del problema 1

5.2. Problema B
Un generador síncrono de 75 M V A, 13,8KV , trifasico con una reactancia síncrona de
XS = 1,35 por unidad y una reactancia no saturada de XS,u = 1,56 por unidad, se conecta
a un sistema externo con reactancia equivalente a XEQ = 0,23 por uhnidad y un voltaje
VEQ = 1,0 por unidad, ambos sobre la base del generador. Esta máquina alcanza su voltaje
de circuito abierto nominal a una corriente de campo de 297 amperes.

5 PROBLEMAS: OPCIONALES 11
Ingenieria Electrónica UNSAAC 2019-II

a) Determine la potencia máxima Pmax (en MW y por unidad) que puede abastecer
a un sistema externo si el voltaje interno del generador se mantiene igual a 1.0 por
unidad.
b) Mediante MatLab, grafique el voltaje terminal del generador, mientras que la
salida de éste varía de cero a Pmas bajo las condiciones del inciso a).
c) Ahora asuma que el generador se equipa con un regulador de voltaje automático que
controla la corriente de campo con el fin de mantener constante el voltaje terminal. Si
el generador se carga a su valor nominal, calcule el ángulo de potencia correspondiente,
el voltaje terminal por unidad y la corriente de campo. Mediante MatLab, grafique Eaf
por uniodad como una función de la potencia por unidad.

5.2.1. Solución a)
A partir de la ecuación

Eaf VEQ
P = (14)
XS + XEQ
Observe que aunque éste e sun generador trifásico, no se reuiqre un factor de 3 debido a
que se realizan los cálculos en el sistema por unidad.
Debido a que la máquina opera con un voltaje terminal aproximado a su valor nominal,
pmax se debe expresar en términos de la reactancia síncrona saturada. Así se tiene
1
Pmax = = 0,633 por unidad = 47,5M W (15)
1,35 + 0,23

5.2.2. Solución b)
A partir de la figura (13), la corriente del terminal del generador se obtiene por

baf − VbEQ
E Eaf ejδ − VEQ ejδ − 1,0
Iba = = = (16)
j(XS + XEQ ) j(XS + XEQ ) j1,58
El voltaje terminal del generador se obtiene por
0,23 jδ−1,0
Vba = VbEQ + jXEQ Iba = 1,0 + (e ) (17)
1,58

Figura 13: Representación de un circuito equivalente de máquina síncrona conectada a un

sistema externo.
La figura (14) es una gráfica de MatLab apropiada. Al parecer, el voltaje terminal varía
de 1,0 a δ = 0o hasta 0,87 a δ = 90o

5 PROBLEMAS: OPCIONALES 12
Ingenieria Electrónica UNSAAC 2019-II

Figura 14: Gráfica de MatLab, voltaje terminal frente a δ para el inciso b).

5.2.3. Solución c)
Con el voltaje terminal mantenido constante a Va = 1,0 por unidad, la potencia puede
expresarse como:

Va VEQ
P = senδt = 4,35senδt (18)
XEQ
donde δt representa el ángulo del voltaje terminal con respecto a VbEQ
Para P = 1,0 por unidad, δt = 13,3o y por tanto Ib es igual a

Va ejδt − VEQ
(19)
o
Iba = = 1,007ej6,65
jXEQ
y

baf = VbEQ + j(XEQ + XS )Iba = 1,78ej62,7o

E (20)
o Eaf = 1,78 por unidad, que corresponde a una corriente de campo de If = 1,78 × 297 =
529 amperes.
El ángulo de potencia correspondiente es 62,7o .

Figura 15: Gráfica de MatLab, Eaf frente a la potencia para el inciso c).

5 PROBLEMAS: OPCIONALES 13
 Ingenieria Electrónica UNSAAC 2019-II

La figura (15) es la gráfica de MatLab apropiada. Puede verse que Eaf varía de 1,0 a
P = 0 hasta 1,78 a P = 1,0
A continuación se cita el código en MaLab

1 clc
2 clear
3 %solucion del inciso b
4 %parametros del sistema
5 Veq = 1.0;
6 Eaf = 1.0;
7 Xeq = .23;
8 Xs = 1.35;
9 %Resuelto para Va al tiempo que delta varia de 0 a 90 grados
10 for n = 1:101
11 delta(n)=(pi /2.) *(n-1) /100;
12 Ia(n) = Eaf *exp(j*delta(n)) - Veq/(j*(Xs + Xeq));
13 Va(n) = abs(Veq + j*Xeq*Ia(n));
14 degrees(n) = 180* delta(n)/pi;
15 end
16
17 %Ahora grfique los reusltados
18 plot (degrees , Va)
19 xlabel ('Angulo de potencia , delta [grados]')
20 ylabel ('Terminal de voltaje [por unidad]')
21 title ('Voltaje terminal vs Angulo de potencia para el inciso (b)')
22 fprintf ('\n\nHit any Key to continue\n')
23 pause
24
25 %solucion para el inciso c
26 %Establecer el voltaje terminal por unidad
27 Vterm = 1.0;
28 for n = 1:101
29 P(n) = (n-1) /100;
30 deltat(n) = asin(P(n)*Xeq/( Vterm*Veq));
31 Ia(n) = (Vterm*exp(j*deltat(n))-Veq)/(j*Xeq);
32 Eaf(n) = abs(Vterm + j*(Xs+Xeq)*Ia(n));
33 end
34
35 %Ahora a graficar los resultados
36
37 plot(P, Eaf)
38 xlabel ('Angulo de potencia [por unidad]')
39 ylabel ('Eaf [por unidad]')
40 title ('Eaf vs. potencia para el inciso (c)')

5 PROBLEMAS: OPCIONALES 14