Esfera rodante y método del ángulo de protección según IEC 62305

Fundamentos Físicos de las Descargas Atmosféricas y Modelo Electrogeométrico (EGM)

La protección contra descargas atmosféricas en infraestructuras complejas se basa en la comprensión rigurosa del fenómeno físico de la ruptura dieléctrica del aire a gran escala. La formación de un rayo se inicia mediante un proceso de ionización gradual conocido como trazador escalonado descendente (downward stepped leader), el cual transporta una densidad de carga considerable desde la nube hasta la superficie terrestre. A medida que la cabeza del trazador descendente se aproxima al suelo, el campo eléctrico macroscópico en las puntas de las estructuras elevadas excede la rigidez dieléctrica crítica del aire impuro o perturbado por la precipitación ( $Ecrit \approx 300 kV/m$ a $3 MV/m$ según las condiciones de humedad y presión local), desencadenando uno o varios trazadores ascendentes conectores (upward connecting leaders).

El punto en el espacio tridimensional donde el trazador descendente y el trazador ascendente se intersecan define la fase de salto final o final jump. En este instante exacto, la impedancia del canal se colapsa y se inicia la descarga de retorno (return stroke), caracterizada por un frente de onda de corriente extremadamente rápido. La distancia euclidiana entre la punta de la estructura que emite el trazador ascendente exitoso y la cabeza del trazador descendente en el momento del salto final se define formalmente como la distancia de impacto o distancia de choque ( $r_s$ ).

El Modelo Electrogeométrico (EGM, Electro-Geometric Model), derivado originalmente por Armstrong, Whitehead y Love, e integrado de forma probabilística en la serie normativa IEC 62305, establece una correlación no lineal entre la amplitud de pico de la corriente de choque ( $I_p$ ) y la distancia de impacto ( $r_s$ ). Físicamente, una mayor carga almacenada en el trazador descendente produce un campo eléctrico más severo a mayor distancia, incrementando la distancia a la cual la rigidez dieléctrica del espacio interelectródico sucumbe. La formulación matemática estándar adoptada para relacionar la distancia de impacto con la corriente de pico es:

r_s = A \cdot I_p^b

Donde, según los coeficientes empíricos validados operacionalmente por la IEC 62305 y la CIGRE para la modelación de estructuras terrestres:

r_s = 10 \cdot I_p^{0.65}

En esta relación, $r_s$ se expresa en metros ( $m$ ) e $I_p$ en kiloamperios ( $kA$ ). El fundamento del diseño del sistema de captación mediante la Esfera Rodante radica en el concepto de corriente mínima crítica ( $Imin$ ). Si un rayo posee una corriente inferior a $Imin$ , su distancia de impacto será menor que el radio de la esfera rodante ( $R$ ). Por lo tanto, para garantizar que una descarga con corriente equivalente a $Imin$ sea interceptada con certeza geométrica absoluta por los captadores y no por la estructura protegida, el radio de la esfera rodante debe fijarse rigurosamente en función de dicha corriente mínima:

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

Por consiguiente, un sistema dimensionado con un radio de esfera rodante menor ofrece un nivel de protección superior, puesto que es capaz de interceptar rayos con impulsos de corriente más pequeños (cuyas distancias de choque son más cortas) que de otro modo penetrarían el volumen de sombra electrogeométrica.

Parámetros Estadísticos e Impulsionales según IEC 62305-1

La norma IEC 62305-1 clasifica los niveles de protección contra rayos (LPL, Lightning Protection Levels) en cuatro categorías (LPL I a LPL IV). Cada nivel define un conjunto estandarizado de parámetros de la corriente de rayo, que representan los valores límite máximos para el diseño de la resistencia mecánica y térmica de los componentes, y los valores límite mínimos para la geometría de intersección de las captaciones.

La forma de onda de la corriente de retorno primaria se modela mediante la función analítica de double-exponential o la ecuación de Heidler, donde el impulso de primera descarga se caracteriza formalmente por una relación de tiempo de frente y tiempo de semi-valor de $10/350\ \mu s$ . La ecuación diferencial de la energía específica transportada por la corriente de rayo $i(t)$ a través de una unidad de resistencia eléctrica se expresa mediante la integral de acción:

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

