Ferrorresonancia en Transformadores y Coordinación gPV en Plantas Automotrices (IEEE 1584 / NFPA 70E)

Fenomenología No Lineal de la Ferrorresonancia en Transformadores de Plantas Automotrices

La ferrorresonancia es un fenómeno de resonancia no lineal que difiere fundamentalmente de la resonancia lineal clásica. En un circuito lineal RLC, la resonancia ocurre a una frecuencia angular única determinada por la interacción de elementos reactivos constantes:

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

En sistemas de distribución de plantas automotrices, caracterizados por densas redes de cables subterráneos de media tensión con alta capacitancia shunt ( $C_0$ ) y transformadores de distribución con núcleos ferromagnéticos saturables, la inductancia magnetizante ( $L_m$ ) no es constante. Esta inductancia es una función no lineal altamente sensible al flujo magnético acoplado ( $\phi$ ), descrita por la relación entre la corriente de magnetización ( $i_m$ ) y el flujo:

i_m(\phi) = a\cdot\phi + b\cdot\phi^n

Donde $a$ y $b$ son coeficientes de ajuste magnético dependientes de la geometría y material del núcleo, y $n$ es un exponente impar (típicamente $n \ge 7$ para aceros al silicio de grano orientado modernos). Esta no linealidad introduce múltiples puntos de operación estables e inestables para una misma frecuencia de excitación industrial (50 Hz o 60 Hz). El comportamiento dinámico del sistema se modela mediante la ecuación diferencial no lineal de Duffing:

\frac{d^2\phi}{dt^2} + \frac{1}{R_p C_0}\frac{d\phi}{dt} + \frac{1}{C_0}i_m(\phi) = \omega E_m \cos(\omega t + \theta)

Donde $R_p$ representa las pérdidas en el cobre y las pérdidas del núcleo (parásitas e histéresis), y $E_m$ es la amplitud de la tensión de excitación. La transición de un estado operativo normal de baja inducción a un estado ferrorresonante de alta inducción ocurre a través de bifurcaciones en el espacio de fases.

Mecanismos de Bifurcación y Modos de Ferrorresonancia

La ferrorresonancia se clasifica según la topología de su órbita en el espacio de fases y el contenido espectral de las variables de estado:

Modo Fundamental (Periodo 1): Las variables de tensión y flujo son periódicas con el mismo periodo $T = 1/f$ de la fuente de excitación. Se caracteriza por sobretensiones severas y distorsiones simétricas de la onda de corriente de magnetización.
Modo Subarmónico: Las órbitas en el espacio de fases se cierran después de un múltiplo entero del periodo de la fuente ( $nT$ , típicamente $3T$ o 1/3 de la frecuencia fundamental). Este modo es común en sistemas con baja amortiguación y puede causar vibraciones mecánicas severas debido a las fuerzas electrodinámicas en los devanados.
Modo Cuasi-periódico: Se presenta cuando la trayectoria en el espacio de fases no es cerrada, sino que envuelve un toroide. El espectro de frecuencia muestra frecuencias no armónicas superpuestas a la fundamental.
Modo Caótico: Caracterizado por una extrema sensibilidad a las condiciones iniciales. Las trayectorias en el espacio de fases se comportan como atractores extraños, con un espectro de frecuencia continuo de ruido de fondo y picos desordenados. Este modo somete al aislamiento del transformador a un estrés térmico y dieléctrico impredecible.

\lambda_{\text{Lyapunov}} > 0

La existencia de un exponente de Lyapunov positivo ( $\lambda_{\text{Lyapunov}}$ ) confirma la presencia del modo caótico, donde el sistema pierde predictibilidad determinista a largo plazo, acelerando la degradación del aislamiento sólido (papel Kraft) y líquido (aceite mineral o éster sintético).

---

Análisis Forense de Fallas por Ferrorresonancia y Transitorios de Conmutación

En el entorno de manufactura robotizada de la industria automotriz, la continuidad de servicio y la calidad de la energía son críticas. Las fallas catastróficas en subestaciones unitarias de media tensión (MT) a menudo son causadas por transitorios de conmutación mal amortiguados que desencadenan ferrorresonancia. Los escenarios típicos de iniciación forense incluyen:

Operación Monofásica de Aparamenta (Fusibles o Seccionadores): La apertura o cierre no simultáneo de los polos de un interruptor o la fusión de un solo fusible de MT debido a una falla aguas abajo deja al transformador alimentado por una o dos fases. Las fases abiertas quedan acopladas capacitivamente a través de las capacidades mutuas de los cables de alimentación ( $Cg$ y $Cm$ ).
Rotura de Conductores de Alimentación: Un cable subterráneo con pantalla metálica a tierra que sufre una falla de aislamiento y posterior corte físico de un conductor crea una topología serie-paralelo ideal para la resonancia entre la capacitancia de la pantalla del tramo desconectado y la reactancia magnetizante no lineal del transformador.
Sistemas de Neutro Aislado o Puesta a Tierra Ineficaz: Los transformadores con conexión estrella-triángulo o triángulo-estrella alimentados por sistemas con neutro flotante son propensos a desplazamientos severos del punto neutro durante fallas a tierra monofásicas, lo que satura los núcleos magnéticos de las fases sanas.

Degradación Térmica y Dieléctrica del Aislamiento

Durante la ferrorresonancia, la corriente de magnetización que fluye por los devanados del transformador puede incrementarse hasta 100 veces su valor nominal de vacío, convirtiéndose en corrientes de inrush continuas y repetitivas. Esto genera pérdidas por efecto Joule localizadas en los devanados de MT:

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

Donde $Rac(f)$ es la resistencia de corriente alterna corregida por el efecto pelicular (skin effect) y proximidad a altas frecuencias armónicas. Simultáneamente, el flujo de dispersión magnética en condiciones de saturación profunda abandona el núcleo de acero silicio y penetra en las estructuras metálicas de soporte, las paredes del tanque y los pernos de sujeción, induciendo corrientes de Foucault masivas que causan puntos calientes (hotspots) locales que superan rápidamente los 200 °C.

El aceite dieléctrico sufre pirólisis térmica acelerada, generando gases inflamables como hidrógeno ( $H_2$ ), etileno ( $C_2H_4$ ) y acetileno ( $C_2H_2$ ), detectables mediante Análisis de Gases Disueltos (DGA) según la norma IEC 60599. El estrés dieléctrico se ve agravado por sobretensiones transitorias con altas tasas de aumento de tensión ( $dv/dt$ ), las cuales dañan la rigidez dieléctrica del aislamiento de las espiras extremas del devanado debido a la distribución no uniforme del potencial capacitivo a lo largo de la bobina.

Degradación de Descargadores de Sobretensión (MOV)

Los descargadores de sobretensión de óxido metálico (MOV) instalados para proteger el transformador contra sobretensiones atmosféricas y de maniobra se convierten en víctimas colaterales durante la ferrorresonancia. Un MOV está diseñado para absorber energía transitoria de corta duración (microsegundos a milisegundos). Sin embargo, la ferrorresonancia sostiene sobretensiones temporales (TOV) durante segundos o minutos. La energía disipada por el varistor excede su capacidad de absorción térmica admisible:

EMOV = \int0^{tfall} vMOV(t) \cdot iMOV(t) \, dt > E_{max\_admisible}

Esto provoca un embalamiento térmico (thermal runaway). Al aumentar la temperatura del varistor de ZnO, su corriente de fuga en estado de conducción normal aumenta de forma exponencial debido a su coeficiente de temperatura negativo, lo que reduce su resistencia óhmica interna hasta su destrucción física, explosión y cortocircuito franco en las barras de MT.

---

Coordinación de Protecciones en Corriente Continua: Fusibles gPV en Sistemas de Almacenamiento y Generación Fotovoltaica Integrada

Las plantas automotrices modernas integran recursos energéticos distribuidos (DER), como plantas fotovoltaicas en cubiertas y sistemas de almacenamiento de energía por baterías (BESS), para alimentar cargas críticas y suavizar los picos de demanda de las líneas de soldadura robotizada. La protección de estos sistemas de corriente continua (DC) requiere una coordinación precisa mediante fusibles de categoría de utilización gPV (diseñados específicamente para protección de sistemas fotovoltaicos según IEC 60269-6).

Física de la Interrupción del Arco Eléctrico en Corriente Continua

A diferencia de los sistemas de corriente alterna (AC), donde la corriente cruza por cero de forma natural dos veces por ciclo, facilitando la extinción del arco eléctrico, en los sistemas DC la corriente no presenta cruces por cero. La extinción del arco eléctrico exige que el dispositivo de protección genere activamente una tensión de arco ( $Varco$ ) sustancialmente superior a la tensión nominal del sistema de corriente continua ( $Vsys$ ):

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

Para que la corriente decrezca ( $di/dt < 0$ ), es imperativo que:

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

