Publicidad

sábado, 3 de diciembre de 2016

Calculo de sistemas electricos de pontencia en MATLAB


Dejo en sus manos los cálculos en matlab mas complejos para realizar cálculos precisos sobre los sistemas eléctricos de potencia,trabajo de años de carrera en el instituto tecnológico de mexicali,disfrútenlo.

encontraras modelos en corta, media, y baja tension en modelos pi & T

CALCULO DE UNA LINEA CORTA
format long
% Analisis de la Linea CORTA
%***********************************
VR=215/sqrt(3)
Z=( 0.111971860750815+35.1506518824855i);
PR=100;
FP=0.9;
FR=sind(25.84);% AQUI VA EL COSENO INVERSO DEL FP
SR=PR/FP;
IR=SR/(sqrt(3)*215)*(FP-FR*i);
%***********************************
 

QR=SR*FR;
A=1;
B=Z;
C=0;
D=A;
T=[A,B;C,D]
R=[VR;IR]
S=T*R;%>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>VS y IS
S=T*R

%% CALCULOS CORRECTOS
VS=S(1,1);
MAGVS=abs(VS)
IS=S(2,1);
MAGIS=abs(IS)
ISCONJ=conj(IS);
Ss=3*VS*ISCONJ
SsMAG=abs(Ss)
PS=real(Ss)

QS=imag(Ss)

IC=IS-IR;

MAGIC=abs(IC);
IRMAG=abs(IR);
%PERDIDAS CU
Pcu=PS-PR
Zreal=real(Z);
PCU=(3*MAGIS^2*(Zreal/2))+(3*IRMAG^2*(Zreal/2))

QT=QS-QR;
Zimag=imag(Z);
QL=(3*MAGIS^2*(Zimag/2))+(3*IRMAG^2*(Zimag/2));
 
 

%PERDIDAS NETAS QRES=QL

% REGULACION
REG=((MAGVS-VR)/VR)*100
%EFICIENCIA
N=(PR/PS)*100
 

LINEA MEDIA EN MODELO T

format long
% Analisis de la Linea Media con modelo T
VR=66/sqrt(3)
% SI ES POR KILOMETRO SACARLO A PIE PARA LLEGAR A Zt

Z=(25+42i);
YC=j*(3.95*10^-4);
PR=20;
FP=.8;
FR=sind(25.84);
SR=20/FP;



QR=SR*FR;
IR=(SR/(sqrt(3)*66))*(FP-FR*i);


A=1+(Z*YC/2);
B=Z+(Z^2*YC/4);
C=YC;
D=A;
T=[A,B;C,D]


R=[VR;IR];
S=T*R;
S=T*R



%% CALCULOS CORRECTOS
VS=S(1,1)

MAGVS=abs(VS)

IS=S(2,1)

MAGIS=abs(IS)

ISCONJ=conj(IS)

Ss=3*VS*ISCONJ

PS=real(Ss)


QS=imag(Ss)


IC=IS-IR


MAGIC=abs(IC)

IRMAG=abs(IR)

%PERDIDAS CU
Pcu=PS-PR

PCU=37.5*(MAGIS^2+IRMAG^2)


QT=QS-QR

QL=63*(MAGIS^2+IRMAG^2)


QC=3*MAGIC^2*0.0004

QN=QL-QC

%PERDIDAS NETAS
QRES=QT+QC

% REGULACION
REG=((MAGVS-VR)/VR)*100

%EFICIENCIA
N=(PR/PS)*100



LINEA MEDIA EN MODELO PI


%%%% LINEA MEDIA PI
format 'long'
VR=66/sqrt(3)    %%kv
Z=(11.25+52.6202i);                 %%%impedancia de la linea entre fase
YC=j*(3.95*10^-4);                %%%%admitancia entre fase
FP=.8                     %%%factor en adelanto
OR=acos(FP)*180/pi         %%%grados de la corriente de resepcion IR 
FR=sind(OR)                %%%%%% factor de potencia de resepcion 
PR=20     %%%Mw
SR=PR/(FP)    %%%%MvA
IR=SR/((sqrt(3)*66))*(FP-FR*i)  %%%kA
IRMAG=abs(IR)
IRANGLE=angle(IR)*180/pi %%%saled negativo el angulo
QR=SR*tan(OR)   %%%%%MvAr



%%%%%%% formula de los parametros electricos 
A=(1+Z*YC/2)
B=Z
C=(1+Z*YC/4)*YC
D=A
T=[A,B;C,D]      %%%%parametros de electricos de la linea 
R=[VR;IR];         %%%%parametros de envio (R)
S=T*R;              %%%%%% parametros de resepcion (S)

VS=S(1,1);   %%%%% kv
VSMAG=abs(VS)
VSANGLE=angle(VS)*180/pi

IS=S(2,1);%%%%%%kA
ISMAG=abs(IS)
ISANGLE=angle(IS)*180/pi

Ss=3*VS*conj(IS) %%Mw
Ss2MAG=abs(Ss)
PS=real(Ss)      %%%%MvA
QS=imag(Ss)      %%%%%MvAr

IC1=VS*(YC/2);  %%%%%corriente capasitiva al principio de la linea 
IC1MAG=abs(IC1)

IC2=VR*YC/2;    %%%%%corriente capasitiva al final de la linea 
IC2MAG=abs(IC2)

Im=IS-IC1;        %%%%%%% corriente de la linea Kv
ImMAG=abs(Im)

Ir=Im-IC2;           %%%%%corriente de resepcion ******** ai dos de resepcion ????????????????? preguntar
IrMAG=abs(Ir)

VZ=VS/Im;              %%%%%%%%%% caida de tenciion de la linea Kv
VZMAG=abs(VZ)

