Perdite per Isteresi e Correnti Parassite nell'Acciaio CRGO

¿Por qué la degradación del aislamiento interlaminar en acero CRGO dispara las pérdidas en vacío P0 hasta un 35%? Desliza este dossier técnico para dominar la d

Ing. Francisco Ramírez

Termodinamica Elettromagnetica e Microstruttura dell'Acciaio al Silicio CRGO

Il nucleo magnetico di un trasformatore di potenza opera come il circuito di accoppiamento di flusso fondamentale tra gli avvolgimenti, lavorando sotto regimi ciclici non lineari di eccitazione elettromagnetica. L'efficienza globale e la densità di potenza del trasformatore sono condizionate in modo determinante dal comportamento microstrutturale e cristallografico del materiale ferromagnetico impiegato. Lo standard industriale predominante per i trasformatori di media, alta e altissima tensione (EHV/UHV) è l'acciaio elettrico a grani orientati laminato a freddo (Cold-Rolled Grain-Oriented Silicon Steel, CRGO), ottimizzato mediante trattamenti termomeccanici per allineare i suoi assi cristallografici di facile magnetizzazione lungo la direzione di laminazione.

L'acciaio CRGO presenta una struttura cristallina cubica a corpo centrato (BCC, reticolo α\alpha-Fe), in cui l'asse cristallografico 100\langle 100 \rangle costituisce la direzione di minima energia di anisotropia magnetocristallina, mentre le direzioni 110\langle 110 \rangle e 111\langle 111 \rangle rappresentano rispettivamente assi di anisotropia intermedia e dura. Mediante una rigorosa sequenza di laminazione a freddo e ricotture di ricristallizzazione secondaria (inibita tramite precipitati di solfuro di manganese, MnSMnS, o nitruro di alluminio, AlNAlN), si sviluppa la cosiddetta tessitura di Goss, caratterizzata dall'orientamento cristallografico {110}001\{110\}\langle 001 \rangle. In questa configurazione, i piani {110}\{110\} risultano paralleli alla superficie della lamiera e le direzioni 001\langle 001 \rangle si orientano parallelamente alla direzione longitudinale di laminazione (RD, Rolling Direction).

Ea=K1(α12α22+α22α32+α32α12)+K2(α12α22α32)E_a = K_1 \left( \alpha_1^2 \alpha_2^2 + \alpha_2^2 \alpha_3^2 + \alpha_3^2 \alpha_1^2 \right) + K_2 \left( \alpha_1^2 \alpha_2^2 \alpha_3^2 \right)

Dove EaE_a è la densità di energia di anisotropia magnetocristallina (J/m3J/m^3), K14.8×104J/m3K_1 \approx 4.8 \times 10^4 J/m^3 e K21.5×104J/m3K_2 \approx 1.5 \times 10^4 J/m^3 sono le costanti di anisotropia del ferro a temperatura ambiente (298K298 K), e α1,α2,α3\alpha_1, \alpha_2, \alpha_3 rappresentano i coseni direttori del vettore di magnetizzazione M\mathbf{M} rispetto agli assi della cella elementare. Il disallineamento angolare medio della tessitura di Goss rispetto alla direzione di laminazione definisce la qualità dell'acciaio: negli acciai CRGO convenzionali la dispersione angolare oscilla tra 55^\circ e 77^\circ, mentre negli acciai ad alta permeabilità (Hi-B) si riduce a intervalli compresi tra 22^\circ e 33^\circ.

L'aggiunta di silicio (SiSi) in concentrazioni nominali dal 3.0%3.0\% al 3.4%3.4\% in peso aumenta la resistività elettrica volumetrica (ρ\rho) da circa 10×108 Ωm10 \times 10^{-8} \ \Omega\cdot m (per il ferro puro) fino a (4550)×108 Ωm(45 - 50) \times 10^{-8} \ \Omega\cdot m. Questo incremento di resistività mitiga drasticamente il trasporto libero di cariche parassite. Tuttavia, concentrazioni superiori al 3.5%3.5\% in peso inducono una severa fragilità meccanica a causa della formazione di fasi ordinate intermetalliche Fe3SiFe_3Si e FeSiFeSi (B2B2 e D03D0_3), impedendo i processi di formatura a freddo su scala industriale. Inoltre, il silicio riduce la costante di magnetostrizione di saturazione longitudinale (λ100\lambda_{100}), minimizzando le deformazioni elastiche indotte dal campo magnetico e la conseguente emissione di rumore acustico.

λ(H)=Δll=32λs((MMs)213)\lambda(H) = \frac{\Delta l}{l} = \frac{3}{2} \lambda_s \left( \left(\frac{M}{M_s}\right)^2 - \frac{1}{3} \right)

La struttura dei domini magnetici (domini di Weiss) nel CRGO è composta da domini principali orientati a 180180^\circ separati da pareti di Bloch, e da domini di chiusura a 9090^\circ in prossimità delle discontinuità superficiali e delle inclusioni. La larghezza caratteristica del dominio (2d02d_0) dipende dall'equilibrio tra l'energia magnetostatica superficiale e l'energia elastica della parete del dominio stesso. L'applicazione di rivestimenti dielettrici in trazione superficiale (come il rivestimento a base di fosfati e silicati, denominato commercialmente Carlite) introduce una trazione biassiale permanente nella lamiera (σt410MPa\sigma_t \approx 4 - 10 MPa), riducendo la larghezza media dei domini a 180180^\circ ed eliminando i domini secondari trasversali di chiusura, con una conseguente drastica diminuzione delle perdite per correnti parassite anomale.

