Sfera rotolante e metodo dell'angolo di protezione secondo IEC 62305

Fondamenti Fisici delle Scariche Atmosferiche e Modello Elettrogeometrico (EGM)

La protezione contro le scariche atmosferiche nelle infrastrutture complesse si basa sulla comprensione rigorosa del fenomeno fisico della scarica dielettrica dell'aria su grande scala. La formazione di un fulmine viene avviata da un processo di ionizzazione graduale noto come leader a scatti discendente (downward stepped leader), il quale trasporta una notevole densità di carica dalla nube verso la superficie terrestre. Man mano che la testa del leader discendente si avvicina al suolo, il campo elettrico macroscopico sulle punte delle strutture elevate supera la rigidità dielettrica critica dell'aria impura o perturbata dalla precipitazione ( $Ecrit \approx 300 kV/m$ a $3 MV/m$ a seconda delle condizioni di umidità e pressione locale), innescando uno o più leader ascendenti di connessione (upward connecting leaders).

Il punto nello spazio tridimensionale in cui il leader discendente e il leader ascendente si intersecano definisce la fase del salto finale o final jump. In questo esatto istante, l'impedenza del canale crolla e ha inizio la scarica di ritorno (return stroke), caratterizzata da un fronte d'onda di corrente estremamente rapido. La distanza euclidea tra la punta della struttura che emette il leader ascendente efficace e la testa del leader discendente al momento del salto finale viene definita formalmente come distanza di impatto o distanza di scarica ( $r_s$ ).

Il Modello Elettrogeometrico (EGM, Electro-Geometric Model), derivato originariamente da Armstrong, Whitehead e Love, e integrato in forma probabilistica nella serie normativa IEC 62305, stabilisce una correlazione non lineare tra l'ampiezza di picco della corrente di scarica ( $I_p$ ) e la distanza di impatto ( $r_s$ ). Dal punto di vista fisico, una maggiore carica immagazzinata nel leader discendente produce un campo elettrico più severo a una distanza maggiore, incrementando la distanza alla quale la rigidità dielettrica dello spazio interelettrodico cede. La formulazione matematica standard adottata per mettere in relazione la distanza di impatto con la corrente di picco è:

r_s = A \cdot I_p^b

Dove, secondo i coefficienti empirici convalidati operativamente dalla IEC 62305 e dalla CIGRE per la modellazione di strutture terrestri:

r_s = 10 \cdot I_p^{0.65}

In questa relazione, $r_s$ è espressa in metri ( $m$ ) e $I_p$ in chiloampere ( $kA$ ). Il fondamento della progettazione del sistema di captazione mediante la Sfera Rotolante risiede nel concetto di corrente minima critica ( $Imin$ ). Se un fulmine possiede una corrente inferiore a $Imin$ , la sua distanza di impatto sarà minore del raggio della sfera rotolante ( $R$ ). Pertanto, per garantire che una scarica con corrente equivalente a $Imin$ sia intercettata con assoluta certezza geometrica dagli elementi captatori e non dalla struttura protetta, il raggio della sfera rotolante deve essere fissato rigorosamente in funzione di tale corrente minima:

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

Di conseguenza, un sistema dimensionato con un raggio della sfera rotolante inferiore offre un livello di protezione superiore, poiché è in grado di intercettare fulmini con impulsi di corrente più piccoli (le cui distanze di scarica sono più brevi) che altrimenti penetrerebbero nel volume di ombra elettrogeometrica.

Parametri Statistici e Impulsivi secondo IEC 62305-1

La norma IEC 62305-1 classifica i livelli di protezione contro i fulmini (LPL, Lightning Protection Levels) in quattro categorie (da LPL I a LPL IV). Ogni livello definisce un insieme standardizzato di parametri della corrente di fulmine, che rappresentano i valori limite massimi per la progettazione della resistenza meccanica e termica dei componenti, e i valori limite minimi per la geometria di intercettazione delle captazioni.

La forma d'onda della corrente di ritorno primaria viene modellata mediante la funzione analitica doppio-esponenziale o l'equazione di Heidler, in cui l'impulso della prima scarica è caratterizzato formalmente da un rapporto tra tempo di fronte e tempo di semivalore pari a $10/350\ \mu s$ . L'equazione differenziale dell'energia specifica trasportata dalla corrente di fulmine $i(t)$ attraverso un'unità di resistenza elettrica è espressa dall'integrale d'azione:

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

