Ferrorisonanza nei Trasformatori e Coordinamento DC gPV negli Impianti Automotive

¿Sobretensiones de 3.2 pu y disparos intempestivos al maniobrar celdas de soldadura robotizada? Desliza este dossier técnico para mitigar ferrorresonancia y coo

Ing. Francisco Ramírez

Fenomenologia Non Lineare della Ferrorisonanza nei Trasformatori degli Impianti Automobilistici

La ferrorisonanza è un fenomeno di risonanza non lineare che differisce fondamentalmente dalla risonanza lineare classica. In un circuito lineare RLC, la risonanza si verifica a una singola frequenza angolare determinata dall'interazione di elementi reattivi costanti:

ω0=1LC\omega_0 = \frac{1}{\sqrt{LC}}

Nei sistemi di distribuzione degli impianti automobilistici, caratterizzati da fitte reti di cavi sotterranei di media tensione con elevata capacità shunt (C0C_0) e trasformatori di distribuzione con nuclei ferromagnetici saturabili, l'induttanza magnetizzante (LmL_m) non è costante. Questa induttanza è una funzione non lineare altamente sensibile al flusso magnetico concatenato (ϕ\phi), descritta dalla relazione tra la corrente di magnetizzazione (imi_m) e il flusso:

im(ϕ)=aϕ+bϕni_m(\phi) = a\cdot\phi + b\cdot\phi^n

Dove aa e bb sono coefficienti di accoppiamento magnetico dipendenti dalla geometria e dal materiale del nucleo, e nn è un esponente dispari (tipicamente n7n \ge 7 per i moderni lamierini magnetici al silicio a grani orientati). Questa non linearità introduce molteplici punti di funzionamento stabili e instabili per la stessa frequenza di eccitazione industriale (50 Hz o 60 Hz). Il comportamento dinamico del sistema viene modellato mediante l'equazione differenziale non lineare di Duffing:

d2ϕdt2+1RpC0dϕdt+1C0im(ϕ)=ωEmcos(ωt+θ)\frac{d^2\phi}{dt^2} + \frac{1}{R_p C_0}\frac{d\phi}{dt} + \frac{1}{C_0}i_m(\phi) = \omega E_m \cos(\omega t + \theta)

Dove RpR_p rappresenta le perdite nel rame e le perdite nel nucleo (parassite e per isteresi), ed EmE_m è l'ampiezza della tensione di eccitazione. La transizione da uno stato operativo normale a bassa induzione a uno stato ferrorisonante ad alta induzione avviene attraverso biforcazioni nello spazio delle fasi.

Meccanismi di Biforcazione e Modi di Ferrorisonanza

La ferrorisonanza viene classificata in base alla topologia della sua orbita nello spazio delle fasi e al contenuto spettrale delle variabili di stato:

  • Modo Fondamentale (Periodo 1): Le variabili di tensione e flusso sono periodiche con lo stesso periodo T=1/fT = 1/f della sorgente di eccitazione. È caratterizzato da gravi sovratensioni e distorsioni simmetriche della forma d'onda della corrente di magnetizzazione.
  • Modo Subarmonico: Le orbite nello spazio delle fasi si chiudono dopo un multiplo intero del periodo della sorgente (nTnT, tipicamente 3T3T o 1/3 della frequenza fondamentale). Questo modo è comune nei sistemi con basso smorzamento e può causare gravi vibrazioni meccaniche dovute alle forze elettrodinamiche negli avvolgimenti.
  • Modo Quasi-Periodico: Si presenta quando la traiettoria nello spazio delle fasi non è chiusa, ma avvolge un toroide. Lo spettro di frequenza mostra frequenze non armoniche sovrapposte alla fondamentale.
  • Modo Caotico: Caratterizzato da un'estrema sensibilità alle condizioni iniziali. Le traiettorie nello spazio delle fasi si comportano come attrattori strani, con uno spettro di frequenza continuo di rumore di fondo e picchi disordinati. Questo modo sottopone l'isolamento del trasformatore a uno stress termico e dielettrico imprevedibile.
λLyapunov>0\lambda_{\text{Lyapunov}} > 0

L'esistenza di un esponente di Lyapunov positivo (λLyapunov>0\lambda_{\text{Lyapunov}} > 0) conferma la presenza del modo caotico, in cui il sistema perde la prevedibilità deterministica a lungo termine, accelerando la degradazione dell'isolamento solido (carta Kraft) e liquido (olio minerale o estere sintetico).

