% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG52n\n');
fprintf('Esquema Linea Doble SR, Bifurcando a X,Y\n');
fprintf('5 Tramos\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Eleccion: kE,L1,L2,L3,L4,L5,Sccs,Sccx,Sccr,Sccy\n');
fprintf(' Kx0,Kr0,Ky0\n');
fprintf('*** Parametro k0=zs0/zr0=0,1-1-10 ***\n');
fprintf('Compensacion 21: Si,No\n');
fprintf('Falta Simple. Para cada Tramo: m=0...1\n');
fprintf('Barra m a elegir entre a,b,c,d,e,f,g,h,q,r\n');
fprintf('Nota: Barras Ficticias, por Linea y Tramo\n');
fprintf('Analisis: Extremos de la Linea en Falta\n');
fprintf('Salida Grafica 1: Extr1,m: Unidad Pref.(pu). k0=0.1-1-10\n');
fprintf(' : Extr2,m: Unidad Pref.(pu). k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr1,m: Unidad Pref.(pua).k0=0.1-1-10\n');
fprintf(' : Extr2,m: Unidad Pref.(pua).k0=0.1-1-10\n');
fprintf('\n');
fprintf('Nota: Extremo Primero(asimilado a Sm) y Segundo(asimilado a Xm) Circulares\n');
% ----------------------------------------------------------------------
%PG52n
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220;
kE=1;
L1=100;
L2=100;
L3=100;
L4=100;
L5=100;
Sccs=5000;
Sccx=5000;
Sccr=5000;
Sccy=5000;
Kx0=1;
Kr0=1;
Ky0=1;
compens21=2;
% ---------------------------------------------------------------
else
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L1=input('Longitud Tramo 1 km: ');
L2=input('Longitud Tramo 2 km: ');
L3=input('Longitud Tramo 3 km: ');
L4=input('Longitud Tramo 4 km: ');
L5=input('Longitud Tramo 5 km: ');
Sccs=input('Sccs(MVA): ');
Sccx=input('Sccx(MVA): ');
Sccr=input('Sccr(MVA): ');
Sccy=input('Sccy(MVA): ');
Kx0=input('Kx0: ');
Kr0=input('Kr0: ');
Ky0=input('Ky0: ');
compens21=input('Compensacion 21: Si=1, No=2 :');
end
% ---------------------------------------------------------------
fprintf('Punto en Falta: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15\n');
fprintf(' : T3inf=16,T5der=17,T5izq=18,T2inf=19,T1inf=20\n');
pfaltam=input('Punto en Falta: ');
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 :');
Rfm=input('Resistencia de Falta en m. Rf(Ohm) = ');
if faltam==7|faltam==8|faltam==9|faltam==11
Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
else
Rgm=0;
end
% ---------------------------------------------------------------
a=0.5;b=0.5;c=0.5;d=0.5;e=0.5;f=0.5;g=0.5;h=0.5;q=0.5;r=0.5;
% ---------------------------------------------------------------
En=Un/sqrt(3);
U=kE*Un;
E=kE*En;
if Un==400
zl1u=0.0346+0.3208i; zl2u=zl1u; zl0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif Un==220
zl1u=0.0829+0.4135i; zl2u=zl1u; zl0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif Un==132
zl1u=0.1342+0.4042i; zl2u=zl1u; zl0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif Un==66
zl1u=0.2156+0.3952i; zl2u=zl1u; zl0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i;
end
k=(zl0u-zl1u)/(3*zl1u);
zs1=(Un^2/Sccs)*i;
zs2=zs1;
zx1=(Un^2/Sccx)*i;
zx2=zx1;
zr1=(Un^2/Sccr)*i;
zr2=zr1;
zy1=(Un^2/Sccy)*i;
zy2=zy1;
zx0=Kx0*zx1; zr0=Kr0*zr1; zy0=Ky0*zy1;
Km=zlM0u/zl1u;
fprintf('\n');
[Minc]=subMinc(1);
fila=1;
k0=0.1;
while k0<11
zs0=k0*zr0;
i=1;
m=0.01;
while m<1
if pfaltam==11 a=m;elseif pfaltam==12 b=m;elseif pfaltam==13 c=m;elseif pfaltam==14 d=m;
elseif pfaltam==15 e=m;elseif pfaltam==16 f=m;elseif pfaltam==17 g=m;
elseif pfaltam==18 h=m;elseif pfaltam==19 q=m;elseif pfaltam==20 r=m;
end
[z0]=subz0(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r);
[z1]=subz1(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r);
[z2]=subz2(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r);
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc,z0);
[Y1]=subY1(Minc,z1);
[Y2]=subY2(Minc,z2);
% ZBUS(0,1,2)
Z0=inv(Y0); Z0(11,11);
Z1=inv(Y1); Z1(11,11);
Z2=inv(Y2); Z2(11,11);
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2]=mag2m(Z0,Z1,Z2,E,faltam,Rfm,Rgm,pfaltam);
% IBUS(0,1,2)
I0=zeros(20,1);
I1=zeros(20,1);
I2=zeros(20,1);
I0(pfaltam,1)=-Im0;
I1(pfaltam,1)=-Im1;
I2(pfaltam,1)=-Im2;
%if pfaltam==11
%I0=[0;0;0;0;0;0;0;0;0;0;-Im0;0;0;0;0;0;0;0;0;0];
%I1=[0;0;0;0;0;0;0;0;0;0;-Im1;0;0;0;0;0;0;0;0;0];
%I2=[0;0;0;0;0;0;0;0;0;0;-Im2;0;0;0;0;0;0;0;0;0];
%end
% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E;
V2=Z2*I2;
% I(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minc,V1,z1);
[Ie2]=subIe2(Minc,V2,z2);
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN EL PUNTO m DE FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
%
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA
Vs0=V0(1);Vs1=V1(1);Vs2=V2(1);
Vx0=V0(2);Vx1=V1(2);Vx2=V2(2);
Vr0=V0(3);Vr1=V1(3);Vr2=V2(3);
Vy0=V0(4);Vy1=V1(4);Vy2=V2(4);
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vxa,Vxb,Vxc]=abc(Vx0,Vx1,Vx2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Vya,Vyb,Vyc]=abc(Vy0,Vy1,Vy2);
Isx0=Ie0(5);Isx1=Ie1(5);Isx2=Ie2(5);
Isy0=Ie0(6);Isy1=Ie1(6);Isy2=Ie2(6);
Ixs0=Ie0(7);Ixs1=Ie1(7);Ixs2=Ie2(7);
Ixr0=Ie0(8);Ixr1=Ie1(8);Ixr2=Ie2(8);
Irx0=Ie0(9);Irx1=Ie1(9);Irx2=Ie2(9);
Iry0=Ie0(10);Iry1=Ie1(10);Iry2=Ie2(10);
Iys0=Ie0(11);Iys1=Ie1(11);Iys2=Ie2(11);
Iyr0=Ie0(12);Iyr1=Ie1(12);Iyr2=Ie2(12);
[Isxa,Isxb,Isxc]=abc(Isx0,Isx1,Isx2);
[Isya,Isyb,Isyc]=abc(Isy0,Isy1,Isy2);
[Ixsa,Ixsb,Ixsc]=abc(Ixs0,Ixs1,Ixs2);
[Ixra,Ixrb,Ixrc]=abc(Ixr0,Ixr1,Ixr2);
[Irxa,Irxb,Irxc]=abc(Irx0,Irx1,Irx2);
[Irya,Iryb,Iryc]=abc(Iry0,Iry1,Iry2);
[Iyra,Iyrb,Iyrc]=abc(Iyr0,Iyr1,Iyr2);
[Iysa,Iysb,Iysc]=abc(Iys0,Iys1,Iys2);
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k,Km,Isy0);
[Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k,Km,Isx0);
[Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=...
PDistComp(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k,Km,Ixr0);
[Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=...
PDistComp(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k,Km,Ixs0);
[Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k,Km,Iry0);
[Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k,Km,Irx0);
[Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=...
PDistComp(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k,Km,Iys0);
[Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=...
PDistComp(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k,Km,Iyr0);
case{2},
% NO COMPENSADA
[Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k);
[Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k);
[Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=...
PDist(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k);
[Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=...
PDist(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k);
[Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=...
PDist(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k);
[Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=...
PDist(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k);
[Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=...
PDist(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k);
[Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=...
PDist(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k);
end
switch(pfaltam)
case{11},
L=m*L1;
La=L1+L4;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,La);
L=L4+(1-m)*L1;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,La);
cadylabel1='alcSmpu';
cadylabel2='alcXmpu';
cadylabel1a='alcSmpua';
cadylabel2a='alcXmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=0;
bp=L1/La;
as=(L4+L1)/La;
bs=-L1/La;
% Salida pu
Lim1=0.8*(L1+L4)/(m*L1);
Lim2=0.8*(L1+L4)/(L4+(1-m)*L1);
case{12},
L=L1+m*L4;
La=L1+L4;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,La);
L=(1-m)*L4;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,La);
cadylabel1='alcSmpu';
cadylabel2='alcXmpu';
cadylabel1a='alcSmpua';
cadylabel2a='alcXmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L1/La;
bp=L4/La;
as=L4/La;
bs=-L4/La;
% Salida pu
Lim1=0.8*(L1+L4)/(L1+m*L4);
Lim2=0.8*(L1+L4)/((1-m)*L4);
case{13},
L=(1-m)*L4;
La=L4+L2+L3;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La);
L=L3+L2+m*L4;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
cadylabel1a='alcXmpua';
cadylabel2a='alcRmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L4/La;
bp=-L4/La;
as=(L3+L2)/La;
bs=L4/La;
% Salida pu
Lim1=0.8*(L4+L2+L3)/((1-m)*L4);
Lim2=0.8*(L4+L2+L3)/(L3+L2+m*L4);
case{14},
L=L4+m*L2;
La=L4+L2+L3;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La);
L=L3+(1-m)*L2;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
cadylabel1a='alcXmpua';
cadylabel2a='alcRmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L4/La;
bp=L2/La;
as=(L3+L2)/La;
bs=-L2/La;
Lim1=0.8*(L4+L2+L3)/(L4+m*L2);
Lim2=0.8*(L4+L2+L3)/(L3+(1-m)*L2);
case{15},
L=L4+L2+m*L3;
La=L4+L2+L3;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La);
L=(1-m)*L3;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
cadylabel1a='alcXmpua';
cadylabel2a='alcRmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=(L4+L2)/La;
bp=L3/La;
as=L3/La;
bs=-L3/La;
Lim1=0.8*(L4+L2+L3)/(L4+L2+m*L3);
Lim2=0.8*(L4+L2+L3)/((1-m)*L3);
case{16},
L=(1-m)*L3;
La=L3+L5;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,La);
L=L5+m*L3;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,La);
cadylabel1='alcRmpu';
cadylabel2='alcYmpu';
cadylabel1a='alcRmpua';
cadylabel2a='alcYmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L3/La;
bp=-L3/La;
as=L5/La;
bs=L3/La;
Lim1=0.8*(L3+L5)/((1-m)*L3);
Lim2=0.8*(L3+L5)/(L5+m*L3);
case{17},
L=L3+m*L5;
La=L3+L5;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,La);
L=(1-m)*L5;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,La);
cadylabel1='alcRmpu';
cadylabel2='alcYmpu';
cadylabel1a='alcRmpua';
cadylabel2a='alcYmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L3/La;
bp=L5/La;
as=L5/La;
bs=-L5/La;
Lim1=0.8*(L3+L5)/(L3+m*L5);
Lim2=0.8*(L3+L5)/((1-m)*L5);
case{18},
L=(1-m)*L5;
La=L5+L2+L1;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La);
L=L1+L2+m*L5;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La);
cadylabel1='alcYmpu';
cadylabel2='alcSmpu';
cadylabel1a='alcYmpua';
cadylabel2a='alcSmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L5/La;
bp=-L5/La;
as=(L1+L2)/La;
bs=L5/La;
Lim1=0.8*(L1+L2+L5)/((1-m)*L5);
Lim2=0.8*(L1+L2+L5)/(L1+L2+m*L5);
case{19},
L=L5+(1-m)*L2;
La=L5+L2+L1;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La);
L=L1+m*L2;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La);
cadylabel1='alcYmpu';
cadylabel2='alcSmpu';
cadylabel1a='alcYmpua';
cadylabel2a='alcSmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=(L5+L2)/La;
bp=-L2/La;
as=L1/La;
bs=L2/La;
Lim1=0.8*(L1+L2+L5)/(L5+(1-m)*L2);
Lim2=0.8*(L1+L2+L5)/(L1+m*L2);
case{20},
L=L5+L2+(1-m)*L1;
La=L5+L2+L1;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La);
L=m*L1;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La);
cadylabel1='alcYmpu';
cadylabel2='alcSmpu';
cadylabel1a='alcYmpua';
cadylabel2a='alcSmpua';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=(L5+L2+L1)/La;
bp=-L1/La;
as=0;
bs=L1/La;
Lim1=0.8*(L1+L2+L5)/(L5+L2+(1-m)*L1);
Lim2=0.8*(L1+L2+L5)/(m*L1);
end
Abscisa(fila,i)=m;
OrdenadaSmAB(fila,i)=alcSmABpu;
OrdenadaSmBC(fila,i)=alcSmBCpu;
OrdenadaSmCA(fila,i)=alcSmCApu;
OrdenadaSmAG(fila,i)=alcSmAGpu;
OrdenadaSmBG(fila,i)=alcSmBGpu;
OrdenadaSmCG(fila,i)=alcSmCGpu;
OrdenadaXmAB(fila,i)=alcXmABpu;
OrdenadaXmBC(fila,i)=alcXmBCpu;
OrdenadaXmCA(fila,i)=alcXmCApu;
OrdenadaXmAG(fila,i)=alcXmAGpu;
OrdenadaXmBG(fila,i)=alcXmBGpu;
OrdenadaXmCG(fila,i)=alcXmCGpu;
OrdenadaSmABa(fila,i)=alcSmABpua;
OrdenadaSmBCa(fila,i)=alcSmBCpua;
OrdenadaSmCAa(fila,i)=alcSmCApua;
OrdenadaSmAGa(fila,i)=alcSmAGpua;
OrdenadaSmBGa(fila,i)=alcSmBGpua;
OrdenadaSmCGa(fila,i)=alcSmCGpua;
OrdenadaXmABa(fila,i)=alcXmABpua;
OrdenadaXmBCa(fila,i)=alcXmBCpua;
OrdenadaXmCAa(fila,i)=alcXmCApua;
OrdenadaXmAGa(fila,i)=alcXmAGpua;
OrdenadaXmBGa(fila,i)=alcXmBGpua;
OrdenadaXmCGa(fila,i)=alcXmCGpua;
% Formar Vectores para LimIzq y LimDer
vectorLim1(1,i)=Lim1;
vectorLim2(1,i)=Lim2;
i=i+1;
m=m+0.01;
end %(while m)
fila=fila+1;
k0=k0*10;
end
if faltam==1
cfaltam='ag';
elemento='AG';
OrdenadaSm=OrdenadaSmAG;
OrdenadaXm=OrdenadaXmAG;
OrdenadaSma=OrdenadaSmAGa;
OrdenadaXma=OrdenadaXmAGa;
elseif faltam==2
cfaltam='bg';
elemento='BG';
OrdenadaSm=OrdenadaSmBG;
OrdenadaXm=OrdenadaXmBG;
OrdenadaSma=OrdenadaSmBGa;
OrdenadaXma=OrdenadaXmBGa;
elseif faltam==3
cfaltam='cg';
elemento='CG';
OrdenadaSm=OrdenadaSmCG;
OrdenadaXm=OrdenadaXmCG;
OrdenadaSma=OrdenadaSmCGa;
OrdenadaXma=OrdenadaXmCGa;
elseif faltam==4
cfaltam='ab';
elemento='AB';
OrdenadaSm=OrdenadaSmAB;
OrdenadaXm=OrdenadaXmAB;
OrdenadaSma=OrdenadaSmABa;
OrdenadaXma=OrdenadaXmABa;
elseif faltam==5
cfaltam='bc';
elemento='BC';
OrdenadaSm=OrdenadaSmBC;
OrdenadaXm=OrdenadaXmBC;
OrdenadaSma=OrdenadaSmBCa;
OrdenadaXma=OrdenadaXmBCa;
elseif faltam==6
cfaltam='ca';
elemento='CA';
OrdenadaSm=OrdenadaSmCA;
OrdenadaXm=OrdenadaXmCA;
OrdenadaSma=OrdenadaSmCAa;
OrdenadaXma=OrdenadaXmCAa;
elseif faltam==7
cfaltam='abg';
elemento='AB';
OrdenadaSm=OrdenadaSmAB;
OrdenadaXm=OrdenadaXmAB;
OrdenadaSma=OrdenadaSmABa;
OrdenadaXma=OrdenadaXmABa;
elseif faltam==8
cfaltam='bcg';
elemento='BC';
OrdenadaSm=OrdenadaSmBC;
OrdenadaXm=OrdenadaXmBC;
OrdenadaSma=OrdenadaSmBCa;
OrdenadaXma=OrdenadaXmBCa;
elseif faltam==9
cfaltam='cag';
elemento='CA';
OrdenadaSm=OrdenadaSmCA;
OrdenadaXm=OrdenadaXmCA;
OrdenadaSma=OrdenadaSmCAa;
OrdenadaXma=OrdenadaXmCAa;
elseif faltam==10
cfaltam='abc';
elemento='BC';
OrdenadaSm=OrdenadaSmBC;
OrdenadaXm=OrdenadaXmBC;
OrdenadaSma=OrdenadaSmBCa;
OrdenadaXma=OrdenadaXmBCa;
elseif faltam==11
cfaltam='abcg';
elemento='BC';
OrdenadaSm=OrdenadaSmBC;
OrdenadaXm=OrdenadaXmBC;
OrdenadaSma=OrdenadaSmBCa;
OrdenadaXma=OrdenadaXmBCa;
end
% ----------------------------------------------------------------------
% SALIDA GRAFICA
% ----------------------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
cadenam0=sprintf('PG52n,%s,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,',...
comp,Un,kE,L1,L2);
cadenam1=sprintf('L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm,Sccs=%-5.0f,\n',...
L3,L4,L5,Sccs);
cadenam1=sprintf('%s %s',cadenam0,cadenam1);
cadenam2=sprintf('Sccx=%-5.0f,Sccr=%-5.0f,Sccy=%-5.0f,',Sccx,Sccr,Sccy);
cadenam3=sprintf('Kx0=%-4.1f,Kr0=%-4.1f,Ky0=%-4.1f,T=%3.0f,%s,Rfm=%-3.0f,',...
Kx0,Kr0,Ky0,pfaltam,cfaltam,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
cadenam=sprintf('%s %s %s',cadenam1,cadenam2,cadenam3);
elseif faltam==7|faltam==8|faltam==9|faltam==11
cadenam4=sprintf('Rgm=%-3.0f,',Rgm);
cadenam=sprintf('%s %s %s %s',cadenam1,cadenam2,cadenam3,cadenam4);
end
cadenam5=sprintf('%s',elemento);
cadenam=sprintf('%s %s',cadenam,cadenam5);
% -----------------------------------------------------------------------
variassalidas=1;
while variassalidas==1
fprintf('Salida Grafica 1: Extr1,m: Unidad Pref.(pu). k0=0.1-1-10\n');
fprintf(' : Extr2,m: Unidad Pref.(pu). k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr1,m: Unidad Pref.(pua).k0=0.1-1-10\n');
fprintf(' : Extr2,m: Unidad Pref.(pua).k0=0.1-1-10\n');
salida=input('Salida: 1,2:');
switch(salida)
case{1},
Abscisa1m=Abscisa(1:1,1:(i-1));
Ordenada1m=OrdenadaSm(1:1,1:(i-1));
Abscisa2m=Abscisa(2:2,1:(i-1));
Ordenada2m=OrdenadaSm(2:2,1:(i-1));
Abscisa3m=Abscisa(3:3,1:(i-1));
Ordenada3m=OrdenadaSm(3:3,1:(i-1));
Abscisa4m=Abscisa(1:1,1:(i-1));
Ordenada4m=OrdenadaXm(1:1,1:(i-1));
Abscisa5m=Abscisa(2:2,1:(i-1));
Ordenada5m=OrdenadaXm(2:2,1:(i-1));
Abscisa6m=Abscisa(3:3,1:(i-1));
Ordenada6m=OrdenadaXm(3:3,1:(i-1));
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
cont=cont+1;
mm=mm+0.01;
end
ejemasy1=input('Subplot 1: Eje +y: ');
ejemenosy1=input('Subplot 1: Eje -y (con signo): ');
ejemasy2=input('Subplot 2: Eje +y: ');
ejemenosy2=input('Subplot 2: Eje -y (con signo): ');
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa3m,vectorLim1,'m',Abscisa1m,Ordenada1m,'r',Abscisa2m,Ordenada2m,'g',...
Abscisa3m,Ordenada3m,'b'),grid on
axis([0 1 ejemenosy1 ejemasy1]);
title(cadenam)
xlabel('m'),ylabel(cadylabel1)
%legend('Ref','Lim','k0=0.1','k0=1','k0=10')
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,vectorLim2,'m',Abscisa4m,Ordenada4m,'r',Abscisa5m,Ordenada5m,'g',...
Abscisa6m,Ordenada6m,'b'),grid on
axis([0 1 ejemenosy2 ejemasy2]);
%title(cadenam)
xlabel('m'),ylabel(cadylabel2)
legend('Ref','Lim','k0=0.1','k0=1','k0=10')
case(2)
Abscisa1m=Abscisa(1:1,1:(i-1));
Ordenada1ma=OrdenadaSma(1:1,1:(i-1));
Abscisa2m=Abscisa(2:2,1:(i-1));
Ordenada2ma=OrdenadaSma(2:2,1:(i-1));
Abscisa3m=Abscisa(3:3,1:(i-1));
Ordenada3ma=OrdenadaSma(3:3,1:(i-1));
Abscisa4m=Abscisa(1:1,1:(i-1));
Ordenada4ma=OrdenadaXma(1:1,1:(i-1));
Abscisa5m=Abscisa(2:2,1:(i-1));
Ordenada5ma=OrdenadaXma(2:2,1:(i-1));
Abscisa6m=Abscisa(3:3,1:(i-1));
Ordenada6ma=OrdenadaXma(3:3,1:(i-1));
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'m',Abscisa0,Ordenadacp,'c',...
Abscisa1m,Ordenada1ma,'r',Abscisa2m,Ordenada2ma,'g',...
Abscisa3m,Ordenada3ma,'b'),grid on
%axis([0 1 0 2]);
title(cadenam)
xlabel('m'),ylabel(cadylabel1a)
%legend('z1','Ref,Cob','k0=0.1','k0=1','k0=10')
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'m',Abscisa0,Ordenadacs,'c',...
Abscisa4m,Ordenada4ma,'r',Abscisa5m,Ordenada5ma,'g',...
Abscisa6m,Ordenada6ma,'b'),grid on
%axis([0 1 0 2]);
%title(cadenam)
xlabel('m'),ylabel(cadylabel2a)
legend('z1','Ref,Cob','k0=0.1','k0=1','k0=10')
end % switch Salida
variassalidas=input('Mas Salidas=1,Acabar=2 :');
end % end del while variassalidas
variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :');
end % end del while variasejecuciones
fprintf(1,'Pulse Enter\n');
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
