UNIVERSIDADE FEDERAL DE SANTA CATARINA

CAMPUS TRINDADE

Pedro Henrique Cardoso Costa

ESTUDO DE PERDAS POR CORRENTES DE FOUCAULT

EM ÍMÃS DE MOTOR SEM ESCOVAS DE ROTOR INTERNO
E ÍMÃS SUPERFICIAIS

Florianópolis
2023
 PEDRO HENRIQUE CARDOSO COSTA

ESTUDO DE PERDAS POR CORRENTES DE FOUCAULT

EM ÍMÃS DE MOTOR SEM ESCOVAS DE ROTOR
INTERNO E ÍMÃS SUPERFICIAIS

Trabalho de Conclusão de Curso sub-

metido à Universidade Federal de
Santa Catarina, como requisito neces-
sário para obtenção do grau de Bacha-
rel em Engenharia Elétrica

Florianópolis, julho de 2023

 Resumo
O trabalho apresenta uma análise abrangente das perdas por correntes de Foucault em
ímãs e no revestimento do rotor de um motor sem escovas. É realizada uma descrição dos
conceitos-chave relacionados ao funcionamento de motores sem escovas, fornecendo uma
visão geral do tema. Em seguida, uma geometria específica de motor é analisada, e as perdas
nos ímãs e no revestimento do rotor são calculadas por meio do método analítico proposto
no trabalho. Além disso, são realizadas simulações utilizando a técnica de elementos
finitos, tanto em uma como em duas dimensões, para obter resultados comparativos. A
análise desses resultados apresenta detalhes importantes do comportamento das perdas
nessa topologia de máquina, bem como as limitações da abordagem de cálculo de perdas
utilizada.

Palavras-chave: Perdas por correntes de Foucault; Motores sem escovas; Ímãs perma-
nentes; Cálculo analítico; Elementos finitos.
 Lista de ilustrações

Figura 1 – Comparação entre diferentes configurações de rotor . . . . . . . . . . . 18

Figura 2 – Comparação entre diferentes configurações de ímãs no rotor . . . . . . 19
Figura 3 – Visão do revestimento do rotor . . . . . . . . . . . . . . . . . . . . . . 20
Figura 4 – Ilustração do segundo quadrante do laço de histerese de um ímã perma-
nente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Figura 5 – Ilustração da modulação do fluxo no entreferro (Bg ) pela abertura das
ranhuras do estator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
Figura 6 – Visão do motor em corte . . . . . . . . . . . . . . . . . . . . . . . . . . 28
Figura 7 – Divisão de regiões para o cálculo de perdas . . . . . . . . . . . . . . . . 31
Figura 8 – Interface do SimcenterT M SPEED . . . . . . . . . . . . . . . . . . . . . 35
Figura 9 – Interface do SimcenterT M MAGNET . . . . . . . . . . . . . . . . . . . 36
Figura 10 – Visão frontal e axial da máquina . . . . . . . . . . . . . . . . . . . . . 37
Figura 11 – Malhas utilizadas para as simulações em elementos finitos em t = 0 ms 38
Figura 12 – Mapa de perdas médias no rotor resultante do MEF3D . . . . . . . . . 42
Figura 13 – Perdas instantâneas em um dos ímãs - MEF2D . . . . . . . . . . . . . 43
Figura 14 – Perdas instantâneas no revestimento - MEF2D . . . . . . . . . . . . . . 43
Figura 15 – Ilustração das dimensões da máquina . . . . . . . . . . . . . . . . . . . 51
 Lista de tabelas

Tabela 1 –
Dimensões da máquina . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Tabela 2 –
Dados da ferrite NeIGT . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Tabela 3 –
Dados do aço 304 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Tabela 4 –
Parâmetros da simulação transiente . . . . . . . . . . . . . . . . . . . . 39
Tabela 5 –
Valor médio das perdas nos ímãs por correntes induzidas para diferentes
valores de resistividade . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
Tabela 6 – Valor médio das perdas no tempo por correntes induzidas . . . . . . . 41
Tabela 7 – Comparação entre os valores de perdas nos ímãs com e sem o revestimento 42
Tabela 8 – Diferença percentual entre métodos de cálculo das perdas por correntes
induzidas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
Lista de abreviaturas e siglas

CC = Corrente contínua

CA = Corrente alternada

RMS = Root-mean-square

MEF = Método de elemento finitos

BLDC = Brushless DC
 Sumário

1 INTRODUÇÃO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.1 Problemática . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2 Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.2.1 Objetivo Geral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.2.2 Objetivos Específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.3 Organização . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

2 MÁQUINA SEM ESCOVAS DE ROTOR INTERNO COM ÍMÃS PER-

MANENTES SUPERFICIAIS . . . . . . . . . . . . . . . . . . . . . . . 17
2.1 Definição . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
[Link] Máquina sem escovas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
[Link] Rotor interno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
[Link] Ímãs superficiais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
[Link] Revestimento do rotor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2 Princípio de funcionamento . . . . . . . . . . . . . . . . . . . . . . . 20
2.3 Características dos ímãs permanentes . . . . . . . . . . . . . . . . 20

3 PERDAS EM MÁQUINAS ELÉTRICAS . . . . . . . . . . . . . . . . . 23

3.1 Perdas por correntes induzidas no rotor de máquinas síncronas . 24
3.1.1 Natureza das correntes induzidas nos ímãs e no revestimento do rotor 26
3.2 Modelo do rotor com 4 regiões . . . . . . . . . . . . . . . . . . . . . 27
3.3 Solução da equação de difusão complexa . . . . . . . . . . . . . . 28
3.3.1 Aplicação nas regiões consideradas . . . . . . . . . . . . . . . . . . . 30

4 MODELAGEM E SIMULAÇÃO . . . . . . . . . . . . . . . . . . . . . 35
4.1 Softwares de simulação . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2 Máquina considerada . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.3 Simulação em Elementos Finitos: Configurações e Parâmetros . 38

5 RESULTADOS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
5.1 Natureza das correntes induzidas . . . . . . . . . . . . . . . . . . . 41
5.2 Perdas calculadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
5.2.1 Distribuição das perdas . . . . . . . . . . . . . . . . . . . . . . . . . . 42
[Link] Perdas instantâneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
5.2.2 Comparação entre os métodos de cálculo . . . . . . . . . . . . . . . . 44
5.3 Estratégias de mitigação . . . . . . . . . . . . . . . . . . . . . . . . . 44
6 CONCLUSÃO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

REFERÊNCIAS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49

APÊNDICE A – INDICAÇÃO DAS DIMENSÕES DA MÁQUINA . . 51

 15

1 Introdução

A eficiência energética é uma questão de grande importância no contexto atual de

crescente preocupação ambiental e busca por soluções sustentáveis [Paris Agreement]. Além
disso, há grande pressão para o aumento da relação custo-benefício devido às exigências
cada vez mais rigorosas por parte de órgãos governamentais. Em particular, os motores
sem escovas têm ganhado popularidade devido à sua alta eficiência, seu design compacto e
seu baixo consumo de energia quando comparados a motores de indução [Han et al. 2018].
Isso os torna amplamente utilizados em diversas aplicações, desde veículos elétricos até
aparelhos domésticos. No entanto, para garantir o desempenho otimizado desses motores,
é crucial compreender e minimizar suas perdas.
O cálculo preciso das perdas por corrente induzidas desempenha um papel funda-
mental no projeto térmico desses motores [Umans 2014]. À medida que a corrente elétrica
circula pelos ímãs, são geradas perdas que se manifestam na forma de calor. A dissipação
eficiente desse calor é essencial para evitar o superaquecimento dos ímãs, o que pode levar
à desmagnetização irreversível e falha fatal na operação. Além disso, altas temperaturas
no rotor podem afetar a lubrificação e confiabilidade dos rolamentos do sistema e limitar a
velocidade rotacional [Han et al. 2018]. Devido à complexidade da distribuição do fluxo
magnético e às não idealidades geométricas, calcular e minimizar essas perdas representam
desafios significativos [Hendershot e Miller 2010].
Nesse contexto, este trabalho de conclusão de curso investiga as perdas por correntes
induzidas nos ímãs do rotor de um motor sem escovas de ímãs permanentes superficiais,
além de analisar a precisão de aproximações analíticas em comparação com métodos
numéricos.
Durante o restante desse trabalho, a configuração de motor utilizada será referenci-
ada como BLDC (do inglês Brushless DC ).