Questa energia specifica è responsabile del riscaldamento adiabatico istantaneo e delle forze elettrodinamiche nei conduttori di discesa. L'aumento di temperatura in un conduttore di sezione trasversale $A$ soggetto a un impulso di fulmine viene calcolato mediante la seguente espressione derivata dal bilancio termico adiabatico:

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

Dove $\alpha$ è il coefficiente di temperatura della resistenza elettrica ( $K ^{-1}$ ), $\rho_0$ è la resistività a temperatura ambiente ( $\Omega\cdot m$ ), e $C_v$ è la capacità calorifica volumetrica del materiale ( $J /( m ^3\cdot K )$ ).

Parametro Tecnico / LPL	LPL I	LPL II	LPL III	LPL IV
Raggio della Sfera Rotolante ( $R$ )	$20 m$	$30 m$	$45 m$	$60 m$
Corrente Minima Intercettata ( $Imin$ )	$3.0 kA$	$5.4 kA$	$10.1 kA$	$15.7 kA$
Corrente Massima di Picco ( $Imax$ , $10/350\ \mu s$ )	$200 kA$	$150 kA$	$100 kA$	$100 kA$
Carica Totale della Scarica ( $Qflash$ )	$300 C$	$225 C$	$150 C$	$150 C$
Energia Specifica ( $W/R$ )	$10.00 MJ /\Omega$	$5.62 MJ /\Omega$	$2.50 MJ /\Omega$	$2.50 MJ /\Omega$
Dimensione della Maglia ( $W \times W$ )	$5 \times 5 m$	$10 \times 10 m$	$15 \times 15 m$	$20 \times 20 m$
Angolo di Protezione ( $\alpha$ ) per $h = 20 m$	$25^\circ$	$36^\circ$	$46^\circ$	$54^\circ$
Probabilità di Intercettazione ( $P_B$ )	$0.99$	$0.98$	$0.95$	$0.80$

Metodo della Sfera Rotolante (Rolling Sphere Method - RSM): Formulazione e Geometria Tridimensionale

Il Metodo della Sfera Rotolante è il metodo analitico universale imposto dalla norma IEC 62305-3 per la determinazione delle zone protette in strutture a geometria complessa, risultando indipendente dall'altezza o dall'inclinazione della superficie. Dal punto di vista geometrico, consiste nel far rotolare una sfera immaginaria di raggio $R$ (associato al corrispondente LPL) sopra e attorno all'intera struttura da proteggere, in tutte le direzioni possibili, fino a quando non entra in contatto con il piano del suolo o con qualsiasi oggetto connesso a terra accessibile.

Tutti i punti della struttura che entrano in contatto diretto con la superficie della sfera rotolante sono considerati esposti a impatti diretti di fulmine e devono essere dotati di elementi captatori (aste di captazione, conduttori aerei o maglie). Il volume al di sotto della superficie tangenziale tracciata dalla sfera tra i punti di captazione rimane completamente protetto all'interno della cosiddetta zona di ombra elettrogeometrica.

Penetrazione della Sfera tra gli Elementi Captatori

Quando la sfera rotolante si appoggia su due captatori paralleli o isolati (come aste di captazione di altezza equivalente $h_p$ o cavi di guardia aerei) separati da una distanza fisica $d$ , la sfera penetra nello spazio compreso tra di essi. Per evitare che la sfera tocchi qualsiasi apparecchiatura o elemento della struttura situato sotto i captatori, è necessario calcolare analiticamente la profondità di penetrazione ( $p$ ).

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

Questa formula richiede intrinsecamente che la distanza tra i captatori soddisfi la condizione limite di continuità geometrica $d \le 2R$ . Se $d > 2R$ , la sfera passa tra i captatori e impatta direttamente il piano orizzontale inferiore o la struttura da proteggere.

L'altezza effettiva dell'inviluppo di protezione ( $h_z$ ) nel punto medio tra due captatori di altezza $h$ distanziati di $d$ è espressa da:

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

Per disposizioni tridimensionali formate da quattro aste captatrici disposte ai vertici di un rettangolo di lati $a$ e $b$ , la distanza diagonale rappresentativa è $ddiag = \sqrt{a^2 + b^2}$ . La penetrazione tridimensionale massima al centro geometrico del rettangolo diventa:

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

Condizione di non penetrazione distruttiva per un oggetto di altezza $hobj$ situato nello spazio intermedio:

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

