Pérdidas Anómalas y Foucault en Acero CRGO

Termodinámica Electromagnética y Microestructura del Acero al Silicio CRGO

El núcleo magnético de un transformador de potencia opera como el circuito de acoplamiento de flujo fundamental entre devanados, operando bajo regímenes cíclicos no lineales de excitación electromagnética. La eficiencia global y la densidad de potencia del transformador están condicionadas de forma determinante por el comportamiento microestructural y cristalográfico del material ferromagnético utilizado. El estándar industrial predominante para transformadores de media, alta y extra alta tensión (EHV/UHV) es el acero eléctrico de grano orientado laminado en frío (Cold-Rolled Grain-Oriented Silicon Steel, CRGO), optimizado mediante tratamientos termomecánicos para alinear sus ejes cristalográficos de fácil magnetización en la dirección de laminación.

El acero CRGO presenta una estructura cristalina cúbica centrada en el cuerpo (BCC, red $\alpha$ -Fe), donde el eje cristalográfico $\langle 100 \rangle$ constituye la dirección de mínima energía de anisotropía magnetocristalina, mientras que las direcciones $\langle 110 \rangle$ y $\langle 111 \rangle$ representan ejes de anisotropía intermedia y dura, respectivamente. Mediante una secuencia rigurosa de laminación en frío y recocidos de recristalización secundaria (inhibida mediante precipitados de sulfuro de manganeso, $MnS$ , o nitruro de aluminio, $AlN$ ), se desarrolla la denominada textura de Goss, caracterizada por la orientación cristalográfica $\{110\}\langle 001 \rangle$ . En esta configuración, los planos $\{110\}$ se sitúan paralelos a la superficie de la lámina y las direcciones $\langle 001 \rangle$ se orientan paralelamente a la dirección longitudinal de laminación (RD, Rolling Direction).

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)

Donde $E_a$ es la densidad de energía de anisotropía magnetocristalina ( $J/m ^3$ ), $K_1 \approx 4.8 \times 10^4 J/m ^3$ y $K_2 \approx 1.5 \times 10^4 J/m ^3$ son las constantes de anisotropía del hierro a temperatura ambiente ( $298 K$ ), y $\alpha_1, \alpha_2, \alpha_3$ representan los cosenos directores del vector de magnetización $\mathbf{M}$ respecto a los ejes de la celda unitaria. La desalineación angular media de la textura de Goss respecto a la dirección de laminación define la calidad del acero: en aceros CRGO convencionales la dispersión angular oscila entre $5^\circ$ y $7^\circ$ , mientras que en aceros de alta permeabilidad (Hi-B) se reduce a rangos de $2^\circ$ a $3^\circ$ .

La adición de silicio ( $Si$ ) en concentraciones nominales del $3.0\%$ al $3.4\%$ en peso incrementa la resistividad volumétrica eléctrica ( $\rho$ ) desde aproximadamente $10 \times 10^{-8} \ \Omega\cdot m$ (para hierro puro) hasta $(45 - 50) \times 10^{-8} \ \Omega\cdot m$ . Este incremento de resistividad mitiga drásticamente el transporte libre de cargas parásitas. Sin embargo, concentraciones superiores al $3.5\%$ en peso inducen fragilidad mecánica severa debido a la formación de fases ordenadas intermetálicas $Fe _3 Si$ y $FeSi$ ( $B2$ y $D0_3$ ), impidiendo los procesos de conformado en frío a escala industrial. Adicionalmente, el silicio reduce la constante de magnetoestricción de saturación longitudinal ( $\lambda_{100}$ ), minimizando las deformaciones elásticas inducidas por el campo magnético y la consecuente emisión de ruido acústico.

\lambda(H) = \frac{\Delta l}{l} = \frac{3}{2} \lambda_s \left( \left(\frac{M}{M_s}\right)^2 - \frac{1}{3} \right)

