clear all
close all
S=30e6;
V=13.8e3;
X0=0.08;
X1=0.2;
X2=0.25;
Zaterramento_ohms=0.571320;
Zbase=((V.^2)/S);
Zn=3*(Zaterramento_ohms/Zbase);
Zparalelo_X2_X0=((j*(X0+Zn)*j*X2)/(j*(X0+Zn)+j*X2));
Ztotal=Zparalelo_X2_X0+(j*X1);
Ibase=(S/(sqrt(3)*V));
Ia1=j/Ztotal;
Ia2=(-Ia1*j*(X0+Zn))/(j*(X2+X0+Zn));
Ia0=(-Ia1*j*(X2))/(j*(X2+X0+Zn));
a=((-1/2)+j*(sqrt(3)/2));
a_quadrado=((-1/2)-j*(sqrt(3)/2));
Ib0=Ia0;
Ib1=Ia1*((-1/2)-j*(sqrt(3)/2));
Ib2=Ia2*((-1/2)+1i*(sqrt(3)/2));
Ic0=Ia0;
Ic1=Ia1*((-1/2)+j*(sqrt(3)/2));
Ic2=Ia2*((-1/2)-j*(sqrt(3)/2));
Ia0_A=Ibase*Ia0;
Ia1_A=Ibase*Ia1;
Ia2_A=Ibase*Ia2;
Ib0_B=Ibase*Ib0;
Ib1_B=Ibase*Ib1;
Ib2_B=Ibase*Ib2;
Ic0_C=Ibase*Ic0;
Ic1_C=Ibase*Ic1;
Ic2_C=Ibase*Ic2;
TDI=[1 1 1; 1 a_quadrado a; 1 a a_quadrado];
Matriz_Ia0_Ia1_Ia2=[Ia0; Ia1; Ia2];
Ia_Ib_Ic_pu=TDI*Matriz_Ia0_Ia1_Ia2;
IA=Ibase*Ia_Ib_Ic_pu(1,1);
IB=Ibase*Ia_Ib_Ic_pu(2,1);
IC=Ibase*Ia_Ib_Ic_pu(3,1);
Icc_pu=3*Ia0;
Icc=Ibase*Icc_pu;
Va0=-1*Ia2*j*X2
Va1=-1*Ia2*j*X2
Va2=-1*Ia2*j*X2
Va0_Va1_Va2=[Va0; Va1; Va2];
Va_Vb_Vc_pu=TDI*Va0_Va1_Va2
VA_fase=(V/(sqrt(3)))*Va_Vb_Vc_pu(1,1)
VB_fase=(V/(sqrt(3)))*Va_Vb_Vc_pu(2,1)
VC_fase=(V/(sqrt(3)))*Va_Vb_Vc_pu(3,1)
VA=(V)*Va_Vb_Vc_pu(1,1)
VB=(V)*Va_Vb_Vc_pu(2,1)
VC=(V)*Va_Vb_Vc_pu(3,1)
V_aterramento=Zn*Ia0*(V/(sqrt(3)))
origem =[0 0];
pIA = [real(IA) imag(IA)];
dp1 = pIA-origem;
pIB = [real(IB) imag(IB)];
dp2 = pIB-origem;
pIC = [real(IC) imag(IC)];
dp3 = pIC-origem;
pIa1 = [real(Ia1_A) imag(Ia1_A)];
dp4 = pIa1-origem;
pIa2 = [real(Ia2_A) imag(Ia2_A)];
dp5 = pIa2-origem;
pIa0 = [real(Ia0_A) imag(Ia0_A)];
origem_Ia0=dp4+dp5;
dp6 = pIa0-origem;
pIb1 = [real(Ib1_B) imag(Ib1_B)];
dp7 = pIb1-origem;
pIb2 = [real(Ib2_B) imag(Ib2_B)];
dp8 = pIb2-origem;
pIb0 = [real(Ib0_B) imag(Ib0_B)];
origem_Ib0=dp7+dp8;
dp9 = pIb0-origem;
pIc1 = [real(Ic1_C) imag(Ic1_C)];
dp10 = pIc1-origem;
pIc2 = [real(Ic2_C) imag(Ic2_C)];
dp11 = pIc2-origem;
pIc0 = [real(Ic0_C) imag(Ic0_C)];
origem_Ic0=dp10+dp11;
dp12 = pIc0-origem;
title('Fasores de corrente e suas sequências');
figure(1)
hold on
grid on
grid minor
quiver(origem(1),origem(2),dp1(1),dp1(2),0)
quiver(origem(1),origem(2),dp2(1),dp2(2),0)
quiver(origem(1),origem(2),dp3(1),dp3(2),0)
quiver(origem(1),origem(2),dp4(1),dp4(2),0)
quiver(dp4(1),dp4(2),dp5(1),dp5(2),0)
quiver(origem_Ia0(1),origem_Ia0(2),dp6(1),dp6(2),0)
quiver(origem(1),origem(2),dp7(1),dp7(2),0)
quiver(dp7(1),dp7(2),dp8(1),dp8(2),0)
quiver(origem_Ib0(1),origem_Ib0(2),dp9(1),dp9(2),0)
quiver(origem(1),origem(2),dp10(1),dp10(2),0)
quiver(pIc1(1),pIc1(2),dp11(1),dp11(2),0)
quiver(origem_Ic0(1),origem_Ic0(2),dp12(1),dp12(2),0)
legend('IA','IB','IC','Ia1','Ia2','Ia0','Ib1','Ib2','Ib0','Ic1','Ic2','Ic0')
legend show
Vb0=Va0;
Vb1=Va1*((-1/2)-j*(sqrt(3)/2));
Vb2=Va2*((-1/2)+j*(sqrt(3)/2));
Vc0=Va0;
Vc1=Va1*((-1/2)+j*(sqrt(3)/2));
Vc2=Va2*((-1/2)-j*(sqrt(3)/2));
Va0_A=(Va0*(V/sqrt(3)))
Va1_A=(Va1*(V/sqrt(3)))
Va2_A=(Va2*(V/sqrt(3)))
Vb0_B=(Vb0*(V/sqrt(3)))
Vb1_B=(Vb1*(V/sqrt(3)))
Vb2_B=(Vb2*(V/sqrt(3)))
Vc0_C=(Vc0*(V/sqrt(3)))
Vc1_C=(Vc1*(V/sqrt(3)))
Vc2_C=(Vc2*(V/sqrt(3)))
pVA = [real(VA_fase) imag(VA_fase)];
dp13 = pVA-origem;
pVB = [real(VB_fase) imag(VB_fase)];
dp14 = pVB-origem;
pVC = [real(VC_fase) imag(VC_fase)];
dp15 = pVC-origem;
pVa1 = [real(Va1_A) imag(Va1_A)];
dp16 = pVa1-origem;
pVa2 = [real(Va2_A) imag(Va2_A)];
dp17 = pVa2-origem;
pVa0 = [real(Va0_A) imag(Va0_A)];
origem_Va0=dp16+dp17;
dp18 = pVa0-origem;
pVb1 = [real(Vb1_B) imag(Vb1_B)];
dp19 = pVb1-origem;
pVb2 = [real(Vb2_B) imag(Vb2_B)];
dp20 = pVb2-origem;
pVb0 = [real(Vb0_B) imag(Vb0_B)];
origem_Vb0=dp19+dp20;
dp21 = pVb0-origem;
pVc1 = [real(Vc1_C) imag(Vc1_C)];
dp22 = pVc1-origem;
pVc2 = [real(Vc2_C) imag(Vc2_C)];
dp23 = pVc2-origem;
pVc0 = [real(Vc0_C) imag(Vc0_C)];
origem_Vc0=dp22+dp23;
dp24 = pVc0-origem;
figure(2)
hold on
title('Fasores de tensão de fase e suas sequências')
grid on
grid minor
quiver(origem(1),origem(2),dp13(1),dp13(2),0)
quiver(origem(1),origem(2),dp14(1),dp14(2),0)
quiver(origem(1),origem(2),dp15(1),dp15(2),0)
quiver(origem(1),origem(2),dp16(1),dp16(2),0)
quiver(dp16(1),dp16(2),dp17(1),dp17(2),0)
quiver(origem_Va0(1),origem_Va0(2),dp18(1),dp18(2),0)
quiver(origem(1),origem(2),dp19(1),dp19(2),0)
quiver(dp19(1),dp19(2),dp20(1),dp20(2),0)
quiver(origem_Vb0(1),origem_Vb0(2),dp21(1),dp21(2),0)
quiver(origem(1),origem(2),dp22(1),dp22(2),0)
quiver(dp22(1),dp22(2),dp23(1),dp23(2),0)
quiver(origem_Vc0(1),origem_Vc0(2),dp24(1),dp24(2),0)
hold off
legend('VA','VB','VC','Va1','Va2','Va0','Vb1','Vb2','Vb0','Vc1','Vc2','Vc0')
legend show