---

Analisi Forense dei Guasti da Ferrorisonanza e Transitori di Commutazione

Nel contesto della produzione robotizzata dell'industria automobilistica, la continuità di servizio e la qualità dell'energia sono critiche. I guasti catastrofici nelle cabine di trasformazione secondarie di media tensione (MT) sono spesso causati da transitori debolmente smorzati che innescano la ferrorisonanza. Gli scenari tipici di innesco forense includono:

  • Funzionamento Monofase delle Apparecchiature di Manovra (Fusibili o Sezionatori): L'apertura o la chiusura non simultanea dei poli di un interruttore o la fusione di un singolo fusibile MT a causa di un guasto a valle lascia il trasformatore alimentato da una o due sole fasi. Le fasi aperte rimangono accoppiate capacitivamente attraverso le capacità mutue dei cavi di alimentazione (CgC_g e CmC_m).
  • Rottura dei Conduttori di Alimentación: Un cavo sotterraneo con schermo metallico a terra che subisce un guasto all'isolamento e la successiva interruzione fisica di un conduttore crea una topologia serie-parallelo ideale per la risonanza tra la capacità dello schermo del tratto scollegato e la reattanza magnetizzante non lineare del trasformatore.
  • Sistemi a Neutro Isolato o con Messa a Terra Inefficace: I trasformatori con collegamento stella-triangolo o triangolo-stella alimentati da sistemi con neutro flottante sono soggetti a gravi spostamenti del punto di neutro durante i guasti a terra monofase, saturando i nuclei magnetici delle fasi sane.

Degradazione Termica e Dielettrica dell'Isolamento

Durante la ferrorisonanza, la corrente de magnetizzazione che fluisce attraverso gli avvolgimenti del trasformatore può aumentare fino a 100 veces il suo valore nominale a vuoto, trasformandosi in correnti di inrush continue e ripetitive. Ciò genera perdite per effetto Joule localizzate negli avvolgimenti MT:

PJoule(t)=Rac(f)Irms2(t)PJoule(t) = Rac(f) \cdot Irms^2(t)

Dove Rac(f)Rac(f) è la resistenza in corrente alternata corretta per l'effetto pelle (skin effect) e di prossimità alle alte frequenze armoniche. Simultaneamente, il flusso di dispersione magnetica in condizioni di profonda saturazione abbandona il nucleo in lamierino al silicio e penetra nelle strutture metalliche di supporto, nelle pareti della cassa e nei bulloni di serraggio, inducendo massicce correnti parassite (correnti di Foucault) che causano punti caldi (hotspot) locali che superano rapidamente i 200 °C.

L'olio dielettrico subisce una pirolisi termica accelerata, generando gas infiammabili come idrogeno (H2H_2), etilene (C2H4C_2H_4) e acetilene (C2H2C_2H_2), rilevabili tramite l'Analisi dei Gas Disciolti (DGA) secondo la norma IEC 60599. Lo stress dielettrico è aggravato da sovratensioni transitorie con elevati gradienti di salita della tensione (dv/dtdv/dt), che compromettono la rigidità dielettrica dell'isolamento delle spire d'estremità dell'avvolgimento a causa della distribuzione non uniforme del potenziale capacitivo lungo la bobina.

Degradazione degli Scaricatori di Sovratensione (MOV)

Gli scaricatori di sovratensione a ossido metallico (MOV) installati per proteggere il trasformatore da sovratensioni atmosferiche e di manovra diventano vittime collaterali durante la ferrorisonanza. Un MOV è progettato per assorbire energia transitoria di breve durata (da microsecondi a millisecondi). Tuttavia, la ferrorisonanza sostiene sovratensioni temporanee (TOV) per secondi o minuti. L'energia dissipata dal varistore supera la sua capacità di assorbimento termico ammissibile:

EMOV=0tfallvMOV(t)iMOV(t)dt>Emax_ammissibileEMOV = \int_0^{tfall} vMOV(t) \cdot iMOV(t) \, dt > E_{max\_ammissibile}

Questo provoca una deriva termica (thermal runaway). All'aumentare della temperatura del varistore di ZnO, la sua corrente di fuga in stato di normale conduzione aumenta in modo esponenziale a causa del suo coefficiente di temperatura negativo, riducendo la sua resistenza ohmica interna fino alla sua distruzione fisica, esplosione e cortocircuito franco sulle sbarre MT.