Tripartizione Classica e Spettro delle Perdite: Il Modello di Bertotti

Sotto eccitazione magnetodinamica ciclica, la dissipazione volumetrica di energia per unità di tempo nel nucleo viene quantificata mediante l'integrale chiuso sul ciclo del ciclo di isteresi dinamico. La formulazione generale delle perdite specifiche di potenza (PtP_t, in W/kgW/kg) è governata dalla teoria statistica delle perdite elettromagnetiche formulata da Giorgio Bertotti, la quale supera le limitazioni delle formulazioni empiriche classiche di Steinmetz scomponendo il fenomeno in tre meccanismi dissipativi fisicamente disaccoppiati:

Ptotal(f,Bp)=Ph(f,Bp)+Pcl(f,Bp)+Pexc(f,Bp)P_{ total }(f, B_p) = P_h(f, B_p) + Pcl(f, B_p) + Pexc(f, B_p)
Ptotal=khfBpα+σπ2d26ρmf2Bp2+kexcρmf1.5Bp1.5P_{ total } = k_h \, f \, B_p^{\alpha} + \frac{\sigma \, \pi^2 \, d^2}{6 \, \rho_m} \, f^2 \, B_p^2 + \frac{kexc}{\rho_m} \, f^{1.5} \, B_p^{1.5}

Dove ff rappresenta la frequenza fondamentale di eccitazione (HzHz), BpB_p è il valore di picco dell'induzione magnetica (TT), σ\sigma è la conducibilità elettrica del materiale (S/mS/m), dd è lo spessore della lamiera magnetica (mm), ρm\rho_m è la densità di massa volumetrica dell'acciaio (kg/m3kg/m^3, tipicamente 7650kg/m37650 kg/m^3), khk_h è il coefficiente di isteresi quasi-statica, α\alpha è l'esponente di Steinmetz (frequentemente 1.6α2.01.6 \le \alpha \le 2.0), e kexckexc è il parametro microscopico delle perdite in eccesso.

Perdite per Isteresi Quasi-Statica

Le perdite per isteresi statica (PhP_h) hanno origine nei processi microscopici di dissipazione termodinamica irreversibile quando le pareti dei domini magnetici si spostano attraverso un mezzo ferromagnetico contenente difetti del reticolo cristallino, quali dislocazioni, vacanze, bordi di grano, tensioni residue e inclusioni non metalliche (ad es. precipitati di SiO2SiO_2 o Fe3CFe_3C). Durante il movimento di una parete di Bloch a 180180^\circ, essa rimane meccanicamente ancorata in questi centri di fissaggio (pinning sites).

Per liberare la parete di dominio dalla buca di potenziale di ancoraggio è necessario un incremento del campo magnetico esterno fino al raggiungimento del campo critico locale. Una volta superata questa barriera energetica, la parete avanza discontinuamente ad alta velocità verso la buca di potenziale successiva mediante un salto microscopico non lineare, fenomeno noto come effetto Barkhausen. L'energia dissipata in questo salto si converte irreversibilmente in fononi (calore) all'interno del reticolo cristallino:

wh=cicloHdB=μ0cicloHdM[J/m3]w_h = \oint_{ ciclo } \mathbf{H} \cdot d\mathbf{B} = \mu_0 \oint_{ ciclo } H \, dM \quad \left[ J/m^3\right]
Ph=fρmcicloHdB=khfBpα[W/kg]P_h = \frac{f}{\rho_m} \oint_{ ciclo } \mathbf{H} \cdot d\mathbf{B} = k_h f B_p^\alpha \quad \left[ W/kg \right]

L'area del ciclo di isteresi quasi-statico (whw_h) dipende esclusivamente dalla struttura metallurgica del materiale ed è indipendente dalla frequenza di variazione temporale del flusso per regimi in cui il tempo di rilassamento del dominio è infinitamente inferiore rispetto al periodo dell'eccitazione. La modellazione matematica di PhP_h in presenza di forme d'onda complesse o saturazione non lineare richiede la formulazione dell'operatore continuo di Preisach:

M(t)=αβμ(α,β)γ^αβ[H(t)]dαdβM(t) = \iint_{\alpha \ge \beta} \mu(\alpha, \beta) \, \hat{\gamma}_{\alpha\beta} [H(t)] \, d\alpha \, d\beta

Dove γ^αβ\hat{\gamma}_{\alpha\beta} è un operatore isteretico elementare bistabile con soglie di commutazione α\alpha (di accensione/innesco) e β\beta (di spegnimento/rilascio), e μ(α,β)\mu(\alpha, \beta) è la funzione di distribuzione di densità di Preisach specifica per il materiale CRGO in prova, identificabile sperimentalmente mediante famiglie di curve di inversione del primo ordine (FORC, First-Order Reversal Curves).

Perdite per Correnti Parassite Classiche (Foucault)

