Esfera rolante e método do ângulo de proteção segundo a IEC 62305

Fundamentos Físicos das Descargas Atmosféricas e Modelo Eletrogeométrico (EGM)

A proteção contra descargas atmosféricas em infraestruturas complexas baseia-se na compreensão rigorosa do fenômeno físico da ruptura dielétrica do ar em grande escala. A formação de um raio inicia-se mediante um processo de ionização gradual conhecido como líder escalonado descendente (downward stepped leader), o qual transporta uma densidade de carga considerável desde a nuvem até a superfície terrestre. À medida que a cabeça do líder descendente se aproxima do solo, o campo elétrico macroscópico nas pontas das estruturas elevadas excede a rigidez dielétrica crítica do ar impuro ou perturbado pela precipitação ( $Ecrit \approx 300 kV/m$ a $3 MV/m$ segundo as condições de umidade e pressão local), desencadeando um ou vários líderes ascendentes conectores (upward connecting leaders).

O ponto no espaço tridimensional onde o líder descendente e o líder ascendente se interseccionam define a fase de salto final ou final jump. Neste instante exato, a impedância do canal colapsa e inicia-se a descarga de retorno (return stroke), caracterizada por uma frente de onda de corrente extremamente rápida. A distância euclidiana entre a ponta da estrutura que emite o líder ascendente bem-sucedido e a cabeça do líder descendente no momento do salto final define-se formalmente como a distância de impacto ou distância de choque ( $r_s$ ).

O Modelo Eletrogeométrico (EGM, Electro-Geometric Model), derivado originalmente por Armstrong, Whitehead e Love, e integrado de forma probabilística na série normativa IEC 62305, estabelece uma correlação não linear entre a amplitude de pico da corrente de choque ( $I_p$ ) e a distância de impacto ( $r_s$ ). Físicamente, uma maior carga armazenada no líder descendente produz um campo elétrico mais severo a uma distância maior, incrementando a distância à qual a rigidez dielétrica do espaço intereletródico sucumbe. A formulação matemática padrão adotada para relacionar a distância de impacto com a corrente de pico é:

r_s = A \cdot I_p^b

Onde, segundo os coeficientes empíricos validados operacionalmente pela IEC 62305 e pelo CIGRE para a modelagem de estruturas terrestres:

r_s = 10 \cdot I_p^{0.65}

Nesta relação, $r_s$ expressa-se em metros ( $m$ ) e $I_p$ em quiloamperes ( $kA$ ). O fundamento do projeto do sistema de captação mediante a Esfera Rolante radica no conceito de corrente mínima crítica ( $Imin$ ). Se um raio possui uma corrente inferior a $Imin$ , sua distância de impacto será menor que o raio da esfera rolante ( $R$ ). Portanto, para garantir que uma descarga com corrente equivalente a $Imin$ seja interceptada com certeza geométrica absoluta pelos captores e não pela estrutura protegida, o raio da esfera rolante deve fixar-se rigorosamente em função daquela corrente mínima:

R = 10 \cdot (Imin)^{0.65} \implies Imin = \left( \frac{R}{10} \right)^{\frac{1}{0.65}}

Por conseguinte, um sistema dimensionado com um raio de esfera rolante menor oferece um nível de proteção superior, visto que é capaz de interceptar raios com impulsos de corrente menores (cujas distâncias de choque são mais curtas) que de outro modo penetrariam o volume de sombra eletrogeométrica.

Parâmetros Estatísticos e Impulsionais segundo a IEC 62305-1

A norma IEC 62305-1 classifica os níveis de proteção contra raios (LPL, Lightning Protection Levels) em quatro categorias (LPL I a LPL IV). Cada nível define um conjunto padronizado de parâmetros da corrente de raio, que representam os valores limite máximos para o projeto da resistência mecânica e térmica dos componentes, e os valores limite mínimos para a geometria de interseção das captações.

