Analisi SFRA per la Rilevazione di Deformazioni degli Avvolgimenti

๐——๐—˜๐—™๐—ข๐—ฅ๐— ๐—”๐—–๐—œ๐—ข๐—ก ๐—œ๐—ก๐—ฉ๐—œ๐—ฆ๐—œ๐—•๐—Ÿ๐—˜ ๐——๐—˜ ๐——๐—˜๐—ฉ๐—”๐—ก๐—”๐——๐—ข๐—ฆ: ๐—˜๐—Ÿ ๐—ฅ๐—œ๐—˜๐—ฆ๐—š๐—ข ๐—ฆ๐—œ๐—Ÿ๐—˜๐—ก๐—–๐—œ๐—ข๐—ฆ๐—ข ๐—ง๐—ฅ๐—”๐—ฆ ๐—จ๐—ก ๐—–๐—ข๐—ฅ๐—ง๐—ข๐—–๐—œ๐—ฅ๐—–๐—จ๐—œ๐—ง๐—ข Un transform

Ing. Francisco Ramรญrez

Fondamenti Fisico-Matematici dell'Analisi della Risposta in Frequenza a Scansione (SFRA)

L'Analisi della Risposta in Frequenza a Scansione (SFRA, dall'inglese Sweep Frequency Response Analysis) รจ la metodologia non distruttiva a piรน alta risoluzione per valutare l'integritร  meccanica e idraulico-strutturale del gruppo nucleo-avvolgimenti nei trasformatori di potenza. La premessa fondamentale del metodo risiede nel fatto che un trasformatore di potenza costituisce una rete passiva complessa a parametri distribuiti RLC (resistenza, induttanza, capacitร  e conduttanza), la cui funzione di trasferimento spettrale รจ biiettiva rispetto alla sua geometria fisica e iperstruttura dielettrica.

Ad alte e medie frequenze, la geometria dell'avvolgimento, le distanze dielettriche inter-disco, lo spazio di isolamento rispetto alla cassa e la posizione relativa tra gli avvolgimenti di alta, bassa tensione e regolazione determinano un'impronta digitale spettrale unica. Qualsiasi alterazione dimensionale impercettibile su scala macroscopica โ€” indotta da forze elettrodinamiche di corto circuito, eventi sismici o sollecitudini meccaniche durante il trasporto โ€” altera la distribuzione locale di induttanze e capacitร , provocando uno spostamento inequivocabile nei poli e negli zeri della funzione di trasferimento del sistema.

Modello a Parametri Distribuiti degli Avvolgimenti

Per modellare analiticamente il comportamento spettrale di un avvolgimento, questo viene suddiviso in infinitesimi o elementi discreti di una rete a scala di tipo Ladder Network. Ogni sezione elementare dxdx lungo l'asse assiale dell'avvolgimento possiede induttanze proprie, induttanze mutue rispetto alle sezioni adiacenti e non adiacenti, capacitร  serie tra spire/dischi, capacitร  verso terra (nucleo e cassa) e perdite dielettriche/ohmiche associate.

Le equazioni differenziali parziali che regolano la propagazione dell'onda di tensione v(x,t)v(x,t) e corrente i(x,t)i(x,t) lungo l'avvolgimento corrispondono alla generalizzazione delle Equazioni del Telegrafo nel dominio della frequenza spettrale:

โˆ’โˆ‚V(x,s)โˆ‚x=(R(x,s)+sL(x,s))I(x,s)โˆ’โˆซ0LwsM(x,ฮพ)I(ฮพ,s)โ€‰dฮพ-\frac{\partial V(x, s)}{\partial x} = \left( R(x, s) + s L(x, s) \right) I(x, s) - \int0^{Lw} s M(x, \xi) I(\xi, s) \, d\xi
โˆ’โˆ‚I(x,s)โˆ‚x=(G(x,s)+sCg(x))V(x,s)โˆ’โˆ‚โˆ‚x(sCs(x)โˆ‚V(x,s)โˆ‚x)-\frac{\partial I(x, s)}{\partial x} = \left( G(x, s) + s Cg(x) \right) V(x, s) - \frac{\partial}{\partial x} \left( s Cs(x) \frac{\partial V(x, s)}{\partial x} \right)

