clc;clear all;format short %valores base SB=100e6; VB3=138e3; VB1=VB3*10/138; VB5=VB3*13.8/138; VB7=VB5; ZB=VB3^2/SB; disp('///////////////// PARTE A \\\\\\\\\\\\\\\\\') disp('DIBUJAR LAS REDES DE SECUENCIA CONSIDERANDO LAS BASES MOSTRADAS') disp('************** Generador 1 ***************') XG1_1=i*.275*SB/27.5e6*(10e3/VB1)^2 XG1_2=i*.165*SB/27.5e6*(10e3/VB1)^2 XG1_0=i*.11*SB/27.5e6*(10e3/VB1)^2 disp('************** Generador 2 ***************') XG2_1=i*.285*SB/20e6*(13.8e3/VB7)^2 XG2_2=i*.175*SB/20e6*(13.8e3/VB7)^2 XG2_0=i*.12*SB/20e6*(13.8e3/VB7)^2 disp('*********** Transformador 1 ***************') XT1=i*.0966*SB/27.6e6*(10e3/VB1)^2 XTP1_012=XT1/5 XTS1_012=4*XT1/5 disp('*********** Transformador 2 ***************') XT2=i*.09625*SB/27.5e6*(13.8e3/VB5)^2 XTP2_012=4*XT2/5 XTS2_012=XT2/5 disp('**************** Linea 1 ******************') XL1=(0.092+i*.484)*69.96/ZB XL1_0=XL1*3 disp('**************** Linea 2 ******************') XL2=(0.088+i*.457)*72/ZB XL2_0=XL2*3 disp('********** Carga 1 *********') SC1=(20+i*10)*1e6/SB; Z1=conj((13.8e3/VB7)^2/SC1) disp('***********SEIN ************') ZSEIN=i*1/conj(12000*1e6/SB) disp('***********impedancias shunt *********') %Las admitancias shunt para este caso han sido transformadas a impedancias %shunt y posteriormente llevadas al sistema pu Xsh34=(1/((i*3.4263*10^-6)*69.96))/ZB; Xsh36=Xsh34; Xsh4=Xsh34 Xsh68=(1/((i*3.0678*10^-6)*72))/ZB; Xsh3=Xsh34/2 %shunt en la barra3 Xsh6=Xsh34*Xsh68/(Xsh34+Xsh68) %shunt en la barra6 Xsh8=Xsh68 %% SECUENCIA POSITIVA Zdatos=[ %ramas R X 0 1 real(XG1_1) imag(XG1_1); 1 3 real(XT1) imag(XT1); 0 2 real(XG1_1) imag(XG1_1); 2 3 real(XT1) imag(XT1); 3 4 real(XL1) imag(XL1); 3 6 real(XL1) imag(XL1); 6 7 real(XT2) imag(XT2); 4 5 real(XT2) imag(XT2); 0 5 real(Z1) imag(Z1); 0 7 real(XG2_1) imag(XG2_1) 0 8 real(ZSEIN) imag(ZSEIN); 6 8 real(XL2) imag(XL2); 0 4 real(Xsh4) imag(Xsh4); 0 3 real(Xsh3) imag(Xsh3); 0 6 real(Xsh6) imag(Xsh6); 0 8 real(Xsh8) imag(Xsh8); ]; Rdatos=Zdatos(:,3); Xdatos=Zdatos(:,4); Z=Rdatos+i*Xdatos;Y=Z.^-1; nrbus=max(max([Zdatos(:,1) Zdatos(:,2)])); nrramas=length(Zdatos);Ybus=zeros(nrbus); %DIAGONAL PRINCIPAL for n=1:nrbus
for x=1:nrramas if Zdatos(x,1)==n||Zdatos(x,2)==n Ybus(n,n)=Ybus(n,n)+Y(x); else if Zdatos(x,1)~=0&&Zdatos(x,2)~=0 Ybus(Zdatos(x,1),Zdatos(x,2))=Y(x);Ybus(Zdatos(x,2),Zdatos(x,1))=Ybus(Zdatos(x,1), Zdatos(x,2)); end end end end disp(' *************************************** secuencia positiva ********************************************') disp(' ********************************************* Z BARRA ***************************************************' ) fprintf('\n') disp(Ybus^-1) fprintf('\n') %% SECUENCIA NEGATIVA Zdatos=[ %ramas R X 0 1 real(XG1_2) imag(XG1_2); 1 3 real(XT1) imag(XT1); 0 2 real(XG1_2) imag(XG1_2); 2 3 real(XT1) imag(XT1); 3 4 real(XL1) imag(XL1); 3 6 real(XL1) imag(XL1); 6 7 real(XT2) imag(XT2); 4 5 real(XT2) imag(XT2); 0 5 real(Z1) imag(Z1); 0 7 real(XG2_2) imag(XG2_2) 0 8 real(ZSEIN) imag(ZSEIN); 6 8 real(XL2) imag(XL2); 0 4 real(Xsh4) imag(Xsh4); 0 3 real(Xsh3) imag(Xsh3); 0 6 real(Xsh6) imag(Xsh6); 0 8 real(Xsh8) imag(Xsh8); ]; Rdatos=Zdatos(:,3); Xdatos=Zdatos(:,4); Z=Rdatos+i*Xdatos;Y=Z.^-1; nrbus=max(max([Zdatos(:,1) Zdatos(:,2)])); nrramas=length(Zdatos);Ybus=zeros(nrbus); %DIAGONAL PRINCIPAL for n=1:nrbus for x=1:nrramas if Zdatos(x,1)==n||Zdatos(x,2)==n Ybus(n,n)=Ybus(n,n)+Y(x); else if Zdatos(x,1)~=0&&Zdatos(x,2)~=0 Ybus(Zdatos(x,1),Zdatos(x,2))=Y(x);Ybus(Zdatos(x,2),Zdatos(x,1))=Ybus(Zdatos(x,1), Zdatos(x,2)); end end end end disp(' *************************************** secuencia negativa ********************************************') disp(' ********************************************* Z BARRA ***************************************************' ) fprintf('\n') disp(Ybus^-1) fprintf('\n') %% SECUENCIA CERO Zdatos=[ %ramas R X 0 1 real(XTS1_012) imag(XTS1_012); 0 1 real(XTS1_012) imag(XTS1_012); 1 2 real(XL1_0) imag(XL1_0);
1 3 real(XL1_0) imag(XL1_0); 0 2 real(XTP2_012) imag(XTP2_012); 0 4 real(XG2_0) imag(XG2_0) 0 5 real(ZSEIN) imag(ZSEIN); 3 5 real(XL2_0) imag(XL2_0); 3 4 real(XT2) imag(XT2); 0 2 real(Xsh4) imag(Xsh4); 0 1 real(Xsh3) imag(Xsh3); 0 3 real(Xsh6) imag(Xsh6); 0 5 real(Xsh8) imag(Xsh8); ]; Rdatos=Zdatos(:,3); Xdatos=Zdatos(:,4); Z=Rdatos+i*Xdatos;Y=Z.^-1; nrbus=max(max([Zdatos(:,1) Zdatos(:,2)])); nrramas=length(Zdatos);Ybus=zeros(nrbus); %DIAGONAL PRINCIPAL for n=1:nrbus for x=1:nrramas if Zdatos(x,1)==n||Zdatos(x,2)==n Ybus(n,n)=Ybus(n,n)+Y(x); else if Zdatos(x,1)~=0&&Zdatos(x,2)~=0 Ybus(Zdatos(x,1),Zdatos(x,2))=Y(x);Ybus(Zdatos(x,2),Zdatos(x,1))=Ybus(Zdatos(x,1), Zdatos(x,2)); end end end end disp(' *************************************** secuencia CERO ********************************************') disp(' ********************************************* Z BARRA ***********************************************') fprintf('\n') disp(Ybus^-1) fprintf('\n') RESULTADOS ///////////////// PARTE A \\\\\\\\\\\\\\\\\
XT1 = 0.0000 + 0.3500i XTP1_012 = 0.0000 + 0.0700i XTS1_012 = 0.0000 + 0.2800i *********** Transformador 2 *************** XT2 = 0.0000 + 0.3500i XTP2_012 = 0.0000 + 0.2800i XTS2_012 = 0.0000 + 0.0700i **************** Linea 1 ****************** XL1 = 0.0338 + 0.1778i XL1_0 = 0.1014 + 0.5334i **************** Linea 2 ****************** XL2 = 0.0333 + 0.1728i XL2_0 =
DIBUJAR LAS REDES DE SECUENCIA CONSIDERANDO LAS BASES MOSTRADAS 0.0998 + 0.5183i ************** Generador 1 *************** ********** Carga 1 ********* XG1_1 = Z1 = 0.0000 + 1.0000i 4.0000 + 2.0000i XG1_2 = ***********SEIN ************ 0.0000 + 0.6000i ZSEIN = XG1_0 = 0.0000 + 0.0083i 0.0000 + 0.4000i ***********impedancias shunt ********* ************** Generador 2 *************** Xsh4 = XG2_1 = 0.0000 -21.9062i 0.0000 + 1.4250i Xsh3 = XG2_2 = 0.0000 -10.9531i 0.0000 + 0.8750i Xsh6 = XG2_0 = 0.0000 -11.4007i 0.0000 + 0.6000i Xsh8 = *********** Transformador 1 *************** 0.0000 -23.7729i
*************************************** secuencia positiva ******************************************** ********************************************* Z BARRA ***************************************************
0.0201 + 0.3846i 0.0201 + 0.1254i 0.0272 + 0.1693i 0.0317 + 0.1654i 0.0413 + 0.1570i 0.0123 + 0.0820i 0.0099 + 0.0658i -0.0001 + 0.0038i 0.0201 + 0.1254i 0.0201 + 0.3846i 0.0272 + 0.1693i 0.0317 + 0.1654i 0.0413 + 0.1570i 0.0123 + 0.0820i 0.0099 + 0.0658i -0.0001 + 0.0038i 0.0272 + 0.1693i 0.0272 + 0.1693i 0.0367 + 0.2285i 0.0428 + 0.2233i 0.0557 + 0.2119i 0.0166 + 0.1107i 0.0133 + 0.0889i -0.0002 + 0.0051i 0.0317 + 0.1654i 0.0317 + 0.1654i 0.0428 + 0.2233i 0.0872 + 0.3915i 0.1093 + 0.3708i 0.0196 + 0.1082i 0.0157 + 0.0869i -0.0000 + 0.0050i 0.0413 + 0.1570i 0.0413 + 0.1570i 0.0557 + 0.2119i 0.1093 + 0.3708i 0.1520 + 0.6862i 0.0259 + 0.1028i 0.0208 + 0.0825i 0.0003 + 0.0048i 0.0123 + 0.0820i 0.0123 + 0.0820i 0.0166 + 0.1107i 0.0196 + 0.1082i 0.0259 + 0.1028i 0.0229 + 0.1398i 0.0184 + 0.1122i -0.0001 + 0.0064i 0.0099 + 0.0658i 0.0099 + 0.0658i 0.0133 + 0.0889i 0.0157 + 0.0869i 0.0208 + 0.0825i 0.0184 + 0.1122i 0.0147 + 0.3711i -0.0001 + 0.0051i -0.0001 + 0.0038i -0.0001 + 0.0038i -0.0002 + 0.0051i -0.0000 + 0.0050i 0.0003 + 0.0048i -0.0001 + 0.0064i -0.0001 + 0.0051i 0.0000 + 0.0082i
*************************************** secuencia negativa ******************************************** ********************************************* Z BARRA ***************************************************
0.0110 + 0.3001i 0.0110 + 0.0791i 0.0174 + 0.1252i 0.0208 + 0.1224i 0.0280 + 0.1164i 0.0074 + 0.0594i 0.0053 + 0.0424i -0.0002 + 0.0027i 0.0110 + 0.0791i 0.0110 + 0.3001i 0.0174 + 0.1252i 0.0208 + 0.1224i 0.0280 + 0.1164i 0.0074 + 0.0594i 0.0053 + 0.0424i -0.0002 + 0.0027i 0.0174 + 0.1252i 0.0174 + 0.1252i 0.0275 + 0.1983i 0.0329 + 0.1939i 0.0443 + 0.1843i 0.0117 + 0.0940i 0.0083 + 0.0672i -0.0002 + 0.0043i 0.0208 + 0.1224i 0.0208 + 0.1224i 0.0329 + 0.1939i 0.0766 + 0.3629i 0.0973 + 0.3441i 0.0143 + 0.0920i 0.0102 + 0.0657i -0.0001 + 0.0042i 0.0280 + 0.1164i 0.0280 + 0.1164i 0.0443 + 0.1843i 0.0973 + 0.3441i 0.1387 + 0.6612i 0.0197 + 0.0875i 0.0141 + 0.0625i 0.0002 + 0.0041i 0.0074 + 0.0594i 0.0074 + 0.0594i 0.0117 + 0.0940i 0.0143 + 0.0920i 0.0197 + 0.0875i 0.0196 + 0.1290i 0.0140 + 0.0921i -0.0002 + 0.0059i 0.0053 + 0.0424i 0.0053 + 0.0424i 0.0083 + 0.0672i 0.0102 + 0.0657i 0.0141 + 0.0625i 0.0140 + 0.0921i 0.0100 + 0.3158i -0.0001 + 0.0042i -0.0002 + 0.0027i -0.0002 + 0.0027i -0.0002 + 0.0043i -0.0001 + 0.0042i 0.0002 + 0.0041i -0.0002 + 0.0059i -0.0001 + 0.0042i 0.0000 + 0.0082i
*************************************** secuencia CERO ******************************************** ********************************************* Z BARRA ***********************************************
0.0038 + 0.1067i -0.0032 + 0.0367i -0.0002 + 0.0422i -0.0001 + 0.0267i -0.0001 + 0.0006i -0.0032 + 0.0367i 0.0094 + 0.1991i -0.0018 + 0.0144i -0.0012 + 0.0091i -0.0001 + 0.0002i -0.0002 + 0.0422i -0.0018 + 0.0144i 0.0311 + 0.2290i 0.0197 + 0.1446i -0.0002 + 0.0036i -0.0001 + 0.0267i -0.0012 + 0.0091i 0.0197 + 0.1446i 0.0124 + 0.3124i -0.0001 + 0.0023i -0.0001 + 0.0006i -0.0001 + 0.0002i -0.0002 + 0.0036i -0.0001 + 0.0023i 0.0000 + 0.0083i