1.1 Problemática

Perdas por correntes de Foucault na região dos ímãs e do revestimento do rotor em

máquinas BLDC são problemáticas, visto a dificuldade em remover o calor resultante em
regiões de difícil acesso do ponto de vista fluidodinâmico [Han et al. 2018]. Além disso, o
cálculo dessas correntes é matematicamente complexo, sendo mais difícil, em geral, validar
esses cálculos do que no caso de modelos magnetoestáticos [Siemens 2022].
16 Capítulo 1. Introdução

1.2 Objetivo
1.2.1 Objetivo Geral
O objetivo principal desse Trabalho de Conclusão de Curso é analisar o compor-
tamento das perdas por corrente de Foucault nos ímãs e no revestimento de um motor
BLDC mediante um método analítico e pelo cálculo por elementos finitos.

1.2.2 Objetivos Específicos

• Simular a topologia de motor escolhida através do MEF;

• Comparar as aproximações analíticas às soluções por MEF;

• Determinar o impacto do revestimento condutor sobre o valor das perdas nos ímãs;

• Revisar a teoria acerca da modelagem de perdas em motores.

1.3 Organização
O presente trabalho está organizado da seguinte forma:
Capítulo 2: os conceitos relacionados às características do tipo de motor abordado
são apresentados, bem como seu princípio de funcionamento e algumas das características-
chave analisadas posteriormente.
Capítulo 3: as principais perdas relacionadas à máquina elétrica estudada são apre-
sentadas, com enfoque nas perdas por correntes induzidas na região dos ímãs permanentes
e do revestimento do rotor; é apresentado também um método para o cálculo dessas perdas
a partir de uma formulação em vetor potencial magnético.
Capítulo 4: descreve-se a máquina utilizada nesse trabalho, sua geometria, alimenta-
ção, e as condições determinadas para a simulação em elementos finitos. São apresentados
também os softwares utilizados para o cálculo das perdas.
Capítulo 5: os resultados dos cálculos, as análises de dados e as simulações realizadas
são apresentados neste capítulo.
 17

2 Máquina sem escovas de rotor interno

com ímãs permanentes superficiais

2.1 Definição
[Link] Máquina sem escovas

Como o nome sugere, esse tipo de máquina se caracteriza pela ausência de conexões
deslizantes associadas à alimentação do rotor, como escovas, comutadores ou anéis coletores.
Ao estabelecer o campo do rotor com ímãs permanentes, elimina-se a necessidade de
alimentação do rotor.
A ausência dessas conexões implica em vantagens e desvantagens em relação às
demais configurações [Umans 2014]. Em primeira análise, percebe-se que o uso de ímãs
permanentes elimina as perdas no cobre associadas à energização do rotor, bem como as
perdas mecânicas nos comutadores; essa é uma grande vantagem para situações em que é
difícil dissipar o calor gerado, como o caso de rotores internos. Além disso, cerca de 90%
da manutenção rotineira em motores CC convencionais ocorre nos comutadores [Gieras e
Wing 1996].
Essas características tornam o motor sem escovas atrativo para aplicações como
discos rígidos, em que os detritos das escovas são indesejáveis; em sopradores e ventiladores,
onde se buscam baixo ruído e controle de velocidade; e em veículos, visto o custo de
manutenção das escovas e a busca por eficiência.
Entretanto, a substituição de enrolamentos por ímãs permanentes implica na perda
de controle do campo magnético gerado no rotor. Assim, o desenvolvimento desse tipo
de motor se popularizou com avanços referentes à tecnologia de ímãs permanentes e
áreas-chave de controle, como transistores de potência, sensores e microeletrônica. Uma vez
que o campo no rotor é relativamente constante, técnicas mais avançadas de alimentação
do estator são necessárias para compensar a perda de controle no campo do rotor; é preciso
inverter a polaridade da alimentação toda vez que ocorre o alinhamento com o polo de um
dos ímãs para garantir torque unidirecional. E mesmo com a melhoria nas características
dos ímãs nas últimas décadas, o custo proibitivo de ímãs permanentes continua sendo um
dos fatores limitantes no projeto desse tipo de máquina [Barcelos 2021].

[Link] Rotor interno

Uma máquina sem escovas pode possuir diferentes topologias. A máquina com
rotor interno é a mais parecida com uma máquina CA síncrona, na qual o rotor com ímãs é
18 Capítulo 2. Máquina sem escovas de rotor interno com ímãs permanentes superficiais

circundado pelo estator. Isso permite uma alta relação de torque/inércia e certa facilidade
na refrigeração, mas apresenta problemas quanto à fixação dos ímãs no rotor; por vezes é
necessário o uso de um revestimento envolvendo os ímãs para evitar que escapem devido à
força centrífuga em altas velocidades.
Outras configurações são possíveis, como um estator central circundado pelo rotor,
que possui os ímãs em sua superfície interna. O rotor externo permite a proteção da carcaça
do rotor contra o desprendimento dos ímãs e, em geral, possui maior raio de entreferro.
Com isto, possibilita maior capacidade de atingir torque elevado [Barcelos 2021]. Por outro
lado, ele apresenta maiores complicações para a refrigeração dos enrolamentos.
Uma configuração comum é um disco com um ímã permanente em forma de arruela
fixado em um dos lados formando o rotor, enquanto o estator tem enrolamentos geralmente
elaborados em circuito impresso. Isso configura uma máquina de rotor axial; ele possui uma
espessura pequena, operação suave e um custo baixo, mas não possui um bom desempenho
para velocidades e torques elevados [Hendershot e Miller 2010].

Figura 1 – Comparação entre diferentes configurações de rotor

Fonte: adaptado de [Barcelos 2021]

2.1. Definição 19

[Link] Ímãs superficiais

Esse trabalho apresenta enfoque na análise de uma máquina com rotor interno.
Nesta configuração é possível encontrar rotores com ímãs fixados na superfície ou internos
ao corpo do rotor.

Figura 2 – Comparação entre diferentes configurações de ímãs no rotor

Fonte: adaptado de [Hendershot e Miller 2010]

Em motores com seis ou mais polos, é possível alcançar altos níveis de fluxo
magnético no entreferro com ímãs de ferrite internos, o que, com rotores de ímãs superficiais,
só seria possível com materiais de terras raras [Hendershot e Miller 2010]. Outro ponto
relevante é maior resistência dos ímãs internos ao desprendimento em altas velocidades
nessa configuração. Entretanto, essa geometria causa uma saliência entre o eixo direto e
o de quadratura, que pode ser indesejada. Além disso, devido à ausência de aberturas,
a capacidade de troca de calor com o ar reduz a densidade de corrente permitida nos
enrolamentos.

[Link] Revestimento do rotor

Em condições de velocidade elevada, as máquinas com ímãs superficiais e rotores

