Perdas por Histerese e Foucault em Aço Silício 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

Termodinâmica Eletromagnética e Microestrutura do Aço Silício CRGO

O núcleo magnético de um transformador de potência opera como o circuito de acoplamento de fluxo fundamental entre enrolamentos, operando sob regimes cíclicos não lineares de excitação eletromagnética. A eficiência global e a densidade de potência do transformador são condicionadas de forma determinante pelo comportamento microestrutural e cristalográfico do material ferromagnético utilizado. O padrão industrial predominante para transformadores de média, alta e extra-alta tensão (EHV/UHV) é o aço elétrico de grão orientado laminado a frio (Cold-Rolled Grain-Oriented Silicon Steel, CRGO), otimizado por meio de tratamentos termomecânicos para alinhar seus eixos cristalográficos de fácil magnetização na direção de laminação.

O aço CRGO apresenta uma estrutura cristalina cúbica de corpo centrado (BCC, rede α\alpha-Fe), onde o eixo cristalográfico 100\langle 100 \rangle constitui a direção de mínima energia de anisotropia magnetocristalina, enquanto as direções 110\langle 110 \rangle e 111\langle 111 \rangle representam eixos de anisotropia intermediária e dura, respectivamente. Por meio de uma sequência rigorosa de laminação a frio e recozimentos de recristalização secundária (inibida mediante precipitados de sulfeto de manganês, MnSMnS, ou nitreto de alumínio, AlNAlN), desenvolve-se a denominada textura de Goss, caracterizada pela orientação cristalográfica {110}001\{110\}\langle 001 \rangle. Nessa configuração, os planos {110}\{110\} situam-se paralelos à superfície da chapa e as direções 001\langle 001 \rangle orientam-se paralelamente à direção longitudinal de laminação (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)

Onde EaE_a é a densidade de energia de anisotropia magnetocristalina (J/m3J/m^3), K14.8×104 J/m3K_1 \approx 4.8 \times 10^4 \ J/m^3 e K21.5×104 J/m3K_2 \approx 1.5 \times 10^4 \ J/m^3 são as constantes de anisotropia do ferro à temperatura ambiente (298 K298 \ K), e α1,α2,α3\alpha_1, \alpha_2, \alpha_3 representam os cossenos diretores do vetor de magnetização M\mathbf{M} em relação aos eixos da célula unitária. O desalinhamento angular médio da textura de Goss em relação à direção de laminação define a qualidade do aço: em aços CRGO convencionais, a dispersão angular oscila entre 55^\circ e 77^\circ, enquanto em aços de alta permeabilidade (Hi-B) reduz-se para faixas de 22^\circ a 33^\circ.

A adição de silício (SiSi) em concentrações nominais de 3.0%3.0\% a 3.4%3.4\% em peso incrementa a resistividade volumétrica elétrica (ρ\rho) desde aproximadamente 10×108 Ωm10 \times 10^{-8} \ \Omega\cdot m (para ferro puro) até (4550)×108 Ωm(45 - 50) \times 10^{-8} \ \Omega\cdot m. Esse incremento de resistividade atenua drasticamente o transporte livre de cargas parasitas. No entanto, concentrações superiores a 3.5%3.5\% em peso induzem fragilidade mecânica severa devido à formação de fases ordenadas intermetálicas Fe3SiFe_3Si e FeSiFeSi (B2B2 e D03D0_3), impedindo os processos de conformação a frio em escala industrial. Adicionalmente, o silício reduz a constante de magnetoestrição de saturação longitudinal (λ100\lambda_{100}), minimizando as deformações elásticas induzidas pelo campo magnético e a consequente emissão de ruído acústico.

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

A estrutura de domínios magnéticos (domínios de Weiss) no CRGO é composta por domínios principais orientados a 180180^\circ separados por paredes de Bloch, e domínios de fechamento a 9090^\circ nas imediações de descontinuidades superficiais e inclusões. A largura de domínio característica (2d02d_0) depende do equilíbrio entre a energia magnetostática superficial e a energia elástica da parede do domínio. A aplicação de revestimentos dielétricos com tensão superficial (como o recobrimento de fosfatos e silicatos, denominado comercialmente Carlite) introduz uma tração biaxial permanente na chapa (σt410 MPa\sigma_t \approx 4 - 10 \ MPa), reduzindo a largura média dos domínios a 180180^\circ e eliminando domínios secundários transversais de fechamento, diminuindo radicalmente as perdas por correntes parasitas anômalas.

