Análise de Redes de distribuição

...

V2 = complex(0,0);

V3_anterior = complex(1,0);

V3 = complex(0,0);

V4_anterior = complex(1,0);

V4 = complex(0,0);

V5_anterior = complex(1,0);

V5 = complex(0,0);

V6_anterior = complex(1,0);

V6 = complex(0,0);

V7_anterior = complex(1,0);

V7 = complex(0,0);

V8_anterior = complex(1,0);

V8 = complex(0,0);

erro=0.001;

while abs(V1 - V1_anterior) > erro

V1 = 1 - ((conj(S_SE_1))/conj(V1_anterior))*(Z_SE_1);

V1_anterior = V1;

end

while abs(V2 - V2_anterior) > erro

V2 = V1 - ((conj(S_1_2))/conj(V2_anterior))*(Z_1_2);

V2_anterior = V2;

end

while abs(V3 - V3_anterior) > erro

V3 = V2 - (conj((S_2_3)))/conj(V3_anterior)*(Z_2_3);

V3_anterior = V3;

end

while abs(V4 - V4_anterior) > erro

V4 = V3 - (conj((S_3_4)))/conj(V4_anterior)*(Z_3_4);

V4_anterior = V4;

end

while abs(V5 - V5_anterior) > erro

V5 = V4 - (conj((S_4_5)))/conj(V5_anterior)*(Z_4_5);

V5_anterior = V5;

end

while abs(V7 - V7_anterior) > erro

V7 = V2 - (conj((S_2_7)))/conj(V7_anterior)*(Z_2_7);

V7_anterior = V7;

end

while abs(V6 - V6_anterior) > erro

V6 = V5 - (conj((S_5_6)))/conj(V6_anterior)*(Z_5_6);

V6_anterior = V6;

end

while abs(V8 - V8_anterior) > erro

V8 = V5 - ((conj(S_5_8)))/conj(V8_anterior)*(Z_5_8);

V8_anterior = V8;

end

%Quedas de tensão no trecho

QuedaSE_1=(abs(V_SE-V1))*100; %Queda do trecho - SE ao ponto 1

Queda1_2=(abs(V1-V2))*100; %Queda do trecho - 1 ao ponto 2

Queda2_3=(abs(V2-V3))*100; %%Queda do trecho - 2 ao ponto 3

Queda3_4=(abs(V3-V4))*100; %Queda do trecho - 3 ao ponto 4

Queda4_5=(abs(V4-V5))*100; %Queda do trecho - 4 ao ponto 5

Queda5_6=(abs(V5-V6))*100; %Queda do trecho - 5 ao ponto 6

Queda2_7=(abs(V2-V7))*100; %Queda do trecho - 2 ao ponto 7

Queda5_8=(abs(V5-V8))*100; %Queda do trecho - 5 ao ponto 8

%Queda de tensão total

Queda_1=(abs(1-V1))*100;

Queda_2=(abs(1-V2))*100;

Queda_3=(abs(1-V3))*100;

Queda_4=(abs(1-V4))*100;

Queda_5=(abs(1-V5))*100;

Queda_6=(abs(1-V6))*100;

Queda_7=(abs(1-V7))*100;

Queda_8=(abs(1-V8))*100;

%Queda de tensão unitária

Quedauni_1=QuedaSE_1/(S_SE_1*L_SE_1);

Quedauni_2=Queda1_2/(S_1_2*L_1_2);

Quedauni_3=Queda2_3/(S_2_3*L_2_3);

Quedauni_4=Queda3_4/(S_3_4*L_3_4);

Quedauni_5=Queda4_5/(S_4_5*L_4_5);

Quedauni_6=Queda5_6/(S_5_6*L_5_6);

Quedauni_7=Queda2_7/(S_2_7*L_2_7);

Quedauni_8=Queda5_8/(S_5_8*L_5_8);

A partir desse script foi possível obter os resultados apresentados na Tabela 1.

Tabela 1 - Resultados do cálculo da queda de tensão pelo método exato

TRECHO

CARGA (MVA)

CONDUTOR

QUEDA DE TENSÃO (%)

Inicial

Final

Distância (km)

Acumulada

Fase

Unitária

Trecho

Total

SE

1

0,9

2,944

4/0

0,2686

0,7116

0,7116

1

2

1,7

2,576

4/0