A forma de onda da corrente de retorno primária é modelada mediante a função analítica de exponencial dupla ou a equação de Heidler, onde o impulso de primeira descarga caracteriza-se formalmente por uma relação de tempo de frente e tempo de semivalor de $10/350\ \mu s$ . A equação diferencial da energia específica transportada pela corrente de raio $i(t)$ através de uma unidade de resistência elétrica expressa-se mediante a integral de ação:

\frac{W}{R} = \int0^{\infty} i^2(t) \, dt

Esta energia específica é a responsável pelo aquecimento adiabático instantâneo e pelas forças eletrodinâmicas nos condutores de descida. A elevação de temperatura em um condutor de seção transversal $A$ sujeito a um impulso de raio calcula-se mediante a seguinte expressão derivada do balanço térmico adiabático:

\Delta T = \frac{1}{\alpha} \left[ \exp \left( \frac{\alpha \cdot \rho_0}{C_v \cdot A^2} \int0^{\infty} i^2(t) \, dt \right) - 1 \right]

Onde $\alpha$ é o coeficiente de temperatura da resistência elétrica ( $K^{-1}$ ), $\rho_0$ é a resistividad à temperatura ambiente ( $\Omega\cdot m$ ), e $C_v$ é a capacidade calorífica volumétrica do material ( $J/(m^3\cdot K)$ ).

Parâmetro Técnico / LPL	LPL I	LPL II	LPL III	LPL IV
Raio da Esfera Rolante ( $R$ )	$20\ m$	$30\ m$	$45\ m$	$60\ m$
Corriente Mínima Interceptada ( $Imin$ )	$3.0\ kA$	$5.4\ kA$	$10.1\ kA$	$15.7\ kA$
Corrente Máxima de Pico ( $Imax$ , $10/350\ \mu s$ )	$200\ kA$	$150\ kA$	$100\ kA$	$100\ kA$
Carga Total da Descarga ( $Qflash$ )	$300\ C$	$225\ C$	$150\ C$	$150\ C$
Energia Específica ( $W/R$ )	$10.00\ MJ/\Omega$	$5.62\ MJ/\Omega$	$2.50\ MJ/\Omega$	$2.50\ MJ/\Omega$
Tamanho de Malha Reticular ( $W \times W$ )	$5 \times 5\ m$	$10 \times 10\ m$	$15 \times 15\ m$	$20 \times 20\ m$
Ângulo de Proteção ( $\alpha$ ) para $h = 20\ m$	$25^\circ$	$36^\circ$	$46^\circ$	$54^\circ$
Probabilidade de Intercepção ( $P_B$ )	$0.99$	$0.98$	$0.95$	$0.80$

Método da Esfera Rolante (Rolling Sphere Method - RSM): Formulação e Geometria Tridimensional

O Método da Esfera Rolante é o método analítico universal imposto pela norma IEC 62305-3 para a determinação das zonas protegidas em estruturas de geometria complexa, sendo independente da altura ou da inclinação da superfície. Geometricamente, consiste em fazer rolar uma esfera imaginária de raio $R$ (associado ao LPL correspondente) sobre e ao redor de toda a estrutura a proteger, em todas as direções possíveis, até que entre em contato com o plano do solo ou com qualquer objeto aterrado acessível.

Todos os pontos da estrutura que entrem em contato direto com a superfície da esfera rolante consideram-se expostos a impactos diretos de raios e devem dotar-se de elementos captores (pontas de captação, condutores aéreos ou malhas). O volume abaixo da superfície tangencial traçada pela esfera entre os pontos de captação fica completamente resguardado na denominada zona de sombra eletrogeométrica.

Penetração da Esfera entre Captores

Quando a esfera rolante apoia-se sobre dois captores paralelos ou isolados (como pontas de captação de altura equivalente $h_p$ ou cabos para-raios aéreos) separados por uma distância física $d$ , a esfera penetra no espaço compreendido entre eles. Para evitar que a esfera toque qualquer equipamento ou elemento da estrutura situado sob os captores, é necessário calcular analiticamente a profundidade de penetração ( $p$ ).