Le perdite classiche per correnti parassite (PclPcl) derivano rigorosamente dall'applicazione macroscopica delle equazioni di Maxwell a un conduttore omogeneo e isotropo continuo di geometria planare finita. Si consideri una lamiera di acciaio di spessore dd nella direzione dell'asse zz (da z=d/2z = -d/2 a z=+d/2z = +d/2), larghezza infinita lungo l'asse yy e lunghezza lungo l'asse xx. Assumendo una permeabilità magnetica effettiva μ\mu, una conducibilità volumetrica σ=1/ρ\sigma = 1/\rho, e un vettore di densità di flusso magnetico parallelo alla superficie: B(t)=Bpsin(ωt)x^\mathbf{B}(t) = B_p \sin(\omega t) \hat{\mathbf{x}}.

×E=Bt    Ey(z,t)z=Bx(z,t)t\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t} \implies \frac{\partial E_y(z,t)}{\partial z} = -\frac{\partial B_x(z,t)}{\partial t}

Ipotizzando che lo spessore dd sia sostanzialmente inferiore alla profondità di penetrazione elettromagnetica (skin depth) δ=2/(ωσμ)\delta = \sqrt{2/(\omega \sigma \mu)}, l'induzione magnetica si distribuisce in modo spazialmente uniforme attraverso la sezione trasversale della lamiera: Bx(z,t)Bpsin(ωt)B_x(z,t) \approx B_p \sin(\omega t). Integrando rispetto a zz a partire dal piano neutro centrale (z=0z = 0, dove Ey=0E_y = 0 per simmetria):

Ey(z,t)=zdBx(t)dt=zωBpcos(ωt)E_y(z,t) = -z \, \frac{d B_x(t)}{dt} = -z \, \omega B_p \cos(\omega t)

La densità di corrente indotta J(z,t)=σE(z,t)\mathbf{J}(z,t) = \sigma \mathbf{E}(z,t) dissipa una potenza istantanea per effetto Joule la cui densità volumetrica locale è p(z,t)=Jy2(z,t)/σ=σEy2(z,t)p(z,t) = J_y^2(z,t) / \sigma = \sigma E_y^2(z,t). Integrando sull'intero spessore della lamiera e mediando nel dominio del tempo durante un periodo T=2π/ωT = 2\pi/\omega:

Pcl,v=1dd/2d/2[1T0Tσ(zωBpcos(ωt))2dt]dz\langle P_{cl,v} \rangle = \frac{1}{d} \int_{-d/2}^{d/2} \left[ \frac{1}{T} \int_0^T \sigma \left( -z \omega B_p \cos(\omega t) \right)^2 dt \right] dz
Pcl,v=σω2Bp22dd/2d/2z2dz=σω2Bp22d[z33]d/2d/2=σω2Bp2d224=σπ2f2d2Bp26[W/m3]\langle P_{cl,v} \rangle = \frac{\sigma \, \omega^2 B_p^2}{2 d} \int_{-d/2}^{d/2} z^2 \, dz = \frac{\sigma \, \omega^2 B_p^2}{2 d} \left[ \frac{z^3}{3} \right]_{-d/2}^{d/2} = \frac{\sigma \, \omega^2 B_p^2 \, d^2}{24} = \frac{\sigma \, \pi^2 f^2 d^2 B_p^2}{6} \quad \left[ W/m^3\right]

Dividendo per la densità di massa volumetrica ρm\rho_m si ottiene la formulazione specifica normalizzata:

Pcl=σπ2d26ρmf2Bp2=π2d26ρρmf2Bp2[W/kg]Pcl = \frac{\sigma \, \pi^2 \, d^2}{6 \, \rho_m} \, f^2 \, B_p^2 = \frac{\pi^2 \, d^2}{6 \, \rho \, \rho_m} \, f^2 \, B_p^2 \quad \left[ W/kg \right]

Questa derivazione evidenzia la dipendenza quadratica diretta dallo spessore dd della lamiera e dalla frequenza fondamentale ff, giustificando la necessità tecnica di ridurre gli spessori commerciali da 0.35mm0.35 mm (M4) a 0.23mm0.23 mm o 0.18mm0.18 mm (gradi laser a bassissime perdite).

Perdite Anomale o in Eccesso (Excess Losses)

Storicamente, la somma diretta delle perdite per isteresi quasi-statica e delle perdite classiche per correnti parassite risultava significativamente inferiore rispetto alle perdite totali misurate sperimentalmente mediante il metodo del wattmetro su telaio di Epstein (Ptotal,mis>Ph+PclP_{ total, mis } > P_h + Pcl). Questa discrepanza energetica è stata formalmente denominata anomalia delle perdite.

Il modello microscopico di Pry e Bean (1958) ha dimostrato che il flusso magnetico non varia uniformemente all'interno della lamiera considerata come un continuo, ma si concentra esclusivamente all'interno dei domini magnetici e si trasferisce attraverso lo spostamento localizzato delle pareti di Bloch a 180180^\circ. La velocità di spostamento di una parete di dominio vw(t)v_w(t) sotto eccitazione sinusoidale induce un gradiente di campo elettrico locale infinitamente più marcato nelle immediate vicinanze della parete rispetto a quanto previsto dalla teoria classica del continuo:

