UNIVERSIDADE FEDERAL DE UBERLÂNDIA

FACULDADE DE ENGENHARIA ELÉTRICA

3º Trabalho de Elementos Finitos

Análise eletrostática de um eletroímã modelado em Lua script e

executado no FEMM
Terceiro Relatório apresentado ao prof. PhD José
Roberto Camacho como parte das atividades avaliativas
da disciplina Elementos Finitos.

Walisson Martins Teixeira

11511EEL024

UBERLÂNDIA – MG

2022

1
Sumário

1. Introdução .............................................................................................................. 3

2. Objetivo ................................................................................................................. 3

3. Desenvolvimento .................................................................................................... 4

3.1 Antes do processamento pelo FEMM:......................................................................... 5

6. (código LUA) ......................................................................................................... 5

3.2 Depois do processamento pelo FEMM: ....................................................................... 7

4. Conclusões ........................................................................................................... 11

5. Referências Bibliográficas ...................................................................................... 11

2
 1. Introdução
O software LUA Script para desenvolver a parte que antecede o processamento em
código, posteriormente executaremos o código no FEMM e chegaremos à resolução do
problema.

A utilização de simulações computacionais apresentam excelentes resultados, além de

possibilitarem a realização de estudos de problemas complexos, que dificilmente poderiam ser
analisados considerando a aplicação das simplificações da teoria de circuitos magnéticos. Um
destes simuladores, o software FEEM, utiliza a técnica dos elementos finitos na resolução de
problemas magnéticos e eletrostáticos, possibilitando ainda a comunicação com outros
softwares de simulações computacionais, como o Octave, Matlab e Mathematica, o que
viabiliza os pós processamento dos resultados obtidos.

O pré-processamento no software FEMM pode ser realizado inserindo-se os diversos objetos

existentes no FEMM via teclado e mouse ou de forma mais automática utilizando-se scripts,
desenvolvidos em linguagem de programação Lua. O segundo método apresenta inúmeras vantagens, já
que a linguagem de programação permite um desenvolvimento automatizado de diversas modelagens,
tornando viável a obtenção de diversos tipos de resultados.

2. Objetivo
O problema utilizado nesta terceira etapa será o mesmo utilizado na etapa anterior, todo
o desenvolvimento e processamento do problema eletrostático em código, será executado no
FEMM até a resolução dele.

O problema é um capacitor que possui placas planas e paralelas, onde todo o

processamento será desenvolvido em Lua.

3
 3. Desenvolvimento
O desenho do capacitor de placas paralelas na Figura 1:

Figura 1

Cálculo da capacitância:

𝑄 𝑆
𝐶= = 𝜀
𝑉 𝑑

Onde:

C = capacitância do capacitor

Q = carga total armazenada

V = diferença de potencial entre as placas

𝜀 = permissividade do material dielétrico

S = área de uma placa

D = distância entre as placas

4
3.1 Antes do processamento pelo FEMM:
Na Figura 2 está definido as coordenadas cartesianas para cada ponto em questão, todo
o problema está sendo analisado em milímetros

Figura 2

Será definido na simulação os valores de borda, designando a fronteira de simulação

dos elementos finitos. O dielétrico escolhido para o capacitor segue sendo o ar e os potenciais
definidos para as placas superior e inferior foram respectivamente de: +100V e -100V.

Abaixo apresenta o código desenvolvido no software LUA para modelar o problema no

FEMM.

6. (código LUA)
clearconsole()

-- Iniciando projeto Lua script

newdocument(1)

messagebox("Capacitor no Lua Script")

5
ei_probdef("millimeters","planar","1e-8",20,30,"Succ. Approx")

ei_saveas("D:\\Elementos Finitos\\capacitor.fee")

-- Criando nós

messagebox("Inserir nos")

ei_addnode(0,0)

ei_addnode(60,0)

ei_addnode(60,101)

ei_addnode(0,101)

ei_addnode(0,50)

ei_addnode(10,50)

ei_addnode(10,51)

ei_addnode(0,51)

-- Criando segmentos

messagebox("Inserir segmentos")

ei_addsegment(0,0,60,0)

ei_addsegment(60,0,60,101)

ei_addsegment(60,101,0,101)

ei_addsegment(0,101,0,0)

ei_addsegment(0,50,10,50)

ei_addsegment(0,51,10,51)

-- Adicionando bloco de propriedades de materiais

ei_addblocklabel(33,82)

-- Importando materiais da biblioteca FEMM

ei_getmaterial("Air")

6
ei_selectlabel(32,82)

ei_setblockprop("Air", 0, 1, "None",0,0, 1)

Representação do problema no FEMM.

Figura 4 – Representação do problema no FEMM

3.2 Depois do processamento pelo FEMM:

A Figura 5 nos mostra que foi criado 6921 nós para este problema, o que coincidiu com
o resultado apresentado na primeira etapa onde o software LUA não foi utilizado.

7
Figura 5 – Processamento

8
 A Figura 6 nos mostra o comportamento da densidade de tensão em volts para este
capacitor.

Figura 6 – Densidade de tensão no capacitor

Na Figura 7 podemos observar o comportamento do campo elétrico, sendo constante

entre as placas do capacitor. Também podemos observar o efeito de borda, situação que
desconsiderada na análise teórica.

9
 Figura 7 – Campo elétrico

Na Figura 8 podemos ver as linhas equipotenciais dentro e fora do capacitor.

Figura 8 – Linhas equipotenciais

Cálculo da capacitância de forma teórica:

𝑄 𝑆
𝐶= = 𝜀
𝑉 𝑑

(20 ∗ 10−3 )2
𝐶 = (8.854 ∗ 10−12 ) ∗
1 ∗ 10−3

𝐶 = 3.5416 ∗ 10−12 𝐹

Cálculo da capacitância por elementos finitos:

O valor da carga encontrado através da simulação foi:

𝑄 = 3.90496 ∗ 10−10 𝐶

É necessário multiplicar o valor da carga por dois, já que no momento anterior a

simulação a placa modelada continha só metade da área da placa do capacitor teórico. Assim

10
sendo, se multiplicarmos a área por dois também iremos dobrar a carga. Nesse contexto temos
que a carga total é:

𝑄 𝑡 = 2 ∗ 3.90496 ∗ 10−10 𝐶

𝑄 𝑡 = 7.80992 ∗ 10−10 𝐶

A diferença de potencial:

𝑉 = 100 − (−100) = 200 𝑉

A capacitância:

𝑄𝑡
𝐶 = = 3.90496 ∗ 10−12 𝐹
𝑉

4. Conclusões
Neste terceiro trabalho continuamos a representação do capacitor para representar a
solução do problema utilizando o FEMM a partir de um script desenvolvido no software
LuaEdit 2010.

Além do mais, a possibilidade de aplicação do método dos elementos finitos no FEMM a partir
do desenvolvimento de códigos em linguagem de programação Lua evidenciam o potencial de
utilização do software FEMM no auxílio à solução de diferentes tipos de problemas
relacionados aos fenômenos eletrostáticos.

5. Referências Bibliográficas
6. [1]https://www.etapa.com.br/etaparesolve/etaparesolve/2015/Fuvest/2Fase_217/resolu
cao/2fuv15f.pdf. Acesso em: 01 de fevereiro de 2022.
7.
8. [2] MEEKER, David Meeker. Finite Element Method Magnetics: OctaveFEMM.
Versão 1.2. Março de 2009.
9. [3] CHAPMAN, Stephen J. Chapman. Fundamentos de Máquinas Elétricas. 5º Edição.
Editora AMGH LTDA, 2013.

11