Tripartição Clássica e Espectro de Perdas: O Modelo de Bertotti

Sob excitação magnetodinâmica cíclica, a dissipação volumétrica de energia por unidade de tempo no núcleo é quantificada por meio da integral fechada sobre o ciclo do laço de histerese dinâmico. A formulação geral de perdas específicas de potência (PtP_t, em W/kgW/kg) é regida pela teoria estatística de perdas eletromagnéticas formulada por Giorgio Bertotti, a qual supera as limitações das formulações empíricas clássicas de Steinmetz ao decompor o fenômeno em três mecanismos dissipativos desacoplados fisicamente:

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}

Onde ff representa a frequência fundamental de excitação (HzHz), BpB_p é o valor de pico da indução magnética (TT), σ\sigma é a condutividade elétrica do material (S/mS/m), dd é a espessura da chapa laminada (mm), ρm\rho_m é a densidade de massa volumétrica do aço (kg/m3kg/m^3, tipicamente 7650 kg/m37650 \ kg/m^3), khk_h é o coeficiente de histerese quase-estática, α\alpha é o expoente de Steinmetz (frequentemente 1.6α2.01.6 \le \alpha \le 2.0), e kexckexc é o parâmetro microscópico de perdas em excesso.

Perdas por Histerese Quase-Estática

As perdas por histerese estática (PhP_h) originam-se nos processos microscópicos de dissipação termodinâmica irreversível quando as paredes dos domínios magnéticos se deslocam através de um meio ferromagnético com defeitos de rede cristalina, tais como discordâncias, vacâncias, contornos de grão, tensões residuais e inclusões não metálicas (e.g., precipitados de SiO2SiO_2 ou Fe3CFe_3C). Durante o movimento de uma parede de Bloch de 180180^\circ, esta fica ancorada mecanicamente nesses centros de ancoramento (pinning sites).

Para liberar a parede de domínio do poço de potencial de ancoramento, faz-se necessário um incremento do campo magnético externo até alcançar o campo crítico local. Uma vez superada essa barreira energética, a parede avança descontinuamente em alta velocidade em direção ao poço de potencial seguinte por meio de um salto microscópico não linear, fenômeno conhecido como efeito Barkhausen. A energia dissipada nesse salto converte-se irreversivelmente em fônons (calor) na rede cristalina:

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

A área do laço de histerese quase-estático (whw_h) depende exclusivamente da estrutura metalúrgica do material e é independente da frequência de variação temporal do fluxo para regimes onde o tempo de relaxação do domínio é infinitamente menor que o período da excitação. A modelagem matemática de PhP_h na presença de formas de onda complexas ou saturação não linear requer a formulação do operador contínuo de 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

Onde γ^αβ\hat{\gamma}_{\alpha\beta} é um operador histerético elementar biestável com limites de comutação α\alpha (de ativação) e β\beta (de desativação), e μ(α,β)\mu(\alpha, \beta) é a função de distribuição de densidade de Preisach específica do material CRGO ensaiado, identificável experimentalmente por meio de famílias de curvas de reversão de primeira ordem (FORC, First-Order Reversal Curves).

Perdas por Correntes Parasitas Clássicas (Foucault)

As perdas clássicas por correntes de Foucault (PclPcl) derivam rigorosamente da aplicação macroscópica das equações de Maxwell a um condutor homogêneo e isotrópico contínuo de geometria planar finita. Considere-se uma chapa de aço de espessura dd na direção do eixo zz (desde z=d/2z = -d/2 até z=+d/2z = +d/2), largura infinita no eixo yy, e comprimento no eixo xx. Assumindo uma permeabilidade magnética efetiva μ\mu, uma condutividade volumétrica σ=1/ρ\sigma = 1/\rho, e um vetor de densidade de fluxo magnético paralelo à superfície: 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}

