% Directorio: P24T % Archivo: ma2mp.m function[Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2,Iq0,Iq1,Iq2,Vq0,Vq1,Vq2]=... ma2mp(Z0,Z1,Z2,E,faltam,faltap,faltaq,Rfm,Rgm,Rfp,Rgp,Rfq,Rgq); Rgm=0;Rgp=0;Rgq=0; a=cos(2*pi/3)+sin(2*pi/3)*i; Mmp=zeros(12); %faltam=input('Falta en m: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :'); switch(faltam) case{1}, %Rfm=input('Resistencia de Falta en m. Rf(Ohm) = '); Mmp(1,1)=1;Mmp(1,2)=-1; Mmp(2,1)=1;Mmp(2,3)=-1; Mmp(3,1)=-3*Rfm;Mmp(3,4)=1;Mmp(3,5)=1;Mmp(3,6)=1; case{2}, %Rfm=input('Resistencia de Falta en m. Rf(Ohm): '); Mmp(1,1)=1;Mmp(1,2)=-a^2; Mmp(2,1)=1;Mmp(2,3)=-a; Mmp(3,1)=-Rfm;Mmp(3,2)=-a^2*Rfm;Mmp(3,3)=-a*Rfm;Mmp(3,4)=1;Mmp(3,5)=a^2;Mmp(3,6)=a; case{3}, %Rfm=input('Resistencia de Falta en m. Rf(Ohm): '); Mmp(1,1)=1;Mmp(1,2)=-a; Mmp(2,1)=1;Mmp(2,3)=-a^2; Mmp(3,1)=-Rfm;Mmp(3,2)=-a*Rfm;Mmp(3,3)=-a^2*Rfm;Mmp(3,4)=1;Mmp(3,5)=a;Mmp(3,6)=a^2; case{4}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): '); Mmp(1,1)=1; Mmp(2,2)=1;Mmp(2,3)=a; Mmp(3,2)=-Rfm;Mmp(3,3)=a*Rfm;Mmp(3,5)=1;Mmp(3,6)=-a; case{5}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): '); Mmp(1,1)=1; Mmp(2,2)=1;Mmp(2,3)=1; Mmp(3,2)=-Rfm;Mmp(3,3)=Rfm;Mmp(3,5)=1;Mmp(3,6)=-1; case{6}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): '); Mmp(1,1)=1; Mmp(2,2)=1;Mmp(2,3)=a^2; Mmp(3,2)=-Rfm;Mmp(3,3)=a^2*Rfm;Mmp(3,5)=1;Mmp(3,6)=-a^2; case{7}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): '); %Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): '); Mmp(1,1)=1;Mmp(1,2)=a;Mmp(1,3)=a^2; Mmp(2,2)=-a^2*Rfm;Mmp(2,3)=Rfm;Mmp(2,5)=a^2;Mmp(2,6)=-1; Mmp(3,1)=-(Rfm+3*Rgm);Mmp(3,2)=-Rfm;Mmp(3,3)=-Rfm; Mmp(3,4)=1;Mmp(3,5)=1;Mmp(3,6)=1; case{8}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): '); %Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): '); Mmp(1,1)=1;Mmp(1,2)=1;Mmp(1,3)=1; Mmp(2,2)=-Rfm;Mmp(2,3)=Rfm;Mmp(2,5)=1;Mmp(2,6)=-1; Mmp(3,1)=-(Rfm+3*Rgm);Mmp(3,2)=-a^2*Rfm;Mmp(3,3)=-a*Rfm; Mmp(3,4)=1;Mmp(3,5)=a^2;Mmp(3,6)=a; case{9}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): '); %Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): '); Mmp(1,1)=1;Mmp(1,2)=a^2;Mmp(1,3)=a; Mmp(2,2)=-Rfm;Mmp(2,3)=a^2*Rfm;Mmp(2,5)=1;Mmp(2,6)=-a^2; Mmp(3,1)=-(Rfm+3*Rgm);Mmp(3,2)=-a*Rfm;Mmp(3,3)=-a^2*Rfm; Mmp(3,4)=1;Mmp(3,5)=a;Mmp(3,6)=a^2; case{10}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = '); Mmp(1,1)=1; Mmp(2,2)=Rfm*a^2;Mmp(2,3)=-Rfm;Mmp(2,5)=-a^2;Mmp(2,6)=1; Mmp(3,2)=-Rfm;Mmp(3,3)=Rfm*a^2;Mmp(3,5)=1;Mmp(3,6)=-a^2; case{11}, %Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = '); %Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): '); Mmp(1,1)=-(Rfm+3*Rgm);Mmp(1,4)=1; Mmp(2,2)=-Rfm;Mmp(2,3)=-Rfm;Mmp(2,5)=1;Mmp(2,6)=1; Mmp(3,2)=-a*Rfm;Mmp(3,3)=-Rfm;Mmp(3,5)=a;Mmp(3,6)=1; end %faltap=input('Falta en p: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :'); switch(faltap) case{1}, %Rfp=input('Resistencia de Falta en p. Rf(Ohm): '); Mmp(4,7)=1;Mmp(4,8)=-1; Mmp(5,7)=1;Mmp(5,9)=-1; Mmp(6,7)=-3*Rfp;Mmp(6,10)=1;Mmp(6,11)=1;Mmp(6,12)=1; case{2}, %Rfp=input('Resistencia de Falta en p. Rf(Ohm): '); Mmp(4,7)=1;Mmp(4,8)=-a^2; Mmp(5,7)=1;Mmp(5,9)=-a; Mmp(6,7)=-Rfp;Mmp(6,8)=-a^2*Rfp;Mmp(6,9)=-a*Rfp;Mmp(6,10)=1;Mmp(6,11)=a^2;Mmp(6,12)=a; case{3}, %Rfp=input('Resistencia de Falta en p. Rf(Ohm): '); Mmp(4,7)=1;Mmp(4,8)=-a; Mmp(5,7)=1;Mmp(5,9)=-a^2; Mmp(6,7)=-Rfp;Mmp(6,8)=-a*Rfp;Mmp(6,9)=-a^2*Rfp;Mmp(6,10)=1;Mmp(6,11)=a;Mmp(6,12)=a^2; case{4}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Mmp(4,7)=1; Mmp(5,8)=1;Mmp(5,9)=a; Mmp(6,8)=-Rfp;Mmp(6,9)=a*Rfp;Mmp(6,11)=1;Mmp(6,12)=-a; case{5}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Mmp(4,7)=1; Mmp(5,8)=1;Mmp(5,9)=1; Mmp(6,8)=-Rfp;Mmp(6,9)=Rfp;Mmp(6,11)=1;Mmp(6,12)=-1; case{6}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Mmp(4,7)=1; Mmp(5,8)=1;Mmp(5,9)=a^2; Mmp(6,8)=-Rfp;Mmp(6,9)=a^2*Rfp;Mmp(6,11)=1;Mmp(6,12)=-a^2; case{7}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); %Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(4,7)=1;Mmp(4,8)=a;Mmp(4,9)=a^2; Mmp(5,8)=-a^2*Rfp;Mmp(5,9)=Rfp;Mmp(5,11)=a^2;Mmp(5,12)=-1; Mmp(6,7)=-(Rfp+3*Rgp);Mmp(6,8)=-Rfp;Mmp(6,9)=-Rfp; Mmp(6,10)=1;Mmp(6,11)=1;Mmp(6,12)=1; case{8}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); %Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(4,7)=1;Mmp(4,8)=1;Mmp(4,9)=1; Mmp(5,8)=-Rfp;Mmp(5,9)=Rfp;Mmp(5,11)=1;Mmp(5,12)=-1; Mmp(6,7)=-(Rfp+3*Rgp);Mmp(6,8)=-a^2*Rfp;Mmp(6,9)=-a*Rfp; Mmp(6,10)=1;Mmp(6,11)=a^2;Mmp(6,12)=a; case{9}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); %Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(4,7)=1;Mmp(4,8)=a^2;Mmp(4,9)=a; Mmp(5,8)=-Rfp;Mmp(5,9)=a^2*Rfp;Mmp(5,11)=1;Mmp(5,12)=-a^2; Mmp(6,7)=-(Rfp+3*Rgp);Mmp(6,8)=-a*Rfp;Mmp(6,9)=-a^2*Rfp; Mmp(6,10)=1;Mmp(6,11)=a;Mmp(6,12)=a^2; case{10}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm) = '); Mmp(4,7)=1; Mmp(5,8)=Rfp*a^2;Mmp(5,9)=-Rfp;Mmp(5,11)=-a^2;Mmp(5,12)=1; Mmp(6,8)=-Rfp;Mmp(6,9)=Rfp*a^2;Mmp(6,11)=1;Mmp(6,12)=-a^2; case{11}, %Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm) = '); %Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(4,7)=-(Rfp+3*Rgp);Mmp(4,10)=1; Mmp(5,8)=-Rfp;Mmp(5,9)=-Rfp;Mmp(5,11)=1;Mmp(5,12)=1; Mmp(6,8)=-a*Rfp;Mmp(6,9)=-Rfp;Mmp(6,11)=a;Mmp(6,12)=1; end %faltaq=input('Falta en q: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :'); switch(faltaq) case{1}, %Rfq=input('Resistencia de Falta en q. Rf(Ohm): '); Mmp(7,13)=1;Mmp(7,14)=-1; Mmp(8,13)=1;Mmp(8,15)=-1; Mmp(9,13)=-3*Rfq;Mmp(9,16)=1;Mmp(9,17)=1;Mmp(9,18)=1; case{2}, %Rfq=input('Resistencia de Falta en p. Rf(Ohm): '); Mmp(7,13)=1;Mmp(7,14)=-a^2; Mmp(8,13)=1;Mmp(8,15)=-a; Mmp(9,13)=-Rfq;Mmp(9,14)=-a^2*Rfq;Mmp(9,15)=-a*Rfq;Mmp(9,16)=1;Mmp(9,17)=a^2;Mmp(9,18)=a; case{3}, %Rfq=input('Resistencia de Falta en p. Rf(Ohm): '); Mmp(7,13)=1;Mmp(7,14)=-a; Mmp(8,13)=1;Mmp(8,15)=-a^2; Mmp(9,13)=-Rfq;Mmp(9,14)=-a*Rfq;Mmp(9,15)=-a^2*Rfq;Mmp(9,16)=1;Mmp(9,17)=a;Mmp(9,18)=a^2; case{4}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Mmp(7,13)=1; Mmp(8,14)=1;Mmp(8,15)=a; Mmp(9,14)=-Rfq;Mmp(9,15)=a*Rfq;Mmp(9,17)=1;Mmp(9,18)=-a; case{5}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Mmp(7,13)=1; Mmp(8,14)=1;Mmp(8,15)=1; Mmp(9,14)=-Rfq;Mmp(9,15)=Rfq;Mmp(9,17)=1;Mmp(9,18)=-1; case{6}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Mmp(7,13)=1; Mmp(8,14)=1;Mmp(8,15)=a^2; Mmp(9,14)=-Rfq;Mmp(9,15)=a^2*Rfq;Mmp(9,17)=1;Mmp(9,18)=-a^2; case{7}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); Rgq=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(7,13)=1;Mmp(7,14)=a;Mmp(7,15)=a^2; Mmp(8,14)=-a^2*Rfq;Mmp(8,15)=Rfq;Mmp(8,17)=a^2;Mmp(8,18)=-1; Mmp(9,13)=-(Rfq+3*Rgq);Mmp(9,14)=-Rfq;Mmp(9,15)=-Rfq; Mmp(9,16)=1;Mmp(9,17)=1;Mmp(9,18)=1; case{8}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); %Rgq=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(7,13)=1;Mmp(7,14)=1;Mmp(7,15)=1; Mmp(8,14)=-Rfq;Mmp(8,15)=Rfp;Mmp(8,17)=1;Mmp(8,18)=-1; Mmp(9,13)=-(Rfp+3*Rgp);Mmp(9,14)=-a^2*Rfp;Mmp(9,15)=-a*Rfp; Mmp(9,16)=1;Mmp(9,17)=a^2;Mmp(9,18)=a; case{9}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm): '); %Rgq=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(7,13)=1;Mmp(7,14)=a^2;Mmp(7,15)=a; Mmp(8,14)=-Rfp;Mmp(8,15)=a^2*Rfp;Mmp(8,17)=1;Mmp(8,18)=-a^2; Mmp(9,13)=-(Rfp+3*Rgp);Mmp(9,14)=-a*Rfp;Mmp(9,15)=-a^2*Rfp; Mmp(9,16)=1;Mmp(9,17)=a;Mmp(9,18)=a^2; case{10}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm) = '); Mmp(7,13)=1; Mmp(8,14)=Rfp*a^2;Mmp(8,15)=-Rfp;Mmp(8,17)=-a^2;Mmp(8,18)=1; Mmp(9,14)=-Rfp;Mmp(9,15)=Rfp*a^2;Mmp(9,17)=1;Mmp(9,18)=-a^2; case{11}, %Rfq=input('Resistencia de Falta,por fase, en p. Rf(Ohm) = '); %Rgq=input('Resistencia de Falta,comun en p. Rg(Ohm): '); Mmp(7,13)=-(Rfp+3*Rgp);Mmp(7,16)=1; Mmp(8,14)=-Rfp;Mmp(8,15)=-Rfp;Mmp(8,17)=1;Mmp(8,18)=1; Mmp(9,14)=-a*Rfp;Mmp(9,15)=-Rfp;Mmp(9,17)=a;Mmp(9,18)=1; end Mmp(10,1)=Z0(1,1); Mmp(10,4)=1; Mmp(10,7)=Z0(1,2); Mmp(10,13)=Z0(1,3); Mmp(11,1)=Z0(2,1); Mmp(11,7)=Z0(2,2); Mmp(11,10)=1; Mmp(11,13)=Z0(2,3); Mmp(12,1)=Z0(3,1); Mmp(12,7)=Z0(3,2); Mmp(12,13)=Z0(3,3); Mmp(12,16)=1; Mmp(13,2)=Z1(1,1); Mmp(13,5)=1; Mmp(13,8)=Z1(1,2); Mmp(13,14)=Z1(1,3); Mmp(14,2)=Z1(2,1); Mmp(14,8)=Z1(2,2); Mmp(14,11)=1; Mmp(14,14)=Z1(2,3); Mmp(15,2)=Z1(3,1); Mmp(15,8)=Z1(3,2); Mmp(15,14)=Z1(3,3); Mmp(15,17)=1; Mmp(16,3)=Z2(1,1); Mmp(16,6)=1; Mmp(16,9)=Z2(1,2); Mmp(16,15)=Z2(1,3); Mmp(17,3)=Z2(2,1); Mmp(17,9)=Z2(2,2); Mmp(17,12)=1; Mmp(17,15)=Z2(2,3); Mmp(18,3)=Z2(3,1); Mmp(18,9)=Z2(3,2); Mmp(18,15)=Z2(3,3); Mmp(18,18)=1; P=[0;0;0;0;0;0;0;0;0;0;0;0;E;E;E;0;0;0]; N=Mmp\P; % ------------------------------------------------------- % Im0=N(1);Im1=N(2);Im2=N(3); Vm0=N(4);Vm1=N(5);Vm2=N(6); Ip0=N(7);Ip1=N(8);Ip2=N(9); Vp0=N(10);Vp1=N(11);Vp2=N(12); Iq0=N(13);Iq1=N(14);Iq2=N(15); Vq0=N(16);Vq1=N(17);Vq2=N(18); 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  >  PGc  >  PG173ic  >  MA2MP.M