Hysterese und Wirbelstromverluste in CRGO-Siliziumstahl

¿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

Elektromagnetische Thermodynamik und Mikrostruktur von kornorientiertem Silizium-Elektroblech (CRGO)

Der Magnetkern eines Leistungstransformators fungiert als grundlegender Flusskopplungskreis zwischen den Wicklungen und arbeitet unter nichtlinearen zyklischen elektromagnetischen Erregungsregimen. Der Gesamtwirkungsgrad und die Leistungsdichte des Transformators werden maßgeblich durch das mikrostrukturelle und kristallographische Verhalten des verwendeten ferromagnetischen Materials bestimmt. Der vorherrschende Industriestandard für Mittel-, Hoch- und Höchstspannungstransformatoren (EHV/UHV) ist kaltgewalztes kornorientiertes Elektroblech (Cold-Rolled Grain-Oriented Silicon Steel, CRGO), das durch thermomechanische Behandlungen optimiert wird, um seine kristallographischen Vorzugsachsen der leichten Magnetisierung parallel zur Walzrichtung auszurichten.

CRGO-Stahl weist eine raumzentrierte kubische Kristallstruktur (BCC, α\alpha-Fe-Gitter) auf, in der die kristallographische 100\langle 100 \rangle-Achse die Richtung minimaler magnetokristalliner Anisotropieenergie darstellt, während die 110\langle 110 \rangle- und 111\langle 111 \rangle-Richtungen jeweils Achsen mittlerer bzw. schwerer Anisotropie abbilden. Durch eine präzise Sequenz von Kaltwalzungen und sekundären Rekristallisationsglühungen (gehemmt durch Ausscheidungen von Mangansulfid, MnSMnS, oder Aluminiumnitrid, AlNAlN) bildet sich die sogenannte Goss-Textur heraus, die durch die kristallographische Orientierung {110}001\{110\}\langle 001 \rangle gekennzeichnet ist. In dieser Konfiguration liegen die {110}\{110\}-Ebenen parallel zur Blechoberfläche und die 001\langle 001 \rangle-Richtungen richten sich exakt parallel zur Längswalzrichtung (RD, Rolling Direction) aus.

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)

Hierbei ist EaE_a die magnetokristalline Anisotropieenergiedichte (J/m3J/m^3), K14,8×104 J/m3K_1 \approx 4{,}8 \times 10^4 \ J/m^3 und K21,5×104 J/m3K_2 \approx 1{,}5 \times 10^4 \ J/m^3 sind die Anisotropiekonstanten von Eisen bei Raumtemperatur (298 K298 \ K), und α1,α2,α3\alpha_1, \alpha_2, \alpha_3 repräsentieren die Richtungskosinus des Magnetisierungsvektors M\mathbf{M} bezüglich der Achsen der Elementarzelle. Die mittlere Winkelabweichung der Goss-Textur von der Walzrichtung definiert die Güteklasse des Stahls: Bei konventionellen CRGO-Stählen liegt die Winkelstreuung zwischen 55^\circ und 77^\circ, während sie bei hochpermeablen Stählen (Hi-B) auf Bereiche von 22^\circ bis 33^\circ reduziert ist.

Die Zulegierung von Silizium (SiSi) in nominalen Konzentrationen von 3,0%3{,}0\% bis 3,4%3{,}4\% Massenanteil erhöht den spezifischen elektrischen Volumenwiderstand (ρ\rho) von etwa 10×108 Ωm10 \times 10^{-8} \ \Omega\cdot m (für Reinstfeisen) auf (4550)×108 Ωm(45 - 50) \times 10^{-8} \ \Omega\cdot m. Diese Erhöhung des spezifischen Widerstands dämpft den freien Transport parasitärer Wirbelströme drastisch. Siliziumkonzentrationen über 3,5%3{,}5\% Massenanteil führen jedoch zu gravierender mechanischer Versprödung infolge der Bildung intermetallischer Ordnungsphasen wie Fe3SiFe_3Si und FeSiFeSi (B2B2 und D03D0_3), was die Kaltumformprozesse im industriellen Maßstab verhindert. Darüber hinaus senkt Silizium die longitudinale Sättigungsmagnetostriktionskonstante (λ100\lambda_{100}), wodurch die durch das Magnetfeld induzierten elastischen Verformungen und die daraus resultierende Schallemission minimiert werden.

λ(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)