Assumindo que a espessura dd é substancialmente menor que a profundidade de penetração eletromagnética (skin depth) δ=2/(ωσμ)\delta = \sqrt{2/(\omega \sigma \mu)}, a indução magnética distribui-se de maneira espacialmente uniforme através da seção transversal da lâmina: Bx(z,t)Bpsin(ωt)B_x(z,t) \approx B_p \sin(\omega t). Integrando em relação a zz a partir do plano neutro central (z=0z = 0, onde Ey=0E_y = 0 por simetria):

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)

A densidade de corrente induzida J(z,t)=σE(z,t)\mathbf{J}(z,t) = \sigma \mathbf{E}(z,t) dissipa uma potência instantânea por efeito Joule cuja densidade volumétrica local é 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 sobre a espessura total da chapa e calculando a média no domínio temporal durante um período 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]

Dividindo pela densidade de massa volumétrica ρm\rho_m, obtém-se a formulação específica normalizada:

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]

Essa dedução expõe a dependência quadrática direta em relação à espessura dd da chapa e à frequência fundamental ff, justificando a necessidade de reduzir os calibres comerciais de 0.35 mm0.35 \ mm (M4) para 0.23 mm0.23 \ mm ou 0.18 mm0.18 \ mm (graus tratados a laser de ultrabaixas perdas).

Perdas Anômalas ou em Excesso (Excess Losses)

Historicamente, a soma direta das perdas por histerese quase-estática e das perdas clássicas de Foucault resultava significativamente inferior às perdas totais medidas experimentalmente pelo método do wattímetro em quadro de Epstein (Ptotal,med>Ph+PclP_{total,med} > P_h + Pcl). Essa discrepância energética foi formalmente denominada anomalia de perdas.

O modelo microscópico de Pry e Bean (1958) demonstrou que o fluxo magnético não varia uniformemente em todo o contínuo da lâmina, mas se concentra exclusivamente dentro dos domínios magnéticos e se transfere por meio do deslocamento localizado das paredes de Bloch de 180180^\circ. A velocidade de deslocamento de uma parede de domínio vw(t)v_w(t) sob excitação senoidal induz um gradiente de campo elétrico local infinitamente mais pronunciado na vizinhança imediata da parede do que o previsto pela teoria clássica do contínuo:

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

Onde LL representa o espaçamento entre paredes de domínio e BsB_s é a indução de saturação. A microcorrente parasita confinada no microvolume adjacente à parede gera um campo de frenagem magnética local. Bertotti generalizou essa física formulando a teoria estatística dos Objetos Magnéticos Ativos (Magnetic Objects, MO), compreendidos como agregados coerentes de paredes de domínio que interagem com as heterogeneidades locais do meio:

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]

Onde G=0.1356G = 0.1356 é uma constante adimensional associada à dissipação do campo de Foucault em meios infinitos, SS é a área da seção transversal da lâmina perpendicular à direção do fluxo magnético, e V0V_0 é um parâmetro intrínseco do material com dimensões de campo magnético (A/mA/m), o qual quantifica a interação estocástica entre os objetos magnéticos e o campo microscópico de ancoramento (internal pinning field distribution).