La estructura de dominios magnéticos (dominios de Weiss) en el CRGO se compone de dominios principales orientados a $180^\circ$ separados por paredes de Bloch, y dominios de cierre a $90^\circ$ en las inmediaciones de discontinuidades superficiales e inclusiones. El ancho de dominio característico ( $2d_0$ ) depende del balance entre la energía magnetostática superficial y la energía elástica de la pared del dominio. La aplicación de revestimientos dieléctricos en tensión superficial (como el recubrimiento de fosfatos y silicatos, denominado comercialmente Carlite) introduce una tracción biaxial permanente en la lámina ( $\sigma_t \approx 4 - 10 MPa$ ), reduciendo el ancho medio de los dominios a $180^\circ$ y eliminando dominios secundarios transversales de cierre, disminuyendo radicalmente las pérdidas por corrientes parásitas anómalas.

Tripartición Clásica y Espectro de Pérdidas: El Modelo de Bertotti

Bajo excitación magnetodinámica cíclica, la disipación volumétrica de energía por unidad de tiempo en el núcleo se cuantifica mediante la integral cerrada sobre el ciclo del lazo de histéresis dinámico. La formulación general de pérdidas específicas de potencia ( $P_t$ , en $W/kg$ ) se rige por la teoría estadística de pérdidas electromagnéticas formulada por Giorgio Bertotti, la cual supera las limitaciones de las formulaciones empíricas clásicas de Steinmetz al descomponer el fenómeno en tres mecanismos disipativos desacoplados físicamente:

P_{ total }(f, B_p) = P_h(f, B_p) + Pcl(f, B_p) + Pexc(f, B_p)

P_{ 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}

Donde $f$ representa la frecuencia fundamental de excitación ( $Hz$ ), $B_p$ es el valor pico de la inducción magnética ( $T$ ), $\sigma$ es la conductividad eléctrica del material ( $S/m$ ), $d$ es el espesor de la chapa laminada ( $m$ ), $\rho_m$ es la densidad de masa volumétrica del acero ( $kg/m ^3$ , típicamente $7650 kg/m ^3$ ), $k_h$ es el coeficiente de histéresis cuasi-estática, $\alpha$ es el exponente de Steinmetz (frecuentemente $1.6 \le \alpha \le 2.0$ ), y $kexc$ es el parámetro microscópico de pérdidas de exceso.

Pérdidas por Histéresis Cuasi-Estática

Las pérdidas por histéresis estática ( $P_h$ ) se originan en los procesos microscópicos de disipación termodinámica irreversible cuando las paredes de los dominios magnéticos se desplazan a través de un medio ferromagnético con defectos de red cristalina, tales como dislocaciones, vacancias, bordes de grano, tensiones residuales e inclusiones no metálicas (e.g., precipitados de $SiO _2$ o $Fe _3 C$ ). Durante el movimiento de una pared de Bloch de $180^\circ$ , ésta queda anclada mecánicamente en estos centros de fijación (pinning sites).

Para liberar la pared de dominio del pozo de potencial de fijación es necesario un incremento del campo magnético exterior hasta alcanzar el campo crítico local. Una vez superada esta barrera energética, la pared avanza discontinuamente a gran velocidad hacia el siguiente pozo de potencial mediante un salto microscópico no lineal, fenómeno conocido como efecto Barkhausen. La energía disipada en este salto se convierte irreversiblemente en fonones (calor) en la red cristalina:

w_h = \oint_{ ciclo } \mathbf{H} \cdot d\mathbf{B} = \mu_0 \oint_{ ciclo } H \, dM \quad \left[ J/m ^3\right]

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]

El área del lazo de histéresis cuasi-estático ( $w_h$ ) depende exclusivamente de la estructura metalúrgica del material y es independiente de la frecuencia de variación temporal del flujo para regímenes donde el tiempo de relajación del dominio es infinitamente menor que el periodo de la excitación. La modelación matemática de $P_h$ en presencia de formas de onda complejas o saturación no lineal requiere la formulación del operador continuo de Preisach:

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

Donde $\hat{\gamma}_{\alpha\beta}$ es un operador histerético elemental biestable con umbrales de conmutación $\alpha$ (de encendido) y $\beta$ (de apagado), y $\mu(\alpha, \beta)$ es la función de distribución de densidad de Preisach específica del material CRGO ensayado, identificable experimentalmente mediante familias de curvas de inversión de primer orden (FORC, First-Order Reversal Curves).