Metodo dell'Angolo di Protezione (Protection Angle Method - PAM) e Metodo della Maglia (Mesh Method)

Il Metodo dell'Angolo di Protezione (PAM) è una derivazione semplificata e bidimensionale del Modello Elettrogeometrico. È geometricamente equivalente al Metodo della Sfera Rotolante unicamente per strutture semplici con confini orizzontali semplici e altezze limitate che non superano il raggio della sfera rotolante corrispondente ( $h \le R$ ).

Variazione Non Lineare dell'Angolo di Protezione

Nel metodo PAM, il volume protetto da un'asta di captazione o da un conduttore aereo assume la forma di un cono o di un prisma triangolare il cui vertice coincide con l'estremità del captatore. L'angolo di protezione ( $\alpha$ ) non è costante; decresce in modo marcatamente non lineare con l'aumento dell'altezza del captatore rispetto al piano di riferimento ( $h$ ). Questa riduzione matematica garantisce che il cono formato rimanga strettamente inscritto all'interno dell'inviluppo circolare descritto dalla sfera rotolante di raggio $R$ .

L'equazione che governa l'angolo di protezione al limite della tangenza stretta con la sfera rotolante è:

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

Il raggio del cerchio di base della protezione generato sul suolo ( $rprot$ ) da un'asta di captazione a un'altezza $h$ viene determinato analiticamente tramite:

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

Sostituendo la dipendenza implicita di $\alpha$ in funzione dell'altezza geometrica e del raggio della sfera rotolante $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}}

Questa espressione dimostra rigorosamente che quando l'altezza del captatore si avvicina al raggio della sfera rotolante ( $h \to R$ ), l'angolo di protezione tende a zero ( $\alpha \to 0^\circ$ ), annullando l'applicabilità del metodo dell'angolo di protezione per altezze $h > R$ .

Metodo della Maglia (Mesh Method)

Per coperture piane di grande estensione o strutture industriali dove l'applicazione di aste o cavi aerei risulta tecnicamente non fattibile, la norma IEC 62305-3 prescrive il Metodo della Maglia. Una rete reticolare di conduttori disposti sulla superficie del tetto fornisce una protezione integrale se sono soddisfatte simultaneamente tre condizioni geometriche vincolanti:

I conduttori di captazione sono installati sui bordi della copertura, sulle sporgenze della struttura e sulle linee di colmo.
Le dimensioni delle maglie (larghezza della griglia $W$ ) non superano i valori limite definiti nella tabella normativa in base al LPL (da $5 \times 5 m$ per LPL I fino a $20 \times 20 m$ per LPL IV).
Nessun impianto metallico sopra la copertura sporge al di fuori del volume protetto dalla rete conduttrice senza un captatore ausiliario.

Analisi Forense dei Guasti nei Sistemi di Captazione e Dielettrici Industriali

I guasti nella protezione contro le scariche atmosferiche sono originati da una sottostima del LPL, da un modellamento elettrogeometrico difettoso o da un calcolo carente delle distanze di isolamento dielettrico. L'analisi forense di tali guasti comprende tre meccanismi fisici principali:

Impatti Laterali su Strutture di Grande Altezza

Nelle strutture la cui altezza supera il raggio della sfera rotolante ( $h > R$ ), la sfera di raggio $R$ può toccare le pareti laterali verticali in punti situati a un'altezza superiore a $hlim = R$ . La corrente delle scariche sulle pareti laterali è solitamente inferiore, ma produce il distacco di murature, la perforazione di facciate ventilate in alluminio e la distruzione di sensori esterni. Il guasto forense risiede nel non aver esteso l'anello di captazione e i conduttori di calata ai bordi e alle sporgenze del terzo superiore della struttura negli edifici di grande altezza (dove $h > 60 m$ secondo la IEC 62305-3).

Scarica Inversa (Back-Flashover) per Induttanza della Calata

Quando la corrente di fulmine $i(t)$ scorre attraverso un conduttore di calata, la caduta di tensione totale lungo il conduttore rispetto alla massa profonda di terra non dipende unicamente dalla resistenza di messa a terra ( $Rpat$ ), ma dall'induttanza distribuita del discesore ( $Lbaj \approx 1\ \mu H/m$ ) a causa dell'elevato valore della derivata temporale della corrente ( $di/dt$ ):

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

