Cálculo de resistência de aterramento na UFSC, disciplina de EEL7081 - Aspectos de Segurança em Engenharia Elétrica

projeto_aterramento.m

clc
close all
clearvars

% function handles
Rel = @(rho,Lh,Dh) rho./(2*pi*Lh).*log(400/2.54*Lh./Dh);
Rem = @(rho,Lh,Dem) 0.183*rho./Lh.* ...
    log10( ((sqrt(Lh.^2+Dem.^2) + Lh).^2 - Dem.^2)./ ...
    (Dem.^2 - (sqrt(Lh.^2+Dem.^2) - Lh).^2));

n = 6;      % numero de eletrodos
Lh = 3;     % comprimento das hastes [metros]
Dem = 3;    % distância horizontal entre eltrodos [metros]
Dh = 3/4;   % diâmetro da haste [polegadas]
rho = 600;  % resistividade aparente do solo [ohm * metro]

% distância horizontal em relação à origem
d_el = linspace(0,(n-1)*Dem,n)';

% matriz de distâncias entre as hastes i-j
D = zeros(n);
for i=1:n
    for j=1:n
        D(i,j) = abs(d_el(i) - d_el(j));
    end
end
D

% resistência de aterramento de um eletrodo vertical
r_el = Rel(rho,Lh,Dh);

% resistência da haste (i) por influencia da haste (j)
R = zeros(n);
for i=1:n
    for j=1:n
        if i == j
            R(i,j) = r_el;
        else
            R(i,j) = Rem(rho,Lh,D(i,j));
        end
    end
end
R

% resistência total em cada eletrodo
Re = sum(R)

% resistência equivalente total
Rne = 1/sum(1./Re)

% coeficiente de redução da resistência
K = Rne/r_el