% ----------------------------------------------------------------------- clc variasejecuciones=1; while variasejecuciones==1 clear % ----------------------------------------------------------------------- fprintf('\n'); fprintf('PG128Mn\n'); fprintf('Esquema Linea Doble SR\n'); fprintf('Unominal AT kV: 220\n'); fprintf('Unominal BT kV: 66\n'); fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n'); fprintf('Eleccion Ctes. del Tf\n'); fprintf('Tf AT/BT/T: Neutro a Tierra en AT y BT\n'); fprintf('Eleccion: kE,LA=LB=L,Sccs,Sccr,Kr0,k0\n'); fprintf('*** Parametro k0=zs0/zr0=0,1-1-10 ***\n'); fprintf('Impedancia de Transferencia Si, No\n'); fprintf('Eleccion: (zsr1,Ksr0)\n'); fprintf('Compensacion 21, Si, No\n'); fprintf('Falta Intercircuito entre Lineas de Diferente Tension\n'); fprintf('m=p=0...1\n'); fprintf('Barra m=1, Barra p=2\n'); %fprintf('Salida Grafica 1: Analisis SAm, Elegir Unidad (pu)\n'); %fprintf(' : Analisis SBp, Elegir Unidad (pu)\n'); %fprintf('Salida Grafica 2: Analisis SAm, Elegir Unidad (pua)\n'); %fprintf(' : Analisis SBp, Elegir Unidad (pua)\n'); fprintf('\n'); % ----------------------------------------------------------------------- %PG128Mc datos=input('datos: standard=1,input=2: '); % ----------------------------------------------------------------------- if datos==1 UA=220; UB=66; kE=1; L=66; Sp=120; Ss=120; St=40; Up=220; Us=66; Ut=33; tccps=12; tccpt=5; tccst=2; Sccs=5000; Sccr=5000; Kr0=1; imptransfer=2; compens21=2; % ------------------------------------------------------------------------ else UA=input('UA(Fase-Fase) kV, poner 220: '); UB=input('UB(Fase-Fase) kV, poner 66: '); Ut=input('Ut(Fase-Fase) kV, poner 33: '); kE=input('kE= '); L=input('Longitud L=LA=LB km: '); Sccs=input('Sccs MVA: '); Sccr=input('Sccr MVA: '); Kr0=input('Fuente R, zr0/zr1: '); Sp=input('Sp MVA, poner 120: '); Ss=input('Ss MVA, poner 120: '); St=input('St MVA, poner 40: '); tccps=input('tccps %, Ref a Sp, poner 12: '); tccpt=input('tccpt %, Ref a St, poner 5: '); tccst=input('tccst %, Ref a St, poner 2: '); imptransfer=input('imptransfer:si=1,no=2:'); compens21=input('Compensacion 21: Si=1, No=2 :'); end % ------------------------------------------------------------------------ %k0=input('k0=zs0/zr0 (Tipicos=0.1-1-10): '); fprintf('Tipo de Falta: \n'); fprintf('Aa=1,Bc=2,Cb=3,Aag=4,Bcg=5,Cbg=6,Aab=7,Bac=8,Cbc=9\n'); fprintf('Aabg=10,Bacg=11,Cbcg=12,Aabc=13,Aabcg=14,ABCabc=15,ABCabcg=16\n'); falta=input('Elegir Tipo de Falta: '); Rf=input('Resistencia de Falta. Rf(Ohm) = '); if falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16 Rg=input('Resistencia comun de Falta. Rg(Ohm) = '); else Rg=0; end % ------------------------------------------------------------------------ if UA==400 zlA1u=0.0346+0.3208i; zlA2u=zlA1u; zlA0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i; elseif UA==220 zlA1u=0.0829+0.4135i; zlA2u=zlA1u; zlA0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i; elseif UA==132 zlA1u=0.1342+0.4042i; zlA2u=zl1u; zlA0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i; elseif UA==66 zlA1u=0.2156+0.3952i; zlA2u=zl1u; zlA0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i; end if UB==400 zlB1u=0.0346+0.3208i; zlB2u=zlB1u; zlB0u=0.1828+1.1482i; elseif UB==220 zlB1u=0.0829+0.4135i; zlB2u=zlB1u; zlB0u=0.2311+1.3085i; elseif UB==132 zlB1u=0.1342+0.4042i; zlB2u=zlB1u; zlB0u=0.2824+1.3718i; elseif UB==66 zlB1u=0.2156+0.3952i; zlB2u=zlB1u; zlB0u=0.3636+1.4278i; end % ------------------------------------------------------------------------ LA=L; LB=L; %EAv=UA*1000/sqrt(3)+0i; EA=1.0; %EBv=UB*1000/sqrt(3)+0i; EB=1.0; E=1.0; E=kE*E; kA=(zlA0u-zlA1u)/(3*zlA1u); kB=(zlB0u-zlB1u)/(3*zlB1u); Sbase=100; % MVA UbaseA=UA; % kV UbaseB=UB; % kV VbaseA=UbaseA/sqrt(3); % kV VbaseB=UbaseB/sqrt(3); % kV zbaseA=UA^2/Sbase; % Ohm zbaseB=UB^2/Sbase; % Ohm IbaseA=(UA/sqrt(3))/zbaseA; % kA IbaseB=(UB/sqrt(3))/zbaseB; %kA KmA=zlM0u/zlA1u; KmA=KmA*(IbaseB/IbaseA); KmB=zlM0u/zlB1u; KmB=KmB*(IbaseA/IbaseB); fprintf('\n'); if imptransfer==1 zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): '); zsr1=zsr1/zbaseA; zsr2=zsr1; Ksr0=input('Ksr0=zsr0/zsr1= '); zsr0=Ksr0*zsr1; end if falta==2 cfalta='Bc';elseif falta==3 cfalta='Cb';elseif falta==5 cfalta='Bcg'; elseif falta==6 cfalta='Cbg';elseif falta==7 cfalta='Aab';elseif falta==8 cfalta='Bac'; elseif falta==9 cfalta='Cbc';elseif falta==10 cfalta='Aabg';elseif falta==11 cfalta='Bacg'; elseif falta==12 cfalta='Cbcg';elseif falta==13 cfalta='Aabc';elseif falta==14 cfalta='Aabcg'; elseif falta==15 cfalta='ABCabc';elseif falta==16 cfalta='ABCabcg'; end % ------------------------------------------------------------------------ fila=1; k0=0.1; while k0<11 i=1; m=0.01; while m<1; p=m; % REDES DE SECUENCIA 0,1,2 pu [zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1]=pusec1(UA,Sccs,Sccr,zlA1u,m,LA,zlB1u,p,LB,zbaseA,zbaseB); [zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2]=pusec2(UA,Sccs,Sccr,zlA2u,m,LA,zlB2u,p,LB,zbaseA,zbaseB); [zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0]=... pusec0(k0,Kr0,zr1,zlA0u,m,LA,zlB0u,p,LB,zlM0u,zbaseA,zbaseB,UbaseA,UbaseB,Sbase); zpTS=1/2*((tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St-(tccst/100)*UA^2/St); zpTS=zpTS/zbaseA; zsTS=1/2*((tccps/100)*UA^2/Sp-(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St); zsTS=zsTS/zbaseA; ztTS=1/2*(-(tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St); ztTS=ztTS/zbaseA; zpTR=zpTS; zsTR=zsTS; ztTR=ztTS; % CALCULO DE LA REDUCCION DE REDES 0,1,2 A LOS PUNTOS m,p DE FALTA switch(imptransfer) case{1}, [Minc0]=subMinc0T(1); [Minc1]=subMinc1T(1); [Minc2]=subMinc2T(1); [z1]=subz1T(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR,zsr1); [z2]=subz2T(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR,zsr2); [z0]=subz0T(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zpTS,zsTS,ztTS,zpTR,zsTR,ztTR,zsr0); case(2) [Minc0]=subMinc0(1); [Minc1]=subMinc1(1); [Minc2]=subMinc2(1); [z1]=subz1(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR); [z2]=subz2(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR); [z0]=subz0(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zpTS,zsTS,ztTS,zpTR,zsTR,ztTR); end % MATRIZ DE ADMITANCIA DE BARRAS [Y0]=subY0(Minc0,z0); [Y1]=subY1(Minc1,z1); [Y2]=subY2(Minc2,z2); % ZBUS(0,1,2) Z0=inv(Y0); Z1=inv(Y1); Z2=inv(Y2); % CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA [Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2]=... maMDmp(Z0,Z1,Z2,E,falta,Rf,Rg,VbaseA,VbaseB,IbaseA,IbaseB); % IBUS(0,1,2) I0=[-Im0;-Ip0;0;0;0;0;0;0]; I1=[-Im1;-Ip1;0;0;0;0]; I2=[-Im2;-Ip2;0;0;0;0]; % VBUS(0,1,2) V0=Z0*I0; V1=Z1*I1+E; V2=Z2*I2; % I(0,1,2) DE LOS ELEMENTOS [Ie0]=subIe0(Minc0,V0,z0); [Ie1]=subIe1(Minc1,V1,z1); [Ie2]=subIe2(Minc2,V2,z2); VSA0=V0(3); VSA1=V1(3); VSA2=V2(3); VRA0=V0(4); VRA1=V1(4); VRA2=V2(4); VSB0=V0(5); VSB1=V1(5); VSB2=V2(5); VRB0=V0(6); VRB1=V1(6); VRB2=V2(6); ISAm0=Ie0(3); ISAm1=Ie1(3); ISAm2=Ie2(3); IRAm0=Ie0(4); IRAm1=Ie1(4); IRAm2=Ie2(4); ISBp0=Ie0(7); ISBp1=Ie1(7); ISBp2=Ie2(7); IRBp0=Ie0(8); IRBp1=Ie1(8); IRBp2=Ie2(8); % CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS PUNTOS m,p DE FALTA [Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); [Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2); [Ima,Imb,Imc]=abc(Im0,Im1,Im2); [Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2); % CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA [VSAa,VSAb,VSAc]=abc(VSA0,VSA1,VSA2); [VRAa,VRAb,VRAc]=abc(VRA0,VRA1,VRA2); [VSBa,VSBb,VSBc]=abc(VSB0,VSB1,VSB2); [VRBa,VRBb,VRBc]=abc(VRB0,VRB1,VRB2); [ISAma,ISAmb,ISAmc]=abc(ISAm0,ISAm1,ISAm2); [ISBpa,ISBpb,ISBpc]=abc(ISBp0,ISBp1,ISBp2); [IRAma,IRAmb,IRAmc]=abc(IRAm0,IRAm1,IRAm2); [IRBpa,IRBpb,IRBpc]=abc(IRBp0,IRBp1,IRBp2); % % PROTECCION DE DISTANCIA (21) switch(compens21) case{1}, % COMPENSADA [SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=... PDistComp(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA,KmA,ISBp0); [RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG]=... PDistComp(VRAa,VRAb,VRAc,IRAma,IRAmb,IRAmc,IRAm0,kA,KmA,IRBp0); [SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=... PDistComp(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB,KmB,ISAm0); [RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG]=... PDistComp(VRBa,VRBb,VRBc,IRBpa,IRBpb,IRBpc,IRBp0,kB,KmB,IRAm0); case{2}, % NO COMPENSADA [SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=... PDist(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA); [RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG]=... PDist(VRAa,VRAb,VRAc,IRAma,IRAmb,IRAmc,IRAm0,kA); [SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=... PDist(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB); [RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG]=... PDist(VRBa,VRBb,VRBc,IRBpa,IRBpb,IRBpc,IRBp0,kB); end % CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION % DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA % [VSApol,ISAmop,AngDirSAm]=PDire(VSA0,ISAm0); % pu [VSBpol,ISBpop,AngDirSBp]=PDire(VSB0,ISBp0); [VRApol,IRAmop,AngDirRAm]=PDire(VRA0,IRAm0); [VRBpol,IRBpop,AngDirRBp]=PDire(VRB0,IRBp0); % % Conversion I(0,1,2)pu a I(0,1,2)kA [IIm0,IIm1,IIm2]=conv3(Im0,Im1,Im2,IbaseA); [IIp0,IIp1,IIp2]=conv3(Ip0,Ip1,Ip2,IbaseB); [IISAm0,IISAm1,IISAm2,IIRAm0,IIRAm1,IIRAm2]=... conv6(ISAm0,ISAm1,ISAm2,IRAm0,IRAm1,IRAm2,IbaseA); [IISBp0,IISBp1,IISBp2,IIRBp0,IIRBp1,IIRBp2]=... conv6(ISBp0,ISBp1,ISBp2,IRBp0,IRBp1,IRBp2,IbaseB); % Conversion V(0,1,2)pu a V(0,1,2)kV [VVm0,VVm1,VVm2,VVSA0,VVSA1,VVSA2,VVRA0,VVRA1,VVRA2]=... conv9(Vm0,Vm1,Vm2,VSA0,VSA1,VSA2,VRA0,VRA1,VRA2,VbaseA); [VVp0,VVp1,VVp2,VVSB0,VVSB1,VVSB2,VVRB0,VVRB1,VVRB2]=... conv9(Vp0,Vp1,Vp2,VSB0,VSB1,VSB2,VRB0,VRB1,VRB2,VbaseB); % Conversion I(a,b,c)pu a I(a,b,c)kA [IIma,IImb,IImc]=conv3(Ima,Imb,Imc,IbaseA); [IIpa,IIpb,IIpc]=conv3(Ipa,Ipb,Ipc,IbaseB); [IISAma,IISAmb,IISAmc,IIRAma,IIRAmb,IIRAmc]=... conv6(ISAma,ISAmb,ISAmc,IRAma,IRAmb,IRAmc,IbaseA); [IISBpa,IISBpb,IISBpc,IIRBpa,IIRBpb,IIRBpc]=... conv6(ISBpa,ISBpb,ISBpc,IRBpa,IRBpb,IRBpc,IbaseB); % Conversion V(a,b,c)pu a V(a,b,c)Volt [VVma,VVmb,VVmc]=conv3(Vma,Vmb,Vmc,VbaseA); [VVpa,VVpb,VVpc]=conv3(Vpa,Vpb,Vpc,VbaseB); [VVSAa,VVSAb,VVSAc,VVRAa,VVRAb,VVRAc]=... conv6(VSAa,VSAb,VSAc,VRAa,VRAb,VRAc,VbaseA); [VVSBa,VVSBb,VVSBc,VVRBa,VVRBb,VVRBc]=... conv6(VSBa,VSBb,VSBc,VRBa,VRBb,VRBc,VbaseB); % 21: Conversion Z(pu) a Z(Ohm) [SAm_ZZAB,SAm_ZZBC,SAm_ZZCA,SAm_ZZAG,SAm_ZZBG,SAm_ZZCG]=... conv6(SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG,zbaseA); [RAm_ZZAB,RAm_ZZBC,RAm_ZZCA,RAm_ZZAG,RAm_ZZBG,RAm_ZZCG]=... conv6(RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG,zbaseA); [SBp_ZZAB,SBp_ZZBC,SBp_ZZCA,SBp_ZZAG,SBp_ZZBG,SBp_ZZCG]=... conv6(SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG,zbaseB); [RBp_ZZAB,RBp_ZZBC,RBp_ZZCA,RBp_ZZAG,RBp_ZZBG,RBp_ZZCG]=... conv6(RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG,zbaseB); % 67N: Conversion pu a kV [VVSApol,VVRApol]=conv2(VSApol,VRApol,VbaseA); [VVSBpol,VVRBpol]=conv2(VSBpol,VRBpol,VbaseB); % 67N: Conversion pu a kA [IISAmop,IIRAmop]=conv2(ISAmop,IRAmop,IbaseA); % kA [IISBpop,IIRBpop]=conv2(ISBpop,IRBpop,IbaseB); alcSAmAB=abs(SAm_ZZAB)/cos(angle(SAm_ZZAB)-angle(zlA1u)); alcSAmBC=abs(SAm_ZZBC)/cos(angle(SAm_ZZBC)-angle(zlA1u)); alcSAmCA=abs(SAm_ZZCA)/cos(angle(SAm_ZZCA)-angle(zlA1u)); alcSAmAG=abs(SAm_ZZAG)/cos(angle(SAm_ZZAG)-angle(zlA1u)); alcSAmBG=abs(SAm_ZZBG)/cos(angle(SAm_ZZBG)-angle(zlA1u)); alcSAmCG=abs(SAm_ZZCG)/cos(angle(SAm_ZZCG)-angle(zlA1u)); alcSBpAB=abs(SBp_ZZAB)/cos(angle(SBp_ZZAB)-angle(zlB1u)); alcSBpBC=abs(SBp_ZZBC)/cos(angle(SBp_ZZBC)-angle(zlB1u)); alcSBpCA=abs(SBp_ZZCA)/cos(angle(SBp_ZZCA)-angle(zlB1u)); alcSBpAG=abs(SBp_ZZAG)/cos(angle(SBp_ZZAG)-angle(zlB1u)); alcSBpBG=abs(SBp_ZZBG)/cos(angle(SBp_ZZBG)-angle(zlB1u)); alcSBpCG=abs(SBp_ZZCG)/cos(angle(SBp_ZZCG)-angle(zlB1u)); alcRAmAB=abs(RAm_ZZAB)/cos(angle(RAm_ZZAB)-angle(zlA1u)); alcRAmBC=abs(RAm_ZZBC)/cos(angle(RAm_ZZBC)-angle(zlA1u)); alcRAmCA=abs(RAm_ZZCA)/cos(angle(RAm_ZZCA)-angle(zlA1u)); alcRAmAG=abs(RAm_ZZAG)/cos(angle(RAm_ZZAG)-angle(zlA1u)); alcRAmBG=abs(RAm_ZZBG)/cos(angle(RAm_ZZBG)-angle(zlA1u)); alcRAmCG=abs(RAm_ZZCG)/cos(angle(RAm_ZZCG)-angle(zlA1u)); alcRBpAB=abs(RBp_ZZAB)/cos(angle(RBp_ZZAB)-angle(zlB1u)); alcRBpBC=abs(RBp_ZZBC)/cos(angle(RBp_ZZBC)-angle(zlB1u)); alcRBpCA=abs(RBp_ZZCA)/cos(angle(RBp_ZZCA)-angle(zlB1u)); alcRBpAG=abs(RBp_ZZAG)/cos(angle(RBp_ZZAG)-angle(zlB1u)); alcRBpBG=abs(RBp_ZZBG)/cos(angle(RBp_ZZBG)-angle(zlB1u)); alcRBpCG=abs(RBp_ZZCG)/cos(angle(RBp_ZZCG)-angle(zlB1u)); % para Salida 1 alcSAmABpu=alcSAmAB/abs(zlA1u*m*L); alcSAmBCpu=alcSAmBC/abs(zlA1u*m*L); alcSAmCApu=alcSAmCA/abs(zlA1u*m*L); alcSAmAGpu=alcSAmAG/abs(zlA1u*m*L); alcSAmBGpu=alcSAmBG/abs(zlA1u*m*L); alcSAmCGpu=alcSAmCG/abs(zlA1u*m*L); alcSBpABpu=alcSBpAB/abs(zlB1u*m*L); alcSBpBCpu=alcSBpBC/abs(zlB1u*m*L); alcSBpCApu=alcSBpCA/abs(zlB1u*m*L); alcSBpAGpu=alcSBpAG/abs(zlB1u*m*L); alcSBpBGpu=alcSBpBG/abs(zlB1u*m*L); alcSBpCGpu=alcSBpCG/abs(zlB1u*m*L); alcRAmABpu=alcRAmAB/abs(zlA1u*(1-m)*L); alcRAmBCpu=alcRAmBC/abs(zlA1u*(1-m)*L); alcRAmCApu=alcRAmCA/abs(zlA1u*(1-m)*L); alcRAmAGpu=alcRAmAG/abs(zlA1u*(1-m)*L); alcRAmBGpu=alcRAmBG/abs(zlA1u*(1-m)*L); alcRAmCGpu=alcRAmCG/abs(zlA1u*(1-m)*L); alcRBpABpu=alcRBpAB/abs(zlB1u*(1-m)*L); alcRBpBCpu=alcRBpBC/abs(zlB1u*(1-m)*L); alcRBpCApu=alcRBpCA/abs(zlB1u*(1-m)*L); alcRBpAGpu=alcRBpAG/abs(zlB1u*(1-m)*L); alcRBpBGpu=alcRBpBG/abs(zlB1u*(1-m)*L); alcRBpCGpu=alcRBpCG/abs(zlB1u*(1-m)*L); cadylabel1='alcSAmpu'; cadylabel2='alcSBppu'; cadylabel3='alcRAmpu'; cadylabel4='alcRBppu'; % para Salida 2 alcSAmABpua=alcSAmAB/abs(zlA1u*L); alcSAmBCpua=alcSAmBC/abs(zlA1u*L); alcSAmCApua=alcSAmCA/abs(zlA1u*L); alcSAmAGpua=alcSAmAG/abs(zlA1u*L); alcSAmBGpua=alcSAmBG/abs(zlA1u*L); alcSAmCGpua=alcSAmCG/abs(zlA1u*L); alcSBpABpua=alcSBpAB/abs(zlB1u*L); alcSBpBCpua=alcSBpBC/abs(zlB1u*L); alcSBpCApua=alcSBpCA/abs(zlB1u*L); alcSBpAGpua=alcSBpAG/abs(zlB1u*L); alcSBpBGpua=alcSBpBG/abs(zlB1u*L); alcSBpCGpua=alcSBpCG/abs(zlB1u*L); alcRAmABpua=alcRAmAB/abs(zlA1u*L); alcRAmBCpua=alcRAmBC/abs(zlA1u*L); alcRAmCApua=alcRAmCA/abs(zlA1u*L); alcRAmAGpua=alcRAmAG/abs(zlA1u*L); alcRAmBGpua=alcRAmBG/abs(zlA1u*L); alcRAmCGpua=alcRAmCG/abs(zlA1u*L); alcRBpABpua=alcRBpAB/abs(zlB1u*L); alcRBpBCpua=alcRBpBC/abs(zlB1u*L); alcRBpCApua=alcRBpCA/abs(zlB1u*L); alcRBpAGpua=alcRBpAG/abs(zlB1u*L); alcRBpBGpua=alcRBpBG/abs(zlB1u*L); alcRBpCGpua=alcRBpCG/abs(zlB1u*L); cadylabel1a='alcSAmpua'; cadylabel2a='alcSBppua'; cadylabel3a='alcRAmpua'; cadylabel4a='alcRBppua'; Abscisa(fila,i)=m; OrdenadaSAmAB(fila,i)=alcSAmABpu; OrdenadaSAmBC(fila,i)=alcSAmBCpu; OrdenadaSAmCA(fila,i)=alcSAmCApu; OrdenadaSAmAG(fila,i)=alcSAmAGpu; OrdenadaSAmBG(fila,i)=alcSAmBGpu; OrdenadaSAmCG(fila,i)=alcSAmCGpu; OrdenadaSBpAB(fila,i)=alcSBpABpu; OrdenadaSBpBC(fila,i)=alcSBpBCpu; OrdenadaSBpCA(fila,i)=alcSBpCApu; OrdenadaSBpAG(fila,i)=alcSBpAGpu; OrdenadaSBpBG(fila,i)=alcSBpBGpu; OrdenadaSBpCG(fila,i)=alcSBpCGpu; OrdenadaRAmAB(fila,i)=alcRAmABpu; OrdenadaRAmBC(fila,i)=alcRAmBCpu; OrdenadaRAmCA(fila,i)=alcRAmCApu; OrdenadaRAmAG(fila,i)=alcRAmAGpu; OrdenadaRAmBG(fila,i)=alcRAmBGpu; OrdenadaRAmCG(fila,i)=alcRAmCGpu; OrdenadaRBpAB(fila,i)=alcRBpABpu; OrdenadaRBpBC(fila,i)=alcRBpBCpu; OrdenadaRBpCA(fila,i)=alcRBpCApu; OrdenadaRBpAG(fila,i)=alcRBpAGpu; OrdenadaRBpBG(fila,i)=alcRBpBGpu; OrdenadaRBpCG(fila,i)=alcRBpCGpu; OrdenadaSAmABa(fila,i)=alcSAmABpua; OrdenadaSAmBCa(fila,i)=alcSAmBCpua; OrdenadaSAmCAa(fila,i)=alcSAmCApua; OrdenadaSAmAGa(fila,i)=alcSAmAGpua; OrdenadaSAmBGa(fila,i)=alcSAmBGpua; OrdenadaSAmCGa(fila,i)=alcSAmCGpua; OrdenadaSBpABa(fila,i)=alcSBpABpua; OrdenadaSBpBCa(fila,i)=alcSBpBCpua; OrdenadaSBpCAa(fila,i)=alcSBpCApua; OrdenadaSBpAGa(fila,i)=alcSBpAGpua; OrdenadaSBpBGa(fila,i)=alcSBpBGpua; OrdenadaSBpCGa(fila,i)=alcSBpCGpua; OrdenadaRAmABa(fila,i)=alcRAmABpua; OrdenadaRAmBCa(fila,i)=alcRAmBCpua; OrdenadaRAmCAa(fila,i)=alcRAmCApua; OrdenadaRAmAGa(fila,i)=alcRAmAGpua; OrdenadaRAmBGa(fila,i)=alcRAmBGpua; OrdenadaRAmCGa(fila,i)=alcRAmCGpua; OrdenadaRBpABa(fila,i)=alcRBpABpua; OrdenadaRBpBCa(fila,i)=alcRBpBCpua; OrdenadaRBpCAa(fila,i)=alcRBpCApua; OrdenadaRBpAGa(fila,i)=alcRBpAGpua; OrdenadaRBpBGa(fila,i)=alcRBpBGpua; OrdenadaRBpCGa(fila,i)=alcRBpCGpua; i=i+1; m=m+0.01; end %(while m) fila=fila+1; k0=k0*10; end % FIN DE BUCLES k0, m %*************************************** if compens21==1 comp='21c'; elseif compens21==2 comp='21nc'; end %*************************************** % Titulo para las Salidas 1,2 %*************************************** cadenamc1=sprintf('PG128Mn,%s,%3.0fkV,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f\n',... comp,UA,UB,kE,L,Sccs,Sccr,Kr0); if imptransfer==1 cadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1)*zbaseA,Ksr0); cadenamc1=sprintf('%s,%s',cadenamc1,cadenamT); end cadenamc2=sprintf('%s,Rf=%-3.0f',cfalta,Rf); if falta==1|falta==2|falta==3|falta==4|falta==7|falta==8|falta==9|falta==13|falta==15 cadenamc=sprintf('%s,%s',cadenamc1,cadenamc2); elseif falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16 cadenamc3=sprintf('Rg=%-3.0f',Rg); cadenamc=sprintf('%s,%s,%s',cadenamc1,cadenamc2,cadenamc3); end % ---------------------------------------------------------------------------------- % SALIDA GRAFICA % ---------------------------------------------------------------------------------- variassalidas=1; while variassalidas==1 fprintf('\n'); fprintf('Salida Grafica 1: Analisis SAm, Elegir Unidad (pu)\n'); fprintf(' : SBp, " \n'); fprintf('Salida Grafica 2: Analisis SAm, Elegir Unidad (pua)\n'); fprintf(' : SBp, " \n'); fprintf('\n'); salida=input('salida: 1, 2 :'); switch(salida) case{1},% **salida 1** : alc...pu mm=0.01; cont=1; while mm<1; Abscisa0(cont)=mm; Ordenada0(cont)=1; Ordenadadi(cont)=0.8/mm; Ordenadadd(cont)=0.8/(1-mm); cont=cont+1; mm=mm+0.01; end extremo=input('extremo: SAm=1, SBp=2, (RAm=3), (RBp=4):'); switch(extremo) case{1},% salida 1 : alc...pu, **extremo Sm** unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:'); switch(unidad) case{1},% salida 1 : alc...pu, extremo SAm, **unidad AB** Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmAB(1:1,1:(i-1)); Ordenada12=OrdenadaSAmAB(2:2,1:(i-1)); Ordenada13=OrdenadaSAmAB(3:3,1:(i-1)); cadenaylabel='alcanceSAmABpu'; case{2},% salida 1 : alc...pu, extremo SAm, unidad BC Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmBC(1:1,1:(i-1)); Ordenada12=OrdenadaSAmBC(2:2,1:(i-1)); Ordenada13=OrdenadaSAmBC(3:3,1:(i-1)); cadenaylabel='alcanceSAmBCpu'; case{3},% salida 1 : alc...pu, extremo SAm, unidad CA Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmCA(1:1,1:(i-1)); Ordenada12=OrdenadaSAmCA(2:2,1:(i-1)); Ordenada13=OrdenadaSAmCA(3:3,1:(i-1)); cadenaylabel='alcanceSAmCApu'; case{4},% salida 1 : alc...pu, extremo SAm, unidad AG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmAG(1:1,1:(i-1)); Ordenada12=OrdenadaSAmAG(2:2,1:(i-1)); Ordenada13=OrdenadaSAmAG(3:3,1:(i-1)); cadenaylabel='alcanceSAmAGpu'; case{5},% salida 1 : alc...pu, extremo SAm, unidad BG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmBG(1:1,1:(i-1)); Ordenada12=OrdenadaSAmBG(2:2,1:(i-1)); Ordenada13=OrdenadaSAmBG(3:3,1:(i-1)); cadenaylabel='alcanceSAmBGpu'; case{6},% salida 1 : alc...pu, extremo SAm, unidad CG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmCG(1:1,1:(i-1)); Ordenada12=OrdenadaSAmCG(2:2,1:(i-1)); Ordenada13=OrdenadaSAmCG(3:3,1:(i-1)); cadenaylabel='alcanceSAmCGpu'; end %end del switch(unidad) case{2},% salida 1 : alc...pu, **extremo SBp** unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:'); switch(unidad) case{1},% salida 1 : alc...pu, extremo SBp, **unidad AB** Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpAB(1:1,1:(i-1)); Ordenada12=OrdenadaSBpAB(2:2,1:(i-1)); Ordenada13=OrdenadaSBpAB(3:3,1:(i-1)); cadenaylabel='alcanceSBpABpu'; case{2},% salida 1 : alc...pu, extremo SBp, unidad BC Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpBC(1:1,1:(i-1)); Ordenada12=OrdenadaSBpBC(2:2,1:(i-1)); Ordenada13=OrdenadaSBpBC(3:3,1:(i-1)); cadenaylabel='alcanceSBpBCpu'; case{3},% salida 1 : alc...pu, extremo SBp, unidad CA Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpCA(1:1,1:(i-1)); Ordenada12=OrdenadaSBpCA(2:2,1:(i-1)); Ordenada13=OrdenadaSBpCA(3:3,1:(i-1)); cadenaylabel='alcanceSBpCApu'; case{4},% salida 1 : alc...pu, extremo SBp, unidad AG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpAG(1:1,1:(i-1)); Ordenada12=OrdenadaSBpAG(2:2,1:(i-1)); Ordenada13=OrdenadaSBpAG(3:3,1:(i-1)); cadenaylabel='alcanceSBpAGpu'; case{5},% salida 1 : alc...pu, extremo SBp, unidad BG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpBG(1:1,1:(i-1)); Ordenada12=OrdenadaSBpBG(2:2,1:(i-1)); Ordenada13=OrdenadaSBpBG(3:3,1:(i-1)); cadenaylabel='alcanceSBpBGpu'; case{6},% salida 1 : alc...pu, extremo SBp, unidad CG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpCG(1:1,1:(i-1)); Ordenada12=OrdenadaSBpCG(2:2,1:(i-1)); Ordenada13=OrdenadaSBpCG(3:3,1:(i-1)); cadenaylabel='alcanceSBpCGpu'; end %end del switch(unidad) end %end del switch(extremo) ejemasy=input('Eje +y: '); ejemenosy=input('Eje -y (con signo): '); plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa,Ordenada11,'r',Abscisa,Ordenada12,'g',... Abscisa,Ordenada13,'b'),grid on axis([0 1 ejemenosy ejemasy]); title(cadenamc) xlabel('m=p'),ylabel(cadenaylabel) legend('Ref','Lim','k0=0.1','k0=1','k0=10') % ********************************************************** case{2},% **salida 2** : alc...pua extremo=input('extremo: SAm=1, SBp=2, (RAm=3), (RBp=4):'); switch(extremo) case{1},% salida 1 : alc...pu, **extremo Sm** mm=0.01; cont=1; while mm<1; Abscisa0(cont)=mm; Ordenada0(cont)=0.8; Ordenada45(cont)=mm; cont=cont+1; mm=mm+0.01; end unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:'); switch(unidad) case{1},% salida 2 : alc...pua, extremo SAm, **unidad AB** Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmABa(1:1,1:(i-1)); Ordenada12=OrdenadaSAmABa(2:2,1:(i-1)); Ordenada13=OrdenadaSAmABa(3:3,1:(i-1)); cadenaylabel='alcanceSAmABpua'; case{2},% salida 2 : alc...pua, extremo SAm, unidad BC Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmBCa(1:1,1:(i-1)); Ordenada12=OrdenadaSAmBCa(2:2,1:(i-1)); Ordenada13=OrdenadaSAmBCa(3:3,1:(i-1)); cadenaylabel='alcanceSAmBCpua'; case{3},% salida 2 : alc...pua, extremo SAm, unidad CA Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmCAa(1:1,1:(i-1)); Ordenada12=OrdenadaSAmCAa(2:2,1:(i-1)); Ordenada13=OrdenadaSAmCAa(3:3,1:(i-1)); cadenaylabel='alcanceSAmCApua'; case{4},% salida 2 : alc...pua, extremo SAm, unidad AG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmAGa(1:1,1:(i-1)); Ordenada12=OrdenadaSAmAGa(2:2,1:(i-1)); Ordenada13=OrdenadaSAmAGa(3:3,1:(i-1)); cadenaylabel='alcanceSAmAGpua'; case{5},% salida 2 : alc...pua, extremo SAm, unidad BG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmBGa(1:1,1:(i-1)); Ordenada12=OrdenadaSAmBGa(2:2,1:(i-1)); Ordenada13=OrdenadaSAmBGa(3:3,1:(i-1)); cadenaylabel='alcanceSAmBGpua'; case{6},% salida 2 : alc...pua, extremo SAm, unidad CG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSAmCGa(1:1,1:(i-1)); Ordenada12=OrdenadaSAmCGa(2:2,1:(i-1)); Ordenada13=OrdenadaSAmCGa(3:3,1:(i-1)); cadenaylabel='alcanceSAmCGpua'; end %end del switch(unidad) case{2},% salida 2 : alc...pua, **extremo SBp** mm=0.01; cont=1; while mm<1; Abscisa0(cont)=mm; Ordenada0(cont)=0.8; Ordenada45(cont)=mm; cont=cont+1; mm=mm+0.01; end unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:'); switch(unidad) case{1},% salida 2 : alc...pua, extremo SBp, **unidad AB** Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpABa(1:1,1:(i-1)); Ordenada12=OrdenadaSBpABa(2:2,1:(i-1)); Ordenada13=OrdenadaSBpABa(3:3,1:(i-1)); cadenaylabel='alcanceSBpABpua'; case{2},% salida 2 : alc...pua, extremo SBp, unidad BC Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpBCa(1:1,1:(i-1)); Ordenada12=OrdenadaSBpBCa(2:2,1:(i-1)); Ordenada13=OrdenadaSBpBCa(3:3,1:(i-1)); cadenaylabel='alcanceSBpBCpua'; case{3},% salida 2 : alc...pua, extremo SBp, unidad CA Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpCAa(1:1,1:(i-1)); Ordenada12=OrdenadaSBpCAa(2:2,1:(i-1)); Ordenada13=OrdenadaSBpCAa(3:3,1:(i-1)); cadenaylabel='alcanceSBpCApua'; case{4},% salida 2 : alc...pua, extremo SBp, unidad AG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpAGa(1:1,1:(i-1)); Ordenada12=OrdenadaSBpAGa(2:2,1:(i-1)); Ordenada13=OrdenadaSBpAGa(3:3,1:(i-1)); cadenaylabel='alcanceSBpAGpua'; case{5},% salida 2 : alc...pua, extremo SBp, unidad BG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpBGa(1:1,1:(i-1)); Ordenada12=OrdenadaSBpBGa(2:2,1:(i-1)); Ordenada13=OrdenadaSBpBGa(3:3,1:(i-1)); cadenaylabel='alcanceSBpBGpua'; case{6},% salida 2 : alc...pua, extremo SBp, unidad CG Abscisa=Abscisa(1:1,1:(i-1)); Ordenada11=OrdenadaSBpCGa(1:1,1:(i-1)); Ordenada12=OrdenadaSBpCGa(2:2,1:(i-1)); Ordenada13=OrdenadaSBpCGa(3:3,1:(i-1)); cadenaylabel='alcanceSBpCGpua'; end %end del switch(unidad) end %end del switch(extremo) ejemasy=input('Eje +y: '); ejemenosy=input('Eje -y (con signo): '); plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenada45,'m',Abscisa,Ordenada11,'r',Abscisa,Ordenada12,'g',... Abscisa,Ordenada13,'b'),grid on axis([0 1 ejemenosy ejemasy]); title(cadenamc) xlabel('m=p'),ylabel(cadenaylabel) legend('Cob','Ref','k0=0.1','k0=1','k0=10') end %end del switch(salida) variassalidas=input('Mas Salidas=1,Acabar=2 :'); end % end del while variassalidas variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :'); fprintf('\n'); end % end del while variasejecuciones fprintf('Pulse Enter\n'); e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.


ANÁLISIS DE ACTUACIÓN DE LAS PROTECCIONES 21 Y 67N ANTE FALTAS SIMPLES, MULTIPLES E INTERCIRCUITO EN LÍNEAS ELÉCTRICAS SOBRE LOS MISMOS APOYOS

: Ruiz Lozano, José Luis
: Ingeniería Industrial
Contenido del proyecto:
Directorio raíz  >  z_Prog  >  PGn  >  PG128Mn  >  FG128Mn.m