Esta energía específica es la responsable del calentamiento adiabático instantáneo y de las fuerzas electrodinámicas en los conductores de bajada. La elevación de temperatura en un conductor de sección transversal $A$ sujeto a un impulso de rayo se calcula mediante la siguiente expresión derivada del balance 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]

Donde $\alpha$ es el coeficiente de temperatura de la resistencia eléctrica ( $K ^{-1}$ ), $\rho_0$ es la resistividad a temperatura ambiente ( $\Omega\cdot m$ ), y $C_v$ es la capacidad calorífica volumétrica del material ( $J /( m ^3\cdot K )$ ).

Parámetro Técnico / LPL	LPL I	LPL II	LPL III	LPL IV
Radio de la Esfera Rodante ( $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$
Corriente Máxima de Pico ( $Imax$ , $10/350\ \mu s$ )	$200 kA$	$150 kA$	$100 kA$	$100 kA$
Carga Total de la Descarga ( $Qflash$ )	$300 C$	$225 C$	$150 C$	$150 C$
Energía Específica ( $W/R$ )	$10.00 MJ /\Omega$	$5.62 MJ /\Omega$	$2.50 MJ /\Omega$	$2.50 MJ /\Omega$
Tamaño de Malla Reticular ( $W \times W$ )	$5 \times 5 m$	$10 \times 10 m$	$15 \times 15 m$	$20 \times 20 m$
Ángulo de Protección ( $\alpha$ ) para $h = 20 m$	$25^\circ$	$36^\circ$	$46^\circ$	$54^\circ$
Probabilidad de Intercepción ( $P_B$ )	$0.99$	$0.98$	$0.95$	$0.80$

Método de la Esfera Rodante (Rolling Sphere Method - RSM): Formulación y Geometría Tridimensional

El Método de la Esfera Rodante es el método analítico universal impuesto por la norma IEC 62305-3 para la determinación de las zonas protegidas en estructuras de geometría compleja, siendo independiente de la altura o de la inclinación de la superficie. Geométricamente, consiste en hacer rodar una esfera imaginaria de radio $R$ (asociado al LPL correspondiente) sobre y alrededor de toda la estructura a proteger, en todas las direcciones posibles, hasta que entre en contacto con el plano del suelo o con cualquier objeto puesto a tierra accesible.

Todos los puntos de la estructura que entren en contacto directo con la superficie de la esfera rodante se consideran expuestos a impactos directos de rayos y deben dotarse de elementos captadores (puntas de captación, conductores aéreos o mallas). El volumen por debajo de la superficie tangencial trazada por la esfera entre los puntos de captación queda completamente resguardado en la denominada zona de sombra electrogeométrica.

Penetración de la Esfera entre Captadores

Cuando la esfera rodante se apoya sobre dos captadores paralelos o aislados (como puntas de captación de altura equivalente $h_p$ o cables de guarda aéreos) separados por una distancia física $d$ , la esfera penetra en el espacio comprendido entre ellos. Para evitar que la esfera toque cualquier equipo o elemento de la estructura situado bajo los captadores, es necesario calcular analíticamente la profundidad de penetración ( $p$ ).

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

Esta fórmula exige intrínsecamente que la distancia entre captadores cumpla con la condición límite de continuidad geométrica $d \le 2R$ . Si $d > 2R$ , la esfera pasa entre los captadores e impacta directamente el plano horizontal inferior o la estructura a proteger.

La altura efectiva del envolvente de protección ( $h_z$ ) en el punto medio entre dos captadores de altura $h$ distanciados por $d$ viene expresada por:

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

Para arreglos tridimensionales formados por cuatro puntas captadoras dispuestas en los vértices de un rectángulo de lados $a$ y $b$ , la distancia diagonal representativa es $ddiag = \sqrt{a^2 + b^2}$ . La penetración tridimensional máxima en el centro geométrico del rectángulo se convierte en:

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

Condición de no penetración destructiva para un objeto de altura $hobj$ ubicado en el espacio intermedio:

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

Método del Ángulo de Protección (Protection Angle Method - PAM) y Método de la Malla (Mesh Method)

El Método del Ángulo de Protección (PAM) es una derivación simplificada y bidimensional del Modelo Electrogeométrico. Es geométricamente equivalente al Método de la Esfera Rodante únicamente para estructuras sencillas con límites horizontales simples y alturas restringidas que no excedan el radio de la esfera rodante correspondiente ( $h \le R$ ).

Variación No Lineal del Ángulo de Protección

En el método PAM, el volumen protegido por una punta de captación o conductor aéreo adopta la forma de un cono o prisma triangular cuyo vértice coincide con el extremo del captador. El ángulo de protección ( $\alpha$ ) no es constante; decrece de forma marcadamente no lineal con el incremento de la altura del captador sobre el plano de referencia ( $h$ ). Esta reducción matemática garantiza que el cono formado permanezca estrictamente inscrito dentro de la envolvente circular descrita por la esfera rodante de radio $R$ .

La ecuación que gobierna el ángulo de protección en el límite de tangencia estricta con la esfera rodante es:

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

El radio del círculo base de protección generado en el suelo ( $rprot$ ) por una punta de captación a una altura $h$ se determina analíticamente mediante:

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

Sustituyendo la dependencia implícita de $\alpha$ en función de la altura geométrica y el radio de la esfera rodante $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 expresión demuestra rigurosamente que cuando la altura del captador se aproxima al radio de la esfera rodante ( $h \to R$ ), el ángulo de protección tiende a cero ( $\alpha \to 0^\circ$ ), cancelando la aplicabilidad del método del ángulo de protección para alturas $h > R$ .

Método de la Malla (Mesh Method)

Para cubiertas planas de gran extensión o estructuras industriales donde la aplicación de puntas o cables aéreos resulta técnicamente inviable, la norma IEC 62305-3 prescribe el Método de la Malla. Una red reticular de conductores dispuestos sobre la superficie del techo proporciona una protección integral si se cumplen simultáneamente tres condiciones geométricas infranqueables:

Los conductores de captación se instalan en los bordes de la cubierta, en las salientes de la estructura y en las líneas de cumbrera.
Las dimensiones de las mallas (ancho de retícula $W$ ) no superan los valores límites definidos en la tabla normativa según el LPL (desde $5 \times 5 m$ para LPL I hasta $20 \times 20 m$ para LPL IV).
Ninguna instalación metálica sobre la cubierta sobresale fuera del volumen protegido por la red conductora sin un captador auxiliar.

Análisis Forense de Fallas en Sistemas de Captación y Dieléctricos Industriales

Las fallas en la protección contra descargas atmosféricas se originan por una subestimación del LPL, un modelado electrogeométrico defectuoso o un cálculo deficiente de las distancias de aislamiento dieléctrico. El análisis forense de estas fallas abarca tres mecanismos físicos principales:

Impactos Laterales en Estructuras de Gran Altura

En estructuras cuya altura excede el radio de la esfera rodante ( $h > R$ ), la esfera de radio $R$ puede tocar las paredes laterales verticales en puntos situados a una altura superior a $hlim = R$ . La corriente de las descargas en las paredes laterales suele ser menor, pero produce el desprendimiento de mampostería, la perforación de fachadas ventiladas de aluminio y la destrucción de sensores exteriores. El fallo forense reside en no extender el anillo de captación y los conductores de bajada a los bordes y salientes del tercio superior de la estructura en edificios de gran altura (donde $h > 60 m$ según IEC 62305-3).

Flameo Inverso (Back-Flashover) por Inductancia del Bajante

Cuando la corriente de rayo $i(t)$ fluye a través de un conductor de bajada, la caída de tensión total a lo largo del conductor respecto a la masa profunda de la tierra no depende únicamente de la resistencia de puesta a tierra ( $Rpat$ ), sino de la inductancia distribuida del bajante ( $Lbaj \approx 1\ \mu H/m$ ) debido al alto valor de la derivada temporal de la corriente ( $di/dt$ ):

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

Para un rayo LPL I de $200 kA$ con un tiempo de frente de $10\ \mu s$ , la derivada inicial promedia $\frac{di}{dt} = 20 kA/ \mu s = 2 \times 10^{10} A/s$ . En un bajante de $20 m$ de longitud ( $Lbaj = 20\ \mu H$ ), el componente inductivo genera un pico de tensión transitoria impulsional de:

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

Si la rigidez dieléctrica del aire ( $\sim 3 MV/m$ ) o de los aislantes sólidos que separan el bajante de las tuberías o cableados internos no es suficiente para soportar esta diferencia de potencial, se produce un arco eléctrico destructivo (flameo inverso o back-flashover), provocando incendios, perforación de tuberías de gas o la inyección directa de sobretensiones brutales en las redes de baja tensión.

Deformación Térmica y Fusión Dieléctrica del Conductor

El estrés térmico en el punto de impacto directo origina la fusión del metal por evaporación puntual. La masa de material fundido ( $mfund$ ) en el punto de entrada de la corriente continua secundaria o del pulso impulsional se relaciona con la carga total $Q$ mediante el potencial de caída catódica/anódica ( $U_{c,a} \approx 10 V$ ):

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

Donde $L_f$ es el calor latente de fusión y $C_p$ el calor específico del conductor. Un diseño forense deficiente omite el espesor mínimo de las chapas metálicas continuas (definido en IEC 62305-3 Tabla 3), provocando la perforación directa de tanques de almacenamiento de hidrocarburos.

Estrategias de Diseño Avanzado y Mitigación de Sobretensiones Transitorias

Para garantizar la integridad estructural y funcional de instalaciones complejas, el diseño del Sistema de Protección contra Rayos (SPCR / LPS) debe implementar topologías aisladas o no aisladas evaluadas rigurosamente mediante la distancia de separación ( $s$ ).

Cálculo de la Distancia de Separación Normativa

La distancia de separación $s$ representa la distancia de seguridad aire/material requerida entre el sistema de captación o de bajada y los componentes metálicos o líneas eléctricas internas, para prevenir el fenómeno de flameo inverso. La formulación general según la norma IEC 62305-3 es:

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

Donde:

$k_i$ : Coeficiente que depende del Nivel de Protección contra Rayos (LPL) seleccionado:
$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 aislamiento eléctrico del medio donde se evalúa la separación:
$k_m = 1.0 \quad ( Aire ), \quad k_m = 0.5 \quad ( Hormigón / Obra de fábrica ), \quad k_m = 0.7 \quad ( Aislante sintético / PRFV )$
$k_c$ : Coeficiente de división de la corriente de rayo entre los conductores de bajada. Para un arreglo espacial de $N$ bajantes interconnectados en un sistema no aislado:
$k_c = \frac{1}{2N} \quad o derivado de la matriz de impedancias de red.$
$l$ : Distancia geométrica lineal a lo largo del conductor de bajada o captador desde el punto donde se evalúa la distancia de separación hasta el punto de equipotencialidad más próximo (el punto de puesta a tierra).

Coordinación de Supresores de Sobretensión (SPD) en Zonas LPZ

La mitigación contra los efectos inducidos electromagnéticamente requiere la zonificación del volumen en Zonas de Protección contra Rayos (LPZ, Lightning Protection Zones) según la norma IEC 62305-4:

LPZ 0A: Zona expuesta a impactos directos y a campos electromagnéticos no atenuados.
LPZ 0B: Zona protegida contra impactos directos mediante el método de la esfera rodante, pero expuesta a campos electromagnéticos completos.
LPZ 1: Zona interior donde la corriente de choque está limitada por la división de corrientes y por la presencia de dispositivos supresores de sobretensión Tipo 1 (ensayados con onda $10/350\ \mu s$ ).
LPZ 2: Zona interior con apantallamiento adicional donde las sobretensiones impulsionales se reducen aún más mediante dispositivos SPD Tipo 2 (ensayados con onda $8/20\ \mu s$ ).

El acoplamiento inductivo en los bucles formados por los cables de señal y de energía dentro de la estructura se rige por la Ley de Faraday-Lenz. La tensión inducida en un bucle de área $Abucle$ sujeto a un campo magnético variable $B(t)$ acoplado a la corriente del bajante $i(t)$ es:

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 ende, la reducción del área del bucle ( $Abucle$ ) mediante el trazado paralelo trenzado o el apantallamiento coaxial blindado resulta tan indispensable como la instalación del propio supresor de sobretensión.

Integración Técnica y Simulación Electromagnética en Vexten Suite

En el ámbito de la ingeniería electromagnética moderna, el software Vexten Suite actúa como el motor central de análisis determinista y estocástico para el modelado en 3D de sistemas de captación contra descargas atmosféricas. Vexten Suite automatiza las metodologías computacionales especificadas por la norma IEC 62305 a través de una integración analítica entre la distribución espacial de la corriente y la respuesta en frecuencia de la red dieléctrica.

Algoritmo de Trazado de Rayos y Cálculo Tridimensional de la Esfera Rodante

El motor de cálculo geométrico de Vexten Suite proyecta numéricamente una superficie discreta hiper-esférica sobre la estructura importada desde formatos BIM/CAD (STEP, IFC). El sistema resuelve la interacción entre la malla de la estructura y la esfera rodante mediante ecuaciones polinómicas de intersección esfera-superficie:

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

Donde $\hat{\mathbf{n}}$ es el vector normal a la superficie de contacto en el espacio tridimensional. El programa realiza un barrido angular infinitesimal ( $d\theta, d\phi$ ) registrando todos los puntos de penetración donde la distancia entre la esfera y la envolvente de la estructura resulta ser menor que cero ( $dreal - R < 0$ ). Los puntos vulnerables son inmediatamente señalados vectorialmente con su respectivo Nivel de Protección contra Rayos (LPL).

Cálculo de Impedancia Impulsional e Integración con Redes de Puesta a Tierra

Vexten Suite integra simultáneamente el análisis de corrientes de cortocircuito en frecuencia industrial (IEC 60909) con el análisis transitorio impulsional de alta frecuencia de la descarga atmosférica (IEC 62305-3). El cálculo de la impedancia transitoria de la red de puesta a tierra ( $Z(t)$ ) no se limita al valor de resistencia en corriente continua ( $Rdc$ ), sino que contempla el comportamiento inductivo de las picas y conductores enterrados a través del modelado de líneas de transmisión distribuidas:

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\}