Die magnetische Domänenstruktur (Weiss-Bezirke) in CRGO besteht aus 180180^\circ-Hauptdomänen, die durch Bloch-Wände getrennt sind, sowie aus 9090^\circ-Abschlussdomänen in der Nähe von Oberflächendiskontinuitäten und Einschlüssen. Die charakteristische Domänenbreite (2d02d_0) resultiert aus dem Gleichgewicht zwischen der magnetostatischen Oberflächenenergie und der elastischen Domänenwandenergie. Das Aufbringen dielektrischer Beschichtungen unter Oberflächenspannung (wie die Phosphat- und Silikatbeschichtung, kommerziell als Carlite bezeichnet) erzeugt eine permanente biaxiale Zugspannung im Blech (σt410 MPa\sigma_t \approx 4 - 10 \ MPa). Dies verringert die mittlere Breite der 180180^\circ-Domänen, eliminiert transversale sekundäre Abschlussdomänen und senkt die anomalen Wirbelstromverluste radikal.

Klassische Verlustdreiteilung und Verlustspektrum: Das Bertotti-Modell

Unter zyklischer magnetodynamischer Erregung wird die volumetrische Energiedissipation pro Zeiteinheit im Kern durch das geschlossene Schleifenintegral der dynamischen Hysteresekurve quantifiziert. Die allgemeine Formulierung der spezifischen Verlustleistung (PtP_t, in W/kgW/kg) folgt der statistischen Theorie elektromagnetischer Verluste nach Giorgio Bertotti. Dieses Modell überwindet die Restriktionen der klassischen empirischen Steinmetz-Gleichungen, indem es das Phänomen in drei physikalisch entkoppelte Dissipationsmechanismen zerlegt:

Ptotal(f,Bp)=Ph(f,Bp)+Pcl(f,Bp)+Pexc(f,Bp)Ptotal(f, B_p) = P_h(f, B_p) + Pcl(f, B_p) + Pexc(f, B_p)
Ptotal=khfBpα+σπ2d26ρmf2Bp2+kexcρmf1.5Bp1.5Ptotal = 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}

Dabei repräsentiert ff die Grundfrequenz der Erregung (HzHz), BpB_p den Spitzenwert der magnetischen Induktion (TT), σ\sigma die elektrische Leitfähigkeit des Materials (S/mS/m), dd die Dicke des lamellierten Blechs (mm), ρm\rho_m die Massendichte des Stahls (kg/m3kg/m^3, typischerweise 7650 kg/m37650 \ kg/m^3), khk_h den quasistatischen Hysteresekoeffizienten, α\alpha den Steinmetz-Exponenten (häufig 1,6α2,01{,}6 \le \alpha \le 2{,}0) und kexckexc den mikroskopischen Überschussverlustparameter.

Quasistatische Hystereseverluste

Statische Hystereseverluste (PhP_h) entstehen durch mikroskopische irreversible thermodynamische Dissipationsprozesse, wenn sich magnetische Domänenwände durch ein ferromagnetisches Medium mit Kristallgitterdefekten wie Versetzungen, Leerstellen, Korngrenzen, Eigenspannungen und nichtmetallischen Einschlüssen (z. B. SiO2SiO_2- oder Fe3CFe_3C-Ausscheidungen) bewegen. Während der Verschiebung einer 180180^\circ-Bloch-Wand verankert sich diese mechanisch an diesen Haftzentren (Pinning Sites).

Um die Domänenwand aus der Verankerungspotenzialmulde zu befreien, ist eine Erhöhung des äußeren Magnetfelds bis zum Erreichen der lokalen kritischen Feldstärke erforderlich. Sobald diese Energiebarriere überwunden ist, springt die Wand unstetig mit hoher Geschwindigkeit in die nächste Potenzialmulde über — ein mikroskopischer, nichtlinearer Sprungprozess, der als Barkhausen-Effekt bekannt ist. Die bei diesem Sprung dissipierte Energie wandelt sich irreversibel in Phononen (Gitterwärme) um:

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

Die Fläche der quasistatischen Hystereseschleife (whw_h) hängt ausschließlich vom metallurgischen Gefüge des Werkstoffs ab und ist unabhängig von der zeitlichen Änderungsrate des Flusses für Betriebsbereiche, in denen die Relaxationszeit der Domänen unendlich viel kleiner ist als die Periode der Erregung. Die mathematische Modellierung von PhP_h bei komplexen Signalformen oder nichtlinearer Sättigung erfordert die Formulierung des kontinuierlichen Preisach-Operators:

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