La arena de sílice de alta pureza que rellena el cuerpo del fusible gPV actúa como un medio de enfriamiento y desionización del canal de plasma del arco. Al fundirse el elemento conductor de plata pura (Ag) bajo condiciones de sobrecorriente, el arco eléctrico se estira y segmenta en los cuellos restrictivos del fusible. La energía del arco funde la sílice circundante, formando una estructura vítrea no conductora llamada fulgurita, que absorbe la energía térmica e interrumpe el flujo de corriente.

Coordinación de Fusibles gPV bajo Condiciones de Cortocircuito Extremas

En el diseño de la protección DC para plantas automotrices con BESS y generación fotovoltaica integrada, la coordinación se complica debido a los dos perfiles de corriente de falla extremos del sistema:

Corrientes de Falla de Baja Magnitud (Típicas de Strings Fotovoltaicos): Los módulos fotovoltaicos son fuentes de corriente limitada por la irradiancia solar. Las corrientes de cortocircuito típicas ( $Isc$ ) son solo de un 10% a un 20% superiores a la corriente de máxima potencia ( $Imp$ ). Los fusibles gPV deben ser capaces de operar de forma segura ante estas corrientes de falla mínimas ( $Imin \approx 1.35 \cdot I_n$ a $2 \cdot I_n$ ) en un tiempo inferior al que comprometa la integridad de los conductores de DC y los módulos.
Corrientes de Falla de Alta Magnitud y Elevada Constante de Tiempo (Típicas de BESS): Los bancos de baterías de ion-litio presentan una resistencia interna extremadamente baja, lo que da lugar a corrientes de cortocircuito presuntas muy elevadas (decenas de kiloamperios) con constantes de tiempo de red muy rápidas ( $\tau = L/R$ ) o lentas según la longitud del embarrado. El fusible gPV debe poseer una capacidad de ruptura extrema (típicamente hasta 100 kA DC) y limitar la energía pasante (integral de Joule, $I^2t$ ).

La integral de Joule de pre-arco ( $I^2t_{pre}$ ) y la integral de Joule de operación total ( $I^2t_{op}$ ) definen los límites de selectividad entre fusibles instalados en serie (por ejemplo, fusible de string vs. fusible de subgrupo vs. fusible principal del inversor):

I^2t_{pre} = \int0^{tpre} i^2(t) \, dt = Constante para un material y geometría de fusible dados

Para garantizar una coordinación selectiva total entre dos fusibles en serie (aguas arriba y aguas abajo), se debe cumplir rigurosamente la siguiente relación de energía en todo el espectro de corrientes de cortocircuito:

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

Si esta condición no se cumple, ambos fusibles se fundirán simultáneamente (pérdida de selectividad), provocando un apagón innecesario de secciones completas de la microrred de la planta automotriz.

---

Tabla Comparativa de Estándares, Parámetros Eléctricos y Límites Críticos

La siguiente tabla detalla los parámetros de diseño y límites operativos normalizados según estándares internacionales para la prevención de ferrorresonancia y la correcta coordinación de protecciones DC gPV.

Parámetro / Fenómeno	Estándar de Referencia	Límite Crítico Operativo	Consecuencia Dieléctrica / Térmica	Estrategia de Mitigación / Diseño
Sobretensión Temporal (TOV) por Ferrorresonancia	IEEE C57.105 / IEC 60076-1	$Vpeak \ge 1.5 p.u.$ por más de 2 segundos.	Falla del aislamiento entre espiras, degradación térmica del papel Kraft, arco interno.	Instalación de resistencias de amortiguamiento en el devanado en triángulo abierto.
Distorsión Armónica Total de Tensión ( $THD_v$ )	IEEE 519 / IEC 61000-4-30	$THD_v > 8.0\%$ en barras de MT (1 kV a 69 kV).	Aumento de pérdidas por histéresis y corrientes parásitas en el núcleo del transformador.	Filtros de armónicos sintonizados y reactores de desintonización en bancos de capacitores.
Capacidad de Ruptura de Fusibles gPV	IEC 60269-6 / UL 248-19	$I_{sc\_presunta} > I_{bc\_fusible}$ (Típicamente < 10 kA a 50 kA).	Ruptura física violenta del fusible, arco persistente, riesgo de incendio en la subestación.	Selección de fusibles con alta capacidad de ruptura certificada ( $\ge 100 kA DC$ ).
Constante de Tiempo de Falla DC ( $\tau$ )	IEC 60269-1 / IEC 60949	$\tau = L/R > 15 ms$ en cortocircuito DC.	Dificultad extrema para extinguir el arco eléctrico; prolongación del tiempo de despeje.	Uso de fusibles gPV ensayados bajo constantes de tiempo severas ( $\tau \approx 10-20 ms$ ).
Factor de Desclasificación Térmica de Cables	IEC 60287 / NEC Art. 310	Temperatura de operación en conductor XLPE $\ge 90 °C$ .	Envejecimiento térmico acelerado del aislamiento de polietileno reticulado, perforación.	Cálculo de ampacidad considerando armónicos mediante el factor de corrección de pérdidas de Vexten.
Rigidez Dieléctrica del Aceite Dieléctrico	IEC 60156 / ASTM D1816	Tensión de ruptura dieléctrica $< 30 kV$ .	Formación de arcos internos destructivos bajo sobretensiones transitorias ordinarias.	Tratamiento de termovacío, desgasificación y regeneración activa del aceite mineral.