Pérdidas por Corrientes Parásitas Clásicas (Foucault)

Las pérdidas clásicas por corrientes de Foucault ( $Pcl$ ) derivan rigurosamente de la aplicación macroscópica de las ecuaciones de Maxwell a un conductor homogéneo e isotrópico continuo de geometría planar finita. Considérese una lámina de acero de espesor $d$ en la dirección del eje $z$ (desde $z = -d/2$ hasta $z = +d/2$ ), ancho infinito en el eje $y$ , y longitud en el eje $x$ . Asumiendo una permeabilidad magnética efectiva $\mu$ , una conductividad volumétrica $\sigma = 1/\rho$ , y un vector de densidad de flujo magnético paralelo a la superficie: $\mathbf{B}(t) = B_p \sin(\omega t) \hat{\mathbf{x}}$ .

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

Asumiendo que el espesor $d$ es sustancialmente menor que la profundidad de penetración electromagnética (skin depth) $\delta = \sqrt{2/(\omega \sigma \mu)}$ , la inducción magnética se distribuye de manera espacialmente uniforme a través de la sección transversal de la lámina: $B_x(z,t) \approx B_p \sin(\omega t)$ . Integrando respecto a $z$ desde el plano neutro central ( $z = 0$ , donde $E_y = 0$ por simetría):

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

La densidad de corriente inducida $\mathbf{J}(z,t) = \sigma \mathbf{E}(z,t)$ disipa una potencia instantánea por efecto Joule cuya densidad volumétrica local es $p(z,t) = J_y^2(z,t) / \sigma = \sigma E_y^2(z,t)$ . Integrando sobre el espesor total de la chapa y promediando en el dominio temporal durante un periodo $T = 2\pi/\omega$ :

\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

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

Dividiendo por la densidad de masa volumétrica $\rho_m$ se obtiene la formulación específica normalizada:

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]

Esta deducción expone la dependencia cuadrática directa respecto al espesor $d$ de la lámina y a la frecuencia fundamental $f$ , justificando la necesidad de reducir los calibres comerciales desde $0.35 mm$ (M4) hasta $0.23 mm$ o $0.18 mm$ (grados láser de ultra-bajas pérdidas).

Pérdidas Anómalas o de Exceso (Excess Losses)

Históricamente, la suma directa de las pérdidas por histéresis cuasi-estática y las pérdidas clásicas de Foucault resultaba significativamente inferior a las pérdidas totales medidas experimentalmente mediante el método del vatímetro en cuadro de Epstein ( $P_{ total, med } > P_h + Pcl$ ). Esta discrepancia energética se denominó formalmente anomalía de pérdidas.

El modelo microscópico de Pry y Bean (1958) demostró que el flujo magnético no varía uniformemente en todo el continuo de la lámina, sino que se concentra exclusivamente dentro de los dominios magnéticos y se transfiere a través del desplazamiento localizado de las paredes de Bloch de $180^\circ$ . La velocidad de desplazamiento de una pared de dominio $v_w(t)$ bajo excitación sinusoidal induce un gradiente de campo eléctrico local infinitamente más pronunciado en la vecindad inmediata de la pared que el previsto por la teoría clásica del continuo:

v_w(t) = \frac{L}{2 B_s} \frac{dB(t)}{dt}

Donde $L$ representa el espaciamiento inter-paredes de dominio y $B_s$ es la inducción de saturación. La micro-corriente parásita confinada en el micro-volumen adyacente a la pared genera un campo de frenado magnético local. Bertotti generalizó esta física formulando la teoría estadística de los Objetos Magnéticos Activos (Magnetic Objects, MO), entendidos como agregados coherentes de paredes de dominio que interactúan con las heterogeneidades locales del medio:

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]

Donde $G = 0.1356$ es una constante adimensional asociada a la disipación del campo de Foucault en medios infinitos, $S$ es el área de la sección transversal de la lámina perpendicular a la dirección del flujo magnético, y $V_0$ es un parámetro intrínseco del material con dimensiones de campo magnético ( $A/m$ ), el cual cuantifica la interacción estocástica entre los objetos magnéticos y el campo microscópico de anclaje (internal pinning field distribution).