vw(t)=L2BsdB(t)dtv_w(t) = \frac{L}{2 B_s} \frac{dB(t)}{dt}

Dove LL rappresenta la spaziatura tra le pareti di dominio e BsB_s è l'induzione di saturazione. La micro-corrente parassita confinata nel micro-volume adiacente alla parete genera un campo di frenatura magnetica locale. Bertotti ha generalizzato questo fenomeno fisico formulando la teoria statistica degli Oggetti Magnetici Attivi (Magnetic Objects, MO), intesi come aggregati coerenti di pareti di dominio che interagiscono con le eterogeneità locali del mezzo:

Pexc=kexcρmf1.5Bp1.5=1ρmσGSV0(fBp)1.5[W/kg]Pexc = \frac{kexc}{\rho_m} \, f^{1.5} \, B_p^{1.5} = \frac{1}{\rho_m} \sqrt{\sigma \, G \, S \, V_0} \, (f \, B_p)^{1.5} \quad \left[ W/kg \right]

Dove G=0.1356G = 0.1356 è una costante adimensionale associata alla dissipazione del campo di Foucault in mezzi infiniti, SS è l'area della sezione trasversale della lamiera perpendicolare alla direzione del flusso magnetico, e V0V_0 è un parametro intrinseco del materiale avente le dimensioni di un campo magnetico (A/mA/m), il quale quantifica l'interazione stocastica tra gli oggetti magnetici e la distribuzione del campo microscopico di ancoraggio (internal pinning field distribution).

Grado CRGO (Classificazione AISI / EN) Spessore Nominale dd (mm) Induzione Massima BB a 800A/m800 A/m (T) Perdite Specifiche P1.7/50P_{1.7/50} (W/kg) Perdite Specifiche P1.7/60P_{1.7/60} (W/kg) Fattore di Impaccamento (\%)
M4 (AISI 35G155 / EN 10107) 0.35 1.82 - 1.84 1.25 - 1.35 1.64 - 1.77 96.0\ge 96.0
M3 (AISI 30G130 / EN 10107) 0.30 1.84 - 1.86 1.05 - 1.15 1.38 - 1.51 95.5\ge 95.5
Hi-B Convenzionale (27Q110) 0.27 1.90 - 1.93 0.95 - 1.05 1.24 - 1.38 95.0\ge 95.0
Hi-B Alta Permeabilità (23Q090) 0.23 1.92 - 1.95 0.80 - 0.90 1.05 - 1.18 94.5\ge 94.5
Hi-B Laser-Scribed (20Q080-L) 0.20 1.93 - 1.96 0.70 - 0.78 0.92 - 1.02 94.0\ge 94.0
Metallo Amorfo (Fe-Si-B Metglas 2605SA1) 0.025 1.56 - 1.58 0.18 - 0.25 0.24 - 0.33 86.0\ge 86.0

Effetti Non Lineari, Distorsione Armonica e Polarizzazione in Corrente Continua (DC Bias)

Nei sistemi elettrici di potenza contemporanei, la diffusione massiva dell'elettronica di potenza (inverter su larga scala per impianti fotovoltaici, collegamenti HVDC, sistemi FACTS e azionamenti a frequenza variabile) introduce regimi di magnetizzazione non sinusoidali e tensioni con elevato contenuto armonico (THDvTHD_v). La risposta non lineare della curva di magnetizzazione dinamica B(H)B(H) trasforma queste sollecitazioni in severe perdite energetiche anomale.

Eccitazione Armonica e Cicli di Isteresi Minori

Quando la forma d'onda di tensione applicata contiene armoniche superiori (h=3,5,7,h = 3, 5, 7, \dots), la derivata temporale dell'induzione magnetica B/t\partial \mathbf{B}/\partial t presenta molteplici passaggi per lo zero nell'arco di un singolo ciclo fondamentale. Fisicamente, questo impone l'inversione transitoria del movimento delle pareti dei domini magnetici prima del completamento della polarizzazione macroscopica totale, originando cicli di isteresi minori (minor hysteresis loops) innestati all'interno del ciclo principale di saturazione.

Pcl,arm=σπ2d26ρmh=1(hf1)2Bp,h2=Pcl,1h=1h2(Bp,hBp,1)2P_{cl, arm } = \frac{\sigma \, \pi^2 \, d^2}{6 \, \rho_m} \sum_{h=1}^{\infty} \left( h \, f_1 \right)^2 B_{p,h}^2 = P_{cl,1} \sum_{h=1}^{\infty} h^2 \left( \frac{B_{p,h}}{B_{p,1}} \right)^2

Sotto tensione puramente sinusoidale, BhVh/(hf1)B_h \propto V_h / (h f_1). Se la tensione presenta una distorsione armonica con componenti VhV_h, il campo armonico indotto si approssima con Bp,hVhhV1Bp,1B_{p,h} \approx \frac{V_h}{h \, V_1} B_{p,1}. Sostituendo nella formulazione delle correnti parassite classiche:

Pcl,total=Pcl,1(1+h=2(VhV1)2)=Pcl,1(1+THDv2)P_{cl, total } = P_{cl,1} \left( 1 + \sum_{h=2}^{\infty} \left( \frac{V_h}{V_1} \right)^2 \right) = P_{cl,1} \left( 1 + THD_v^2 \right)

