# data[].attributes.attachment_link.readable_text
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%          UNIVERSIDAD T��CNICA DE COTOPAXI (SEDE LATACUNGA)               
%         FACULTAD DE CIENCIAS DE LA INGENIERIA Y APLICADAS               
%                                                                         
%                          INTEGRANTES:
%                       WIDYNSON CHARIGUAMAN
%                       JORGE CHICAIZA
%                       KEVIN MAIQUIZA
%                       HUGO YUGSI                                   
%                             TEMA:                                         
%               EJERCICIOS DE MAQUINAS SINCRONICAS
% INICIO DEL CODIGO

clc
clear all

% ITEM DEL EJERCICIO

disp('2.- Las reactancias de un generador s��ncrono de polos salientes son �������� = 1 �������� ���� �������� = 0.6 �������� ');
disp('La resistencia de armadura puede despreciarse. Calcular el voltaje generado cuando la')
disp('m��quina entrega su potencia nominal a un �������� = 0.8 atrasado, y voltaje terminal nominal')
disp('(���� = 1���0 �������� = 1���cos���1 ��������).')

% VALORES PARA CALCULAR

disp('INGRESE LOS VALORES QUE TIENE.')
xd=input('INGRESE REACTANCIA DE EJE DIRECTO : ')
xq=input('INGRESE REACTANCIA DE EJE DE CUADRATURA:')
fp=input('INGRESE FACTOR DE POTENCIA:')
V=input('INGRESE EL VOLTAJE NOMINAL:')
I=input('INGRESE LA CORRIENTE:')

%CALCULOS
angV1=0;
angI1=-acosd(fp);
fprintf('I:     %.2f    <    %.2f  \n',I,angI1);

angV = deg2rad(angV1) ;
angI = deg2rad(angI1) ;

[Vx,Vy] = pol2cart(angV,V) ;
[Ix,Iy] = pol2cart(angI,I); 

sig=-1;

Va= complex(Vx,Vy);
Ia=complex(Ix,Iy);

Vax=real(Va);
Vay=imag(Va);

Iax=real(Ia);
Iay=imag(Ia);

Xq=xq*j;
Xd=xd*j;

E= Va+(Xq*(Ia));
Ex=real(E);
Ey=imag(E);
[angE,E1] = cart2pol(Ex,Ey);

angE=rad2deg(angE);

fprintf('E:     %.2f    <     %.2f  \n',E1,angE);

Id=I*(sind(sig*(angI1-angE)));
angId=angE-90;

fprintf('Id:    %.2f    <    %.2f  \n',Id,angId);

Iq=I*(cosd(angI1-angE));
angIq=angE;

fprintf('Iq:    %.2f    <    %.2f  \n',Iq,angIq);

angIdr = deg2rad(angId);
[Idx,Idy] = pol2cart(angIdr,Id); 

Id1=complex(Idx,Idy);


angIqr = deg2rad(angIq); 
[Iqx,Iqy] = pol2cart(angIqr,Iq) ;

Iq1=complex(Iqx,Iqy);

Eaf=1+(Xd*Id1)+(Xq*Iq1);
Eafx=real(Eaf);
Eafy=imag(Eaf);

Eaf1=sqrt((Eafx^2)+(Eafy^2));
angEaff=atand(Eafy/Eafx);

fprintf('Eaf:   %.2f    <   %.2f  \n',Eaf1,angEaff);

# data[].attributes.attachment_link.readable_text
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%          UNIVERSIDAD TÉCNICA DE COTOPAXI (SEDE LATACUNGA)               
%         FACULTAD DE CIENCIAS DE LA INGENIERIA Y APLICADAS               
%                                                                         
%                          INTEGRANTES:
%                       WIDYNSON CHARIGUAMAN
%                       JORGE CHICAIZA
%                       KEVIN MAIQUIZA
%                       HUGO YUGSI                                   
%                             TEMA:                                         
% EJERCICIOS DE MAQUINAS SINCRONICAS DEL TALLER 
% INICIO DEL CODIGO


% ITEM DEL EJERCICIO

disp('1.-Un generador síncrono de 75 𝑀𝑉𝐴, 13.8 𝑘𝑉 con reactancia síncrona NO saturada ');
disp('𝑋𝑠 = 1.35 𝑝𝑢, se conecta a un sistema cuya reactancia equivalente es 𝑋𝐸𝑄 = 0.23 y')
disp('𝑉𝐸𝑄 = 1.')
disp('Alcanza el voltaje de circuito abierto nominal a una corriente de campo de 297 𝐴.')
disp('-Determinar la potencia máxima(en pu y en MW) que puede suplir al sistema externo, si el voltaje')
disp(' interno se mantiene constante a 1 pusi el voltaje interno se mantiene constante a 1 𝑝𝑢.')

% VALORES PARA CALCULAR