internos precisam tomar medidas quanto ao desprendimento dos ímãs, sendo muitas vezes
empregado um revestimento no rotor, permitindo a fixação mecânica [Han et al. 2018].
Esse revestimento pode ser feito de material condutor, ocasionando perdas por correntes
induzidas. Uma visão do rotor, dos ímãs e do revestimento é apresentada na figura 3.
Embora materiais não-condutores sejam por vezes considerados para a construção
do revestimento, eles possuem baixa condutividade térmica e dificultam a dissipação de
20 Capítulo 2. Máquina sem escovas de rotor interno com ímãs permanentes superficiais

Figura 3 – Visão do revestimento do rotor

Fonte: autoria própria

calor no rotor, resultando em um aumento de temperatura. Além disso, o coeficiente

de expansão de revestimentos não-condutores é, em geral, menor do que o da corroa
do rotor e dos ímãs. Isso leva a um grande estresse mecânico sobre a estrutura. Esses
materiais também costumam ser feitos de materiais compósitos que complicam os aspectos
construtivos do rotor. Essas características fazem com que materiais não-condutores sejam
uma escolha inadequada para a maioria das aplicações.

2.2 Princípio de funcionamento

Como no motor sem escovas não há contato entre a alimentação do estator e o rotor,
a transferência de energia se dá completamente através da conversão da energia elétrica do
estator em energia magnética no entreferro, que então é transformada em energia mecânica
para movimentar o rotor através da força magnética proveniente da interação com o campo
dos ímãs permanentes. Em um cenário ideal, toda a energia elétrica fornecida ao motor
seria convertida em energia mecânica. No entanto, em dispositivos reais, o processo de
conversão eletromagnética de energia inevitavelmente resulta em perdas no núcleo de aço,
bem como perdas no cobre e perdas mecânicas [Umans 2014].

2.3 Características dos ímãs permanentes

Os ímãs permanentes são objetos que possuem propriedades magnéticas duradouras,
ou seja, mantêm um campo magnético mesmo na ausência de uma fonte externa de energia.
Eles são materiais chamados magnéticos duros e suas características resultam da orientação
ordenada nos movimentos dos elétrons ao redor de seus núcleos.
Em um nível microscópico existem pequenas regiões no material onde os átomos
estão alinhados magneticamente, chamadas Domínios de Weiss, e separadas entre si por
2.3. Características dos ímãs permanentes 21

Paredes de Bloch, geralmente posicionadas nas imperfeições e impurezas do material [Jiles

e Atherton 1986]. Cada domínio possui um campo magnético próprio, mas naturalmente
eles estão orientados aleatoriamente, resultando em um campo magnético líquido igual a
zero para o material como um todo.
Para adquirirem direções definidas e formarem os domínios, é necessário submeter
o material à ação de um campo magnético externo suficientemente forte, em um processo
conhecido como magnetização, que consiste na expansão dos domínios com direção igual
à do campo externo e encolhimento daqueles com direção contrária. Após este processo,
mesmo que se retire o campo magnético externo, o material apresentará uma magnetização
residual significativa devido aos momentos magnéticos estarem alinhados em uma direção
bem definida [Barcelos 2021].
Diferentes materiais magnéticos podem apresentar características distintas. Alguns
materiais, como o ferro, o cobalto e o níquel, são naturalmente magnéticos e podem ser
permanentemente magnetizados. Ligas de terras raras também são amplamente utilizados
em ímãs permanentes de alto desempenho. As principais características de comparação
entre os materiais são:

• Indução remanente (Br ): indução magnética do ímã magnetizado, quando há ausência

de um campo externo;

• Coercitividade (Hcj ): força magnetizante oposta necessária para eliminar a remanên-

cia;

• Coercitividade intrínseca (Hci ): força magnetizante que reduz a polarização intrínseca

a zero;

• Permeabilidade relativa (µr ): gradiente da curva B/H no ponto de remanência relativo

à permeabilidade magnética no vácuo (µo );

• Magnetização no joelho da curva (Hk ): valor do campo magnético no segundo

quadrante abaixo do qual o ímã sofre desmagnetização irreversível;

• Produto BH máximo (BHmax): medida da sua energia magnética máxima armaze-

nada.

Essas características são obtidas no segundo quadrante da curva de histerese de um ímã

permanente, e indicadas na figura 4.
22 Capítulo 2. Máquina sem escovas de rotor interno com ímãs permanentes superficiais

Figura 4 – Ilustração do segundo quadrante do laço de histerese de um ímã permanente

Fonte: autoria própria

 23

3 Perdas em máquinas elétricas

No projeto de máquinas elétricas, um aspecto de extrema importância é o dimensio-

namento adequado das perdas. As perdas desempenham um papel importante na eficiência
energética, no desempenho térmico e na vida útil desses dispositivos; o dimensionamento
preciso das perdas permite um projeto eficiente e confiável.
A seguir serão discutidas a natureza das perdas no rotor da máquina analisada, as
aproximações e a divisão de regiões considerada para o cálculo das perdas e a resolução
do problema através da solução da equação de difusão complexa aplicada para estruturas
cilíndricas de várias camadas [Miller e Lawrenson 1976]. Nesse trabalho são abordados
mais detalhadamente os mecanismos de ocorrência das perdas por correntes induzidas
nos ímãs e no revestimento do rotor, suas implicações no desempenho das máquinas e
possíveis estratégias utilizadas para minimizar essas perdas.
As perdas em máquinas elétricas podem ocorrer em várias formas, sendo os princi-
pais tipos de perdas as perdas no cobre, perdas mecânicas e perdas no ferro. As perdas no
cobre estão relacionadas à resistência elétrica dos enrolamentos e a correntes que passam
por eles. Já as perdas no ferro estão associadas ao processo de magnetização e desmagneti-
zação do material magnético utilizado nas partes fixas da máquina. As perdas mecânicas,
por sua vez, estão relacionadas ao atrito e às vibrações decorrentes do movimento mecânico
da máquina.

Perdas no cobre
As perdas no cobre referem-se à dissipação de energia devido à resistência elétrica
do condutor, ocorrendo principalmente nas bobinas do enrolamento do estator e do rotor.
Na topologia de motor estudada não há preocupação com as perdas no cobre do rotor, pois
o campo magnético é mantido pelos ímãs permanentes e não há alimentação de correntes.
Existem várias estratégias para mitigar as perdas no cobre. Uma delas é o uso de
condutores de maior condutividade, como cobre puro ou ligas condutoras especiais. Isso
reduz a resistência do condutor e, consequentemente, as perdas no cobre. Além disso, a
otimização do design do motor, como a seleção adequada do tamanho e da geometria das
bobinas, pode minimizar as perdas no cobre. O uso de técnicas avançadas de resfriamento,
como a utilização de líquidos refrigerantes ou ventilação forçada, também pode ajudar a
reduzir as perdas no cobre, melhorando a dissipação térmica nas bobinas.
Em um motor com ímãs permanentes elimina-se uma parte dos condutores utilizados,
e assim as perdas no cobre. Essa característica pode ser especialmente atrativa em contextos
de baixa dissipação do calor gerado.
24 Capítulo 3. Perdas em máquinas elétricas

Perdas mecânicas
As perdas mecânicas em máquinas elétricas referem-se à dissipação de energia
devido a atritos e resistências mecânicas durante o funcionamento do motor. Essas perdas
são influenciadas por vários fatores, como a fricção entre os componentes rotativos, a
resistência do rolamento e a perda por ventilação. Um dos atrativos de máquinas sem
escovas é a diminuição de componentes em atrito durante a operação do sistema.