Tuttavia, le perdite in eccesso non seguono questo disaccoppiamento lineare. L'interazione tra frequenze multiple accelera la velocità locale delle pareti di Bloch secondo un integrale temporale della derivata del flusso:

Pexc,nonsin=kexcρm1T0TdB(t)dt1.5dtP_{exc, non-sin } = \frac{kexc}{\rho_m} \frac{1}{T} \int_0^T \left| \frac{d B(t)}{dt} \right|^{1.5} dt

Magnetizzazione con Polarizzazione in Corrente Continua (DC Bias)

L'iniezione di corrente continua nelle fasi del trasformatore — causata da correnti geomagnetizzanti indotte (GIC), topologie di inverter ad accoppiamento diretto o guasti asimmetrici monopolari su linee di trasmissione HVDC con ritorno a terra — produce una traslazione del ciclo dinamico lungo l'asse del campo magnetizzante (Hdc=NIdc/lfeHdc = N Idc / lfe).

B(t)=Bdc+ΔB(t)=μrev(Hdc)Hdc+B^sin(ωt)B(t) = Bdc + \Delta B(t) = \mu_{rev}(Hdc) \cdot Hdc + \hat{B} \sin(\omega t)

Questo spostamento asimmetrico spinge il nucleo verso il ginocchio di saturazione durante un intero semiciclo (saturazione semi-ciclica). Nella regione di saturazione (B>1.9TB > 1.9 T nel CRGO), la permeabilità differenziale μd=dB/dH\mu_d = dB/dH decade drasticamente verso valori prossimi a quelli della permeabilità del vuoto (μ0=4π×107H/m\mu_0 = 4\pi \times 10^{-7} H/m). Ciò provoca massicci picchi di corrente magnetizzante asimmetrica (imag(t)imag(t) con fattori di cresta superiori a 10), trasformando il campo magnetico da un vettore confinato nel nucleo laminato in un flusso di dispersione tridimensionale di grande ampiezza che concatena gli elementi strutturali adiacenti (piastre di serraggio, cassa del trasformatore, tiranti e schermi magnetici).

Qdc=ωimag(t)dλ(t)ω0TH(t)dBdtdtQdc = \omega \oint imag(t) \, d\lambda(t) \propto \omega \int_0^T H(t) \frac{dB}{dt} \, dt

Il consumo di potenza reattiva QQ aumenta di ordini di grandezza, provocando severe cadute di tensione ai morsetti, la comparsa di armoniche pari e dispari di basso ordine (in particolare h=2h = 2 e h=3h = 3), e surriscaldamenti termici critici dovuti a correnti indotte nelle parti metalliche strutturali.

Analisi Forense dei Guasti Termici e Degradazione Dielettrica Indotta dal Nucleo

Le patologie elettromagnetiche e la degradazione termica che si sviluppano nel nucleo costituiscono un meccanismo di guasto critico nei trasformatori di potenza di grande taglia. A differenza dei guasti meccanici negli avvolgimenti causati dalle sollecitazioni radiali e assiali di cortocircuito, i guasti al nucleo si manifestano tipicamente come processi di degradazione elettrochimica e termica progressiva che evolvono lungo estesi periodi di servizio.

Degradazione dell'Isolamento Interlaminare (Rivestimento Tipo Carlite)

L'isolamento interlaminare superficiale dei lamierini magnetici (spessore 1.53.0 μm\sim 1.5 - 3.0 \ \mu m, composto da una base di forsterite Mg2SiO4Mg_2SiO_4 con strato superiore di fosfato di alluminio AlPO4AlPO_4 e acido fosforico) è progettato per sopportare tensioni interlaminari di pochi volt (da 1.5V1.5 V a 5.0V5.0 V efficaci secondo la norma IEC 60404-6). Tuttavia, sotto sollecitazioni termiche locali continue (T>150CT > 150^\circ C), fatica vibrazionale da magnetostrizione (100/120Hz100/120 Hz e armoniche acustiche) o eccessive tensioni meccaniche da serraggio durante l'assemblaggio dei gioghi, il rivestimento dielettrico si frammenta e si polverizza meccanicamente.

La perdita dell'isolamento galvanico tra lamierini contigui provoca micro-saldature da arco o contatti galvanici diretti. Se si stabiliscono due o più punti di contatto nello stesso pacco di laminazione, si forma una spira in cortocircuito di ampia sezione che concatena il flusso magnetico principale del nucleo. La forza elettromotrice indotta nella spira chiusa risulta:

Eloop=SBtdS=ωBpAloopcos(ωt)\mathcal{E}_{loop} = -\oint_{\partial S} \frac{\partial \mathbf{B}}{\partial t} \cdot d\mathbf{S} = -\omega \, B_p \, Aloop \cos(\omega t)
Ifault=EloopRloop2+(ωLloop)2EloopRloopIfault = \frac{\mathcal{E}_{loop}}{\sqrt{Rloop^2 + (\omega Lloop)^2}} \approx \frac{\mathcal{E}_{loop}}{Rloop}