disp('INGRESE LOS VALORES QUE TIENE.')
Sg=input('INGRESE LA POTENCIA APARENTE: ')
xs=input('INGRESE LA REACTANCIA SINCRONA:')
xeq=input('INGRESE LA REACTANCIA EQUIVALENTE :')
veq=input('INGRESE EL VOLTAJE EQUIVALENTEW:')
Eaf=input('INGRESE LA TENSION :')

%CALCULOS

Pmax=Eaf*veq/(xs)

fprintf('Pmax  %.1f \n',Pmax)


# data[].attributes.attachment_link.readable_text
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%          UNIVERSIDAD TÉCNICA DE COTOPAXI (SEDE LATACUNGA)               
%         FACULTAD DE CIENCIAS DE LA INGENIERIA Y APLICADAS               
%                                                                         
%                          INTEGRANTES:
%                       WIDYNSON CHARIGUAMAN
%                       JORGE CHICAIZA
%                       KEVIN MAIQUIZA
%                       HUGO YUGSI                                   
%                             TEMA:                                         
%               EJERCICIOS DE MAQUINAS SINCRONICAS
% INICIO DEL CODIGO

clc
clear all
% ITEM DEL EJERCICIO

disp('3.-Un motor síncrono de polos salientes de 1490 𝑘𝑊, 2300 𝑉, tiene como reactancias:');
disp('𝑋𝑑 = 1.95 Ω 𝑦 𝑋𝑞 = 1.4 Ω.')
disp('Despreciando todas las pérdidas, y considerando que está siendo alimentado desde un bus')
disp('infinito, y su voltaje de campo se mantiene constante a un valor que resultará')
disp('en un 𝑓𝑝 = 1.')
disp('-Determinar el ángulo 𝛿 que corresponde a la máxima potencia mecánica.')
disp('-Calcular la máxima potencia mecánica que puede entregar')

% VALORES PARA CALCULAR

disp('INGRESE LOS VALORES QUE TIENE.')
xd=input('INGRESE REACTANCIA DE EJE DIRECTO : ')
xq=input('INGRESE REACTANCIA DE EJE DE CUADRATURA:')
P=input('INGRESE LA POTENCIA ACTIVA:')
V=input('INGRESE EL VOLTAJE NOMINAL:')

%CALCULOS

Va=V/(sqrt(3));

Ia=(P/1000)/(3*(Va/1000));

Xq=xq*j;
Xd=xd*j;

E=Va-(Xq*Ia);
Ex=real(E);
Ey=imag(E);
E1=sqrt((Ex^2)+(Ey^2));
angEg=atand(Ey/Ex);


fprintf('ang: %.1f \n',angEg)


Id=Ia*sind(-angEg);

fprintf('Id:  %.1f \n',Id)


Iq=Ia*cosd(-angEg);

fprintf('Iq:  %.1f \n',Iq)

Eaf=Va-(Xd*Id)-(Xq*Iq);
Eafx=real(Eaf);
Eafy=imag(Eaf);

Eafm=sqrt((Eafx^2)+(Eafy^2));
angEr=atand(Ey/Ex);

fprintf('Eaf:   %.1f \n',Eafm)

Px=(Eafm*Va/xd)*sin(angEg);
Py=(Va^2*((xd-xq)/(2*xd*xq)))*sin(2*angEg);

fprintf('P:     %.1f        %.1f  \n',Px,Py)

# data[].attributes.attachment_link.readable_text
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%          UNIVERSIDAD TÉCNICA DE COTOPAXI (SEDE LATACUNGA)               
%         FACULTAD DE CIENCIAS DE LA INGENIERIA Y APLICADAS               
%                                                                         
%                          INTEGRANTES:
%                       WIDYNSON CHARIGUAMAN
%                       JORGE CHICAIZA
%                       KEVIN MAIQUIZA
%                       HUGO YUGSI                                   
%                             TEMA:                                         
%               EJERCICIOS DE MAQUINAS SINCRONICAS
% INICIO DEL CODIGO

clc
clear all

% ITEM DEL EJERCICIO

disp('2.- Las reactancias de un generador síncrono de polos salientes son 𝑋𝑑 = 1 𝑝𝑢 𝑦 𝑋𝑞 = 0.6 𝑝𝑢 ');
disp('La resistencia de armadura puede despreciarse. Calcular el voltaje generado cuando la')
disp('máquina entrega su potencia nominal a un 𝑓𝑝 = 0.8 atrasado, y voltaje terminal nominal')
disp('(𝑉 = 1∠0 𝐼𝑎 = 1∠cos−1 𝑓𝑝).')

% VALORES PARA CALCULAR

disp('INGRESE LOS VALORES QUE TIENE.')
xd=input('INGRESE REACTANCIA DE EJE DIRECTO : ')
xq=input('INGRESE REACTANCIA DE EJE DE CUADRATURA:')
fp=input('INGRESE FACTOR DE POTENCIA:')
V=input('INGRESE EL VOLTAJE NOMINAL:')
I=input('INGRESE LA CORRIENTE:')

