(2021) CapítuloDeLivroEditoraPascal SomenteOCapítulo10
(2021) CapítuloDeLivroEditoraPascal SomenteOCapítulo10
(2021) CapítuloDeLivroEditoraPascal SomenteOCapítulo10
MODELAGEM MATEMÁTICA
E ANÁLISE DE UM MOTOR A
RELUTÂNCIA VARIÁVEL 8/6
ATRAVÉS DA APLICAÇÃO DO MÉTODO
DE ELEMENTOS FINITOS
E
ste trabalho tem o propósito de contribuir com o estudo de máquinas a relutân-
cia variável através da aplicação do método de elementos finitos. Por meio da
utilização do software FEMM é analisada a distribuição de fluxo magnético nas
lâminas do motor e é obtido o comportamento da indutância ao variar a corrente
e o ângulo da posição do rotor. Tomando-se como parâmetros de entrada os dados
obtidos da análise magnetostática no FEMM, é simulado o sistema de acionamento
do motor através da utilização da ferramenta Simulink, sendo possível verificar o
comportamento dinâmico da máquina ao se obter a curva do torque mecânico da
mesma.
Abstract
T
his paper has the purpose of contributing to the study of switched reluctan-
ce machines through the application of the finite element method. Using the
software FEMM, the magnetic flux distribution on the laminated sheets of the
motor is analyzed and the behavior of the inductance is obtained by varying the
current and the angle of the rotor position. Taking as input parameters the data ob-
tained from the magnetostatic analysis of FEMM, the motor drive system is simu-
lated through the use of the Simulink tool, and it is possible to verify the dynamic
behavior of the machine by obtaining the mechanical torque curve of it.
Engenharia 4.0
A era da produção inteligente - Vol. 06
148
1. INTRODUÇÃO
Embora fosse conhecida desde o início dos estudos das máquinas elétricas
e suas primeiras aplicações fossem registradas no século XIX, a Máquina a Re-
lutância Variável (MaRV) teve um desenvolvimento tardio em relação aos outros
tipos de máquinas e seu potencial foi pouco explorado, isso ocorreu principalmente
devido às dificuldades envolvendo seu acionamento (MILLER, 1993). No entanto,
a partir de 1960, com o avanço da tecnologia dos dispositivos semicondutores de
chaveamento como o MOSFET de potência e o IGBT, advindos dos estudos da fís-
ica do estado sólido, foram desenvolvidos diversos tipos de circuitos eletrônicos de
potência que viabilizaram o acionamento eficiente das MaRVs (ANDERSON, 2001).
Atualmente, com a evolução dos microcontroladores e das técnicas de controle é
possível construir MaRVs sofisticadas com acionamentos de alto desempenho.
Dessa forma, este trabalho tem como um dos objetivos realizar a análise mag-
netostática de um MRV com topologia 8/6 (8 polos no estator e 6 polos no rotor)
por meio da aplicação do MEF pelo software FEMM, através do qual é possível re-
alizar o levantamento das curvas características de fluxo magnético e indutância
da máquina em função do ângulo da posição do rotor. Adicionalmente, tem-se o
Engenharia 4.0
A era da produção inteligente - Vol. 06
150
A Figura 1 foi obtida a partir da modelagem do MRV 8/6 feita no software Auto-
CAD com base nas dimensões físicas das lâminas do estator e do rotor da máquina
que foram obtidas seguindo-se os procedimentos de cálculos desenvolvidos em
(KRISHNAN, 2001) e (MILLER, 2001), que abordam projetos de máquinas a relu-
tância variável. Estes cálculos estão detalhados em (MORAES FILHO, 2017), que
foi a dissertação utilizada como referência para este trabalho.
di
R f i + L(θ , i ) + e
v=
(1)
dt
λ (θ , i ) = L(θ , i )i
(2)
∂L(θ , i )
e = iω
m (3)
∂θ
di ∂L(θ , i )
R f i + L(θ , i ) + iωm
v=
(4)
dt ∂θ
pe = vi
di ∂L(θ , i )
R f i 2 + iL(θ , i ) + i 2ωm
pe = (5)
dt ∂θ
d 1 1 ∂L(θ , i )
R f i 2 + L(θ , i )i 2 + i 2
pe =
ωm (6)
dt 2 2 ∂θ
1 ∂L(θ , i )
pent = i 2
ωm (7)
2 ∂θ
pent 1 2 ∂L(θ , i )
T=
emg =
i (8)
ωm 2 ∂θ
Como todo o equacionamento realizado até o momento foi feito para somente
uma fase da máquina, tem-se que para o MRV 8/6 tetrafásico a expressão para o
cálculo do torque eletromagnético total será composta pela somatória dos torques
individuais de cada fase (a, b, c e d), calculado por (9):
Engenharia 4.0
A era da produção inteligente - Vol. 06
152
d ωm
Tmec= Temg −total − Dωm − J
(10)
dt
va ia ia
v i
b b ib
vc i i
=
R c + L c
vd id id
(11)
Tmec ωm ω
m
0 θ θ
ra 0 0 0 0 0
0 r 0 0 0 0
b
0 0 rc 0 0 0
R=
(12)
0 0 0 rd 0 0
m n p q − D 0
0 0 0 0 −1 0
sendo as variáveis m, n, p e q calculadas através de (13), (14), (15) e (16):
1 ∂L (θ , ia )
m = ia a
2 ∂θ (13)
1 ∂L (θ , ib )
n = ib b
(14)
2 ∂θ
1 ∂L (θ , ic )
p = ic c
2 ∂θ (15)
1 ∂L (θ , id )
q = id d
(16)
2 ∂θ
[ I] [ L] [V ] − [ L] [ R][ I ]
= −1
(18)
−1
Feito isso, os materiais de cada parte do MRV foram definidos. Para todo o
entreferro foi utilizado o ar, no eixo da máquina foi utilizado o aço 1020, nas bobi-
nas foram utilizados fios de cobre AWG 14 e nas lâminas do rotor e do estator foi
utilizado o aço E185. Estando definido os materiais, foram designados os circuitos
elétricos a serem analisados para solução do problema. Como MRV é 8/6 (tetrafási-
co), foram definidos 4 circuitos (fases A, B, C e D). O número de espiras envolvidas
em cada polo do motor foi igual a 55, ou seja, 110 voltas por fase, conforme foi
projetado em (MORAES FILHO, 2017).
Engenharia 4.0
A era da produção inteligente - Vol. 06
154
fronteira. Como a relutância do material ferromagnético utilizado (E185) nas lâmi-
nas é muito menor do que a relutância do ar, foi definido sem perda de precisão, a
fronteira do problema como sendo a circunferência externa da lâmina do estator.
Figura 4 – Distribuição do fluxo magnético na estrutura do MRV 8/6 para a posição de 15º
Figura 5 – Curvas de fluxo magnético versus ângulo da posição do rotor para alguns valores de corrente
fixados
Fonte: Autoria própria
Engenharia 4.0
A era da produção inteligente - Vol. 06
156
Figura 6 – Perfil em 3D do fluxo magnético em função da corrente e do ângulo da posição do rotor
Fonte: Autoria própria
Figura 7 – Curvas de indutância versus ângulo da posição do rotor para alguns valores de corrente fixados
Fonte: Autoria própria
Engenharia 4.0
A era da produção inteligente - Vol. 06
158
Figura 9 – Circuito elétrico do conversor meia ponte assimétrica
Fonte: Moraes Filho (2017)
Engenharia 4.0
A era da produção inteligente - Vol. 06
160
5. CONCLUSÕES
Neste trabalho, foi feita inicialmente uma breve apresentação sobre o princípio
de funcionamento da MaRV. Em seguida, foi realizada a modelagem matemática
da máquina obtendo-se as equações que descrevem o comportamento elétrico e
mecânico da mesma. Feito isso, com base nas medidas determinadas para as lâ-
minas do rotor e estator foi modelado em 2D o motor no software AutoCAD. Tendo
a geometria da máquina, essa foi exportada para o software FEMM, no qual foram
realizadas as simulações via método dos elementos finitos que possibilitou que
fossem levantadas as curvas de indutância e de fluxo magnético de uma fase em
função do ângulo da posição do rotor.
Todas essas etapas foram cumpridas com sucesso e os dados coletados, que
se mostraram coerentes, servem de validação para o projeto do motor. Neste con-
texto, a fim de comprovar a coerência entre as curvas obtidas via elementos finitos
e os dados reais do motor durante seu funcionamento na prática, recomenda-se o
desenvolvimento de uma metodologia de ensaio para realizar este levantamento
de forma experimental.
Referências
KRISHNAN, R. Switched Reluctance Motor Drives: Modeling, Simulation, Analysis, Design, and
Applications, Boca Raton, FL: CRC Press, ISBN-13: 978-0849308383, 2001.
MEEKER, D. Finite Element Method Magnetics: Version 4.2: User’s Manual, 2015. Disponível em:
http://www.femm.info/Archives/doc/manual42.pdf Acesso em: 9 de julho de 2018.
MILLER, T. J. E. Electronic Control of Switched Reluctance Machines, Newnes Power Engineering Se-
ries, ISBN-13: 978-0080505244, Oxford, Reino Unido, 2001.
MILLER, T. J. E. Switched Reluctance Motors and Their Control, Magna Physics Publishing and Claren-
don Press, Oxford, Reino Unido, 1993.
UMANS, S. D., Máquinas Elétricas de Fitzgerald e Kingsley, 7. ed., AMGH Editora Ltda, ISBN-10:
8580553733, 2014.
Engenharia 4.0
A era da produção inteligente - Vol. 06
162