---

Estrategias de Mitigación y Diseño de Ingeniería de Detalle

Para evitar de forma proactiva la aparición de ferrorresonancia en las subestaciones de la planta automotriz, se deben implementar soluciones de ingeniería de detalle basadas en el análisis de parámetros físicos.

Diseño de Circuitos de Amortiguamiento (Damping Loops)

Una de las soluciones más efectivas para mitigar la ferrorresonancia en transformadores de medida de tensión (VTs) y transformadores de distribución conectados en estrella con neutro a tierra consiste en la instalación de una resistencia de amortiguamiento ( $R_d$ ) conectada en los terminales de un devanado secundario auxiliar conectado en triángulo abierto (esquema de neutro roto):

R_d \le \frac{1}{3 \cdot \omega \cdot C_0}

Donde $C_0$ representa la capacitancia homopolar por fase de la red de cables de media tensión conectada al transformador. Esta resistencia actúa disipando la energía de las oscilaciones ferrorresonantes cuando aparece una tensión de secuencia cero debido a una asimetría o transitorio de conmutación. La potencia térmica nominal ( $PRd$ ) de la resistencia de amortiguamiento debe calcularse para soportar la tensión de cortocircuito homopolar temporal sin quemarse:

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

Criterios de Conmutación Tripolar Simultánea

Se debe prohibir estrictamente el uso de fusibles de media tensión o seccionadores unipolares aguas arriba de transformadores de distribución propensos a la ferrorresonancia. En su lugar, se deben especificar interruptores automáticos de corte tripolar simultáneo (SF6 o vacío) controlados por relés de protección numéricos de alta velocidad. Esto asegura que las tres fases se energicen o desenergicen con un desfase temporal máximo entre polos inferior a 2 milisegundos:

\Delta tpolos \le 2 ms

Esta simultaneidad impide la formación de circuitos resonantes serie-paralelo asimétricos que caracterizan a las fases abiertas.

Sizing y Coordinación de Fusibles gPV en el Lado DC

Para el dimensionamiento riguroso de los fusibles gPV en los inversores y strings fotovoltaicos de la planta automotriz, se deben seguir las siguientes ecuaciones de coordinación:

La tensión nominal del fusible ( $V_{n\_fuse}$ ) debe corregirse según la temperatura ambiente de la subestación y la tensión máxima de vacío del string fotovoltaico ( $Voc$ ) a la temperatura mínima de diseño del sitio ( $Tmin$ ):

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

Donde $\beta_{voc}$ es el coeficiente de temperatura de la tensión de circuito abierto del panel fotovoltaico ( $\%/ °C$ ).

La corriente nominal del fusible ( $I_{n\_fuse}$ ) debe considerar la desclasificación por temperatura ambiente ( $Ktemp$ ) y por altitud ( $Kalt$ ) si la planta se ubica a más de 2000 metros sobre el nivel del mar:

I_{n\_fuse} \ge \frac{1.4 \cdot I_{sc\_string}}{Ktemp \cdot Kalt}

El factor de seguridad de 1.4 garantiza que el fusible no sufra envejecimiento prematuro debido a corrientes cíclicas provocadas por las fluctuaciones de irradiancia solar y las demandas pulsantes de energía del sistema de soldadura robotizada de la planta.

---

Modelado y Simulación Avanzada mediante Vexten Suite