Wobei γ^αβ\hat{\gamma}_{\alpha\beta} ein elementarer bistabiler Hystereseoperator mit Schaltschwellen α\alpha (Einschalten) und β\beta (Ausschalten) ist und μ(α,β)\mu(\alpha, \beta) die spezifische Preisach-Dichteverteilungsfunktion des untersuchten CRGO-Materials darstellt, die experimentell über Familien von Umkehrkurven erster Ordnung (FORC, First-Order Reversal Curves) identifiziert werden kann.

Klassische Wirbelstromverluste (Foucault)

Klassische Wirbelstromverluste (PclPcl) leiten sich streng aus der makroskopischen Anwendung der Maxwell-Gleichungen auf einen homogenen, isotropen, kontinuierlichen Leiter planarer Geometrie endlicher Dicke ab. Betrachtet wird ein Stahlblech der Dicke dd in Richtung der zz-Achse (von z=d/2z = -d/2 bis z=+d/2z = +d/2), unendlicher Ausdehnung in der yy-Achse und Länge in der xx-Achse. Angenommen werden eine effektive magnetische Permeabilität μ\mu, eine elektrische Leitfähigkeit σ=1/ρ\sigma = 1/\rho und ein zur Oberfläche paralleler magnetischer Flussdichtevektor: 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}

Unter der Annahme, dass die Blechdicke dd wesentlich geringer ist als die elektromagnetische Skintiefe δ=2/(ωσμ)\delta = \sqrt{2/(\omega \sigma \mu)}, verteilt sich die magnetische Induktion räumlich homogen über den Blechquerschnitt: Bx(z,t)Bpsin(ωt)B_x(z,t) \approx B_p \sin(\omega t). Durch Integration über zz von der zentralen Neutralebene (z=0z = 0, wo Ey=0E_y = 0 aus Symmetriegründen gilt):

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)

Die induzierte Stromdichte J(z,t)=σE(z,t)\mathbf{J}(z,t) = \sigma \mathbf{E}(z,t) dissipiert eine Momentanleistung durch Joulesche Erwärmung mit der lokalen Raumleistungsdichte p(z,t)=Jy2(z,t)/σ=σEy2(z,t)p(z,t) = J_y^2(z,t) / \sigma = \sigma E_y^2(z,t). Durch Integration über die gesamte Blechdicke und zeitliche Mittelung über eine Periode 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]

Division durch die Massendichte ρm\rho_m liefert die normierte spezifische Formulierung:

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]

Diese Herleitung verdeutlicht die direkt quadratische Abhängigkeit von der Blechdicke dd und der Grundfrequenz ff. Dies begründet die Notwendigkeit der Reduzierung industrieller Nenndicken von 0,35 mm0{,}35 \ mm (M4) auf 0,23 mm0{,}23 \ mm oder 0,18 mm0{,}18 \ mm (laserbehandelte Ultra-Low-Loss-Güten).

Anomale Verluste bzw. Überschussverluste (Excess Losses)

In der Praxis erwies sich die direkte Summe aus quasistatischen Hystereseverlusten und klassischen Wirbelstromverlusten als signifikant niedriger als die im Epstein-Rahmen wattmetrisch gemessenen Gesamtverluste (Ptotal,gemessen>Ph+PclP_{total,gemessen} > P_h + Pcl). Diese energetische Diskrepanz wurde historisch als Verlustanomalie bezeichnet.

Das mikroskopische Modell von Pry und Bean (1958) wies nach, dass sich der magnetische Fluss nicht gleichmäßig über das gesamte Blechkontinuum ändert, sondern sich innerhalb der magnetischen Domänen konzentriert und durch die lokale Verschiebung von 180180^\circ-Bloch-Wänden übertragen wird. Die Verschiebungsgeschwindigkeit einer Domänenwand vw(t)v_w(t) unter sinusförmiger Erregung induziert in der unmittelbaren Wandumgebung einen lokal wesentlich steileren elektrischen Feldgradienten als von der klassischen Kontinuumstheorie vorhergesagt:

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