Per un fulmine LPL I da $200 kA$ con un tempo di fronte di $10\ \mu s$ , la derivata iniziale media è $\frac{di}{dt} = 20 kA/ \mu s = 2 \times 10^{10} A/s$ . In una calata di $20 m$ di lunghezza ( $Lbaj = 20\ \mu H$ ), la componente induttiva genera un picco di tensione transitoria impulsiva pari a:

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

Se la rigidità dielettrica dell'aria ( $\sim 3 MV/m$ ) o degli isolanti solidi che separano la calata dalle tubazioni o dai cablaggi interni non è sufficiente a sopportare questa differenza di potenziale, si innesca un arco elettrico distruttivo (scarica inversa o back-flashover), provocando incendi, perforazione di tubazioni del gas o l'iniezione diretta di sovratensioni brutali nelle reti a bassa tensione.

Deformazione Termica e Fusione Dielettrica del Conduttore

Lo stress termico nel punto di impatto diretto origina la fusione del metallo per evaporazione puntuale. La massa di materiale fuso ( $mfund$ ) nel punto di ingresso della corrente continua secondaria o dell'impulso di corrente è legata alla carica totale $Q$ mediante il potenziale di caduta catodica/anodica ( $U_{c,a} \approx 10 V$ ):

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

Dove $L_f$ è il calore latente di fusione e $C_p$ il calore specifico del conduttore. Una progettazione forense insufficiente omette lo spessore minimo delle lamiere metalliche continue (definito nella Tabella 3 della IEC 62305-3), provocando la perforazione diretta di serbatoi di stoccaggio di idrocarburi.

Strategie di Progettazione Avanzata e Mitigazione delle Sovratensioni Transitorie

Per garantire l'integrità strutturale e funzionale di impianti complessi, la progettazione del Sistema di Protezione contro i Fulmini (LPS) deve implementare topologie isolate o non isolate valutate rigorosamente mediante la distanza di separazione ( $s$ ).

Calcolo della Distanza di Separación Normativa

La distanza di separazione $s$ rappresenta la distanza di sicurezza in aria/materiale richiesta tra il sistema di captazione o di calata e i componenti metallici o le linee elettriche interne, al fine di prevenire il fenomeno della scarica inversa. La formulazione generale secondo la norma IEC 62305-3 è:

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

Dove:

$k_i$ : Coefficiente che dipende dal Livello di Protezione contro i Fulmini (LPL) selezionato:
$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$ : Coefficiente di isolamento elettrico del mezzo in cui si valuta la separazione:
$k_m = 1.0 \quad ( Aria ), \quad k_m = 0.5 \quad ( Calcestruzzo / Muratura ), \quad k_m = 0.7 \quad ( Isolante sintetico / VTR )$
$k_c$ : Coefficiente di ripartizione della corrente di fulmine tra i conduttori di calata. Per una disposizione spaziale di $N$ calate interconnesse in un sistema non isolato:
$k_c = \frac{1}{2N} \quad o derivato dalla matrice delle impedenze di rete.$
$l$ : Distanza geometrica lineare lungo il conduttore di calata o il captatore dal punto in cui si valuta la distanza di separazione fino al punto di equipotenzialità più vicino (il punto di messa a terra).

Coordinamento degli Scaricatori di Sovratensione (SPD) nelle Zone LPZ

La mitigazione degli effetti indotti elettromagneticamente richiede la suddivisione del volume in Zone di Protezione contro i Fulmini (LPZ, Lightning Protection Zones) secondo la norma IEC 62305-4:

LPZ 0A: Zona esposta a impatti diretti e a campi elettromagnetici non attenuati.
LPZ 0B: Zona protetta contro gli impatti diretti tramite il metodo della sfera rotolante, ma esposta a campi elettromagnetici completi.
LPZ 1: Zona interna in cui la corrente d'impulso è limitata dalla ripartizione delle correnti e dalla presenza di dispositivi scaricatori di sovratensione di Tipo 1 (testati con onda $10/350\ \mu s$ ).
LPZ 2: Zona interna con schermatura aggiuntiva in cui le sovratensioni impulsive sono ulteriormente ridotte tramite dispositivi SPD di Tipo 2 (testati con onda $8/20\ \mu s$ ).

L'accoppiamento induttivo nei spiri formati dai cavi di segnale e di potenza all'interno della struttura è regolato dalla legge di Faraday-Lenz. La tensione indotta in una spira di area $Abucle$ soggetta a un campo magnetico variabile $B(t)$ accoppiato alla corrente del discesore $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}