p = R - \sqrt{R^2 - \left( \frac{d}{2} \right)^2}

Esta fórmula exige intrinsecamente que a distância entre captores cumpra com a condição limite de continuidade geométrica $d \le 2R$ . Se $d > 2R$ , a esfera passa entre os captores e impacta diretamente o plano horizontal inferior ou a estrutura a proteger.

A altura efetiva do invólucro de proteção ( $h_z$ ) no ponto médio entre dois captores de altura $h$ distanciados por $d$ vem expressa por:

h_z = h - p = h - R + \sqrt{R^2 - \left( \frac{d}{2} \right)^2}

Para arranjos tridimensionais formados por quatro pontas captoras dispostas nos vértices de um retângulo de lados $a$ e $b$ , a distância diagonal representativa é $ddiag = \sqrt{a^2 + b^2}$ . A penetração tridimensional máxima no centro geométrico do retângulo converte-se em:

p3D = R - \sqrt{R^2 - \frac{a^2 + b^2}{4}}

Condição de não penetração destrutiva para um objeto de altura $hobj$ localizado no espaço intermediário:

hobj \le h - p3D \implies hobj \le h - R + \sqrt{R^2 - \frac{a^2 + b^2}{4}}

Método do Ângulo de Proteção (Protection Angle Method - PAM) e Método da Malha (Mesh Method)

O Método do Ângulo de Proteção (PAM) é uma derivação simplificada e bidimensional do Modelo Eletrogeométrico. É geometricamente equivalente ao Método da Esfera Rolante unicamente para estruturas singelas com limites horizontais simples e alturas restringidas que não excedam o raio da esfera rolante correspondente ( $h \le R$ ).

Variação Não Linear do Ângulo de Proteção

No método PAM, o volume protegido por uma ponta de captação ou condutor aéreo adota a forma de um cone ou prisma triangular cujo vértice coincide com a extremidade do captor. O ângulo de proteção ( $\alpha$ ) não é constante; decresce de forma acentuadamente não linear com o incremento da altura do captor sobre o plano de referência ( $h$ ). Esta redução matemática garante que o cone formado permaneça estritamente inscrito dentro do invólucro circular descrito pela esfera rolante de raio $R$ .

A equação que governa o ângulo de proteção no limite de tangência estrita com a esfera rolante é:

\alpha(h) = \arcsin \left( 1 - \frac{h}{R} \right)

O raio do círculo base de proteção gerado no solo ( $rprot$ ) por uma ponta de captação a uma altura $h$ determina-se analiticamente mediante:

rprot = h \cdot \tan[\alpha(h)]

Substituindo a dependência implícita de $\alpha$ em função da altura geométrica e do raio da esfera rolante $R$ :

rprot = h \cdot \tan \left[ \arcsin \left( 1 - \frac{h}{R} \right) \right] = h \cdot \frac{1 - \frac{h}{R}}{\sqrt{1 - \left(1 - \frac{h}{R}\right)^2}} = h \cdot \frac{R - h}{\sqrt{2Rh - h^2}}

Esta expressão demonstra rigorosamente que quando a altura do captor se aproxima do raio da esfera rolante ( $h \to R$ ), o ângulo de proteção tende a zero ( $\alpha \to 0^\circ$ ), cancelando a aplicabilidade do método do ângulo de proteção para alturas $h > R$ .

Método da Malha (Mesh Method)

Para coberturas planas de grande extensão ou estruturas industriais onde a aplicação de pontas ou cabos aéreos resulta tecnicamente inviável, a norma IEC 62305-3 prescreve o Método da Malha. Uma rede reticular de condutores dispostos sobre a superfície do telhado proporciona uma proteção integral se cumprirem-se simultaneamente três condições geométricas intransponíveis:

Os condutores de captação instalam-se nas bordas da cobertura, nas saliências da estrutura e nas linhas de cumeeira.
As dimensões das malhas (largura de retícula $W$ ) não superam os valores limite definidos na tabela normativa segundo o LPL (desde $5 \times 5\ m$ para LPL I até $20 \times 20\ m$ para LPL IV).
Nenhuma instalação metálica sobre a cobertura sobressai fora do volume protegido pela rede condutora sem um captor auxiliar.

Análise Forense de Falhas em Sistemas de Captação e Dielétricos Industriais

As falhas na proteção contra descargas atmosféricas originam-se por uma subestimação do LPL, uma modelagem eletrogeométrica defeituosa ou um cálculo deficiente das distâncias de isolamento dielétrico. A análise forense destas falhas abrange três mecanismos físicos principais:

Impactos Laterais em Estruturas de Grande Altura

Em estruturas cuja altura excede o raio da esfera rolante ( $h > R$ ), a esfera de raio $R$ pode tocar as paredes laterais verticais em pontos situados a uma altura superior a $hlim = R$ . A corrente das descargas nas paredes laterais costuma ser menor, mas produz o desprendimento de alvenaria, a perfuração de fachadas ventiladas de alumínio e a destruição de sensores externos. A falha forense reside em não estender o anel de captação e os condutores de descida às bordas e saliências do terço superior da estrutura em edifícios de grande altura (onde $h > 60\ m$ segundo a IEC 62305-3).

Surto Induzido por Inversão de Potencial (Back-Flashover) devido à Indutância da Descida

Quando a corrente de raio $i(t)$ flui através de um condutor de descida, a queda de tensão total ao longo do condutor em relação à massa profunda da terra não depende unicamente da resistência de aterramento ( $Rpat$ ), mas da indutância distribuída da descida ( $Lbaj \approx 1\ \mu H/m$ ) devido ao alto valor da derivada temporal da corrente ( $di/dt$ ):

V(t) = i(t) \cdot Rpat + Lbaj \cdot \frac{di(t)}{dt}

Para um raio LPL I de $200\ kA$ com um tempo de frente de $10\ \mu s$ , a derivada inicial média $\frac{di}{dt} = 20\ kA/\mu s = 2 \times 10^{10}\ A/s$ . Em uma descida de $20\ m$ de comprimento ( $Lbaj = 20\ \mu H$ ), o componente indutivo gera um pico de tensão transitória impulsional de:

Vind = 20 \times 10^{-6}\ H \cdot 2 \times 10^{10}\ A/s = 400\ kV

Se a rigidez dielétrica do ar ( $\sim 3\ MV/m$ ) ou dos isolantes sólidos que separam a descida das tubulações ou fiações internas não for suficiente para suportar esta diferença de potencial, produz-se um arco elétrico destrutivo (flameo inverso ou back-flashover), provocando incêndios, perfuração de tubulações de gás ou a injeção direta de sobretensões brutais nas redes de baixa tensão.

Deformação Térmica e Fusão Dielétrica do Condutor

O estresse térmico no ponto de impacto direto origina a fusão do metal por evaporação pontual. A massa de material fundido ( $mfund$ ) no ponto de entrada da corrente contínua secundária ou do pulso impulsional relaciona-se com a carga total $Q$ mediante o potencial de queda catódica/anódica ( $U_{c,a} \approx 10\ V$ ):

mfund = \frac{U_{c,a} \cdot Q}{L_f + C_p(Tfund - Tamb)}

Onde $L_f$ é o calor latente de fusão e $C_p$ o calor específico do condutor. Um projeto forense deficiente omite a espessura mínima das chapas metálicas contínuas (definida na IEC 62305-3 Tabela 3), provocando a perfuração direta de tanques de armazenamento de hidrocarbonetos.

Estratégias de Projeto Avançado e Mitigação de Sobretensões Transitórias