Perdas no ferro
As perdas no ferro são uma forma de perda energética que ocorre devido às
propriedades magnéticas do material. Elas são causadas principalmente por dois fenômenos:
histerese magnética e correntes parasitas.
A histerese magnética é um comportamento intrínseco do material ferromagnético
quando submetido a variações de campo magnético. Esse processo resulta em uma perda
de energia devido à resistência do material e à reversão magnética a medida que o campo
aplicado percorre o ciclo de histerese. As correntes parasitas, também conhecidas como
correntes de Foucault, são correntes induzidas no núcleo de ferro pela variação do campo
magnético. Essas correntes circulam no material ferroso e causam perdas adicionais devido
à resistência elétrica do núcleo.
Para mitigar as perdas no ferro, são utilizadas técnicas como a seleção de materiais
magnéticos com propriedades otimizadas, o design adequado do núcleo da máquina e
a aplicação de revestimentos isolantes para reduzir as correntes parasitas. Além disso,
simulações computacionais e análises de elementos finitos são empregadas para prever e
otimizar o comportamento magnético e as perdas no ferro durante o projeto.

3.1 Perdas por correntes induzidas no rotor de máquinas

síncronas
Em uma máquina ideal, o campo gira em sincronismo com o rotor e a densidade de
fluxo é invariante no tempo em toda a seção transversal do rotor, visto que o movimento
relativo é nulo. Não há tendência para correntes parasitas fluírem em qualquer lugar do
rotor. Essas condições ideais consideram velocidade constante; uma máquina sem ranhuras;
enrolamentos com distribuição senoidal; e correntes senoidais balanceadas. Na prática, as
correntes parasitas são induzidas por imperfeições ou desvios da máquina síncrona ideal.
Harmônicos de espaço devido à distribuição dos condutores do estator e harmônicos
de tempo na forma de onda da corrente do estator produzem componentes de campo
assíncronos que giram em relação ao rotor. Em geral, os harmônicos de espaço e tempo se
combinam para produzir um número potencialmente grande de componentes assíncronos.
3.1. Perdas por correntes induzidas no rotor de máquinas síncronas 25

Por exemplo, a troca das correntes de uma fase para outra se comporta como um degrau de
tensão nos terminais do estator, induzindo uma força eletromotriz no rotor. As correntes
parasitas resultantes são transitórias, mas excitadas regularmente, resultando em uma
perda média em regime permanente.
Além das possíveis harmônicas nas formas de onda das correntes, qualquer desequi-
líbrio entre as fases pode produzir um componente de sequência negativa do campo de
reação da armadura, mesmo quando as correntes de fase são senoidais. Essa componente
de sequência negativa gira na direção contrária ao rotor e induz uma força eletromotriz
com o dobro da frequência fundamental.
Outra fonte de harmônicas decorre da frequência de passagem das ranhuras do
estator devido ao giro do rotor, que pode modular a relutância geral do circuito magnético.
Esses harmônicos de relutância produzem torque de dentação (do inglês cogging torque) e
geram correntes parasitas no rotor.
As aberturas das ranhuras do estator modulam também a distribuição do fluxo do
entreferro, criando “vales” que giram para trás na velocidade síncrona em relação ao rotor.
Esses picos produzem conjuntos de correntes parasitas que permanecem estacionárias no
espaço, próximas às “sombras” no rotor das aberturas das ranhuras. Esse mecanismo de
perdas é representado na figura 5.

Figura 5 – Ilustração da modulação do fluxo no entreferro (Bg ) pela abertura das ranhuras
do estator

Fonte: autoria própria

26 Capítulo 3. Perdas em máquinas elétricas

3.1.1 Natureza das correntes induzidas nos ímãs e no revestimento do

rotor
As correntes induzidas nos ímãs e em um revestimento condutor são normalmente
causadas pela variação de fluxo imposta sobre esses componentes. O processo é o mesmo
que a geração de correntes induzidas clássica, como descrito em [Bastos 2004] ou [Stoll
1974].
Comumente, se analisa o comportamento das correntes induzidas a partir da relação
entre as grandezas da geometria e a penetração do campo representada pela profundidade
pelicular (δ). Classicamente, a profundidade pelicular é inversamente proporcional à raiz
da frequência angular do campo (ω), à raiz da permeabilidade magnética do material (µ)
e à raiz da condutividade elétrica (σ):
r
2
δ=
ωµσ
Isso leva à identificação de duas situações: o caso em que as grandezas avaliadas são
comparáveis ou pequenas em relação à profundidade pelicular e o caso em que são maiores.
No primeiro caso, pode-se dizer que o campo variante penetra toda a extensão da
geometria analisada e, normalmente, tem natureza de baixa frequência. Nessa situação, as
correntes são limitadas pela resistividade do material, dado que as perdas são calculadas a
partir de
E
W =
R
em que E é o campo elétrico, proporcional ao quadrado da frequência, e R é uma resistência
constante, dependente somente do material. Nessa situação, as correntes induzidas não
possuem impacto no campo que as gerou e diminuem conforme o aumento da resistividade
do material.
Conforme a frequência do campo aumenta, o efeito pelicular toma precedência no
comportamento do fenômeno, e as correntes induzidas se concentram em áreas menores
na superfície da geometria. As correntes de Foucault aumentam em magnitude e há
uma mudança em suas fases, de forma que em altas frequências formam uma espécie de
blindagem que impede a penetração do campo no material condutor.
Ao mesmo tempo, o aumento da frequência aumenta a resistência do caminho
das correntes induzidas, devido à diminuição da profundidade pelicular e consequente
diminuição da área efetiva para circulação das correntes. Isto é, a resistência é proporcional
ao inverso da profundidade pelicular. Podemos analisar as perdas a partir de

W = J2 R
√
em que J se mantém constante e R é proporcional a ω. Nesse contexto, o comportamento
das perdas assume uma natureza diferente: as correntes induzidas produzem um campo
3.2. Modelo do rotor com 4 regiões 27

oposto ao campo original e as perdas aumentam conforme o aumento da resistividade do

material. Nesse caso, pode-se dizer que as correntes são limitadas pela indutância.
Na maioria das situações um aumento de resistividade significa uma diminuição
das perdas por correntes induzidas, mas esse nem sempre é o caso. Em situações nas quais
o componente age como uma blindagem eletromagnética, a baixa resistividade tenderia
a tornar o comportamento ainda mais como o de uma blindagem perfeita. Visto como
esse comportamento pode afetar o projeto magnético da máquina de várias maneiras, é
importante analisar a natureza dessas correntes para se situar corretamente na escolha de
materiais e geometrias.
A análise mostrada é uma simplificação conveniente, embora possa ser válida
em vários contextos; o comportamento das correntes induzidas é complexo e envolve
mais variáveis que a profundidade pelicular ou uma grandeza de referência na direção de
penetração do campo [Stoll e Hammond 1965]. Por exemplo, a natureza do comportamento
pode depender também do espaçamento entre os polos ou da largura da geometria medida
em um plano normal à direção de fluxo.

3.2 Modelo do rotor com 4 regiões

O estudo dessas perdas de uma forma analítica é facilitada pela simplificação do
sistema em um problema bi-dimensional. Esse sistema consiste em um plano que atravessa
a máquina na direção perpendicular ao seu eixo e considera quatro regiões: o entreferro, o
revestimento do rotor, os ímãs e a coroa do rotor. O corte do motor considerado pode ser
visto na figura 6, com a demarcação dos raios para cada região.
Na imagem estão indicadas as diferentes regiões do motor, em que rh é o raio
externo da coroa do estator; r1 é o raio externo dos ímãs; r2 é o raio externo do revestimento
do rotor e rs é o raio externo do entreferro.
As perdas nessa geometria podem ser analisadas através da solução do campo
eletromagnético gerado pelas harmônicas temporais e espaciais. Assumindo que essas
harmônicas não afetam as correntes no estator, de modo a separar as análises de circuito
e de campo, e considerando também cada região como um cilindro contínuo, pode-se
desenvolver a solução através da equação de difusão complexa aplicada a estruturas
cilíndricas de várias camadas [Miller e Lawrenson 1976].
28 Capítulo 3. Perdas em máquinas elétricas