Hierbei steht LL für den Abstand zwischen benachbarten Domänenwänden und BsB_s für die Sättigungsinduktion. Der im mikroskopischen Volumen neben der Wand konzentrierte lokale Mikrowirbelstrom erzeugt ein magnetisches Gegenbremsfeld. Bertotti verallgemeinerte diese physikalischen Zusammenhänge durch die statistische Theorie der Aktiven Magnetischen Objekte (Magnetic Objects, MO), verstanden als kohärente Gruppen von Domänenwänden, die mit lokalen Inhomogenitäten des Mediums wechselwirken:

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]

Dabei ist G=0,1356G = 0{,}1356 eine dimensionslose Konstante der Wirbelfelddissipation in unendlichen Medien, SS die Querschnittsfläche des Blechs senkrecht zur Flussrichtung und V0V_0 ein intrinsischer Materialparameter mit der Einheit der magnetischen Feldstärke (A/mA/m), der die stochastische Wechselwirkung zwischen den magnetischen Objekten und dem mikroskopischen Verankerungsfeld (Internal Pinning Field Distribution) quantifiziert.

CRGO-Güteklasse (Klassifizierung AISI / EN) Nenndicke dd (mm) Max. Induktion BB bei 800 A/m800 \ A/m (T) Spezifische Verluste P1.7/50P_{1.7/50} (W/kg) Spezifische Verluste P1.7/60P_{1.7/60} (W/kg) Stapelfaktor (%)
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 Konventionell (27Q110) 0,27 1,90 - 1,93 0,95 - 1,05 1,24 - 1,38 95,0\ge 95{,}0
Hi-B Hohe Permeabilität (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
Amorphes Metall (Fe-Si-B Metglas 2605SA1) 0,025 1,56 - 1,58 0,18 - 0,25 0,24 - 0,33 86,0\ge 86{,}0

Nichtlineare Effekte, Oberschwingungsverzerrung und Gleichstromvormagnetisierung (DC-Bias)

In modernen Energiesystemen führt der massive Einsatz von Leistungselektronik (Großwechselrichter für Photovoltaikparks, HGÜ-Verbindungen, FACTS-Systeme und Frequenzumrichter) zu nicht-sinusförmigen Magnetisierungsverläufen und Spannungen mit erheblichem Oberschwingungsgehalt (THDvTHD_v). Die nichtlineare Charakteristik der dynamischen Magnetisierungskurve B(H)B(H) transformiert diese Anregungen in massive anomale Energieverluste.

Harmonische Erregung und untergeordnete Hystereseschleifen (Minor Loops)

Enthält die anliegende Spannungsform Oberschwingungen höherer Ordnung (h=3,5,7,h = 3, 5, 7, \dots), weist die Zeitableitung der magnetischen Flussdichte B/t\partial \mathbf{B}/\partial t innerhalb eines Grundschwingungszyklus mehrere Nulldurchgänge auf. Physikalisch erzwingt dies eine transiente Richtungsumkehr der Domänenwandbewegung vor dem Erreichen der makroskopischen Vollsättigung, wodurch untergeordnete Hystereseschleifen (Minor Hysteresis Loops) innerhalb der Hauptsättigungsschleife entstehen.

Pcl,harm=σπ2d26ρmh=1(hf1)2Bp,h2=Pcl,1h=1h2(Bp,hBp,1)2P_{cl,harm} = \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

Unter rein sinusförmiger Spannung gilt BhVh/(hf1)B_h \propto V_h / (h f_1). Weist die Spannung eine Oberschwingungsverzerrung mit den Amplituden VhV_h auf, nähert sich das harmonische Induktionsfeld über Bp,hVhhV1Bp,1B_{p,h} \approx \frac{V_h}{h \, V_1} B_{p,1} an. Einsetzen in die Formulierung der klassischen Wirbelströme ergibt:

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)

Die Überschussverluste folgen dieser linearen Entkopplung jedoch nicht. Die Überlagerung mehrerer Frequenzen beschleunigt die lokale Geschwindigkeit der Bloch-Wände gemäß einem Integral über die Flusszeitableitung:

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

Magnetisierung mit Gleichstromüberlagerung (DC-Bias)

Die Einspeisung von Gleichströmen in Transformatorwicklungen — verursacht durch geomagnetisch induzierte Ströme (GIC), direkt gekoppelte Wechselrichtertopologien oder asymmetrische Monopolfehler in HGÜ-Übertragungsleitungen mit Erdrückleitung — bewirkt eine Verschiebung des dynamischen Hysteresezyklus auf der Feldstärkeachse (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)

