Tension Prefalla.docx

%% tension prefalla v5=1.0; Icarga=v5/(xz1); vpfalla=v5+Icarga*xt2; disp('***************TENSION PRE-FALLA (en la barra 3)*************') fprintf(' Vprefalla = %8.4f [pu] %8.4f° \n',abs(vpfalla),radtodeg(angle(vpfalla))) %% corrientes de secuencia a=-0.5+sqrt(3)/2*i; A=[1 1 1;1 a^2 a;1 a a^2]; Ia1=vpfalla/(ZBpo(4,4)+(ZBne(4,4)*(ZBce(2,2))/(ZBne(4,4)+ZBce(2,2)))); Ia0=-Ia1*(ZBne(4,4))/(ZBne(4,4)+ZBce(2,2)); Ia2=-Ia1*(ZBce(2,2))/(ZBne(4,4)+ZBce(2,2)); disp('*****************CORRIENTES DE SECUENCIA fprintf(' Ia0 = %8.4f [pu] %8.4f [deg] fprintf(' Ia1 = %8.4f [pu] %8.4f [deg] fprintf(' Ia2 = %8.4f [pu] %8.4f [deg]

(base fase a)**************') \n',abs(Ia0),radtodeg(angle(Ia0))) \n',abs(Ia1),radtodeg(angle(Ia1))) \n',abs(Ia2),radtodeg(angle(Ia2)))

disp('************CORRIENTES DE FASE **************') Iabc=A*[Ia0;Ia1;Ia2]; Ia=Iabc(1); Ib=Iabc(2); Ic=Iabc(3); fprintf(' fprintf(' fprintf('

Ia Ib Ic

= %8.4f [pu] = %8.4f [pu] = %8.4f [pu]

%8.4f [deg] \n',abs(Ia),radtodeg(angle(Ia))) %8.4f [deg] \n',abs(Ib),radtodeg(angle(Ib))) %8.4f [deg] \n',abs(Ic),radtodeg(angle(Ic)))

disp('**************CORRIENTES DE FALLA ***************') Ifalla=3*Ia0; fprintf(' IFALLA = %8.4f [pu] %8.4f [deg] \n',abs(Ifalla),radtodeg(angle(Ifalla)))

***************TENSION PRE-FALLA (en la barra 3)************* Vprefalla = 1.0374 [pu] 3.8692° *****************CORRIENTES DE SECUENCIA (base fase a)************** Ia0

= 1.2760 [pu] 101.5821 [deg]

Ia1

= 1.9557 [pu] -75.1974 [deg]

Ia2

= 0.6855 [pu] 110.8046 [deg]

************CORRIENTES DE FASE ************** Ia

= 0.0000 [pu] 0.0000 [deg]

Ib

= 2.8558 [pu] 154.4573 [deg]

Ic

= 3.1004 [pu] 54.3237 [deg]

**************CORRIENTES DE FALLA *************** IFALLA

= 3.8279 [pu] 101.5821 [deg]