Di conseguenza, la riduzione dell'area della spira ( $Abucle$ ) mediante il tracciato parallelo intrecciato o la schermatura coassiale blindata risulta indispensabile tanto quanto l'installazione dello scaricatore di sovratensione stesso.

Integrazione Tecnica e Simulazione Elettromagnetica in Vexten Suite

Nell'ambito dell'ingegneria elettromagnetica moderna, il software Vexten Suite agisce come motore centrale di analisi deterministica e stocastica per la modellazione in 3D dei sistemi di captazione contro le scariche atmosferiche. Vexten Suite automatizza le metodologie computazionali specificate dalla norma IEC 62305 attraverso un'integrazione analitica tra la distribuzione spaziale della corrente e la risposta in frequenza della rete dielettrica.

Algoritmo di Ray Tracing e Calcolo Tridimensionale della Sfera Rotolante

Il motore di calcolo geometrico di Vexten Suite proietta numericamente una superficie discreta iper-sferica sulla struttura importata da formati BIM/CAD (STEP, IFC). Il sistema risolve l'interazione tra la maglia della struttura e la sfera rotolante mediante equazioni polinomiali di intersezione sfera-superficie:

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

Dove $\hat{\mathbf{n}}$ è il vettore normale alla superficie di contatto nello spazio tridimensionale. Il programma esegue una scansione angolare infinitesima ( $d\theta, d\phi$ ) registrando tutti i punti di penetrazione in cui la distanza tra la sfera e l'inviluppo della struttura risulta inferiore a zero ( $dreal - R < 0$ ). I punti vulnerabili vengono immediatamente segnalati in modo vettoriale con il rispettivo Livello di Protezione contro i Fulmini (LPL).

Calcolo dell'Impedenza Impulsiva e Integrazione con le Reti di Dispersione a Terra

Vexten Suite integra simultaneamente l'analisi delle correnti di cortocircuito a frequenza industriale (IEC 60909) con l'analisi transitoria impulsiva ad alta frequenza della scarica atmosferica (IEC 62305-3). Il calcolo dell'impedenza transitoria della rete di terra ( $Z(t)$ ) non si limita al valore di resistenza in corrente continua ( $Rdc$ ), ma contempla il comportamento induttivo dei dispersori e dei conduttori interrati attraverso la modellazione a linee di trasmissione distribuite:

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

Dove $l_e$ è la lunghezza effettiva critica dell'elettrodo di terra sotto impulsi ad alta frequenza. Se la lunghezza fisica del conduttore interrato supera $l_e$ , qualsiasi lunghezza aggiuntiva non contribuisce alla riduzione della sovratensione iniziale della scarica. Vexten Suite calcola la lunghezza effettiva critica mediante la relazione empirica:

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

Dove $\rho_{\text{suelo}}$ rappresenta la resistività del terreno ( $\Omega\cdot m$ ) e $T_1$ il tempo di fronte dell'impulso di corrente ( $\mu s$ ).

Flusso Integrado di Progettazione e Validazione in Vexten Suite

Valutazione Automatica del Rischio (IEC 62305-2): Determinazione della frequenza di impatti annuali ( $N_d$ ) e calcolo delle perdite ammissibili (da $R_1$ a $R_4$ ) per fissare computazionalmente il livello di protezione richiesto (da LPL I a IV).
Generazione Geometrica della Zona Protetta: Esecuzione del modulo della Sfera Rotolante in 3D. Il software determina la posizione precisa delle aste captatrici, l'altezza dei pali o il tracciato delle reti a maglia sui tetti per evitare punti diretti di collisione.
Valutazione delle Distanze di Separazione ( $s$ ): Mappatura nodale in 3D della caduta di tensione su ciascuna calata e verifica automatica della regola di separazione isolante rispetto a passacavi, tubazioni e condotti di processo.
Coordinamento Dinamico degli SPD: Generazione delle specifiche degli scaricatori di sovratensione di Tipo 1, Tipo 2 e Tipo 3 richiesti in corrispondenza di ciascun confine di zona LPZ, garantendo il coordinamento dell'energia basato sulla dissipazione del parametro di energia specifica ( $W/R$ ).

Grazie al rigore matematico dell'algoritmo di Vexten Suite, il progettista ottiene un modello elettromagnetico deterministico che invalida il processo per tentativi, garantendo una protezione integrale contro i severi transitori dielettrici, meccanici e termici provocati dalle scariche atmosferiche secondo i più elevati standard internazionali della serie IEC 62305.