Diese asymmetrische Verschiebung treibt den Kern während einer vollen Halbwelle tief in das Sättigungsknie (halbzyklische Sättigung). Im Sättigungsbereich (B>1,9 TB > 1{,}9 \ T bei CRGO) bricht die differentielle Permeabilität μd=dB/dH\mu_d = dB/dH drastisch auf Werte nahe der Vakuumpermeabilität (μ0=4π×107 H/m\mu_0 = 4\pi \times 10^{-7} \ H/m) ein. Dies führt zu massiven asymmetrischen Magnetisierungsstromspitzen (imag(t)imag(t) mit Scheitelfaktoren von über 10). Das Magnetfeld wird von einem im Blechpaket geführten Vektor in einen hochenergetischen dreidimensionalen Streufluss transformiert, der benachbarte Strukturteile durchdringt (Presskonstruktionen, Kesselwandungen, Zugstangen und magnetische Schirmbleche).

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

Der Blindleistungsbedarf QQ steigt um Größenordnungen an, was zu schweren Klemmenspannungseinbrüchen, gerad- und ungeradzahligen Harmonischen niedriger Ordnung (insbesondere h=2h = 2 und h=3h = 3) sowie katastrophalen thermischen Überhitzungen durch induzierte Wirbelströme in metallischen Konstruktionsteilen führt.

Forensische Analyse thermischer Fehler und kerninduzierter dielektrischer Degradation

Elektromagnetische Fehlfunktionen und thermische Degradationsprozesse im Transformatorkern stellen einen kritischen Ausfallmechanismus bei Hochleistungstransformatoren dar. Im Gegensatz zu mechanischen Wicklungsfehlern infolge radialer und axialer Kurzschlusskräfte manifestieren sich Kernfehler typischerweise als schleichende elektrochemische und thermische Zersetzungsprozesse über lange Betriebszeiträume.

Degradation der Zwischenlagenisolation (Carlite-Beschichtung)

Die anorganische Oberflächenisolierung der Elektrobleche (Schichtdicke 1,53,0 μm\sim 1{,}5 - 3{,}0 \ \mu m, bestehend aus einer Forsterit-Grundschicht Mg2SiO4Mg_2SiO_4 mit einer Deckschicht aus Aluminiumphosphat AlPO4AlPO_4 und Phosphorsäure) ist für Zwischenlagenspannungen von wenigen Volt ausgelegt (1,5 V1{,}5 \ V bis 5,0 V5{,}0 \ V effektiv nach IEC 60404-6). Unter kontinuierlicher lokaler thermischer Überlastung (T>150CT > 150^\circ C), vibrationsbedingter Ermüdung durch Magnetostriktion (100/120 Hz100/120 \ Hz und akustische Harmonische) oder übermäßigen mechanischen Druckspannungen bei der Jochpressung fragmentiert und pulverisiert diese dielektrische Schicht mechanisch.

Der Verlust der galvanischen Trennung benachbarter Bleche erzeugt Mikroschweißungen durch Lichtbogenbildung oder direkten galvanischen Kontakt. Bilden sich zwei oder mehr Kontaktpunkte innerhalb desselben Blechpakets aus, entsteht eine großflächige Kurzschlusswindung, die den Hauptfluss des Kerns umschlingt. Die im geschlossenen Kreis induzierte elektromagnetische Kraft berechnet sich zu:

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}

Da der ohmsche Schleifenwiderstand RloopRloop im Milliohmbereich (mΩm\Omega) liegt, erreicht der Fehlerstrom IfaultIfault im stationären Zustand Amplituden von hunderten bis tausenden Ampere. Dies führt über den Jouleschen Effekt zu massiver lokaler Wärmedissipation und Hot-Spots mit Temperaturen, die schnell 500C800C500^\circ C - 800^\circ C übersteigen.

Degradation des Isolieröls und DGA-Forensik (Gelöste Gasanalyse)

Die extreme Wärmeübertragung vom Kern-Hot-Spot auf die umgebende Isolierflüssigkeit (naphthenisches oder paraffinisches Mineralöl sowie synthetische Ester) initiiert die Pyrolyse und thermokatalytische Spaltung der Kohlenwasserstoffketten (CnH2n+2C_n H_{2n+2}). Die kovalenten Kohlenstoff-Wasserstoff- (CHC-H, Bindungsenergie 413 kJ/mol\sim 413 \ kJ/mol) und Kohlenstoff-Kohlenstoff-Bindungen (CCC-C, 348 kJ/mol\sim 348 \ kJ/mol; C=CC=C, 614 kJ/mol\sim 614 \ kJ/mol; CCC \equiv C, 839 kJ/mol\sim 839 \ kJ/mol) brechen in Abhängigkeit vom lokalen thermischen Kontaktgradienten auf.