---

Coordinamento delle Protezioni in Corriente Continua: Fusibili gPV nei Sistemi di Accumulo e Generazione Fotovoltaica Integrata

I moderni impianti automobilistici integrano risorse energetiche distribuite (DER), come impianti fotovoltaici su tetto e sistemi di accumulo dell'energia a batteria (BESS), per alimentare carichi critici e livellare i picchi di domanda delle linee di saldatura robotizzate. La protezione di questi sistemi in corrente continua (DC) richiede un coordinamento preciso mediante fusibili con categoria di utilizzo gPV (progettati specificamente per la protezione di sistemi fotovoltaicos secondo la norma IEC 60269-6).

Fisica dell'Interruzione dell'Arco Elettrico in Corrente Continua

A differenza dei sistemi in corrente alternata (AC), dove la corrente attraversa lo zero in modo naturale due volte per ciclo, facilitando l'estinzione dell'arco elettrico, nei sistemi DC la corrente non presenta passaggi per lo zero. L'estinzione dell'arco elettrico richiede che il dispositivo di protezione generi attivamente una tensione d'arco (VarcoVarco) sostanzialmente superiore alla tensione nominale del sistema in corrente continua (VsysVsys):

Varco(t)=Vsys+Ldidt+Ri(t)Varco(t) = Vsys + L\frac{di}{dt} + R\cdot i(t)

Affinché la corrente decresca (di/dt<0di/dt < 0), è imperativo che:

Varco(t)>VsysRi(t)Varco(t) > Vsys - R\cdot i(t)

La sabbia di quarzo ad alta purezza che riempie il corpo del fusibile gPV funge da mezzo di raffreddamento e deionizzazione del canale di plasma dell'arco. Quando l'elemento conduttore in argento puro (Ag) si fonde in condizioni di sovracorrente, l'arco elettrico si allunga e si segmenta nei restringimenti (colli) del fusibile. L'energia dell'arco fonde la silice circostante, formando una struttura vetrosa non conduttiva chiamata fulgurite, che assorbe l'energia termica e interrompe il flusso di corrente.

Coordinamento dei Fusibili gPV in Condizioni de Cortocircuito Estreme

Nella progettazione della protezione DC per impianti automobilistici con BESS e generazione fotovoltaica integrata, il coordinamento si complica a causa dei due profili di corrente di guasto estremi del sistema:

  • Corrienti di Guasto di Bassa Entità (Tipiche delle Stringhe Fotovoltaiche): I moduli fotovoltaici sono sorgenti di corrente limitata dall'irraggiamento solare. Le correnti di cortocircuito tipiche (IscIsc) sono solo dal 10% al 20% superiori alla corrente di massima potenza (ImpImp). I fusibili gPV devono essere in grado di intervenire in modo sicuro a fronte di queste correnti di guasto minime (Imin1.35InImin \approx 1.35 \cdot I_n fino a 2In2 \cdot I_n) in un tempo inferiore a quello che comprometterebbe l'integrità dei conduttori DC e dei moduli.
  • Corrienti de Guasto di Elevata Entità e Alta Costante di Tempo (Tipiche dei BESS): I banchi di batterie agli ioni di litio presentano una resistenza interna estremamente bassa, il che dà origine a correnti di cortocircuito presunte molto elevate (decine di kiloampere) con costanti di tempo di rete molto rapide (τ=L/R\tau = L/R) o lente a seconda della lunghezza delle sbarre. Il fusibile gPV deve possedere un potere di interruzione estremo (tipicamente fino a 100 kA DC) e limitare l'energia passante (integrale di Joule, I2tI^2t).