Figura 6 – Visão do motor em corte

Fonte: autoria própria

3.3 Solução da equação de difusão complexa

A partir da formulação do potencial vetorial magnético [Bastos 2004]

∇×A=B (3.1)

aplica-se o rotacional em ambos os lados

∇×∇×A=∇×B

Sendo J a densidade de corrente

∇×H=J

e considerando uma relação linear de permeabilidade B = µH, tem-se na equação anterior

∇ × ∇ × A = µJ

No contexto de correntes induzidas, podemos escrever J em função do potencial vetorial

magnético. A partir da lei de Faraday em sua forma vetorial
∂B
∇×E=−
∂t
3.3. Solução da equação de difusão complexa 29

substitui-se o potencial vetorial magnético

∂∇ × A
∇×E+ =0
 ∂t 
∂A
∇× E+ =0
∂t

a partir do que pode-se definir um potencial escalar Φ:

∂A
E=− − ∇Φ
∂t
pois ∇ × ∇Φ é nulo. Sendo o gradiente desse potencial escalar nulo no contexto das
correntes induzidas em um material condutor, vale

∂A
J = σE = −σ (3.2)
∂t
e assim
∂A
∇ × ∇ × A = µσ
∂t
A identidade vetorial fornece [Bastos 2004]:

∇ × ∇ × A = ∇(∇ · A) − ∇2 A = −∇2 A

dado a consideração de que ∇ · A = 0, o que é normalmente assumido para esse tipo de

análise. Portanto
∂A
∇2 A = µσ
∂t
Dada a geometria do rotor, é adequado o tratamento do problema em coordenadas
cilíndricas com origem no centro do rotor. Aproximando o comportamento de A por uma
única componente na direção z, a equação é simplificada para o que é conhecido como
equação da difusão em duas dimensões.

∂ 2 A 1 ∂A 1 ∂ 2A ∂A
2
+ + 2 2
= µσ
∂r r ∂r r ∂θ ∂t
Uma solução é encontrada através da separação de variáveis, em que um campo rotativo
leva a
A(r, θ, t) = A(r)ej(ωt−pθ)

onde ω é a frequência angular considerada e p é o número de par de polos magnéticos do

rotor. Substituindo a solução na equação da difusão, tem-se

d2 A 1 dA p2 A
+ − 2 = jwµσA
dr2 r dr r

d2 A 1 dA p2
 
+ − + jωµσ A = 0 (3.3)
dr2 r dr r2
30 Capítulo 3. Perdas em máquinas elétricas

Essa é uma variação da equação de Bessel, com solução dada por [Hendershot e Miller
2010]
A = c1 Ip (j 1/2 r/d) + c2 Kp (j 1/2 r/d) (3.4)
1
em que d = √ωµσ , ou ainda uma fração da profundidade pelicular, d = √12 δ; Ip e Kp são
funções de Bessel modificadas do primeiro e segundo tipo, respectivamente [Hendershot e
Miller 2010]; c1 e c2 são constantes complexas.
Em regiões não condutoras σ = 0, e a equação 3.3 se reduz a uma equação de Laplace

d2 A 1 dA p2
+ − 2A = 0
dr2 r dr r
com solução [Bastos 2004]
A = c3 rp + c4 r−p (3.5)

Além da solução para o potencial vetorial magnético A, é preciso relacionar as

soluções aos campos usados em análises práticas, isso é, B, H, E e J. É importante notar
que as condições de contorno são normalmente descritas em termos de B e H. Pela
definição do potencial vetorial magnético em 3.1

1 ∂A ∂A
Br = ; Bθ = µo Hθ = − (3.6)
r ∂θ ∂r
A relação com o campo elétrico é dada por

∂A
E=− = −jωA
∂t
e, em regiões condutoras, a corrente induzida é 3.2, que pode ser escrita como

J = −jωσA

Isso posto, as perdas na região cilíndrica podem ser calculadas a partir da integração
2
de J /σ no volume. Desenvolvendo essa integração, tem-se a densidade de potência por
unidade de comprimento ao longo do eixo em um cilindro contínuo com raios r1 e r2 .
Z r2
2π 2
W = ω σ AA∗ rdr W/mz (3.7)
2 r1

3.3.1 Aplicação nas regiões consideradas

A máquina é dividida em 4 regiões, como delimitado na seção 3.2. A divisão das
regiões é novamente representada na figura 7.
Para a região do entreferro, a solução do potencial vetorial magnético se comporta
como visto para regiões não condutoras, em 3.5. A partir das relações apresentadas em 3.6
3.3. Solução da equação de difusão complexa 31

Figura 7 – Divisão de regiões para o cálculo de perdas

Fonte: autoria própria

deduzem-se os demais campos.

A = c1 rp + c2 r−p
p
Br = −j (c1 rp + c2 r−p )
r
p
Hθ = − (c1 rp−1 + c2 r−p−1 )
µo

Considerando a região do revestimento como condutora, vale a solução obtida em 3.4.

A = c3 Ip (z) + c4 Kp (z)
p
Br = −j (c3 Ip (z) + c4 Kp (z))
r
j 1/2
Hθ = − (c3 Ip 0 (z) + c4 Kp 0 (z))
µR dR

sendo z = j 1/2 drR ; e dR = √ωµ1r σr = √δs2 ; µr , σr e δs a permeabilidade magnética em H/m, a

condutividade em S/m e a profundidade pelicular do revestimento, respectivamente.
Considerando os ímãs também como condutores

A = c5 Ip (z) + c6 Kp (z)
p
Br = −j (c5 Ip (z) + c6 Kp (z))
r
j 1/2
Hθ = − (c5 Ip 0 (z) + c6 Kp 0 (z))
µi di
32 Capítulo 3. Perdas em máquinas elétricas

sendo z = j 1/2 dri , e di = √ωµ1 i σi = √δi2 ; µi , σi e δi a permeabilidade magnética em H/m, a

condutividade em S/m e a profundidade pelicular dos ímãs, respectivamente.
Por fim, sendo a região da coroa do rotor condutora,

A = c7 Ip (z)
p
Br = −j c7 Ip (z)
r
j 1/2
Hθ = − c7 Ip 0 (z)
µe de

sendo z = j 1/2 dre , de = √ωµ1e σe = √δe2 e µe , σe e δe a permeabilidade magnética em H/m, a

condutividade em S/m e a profundidade pelicular da coroa do rotor, respectivamente. O
termo Kp (z) não está presente, pois a solução deve ser finita em r = 0, e Kp tende ao
infinito quando r tende a zero.
Por fim, os coeficientes complexos ck devem ser resolvidos através da solução do
sistema algébrico obtido a partir das condições de contorno do problema.

Condições de contorno

Uma simulação completa da máquina e seu acionamento é necessária para obter

as harmônicas no tempo das formas de onda de corrente e uma análise harmônica do
enrolamento é necessária para obter as harmônicas espaciais da distribuição de ampere-
condutor, mas isso foge do escopo desse trabalho, sendo considerado um dado para a
análise.
Em r = rs tem-se um material de permeabilidade infinita em que se considera
a presença de uma folha de corrente K. Isso significa que a região do estator é uma
superfície em que o campo tangencial é definido pela folha de corrente representando todos
os ampere-condutores do estator. Com isso, a condição de contorno se dá por:

c1 rsp + c2 rs−p = K

