Clc1imp.docx

  • Uploaded by: Fernando Medina Garcia
  • 0
  • 0
  • December 2019
  • PDF

This document was uploaded by user and they confirmed that they have the permission to share it. If you are author or own the copyright of this book, please report to us by using this DMCA report form. Report DMCA


Overview

Download & View Clc1imp.docx as PDF for free.

More details

  • Words: 1,062
  • Pages: 3
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

More Documents from "Fernando Medina Garcia"