Grau CRGO (Classificação AISI / EN) Espessura Nominal dd (mm) Indução Máxima BB a 800 A/m800 \ A/m (T) Perdas Específicas P1.7/50P_{1.7/50} (W/kg) Perdas Específicas P1.7/60P_{1.7/60} (W/kg) Fator de Empilhamento (%)
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 Convencional (27Q110) 0.27 1.90 - 1.93 0.95 - 1.05 1.24 - 1.38 95.0\ge 95.0
Hi-B Alta Permeabilidade (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
Metal 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

Efeitos Não Lineares, Distorção Harmônica e Componente Contínua (DC Bias)

Em sistemas de potência contemporâneos, a disseminação massiva da eletrônica de potência (inversores de grande escala para parques fotovoltaicos, elos HVDC, sistemas FACTS e acionamentos de frequência variável) introduz regimes de magnetização não senoidais e tensões com elevado conteúdo harmônico (THDvTHD_v). A resposta não linear da curva de magnetização dinâmica B(H)B(H) transforma essas excitações em perdas energéticas anômalas severas.

Excitação Harmônica e Laços de Histerese Menores

Quando a forma de onda de tensão aplicada contém harmônicos superiores (h=3,5,7,h = 3, 5, 7, \dots), a derivada temporal da indução magnética B/t\partial \mathbf{B}/\partial t apresenta múltiplos cruzamentos por zero no decorrer de um ciclo fundamental. Fisicamente, isso força a reversão transitória do movimento das paredes de domínio magnético antes de completar a polarização macroscópica total, originando laços menores de histerese (minor hysteresis loops) aninhados dentro do laço principal de saturação.

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

Sob tensão senoidal pura, BhVh/(hf1)B_h \propto V_h / (h f_1). Se a tensão contém distorção harmônica com componentes VhV_h, o campo induzido harmônico é aproximado por Bp,hVhhV1Bp,1B_{p,h} \approx \frac{V_h}{h \, V_1} B_{p,1}. Substituindo na formulação de correntes clássicas:

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)

No entanto, as perdas em excesso não seguem esse desacoplamento linear. A interação de múltiplas frequências acelera a velocidade local das paredes de Bloch de acordo com uma integral temporal da derivada de fluxo:

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

Magnetização com Deslocamento de Corrente Contínua (DC Bias)

A injeção de corrente contínua nas fases do transformador — produzida por correntes geomagneticamente induzidas (GIC), topologias de inversores de acoplamento direto ou faltas assimétricas monopolares de linhas de transmissão HVDC com retorno por terra — produz uma translação do ciclo dinâmico sobre o eixo do campo magnetizante (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)

Esse deslocamento assimétrico força o núcleo em direção ao joelho de saturação durante um semiciclo completo (saturação semicíclica). Na região de saturação (B>1.9 TB > 1.9 \ T em CRGO), a permeabilidade diferencial μd=dB/dH\mu_d = dB/dH colapsa drasticamente para valores próximos à permeabilidade do váCuO (μ0=4π×107 H/m\mu_0 = 4\pi \times 10^{-7} \ H/m). Isso induz picos massivos de corrente magnetizante assimétrica (imag(t)imag(t) com fatores de crista superiores a 10), transformando o campo magnético de um vetor confinado dentro do núcleo laminado em um fluxo de dispersão tridimensional de grande magnitude que enlaça elementos estruturais adjacentes (placas de aperto, tanques de transformadores, tirantes de amarração e blindagens magnéticas).

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

O consumo de potência reativa QQ eleva-se em ordens de magnitude, originando severas quedas de tensão nos terminais, harmônicos pares e ímpares de baixa ordem (especialmente h=2h = 2 e h=3h = 3) e aquecimentos térmicos catastróficos por correntes induzidas em peças metálicas estruturais.

Análise Forense de Falhas Térmicas e Degradação Dielétrica Induzida pelo Núcleo

As patologias eletromagnéticas e de degradação térmica originadas no núcleo constituem um mecanismo de falha crítico em transformadores de potência de grande porte. Diferentemente das falhas mecânicas nos enrolamentos induzidas por esforços de curto-circuito radiais e axiais, as falhas de núcleo costumam manifestar-se como processos de degradação eletroquímica e térmica progressiva ao longo de extensos períodos de operação.

Degradação do Isolamento Interlaminar (Revestimento Tipo Carlite)

O isolamento interlaminar superficial das chapas magnéticas (espessura 1.53.0 μm\sim 1.5 - 3.0 \ \mu m, composto por uma base de forsterita Mg2SiO4Mg_2SiO_4 com sobrecamada de fosfato de alumínio AlPO4AlPO_4 e ácido fosfórico) é projetado para suportar tensões interlaminares de alguns poucos volts (1.5 V1.5 \ V a 5.0 V5.0 \ V eficazes conforme a norma IEC 60404-6). Entretanto, sob solicitações térmicas locais contínuas (T>150CT > 150^\circ C), fadiga vibracional por magnetoestrição (100/120 Hz100/120 \ Hz e harmônicos acústicos) ou tensões mecânicas excessivas por sobreaperto durante a montagem dos jugos, o revestimento dielétrico fragmenta-se e pulveriza-se mecanicamente.