Poiché la resistenza dell'anello ohmico RloopRloop è nell'ordine dei milliohm (mΩm\Omega), la corrente di guasto IfaultIfault raggiunge ampiezze che vanno da centinaia a migliaia di ampere in regime stazionario, dissipando calore localizzato per effetto Joule che genera punti caldi (hot-spots) con temperature che superano rapidamente 500C800C500^\circ C - 800^\circ C.

Degradazione dell'Olio Dielettrico e Analisi Forense mediante DGA

Il trasferimento dell'estremo calore generato dal punto caldo del nucleo al fluido dielettrico circostante (olio minerale naftenico, paraffinico o esteri sintetici) induce la pirolisi e la scissione termocatalitica delle catene molecolari di idrocarburi (CnH2n+2C_n H_{2n+2}). I legami covalenti carbonio-idrogeno (CHC-H, energia di legame 413kJ/mol\sim 413 kJ/mol) e carbonio-carbonio (CCC-C, 348kJ/mol\sim 348 kJ/mol; C=CC=C, 614kJ/mol\sim 614 kJ/mol; CCC \equiv C, 839kJ/mol\sim 839 kJ/mol) si rompono in funzione del gradiente termico di contatto.

In conformità agli standard internazionali IEC 60599 e IEEE C57.104, il profilo dei gas disciolti nell'olio (Dissolved Gas Analysis, DGA) consente l'identificazione termodinamica precisa dell'intervallo di temperatura del guasto nel nucleo:

Intervallo Termico del Guasto Gas Combustibile Dominante Rapporti Diagnostici (IEC 60599 / Duval) Meccanismo Fisico nel Nucleo Azione Critica Richiesta
Guasto Termico T<300CT < 300^\circ C (T1) Metano (CH4CH₄) / Etano (C2H6C₂H₆) C2H2C2H4<0.1\frac{C₂H₂}{C₂H₄} < 0.1, 0.1<CH4H2<10.1 < \frac{CH₄}{H₂} < 1, C2H4C2H6<1\frac{C₂H₄}{C₂H₆} < 1 Riscaldamento diffuso dei gioghi, degradazione iniziale delle vernici interlaminari o della pittura di bloccaggio. Monitoraggio cromatografico periodico mensile; verifica della presenza di armoniche in rete.
Guasto Termico 300C<T<700C300^\circ C < T < 700^\circ C (T2) Etilene (C2H4C₂H₄) C2H2C2H4<0.1\frac{C₂H₂}{C₂H₄} < 0.1, CH4H2>1\frac{CH₄}{H₂} > 1, 1<C2H4C2H6<41 < \frac{C₂H₄}{C₂H₆} < 4 Spire in cortocircuito tra molteplici lamierini di acciaio CRGO. Severo degrado del rivestimento inorganico. Riduzione del carico al 70%70\%; esecuzione della prova di perdite a vuoto e verifica della resistenza di isolamento del nucleo.
Guasto Termico T>700CT > 700^\circ C (T3) Etilene (C2H4C₂H₄) con tracce di Acetilene (C2H2C₂H₂) C2H2C2H4<0.15\frac{C₂H₂}{C₂H₄} < 0.15, CH4H2>1\frac{CH₄}{H₂} > 1, C2H4C2H6>4\frac{C₂H₄}{C₂H₆} > 4 Fusione localizzata del pacco magnetico (core burning); cortocircuito tra lamierini e struttura di serraggio. Fuori servizio immediato forzato. Ispezione interna endoscopica e ricostruzione del nucleo.
Scariche Parziali / Scintillio nel Nucleo Idrogeno (H2H₂) dominante con CH4CH₄ CH4H2<0.1\frac{CH₄}{H₂} < 0.1, H2H₂ oltre i limiti basali (> 100 ppm) Flottazione galvanica del nucleo (perdita della terra fisica singola), che genera scariche capacitive verso la cassa. Verifica del circuito di messa a terra esterna del nucleo e del neutro mediante megaohmmetro a 2.5kV2.5 kV.

Guasto per Messa a Terra Multipla (Correnti Circolanti da Anello di Terra)

Per prescrizione normativa, il nucleo di un trasformatore deve essere collegato a terra galvanicamente in un solo punto fisico mediante un piatto di rame collegato a un apposito isolatore passante dedicato sul coperchio della cassa. Questa messa a terra singola previene l'innalzamento del potenziale elettrostatico flottante indotto dall'accoppiamento capacitivo con gli avvolgimenti di alta tensione (Vind=Cwc/(Cwc+Ccg)VwindingV_{ ind } = C_{w-c}/(C_{w-c} + C_{c-g}) \cdot Vwinding).

Qualora, a causa del degrado dell'isolamento dei tiranti (tie-rods), del cortocircuito provocato da depositi metallici sul fondo della cassa o di errori di montaggio, si verifichi una seconda connessione a terra accidentale all'estremità opposta del nucleo, si crea un circuito chiuso di ampia area fisica che concatena il flusso magnetico principale:

Vloop=ddtareatraleterreBdSV_{ loop } = -\frac{d}{dt} \iint_{ area tra le terre } \mathbf{B} \cdot d\mathbf{S}