Em r = r2 há a interface entre o entreferro e o revestimento, em que Hθ e Br são contínuos.

Dessa forma
p j 1/2
− (c1 r2p−1 + c2 r2−p−1 ) = − (c3 Ip 0 (z2R ) + c4 Kp 0 (z2R ))
µo µR dR
e
p p
(c1 r2p + c2 r2−p ) = −j (c3 Ip (z2R ) + c4 Kp (z2R ))
−j
r2 r2
1/2
onde z2R = j r2 /dR . De maneira similar na interface entre o revestimento e os ímãs, de
maneira que em r = r1 tem-se
j 1/2 j 1/2
− (c3 Ip 0 (z1R ) + c4 Kp 0 (z1R )) = − (c5 Ip 0 (z1i ) + c6 Kp 0 (z1i ))
µR dR µi di
3.3. Solução da equação de difusão complexa 33

e
p p
−j (c3 Ip (z1R ) + c4 Kp (z1R )) = −j (c5 Ip (z1i ) + c6 Kp (z1i ))
r1 r1
em que z1R = j 1/2 r1 /dR e z1i = j 1/2 r1 /di .
Na interface restante em r = rh , entre os ímãs e a coroa do rotor,

j 1/2 j 1/2
− (c5 Ip 0 (zhi ) + c6 Kp 0 (zhi )) = − c7 Ip 0 (zhe )
µi di µe de
e
p p
−j (c5 Ip (zhi ) + c6 Kp (zhi )) = −j c7 Ip (zhe )
rh rh
em que zhi = j 1/2 rh /di e zhe = j 1/2 rh /de .

Resolução do sistema

A partir disso, o problema pode ser escrito como o sistema abaixo.

a1 c1 + a2 c2 = µo K
a3 c1 + a4 c2 + a5 c3 + a6 c4 = 0
a7 c1 + a8 c2 + a9 c3 + a10 c4 = 0
a11 c3 + a12 c4 + a13 c5 + a14 c6 = 0
a15 c3 + a16 c4 + a17 c5 + a18 c6 = 0
a19 c5 + a20 c6 + a21 c7 = 0
a22 c5 + a23 c6 + a24 c7 = 0

em que

a1 = prsp−1 ; a2 = −prs−p−1 ; a3 = −pr2p−1 ;

j 1/2 µo 0 j 1/2 µo
a4 = pr2−p−1 ; a5 = Ip (z2R ); a6 = Kp 0 (z2R );
dR µR dR µR
a7 = r2p ; a8 = r2−p ; a9 = −Ip (z2 R);
j 1/2 0 j 1/2
a10 = −Kp (z2 R); a11 = − Ip (z1R ); a12 = − Kp 0 (z1R ));
µR dR µR dR
j 1/2 0 j 1/2 p
a13 = − Ip (z1i ); a14 =− Kp 0 (z1i ); a15 = −j Ip (z1R );
di di r1
p
a16 = −j Kp (z1R )); a17 = −Ip (z1i ); a18 = −Kp (z1i );
r1
1/2
j j 1/2 j 1/2 µi 0
a19 =− Ip 0 (zhi ); a20 = − Kp 0 (zhi ); a21 = − Ip (zhe );
di di de µe
a22 = Ip (zhi ); a23 = Kp (zhi ); a24 = −Ip (zhe );
34 Capítulo 3. Perdas em máquinas elétricas

A resolução do sistema é obtida através da álgebra linear, e os resultados são

expressos com a ajuda dos seguintes coeficientes auxiliares:

b1 = a19 − a21 a22 /a24 ; b2 = a20 − a21 a23 /a24 ; b3 = −b1 /b2 ;
b4 = a13 + a14 b3 b5 = a17 + a18 b3 ; b6 = a11 − b4 a15 /b5 ;
b7 = a12 − b4 a16 /b5 ; b8 = a9 − b6 a10 /b7 ; b9 = a5 − b6 a6 /b7 ;
b10 = a3 − b9 a7 /b8 ; b11 = a4 − b9 a8 /b8 ; b12 = a1 − b10 a1 /b11 .

e por fim
µo K b10 a7 a8
c1 = ; c2 = − c1 ; c3 = − c1 − c2 ;
b12 b11 b8 b8
b6 a15 a16
c4 = − c3 ; c5 = − c3 − c4 ; c6 = −b3 c5 ;
b7 b5 b5
a22 a23
c7 = − c5 − c6 .
a24 a24

Mediante a solução apresentada, as perdas podem ser calculadas via 3.7.

 35

4 Modelagem e simulação

A computação dos valores de perdas envolve a definição precisa da geometria

e condições de alimentação envolvidos, bem como a aplicação dos métodos numéricos
pertinentes. Nesse trabalho, utilizaram-se os softwares comerciais SimcenterT M SPEED
e SimcenterT M MAGNET para a definição da geometria em formato digital; obtenção
das harmônicas no tempo das formas de onda de corrente e harmônicas espaciais da
distribuição de ampere-condutor para a aplicação do modelo de perdas apresentado na
seção 3.3; a definição e refino de malhas; condições de contorno; análises transientes por
elementos finitos em 2D e 3D; e análise de resultados.
Abaixo, são comentadas as características dos softwares utilizados, os detalhes da
geometria da máquina e as características utilizadas para o cálculo de perdas e para as
simulações transientes.

4.1 Softwares de simulação

SimcenterT M SPEED

O software Simcenter SPEED é um software científico para solução numérica de

equações íntegro-diferenciais.

Figura 8 – Interface do SimcenterT M SPEED

Fonte: autoria própria

36 Capítulo 4. Modelagem e simulação

SimcenterT M MAGNET
O Simcenter MAGNET é um gerador de malha de elementos finitos com um criador
de geometria e pós-processador aberto ao acoplamento de problemas físicos (eletromagné-
tico, térmico, etc.) com resolução por métodos numéricos. A interface do programa pode
ser vista na figura 3.

Figura 9 – Interface do SimcenterT M MAGNET

Fonte: autoria própria

4.2 Máquina considerada

O modelo utilizado para a realização das simulações é um motor trifásico de 40
kW, sem escovas, com ímãs permanentes superficiais, 4 polos, rotor interno e alimentação
senoidal. As dimensões da máquina são indicadas na tabela 1. No apêndice A há uma
representação dessas grandezas na máquina em corte.

Tabela 1 – Dimensões da máquina

RADSH 20 mm
R1 50 mm
R2 100 mm
GAP 3 mm
LM 11 mm
RVS 1 mm
TGD 2 mm
TWS 10 mm
SOang 20o
Comprimento do pacote 100 mm

Na figura 10 é apresentada a visão da máquina como um todo.

4.2. Máquina considerada 37

Figura 10 – Visão frontal e axial da máquina

Fonte: autoria própria

Materiais
Os materiais utilizados são parte da biblioteca padrão dos softwares. O ímã é feito
da ferrite NeIGT, enquanto o revestimento utiliza o aço 304.

Tabela 2 – Dados da ferrite NeIGT

Temperatura (o C) Permeabilidade Relativa Coercitividade (A/m)

-60 1,07352 -344209
-20 1,07361 -320662
20 1,07211 -297016
60 1,07294 -273004
100 1,07285 -249337
Resistividade Elétrica (Ω.m) 1, 436 ∗ 10−6

Tabela 3 – Dados do aço 304

Permeabilidade relativa 1,008

Resistividade Elétrica (Ω.m) 7, 2 ∗ 10−7

Além disso, não foi considerada nenhuma espécie de laminação nesses componentes.
Eles são tratados como condutores sólidos para o cálculo de correntes induzidas.