A perda do isolamento galvânico entre chapas contíguas gera microssoldas por arco ou contato galvânico direto. Se dois ou mais pontos de contato forem estabelecidos no mesmo pacote de laminação, forma-se uma espira em curto-circuito de grande seção transversal que enlaça o fluxo magnético principal do núcleo. A força eletromotriz induzida no circuito fechado:

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}

Dado que a resistência do circuito ôhmico RloopRloop é da ordem de miliohms (mΩm\Omega), a corrente de falha IfaultIfault atinge amplitudes de centenas a milhares de ampères em regime permanente, dissipando calor localizado por efeito Joule que gera pontos quentes (hot-spots) com temperaturas que ultrapassam rapidamente 500C800C500^\circ C - 800^\circ C.

Degradação do Óleo Dielétrico e Análise Forense por DGA

A transferência de calor extremo do ponto quente do núcleo para o fluido dielétrico circundante (óleo mineral naftênico, parafínico ou ésteres sintéticos) induz a pirólise e a cisão termocatalítica das cadeias moleculares de hidrocarbonetos (CnH2n+2C_n H_{2n+2}). As ligações covalentes carbono-hidrogênio (CHC-H, energia de ligação 413 kJ/mol\sim 413 \ kJ/mol) e carbono-carbono (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) rompem-se em função do gradiente térmico de contato.

De acordo com os padrões internacionais IEC 60599 e IEEE C57.104, o perfil de gases dissolvidos no óleo (Dissolved Gas Analysis, DGA) permite a identificação termodinâmica precisa da faixa de temperatura da falha no núcleo:

Faixa Térmica de Falha Gás Combustível Dominante Razões de Diagnóstico (IEC 60599 / Duval) Mecanismo Físico no Núcleo Ação Crítica Requerida
Falha Térmica 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 Aquecimento difuso de jugos, degradação inicial de vernizes interlaminares ou tinta de travamento. Monitoramento cromatográfico periódico mensal; verificação de harmônicos na rede.
Falha Térmica 300C<T<700C300^\circ C < T < 700^\circ C (T2) Etileno (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 Espiras em curto-circuito entre múltiplas lâminas de aço CRGO. Degradação severa do revestimento inorgânico. Redução de carga para 70%70\%; ensaio de perdas a vazio e resistência de isolamento do núcleo.
Falha Térmica T>700CT > 700^\circ C (T3) Etileno (C2H4C₂H₄) com traços de Acetileno (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 Fusão localizada do pacote magnético (core burning); curto-circuito entre lâminas e estrutura de aperto. Desconexão imediata forçada. Inspeção interna endoscópica e reconstrução do núcleo.
Descargas Parciais / Centelhamento no Núcleo Hidrogênio (H2H₂) dominante com CH4CH₄ CH4H2<0.1\frac{CH₄}{H₂} < 0.1, H2H₂ excede limites basais (> 100 ppm) Flutuação galvânica do núcleo (perda do aterramento físico único), produzindo descargas capacitivas para o tanque. Verificação do circuito de aterramento externo do núcleo e neutro por meio de megômetro a 2.5 kV2.5 \ kV.

Falha por Múltiplos Aterramentos (Correntes Circulantes em Laço de Terra)

Por diretriz de projeto normativo, o núcleo de um transformador deve ser aterrado galvanicamente em um único ponto físico por meio de uma barra de cobre conectada a uma bucha dedicada na tampa do tanque. Esse aterramento unificado previne a elevação do potencial eletrostático flutuante induzido por acoplamento capacitivo com os enrolamentos de alta tensão (Vind=Cwc/(Cwc+Ccg)VwindingVind = C_{w-c}/(C_{w-c} + C_{c-g}) \cdot Vwinding).

Se, devido à degradação dos isolamentos dos tirantes de fixação (tie-rods), pontes condutoras formadas por sedimentos metálicos no fundo do tanque ou falhas de montagem, for gerada uma segunda conexão à terra inadvertida na extremidade oposta do núcleo, forma-se um circuito fechado de grande área física que enlaça o fluxo magnético principal:

Vloop=ddt aˊrea entre terras BdSVloop = -\frac{d}{dt} \iint_{\text{ área entre terras }} \mathbf{B} \cdot d\mathbf{S}

Essa tensão induz uma corrente de circulação contínua pela estrutura do tanque e pelas lâminas de aterramento do núcleo, variando tipicamente entre 10 A10 \ A e mais de 150 A150 \ A, queimando o isolamento dos parafusos passantes, carbonizando o óleo adjacente e gerando concentrações alarmantes de etileno (C2H4C₂H₄) e metano (CH4CH₄) sem que as proteções elétricas de sobrecorrente ou diferencial (87T87T) detectem qualquer anomalia nas correntes de linha.

Técnicas Avançadas de Mitigação no Projeto e na Fabricação

A otimização do desempenho eletromagnético dos núcleos magnéticos requer uma abordagem multidisciplinar que abrange a modificação cristalográfica do material, a microengenharia de domínios por radiação coerente e a otimização topológica tridimensional das juntas nos jugos.

Juntas de Jugo Tipo Step-Lap a 4545^\circ

A montagem geométrica das colunas e dos jugos em núcleos empilhados introduz descontinuidades no circuito magnético. Nas junções tradicionais de topo com sobreposição simples a 9090^\circ (Butt-Lap joints), o fluxo magnético é forçado a cruzar entre lâminas perpendicularmente à direção de laminação (na direção cristalográfica dura 110\langle 110 \rangle), provocando uma refração severa das linhas de fluxo, um aumento drástico da relutância magnética interfacial e saturações locais violentas nas pontas de corte.

A tecnologia moderna de corte e empilhamento emprega juntas em esquadria a 4545^\circ com defasagem progressiva de lâminas denominada Step-Lap (geralmente configurações de 5 a 7 degraus por ciclo). O escalonamento distribui o entreferro de corte longitudinalmente ao longo de um gradiente axial:

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

Ao defasar as juntas por uma distância espacial Δx3.57.0 mm\Delta x \approx 3.5 - 7.0 \ mm entre lâminas contíguas, o fluxo magnético não é obrigado a saltar perpendicularmente através de um único entreferro de alta relutância para a lâmina adjacente, mas se distribui progressivamente ao longo de uma seção geométrica efetiva equivalente muito maior:

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

A implementação do projeto Step-Lap reduz as perdas totais de junção nos jugos em até 25%30%25\% - 30\%, mitiga a potência reativa de excitação em 40%40\% e reduz os níveis de emissão acústica do transformador em 3a6 dB(A)3 a 6 \ dB(A) em comparação com montagens convencionais.

Refinamento de Domínios Magnéticos por Feixe Laser (Laser Scribing)

Em aços CRGO de alta permeabilidade (Hi-B), o tamanho dos grãos cristalográficos é grande (1030 mm\sim 10 - 30 \ mm), o que favorece a orientação Goss quase perfeita (3\le 3^\circ), mas resulta contraproducente para o comportamento dinâmico de perdas, visto que o espaçamento entre paredes de Bloch de 180180^\circ (2d02d_0) cresce proporcionalmente à raiz quadrada do tamanho de grão:

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

Uma largura de domínio elevada eleva drasticamente a velocidade instantânea de deslocamento da parede vw(t)v_w(t), o que dispara as perdas anômalas de Foucault (PexcPexc).

A técnica de refinamento de domínios por laser (Domain Refinement ou Laser Scribing) projeta um feixe focalizado de laser contínuo (Nd:YAG ou fibra de alta potência) transversalmente à direção de laminação em intervalos espaciais regulares (p48 mmp \approx 4 - 8 \ mm). A energia térmica do pulso laser induz um aquecimento ultrarrápido e localizado sem ablação do material, gerando um microcampo de tensões elásticas de compressão permanente subsuperficiais decorrentes do gradiente térmico de resfriamento:

σ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

Essas faixas periódicas de tensão interna atuam como barreiras de energia magnetoelástica artificiais que fragmentam os domínios magnéticos largos de 180180^\circ, reduzindo seu espaçamento médio 2d02d_0 à metade ou a um terço de sua dimensão original. Como resultado direto:

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 Ptotal \approx -10\% a -15\%

É imperativo distinguir entre o tratamento a laser de tipo não termicamente resistente (o qual perde o efeito de refinamento após recozimentos de alívio de tensões a T>500CT > 500^\circ C exigidos em núcleos enrolados) e os métodos mecânicos de gravação química/deformação plástica controlada com plasma (recozíveis), projetados para reter o ancoramento de microdomínios após tratamentos térmicos a 800C800^\circ C.

Implementação e Modelagem Computacional com Vexten Suite

A quantificação rigorosa das perdas no núcleo e seu acoplamento multifísico com a rede elétrica são resolvidos operacionalmente no Vexten Suite integrando os módulos de fluxo de potência com harmônicos não lineares (conforme IEEE 519 / IEEE C57.110), cálculo térmico dinâmico acoplado (IEC 60076-7) e simulação de curto-circuito / regimes transitórios magnéticos (IEC 60909 / IEEE 141).

Algoritmo de Correção por Distorção Harmônica e Fator K

Na presença de espectros harmônicos severos provenientes de cargas não lineares ou geradores renováveis, o módulo Vexten PowerFlow & Harmonics executa uma partição de perdas dinâmica recalculando o fator de perdas do núcleo de acordo com o tensor de coeficientes espectrais normalizados:

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}

O software vincula esse fator à redução de capacidade nominal (derating) do transformador por meio do cálculo do Fator KK e do Fator Harmônico de Perdas de Foucault (FHLFHL), ajustando a corrente admissível máxima para garantir que a temperatura do ponto mais quente do enrolamento e do núcleo não exceda os limites da classe térmica de 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)}}

