testecampo

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

%Cálculo do Perfil do Campo Elétrico
h=20; %Altura do perfil
x=[0:0.001:50]; %matriz com a distância horinzontal do perfil
hf=18.85; %altura da fase
df=7.5; %distância entre fases 
Q1=(2.476+1i*0.)*10^(-6); %definição das cargas (em C/m)
Q2=(-1.2704+1i*2.2)*10^(-6); 
Q3=(-1.2704+1i*-2.148)*10^(-6);
%cálculo dos ângulos
theta1=atan((x+df)/(hf-h));
thetap1=atan((x+df)/(hf+h));
theta2=atan((x)/(hf-h));
thetap2=atan((x)/(hf+h));
theta3=atan((x-df)/(hf-h));
thetap3=atan((x-df)/(hf+h));
%cálculo das distâncias
d1=sqrt((hf-h)^2+(df+x).^2);
dp1=sqrt((hf+h)^2+(df+x).^2);
d2=sqrt((hf-h)^2+(x).^2);
dp2=sqrt((hf+h)^2+(x).^2);
d3=sqrt((hf-h)^2+(x-df).^2);
dp3=sqrt((hf+h)^2+(x-df).^2);
%definição das imagens
Qp1=Q1; 
Qp2=Q2; 
Qp3=Q3; 
%cálculo da constante k para o cálculo do campo elétrico
k=2*pi*8.85*10^(-12);
%cálculo dos campos onde o índice p indica as imagens (Qprimo)
E1=Q1./(k*d1);
Ep1=Qp1./(k*dp1);
E2=Q2./(k*d2);
Ep2=Qp2./(k*dp2);
E3=Q3./(k*d3);
Ep3=Qp3./(k*dp3);
%cálculo das componentes dos campos
E1x=-E1.*sin(theta1);
E1y=E1.*cos(theta1);
Ep1x=Ep1.*sin(thetap1);
Ep1y=Ep1.*cos(thetap1);
E2x=-E2.*sin(theta2);
E2y=E2.*cos(theta2);
Ep2x=Ep2.*sin(thetap2);
Ep2y=Ep2.*cos(thetap2);
E3x=-E3.*sin(theta3);
E3y=E3.*cos(theta3);
Ep3x=Ep3.*sin(thetap3);
Ep3y=Ep3.*cos(thetap3);
%cálculo da resultante das componentes 
Ex=E1x+Ep1x+E2x+Ep2x+E3x+Ep3x;
Ey=E1y+Ep1y+E2y+Ep2y+E3y+Ep3y;
%cálculo do módulo do campo - o fator 0.01 é para converter em kV
E=0.001*sqrt((abs(Ex)).^2+(abs(Ey).^2));
plot(x,E,'r'); grid
xlabel('Distância da fase central em metros');
ylabel('Campo elétrico (KV/m)');
title('Perfil do campo elétrico para uma altura h=20m');