Grado CRGO (Clasificación AISI / EN)	Espesor Nominal $d$ (mm)	Inducción Máxima $B$ a $800 A/m$ (T)	Pérdidas Específicas $P_{1.7/50}$ (W/kg)	Pérdidas Específicas $P_{1.7/60}$ (W/kg)	Factor de Apilamiento (\%)
M4 (AISI 35G155 / EN 10107)	0.35	1.82 - 1.84	1.25 - 1.35	1.64 - 1.77	$\ge 96.0$
M3 (AISI 30G130 / EN 10107)	0.30	1.84 - 1.86	1.05 - 1.15	1.38 - 1.51	$\ge 95.5$
Hi-B Convencional (27Q110)	0.27	1.90 - 1.93	0.95 - 1.05	1.24 - 1.38	$\ge 95.0$
Hi-B Alta Permeabilidad (23Q090)	0.23	1.92 - 1.95	0.80 - 0.90	1.05 - 1.18	$\ge 94.5$
Hi-B Laser-Scribed (20Q080-L)	0.20	1.93 - 1.96	0.70 - 0.78	0.92 - 1.02	$\ge 94.0$
Metal Amorfo (Fe-Si-B Metglas 2605SA1)	0.025	1.56 - 1.58	0.18 - 0.25	0.24 - 0.33	$\ge 86.0$

Efectos No Lineales, Distorsión Armónica y Componente Continua (DC Bias)

En sistemas de potencia contemporáneos, la propagación masiva de electrónica de potencia (inversores de gran escala para parques fotovoltaicos, enlaces HVDC, sistemas FACTS y accionamientos de frecuencia variable) introduce regímenes de magnetización no sinusoidales y tensiones con alto contenido armónico ( $THD _v$ ). La respuesta no lineal de la curva de imanación dinámica $B(H)$ transforma estas excitaciones en pérdidas energéticas anómalas severas.

Excitación Armónica y Bucles de Histéresis Menores

Cuando la forma de onda de tensión aplicada contiene armónicos superiores ( $h = 3, 5, 7, \dots$ ), la derivada temporal de la inducción magnética $\partial \mathbf{B}/\partial t$ presenta múltiples cruces por cero en el transcurso de un ciclo fundamental. Físicamente, esto fuerza la inversión transitoria del movimiento de las paredes de dominio magnético antes de completar la polarización macroscópica total, originando bucles de histéresis menores (minor hysteresis loops) anidados dentro del lazo principal de saturación.

P_{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

Bajo tensión sinusoidal pura, $B_h \propto V_h / (h f_1)$ . Si la tensión contiene una distorsión armónica con componentes $V_h$ , el campo inducido armónico se aproxima por $B_{p,h} \approx \frac{V_h}{h \, V_1} B_{p,1}$ . Sustituyendo en la formulación de corrientes clásicas:

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)

Sin embargo, las pérdidas por exceso no siguen este desacoplamiento lineal. La interacción de múltiples frecuencias acelera la velocidad local de las paredes de Bloch de acuerdo con una integral temporal de derivada de flujo:

P_{exc, non-sin } = \frac{kexc}{\rho_m} \frac{1}{T} \int_0^T \left| \frac{d B(t)}{dt} \right|^{1.5} dt

Magnetización con Desplazamiento de Corriente Continua (DC Bias)

La inyección de corriente continua en las fases del transformador —producida por corrientes geomagnéticamente inducidas (GIC), topologías de inversores de acoplamiento directo o fallas asimétricas monopolares de líneas de transmisión HVDC con retorno por tierra— produce una traslación del ciclo dinámico sobre el eje del campo magnetizante ( $Hdc = N Idc / lfe$ ).

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