L'integrale di Joule di pre-arco (I2tpreI^2t_{pre}) e l'integrale di Joule di funzionamento totale (I2topI^2t_{op}) definiscono i limiti di selettività tra fusibili installati in serie (ad esempio, fusibile di stringa vs. fusibile di sotto-gruppo vs. fusibile principale dell'inverter):

I2tpre=0tprei2(t)dt=CostanteperundatomaterialeegeometriadelfusibileI^2t_{pre} = \int_0^{tpre} i^2(t) \, dt = Costante per un dato materiale e geometria del fusibile

Per garantire un coordinamento selettivo totale tra due fusibili in serie (a monte e a valle), deve essere rigorosamente rispettata la seguente relazione di energia su tutto lo spettro delle correnti di cortocircuito:

I2top_downstream<I2tpre_upstreamI^2t_{op\_downstream} < I^2t_{pre\_upstream}

Se questa condizione non viene soddisfatta, entrambi i fusibili si fonderanno simultaneamente (perdita di selettività), causando un disservizio non necessario in intere sezioni della microrete dell'impianto automobilistico.

---

Tabella Comparativa di Standard, Parametri Elettrici e Limiti Critici

La seguente tabella detalla i parametri di progettazione e i limiti operativi normalizzati secondo gli standard internazionali per la prevenzione della ferrorisonanza e il corretto coordinamento delle protezioni DC gPV.

Parametro / Fenomeno Standard di Riferimento Limite Critico Operativo Conseguenza Dielettrica / Termica Strategia di Mitigazione / Progettazione
Sovratensione Temporanea (TOV) da Ferrorisonanza IEEE C57.105 / IEC 60076-1 V^1.5p.u.\hat{V} \ge 1.5 p.u. per più di 2 secondi. Cedimento dell'isolamento tra le spire, degradazione termica della carta Kraft, arco interno. Installazione di resistenze di smorzamento sull'avvolgimento a triangolo aperto.
Distorsione Armonica Totale di Tensione (THDvTHD_v) IEEE 519 / IEC 61000-4-30 THDv>8.0%THD_v > 8.0\% sulle sbarre MT (da 1 kV a 69 kV). Aumento delle perdite per isteresi e delle correnti parassite nel nucleo del trasformatore. Filtri armonici sintonizzati e reattori di disaccordo nei banchi di condensatori.
Potere di Interruzione dei Fusibili gPV IEC 60269-6 / UL 248-19 Isc_presunta>Ibc_fusibileI_{sc\_presunta} > I_{bc\_fusibile} (Tipicamente < 10 kA o 50 kA). Rottura fisica violenta del fusibile, arco persistente, rischio di incendio nella cabina di trasformazione. Selezione di fusibili con elevato potere di interruzione certificato (100kADC\ge 100 kA DC).
Costante di Tempo di Guasto DC (τ\tau) IEC 60269-1 / IEC 60949 τ=L/R>15ms\tau = L/R > 15 ms in cortocircuito DC. Estrema difficoltà nell'estinguere l'arco elettrico; prolungamento del tempo di eliminazione del guasto. Utilizzo di fusibili gPV testati con costanti di tempo severe (τ1020ms\tau \approx 10-20 ms).
Fattore di Declassamento Termico dei Cavi IEC 60287 / NEC Art. 310 Temperatura di funzionamento sul conduttore XLPE 90°C\ge 90 °C. Invecchiamento termico accelerato dell'isolamento in polietilene reticolato, perforazione. Calcolo della portata considerando le armoniche mediante il fattore di correzione delle perdite di Vexten.
Rigidità Dielettrica dell'Olio Isolante IEC 60156 / ASTM D1816 Tensione di scarica dielettrica <30kV< 30 kV. Formazione di archi interni distruttivi sotto sovratensioni transitorie ordinarie. Tratamiento di termovuoto, degasaggio e rigenerazione attiva dell'olio minerale.
---

Strategie di Mitigazione e Progettazione di Ingegneria di Dettaglio

Per evitare in modo proattivo l'insorgenza della ferrorisonanza nelle cabine dell'impianto automobilistico, devono essere implementate soluzioni di ingegneria di dettaglio basate sull'analisi dei parametri fisici.

Progettazione di Circuiti di Smorzamento (Damping Loops)

Una delle soluzioni più efficaci per mitigare la ferrorisonanza nei trasformatori di misura di tensione (TV) e nei trasformatori di distribuzione collegati a stella con neutro a terra consiste nell'installazione di una resistenza di smorzamento (RdR_d) collegata ai terminali di un avvolgimento secondario ausiliario collegato a triangolo aperto (schema a triangolo aperto o neutro rotto):

Rd13ωC0R_d \le \frac{1}{3 \cdot \omega \cdot C_0}

Dove C0C_0 rappresenta la capacità omopolare per fase della rete di cavi di media tensione collegata al trasformatore. Questa resistenza agisce dissipando l'energia delle oscillazioni ferrorisonanti quando si presenta una tensione di sequenza zero dovuta a un'asimmetria o a un transitorio di commutazione. La potenza termica nominale (PRdPRd) della resistenza di smorzamento deve essere calcolata per sopportare la tensione di cortocircuito omopolare temporanea senza bruciarsi:

PRd=(3Vsec_nominal)2RdtsupportotcicloPRd = \frac{(3 \cdot V_{sec\_nominal})^2}{R_d} \cdot \frac{tsupporto}{tciclo}

Criteri di Commutazione Tripolare Simultanea

Deve essere rigorosamente vietato l'uso di fusibili di media tensione o sezionatori unipolari a monte di trasformatori di distribuzione soggetti a ferrorisonanza. Al loro posto, devono essere specificati interruttori automatici a tripla interruzione simultanea (in SF6 o sottovuoto) controllati da relè di protezione numerici ad alta velocità. Ciò garantisce che le tre fasi vengano energizzate o disenergizzate con un disallineamento temporale massimo tra i poli inferiore a 2 millisecondi:

Δtpoli2ms\Delta tpoli \le 2 ms

Questa simultaneità impedisce la formazione di circuiti risonanti serie-parallelo asimmetrici che caratterizzano le fasi aperte.

Dimensionamento e Coordinamento dei Fusibili gPV sul Lato DC

Per il dimensionamiento rigoroso dei fusibili gPV negli inverter e nelle stringhe fotovoltaiche dell'impianto automobilistico, si devono seguire le seguenti equazioni di coordinamento:

La tensione nominale del fusibile (Vn_fuseV_{n\_fuse}) deve essere corretta in base alla temperatura ambiente della cabina e alla tensione massima a vuoto della stringa fotovoltaica (VocVoc) alla temperatura minima di progetto del sito (TminTmin):

Vn_fuseVoc_string[1+βvoc(Tmin25)]V_{n\_fuse} \ge V_{oc\_string} \cdot \left[1 + \beta_{voc} \cdot (Tmin - 25)\right]

Dove βvoc\beta_{voc} è il coefficiente di temperatura della tensione a vuoto del pannello fotovoltaico (%/°C\%/ °C).

La corrente nominale del fusibile (In_fuseI_{n\_fuse}) deve considerare il declassamento per temperatura ambiente (KtempKtemp) e per altitudine (KaltKalt) se l'impianto si trova a più di 2000 metri sopra il livello del mare:

In_fuse1.4Isc_stringKtempKaltI_{n\_fuse} \ge \frac{1.4 \cdot I_{sc\_string}}{Ktemp \cdot Kalt}

Il fattore di sicurezza di 1,4 garantisce che il fusibile non subisca un invecchiamento precoce dovuto a correnti cicliche causate dalle fluttuazioni dell'irraggiamento solare e dalle richieste pulsanti di energia del sistema di saldatura robotizzata dell'impianto.

---

Modellazione e Simulazione Avanzata tramite Vexten Suite

La suite di ingegneria elettrica Vexten Suite è fondamentale per prevedere, analizzare e risolvere i problemi di ferrorisonanza e coordinamento delle protezioni negli impianti industriali complessi. Di seguito viene presentata la metodologia di modellazione e analisi utilizzando i moduli avanzati di questo strumento.

Modulo di Cortocircuito (IEC 60909 / IEEE 141)

Per determinare le correnti di guasto massime e minime nei punti di accoppiamento comune (PCC) e sulle sbarre di distribuzione in corrente continua della microrete dell'impianto automobilistico, si utilizza il motore di calcolo del cortocircuito di Vexten Suite. Il software calcola simultaneamente i parametri di cortocircuito secondo gli standard internazionali:

  • Corrente di cortocircuito di picco (ipi_p): Cruciale per determinare gli sforzi meccanici sulle sbarre e il potere di chiusura delle apparecchiature di manovra.
  • Corrente di cortocircuito termica equivalente (IthIth): Definisce lo stress termico sopportato dai conduttori e dai fusibili durante il tempo di eliminazione del guasto.

Vexten Suite modella in modo preciso l'impedenza subtransitoria dei generatori e il contributo dei motori asincroni e sincroni della linea di assemblaggio automobilistica, applicando la matrice delle impedenze di nodo [Zbus][Zbus]:

Ifcc=[Zbus]1Vpre_guastoIfcc = [Zbus]^{-1} \cdot V_{pre\_guasto}

Nel dominio della corrente continua, il software risolve le equazioni differenziali di scarica dei banchi di batterie (BESS) e dei condensatori del circuito intermedio DC (DC-link) dei variatori di frequenza (VFD), consentendo di parametrizzare la resistenza interna non lineare delle celle al litio e l'induttanza della rete di sbarre.

Modulo di Dimensionamento dei Cavi e Declassamento Armonico (IEC 60287 / NEC 310)

Le linee di produzione robotizzate contengono centinaia di servomotori e sistemi di saldatura a resistenza che iniettano massicce correnti armoniche nella rete di distribuzione di bassa e media tensione. Il modulo di dimensionamento dei cavi di Vexten Suite implementa gli algoritmi iterativi della norma IEC 60287 per il calcolo della temperatura a regime del conduttore.

Il software esegue un declassamento armonico avanzato considerando l'incremento della resistenza in corrente alternata (RacRac) dovuto all'effetto pelle (skin effect) e all'effetto di prossimità in presenza di armoniche ad alta frequenza:

Rac(h)=Rdc[1+ys(h)+yp(h)]Rac(h) = Rdc \cdot \left[1 + y_s(h) + y_p(h)\right]

Dove ys(h)y_s(h) e yp(h)y_p(h) sono rispettivamente i fattori di effetto pelle e di prossimità per l'armonica di ordine hh. Vexten Suite calcola automaticamente il fattore di declassamento armonico (KharmonicKharmonic):

Kharmonic=1+h=1H(IhI1)21+h=1H[(IhI1)2Rac(h)Rac(1)]Kharmonic = \sqrt{\frac{1 + \sum_{h=1}^{H} \left(\frac{I_h}{I_1}\right)^2}{1 + \sum_{h=1}^{H} \left[\left(\frac{I_h}{I_1}\right)^2 \cdot \frac{Rac(h)}{Rac(1)}\right]}}

Questo fattore corregge la portata (ampacity) ammissibile del conduttore in rame con isolamento in XLPE dei circuiti di alimentazione dei trasformatori, prevenendo surriscaldamenti che potrebbero degradare l'isolamento e fungere da innesco per guasti a terra in grado di scatenare la ferrorisonanza.

Modulo di Rifasamento e Mitigazione della Risonanza

Per evitare la risonanza armonica parallela tra i banchi di condensatori di rifasamento e l'induttanza di cortocircuito del trasformatore dell'impianto automobilistico, Vexten Suite consente di progettare filtri armonici disaccordati (detuned filters). Il software calcola la frequenza di accordo del reattore di blocco (frf_r) richiesta per spostare il punto di risonanza parallela al di fuori delle frequenze delle armoniche dominanti inyectate dai raddrizzatori dei robot (tipicamente la 5ª, 7ª, 11ª e 13ª armonica):

fr=f11pf_r = f_1 \cdot \sqrt{\frac{1}{p}}

Dove pp è il fattore di disaccordo o sintonizzazione (tipicamente p=7%p = 7\% per una frequenza di accordo di 189 Hz in reti a 50 Hz, o p=14%p = 14\% per evitare l'amplificazione della 3ª armonica). L'ambiente grafico di simulazione di Vexten Suite fornisce le curve di risposta in frequenza dell'impedenza del sistema vista da qualsiasi nodo dell'impianto:

Znode(f)=Vnode(f)Iinject(f)Znode(f) = \frac{Vnode(f)}{Iinject(f)}

Questo consente di convalidare visivamente e analiticamente l'assenza di picchi di alta impedenza alle frequenze critiche di commutazione, eliminando la possibilità di sovratensioni ferrorisonanti transitorie causate dalla manovra dei banchi di condensatori.

Attraverso il suo motore di simulazione dei transitori elettromagnetici (EMT) integrato, Vexten Suite consente di modellare l'isteresi magnetica del nucleo del trasformatore mediante il modello di Jiles-Atherton. Questo fornisce agli ingegneri di progettazione automobilistica la capacità di prevedere l'insorgenza della ferrorisonanza sotto diversi scenari di manovra degli interruttori, garantendo la robustezza operativa dell'infrastruttura elettrica dell'impianto prima della sua costruzione o modifica sul campo.