Dove:

  • s=ฯƒ+jฯ‰=j2ฯ€fs = \sigma + j\omega = j2\pi f รจ la frequenza complessa di Laplace.
  • R(x,s)R(x, s) รจ la resistenza distribuita dipendente dalla frequenza per effetto pelle (skin effect) ed effetto di prossimitร  nei conduttori (RโˆfR \propto \sqrt{f}).
  • L(x,s)L(x, s) รจ l'induttanza propria distribuita per unitร  di lunghezza.
  • M(x,ฮพ)M(x, \xi) รจ il nucleo dell'operatore integrale di induttanza mutua tra la posizione assiale xx e ฮพ\xi.
  • Cg(x)Cg(x) รจ la capacitร  shunt o distribuita verso terra (tra l'avvolgimento e il nucleo magnetico messo a terra o la parete della cassa).
  • Cs(x)Cs(x) รจ la capacitร  serie distribuita (capacitร  tra spire consecutive e capacitร  tra dischi adiacenti).
  • G(x,s)G(x, s) รจ la conduttanza trasversale dell'isolamento dielettrico cellulosa/olio (G=ฯ‰CgtanโกฮดG = \omega Cg \tan\delta).
  • LwLw rappresenta la lunghezza totale equivalente del conduttore dell'avvolgimento.

Formulazione della Funzione di Trasferimento

La prova SFRA inietta un segnale sinusoidale di ampiezza costante Vin(f)Vin(f) (tipicamente 10 V picco-picco) in un terminale dell'avvolgimento con una scansione logaritmica di frequenza che spazia solitamente da 10Hz10 Hz fino a 2MHz2 MHz, misurando al contempo la risposta in tensione Vout(f)Vout(f) su un altro terminale, terminando entrambi i canali con un'impedenza di riferimento standard di 50ย ฮฉ50\ \Omega.

La risposta viene comunemente espressa in scala decibel di attenuazione o guadagno di tensione:

H(f)dB=20logโก10โˆฃVout(f)Vin(f)โˆฃH(f)_{ dB } = 20 \log10 \left| \frac{Vout(f)}{Vin(f)} \right|

In alternativa, in termini di risposta di fase spettrale ฯ•(f)\phi(f):

ฯ•(f)=argโก(Vout(f)Vin(f))=arctanโก(Im(H(f))Re(H(f)))\phi(f) = \arg\left( \frac{Vout(f)}{Vin(f)} \right) = \arctan \left( \frac{ Im (H(f))}{ Re (H(f))} \right)

Nel dominio operativo di Laplace, l'impedenza di ingresso Zin(s)Zin(s) o la funzione di trasferimento di tensione H(s)H(s) puรฒ essere espressa rigorosamente mediante lo sviluppo in poli (pjp_j) e zeri (ziz_i):

H(s)=Kโ‹…โˆi=1m(sโˆ’zi)โˆj=1n(sโˆ’pj)=Kโ‹…sm+amโˆ’1smโˆ’1+โ‹ฏ+a1s+a0sn+bnโˆ’1snโˆ’1+โ‹ฏ+b1s+b0(nโ‰ฅm)H(s) = K \cdot \frac{\prod_{i=1}^{m} (s - z_i)}{\prod_{j=1}^{n} (s - p_j)} = K \cdot \frac{s^m + a_{m-1}s^{m-1} + \dots + a_1 s + a_0}{s^n + b_{n-1}s^{n-1} + \dots + b_1 s + b_0} \quad (n \ge m)

Le frequenze di risonanza (massimi locali, poli) e antirisonanza (minimi locali, zeri) sono governate dalle radici del polinomio caratteristico. Per una sezione elementare di avvolgimento discreto di ordine NN, le frequenze risonanti primarie ฯ‰k\omega_k rimangono parametrizzate analiticamente da:

ฯ‰kโ‰ˆkฯ€Leq(Cg+k2ฯ€2Cs)\omega_k \approx \frac{k \pi}{\sqrt{Leq \left( Cg + k^2 \pi^2 Cs \right)}}

Questa relazione matematica dimostra in modo indiscutibile l'estrema sensibilitร  del metodo: un'alterazione geometrica millimetrica che modifichi la spaziatura tra i dischi altera la capacitร  serie CsC_s, il che sposta direttamente la frequenza propria ฯ‰k\omega_k nelle regioni di media e alta frequenza dello spettro.

Segmentazione per Bande di Frequenza e Diagnostica Sottostrutturale

L'intero spettro della funzione di trasferimento SFRA (10Hz10 Hz a 2MHz2 MHz) non risponde in modo omogeneo a tutti i componenti fisici del trasformatore. A causa dei diversi comportamenti della reattanza induttiva (XL=2ฯ€fLX_L = 2\pi f L) e della reattanza capacitiva (XC=12ฯ€fCX_C = \frac{1}{2\pi f C}), la risposta in frequenza si suddivide analiticamente in quattro bande funzionali distinte, ciascuna dominata da uno specifico sottinsieme strutturale dell'apparecchiatura.

Banda di Molto Bassa e Bassa Frequenza (< 2 kHz): Dominio Magnetizzante del Nucleo

A frequenze inferiori a 2kHz2 kHz, la reattanza capacitiva degli isolamenti รจ estremamente elevata (XCโ†’โˆžX_C \to \infty), comportandosi come un circuito aperto. Di conseguenza, la corrente iniettata scorre prevalentemente attraverso l'induttanza magnetizzante non lineare del nucleo magnetico (LmL_m) e la resistenza degli avvolgimenti (RdcRdc).

La funzione di trasferimento in questa zona รจ dominata da una prima antirisonanza profonda seguita da una risonanza corrispondente al circuito equivalente costituito dall'induttanza magnetizzante del nucleo in parallelo con la capacitร  equivalente verso terra:

fcore=12ฯ€Lm(Cg+Cinter)f_{ core } = \frac{1}{2\pi \sqrt{L_m \left( C_g + Cinter \right)}}

Qualsiasi fenomeno che alteri la permeabilitร  magnetica del nucleo ฮผr\mu_r, la continuitร  dei lamierini in acciaio al silicio, lo stato della messa a terra del nucleo o la presenza di magnetismo residuo a seguito di un guasto provocherร  un cambiamento massivo nell'ampiezza e nella posizione spettrale di questa banda bassa.

Banda di Media-Bassa Frequenza (2 kHz a 100 kHz): Interazione Inter-Avvolgimento e Strutture Globali

In questo intervallo spettrale, la reattanza dell'induttanza magnetizzante aumenta sostanzialmente e la capacitร  inter-avvolgimento (CinterCinter) inizia ad offrire un percorso a bassa impedenza. L'impronta in questa banda riflette l'interazione mutua tra gli avvolgimenti concentrici di Alta Tensione (AT), Bassa Tensione (BT) e Avvolgimenti di Regolazione (Taps).

Il comportamento รจ regolato dall'induttanza di dispersione a corto circuito (LkL_k) e dalla capacitร  inter-avvolgimento CinterCinter. Le variazioni in questa zona rivelano l'inclinazione globale degli avvolgimenti, spostamenti assiali globali dell'intero gruppo o problemi nelle connessioni interne e nei commutatori di prese sotto carico (OLTC).

Banda di Media-Alta Frequenza (100 kHz a 500 kHz): Struttura Propria dell'Avvolgimento e Geometria Assiale/Radiale

A queste frequenze, l'induttanza di dispersione LkL_k e le capacitร  serie tra dischi o spire (CsC_s) formano una complessa rete multifase di poli e zeri. L'acciaio del nucleo magnetico agisce praticamente come uno schermo al flusso magnetico a causa delle correnti parassite (correnti di Eddy) indotte sulla superficie del laminato, riducendo sostanzialmente la permeabilitร  effettiva del nucleo a ฮผrโ‰ˆ1\mu_r \approx 1.

Pertanto, la risposta dipende quasi esclusivamente dalla geometria interna dell'avvolgimento stesso. Spostamenti locali dei dischi, deformazione per instabilitร  radiale (radial buckling), inclinazione delle spire (tilting) o il cedimento dei tasselli di pressatura assiale provocano variazioni dirette nella capacitร  locale tra dischi Cs(x)C_s(x) e nell'induttanza di dispersione incrementale ฮ”Lk\Delta L_k, manifestandosi come uno spostamento orizzontale (in frequenza) e verticale (in attenuazione) di molteplici picchi risonanti.

Banda di Alta Frequenza (500 kHz a >2 MHz): Cablaggio di Uscita, Connessioni e Isolatori Passanti

All'estremitร  superiore dello spettro (> 500 kHz), l'attenuazione รจ dominata dalle capacitร  parassite di minore entitร , dai conduttori di interconnessione interna (leads), dagli isolatori passanti (bushing) e dall'induttanza propria delle barre flessibili di connessione.

Qualsiasi alterazione osservabile in questa zona che non sia accompagnata da variazioni a frequenze piรน basse รจ solitamente associata a variazioni nella messa a terra dell'apparecchiatura di misura, al degrado delle connessioni dei bushing, a variazioni nella posizione fisica dei conduttori di discesa del trasformatore o a guasti negli schermi di blindaggio elettrostatico.

Analisi Forense dei Modi di Guasto Meccanico e Trasformazioni dello Spettro

Le forze elettromagnetico-dinamiche o forze di Lorentz che agiscono sui conduttori di un trasformatore durante un evento di corto circuito esterno sono proporzionali al quadrato della corrente istantanea (FโˆIsc2F \propto Isc^2). Queste forze possiedono due componenti vettoriali fondamentali: la forza radiale (FrF_r) e la forza assiale (FzF_z). La manifestazione di queste forze imprime deformazioni plastiche irreversibili nella struttura meccanica dell'avvolgimento, le quali si traducono in impronte spettrali SFRA altamente specifiche.

Instabilitร  Radiale (Radial Buckling)

Le forze radiali dirette verso l'interno agiscono comprimendo meccanicamente l'avvolgimento interno (frequentemente l'avvolgimento di bassa tensione che circonda direttamente il nucleo). Se la forza radiale supera il limite critico di deformazione plastica o instabilitร  del conduttore di rame, si genera il fenomeno noto come Buckling.

Esistono due tipologie meccaniche principali:

  • Buckling Libero (Free Buckling): L'avvolgimento si deforma adottando un profilo lobulare privo di punti fissi, formando creste e valli lungo la sua circonferenza.
  • Buckling Forzato (Forced Buckling): L'avvolgimento collassa internamente tra i listelli di supporto o i tasselli assiali della colonna del nucleo, generando deformazioni periodiche schiacciate.

Effetto sullo Spettro SFRA: La deformazione per buckling riduce la distanza fisica radiale tra l'avvolgimento interno e il nucleo, incrementando sostanzialmente la capacitร  verso terra (Cgโ†‘C_g \uparrow). Inoltre, si altera la simmetria dello spazio inter-avvolgimento, modificando la capacitร  mutua (CinterCinter). Sul tracciato SFRA, questo viene diagnosticato mediante uno spostamento sistematico verso frequenze piรน basse dei picchi di risonanza nella banda da 100 kHz a 500 kHz, accompagnato da un'alterazione nell'ampiezza delle valli antirisonanti.

Spostamento Assiale e Disallineamento (Axial Displacement & Tilting)

Le forze assiali si originano dall'asimmetria del flusso magnetico di dispersione alle estremitร  superiore e inferiore dell'avvolgimento. Se gli avvolgimenti di AT e BT non sono perfettamente allineati dal punto di vista magnetico, si genera una componente di forza assiale distruttiva che tende a spingere l'avvolgimento interno verso il basso e quello esterno verso l'alto.

Se la pressione assiale esercitata supera il precarico del sistema di pressatura meccanica (anelli di tenuta e file di tasselli), si verificano due modi di guasto:

  1. Spostamento Assiale Globale: L'intero avvolgimento scivola assialmente lungo la colonna del nucleo.
  2. Inclinazione dei Dischi (Disk Tilting): I dischi dell'avvolgimento adiacenti ruotano tangenzialmente o si inclinano a "zigzag" a causa del cedimento dei blocchi di isolamento inter-disco.

Effetto sullo Spettro SFRA: Lo spostamento assiale modifica drasticamente l'accoppiamento magnetico mutuo e genera una variazione nell'induttanza di dispersione a corto circuito (Lkโ†‘L_k \uparrow). Nell'impronta SFRA, lo spostamento assiale provoca un'alterazione pronunciata nella banda di media-bassa frequenza (2 kHz - 100 kHz), dove si osserva un notevole spostamento della prima valle di risonanza inter-avvolgimento verso sinistra (frequenze piรน basse) e variazioni di ampiezza superiori a 3-6 dB nei picchi intermedi.

Cortocircuito tra le Spire o Cedimento Dielettrico Inter-Disco

Un guasto dielettrico derivante dal degrado della carta isolante Kraft o da una sovratensione transitoria dovuta a manovra o scarica atmosferica puรฒ mettere in corto circuito due o piรน spire consecutive. Un corto circuito tra le spire annulla il flusso magnetico all'interno della spira cortocircuitata a causa della legge di Lenz, generando un'enorme corrente circolante interna che maschera e annulla l'induttanza di quel segmento dell'avvolgimento.

Effetto sullo Spettro SFRA: L'annullamento dell'induttanza equivalente provoca una drammatica perdita della capacitร  di immagazzinamento dell'energia magnetica. Nella risposta SFRA, l'impronta a bassa e media frequenza elimina completamente i picchi risonanti caratteristici dell'avvolgimento interessato, trasformando la curva della funzione di trasferimento in un profilo piatto ad alta attenuazione nelle bande di bassa e media frequenza, simile alla risposta di un circuito fortemente smorzato o in corto circuito permanente.

Algoritmi di Confronto Statistico e Criteri Normativi

L'interpretazione della prova SFRA si basa principalmente sul confronto comparativo delle impronte spettrali mediante tre filosofie metodologiche standard:

  • Confronto Baseline (Fingerprint vs. Attuale): Confronto diretto con l'impronta digitale precedente misurata in fabbrica o prima della messa in servizio. รˆ il metodo che offre la massima sensibilitร  e certezza diagnostica.
  • Confronto Fase a Fase (Phase-to-Phase): Confronto tra le fasi simmetriche della stessa macchina (es. Fase A vs. Fase C) misurate nella medesima configurazione.
  • Confronto con Trasformatore Gemello (Sister Transformer): Confronto con un'unitร  avente un progetto meccanicamente identico e realizzata dallo stesso costruttore.

Per eliminare la soggettivitร  nell'ispezione visiva delle curve di Bode, le norme internazionali come IEEE C57.149, IEC 60076-18 e le raccomandazioni del CIGRE WG A2.26 stabiliscono l'uso di metriche analitiche avanzate, tra cui spiccano il Coefficiente di Correlazione Crociata (RxyRxy o CCCC) e la Differenza Logaritmica dell'Errore Quadratico Medio.

Coefficiente di Correlazione Crociata (RxyRxy)

Dati due vettori di risposta in frequenza espressi in decibel, Xk=Hbase(fk)X_k = Hbase(f_k) e Yk=Hactual(fk)Y_k = Hactual(f_k), valutati discrezionalmente su un numero NN di punti all'interno di sotto-bande di frequenza predefinite:

Rxy=โˆ‘k=1N(Xkโˆ’Xห‰)(Ykโˆ’Yห‰)โˆ‘k=1N(Xkโˆ’Xห‰)2โ‹…โˆ‘k=1N(Ykโˆ’Yห‰)2Rxy = \frac{\sum_{k=1}^{N} \left( X_k - \bar{X} \right) \left( Y_k - \bar{Y} \right)}{\sqrt{\sum_{k=1}^{N} \left( X_k - \bar{X} \right)^2 \cdot \sum_{k=1}^{N} \left( Y_k - \bar{Y} \right)^2}}

Dove Xห‰\bar{X} e Yห‰\bar{Y} corrispondono alle medie aritmetiche delle risposte in ampiezza all'interno della sotto-banda in esame:

Xห‰=1Nโˆ‘k=1NXk,Yห‰=1Nโˆ‘k=1NYk\bar{X} = \frac{1}{N} \sum_{k=1}^{N} X_k, \quad \bar{Y} = \frac{1}{N} \sum_{k=1}^{N} Y_k

Metrica di Distanza Assoluta dello Spettro (DabsDabs)

Per quantificare le deviazioni in ampiezza senza subire distorsioni eccessive dovute a purissimi sfasamenti o scostamenti spettrali di picchi stretti, si utilizza la Distanza Assoluta Spettrale Media:

Dabs(X,Y)=1Nโˆ‘k=1NโˆฃXkโˆ’YkโˆฃDabs(X, Y) = \frac{1}{N} \sum_{k=1}^{N} \left| X_k - Y_k \right|

Criteri di Decisione Diagnostica

In conformitร  con la guida DL/T 911-2016 e gli standard IEEE/IEC, le soglie di severitร  strutturate per il coefficiente di correlazione crociata nella banda critica dell'avvolgimento (100kHzโˆ’500kHz100 kHz - 500 kHz) si classificano secondo la seguente scala rigorosa:

StatoDiagnostico={Normal/Intatto,siRxyโ‰ฅ0.998DeformazioneLeve/Marginale,si0.990โ‰คRxy<0.998DeformacioหŠnModerada,si0.950โ‰คRxy<0.990DeformacioหŠnSevera/FalloMecaหŠnico,siRxy<0.950Stato Diagnostico = \begin{cases} Normal / Intatto , & si Rxy \ge 0.998 \\ Deformazione Leve / Marginale , & si 0.990 \le Rxy < 0.998 \\ Deformaciรณn Moderada , & si 0.950 \le Rxy < 0.990 \\ Deformaciรณn Severa / Fallo Mecรกnico , & si Rxy < 0.950 \end{cases}

Matrice Comparativa di Fenomenologie, Criteri Normativi e Impatto Diagnostico

La seguente tabella consolida in modo esaustivo le impronte spettrali, le variazioni dei parametri elettrici equivalenti, gli standard applicabili, i limiti metrici e le conseguenze operative per i modi di guasto rilevabili tramite SFRA:

Modo di Guasto Meccanico / Strutturale Banda di Frequenza Interessata Variazione Parametrica RLC Equivalente Comportamento dell'Impronta Spettrale SFRA Criterio Metrico Normativo (IEC 60076-18 / IEEE C57.149) Impatto Dielettrico o Operativo Critico
Instabilitร  Radiale (Radial Buckling) Media-Alta Frequenza
(100 kHz โ€“ 500 kHz)
Cgโ†‘C_g \uparrow (10% - 30%)
Cinterโ†“Cinter \downarrow
LkvariabileL_k variabile
Spostamento orizzontale dei picchi risonanti verso frequenze piรน basse. Variazione di ampiezza nei picchi. Rxy<0.95Rxy < 0.95 nella banda media-alta.
Dabs>3.5dBDabs > 3.5 dB.
Riduzione delle distanze di isolamento in olio. Cedimento dielettrico inter-avvolgimento per scarica parziale.
Inclinaciรณn / Spostamento Assiale dei Dischi Media-Bassa a Media-Alta
(20 kHz โ€“ 300 kHz)
Csโ†‘C_s \uparrow o โ†“\downarrow
Lkโ†‘L_k \uparrow (1% - 5%)
MalterataM alterata
Modifica severa nell'ampiezza di picchi/valli. Comparsa di nuove coppie di poli e zeri parassiti. Rxy<0.98Rxy < 0.98 nella banda media.
ฮ”Lk>2%\Delta L_k > 2\% nella prova di reattanza di dispersione.
Cesoia dei tasselli di isolamento chiave. Perdita totale del precarico di serraggio meccanico assiale.
Perdita della Messa a Terra del Nucleo Bassa Frequenza
(10 Hz โ€“ 1 kHz)
CgflotanteC_g flotante
Gโ†’โˆž(parassito)G \to \infty (parassito)
Scomparsa della valle caratteristica a bassa frequenza. Spostamento verticale verso l'alto del tracciato. Deviazione estrema con Rxy<0.90Rxy < 0.90 nella sotto-banda bassa (< 1 kHz). Generazione di scariche distruttive per potenziale flottante. Gassificazione massiva dell'olio dielettrico (DGA).
Magnetismo Residuo nel Nucleo Molto Bassa Frequenza
(10 Hz โ€“ 500 Hz)
ฮผrโ†“(saturazionelocale)\mu_r \downarrow (saturazione locale)
Lmโ†“L_m \downarrow
Spostamento del primo picco di risonanza verso frequenze piรน elevate. L'impronta ritorna allo stato originale dopo la smagnetizzazione. Rxy<0.95Rxy < 0.95 nella banda < 1 kHz, ma normalizzazione assoluta dopo ciclo di smagnetizzazione DC. Correnti di inrush severe durante l'eccitazione. Scatto intempestivo delle protezioni differenziali (87T).
Cortocircuito tra le Spire Bassa, Media e Alta
(10 Hz โ€“ 500 kHz)
Lmโ†’0L_m \to 0
Lkโ†“drammaticoL_k \downarrow drammatico
Rโ†‘(perdite)R \uparrow (perdite)
Crollo totale dell'impronta spettrale. La curva perde i picchi e si trasforma in una risposta ad attenuazione piatta. Rxy<0.80Rxy < 0.80 in piรน bande simultaneamente. Induttanza magnetizzante distrutta. Guasto catastrofico immediato. Scatto della protezione Buchholz (63) e di sovracorrente istantanea (50).
Spostamento delle Connessioni e Inclinazione dei Leads Alta Frequenza
(500 kHz โ€“ 2 MHz)
LleadalterataLlead alterata
CstrayalterataCstray alterata
Spostamento delle creste esclusivamente nella regione superiore a 500 kHz. Stabilitร  assoluta a bassa e media frequenza. Rxy<0.95Rxy < 0.95 esclusivo per f>500kHzf > 500 kHz. Sotto-bande inferiori mantenute con Rxy>0.998Rxy > 0.998. Rischio minore di guasto immediato, salvo possibile riduzione delle distanze dielettriche verso la massa della cassa.

Strategie di Mitigazione Meccanica e Modellazione Quantitativa delle Forze Elettrodinamiche

La prevenzione del danno strutturale rilevato mediante SFRA richiede una modellazione rigorosa delle forze elettrodinamiche che agiscono sulla geometria degli avvolgimenti durante guasti di corto circuito trifase e bifase a terra. Un progetto robusto deve garantire che le sollecitudini dinamiche non superino i limiti elastici del rame o delle leghe di rame drogato con argento (Cu-Ag).

Modellazione delle Forze Radiali di Cortocircuito

La densitร  di forza radiale per unitร  di lunghezza FrF_r su un avvolgimento circolare sottoposto a una corrente massima istantanea di corto circuito asimmetrico ipeakipeak si calcola mediante:

Fr=ฮผ0โ‹…(Nโ‹…ipeak)2โ‹…ฯ€โ‹…Dm2โ‹…heqF_r = \frac{\mu_0 \cdot (N \cdot ipeak)^2 \cdot \pi \cdot Dm}{2 \cdot heq}

Dove:

  • ฮผ0=4ฯ€ร—10โˆ’7H/m\mu_0 = 4\pi \times 10^{-7} H/m รจ la permeabilitร  magnetica del vuoto.
  • NN รจ il numero totale di spire dell'avvolgimento.
  • ipeakipeak รจ la corrente di picco di corto circuito (inclusa la componente continua di asimmetria Kpeakโ‹…2IscKpeak \cdot \sqrt{2} Isc).
  • DmD_m รจ il diametro medio dell'avvolgimento.
  • heqheq รจ l'altezza magnetica equivalente dell'avvolgimento considerando i campi di dispersione di estremitร .

La sollecitudine meccanica di trazione/compressione media di tensione tangenziale (Hoop Stress, ฯƒฮธ\sigma_{\theta}) indotta sulla sezione trasversale del conduttore AcA_c รจ data analiticamente da:

ฯƒฮธ=Fr2ฯ€โ‹…Nโ‹…Ac=ฮผ0โ‹…Nโ‹…ipeak2โ‹…Dm4โ‹…heqโ‹…Ac\sigma_{\theta} = \frac{F_r}{2\pi \cdot N \cdot A_c} = \frac{\mu_0 \cdot N \cdot ipeak^2 \cdot D_m}{4 \cdot heq \cdot A_c}

Per evitare l'instabilitร  permanente (buckling), la condizione di progetto rigorosa impone che:

ฯƒฮธโ‰คฯƒammissibile=ฯƒ0.2FS\sigma_{\theta} \le \sigma_{ ammissibile } = \frac{\sigma_{0.2}}{ FS }

Dove ฯƒ0.2\sigma_{0.2} รจ il limite elastico convenzionale allo 0.2% di deformazione del materiale conduttore (tipicamente 120โˆ’220MPa120 - 220 MPa per rame incrudito) e FSFS รจ il fattore di sicurezza statico/dinamico (FSโ‰ฅ1.5FS \ge 1.5).

Calcolo delle Forze Assiali di Compressione e Pressatura

La forza assiale totale di compressione FzF_z che agisce tendendo a far collassare l'avvolgimento verso il suo centro magnetico medio รจ descritta dall'integrazione del prodotto del flusso radiale di dispersione Br(z)B_r(z) per la corrente totale:

Fz(z)=2ฯ€โˆซ0zDmโ‹…Ndirโ‹…ipeakโ‹…Br(ฮพ)โ€‰dฮพF_z(z) = 2\pi \int0^{z} D_m \cdot Ndir \cdot ipeak \cdot B_r(\xi) \, d\xi

La sollecitudine di schiacciamento assiale sui tasselli di isolamento cellulosico di supporto non deve superare la resistenza alla compressione dinamica a caldo del presspan (Pressboard, tipicamente ฯƒcompโ‰ค40โˆ’70MPa\sigma_{\text{comp}} \le 40 - 70 MPa).

Tecniche Costruttive Moderne di Mitigaciรณn

  • Conduttori Trasposti Continui con Resina Epossidica (Bonded CTC): Utilizzo di cavi CTC essiccati al forno con vernice termoadesiva B-stage. Questo unisce meccanicamente tutte le singole piattine di rame all'interno del conduttore in un'unica trave rigida monolitica, incrementando il modulo di resistenza della sezione e aumentando la resistenza all'instabilitร  radiale di un fattore da 3 a 5 rispetto al conduttore singolo.
  • Pressatura Assiale Dinamica Serrata da Molle o Cilindri Idraulici: Installazione di dispositivi di applicazione di un precarico costante per compensare automaticamente lo scorrimento viscoso (creep) della carta e dell'isolamento in legno lamellare lungo la vita utile dell'apparecchiatura, garantendo che l'avvolgimento non perda mai la forza di serraggio assiale minima raccomandata (Fprecarga>Fz,dinamicoFprecarga > F_{z, dinamico}).

Integrazione Avanzata con Vexten Suite per l'Analisi SFRA e l'Integritร  Strutturale

Nell'ecosistema analitico di Vexten Suite, il modulo specializzato di diagnostica SFRA opera come motore centrale per l'analisi forense e la previsione dello stato di salute strutturale degli asset di trasformazione di potenza. Il software integra l'elaborazione analitica dei dati spettrali con simulazioni deterministiche del sistema elettrico.

Correlazione dei Guasti di Cortocircuito tramite IEC 60909 e IEEE 141

Il motore di calcolo dei cortocircuiti di Vexten Suite (modulo Vexten Fault Engine) simula i transitori elettromagnetici di guasto nelle reti di trasmissione e distribuzione secondo le metodologie IEC 60909 e IEEE 141. Quando si verifica un guasto in prossimitร  della sottostazione, Vexten Suite calcola con precisione la corrente di picco asimmetrica istantanea (ipi_p) che ha attraversato il trasformatore:

ip=ฮบโ‹…2โ‹…Ikโ€ฒโ€ฒi_p = \kappa \cdot \sqrt{2} \cdot I_k''

Dove ฮบ\kappa รจ il fattore di picco in funzione del rapporto R/XR/X del sistema nel punto di guasto e Ikโ€ฒโ€ฒI_k'' รจ la corrente iniziale di corto circuito simmetrica. A partire da questa corrente calcolata, Vexten Suite determina la mappa tridimensionale delle sollecitudini meccaniche dinamiche (Fr,FzF_r, F_z) subite dall'avvolgimento. Se i limiti di fatica accumulata del conduttore vengono superati, il sistema avvisa automaticamente l'ingegnere di manutenzione, emettendo un ordine di lavoro prioritario per l'esecuzione di una prova SFRA in loco.

Algoritmo Vexten di Impronta Spettrale e Decomposizione Sotto-Banda Automatica

Il pacchetto software Vexten SFRA Analyzer automatizza l'importazione di file di tracciato standardizzati XML/ASCII (formato CIGRE / IEEE C57.149). Esegue nativamente i seguenti processi matematici avanzati:

  1. Allineamento Spettrale mediante Spline Cubica: Interpolazione logaritmica della risposta per armonizzare punti di frequenza discontinui acquisiti con diverse apparecchiature di prova (es. Omicron, Doble, Megger).
  2. Decomposizione Adattativa in Quattro Sotto-Bande Standard: Segmentazione automatica del vettore spettrale nelle regioni fisiche del nucleo, interazione inter-avvolgimento, struttura dell'avvolgimento e connessioni.
  3. Matrice del Coefficiente di Correlazione Crociata Vettoriale e Metrica SDM (Spectrum Distance Metrics): Esecuzione parallela dell'analisi statistica multivariabile tra l'impronta originale e le misurazioni di campo recenti.
CCVexten=[RxyBajaRxyMedโˆ’BajaRxyMedโˆ’AltaRxyAlta]โ€…โ€ŠโŸนโ€…โ€ŠMotorediInferenzaDetezionediBuckling/Tilting\mathbf{CC}_{ Vexten } = \begin{bmatrix} Rxy^{ Baja } \\ Rxy^{ Med-Baja } \\ Rxy^{ Med-Alta } \\ Rxy^{ Alta } \end{bmatrix} \implies Motore di Inferenza Detezione di Buckling / Tilting

Se la metrica calcolata da Vexten Suite nella banda media-alta scende al di sotto di 0.980.98, il modulo attiva allarmi precoci e integra i risultati con le prove di capacitร  e fattore di potenza dell'isolamento (C1/C2) e con la prova di induttanza di dispersione a corto circuito (Short-Circuit Impedance test secondo IEC 60076-5), isolando l'esatto modo di guasto e stimando la percentuale di vita utile meccanica residua del trasformatore.