Este desplazamiento asimétrico fuerza el núcleo hacia el codo de saturación durante un semiciclo completo (saturación semi-cíclica). En la región de saturación ( $B > 1.9 T$ en CRGO), la permeabilidad diferencial $\mu_d = dB/dH$ colapsa drásticamente hacia valores cercanos a la permeabilidad del vacío ( $\mu_0 = 4\pi \times 10^{-7} H/m$ ). Esto induce picos masivos de corriente magnetizante asimétrica ( $imag(t)$ con factores de cresta superiores a 10), transformando el campo magnético de un vector confinado dentro del núcleo laminado en un flujo de dispersión tridimensional de gran magnitud que concatena elementos estructurales adyacentes (placas de apriete, cubas de transformadores, pernos de amarre y blindajes magnéticos).

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

El consumo de potencia reactiva $Q$ se eleva en órdenes de magnitud, originando severas caídas de tensión en bornes, armónicos pares e impares de bajo orden (especialmente $h = 2$ y $h = 3$ ), y calentamientos térmicos catastróficos por corrientes inducidas en piezas metálicas estructurales.

Análisis Forense de Fallas Térmicas y Degradación Dieléctrica Inducida por el Núcleo

Las patologías electromagnéticas y de degradación térmica originadas en el núcleo constituyen un mecanismo de fallo crítico en transformadores de potencia de gran capacidad. A diferencia de las fallas mecánicas en bobinados inducidas por esfuerzos de cortocircuito radiales y axiales, las fallas de núcleo suelen manifestarse como procesos de degradación electroquímica y térmica progresiva a lo largo de extensos periodos de servicio.

Degradación del Aislamiento Interlaminar (Revestimiento Tipo Carlite)

El aislamiento interlaminar superficial de las chapas magnéticas (espesor $\sim 1.5 - 3.0 \ \mu m$ , compuesto de una base de forsterita $Mg _2 SiO _4$ con sobrecapa de fosfato de aluminio $AlPO _4$ y ácido fosfórico) está diseñado para soportar tensiones interlaminares de unos pocos voltios ( $1.5 V$ a $5.0 V$ eficaces según norma IEC 60404-6). Sin embargo, bajo esfuerzos térmicos locales continuos ( $T > 150^\circ C$ ), fatiga vibracional por magnetoestricción ( $100/120 Hz$ y armónicos acústicos) o tensiones mecánicas excesivas por sobreapriete durante el ensamblaje de los yugos, el revestimiento dieléctrico se fragmenta y pulveriza mecánicamente.

La pérdida del aislamiento galvánico entre chapas contiguas genera micro-soldaduras por arco o contacto galvánico directo. Si se establecen dos o más puntos de contacto en el mismo paquete de laminación, se forma una espira en cortocircuito de gran sección que concatena el flujo magnético principal del núcleo. La fuerza electromotriz inducida en el bucle cerrado:

\mathcal{E}_{loop} = -\oint_{\partial S} \frac{\partial \mathbf{B}}{\partial t} \cdot d\mathbf{S} = -\omega \, B_p \, Aloop \cos(\omega t)

Ifault = \frac{\mathcal{E}_{loop}}{\sqrt{Rloop^2 + (\omega Lloop)^2}} \approx \frac{\mathcal{E}_{loop}}{Rloop}

Dado que la resistencia del bucle óhmico $Rloop$ es del orden de los miliohmios ( $m \Omega$ ), la corriente de falla $Ifault$ alcanza amplitudes de cientos a miles de amperios en régimen estacionario, disipando calor localizado mediante efecto Joule que genera puntos calientes (hot-spots) con temperaturas que superan rápidamente los $500^\circ C - 800^\circ C$ .

Degradación del Aceite Dieléctrico y Análisis Forense por DGA

La transferencia de calor extremo desde el punto caliente del núcleo hacia el fluido dieléctrico circundante (aceite mineral nafténico, parafínico o ésteres sintéticos) induce la pirólisis y escisión termocatalítica de las cadenas moleculares de hidrocarburos ( $C _n H _{2n+2}$ ). Los enlaces covalentes carbono-hidrógeno ( $C-H$ , energía de enlace $\sim 413 kJ/mol$ ) y carbono-carbono ( $C-C$ , $\sim 348 kJ/mol$ ; $C=C$ , $\sim 614 kJ/mol$ ; $C \equiv C$ , $\sim 839 kJ/mol$ ) se rompen en función del gradiente térmico de contacto.

