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
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 -Fe), onde o eixo cristalográfico constitui a direção de mínima energia de anisotropia magnetocristalina, enquanto as direções e 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, , ou nitreto de alumínio, ), desenvolve-se a denominada textura de Goss, caracterizada pela orientação cristalográfica . Nessa configuração, os planos situam-se paralelos à superfície da chapa e as direções orientam-se paralelamente à direção longitudinal de laminação (RD, Rolling Direction).
Onde é a densidade de energia de anisotropia magnetocristalina (), e são as constantes de anisotropia do ferro à temperatura ambiente (), e representam os cossenos diretores do vetor de magnetização 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 e , enquanto em aços de alta permeabilidade (Hi-B) reduz-se para faixas de a .
A adição de silício () em concentrações nominais de a em peso incrementa a resistividade volumétrica elétrica () desde aproximadamente (para ferro puro) até . Esse incremento de resistividade atenua drasticamente o transporte livre de cargas parasitas. No entanto, concentrações superiores a em peso induzem fragilidade mecânica severa devido à formação de fases ordenadas intermetálicas e ( e ), 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 (), minimizando as deformações elásticas induzidas pelo campo magnético e a consequente emissão de ruído acústico.
A estrutura de domínios magnéticos (domínios de Weiss) no CRGO é composta por domínios principais orientados a separados por paredes de Bloch, e domínios de fechamento a nas imediações de descontinuidades superficiais e inclusões. A largura de domínio característica () 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 (), reduzindo a largura média dos domínios a 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 (, em ) é 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:
Onde representa a frequência fundamental de excitação (), é o valor de pico da indução magnética (), é a condutividade elétrica do material (), é a espessura da chapa laminada (), é a densidade de massa volumétrica do aço (, tipicamente ), é o coeficiente de histerese quase-estática, é o expoente de Steinmetz (frequentemente ), e é o parâmetro microscópico de perdas em excesso.
Perdas por Histerese Quase-Estática
As perdas por histerese estática () 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 ou ). Durante o movimento de uma parede de Bloch de , 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:
A área do laço de histerese quase-estático () 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 na presença de formas de onda complexas ou saturação não linear requer a formulação do operador contínuo de Preisach:
Onde é um operador histerético elementar biestável com limites de comutação (de ativação) e (de desativação), e é 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 () 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 na direção do eixo (desde até ), largura infinita no eixo , e comprimento no eixo . Assumindo uma permeabilidade magnética efetiva , uma condutividade volumétrica , e um vetor de densidade de fluxo magnético paralelo à superfície: .
Assumindo que a espessura é substancialmente menor que a profundidade de penetração eletromagnética (skin depth) , a indução magnética distribui-se de maneira espacialmente uniforme através da seção transversal da lâmina: . Integrando em relação a a partir do plano neutro central (, onde por simetria):
A densidade de corrente induzida dissipa uma potência instantânea por efeito Joule cuja densidade volumétrica local é . Integrando sobre a espessura total da chapa e calculando a média no domínio temporal durante um período :
Dividindo pela densidade de massa volumétrica , obtém-se a formulação específica normalizada:
Essa dedução expõe a dependência quadrática direta em relação à espessura da chapa e à frequência fundamental , justificando a necessidade de reduzir os calibres comerciais de (M4) para ou (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 (). 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 . A velocidade de deslocamento de uma parede de domínio 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:
Onde representa o espaçamento entre paredes de domínio e é 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:
Onde é uma constante adimensional associada à dissipação do campo de Foucault em meios infinitos, é a área da seção transversal da lâmina perpendicular à direção do fluxo magnético, e é um parâmetro intrínseco do material com dimensões de campo magnético (), 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 (mm) | Indução Máxima a (T) | Perdas Específicas (W/kg) | Perdas Específicas (W/kg) | Fator de Empilhamento (%) |
|---|---|---|---|---|---|
| M4 (AISI 35G155 / EN 10107) | 0.35 | 1.82 - 1.84 | 1.25 - 1.35 | 1.64 - 1.77 | |
| M3 (AISI 30G130 / EN 10107) | 0.30 | 1.84 - 1.86 | 1.05 - 1.15 | 1.38 - 1.51 | |
| Hi-B Convencional (27Q110) | 0.27 | 1.90 - 1.93 | 0.95 - 1.05 | 1.24 - 1.38 | |
| Hi-B Alta Permeabilidade (23Q090) | 0.23 | 1.92 - 1.95 | 0.80 - 0.90 | 1.05 - 1.18 | |
| Hi-B Laser-Scribed (20Q080-L) | 0.20 | 1.93 - 1.96 | 0.70 - 0.78 | 0.92 - 1.02 | |
| Metal Amorfo (Fe-Si-B Metglas 2605SA1) | 0.025 | 1.56 - 1.58 | 0.18 - 0.25 | 0.24 - 0.33 |
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 (). A resposta não linear da curva de magnetização dinâmica 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 (), a derivada temporal da indução magnética 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.
Sob tensão senoidal pura, . Se a tensão contém distorção harmônica com componentes , o campo induzido harmônico é aproximado por . Substituindo na formulação de correntes clássicas:
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:
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 ().
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 ( em CRGO), a permeabilidade diferencial colapsa drasticamente para valores próximos à permeabilidade do váCuO (). Isso induz picos massivos de corrente magnetizante assimétrica ( 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).
O consumo de potência reativa eleva-se em ordens de magnitude, originando severas quedas de tensão nos terminais, harmônicos pares e ímpares de baixa ordem (especialmente e ) 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 , composto por uma base de forsterita com sobrecamada de fosfato de alumínio e ácido fosfórico) é projetado para suportar tensões interlaminares de alguns poucos volts ( a eficazes conforme a norma IEC 60404-6). Entretanto, sob solicitações térmicas locais contínuas (), fadiga vibracional por magnetoestrição ( 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:
Dado que a resistência do circuito ôhmico é da ordem de miliohms (), a corrente de falha 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 .
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 (). As ligações covalentes carbono-hidrogênio (, energia de ligação ) e carbono-carbono (, ; , ; , ) 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 (T1) | Metano () / Etano () | , , | 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 (T2) | Etileno () | , , | 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 ; ensaio de perdas a vazio e resistência de isolamento do núcleo. |
| Falha Térmica (T3) | Etileno () com traços de Acetileno () | , , | 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 () dominante com | , 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 . |
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 ().
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:
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 e mais de , queimando o isolamento dos parafusos passantes, carbonizando o óleo adjacente e gerando concentrações alarmantes de etileno () e metano () sem que as proteções elétricas de sobrecorrente ou diferencial () 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
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 (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 ), 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 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:
Ao defasar as juntas por uma distância espacial 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:
A implementação do projeto Step-Lap reduz as perdas totais de junção nos jugos em até , mitiga a potência reativa de excitação em e reduz os níveis de emissão acústica do transformador em 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 (), o que favorece a orientação Goss quase perfeita (), mas resulta contraproducente para o comportamento dinâmico de perdas, visto que o espaçamento entre paredes de Bloch de () cresce proporcionalmente à raiz quadrada do tamanho de grão:
Uma largura de domínio elevada eleva drasticamente a velocidade instantânea de deslocamento da parede , o que dispara as perdas anômalas de Foucault ().
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 (). 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:
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 , reduzindo seu espaçamento médio à metade ou a um terço de sua dimensão original. Como resultado direto:
É 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 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 .
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:
O software vincula esse fator à redução de capacidade nominal (derating) do transformador por meio do cálculo do Fator e do Fator Harmônico de Perdas de Foucault (), 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:
Procedimento de Modelagem Multifísica no Vexten Suite
- 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 () e a componente contínua residual () proveniente do fluxo de potência acoplado com fontes renováveis.
- 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 , coeficiente de histerese , espessura , densidade e o parâmetro de Bertotti .
- Configuração da Topologia de Junta: Selecionar nos parâmetros construtivos o tipo de junta nos jugos: Butt-Lap , Step-Lap de 5 degraus ou Step-Lap de 7 degraus. O software introduz automaticamente o multiplicador de dispersão interfacial de relutância ().
- Execução da Análise Térmica Transitória Acoplada: O solver acopla as perdas eletromagnéticas totais 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.
- 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 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 (, , ), 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.