Pcu=PS-PR              %%%%%%%%%%Mw
PCU=3*(IS^2)*(10);           %%%%%el 10 de la Z +
PCUMGA=abs(PCU)

Qt=QS-QR   %%%%%%%perdidas reactivas de la linea MvAr
QtMGA=abs(Qt)

JX=imag(Z);%%%%parte imaginaria de la linea Z
QL=3*(Im^2)*JX; %%%%%%%perdidas reactivas inductivas comprovacion MvAr 
QLMGA=abs(QL)

QC=3*(IC1^2+IC2^2)/(YC/2);    %%%%%%%perdidas reactivas capasitivas comprovacion MvAr 
QCMGA=abs(QC)

Qnet=QL-QC; %%%%%%%es la diferencia de reactancias
QnetMGA=abs(Qnet)
%***Calc Regulacion
Reg=((VS-VR)/VR)*100;
RegMGA=abs(Reg)

%Eficiencia
N=(PR/PS)*100






LINEA LARGA EN MODELO T

%LINEA LARGA MODELO T

VR=220/sqrt(3)
LONG=386.16
Z=(.174+.487i)*LONG
YC=j*.000003293*LONG
PR=75;
FP=1;
FR=0;
SR=PR/FP;
HZ=60;
L=386.16
IR=SR/((sqrt(3)*220)*(FP-FR*i))
LANDA=sqrt(Z*YC)
MAGLANDA=abs(LANDA)
ANGLANDA=angle(LANDA)*180/pi
Z0=sqrt(Z/YC)
MAGZ0=abs(Z0)
ANGZ0=angle(Z0)*180/pi
Y0=1/Z0
MAGY0=abs(Y0)
ANGY0=angle(Y0)*180/pi
% el angulo es el mismo de y0 y z0 
A=cosh(LANDA)
B=(Y0*sinh(LANDA))*(Z0)
C=Y0
D=A
T=[A,B;C,D]
R=[VR;IR]
S=T*R
VS=S(1,1)
MAGVS=abs(VS)
ANGVS=angle(VS)*180/pi
IS=S(2,1)
MAGIS=abs(IS)
ANGIS=angle(IS)*180/pi
% TRABAJANDO CON EL MODELO PI
ZE=B
RE=real(ZE)
XE=imag(ZE)
YA=((A-1)*2)/B
YE=0.001297504217382i
YEM=abs(YE)
IC1=VS*(YE/2)
MAGIC1=abs(IC1)
IC2=VR*(YE/2)
MAGIC2=abs(IC2)
IT=IS-IC1
ITM=abs(IT)
ITX=IR+IC2
% calculo de perdidas
CONJIS=conj(IS)
Ss=3*VS*CONJIS
PS=real(Ss)
QS=imag(Ss)
% la fuente recive reactivos del sistema, la liena tiene una compeansacion capasitiva muy fuerte, en lineas largas es normal
Pcu=PS-PR
PCU=3*RE*(ITM^2)
% perdidas rectivas
QL=3*XE*(ITM^2)
QC=6*(MAGIC1^2+MAGIC2^2)/YEM
QR=0
QT=QL-QS
% ES SUMA PERO QS NEGATIVO
REG=((MAGVS-VR)/VR)*100
N=(PR/PS)*100





 
LINEA LARGA EN MODELO PI


%LINEA LARGA MODELO PI 
VR=220/sqrt(3)
LONG=350
Z=(.174+.487i)*LONG
YC=j*.000003293*LONG
PR=75;
FP=1;
FR=0;
SR=PR/FP;
HZ=60;
L=350
IR=SR/((sqrt(3)*220)*(FP-FR*i))
MAGIR=abs(IR)
LANDA=sqrt(Z*YC)
MAGLANDA=abs(LANDA)
ANGLANDA=angle(LANDA)*180/pi
Z0=sqrt(Z/YC)
MAGZ0=abs(Z0)
ANGZ0=angle(Z0)*180/pi
Y0=1/Z0
MAGY0=abs(Y0)
ANGY0=angle(Y0)*180/pi
% el angulo es el mismo de y0 y z0 
A=cosh(LANDA)
B=Z0*sinh(LANDA)
C=Y0*sinh(LANDA)
D=A
T=[A,B;C,D]
R=[VR;IR]
S=T*R
VS=S(1,1)
MAGVS=abs(VS)
ANGVS=angle(VS)*180/pi
IS=S(2,1)
MAGIS=abs(IS)
ANGIS=angle(IS)*180/pi
% TRABAJANDO CON EL MODELO PI
ZE=B
RE=real(ZE)
XE=imag(ZE)
YA=((A-1)*2)/B
YE=0.001297504217382i
YEM=abs(YE)
IC1=VS*(YE/2)
MAGIC1=abs(IC1)
IC2=VR*(YE/2)
MAGIC2=abs(IC2)
IT=IS-IC1
ITM=abs(IT)
ITX=IR+IC2
% calculo de perdidas
CONJIS=conj(IS)
Ss=3*VS*CONJIS
PS=real(Ss)
QS=imag(Ss)
% la fuente recive reactivos del sistema, la liena tiene una compeansacion capasitiva muy fuerte, en lineas largas es normal
Pcu=PS-PR
PCU=3*RE*(ITM^2)
% perdidas rectivas
QL=3*XE*(ITM^2)
QC=6*(MAGIC1^2+MAGIC2^2)/YEM
QR=0
QT=QL-QS
% ES SUMA PERO QS NEGATIVO
REG=((MAGVS-VR)/VR)*100
N=(PR/PS)*100


Los archivos originales dispongo de ellos si los deseas puedes localizarme y te enviare una copia.

No hay comentarios.:

Publicar un comentario

Publicidad