De acuerdo con los estándares internacionales IEC 60599 e IEEE C57.104, el perfil de gases disueltos en el aceite (Dissolved Gas Analysis, DGA) permite la identificación termodinámica precisa del rango de temperatura de la falla en el núcleo:

Rango Térmico de Falla	Gas Combustible Dominante	Relaciones de Diagnóstico (IEC 60599 / Duval)	Mecanismo Físico en el Núcleo	Acción Crítica Requerida
Falla Térmica $T < 300^\circ C$ (T1)	Metano ( $CH _4$ ) / Etano ( $C _2 H _6$ )	$\frac{ C _2 H _2}{ C _2 H _4} < 0.1$ , $0.1 < \frac{ CH _4}{ H _2} < 1$ , $\frac{ C _2 H _4}{ C _2 H _6} < 1$	Calentamiento difuso de yugos, degradación inicial de barnices interlaminares o pintura de bloqueo.	Monitoreo cromatográfico periódico mensual; verificación de armónicos en la red.
Falla Térmica $300^\circ C < T < 700^\circ C$ (T2)	Etileno ( $C _2 H _4$ )	$\frac{ C _2 H _2}{ C _2 H _4} < 0.1$ , $\frac{ CH _4}{ H _2} > 1$ , $1 < \frac{ C _2 H _4}{ C _2 H _6} < 4$	Espiras en cortocircuito entre múltiples láminas de acero CRGO. Degradación severa del recubrimiento inorgánico.	Reducción de carga al $70\%$ ; ensayo de pérdidas en vacío y resistencia de aislamiento del núcleo.
Falla Térmica $T > 700^\circ C$ (T3)	Etileno ( $C _2 H _4$ ) con trazas de Acetileno ( $C _2 H _2$ )	$\frac{ C _2 H _2}{ C _2 H _4} < 0.15$ , $\frac{ CH _4}{ H _2} > 1$ , $\frac{ C _2 H _4}{ C _2 H _6} > 4$	Fusión localizada del paquete magnético (core burning); cortocircuito entre láminas y estructura de apriete.	Desconexión inmediata forzada. Inspección interna endoscópica y reconstrucción de núcleo.
Descargas Parciales / Chispas en Núcleo	Hidrógeno ( $H _2$ ) dominante con $CH _4$	$\frac{ CH _4}{ H _2} < 0.1$ , $H _2$ excede límites basales (> 100 ppm)	Flotación galvánica del núcleo (pérdida de tierra física única), produciendo descargas capacitivas a la cuba.	Verificación del circuito de puesta a tierra externa del núcleo y neutro mediante megóhmetro a $2.5 kV$ .

Falla por Múltiples Puestas a Tierra (Ground Loop Corrientes Circulantes)

Por diseño normativo, el núcleo de un transformador debe estar puesto a tierra galvánicamente en un único punto físico mediante una pletina de cobre conectada a un aislador pasatapas dedicado en la tapa del tanque. Esta puesta a tierra unificada previene la elevación del potencial electrostático flotante inducido por acoplamiento capacitivo con los devanados de alta tensión ( $V_{ ind } = C_{w-c}/(C_{w-c} + C_{c-g}) \cdot Vwinding$ ).

Si por degradación de aislamientos de pernos de sujeción (tie-rods), puenteo por sedimentos metálicos en el fondo del tanque o errores de montaje se genera una segunda conexión a tierra inadvertida en el extremo opuesto del núcleo, se forma un circuito cerrado de gran área física que concatena el flujo magnético principal:

V_{ loop } = -\frac{d}{dt} \iint_{\text{ área entre tierras }} \mathbf{B} \cdot d\mathbf{S}