Gemäß den internationalen Normen IEC 60599 und IEEE C57.104 erlaubt das Profil der im Öl gelösten Gase (Dissolved Gas Analysis, DGA) die präzise thermodynamische Zuordnung des Fehlertemperaturbereichs im Transformatorkern:

Thermischer Fehlerbereich Dominantes Brenngas Diagnose-Gasverhältnisse (IEC 60599 / Duval) Physikalischer Fehlermechanismus im Kern Erforderliche betriebliche Maßnahme
Thermischer Fehler T<300CT < 300^\circ C (T1) Methan (CH4CH₄) / Ethan (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 Diffuse Jochüberhitzung, beginnende Zersetzung von Zwischenlagenlacken oder Kernklebern. Monatliche chromatographische DGA-Überwachung; Überprüfung von Netzoberschwingungen.
Thermischer Fehler 300C<T<700C300^\circ C < T < 700^\circ C (T2) Ethylen (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 Kurzschlusswindungen über mehrere CRGO-Bleche. Schwere Zersetzung der anorganischen Isolation. Lastreduzierung auf 70%70\%; Messung der Leerlaufverluste und Kernisolationswiderstandsprüfung.
Thermischer Fehler T>700CT > 700^\circ C (T3) Ethylen (C2H4C₂H₄) mit Acetylenspuren (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 Lokales Kernschmelzen (Core Burning); Massekurzschluss zwischen Blechpaket und Pressbalken. Sofortige Notabschaltung. Endoskopische Kesselinspektion und Kernrekonstruktion.
Teilentladungen / Funkenbildung im Kern Wasserstoff (H2H₂) dominant mit CH4CH₄ CH4H2<0,1\frac{CH₄}{H₂} < 0{,}1, H2H₂ überschreitet Basisgrenzwerte (> 100 ppm) Potenzialfreier Kern (Verlust der Einpunkterdung), kapazitive Entladungen gegen die Kesselwand. Überprüfung des externen Erdungskreises für Kern und Schirmung mittels 2,5kV2{,}5 - kV-Isolationsmessung.

Fehler durch Mehrfacherderdung (Erdschleifen-Zirkulationsströme)

Normativ vorgeschrieben darf der Transformatorkern ausschließlich an einem einzigen definierten Punkt galvanisch geerdet werden, realisiert über ein Kupferband, das an eine isolierte Erddurchführung auf dem Kesseldeckel geführt ist. Diese Einpunkterdung verhindert den Aufbau statischer Schwebepotenziale durch kapazitive Einkopplung von den Oberspannungswicklungen (Vind=Cwc/(Cwc+Ccg)VwindingVind = C_{w-c}/(C_{w-c} + C_{c-g}) \cdot Vwinding).

Kommt es durch Isolationsversagen an Zugankerbolzen (Tie-Rods), metallische Schlammablagerungen am Kesselboden oder Montagefehler zu einer unbeabsichtigten zweiten Erdung am entgegengesetzten Kernende, schließt sich eine großflächige Leiterschleife, die den magnetischen Hauptfluss umfasst:

Vloop=ddtFla¨che zwischen ErdpunktenBdSVloop = -\frac{d}{dt} \iint_{Fl\ddot{a}che\ zwischen\ Erdpunkten} \mathbf{B} \cdot d\mathbf{S}

Diese Schleifenspannung treibt einen kontinuierlichen Zirkulationsstrom durch die Kesselwand und die Erdungsschiene des Kerns, der typischerweise Amplituden zwischen 10 A10 \ A und über 150 A150 \ A erreicht. Dies führt zum thermischen Verbrennen der Isolierhülsen der Kernbolzen, zur Verkokung des umgebenden Öls und zu alarmierenden Bildungsraten von Ethylen (C2H4C₂H₄) und Methan (CH4CH₄), ohne dass der Überstrom- oder Differentialschutz (87T87T) eine Anomalie in den Netzströmen erkennt.

Fortschrittliche Methoden zur Verlustminderung in Design und Fertigung

Die Optimierung der elektromagnetischen Güte von Transformatorkernen erfordert einen multidisziplinären Ansatz, der werkstoffliche Kristallgitteroptimierungen, laserinduzierte Domänenmikrostrukturierung sowie die dreidimensionale topologische Auslegung der Jochschichtungen umfasst.

Step-Lap-Jochverbindungen unter 4545^\circ

Die geometrische Schichtung von Schenkeln und Jochen in Blechpaketen erzeugt unvermeidliche Diskontinuitäten im Magnetkreis. Bei traditionellen stumpfen 9090^\circ-Überlappungsverbindungen (Butt-Lap Joints) muss der magnetische Fluss senkrecht zur Walzrichtung (in der harten kristallographischen 110\langle 110 \rangle-Richtung) zwischen benachbarten Blechen übertreten. Dies erzwingt eine starke Brechung der Flusslinien, erhöht den Grenzflächenreluktanzsprung drastisch und erzeugt lokale Sättigungsspitzen an den Schnittkanten.

Moderne automatisierte Schneid- und Schichtverfahren setzen auf Gehrungsschnitte unter 4545^\circ mit kontinuierlich versetzten Blechlagen, bezeichnet als Step-Lap-Schichtung (typischerweise 5 bis 7 Stufen pro Versatzzyklus). Diese Abstufung verteilt den fertigungstechnischen Trennspalt axial über eine vorgegebene Schichtlänge:

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

Durch den stufenweisen Lagenversatz um Δx3,57,0 mm\Delta x \approx 3{,}5 - 7{,}0 \ mm springt der Fluss nicht mehr orthogonal über einen einzigen hochreluktanten Stoßspalt in das Nachbarblech, sondern verteilt sich progressiv über eine wesentlich größere effektive geometrische Koppelfläche:

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

Die Step-Lap-Technologie senkt die Gesamtverbindungsverluste an den Jochen um bis zu 25%30%25\% - 30\%, reduziert den Magnetisierungsblindleistungsbedarf um 40%40\% und vermindert die Schallemission des Transformators um 3bis6 dB(A)3 bis 6 \ dB(A) gegenüber konventionellen Schichtungen.

Laserinduzierte Domänenverfeinerung (Laser Scribing)

In hochpermeablen CRGO-Güten (Hi-B) weisen die Kristallkörner eine beträchtliche Größe auf (1030 mm\sim 10 - 30 \ mm). Dies begünstigt zwar eine nahezu perfekte Goss-Ausrichtung (3\le 3^\circ), wirkt sich jedoch nachteilig auf die dynamischen Wirbelstromverluste aus, da der Abstand zwischen den 180180^\circ-Bloch-Wänden (2d02d_0) proportional mit der Quadratwurzel der Korngröße anwächst:

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

Ein großer Domänenabstand erhöht die lokale Wandverschiebungsgeschwindigkeit vw(t)v_w(t) stark, was die anomalen Überschussverluste (PexcPexc) drastisch ansteigen lässt.

Beim Verfahren der Laser-Domänenverfeinerung (Laser Scribing) tastet ein hochfokussierter Dauerstrichlaserstrahl (Nd:YAG- oder Hochleistungsfaserlaser) das Blech quer zur Walzrichtung in regelmäßigen Abständen (p48 mmp \approx 4 - 8 \ mm) ab. Der thermische Energieeintrag bewirkt eine extrem schnelle lokale Erwärmung ohne Materialabtrag, wodurch infolge des steilen Abkühlungsgradienten ein permanentes elastisches Druckeigenspannungsfeld unter der Oberfläche erzeugt wird:

σ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} \iint E \, \alpha_{th} \, \Delta T \, dy \, dz

Diese periodischen internen Spannungsstreifen fungieren als künstliche magnetoelastische Energiebarrieren, die breite 180180^\circ-Domänen aufspalten und den mittleren Domänenwandabstand 2d02d_0 auf die Hälfte oder ein Drittel reduzieren. Daraus resultiert direkt:

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

Hierbei ist zwingend zwischen wärmeunbeständigen Laserbehandlungen (die den Verfeinerungseffekt bei Entspannungsglühungen über 500C500^\circ C, wie bei Schnittbandkernen erforderlich, verlieren) und wärmebeständigen Verfahren (z. B. plasmastrahlinduzierte plastische Mikroverformung oder Ätzrillen) zu unterscheiden, die ihre Mikrokernverankerung auch nach Glühungen bei 800C800^\circ C vollständig beibehalten.

Implementierung und Rechenmodellierung in der Vexten Suite

Die rigorose Quantifizierung von Kernverlusten und deren multiphysikalische Rückkopplung mit dem Netzbetrieb wird in der Vexten Suite über integrierte Module gelöst: nichtlineare Oberschwingungs-Leistungsflussberechnung (nach IEEE 519 / IEEE C57.110), dynamische thermische Netzkalkulation (IEC 60076-7) sowie transiente elektromagnetische Kurzschluss- und Sättigungsanalyse (IEC 60909 / IEEE 141).

Korrekturalgorithmus für Oberschwingungsverzerrung und K-Faktor

Bei signifikanten harmonischen Oberschwingungsspektren aus nichtlinearen Lasten oder regenerativen Erzeugern führt das Modul Vexten PowerFlow & Harmonics eine dynamische Verlusttrennung durch, indem der Kernverlustfaktor anhand des normierten Spektralkoeffiziententensors neu berechnet wird:

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}

Die Software verknüpft diesen Faktor mit der thermischen Leistungsreduktion (Derating) des Transformators über die Bestimmung des K-Faktors sowie des Oberschwingungsverlustfaktors für Wirbelströme (FHLFHL). Der maximal zulässige Betriebsstrom wird so eingeregelt, dass weder der Wicklungs-Hot-Spot noch die Kernkerntemperatur die Grenzwerte der Isolationsklasse überschreiten:

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

Ablauf der multiphysikalischen Modellierung in der Vexten Suite

  1. Definition des magnetischen Erregungsspektrums: Import von Feldmessdaten oder normierten IEEE-Spannungsharmonischen (VhV_h) sowie des quasistationären Gleichstromanteils (IdcIdc) aus dem netzgekoppelten Leistungsfluss im Modul Netzqualität.
  2. Zuweisung des Magnetwerkstoffs in der Werkstoffdatenbank: Auswahl der spezifischen CRGO-Blechgüte in den Maschineneigenschaften (z. B. M3, 23Q090 oder 20Q080-L mit Laser-Scribing). Die Vexten Suite lädt automatisch die kalibrierten mikrostrukturellen Kennwerte: Leitfähigkeit σ\sigma, Hysteresekoeffizient khk_h, Nenndicke dd, Dichte ρm\rho_m und den mikroskopischen Bertotti-Parameter V0V_0.
  3. Konfiguration der Jochschichtungstopologie: Parametrierung der Schichtgeometrie: konventionelle 9090^\circ-Stumpfschichtung, 5-stufiges oder 7-stufiges Step-Lap-Verfahren. Das System beaufschlagt die Schnittzonen automatisch mit dem entsprechenden Reluktanzüberhöhungsfaktor (Kjoint\mathcal{K}_{joint}).
  4. Durchführung der transienten gekoppelten thermischen Analyse: Der Berechnungskern koppelt die elektromagnetischen Gesamtverluste Ptotal(x,y,z)Ptotal(x,y,z) als räumliche Wärmequellendichte in den thermischen Netzwerklöser gemäß IEC 60076-7 ein. Hieraus wird das dreidimensionale Temperaturprofil über die oberen und unteren Joche sowie die Kernschenkel ermittelt.
  5. Grenzwertvalidierung und präventive forensische Diagnose: Automatischer Abgleich der berechneten Temperaturgradienten mit den Gefahrenschwellenwerten nach IEEE C57.104. Überschreitet die Dauertemperatur im Blechpaket 130C130^\circ C, generiert das System Warnmeldungen zur Leistungsreduktion (Derating) und kalkuliert die theoretischen Gasbildungsraten (C2H4C₂H₄, CH4CH₄, H2H₂), was die Implementierung hochpräziser zustandsorientierter Instandhaltungsstrategien (CBM) ermöglicht.

Auf Basis dieser integrativen Berechnungsplattform können Entwicklungs-, Inbetriebnahme- und Betriebsingenieure die thermische Alterung des Transformatorkerns exakt prognostizieren, das Risiko katastrophaler Schäden durch Zirkulationsströme eliminieren und die Auswahl moderner Magnetwerkstoffe für Transformatoren mit höchsten Wirkungsgradanforderungen fundiert optimieren.