Questa tensione induce una corrente di circolazione permanente attraverso la struttura della cassa e i lamierini di terra del nucleo, tipicamente compresa tra 10A10 A e oltre 150A150 A, provocando la bruciatura dell'isolamento dei tiranti passanti, la carbonizzazione dell'olio adiacente e concentrazioni allarmanti di etilene (C2H4C₂H₄) e metano (CH4CH₄), senza che le protezioni elettriche di massima corrente o differenziali (87T87 T) rilevino alcuna anomalia nelle correnti di linea.

Tecniche Avanzate di Mitigazione nella Progettazione e Fabbricazione

L'ottimizzazione delle prestazioni elettromagnetiche dei nuclei magnetici richiede un approccio multidisciplinare che comprende la modifica cristallografica del materiale, la microingegneria dei domini tramite radiazione coerente e l'ottimizzazione topologica tridimensionale delle giunzioni nei gioghi.

Giunzioni di Giogo di Tipo Step-Lap a 4545^\circ

L'assemblaggio geometrico delle colonne e dei gioghi nei nuclei a pacco introduce discontinuità nel circuito magnetico. Nelle tradizionali giunzioni testa a testa a sovrapposizione semplice a 9090^\circ (Butt-Lap joints), il flusso magnetico è costretto ad attraversare i lamierini perpendicolarmente alla direzione di laminazione (lungo la direzione cristallografica dura 110\langle 110 \rangle), causando una severa rifrazione delle linee di flusso, un drastico aumento della riluttanza magnetica d'interfaccia e violente saturazioni locali sui bordi di taglio.

La moderna tecnologia di taglio e impilamento impiega giunzioni a 4545^\circ con sfalsamento progressivo dei lamierini denominato Step-Lap (generalmente configurazioni da 5 a 7 gradini per ciclo). La gradinatura distribuisce il traferro di taglio longitudinalmente lungo un gradiente assiale:

Rgap=gμ0Ajoint\mathcal{R}_{gap} = \frac{g}{\mu_0 \cdot Ajoint}

Sfalsando i giunti di una distanza spaziale Δx3.57.0mm\Delta x \approx 3.5 - 7.0 mm tra lamierini contigui, il flusso magnetico non è costretto a saltare perpendicolarmente attraverso un unico traferro ad alta riluttanza verso il lamierino adiacente, ma si distribuisce progressivamente lungo una sezione geometrica effettiva equivalente notevolmente maggiore:

Aeff=Agap+k=1NstepskΔxwsheetAeff = Agap + \sum_{k=1}^{Nsteps} k \cdot \Delta x \cdot wsheet

L'adozione della configurazione Step-Lap riduce le perdite totali di giunzione nei gioghi fino al 25%30%25\% - 30\%, mitiga la potenza reattiva di magnetizzazione del 40%40\% e diminuisce i livelli di emissione acustica del trasformatore da 3a6dB(A)3 a 6 dB(A) rispetto agli assemblaggi convenzionali.

Affinamento dei Domini Magnetici mediante Raggio Laser (Laser Scribing)

Negli acciai CRGO ad alta permeabilità (Hi-B), la dimensione dei grani cristallografici è notevole (1030mm\sim 10 - 30 mm), il che favorisce un orientamento di Goss quasi perfetto (3\le 3^\circ), ma risulta controproducente per il comportamento dinamico delle perdite poiché la spaziatura tra le pareti di Bloch a 180180^\circ (2d02d_0) cresce proporzionalmente alla radice quadrata della dimensione del grano:

2d0=γwLgrain1.38μ0Ms22d_0 = \sqrt{\frac{\gamma_w \cdot Lgrain}{1.38 \cdot \mu_0 \cdot M_s^2}}

Un'elevata larghezza dei domini incrementa drasticamente la velocità istantanea di spostamento della parete vw(t)v_w(t), facendo impennare le perdite anomale per correnti parassite (PexcPexc).

La tecnica di affinamento dei domini mediante laser (Domain Refinement o Laser Scribing) proietta un raggio focalizzato di laser continuo (Nd:YAG o a fibra ad alta potenza) trasversalmente alla direzione di laminazione a intervalli spaziali regolari (p48mmp \approx 4 - 8 mm). L'energia termica dell'impulso laser induce un riscaldamento ultrarapido e localizzato senza asportazione di materiale, generando un micro-campo di tensioni elastiche di compressione permanente sub-superficiale dovuto al gradiente termico di raffreddamento:

σresidual(y,z)=EαthΔT(y,z)+1AEαthΔTdydz\sigma_{ residual }(y,z) = -E \, \alpha_{th} \, \Delta T(y,z) + \frac{1}{A} \int \int E \, \alpha_{th} \, \Delta T \, dy \, dz

Queste bande periodiche di tensione interna agiscono come barriere artificiali di energia magnetoelastica che frammentano i larghi domini magnetici a 180180^\circ, riducendo la loro spaziatura media 2d02d_0 alla metà o a un terzo della dimensione originale. Come conseguenza diretta:

Pexc,lasered=Pexc,base(d0,laseredd0,base)β    ΔPtotal10%a15%P_{exc, lasered } = P_{exc, base } \left( \frac{d_{0, lasered }}{d_{0, base }} \right)^\beta \implies \Delta P_{ total } \approx -10\% a -15\%