Donde $l_e$ es la longitud efectiva crítica del electrodo de puesta a tierra bajo impulsos de alta frecuencia. Si la longitud física del conductor enterrado excede $l_e$ , cualquier longitud adicional no contribuye a la reducción de la sobretensión inicial de la descarga. Vexten Suite calcula la longitud efectiva crítica mediante la relación empírica:

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

Donde $\rho_{\text{suelo}}$ representa la resistividad del terreno ( $\Omega\cdot m$ ) y $T_1$ el tiempo de frente del impulso de la corriente ( $\mu s$ ).

Flujo Integrado de Diseño y Validación en Vexten Suite

Evaluación Automática del Riesgo (IEC 62305-2): Determinación de la frecuencia de impactos anuales ( $N_d$ ) y cálculo de las pérdidas admisibles ( $R_1$ a $R_4$ ) para fijar computacionalmente el nivel de protección requerido (LPL I a IV).
Generación Geométrica de la Zona Protegida: Ejecución del módulo de la Esfera Rodante en 3D. El software determina la posición precisa de las puntas captadoras, altura de mástiles o trazado de redes de malla en techos para evitar puntos directos de colisión.
Evaluación de Distancias de Separación ( $s$ ): Mapeo nodal en 3D de la caída de tensión en cada bajante y verificación automática de la regla de separación aislante respecto a bandejas de cables, tuberías y ductos de proceso.
Coordinación Dinámica de SPDs: Generación de las especificaciones de supresores de sobretensión Tipo 1, Tipo 2 y Tipo 3 requeridos en cada frontera de zona LPZ, garantizando la coordinación de energía basada en la disipación del parámetro de energía específica ( $W/R$ ).

Gracias al rigor matemático del algoritmo de Vexten Suite, el diseñador obtiene un modelo electromagnético determinista que invalida la prueba y error, garantizando una protección integral contra los severos transitorios dieléctricos, mecánicos y térmicos provocados por las descargas atmosféricas según los más elevados estándares internacionales de la serie IEC 62305.