Controle
Foi considerado um controle senoidal com tensão de alimentação Vs = 200 V , ponto
de operação da corrente Isp = 200 A e velocidade de rotação de 6000 RPM.
38 Capítulo 4. Modelagem e simulação

4.3 Simulação em Elementos Finitos: Configurações e Pa-

râmetros
Para a simulação em elementos finitos, a malha foi gerada automaticamente, com
maior refinamento nos ímãs e no revestimento do rotor; as malhas são recalculados em
cada passo de tempo. A figura 11 apresenta as malhas usadas para o cálculo em t = 0 ms,
desconsiderando o envoltório de ar. A região que envolve o rotor foi definida como uma
região de rotação, permitindo o movimento durante as simulações transientes.

Figura 11 – Malhas utilizadas para as simulações em elementos finitos em t = 0 ms

Fonte: autoria própria

A tabela 4 apresenta os parâmetros utilizados para a simulação transiente em duas

e três dimensões. Em ambos os casos, o cálculo das perdas é realizado sobre os últimos 10
ms para garantir a contabilização durante operação estável.
4.3. Simulação em Elementos Finitos: Configurações e Parâmetros 39

Tabela 4 – Parâmetros da simulação transiente

Parâmetro 2D 3D
Tempo total 15 ms 20 ms
Passo de tempo 0,1 ms 0,42 ms
Critério de convergência 0,01% 0,01%
Ordem polinomial 2 2

Em ambos os casos, foi escolhido um tempo de simulação maior do que o período

de 2 ciclos elétricos para eliminar quaisquer efeitos transientes do início do movimento.
Para a simulação 3D, esse tempo foi um pouco maior, devido à maior complexidade os
efeitos envolvidos e o menor passo de tempo utilizado.
Mediante estimativas pelos valores de frequência analisados, bem como ajustes
iterativos, o passo de tempo de 0,1 ms se mostrou um compromisso válido entre velocidade
de simulação e consistência das respostas. No caso tridimensional foi utilizado um passo
maior para diminuir o tempo de simulação, visto o grande requerimento computacional de
geometrias mais complexas.
O critério de convergência de 0,1 ms é o valor padrão utilizado para esse tipo de
análise e mostrou bom desempenho em simulação iniciais. Portanto, manteve-se o valor.
No que diz respeito à ordem polinomial, elementos de segunda e terceira ordem
apresentaram pouca diferença entre si na magnitude e distribuição dos campos calculados
em uma simulação magnetoestática. Uma ordem polinomial igual a 1, por outro lado,
levou a uma representação grosseira dos campos. Por esses motivos, e levando em conta
a necessidade de reduzir ao máximo o custo computacional da solução, utilizaram-se
polinômios de segunda ordem.
 41

5 Resultados

5.1 Natureza das correntes induzidas

Uma análise utilizando o MEF2D com variação da resistividade dos ímãs e do
revestimento permite a constatação que o fator limitante para o surgimento de correntes
nos ímãs, nas principais frequências harmônicas consideradas, é a sua resistência elétrica.
Ou seja, o aumento da sua resistividade leva a quedas expressivas no valor das correntes
induzidas, e consequentemente das perdas na região.
Tabela 5 – Valor médio das perdas nos ímãs por correntes induzidas para diferentes valores
de resistividade

Resistividade (Ω.m) Perdas

1,437*10−7 72,8 W
1,437*10−6 10,52 W

5.2 Perdas calculadas

As perdas foram calculadas considerando a alimentação senoidal, de forma que o
fluxo magnético existente corresponde à variação temporal causada pela alimentação das
bobinas, e no cálculo por métodos numéricos devido à interação dessa com o campo dos
ímãs em movimento.
Os resultados numéricos das perdas médias, em watts, obtidos pela resolução
analítica, por MEF-2D e 3D são apresentados na tabela 6.
Tabela 6 – Valor médio das perdas no tempo por correntes induzidas

Região Analítico 2D 3D
Revestimento 88,3 W 153,6 W 51,7 W
Ímãs 60,9 W 72,8 W 46,41 W

Observa-se que as perdas são mais acentuadas no revestimento do rotor, tendo

em vista sua condutividade e a proximidade com o entreferro. Entretanto, vale notar que
o revestimento, além da fixação mecânica, atua como uma barreira de fluxo, atenuando
a penetração das correntes induzidas nos ímãs ao custo de aumentá-las em si. Para
comparação, o valor médio das perdas nos ímãs para uma resolução em MEF2D sem o
revestimento é apresentado na tabela 7.
Observa-se que o valor das perdas nos ímãs são maiores sem o revestimento, como
esperado.
42 Capítulo 5. Resultados

Tabela 7 – Comparação entre os valores de perdas nos ímãs com e sem o revestimento

Presença do Revestimento Perdas médias nos ímãs

Sim 72,8 W
Não 80 W

5.2.1 Distribuição das perdas

Figura 12 – Mapa de perdas médias no rotor resultante do MEF3D

Fonte: autoria própria

Na figura 12 é possível perceber que as perdas se concentram nas regiões centrais

dos ímãs, bem como nas regiões de quina. Isso pode ser explicado pelo fato de que o
campo magnético é mais intenso próximo ao entreferro, onde ocorre a interação com o
campo magnético gerado pelas bobinas do estator. A maior densidade de fluxo magnético
na região central resulta em maiores variações de fluxo magnético durante o ciclo de
operação, o que aumenta a indução de correntes de Foucault e, consequentemente, as
perdas. A distribuição de corrente nas seções transversais do ímã também pode influenciar
a concentração das perdas na região central.
Quanto às bordas, é esperado uma maior variação do fluxo magnético na região,
viso a modulação do fluxo devido à separação entre os ímãs, de maneira semelhante ao
fenômeno de modulação do fluxo pela ranhura, ilustrado em 5.

[Link] Perdas instantâneas

Podemos analisar o valor instantâneo das perdas para compreender fenômenos

relacionados ao período de repetição das perdas e à operação em regime permanente do
5.2. Perdas calculadas 43

sistema. A análise do comportamento apresentado nas figuras 13 e 14 demonstra que o

cálculo das médias de perdas foi, de fato, realizado em um período com regime estável
de operação. Além disso, é possível observar a repetição do “ciclo de perdas” em uma
frequência próxima à frequência de rotação mecânica do rotor.

Figura 13 – Perdas instantâneas em um dos ímãs - MEF2D

Fonte: autoria própria

Figura 14 – Perdas instantâneas no revestimento - MEF2D

Fonte: autoria própria

44 Capítulo 5. Resultados

5.2.2 Comparação entre os métodos de cálculo

A diferença percentual do método analítico utilizado e o valor da simulação em
duas e três dimensões é apresenta na tabela 8.

Região Dif. 2D Dif. 3D

Revestimento 42,5% -70,79%
Ímãs 16,3% -31,22%

Tabela 8 – Diferença percentual entre métodos de cálculo das perdas por correntes induzi-
das

Existem algumas conjecturas que podem explicar por que o cálculo de perdas por
elementos finitos pode resultar em um valor total significativamente diferente do que uma
aproximação analítica.
A simulação por elementos finitos considera detalhes geométricos complexos, como
o formato exato das peças e a presença de bordas, cantos e descontinuidades. Esse tipo de
comportamento é aproximado, na solução apresentada, com um cilindro contínuo. Esses
detalhes podem levar a um aumento das perdas devido a efeitos como o aumento das
correntes parasitas e a concentração do campo magnético.
Em outro ponto de vista, a alta complexidade e o detalhamento necessários para a
modelagem da geometria também tornam a simulação em três dimensões extremamente
suscetível a erros, sejam de origem da modelagem, de erros numéricos ou de gargalos
nas propriedades dos materiais; de forma que pequenos detalhes podem invalidar a
confiabilidade da resolução.
Além disso, o cálculo analítico considera uma resposta linear do material, enquanto
a simulação por elementos finitos considera o comportamento não linear dos materiais
magnéticos, como a curva de magnetização não linear do aço elétrico. Isso pode resultar
em perdas adicionais que não são consideradas na aproximação analítica.
Uma melhoria significativa na aproximação do comportamento pelas equações
utilizadas seria a consideração de seccionamento no cilindro que representa os ímãs,
caracterizando uma barreira no caminho das correntes entre as regiões de um mesmo
cilindro.