Procedimento de Modelagem Multifísica no Vexten Suite

  1. Definição do Espectro de Excitação Magnética: No módulo de qualidade de energia, importar as medições de campo ou espectros típicos IEEE de harmônicos de tensão (VhV_h) e a componente contínua residual (IdcIdc) proveniente do fluxo de potência acoplado com fontes renováveis.
  2. Atribuição de Material Magnético na Base de Dados: Selecionar nas propriedades da máquina o grau específico de aço CRGO (e.g., M3, 23Q090 ou 20Q080-L com refinamento a laser). O Vexten Suite carrega automaticamente os parâmetros microscópicos calibrados: condutividade σ\sigma, coeficiente de histerese khk_h, espessura dd, densidade ρm\rho_m e o parâmetro de Bertotti V0V_0.
  3. Configuração da Topologia de Junta: Selecionar nos parâmetros construtivos o tipo de junta nos jugos: Butt-Lap 9090^\circ, Step-Lap de 5 degraus ou Step-Lap de 7 degraus. O software introduz automaticamente o multiplicador de dispersão interfacial de relutância (Kjoint\mathcal{K}_{joint}).
  4. Execução da Análise Térmica Transitória Acoplada: O solver acopla as perdas eletromagnéticas totais Ptotal(x,y,z)Ptotal(x,y,z) como fonte de calor volumétrica no solucionador de redes térmicas conforme a IEC 60076-7, determinando a distribuição tridimensional de temperaturas nos jugos superiores, inferiores e colunas do transformador.
  5. Validação de Limites e Diagnóstico Forense Preventivo: O sistema compara os gradientes de temperatura local e densidades de perdas com os limites da norma IEEE C57.104. Se a temperatura no núcleo exceder 130C130^\circ C contínuos, ativam-se alertas automáticos de derating de potência e calcula-se a taxa teórica estimada de desprendimento de gases de falha (C2H4C₂H₄, CH4CH₄, H2H₂), permitindo estabelecer planos de manutenção baseada na condição (CBM) de alta precisão física.

Por meio dessa arquitetura computacional abrangente, engenheiros de projeto, comissionamento e confiabilidade podem prever com exatidão o envelhecimento térmico do pacote magnético, eliminar riscos de falhas catastróficas por correntes circulantes interlaminares e otimizar a seleção de materiais magnéticos avançados para transformadores de altíssima eficiência energética.