Esta tensión induce una corriente de circulación continua por la estructura del tanque y las chapas de tierra del núcleo que típicamente varía entre $10 A$ y más de $150 A$ , quemando el aislamiento de los pernos pasantes, carbonizando el aceite adyacente y generando concentraciones alarmantes de etileno ( $C _2 H _4$ ) y metano ( $CH _4$ ) sin que las protecciones eléctricas de sobrecorriente o diferencial ( $87 T$ ) detecten anomalía alguna en las corrientes de línea.

Técnicas Avanzadas de Mitigación en el Diseño y Fabricación

La optimización del rendimiento electromagnético de los núcleos magnéticos requiere una aproximación multidisciplinaria que abarca la modificación cristalográfica del material, la microingeniería de dominios por radiación coherente y la optimización topológica tridimensional de las uniones en los yugos.

Uniones de Yugo Tipo Step-Lap a $45^\circ$

El ensamblaje geométrico de las columnas y los yugos en núcleos apilados introduce discontinuidades en el circuito magnético. En las uniones tradicionales a tope con solape simple a $90^\circ$ (Butt-Lap joints), el flujo magnético se ve forzado a cruzar entre láminas perpendicularmente a la dirección de laminación (en la dirección cristalográfica dura $\langle 110 \rangle$ ), provocando una refracción severa de líneas de flujo, un incremento drástico de la reluctancia magnética interfacial y saturaciones locales violentas en las puntas de corte.

La tecnología moderna de corte y apilamiento utiliza uniones a inglete a $45^\circ$ con desfase progresivo de láminas denominado Step-Lap (generalmente configuraciones de 5 a 7 escalones por ciclo). El escalonamiento distribuye el entrehierro de corte longitudinalmente a lo largo de un gradiente axial:

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

Al desfasar las juntas una distancia espacial $\Delta x \approx 3.5 - 7.0 mm$ entre láminas contiguas, el flujo magnético no se ve obligado a saltar perpendicularmente a través de un único entrehierro de alta reluctancia hacia la lámina adyacente, sino que se reparte progresivamente a lo largo de una sección geométrica efectiva equivalente mucho mayor:

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

La implementación del diseño Step-Lap reduce las pérdidas totales de unión en los yugos hasta en un $25\% - 30\%$ , mitiga la potencia reactiva de excitación en un $40\%$ y reduce los niveles de emisión acústica del transformador en $3 a 6 dB(A)$ en comparación con ensamblajes convencionales.

Refinamiento de Dominios Magnéticos por Rayo Láser (Laser Scribing)

En aceros CRGO de alta permeabilidad (Hi-B), el tamaño de los granos cristalográficos es grande ( $\sim 10 - 30 mm$ ), lo cual favorece la orientación Goss casi perfecta ( $\le 3^\circ$ ), pero resulta contraproducente para el comportamiento dinámico de pérdidas debido a que el espaciamiento entre paredes de Bloch de $180^\circ$ ( $2d_0$ ) crece proporcionalmente con la raíz cuadrada del tamaño de grano:

2d_0 = \sqrt{\frac{\gamma_w \cdot Lgrain}{1.38 \cdot \mu_0 \cdot M_s^2}}

Un ancho de dominio elevado incrementa drásticamente la velocidad instantánea de desplazamiento de la pared $v_w(t)$ , lo que dispara las pérdidas anómalas de Foucault ( $Pexc$ ).

La técnica de refinamiento de dominios por láser (Domain Refinement o Laser Scribing) proyecta un haz focalizado de láser continuo (Nd:YAG o fibra de alta potencia) transversalmente a la dirección de laminación a intervalos espaciales regulares ( $p \approx 4 - 8 mm$ ). La energía térmica del pulso láser induce un calentamiento ultrarrápido y localizado sin ablación del material, generando un micro-campo de tensiones elásticas de compresión permanente sub-superficiales debido al gradiente térmico de enfriamiento:

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

Estas bandas periódicas de tensión interna actúan como barreras de energía magnetoelástica artificiales que fragmentan los dominios magnéticos anchos de $180^\circ$ , reduciendo su espaciamiento medio $2d_0$ a la mitad o un tercio de su dimensión original. Como resultado directo:

P_{exc, lasered } = P_{exc, base } \left( \frac{d_{0, lasered }}{d_{0, base }} \right)^\beta \implies \Delta P_{ total } \approx -10\% a -15\%