5.3 Estratégias de mitigação

Nota-se que os valores obtidos para as perdas representam uma porção significativa
da potência da máquina. Embora seja um resultado esperado, visto que a geometria e
controle escolhido não possuem em vista um motor otimizado, algumas estratégias podem
ser empregadas para a redução do valor dessas perdas. Duas abordagens comuns são
5.3. Estratégias de mitigação 45

o aumento da resistividade dos materiais utilizados e o seccionamento dos ímãs e do

revestimento ao longo do eixo.
Como apresentado na seção 5.1, as correntes de Foucault analisadas são limitadas
pela resistividade dos materiais empregados. Portanto, a escolha de materiais com maior
resistividade elétrica, como aços elétricos especiais ou ligas com maior teor de resistividade,
são uma forma efetiva da redução dessas perdas [Han et al. 2018].
O seccionamento dos ímãs e do revestimento ao longo do eixo é outra estratégia
eficaz para mitigar as perdas [Siemens 2022] [Hendershot e Miller 2010]. Ao dividir os ímãs
e o revestimento em seções ao longo do eixo do rotor, o caminho percorrido pelas correntes
parasitas é interrompido. Essa interrupção impede que as correntes circulem livremente e
reduz consideravelmente as perdas por correntes de Foucault [Bastos 2004]. Além disso, o
seccionamento também pode ajudar a minimizar o efeito da saturação magnética.
A escolha das estratégias de mitigação das perdas nos ímãs e no revestimento do
rotor depende das características específicas do motor e dos requisitos de desempenho.
Cada abordagem deve ser avaliada em termos de custo, complexidade de fabricação e
impacto nas propriedades magnéticas do motor.
 47

6 Conclusão

Com base nas análises realizadas neste estudo sobre as perdas por correntes de
Foucault em ímãs de motor sem escovas de rotor interno e ímãs superficiais, foram
identificadas diferenças significativas entre o modelo analítico proposto e as simulações em
elementos finitos.
O modelo analítico apresentado demonstrou ser uma ferramenta valiosa para estimar
as perdas de forma rápida, mas esbarra em limitações. Ao comparar esses resultados com as
simulações em elementos finitos, observou-se que as simulações forneceram uma avaliação
mais detalhada das perdas, considerando a complexidade da geometria do motor e a
distribuição espacial das correntes de Foucault.
Essas diferenças ressaltam a importância de utilizar abordagens complementares
na análise das perdas. O modelo analítico pode ser uma opção viável para estimativas
rápidas e simplificadas, especialmente durante a fase de projeto inicial. Por outro lado,
as simulações em elementos finitos permitem uma análise mais precisa e abrangente das
perdas, considerando aspectos geométricos e físicos mais complexos.
Durante o desenvolvimento deste estudo, foi possível obter um maior entendimento
sobre as áreas nas quais as perdas por correntes de Foucault se concentram no rotor de
um motor sem escovas. A análise em elementos finitos e as simulações permitiram uma
visualização detalhada das regiões onde ocorrem as correntes induzidas e, consequentemente,
onde as perdas são mais significativas.
Além disso, foi possível investigar a natureza cíclica das perdas instantâneas por
correntes de Foucault que ocorrem nos ímãs permanentes superficiais durante a operação do
motor sob regime de alimentação senoidal. Essas perdas apresentam um padrão periódico,
seguindo o ciclo da corrente elétrica aplicada. Outro aspecto relevante observado neste
estudo é que as correntes induzidas são limitadas pela resistência do material e não pela
sua indutância. Isso significa que o aumento na resistividade dos materiais utilizados no
rotor pode ser uma estratégia eficaz para reduzir as perdas por correntes de Foucault.
Em suma, este estudo contribuiu para o avanço do conhecimento sobre as perdas
por correntes de Foucault em motores sem escovas, ressaltando a importância de utilizar
abordagens analíticas e de simulação em elementos finitos complementarmente. Essas
ferramentas são essenciais para o projeto e otimização de motores eficientes e de alto
desempenho, levando a avanços significativos na indústria de motores elétricos.
 49

Referências

BARCELOS, P. C. S. Tese de mestrado, PROJETO E OTIMIZAÇÃO DE UM

MOTOR BLDC DE ROTOR INTERNO COM ÍMÃS CERÂMICOS NA SUPERFÍCIE.
Florianópolis, SC: [s.n.], 2021. Citado 3 vezes nas páginas 17, 18 e 21.

BASTOS, J. Eletromagnetismo para engenharia: estática e quase-estática. [S.l.]: Ed.

da UFSC, 2004. (Serie Didatica Universidade Federal de Santa Catarina). ISBN
9788532803061. Citado 5 vezes nas páginas 26, 28, 29, 30 e 45.

GIERAS, J.; WING, M. Permanent Magnet Motor Technology: Design and Applications.
[S.l.]: Marcel Dekker, 1996. Citado na página 17.

HAN, T. et al. Investigation of various rotor retaining sleeve structures in high-speed pm

brushless motors. In: 2018 21st International Conference on Electrical Machines and
Systems (ICEMS). [S.l.: s.n.], 2018. p. 109–114. Citado 3 vezes nas páginas 15, 19 e 45.

HENDERSHOT, J.; MILLER, T. Design of Brushless Permanent-magnet Machines. [S.l.]:

Motor Design Books, 2010. ISBN 9780984068708. Citado 5 vezes nas páginas 15, 18, 19,
30 e 45.

JILES, D.; ATHERTON, D. Theory of ferromagnetic hysteresis. Journal of Magnetism

and Magnetic Materials, v. 61, n. 1, p. 48–60, 1986. ISSN 0304-8853. Disponível em:
<[Link] Citado na
página 21.

MILLER, T.; LAWRENSON. Penetration of transient magnetic fields through conducting

cylindrical structures, with particular reference to superconducting a.c. machines. Proc.
IEE, v. 123, n. 5, p. 437–443, 1976. Citado 2 vezes nas páginas 23 e 27.

PARIS Agreement. Disponível em: <[Link]

TREATY&mtdsg_no=XXVII-7-d&chapter=27&clang=_en>. Citado na página 15.

SIEMENS. Eddy-current Losses in Magnets and Rotor Cans. [S.l.], 2022. Trabalho não
publicado. Citado 2 vezes nas páginas 15 e 45.

STOLL, R. The Analysis of Eddy Currents. [S.l.]: Clarendon Press, 1974. (International
Series of Monographs in Electrical Engineering). ISBN 9780198593119. Citado na página
26.

STOLL, R. L.; HAMMOND, P. Calculation of the magnetic field of rotating machines.

part 4: Approximate determination of the field and the losses associated with eddy
currents in conducting surfaces. In: . [S.l.: s.n.], 1965. Citado na página 27.

UMANS, S. Máquinas Elétricas de Fitzgerald e Kingsley - [Link]. [S.l.]: AMGH Editora,

2014. ISBN 9788580553741. Citado 3 vezes nas páginas 15, 17 e 20.
 51

APÊNDICE A – Indicação das

dimensões da máquina

Figura 15 – Ilustração das dimensões da máquina

Fonte: autoria própria