Para garantir a integridade estrutural e funcional de instalações complexas, o projeto do Sistema de Proteção contra Raios (SPCR / LPS) deve implementar topologias isoladas ou não isoladas avaliadas rigorosamente mediante a distância de separação ( $s$ ).

Cálculo da Distância de Separação Normativa

A distância de separação $s$ representa a distância de segurança ar/material requerida entre o sistema de captação ou de descida e os componentes metálicos ou linhas elétricas internas, para prevenir o fenômeno de back-flashover. A formulação geral segundo a norma IEC 62305-3 é:

s = k_i \cdot \frac{k_c}{k_m} \cdot l

Onde:

$k_i$ : Coeficiente que depende do Nivel de Proteção contra Raios (LPL) selecionado:
$k_i = 0.08 \quad (LPL\ I), \quad k_i = 0.06 \quad (LPL\ II), \quad k_i = 0.04 \quad (LPL\ III\ /\ IV)$
$k_m$ : Coeficiente de isolamento elétrico do meio onde se avalia a separação:
$k_m = 1.0 \quad (Ar), \quad k_m = 0.5 \quad (Concreto\ /\ Obra\ de\ alvenaria), \quad k_m = 0.7 \quad (Isolante\ sintético\ /\ PRFV)$
$k_c$ : Coeficiente de divisão da corrente de raio entre os condutores de descida. Para um arranjo espacial de $N$ descidas interconectadas em um sistema não isolado:
$k_c = \frac{1}{2N} \quad ou\ derivado\ da\ matriz\ de\ impedâncias\ de\ rede.$
$l$ : Distância geométrica linear ao longo do condutor de descida ou captador desde o ponto onde se avalia a distância de separação até o ponto de equipotencialidade mais próximo (o ponto de aterramento).

Coordenação de Dispositivos de Proteção contra Surtos (DPS) em Zonas LPZ

A mitigação contra os efeitos induzidos eletromagneticamente requer a zonificação do volume em Zonas de Proteção contra Raios (LPZ, Lightning Protection Zones) segundo a norma IEC 62305-4:

LPZ 0A: Zona expuesta a impactos diretos e a campos eletromagnéticos não atenuados.
LPZ 0B: Zona protegida contra impactos diretos mediante o método da esfera rolante, mas exposta a campos eletromagnéticos completos.
LPZ 1: Zona interior onde a corrente de choque está limitada pela divisão de correntes e pela presença de dispositivos de proteção contra surtos Tipo 1 (ensaiados com onda $10/350\ \mu s$ ).
LPZ 2: Zona interior com blindagem adicional onde as sobretensões impulsionais reduzem-se ainda mais mediante dispositivos DPS Tipo 2 (ensaiados com onda $8/20\ \mu s$ ).

O acoplamento indutivo nos loops formados pelos cabos de sinal e de energia dentro da estrutura rege-se pela Lei de Faraday-Lenz. A tensão induzida em um loop de área $Abucle$ sujeito a um campo magnético variável $B(t)$ acoplado à corrente da descida $i(t)$ é:

Vind(t) = -\frac{d\Phi}{dt} = -\frac{d}{dt} \int \mathbf{B}(t) \cdot d\mathbf{A} \approx \frac{\mu_0 \cdot Abucle}{2\pi \cdot d_p} \cdot \frac{di(t)}{dt}

Por conseguinte, a redução da área do loop ( $Abucle$ ) mediante o traçado paralelo trançado ou a blindagem coaxial resulta tão indispensável quanto a instalação do próprio dispositivo de proteção contra surtos.

Integração Técnica e Simulação Eletromagnética no Vexten Suite

No âmbito da engenharia eletromagnética moderna, o software Vexten Suite atua como o motor central de análise determinística e estocástica para a modelagem em 3D de sistemas de captação contra descargas atmosféricas. O Vexten Suite automatiza as metodologias computacionais especificadas pela norma IEC 62305 através de uma integração analítica entre a distribuição espacial da corrente e a resposta em frequência da rede dielétrica.