È imperativo distinguere tra il trattamento laser di tipo non termicamente resistente (che perde l'effetto di affinamento dopo ricotture di distensione a T>500CT > 500^\circ C richieste per i nuclei avvolti) e i metodi meccanici di incisione chimica/deformazione plastica controllata al plasma (ricuocibili), progettati per mantenere l'ancoraggio dei microdomini anche a seguito di trattamenti termici a 800C800^\circ C.

Implementazione e Modellazione Computazionale con Vexten Suite

La quantificazione rigorosa delle perdite nel nucleo e il loro accoppiamento multifisico con la rete elettrica vengono risolti operativamente in Vexten Suite integrando i moduli di calcolo dei flussi armonici non lineari (secondo IEEE 519 / IEEE C57.110), il calcolo termico dinamico accoppiato (IEC 60076-7) e la simulazione dei transitori di cortocircuito e magnetici (IEC 60909 / IEEE 141).

Algoritmo di Correzione per Distorsione Armonica e Fattore K

In presenza di spettri armonici gravosi originati da carichi non lineari o generatori da fonti rinnovabili, il modulo Vexten PowerFlow & Harmonics esegue una ripartizione dinamica delle perdite ricalcolando il fattore di perdita del nucleo in base al tensore dei coefficienti spettrali normalizzati:

FNLLoss=Pcore,distPcore,sin=khh=1Hmax(VhhV1)α+kclh=1Hmax(VhV1)2+kexch=1Hmaxh0.5(VhV1)1.5F_{NL-Loss} = \frac{P_{core, dist }}{P_{core, sin }} = k_h \sum_{h=1}^{Hmax} \left(\frac{V_h}{h V_1}\right)^\alpha + kcl \sum_{h=1}^{Hmax} \left(\frac{V_h}{V_1}\right)^2 + kexc \sum_{h=1}^{Hmax} h^{0.5} \left(\frac{V_h}{V_1}\right)^{1.5}

Il software collega questo fattore al declassamento nominale (derating) del trasformatore mediante il calcolo del Fattore KK e del Fattore Armonico di Perdita per Correnti Parassite (FHLFHL), regolando la massima corrente ammissibile per garantire che la temperatura del punto più caldo dell'avvolgimento e del nucleo non superi i limiti imposti dalla classe termica di isolamento:

Imax,pu=Ptotal,nomPcore,dist(Vthd)Pdc,nom(1+FHLPec,nompu+FHLOSLPosl,nompu)I_{max,pu} = \sqrt{\frac{P_{total,nom} - P_{core, dist }(Vthd)}{P_{dc,nom} \left( 1 + FHL \cdot P_{ec,nom-pu} + F_{HL-OSL} \cdot P_{osl,nom-pu} \right)}}

Procedura di Modellazione Multifisica in Vexten Suite

  1. Definizione dello Spettro di Eccitazione Magnetica: Nel modulo Power Quality, importare le misurazioni di campo o gli spettri tipici IEEE delle armoniche di tensione (VhV_h) e la componente continua residua (IdcIdc) derivante dal flusso di carico accoppiato con impianti a fonti rinnovabili.
  2. Assegnazione del Materiale Magnetico nel Database: Selezionare nelle proprietà della macchina il grado specifico di acciaio CRGO (es. M3, 23Q090 o 20Q080-L con affinamento laser). Vexten Suite carica automaticamente i parametri microscopici calibrati: conducibilità σ\sigma, coefficiente di isteresi khk_h, spessore dd, densità ρm\rho_m e il parametro di Bertotti V0V_0.
  3. Configurazione della Topologia di Giunzione: Selezionare nei parametri costruttivi la tipologia di giunzione dei gioghi: Butt-Lap a 9090^\circ, Step-Lap a 5 gradini o Step-Lap a 7 gradini. Il software introduce automaticamente il moltiplicatore di dispersione della riluttanza d'interfaccia (Kjoint\mathcal{K}_{joint}).
  4. Esecuzione dell'Analisi Termica Transitoria Accoppiata: Il solutore accoppia le perdite elettromagnetiche totali Ptotal(x,y,z)P_{ total }(x,y,z) come sorgente volumetrica di calore all'interno del risolutore di reti termiche conforme a IEC 60076-7, determinando la distribuzione termica tridimensionale nei gioghi superiori, inferiori e nelle colonne del trasformatore.
  5. Validazione dei Limiti e Diagnosi Forense Preventiva: Il sistema confronta i gradienti di temperatura locale e le densità di perdita rispetto alle soglie stabilite dallo standard IEEE C57.104. Qualora la temperatura nel nucleo superi continuativamente i 130C130^\circ C, si attivano automaticamente allarmi di declassamento della potenza e viene calcolata la velocità teorica stimata di generazione dei gas di guasto (C2H4C₂H₄, CH4CH₄, H2H₂), consentendo la definizione di piani di manutenzione predittiva basata sulle condizioni reali (CBM) di elevata precisione fisica.

Grazie a questa architettura computazionale integrata, i progettisti, gli ingegneri di collaudo e gli specialisti di affidabilità possono prevedere con esattezza l'invecchiamento termico del circuito magnetico, eliminare il rischio di guasti catastrofici causati da correnti parassite interlaminari e ottimizzare la scelta di materiali magnetici avanzati destinati a trasformatori a rendimento energetico ultra-elevato.