La suite de ingeniería eléctrica Vexten Suite es fundamental para predecir, analizar y resolver problemas de ferrorresonancia y coordinación de protecciones en plantas industriales complejas. A continuación, se presenta la metodología de modelado y análisis utilizando los módulos avanzados de esta herramienta.

Módulo de Cortocircuito (IEC 60909 / IEEE 141)

Para determinar las corrientes de falla máximas y mínimas en los puntos de acoplamiento común (PCC) y en las barras de distribución de corriente continua de la microrred de la planta automotriz, se utiliza el motor de cálculo de cortocircuito de Vexten Suite. El software calcula de forma simultánea los parámetros de cortocircuito según las normas internacionales:

Corriente de cortocircuito pico ( $i_p$ ): Crucial para determinar los esfuerzos mecánicos en los embarrados y la capacidad de cierre de la aparamenta.
Corriente de cortocircuito térmica equivalente ( $Ith$ ): Define el estrés térmico soportado por los conductores y fusibles durante el tiempo de despeje de la falla.

Vexten Suite modela de manera precisa la impedancia subtransitoria de los generadores y la contribución de los motores de inducción y síncronos de la línea de ensamblaje automotriz, aplicando la matriz de impedancias de barra $[Zbus]$ :

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

En el dominio de corriente continua, el software resuelve las ecuaciones diferenciales de descarga de los bancos de baterías (BESS) y capacitores de enlace DC de los variadores de frecuencia (VFD), permitiendo parametrizar la resistencia interna no lineal de las celdas de litio y la inductancia de la red de barras.

Módulo de Dimensionamiento de Cables y Desclasificación Armónica (IEC 60287 / NEC 310)

Las líneas de producción robotizadas contienen cientos de servomotores y sistemas de soldadura por resistencia que inyectan corrientes armónicas masivas en la red de distribución de baja y media tensión. El módulo de dimensionamiento de cables de Vexten Suite implementa los algoritmos iterativos de la norma IEC 60287 para el cálculo de la temperatura de estado estable del conductor.

El software realiza una desclasificación armónica avanzada considerando el incremento de la resistencia en corriente alterna ( $Rac$ ) debido al efecto pelicular (skin effect) y al efecto proximidad en presencia de armónicos de alta frecuencia:

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

Donde $y_s(h)$ e $y_p(h)$ son los factores de efecto pelicular y proximidad para el armónico de orden $h$ , respectivamente. Vexten Suite calcula automáticamente el factor de desclasificación armónica ( $Kharmonic$ ):

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

Este factor corrige la ampacidad admisible del conductor de cobre con aislamiento XLPE de los circuitos alimentadores de los transformadores, previniendo sobrecalentamientos que degraden el aislamiento y actúen como iniciadores de fallas a tierra que desencadenen ferrorresonancia.

Módulo de Corrección del Factor de Potencia y Mitigación de Resonancia

Para evitar la resonancia armónica paralela entre los bancos de capacitores de corrección del factor de potencia y la inductancia de cortocircuito del transformador de la planta automotriz, Vexten Suite permite diseñar filtros de armónicos desintonizados (detuned filters). El software calcula la frecuencia de sintonización del reactor de bloqueo ( $f_r$ ) requerida para desplazar el punto de resonancia paralela fuera de las frecuencias de los armónicos dominantes inyectados por los rectificadores de los robots (típicamente los armónicos 5º, 7º, 11º y 13º):

f_r = f_1 \cdot \sqrt{\frac{1}{p}}

Donde $p$ es el factor de desintonización (típicamente $p = 7\%$ para una frecuencia de sintonización de 189 Hz en redes de 50 Hz, o $p = 14\%$ para evitar la amplificación del 3º armónico). El entorno gráfico de simulación de Vexten Suite proporciona las curvas de respuesta en frecuencia de la impedancia del sistema vista desde cualquier nodo de la planta:

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

Esto permite validar visual y analíticamente que no existan picos de alta impedancia en las frecuencias críticas de conmutación, eliminando la posibilidad de sobretensiones ferrorresonantes transitorias provocadas por la conmutación de los bancos de capacitores.

A través de su motor de simulación de transitorios electromagnéticos (EMT) integrado, Vexten Suite permite modelar la histéresis magnética del núcleo del transformador mediante el modelo de Jiles-Atherton. Esto proporciona a los ingenieros de diseño automotriz la capacidad de predecir la aparición de ferrorresonancia bajo diferentes escenarios de maniobra de disyuntores, asegurando la robustez operativa de la infraestructura eléctrica de la planta antes de su construcción o modificación en campo.