%CALCULOS
angV1=0;
angI1=-acosd(fp);
fprintf('I:     %.2f    <    %.2f  \n',I,angI1);

angV = deg2rad(angV1) ;
angI = deg2rad(angI1) ;

[Vx,Vy] = pol2cart(angV,V) ;
[Ix,Iy] = pol2cart(angI,I); 

sig=-1;

Va= complex(Vx,Vy);
Ia=complex(Ix,Iy);

Vax=real(Va);
Vay=imag(Va);

Iax=real(Ia);
Iay=imag(Ia);

Xq=xq*j;
Xd=xd*j;

E= Va+(Xq*(Ia));
Ex=real(E);
Ey=imag(E);
[angE,E1] = cart2pol(Ex,Ey);

angE=rad2deg(angE);

fprintf('E:     %.2f    <     %.2f  \n',E1,angE);

Id=I*(sind(sig*(angI1-angE)));
angId=angE-90;

fprintf('Id:    %.2f    <    %.2f  \n',Id,angId);

Iq=I*(cosd(angI1-angE));
angIq=angE;

fprintf('Iq:    %.2f    <    %.2f  \n',Iq,angIq);

angIdr = deg2rad(angId);
[Idx,Idy] = pol2cart(angIdr,Id); 

Id1=complex(Idx,Idy);


angIqr = deg2rad(angIq); 
[Iqx,Iqy] = pol2cart(angIqr,Iq) ;

Iq1=complex(Iqx,Iqy);

Eaf=1+(Xd*Id1)+(Xq*Iq1);
Eafx=real(Eaf);
Eafy=imag(Eaf);

Eaf1=sqrt((Eafx^2)+(Eafy^2));
angEaff=atand(Eafy/Eafx);

fprintf('Eaf:   %.2f    <   %.2f  \n',Eaf1,angEaff);

# data[].attributes.attachment_link.readable_text
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%          UNIVERSIDAD T��CNICA DE COTOPAXI (SEDE LATACUNGA)               
%         FACULTAD DE CIENCIAS DE LA INGENIERIA Y APLICADAS               
%                                                                         
%                          INTEGRANTES:
%                       WIDYNSON CHARIGUAMAN
%                       JORGE CHICAIZA
%                       KEVIN MAIQUIZA
%                       HUGO YUGSI                                   
%                            TEMA                                         
% MEDIANTE EL MODELO MATEM��TICO DE UN GENERADOR S��NCRON0 REALIZAR UN 
% EJEMPLO DE CALCULOS PARA DETERMINADOS VALORES DEL GENERADOR 

% INICIO DEL CODIGO
clc
clear all
disp('EJERCICIOS DE GENERADOR SINCRONO');
selet=input('INGRESE UNA OPCION: ') % AL CORRER EL CODIGO EN EL COMMAND WINDOW INSERTA "1" O EN RUN
if(selet==1)
    % INGRESAR VALORES DETERMINADOS
    
    Vl=input('INGRESE VOLTAJE DE LINEA: ')
    S=input('INGRESE POTENCIA APARENTE: ')
    Ra=input('INGRESE RESISTENCIA DEL INDUCIDO: ')
    Xs=input('INGRESE LA REACTANCIA SINCRONA: ')
    F=input('INGRESE LA FRECUENCIA: ')
    p=input('INGRESE N��MERO DE POLOS: ')
    Fp=input('FACTOR DE POTENCIA: ')
    
% FORMULAS DEL GENERADOR SINCRONICO

% TENSI��N POR FASE 
  Vo=Vl/sqrt(3);
  
% ��NGULO DE DESFASE
  Ang=acos(Fp);
  
% CORRIENTE DEL INDUCIDO
  Ia=S/(sqrt(3)*Vl) ;
  
%T��NSION INTERNA DE LA M��QUINA
  Ea=Vo+Xs*Ia+Ra*Ia;
  
%POTENCIA ACTIVA
  P=sqrt(3)*Vo*Ia*Fp;
  
%POTENCIA REACTIVA
  Q=sqrt(3)*Vo*Ia*sin(Ang);
  
%VELOCIDAD DE SINCRONISMO
  V= (120*F)/p;
  
%VALORES CALCULADOS
  disp('LA CORRIENTE DEL INDUCIDO:')
    dIa=['Ia=', num2str(Ia),'A'];
  disp(dIa)
  
    disp('EL TENSI��N INTERNO DE LA M��QUINA:')
  dEa=['Ea=', num2str(Ea),'V'];
  disp(dEa)
  
  disp('VELOCIDAD DE SINCRONISMO:') 
    dV=['V=', num2str(V),'rpm'];
  disp(dV)
  
   disp('LA POTENCIA ACTIVA:')
    dP=['P=', num2str(P),'W'];
  disp(dP)
  
  disp('LA POTENCIA REACTIVA:')
    dQ=['Q=', num2str(Q),'VAr'];
  disp(dQ)
 
 end
    