Es imperativo distinguir entre el tratamiento láser de tipo no térmicamente resistente (el cual pierde el efecto de refinamiento tras recocidos de alivio de tensiones a $T > 500^\circ C$ requeridos en núcleos enrollados) y los métodos mecánicos de grabado químico/deformación plástica controlada con plasma (recocibles), diseñados para retener el anclaje de microdominios tras tratamientos térmicos a $800^\circ C$ .

Implementación y Modelado Computacional con Vexten Suite

La cuantificación rigurosa de las pérdidas en el núcleo y su acoplamiento multifísico con la red eléctrica se resuelve operativamente en Vexten Suite integrando los módulos de flujo de armónicos no lineales (según IEEE 519 / IEEE C57.110), cálculo térmico dinámico acoplado (IEC 60076-7) y simulación de cortocircuito / regímenes transitorios magnéticos (IEC 60909 / IEEE 141).

Algoritmo de Corrección por Distorsión Armónica y Factor K

En presencia de espectros armónicos severos procedentes de cargas no lineales o generadores renovables, el módulo Vexten PowerFlow & Harmonics ejecuta una partición de pérdidas dinámica recalculando el factor de pérdidas del núcleo de acuerdo con el tensor de coeficientes espectrales normalizados:

F_{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}

El software vincula este factor con la reducción de capacidad nominal (derating) del transformador mediante el cálculo del Factor $K$ y el Factor Armónico de Pérdidas de Foucault ( $FHL$ ), ajustando la corriente admisible máxima para garantizar que la temperatura del punto más caliente del devanado y del núcleo no exceda los límites de la clase térmica de aislamiento:

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

Procedimiento de Modelado Multifísico en Vexten Suite

Definición del Espectro de Excitación Magnética: En el módulo de calidad de potencia, importar las mediciones de campo o espectros típicos IEEE de armónicos de tensión ( $V_h$ ) y la componente continua residual ( $Idc$ ) proveniente del flujo de potencia acoplado con fuentes renovables.
Asignación de Material Magnético en la Base de Datos: Seleccionar en las propiedades de máquina el grado específico de acero CRGO (e.g., M3, 23Q090 o 20Q080-L con refinamiento láser). Vexten Suite carga automáticamente los parámetros microscópicos calibrados: conductividad $\sigma$ , coeficiente de histéresis $k_h$ , espesor $d$ , densidad $\rho_m$ y el parámetro de Bertotti $V_0$ .
Configuración de Topología de Unión: Seleccionar en los parámetros constructivos el tipo de unión en yugos: Butt-Lap $90^\circ$ , Step-Lap de 5 escalones o Step-Lap de 7 escalones. El software introduce automáticamente el multiplicador de dispersión interfacial de reluctancia ( $\mathcal{K}_{joint}$ ).
Ejecución del Análisis Térmico Transitorio Acoplado: El solver acopla las pérdidas electromagnéticas totales $P_{ total }(x,y,z)$ como fuente de calor volumétrica en el solucionador de redes térmicas conforme a IEC 60076-7, determinando la distribución tridimensional de temperaturas en los yugos superiores, inferiores y columnas del transformador.
Validación de Límites y Diagnóstico Forense Preventivo: El sistema compara los gradientes de temperatura local y densidades de pérdidas frente a los umbrales de la norma IEEE C57.104. Si la temperatura en el núcleo excede los $130^\circ C$ continuos, se activan alertas automáticas de derating de potencia y se calcula la tasa teórica estimada de desprendimiento de gases de falla ( $C _2 H _4$ , $CH _4$ , $H _2$ ), permitiendo establecer planes de mantenimiento basados en condición (CBM) de alta precisión física.

Mediante esta arquitectura computacional integral, los ingenieros de diseño, comisionamiento y confiabilidad pueden predecir con exactitud el envejecimiento térmico del paquete magnético, eliminar riesgos de fallas catastróficas por corrientes circulantes interlaminares y optimizar la selección de materiales magnéticos avanzados para transformadores de ultra-alta eficiencia energética.