Algoritmo de Traçado de Raios e Cálculo Tridimensional da Esfera Rolante

O motor de cálculo geométrico do Vexten Suite projeta numericamente uma superfície discreta hiperesférica sobre a estrutura importada de formatos BIM/CAD (STEP, IFC). O sistema resolve a interação entre a malha da estrutura e a esfera rolante mediante equações polinomiais de interseção esfera-superfície:

\mathbf{P}_{centro} = \mathbf{P}_{captador} + R \cdot \hat{\mathbf{n}}

Onde $\hat{\mathbf{n}}$ é o vetor normal à superfície de contato no espaço tridimensional. O programa realiza uma varredura angular infinitesimal ( $d\theta, d\phi$ ) registrando todos os pontos de penetração onde a distância entre a esfera e o invólucro da estrutura resulta ser menor que zero ( $dreal - R < 0$ ). Os pontos vulneráveis são imediatamente assinalados vetorialmente com seu respectivo Nível de Proteção contra Raios (LPL).

Cálculo de Impedância Impulsional e Integração com Redes de Aterramento

O Vexten Suite integra simultaneamente a análise de correntes de curto-circuito em frequência industrial (IEC 60909) com a análise transitória impulsional de alta frequência da descarga atmosférica (IEC 62305-3). O cálculo da impedância transitória da rede de aterramento ( $Z(t)$ ) não se limita ao valor de resistência em corrente contínua ( $Rdc$ ), mas contempla o comportamento indutivo das hastes e condutores enterrados através da modelagem de linhas de transmissão distribuídas:

Z(t) = \mathcal{F}^{-1} \left\{ \sqrt{\frac{R + j\omega L}{G + j\omega C}} \cdot \coth\left( \sqrt{(R + j\omega L)(G + j\omega C)} \cdot l_e \right) \right\}

Onde $l_e$ é o comprimento efetivo crítico do eletrodo de aterramento sob impulsos de alta frequência. Se o comprimento físico do condutor enterrado excede $l_e$ , qualquer comprimento adicional não contribui para a redução da sobretensão inicial da descarga. O Vexten Suite calcula o comprimento efetivo crítico mediante a relação empírica:

l_e = 1.6 \cdot \sqrt{\rho_{\text{solo}} \cdot T_1}

Onde $\rho_{\text{solo}}$ representa a resistividade do terreno ( $\Omega\cdot m$ ) e $T_1$ o tempo de frente do impulso de corrente ( $\mu s$ ).

Fluxo Integrado de Projeto e Validação no Vexten Suite

Avaliação Automática do Risco (IEC 62305-2): Determinação da frequência de impactos anuais ( $N_d$ ) e cálculo das perdas admissíveis ( $R_1$ a $R_4$ ) para fixar computacionalmente o nível de proteção requerido (LPL I a IV).
Geração Geométrica da Zona Protegida: Execução do módulo da Esfera Rolante em 3D. O software determina a posição precisa das pontas captoras, altura de mastros ou traçado de redes de malha em telhados para evitar pontos diretos de colisão.
Avaliação de Distâncias de Separação ( $s$ ): Mapeamento nodal em 3D da queda de tensão em cada descida e verificação automática da regra de separação isolante em relação a eletrocalhas, tubulações e dutos de processo.
Coordenação Dinâmica de DPS: Geração das especificações de protetores de surto Tipo 1, Tipo 2 e Tipo 3 requeridos em cada fronteira de zona LPZ, garantindo a coordenação de energia baseada na dissipação do parâmetro de energia específica ( $W/R$ ).

Graças ao rigor matemático do algoritmo do Vexten Suite, o projetista obtém um modelo eletromagnético determinístico que invalida a tentativa e erro, garantindo uma proteção integral contra os severos transitórios dielétricos, mecânicos e térmicos provocados pelas descargas atmosféricas segundo os mais elevados padrões internacionais da série IEC 62305.