Peel1924 D
Peel1924 D
Peel1924 D
Florianópolis
2019
Ficha de identificação da obra elaborada pelo autor,
através do Programa de Geração Automática da Biblioteca Universitária da UFSC.
Tedesco, Duani
Projeto ótimo de um gerador síncrono de polos salientes /
Duani Tedesco ; orientador, Patrick Kuo-Peng,
coorientador, Jean Vianei Leite, 2019.
102 p.
Inclui referências.
O presente trabalho em nível de mestrado foi avaliado e aprovado por banca examinadora
composta pelos seguintes membros:
Certificamos que esta é a versão original e final do trabalho de conclusão que foi julgado
adequado para obtenção do título de mestre em Engenharia Elétrica
Assinado de forma digital
Bartolomeu por Bartolomeu Ferreira
Ferreira Uchoa Uchoa Filho:47636211491
Filho:47636211491 Dados: 2019.11.18
10:58:21 -03'00'
____________________________
Prof. Dr.(a) Bartolomeu Ferreira Uchôa Filho
Coordenador(a) do Programa
Assinado de forma digital por
Patrick Kuo Patrick Kuo Peng:00352002948
Peng:00352002948 Dados: 2019.11.14 15:50:23
____________________________
-03'00'
Salient pole synchronous generators (GSPS – in portuguese) are non-serial equipment and have
high added value. This work aims to develop an algorithm to size a generator, determine its
performance and finally optimize it. An analytical model for the machine was developed and
validated by comparing it with numerical simulation based on the finite element method (FEM).
The optimization is based on multiobjective genetic algorithms, being the objectives
maximization of the yield and minimization of the machine cost. Restrictions were also
imposed to determine the generator performance range. To calculate the cost of the machine, a
cost factor variable which takes into account the influence of ferromagnetic material and copper
in the calculation of the machine was defined. At the end of the algorithm execution, a graph is
presented that lists all the solutions found and the non-dominated solutions are highlighted.
From the non-dominated solutions the simulation is made and the best solution for the presented
problem is chosen. Finally, the initial geometry is compared with the non-dominated optimized
geometry.
1 INTRODUÇÃO.......................................................................................................... 15
1.1 OBJETIVOS ............................................................................................................... 17
4.1.2 Seleção........................................................................................................................ 53
6 CONCLUSÃO............................................................................................................ 76
6.1 MELHORIAS E SUGESTÕES DE TRABALHOS FUTUROS ............................... 77
REFERÊNCIAS ....................................................................................................................... 79
1 INTRODUÇÃO
1.1 OBJETIVOS
a) b)
Fonte: O autor.
Fonte: [28].
22
Uma espira pode ser feita com um ou mais condutores em paralelo. A utilização
de condutores em paralelo facilita a formatação da bobina. Na Figura (2-3a) é mostrada
uma ranhura com duas bobinas, característica de enrolamento de camada dupla. Em
detalhe na Figura (2-3b) é mostrado o corte transversal de uma bobina que mostra um
total de 7 espiras, e cada espira é composta por dois condutores de altura ( ℎ ) e base
( ).
O espaço mostrado entre a bobina inferior e a bobina superior e também a
quantidade de isolante da bobina até as paredes da ranhura são calculados em função da
tensão do gerador e tem a função de prevenir fugas de corrente de uma bobina para a
outra e da bobina para a massa-terra [29].
a) b) c)
Fonte: O autor.
23
A bobina vista de cima é mostrada na Figura (2-4), nesta figura é possível ver as
regiões da cabeça da bobina, a parte reta e o passo de bobinagem ( ) que será calculado
nos próximos capítulos.
Fonte: O autor.
Fonte: O autor.
24
Fonte: O autor.
Fonte: O autor.
26
= = = (2-1)
∗
= � = , [Nm] (2-2)
6
Tabela 2-4 Valores típicos de tensão tangencial para máquinas síncronas de polos
salientes refrigerados a ar
Parâmetro Sigla Valor
Tensão tangencial mínima 21000 [Pa]
Tensão tangencial média 33500 [Pa]
Tensão tangencial máxima 48000 [Pa]
Fonte: [26].
= = 1,5553[ ] (2-3)
� �
= = √ = √ = , (2-4)
Tendo em vista que o volume do rotor está relacionado com o diâmetro e com o
comprimento do pacote conforme Equação (2-5), é possível chegar ao valor do diâmetro
e ao comprimento efetivo do pacote, utilizando as equações (2-6) e (2-7) [27].
�
= × × (2-5)
=( ) =1963[mm] (2-6)
�
= = 514[mm] (2-7)
�
� = = , [mm] (2-8)
ℎ
� = �
≈ , − , (2-9)
optou-se por trabalhar com � médio e variá-lo durante a otimização, portanto o valor de
ℎ é dado pela Equação (2-10) [8].
ℎ =� � = [mm] (2.10)
Fonte: O autor.
30
= , � =1 [mm] (2-11)
ℎ
ℎ = ∴ℎ = ℎ� ≈ [ ] (2-12)
�
ℎ ℎ = , � = [mm] (2-13)
,
= �√ �
= (2-15)
Inicialmente foi optado por utilizar o número de bobinas por polo e por fase (q)
igual a 3, e variá-lo durante a otimização.
Com o número de bobinas por polo e por fase é possível determinar a quantidade
de ranhuras ( ) no estator da máquina [8].
= = . (2-16)
� = = , [pu] (2-17)
− �
= . − ,
= [mm] (2-18)
�
� + �
� = = , [ ] (2-19)
pela empresa Automatic. Desta forma, o número total de canais de ventilação para este
gerador é igual a 9. Considerando os canais de ventilação, o comprimento total do pacote
do estator ( ) é mostrado na Equação (2-20).
= + = [mm]. (2-20)
�
= = , [ ] (2-21)
� �
Para determinar a relação entre base e altura das ranhuras é feito o cálculo do
passo de cada ranhura [8].
� + �
� = = ,53[mm] (2.22)
=� � = [mm] (2.23)
ℎ = = [ ] (2.24)
gl �
ℎys = = 9 [mm] (2.25)
ys �
33
= + � = [mm] (2.26)
= + (ℎ +ℎ )≈ [mm] (2.27)
= − = (2.28)
Fonte: O autor.
34
Passo de enrolamento 8
Altura do corpo do polo ℎ 87 [mm]
Altura da sapata do polo ℎ ℎ 35 [mm]
Base do corpo polar 153 [mm]
Base da sapata polar ℎ 244 [mm]
Número de espiras por bobina � 4
Número de caminhos paralelos 1
Largura da ranhura 15 [mm]
Altura da ranhura ℎ 72 [mm]
Fonte: O autor.
�= (3-1)
√ �
37
Para o cálculo do fluxo magnético necessário por polo são utilizados o valor da
tensão nominal em cada fase ( ), fator de bobinagem ( ), a frequência de operação da
máquina (f ) e o número de espiras por caminho e por fase ( ).
Com o valor do fluxo magnético necessário é possível determinar a indução no
circuito magnético da máquina.
O fluxo magnético gerado por polo descrito na Equação (3.1) deve percorrer o
circuito magnético conforme apresentado na Figura (3-1).
As linhas de distribuição do fluxo magnético são mostradas na Figura (3-2). O
fluxo magnético produzido por polo deve ser igual a �, o valor da indução magnética é
calculado conforme Equação (3-2). O valor da seção transversal à direção do fluxo
magnético é representado pela variável s [3]. O método de cálculo da seção é mostrado
no Apêndice 1.
�
�=∫⃗ ∴ = (3-2)
Fonte: O autor.
38
=� ∴ = (3-3)
Fonte: O autor.
Quando comparada a Tabela (3-2) com a Tabela (2-3), que previa as induções
no gerador é possível perceber que houve uma significativa redução nos valores de
indução. Isso ocorreu pois o número total de espiras por fase depende da quantidade de
espiras por bobina. Para que fosse possível calcular o valor de , foi levado em
consideração que a quantidade de espiras em uma bobina deveria resultar sempre em um
valor igual ou maior de espiras em série por fase que o valor de , encontrado na
Equação (2-15). Nos cálculos realizados, a quantidade de espiras em série por fase ficou
1,32 vezes maior que .
Com o valor da indução e levando em consideração a permeabilidade linear é
calculado o campo magnético, mostrado na Tabela (3-3)
=∮⃗ = (3-5)
+ +�
= +
(3-6)
= ( ) (3-7)
Fonte: [34].
42
3.2.3.2 Reatâncias
= + (3-8)
=( )+ (3-9)
� �
= � +
(3-10)
�
� �
= +� (3-11)
� �
� � � �
= +� + (3-12)
� � � �
= + (3-13)
43
√ �
= �
(3-14)
= |�| + � (3-15)
= |�| + � (3-16)
=√ + (3-17)
FMMCarg Ncp
= (3-19)
pol
� = × (3-20)
�� �
− +
� � �
= � �� (3-20)
+ − + − −
� �
Fonte: O autor.
Fonte: O autor.
47
σd dB
� = ∫ dt [W/kg] (3-21)
12 dt
Figura 3-6 Campo magnético criado pela circulação de corrente induzida no ferro
magnético.
48
Fonte: O autor.
3.3.2.3 Perdas excedentes
dB 1,5
�ex = √� ∫ | dt | dt [W/kg] (3-22)
Figura 3-7 Campo magnético criado pelo descolamento das paredes dos
domínios magnéticos
Fonte: O autor.
49
O cálculo das perdas totais no ferro silício pode ser feito obtendo-se uma curva
com a soma do campo magnético originado pela histerese, pelas correntes de Foucault e
perdas excedentes versus a indução magnética, conforme Figura (3-8). A densidade
volumétrica da energia perdida por ciclo pode ser obtida pela integral da área da curva
BH resultante.
Figura 3-8 Curva BH contendo os efeitos das correntes induzidas, histerese e perdas
excedentes.
Fonte: O autor.
O cálculo das perdas foi realizado somente no estator da máquina pois foi
considerado que o fluxo magnético no rotor é contínuo. O circuito magnético do estator
foi separado em dente e coroa, foi utilizada a indução de pico na seção e variado a mesma
senoidalmente. O cálculo da integral da área do laço BH mostrado na Figura (3-8) resulta
na densidade de perda de origem magnética em [J/m³] logo, para o cálculo das perdas
50
totais foi necessário calcular o volume dos dentes e da coroa do estator. Este cálculo
encontra-se no Apêndice 2.
Com o cálculo das perdas magnéticas é possível calcular o rendimento da
máquina. As perdas mecânicas que compõem o cálculo do rendimento foram estimadas
em 3 kW levando em consideração ensaios em máquinas deste porte pela empresa
Automatic. As perdas totais podem ser calculadas conforme Equação (3-23).
Com o valor das perdas totais calculadas é possível então calcular o rendimento
do gerador (Rend), conforme Equação (3-24).
+ +( + )
= = , (3-25)
51
4.1.2 Seleção
4.1.3 Crossover
4.1.4 Mutação
4.1.5 Elitismo
Segundo [2], dada uma população gerada aleatoriamente cada indivíduo será
analisado e uma cópia do melhor ponto, ou seja, do ponto no qual se obteve o melhor
desempenho é guardada em uma posição extra do algoritmo. Em seguida os operadores
genéticos atuam sobre os primeiros pontos e comparam com a cópia guardada. Caso
algum dos indivíduos da nova geração possuir uma função fitness menor que a cópia
guardada, esse valor será substituído.
Após a execução dos processos anteriores, surge uma nova população. Essa
população passará pelo mesmo processo. Inicialmente será testada e caso não atinja a
condição de parada do algoritmo o processo se inicia novamente.
A maioria dos problemas reais possui mais de um objetivo. Por isso, as soluções
multiobjetivo apresentam diversas soluções que se encaixam dentro das restrições
impostas ao problema. Para determinar qual a melhor solução, devem ser definidos as
prioridades e os pesos entre os diversos objetivos de interesse, assim encontrando uma
solução ótima baseada nas informações fornecidas. Também existe a possibilidade de não
56
Fonte: [1].
4.2 ALGORITMO IMPLEMENTADO
Fonte: O autor.
58
Fonte: O autor.
Indução na coroa do
0,8270 0,7991 -3,37
estator
Indução no corpo do
1,10 1,08 -1,818
polo
Indução na coroa do
0,5288 0,55 4,01
rotor
Fluxo máximo por polo
� 0,087 0,085 -2,29
no entreferro
Fonte: O autor.
O cálculo de campos foi realizado por meio do programa FEMM, cujo código
encontra-se no Apêndice 4.
Como condição de contorno, foi aplicada uma caixa de ar quadrada fora do
gerador, na borda desta caixa foi aplicado o vetor potencial magnético nulo no segmento.
Foram utilizados elementos de malha menores no entreferro, a malha de um corte do
gerador é mostrada na Figura (5-1).
Como o programa foi escrito de forma a gerar o desenho automaticamente, não
foi possível aplicar as condições de periodicidade para se trabalhar com seções menores
da máquina para a simulação.
Na Figura (5-2) são mostradas as induções magnéticas encontradas na
simulação do gerador. Essa imagem ilustra a distribuição de indução e os valores médios
são os mostrados na Tabela (5-1).
66
Tabela 5-3 Comparação entre os valores das variáveis objetivo e de desempenho do GSPS
da geometria inicial versus a geometria encontrada na solução 1.
Sigla Valor Valor Valor encontrado Valor encontrado
mínimo máximo antes da otimização após a otimização
Objetivo Objetivo 51,81 [%] 98,03 [%]
FC Objetivo Objetivo 240,61 296,509
� 0,6 [T] 0,85 [T] 0,6292 [T] 0,7823 [T]
1,2 [T] 1,45 [T] 0,8049 [T] 1,41 [T]
1,2 [T] 1,55 [T] 0,8091 [T] 1,54 [T]
1 [T] 1,55 [T] 1,046 [T] 1,511 [T]
0,6 [T] 1,45 [T] 0,5152 [T] 1,44 [T]
� 0,8 [pu] 1,6 [pu] 3,77 [pu] 1,53 [pu]
� 3 [A/mm²] 5 [A/mm²] -- 4,5 [A/mm²]
� 3 [A/mm²] 5 [A/mm²] -- 4,46 [A/mm²]
Fonte: O autor.
5.2.2 Solução 2
Altura do condutor
ℎ 1 3
enrolamento rotor
Altura do condutor
1 4
enrolamento rotor
Fonte: O autor.
Na Tabela (5-5) são mostrados os valores das variáveis de desempenho
comparando a segunda solução com a máquina calculada no capítulo 2.
6 CONCLUSÃO
que ele possa apresentar respostas com maior fidelidade de resultados, que abranjam
melhor a forma construtiva real da máquina, é necessária a inclusão das barras de
amortecimento nos polos do rotor, realizar o cálculo da distorção harmônica e introduzir
a modelagem do efeito pelicular nos condutores da armadura.
chega ao valor informado pelo usuário. Desta forma, deve-se estudar critérios que
identifiquem quando o algoritmo alcançou a melhor solução.
79
REFERÊNCIAS
[7] BOLDEA, Ion; TUTELEA, Lucian. Electric machines : steady state, transients,
and design with MATLAB. 1. ed. New York: CRC Press, 2010. 800 p.
[8] BOLDEA, I. Synchronous Generators. Second. ed. London: CRC Press, 2016.
487p.
[13] CARLSON, R.; WURTZ, F. Optimal design of a set permanent magnet generators
with the same cross-section. In: 15º. SBMO. Simpósio Brasileiro de Microondas e
Optoeletrônica e 10º. CBMag - Congresso Brasileiro de Eletromagnetismo
(MOMAG 2012), 2012, João Pessoa/Paraíba. Anais de Congresso, 2012. v. 1. p. 1-6.
[21] MARTIGNONI, Alfonso. Máquinas de Corrente Alternada. 7. ed. Porto Alegre RS:
Editora Globo, 1920.
[23] MARTIN, Juan Corrales. Calculo Industrial de Máquinas Electricas: Tomo II.
[S.l.]: U.P.B, 1975.
[24] MEEKER, David. Finite Element Method Magnetics - Version 4.2 - User’s
Manual. 1. ed.[S.l.], 2015.
[26] Pedrosa Silva, Rodrigo & Li, Min & Rahman, Tanvir & Lowther, David. (2017).
Surrogate-Based MOEA/D for Electric Motor Design with Scarce Function Evaluations.
IEEE Transactions on Magnetics. PP. 1-1. 10.1109/TMAG.2017.2661987.
[33] SEARS, F.; ZEMANSKY, Mark W.; YOUNG, Hugh D. Física 3 – Eletricidade e
Magnetismo. 2. ed. [S.l.]: LTC Editora, 1984. 771 p.
Onde:
= Área do polo;
= Largura da base do polo;
= Coeficiente de empilhamento do rotor;
= Comprimento do pacote total do gerador;
bd 01 = D1 − b01
N1
bd 1 =
( bd 01 + bd11 + bd 21 )
3
Ad 1 =
N1
kemp L2bd 1
2p
Onde:
Ad 1 = Área do dente;
hc1 =
( D − ( D + 2h ) )
e1 1 n1
Onde:
Ac1 = Área da coroa do estator;
( D − 2(h ) − g + D
4. Área da coroa do rotor
+ hpk )
hc1 =
1
2
sch ic 2
i
2
Ac 2 = kemp hc 2 L2
Onde:
Ac 2 = Área da coroa do rotor;
kemp ( D1 − 2 gi ) Li
Agi =
2 pkc1
Onde:
Agi = Área do entreferro;
= − + ℎ
= + ℎ − −
Onde:
− −
= � ℎ ℎ
Onde:
= Massa do cobre do rotor;
= Largura média de uma espira do rotor;
� = Número de espiras por polo;
ℎ = Altura do condutor;
= Base do condutor;
ℎ = Quantidade de condutores em paralelo na altura da espira;
= Quantidade de condutores em paralelo na base da espira;
= Massa especifica do cobre;
= Número de pares de polos no rotor;
− −
= �ℎ ℎ
Onde:
= + +
−
=ℎ
−
=ℎ ℎ ℎ
� − � −
=
Onde:
= +
Onde:
%% Simulação
mi_addnode((-(Di3/2+Di2/2)/2),bpko/2);
mi_addnode((-(Di3/2+Di2/2)/2),-bpko/2);
mi_addnode((-(Di3/2+Di2/2)/2),bpko/2);%
mi_addnode((-(Di3/2+Di2/2)/2),-bpko/2);%
mi_addnode(-Di3/2-hpk,bpko/2);
mi_addnode(-Di3/2-hpk,-bpko/2);
mi_addnode(-Di3/2-hpk,bpko/2);
mi_addnode(-Di3/2-hpk,-bpko/2);
mi_addnode(-Di3/2-hpk,-bpko/2-dif);
mi_addnode(-Di3/2-hpk,bpko/2+dif);
mi_addsegment(-D1/2,-bpko/2-dif,-Di3/2-hpk,-bpko/2-dif);
mi_addsegment(-D1/2,+bpko/2+dif,-Di3/2-hpk,+bpko/2+dif);
mi_selectnode(-De2./2+Rp,-Rp); %Deletando nós que sobram
mi_selectnode(-De2./2+Rp,+Rp); %Deletando nós que sobram
mi_deleteselected;
mi_clearselected;
mi_drawrectangle(De1./1.5,De1./1.5,-De1./1.5,-De1./1.5);
mi_zoomnatural
mi_drawarc([-Di3./2 0;Di3./2 0], 180, 1); %Desenhando arco diâmetro externo do rotor
mi_selectarcsegment(Di3./2,0); %Copiando e rotacionando arco diâmetro externo do
rotor
mi_copyrotate2(0, 0, 180, 2, 3 );
mi_clearselected;
%
mi_selectnode(-Di3/2-hpk,bpko/2);
mi_setnodeprop('Polo',2);
mi_clearselected;
%
mi_selectnode(-Di3/2-hpk,-bpko/2);
mi_setnodeprop('Polo',2);
mi_clearselected;
92
%
mi_selectnode(-Di3/2-hpk,bpko/2);
mi_setnodeprop('Polo',2);
mi_clearselected;
%
mi_selectnode(-Di3/2-hpk,-bpko/2);
mi_setnodeprop('Polo',2);
mi_clearselected;
%
mi_selectnode(-Di3/2-hpk,-bpko/2-dif);
mi_setnodeprop('Polo',2);
mi_clearselected;
%
mi_selectnode(-Di3/2-hpk,bpko/2+dif);
mi_setnodeprop('Polo',2);
mi_clearselected;
%
mi_addsegment(-Di3/2-hpk,bpko/2,-Di3/2-hpk,bpko/2+dif);
mi_addsegment(-Di3/2-hpk,-bpko/2,-Di3/2-hpk,-bpko/2-dif);
%
mi_addsegment((-(Di3/2+Di2/2)/2),-bpko/2,-Di3/2-hpk,-bpko/2);
mi_addsegment((-(Di3/2+Di2/2)/2),bpko/2,-Di3/2-hpk,bpko/2);
%
mi_selectnode((-Di3/2-Di2/2)/2,bpko/2);%Deletando nós que sobram
mi_deleteselected;
mi_clearselected;
%
mi_selectnode((-Di3/2-Di2/2)/2,-bpko/2);%Deletando nós que sobram
mi_deleteselected;
mi_clearselected;
%
mi_selectnode(-D1/2,-bpko/2-dif);%Deletando nós que sobram
93
mi_deleteselected;
mi_clearselected;
%
mi_selectnode(-D1/2,+bpko/2+dif);%Deletando nós que sobram
mi_deleteselected;
mi_clearselected;
%
mi_selectsegment(-Di3/2-hpk,bpko/2); %Polo
mi_copyrotate2(0, 0,np, Pol, 1 );
mi_clearselected;
%
mi_selectsegment(-Di3/2-hpk,-bpko/2);%Polo
mi_copyrotate2(0, 0,np, Pol, 1 )
mi_clearselected
%
mi_selectsegment(-De2/2,bsch/2); %Polo
mi_copyrotate2(0, 0,np, Pol, 1 );
mi_clearselected;
%
mi_selectsegment(-De2/2,-bsch/2);%Polo
mi_copyrotate2(0, 0,np, Pol, 1 );
mi_clearselected;
%
mi_selectarcsegment(-De2/2,0);%Polo
mi_copyrotate2(0, 0,np, Pol, 3 );
mi_clearselected;
%
mi_selectsegment(-Di3/2,-bpko/2);%Polo
mi_copyrotate2(0, 0,np, Pol, 1 );
mi_clearselected;
%
94
mi_selectsegment(-Di3/2,bpko/2);%Polo
mi_copyrotate2(0, 0,np, Pol, 1 );
mi_clearselected;
%
mi_drawarc([-Di2./2 0;Di2./2 0], 180, 1);%Desenhando arco diâmetro interno do rotor
mi_selectarcsegment(Di2./2,0);%Copiando e rotacionando arco diâmetro interno do
rotor
mi_copyrotate2(0, 0, 180, 2, 3 );
mi_clearselected;
% P1
mi_addnode(hpk+Di3/2-4,bpko/2+0.4);
mi_addnode(hpk+Di3/2-4,bpko/2+0.4+Lfio);
mi_addnode(hpk+Di3/2-4-Hfio,bpko/2+0.4);
mi_addnode(hpk+Di3/2-4-Hfio,bpko/2+0.4+Lfio);
%
mi_addsegment(hpk+Di3/2-4,bpko/2+0.4,hpk+Di3/2-4,bpko/2+0.4+Lfio);
mi_addsegment(hpk+Di3/2-4-Hfio,bpko/2+0.4,hpk+Di3/2-4-Hfio,bpko/2+0.4+Lfio);
mi_addsegment(hpk+Di3/2-4,bpko/2+0.4+Lfio,hpk+Di3/2-4-Hfio,bpko/2+0.4+Lfio);
mi_addsegment(hpk+Di3/2-4-Hfio,bpko/2+0.4,hpk+Di3/2-4,bpko/2+0.4);
% P2
mi_addnode(hpk+Di3/2-4,-bpko/2-0.4);%
mi_addnode(hpk+Di3/2-4,-bpko/2+0.4-Lfio);%
mi_addnode(hpk+Di3/2-4-Hfio,-bpko/2-0.4);%
mi_addnode(hpk+Di3/2-4-Hfio,-bpko/2-0.4-Lfio);%
%
mi_addsegment(hpk+Di3/2-4,-bpko/2-0.4,hpk+Di3/2-4,-bpko/2-0.4-Lfio);
mi_addsegment(hpk+Di3/2-4-Hfio,-bpko/2-0.4,hpk+Di3/2-4-Hfio,-bpko/2-0.4-Lfio);
mi_addsegment(hpk+Di3/2-4,-bpko/2-0.4-Lfio,hpk+Di3/2-4-Hfio,-bpko/2-0.4-Lfio);
mi_addsegment(hpk+Di3/2-4-Hfio,-bpko/2-0.4,hpk+Di3/2-4,-bpko/2-0.4);
%
mi_selectsegment(hpk+Di3/2-4,bpko/2+0.4);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
95
mi_clearselected;
%
mi_selectsegment(hpk/2+Di3/2,bpko/2+0.4+Lfio);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
mi_selectsegment(hpk/2+Di3/2,bpko/2+0.4);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
mi_selectsegment(hpk+Di3/2-4-Hfio,bpko/2+0.4);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
mi_selectsegment(hpk+Di3/2-4,-bpko/2-0.4);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
%
mi_selectsegment(hpk/2+Di3/2,-bpko/2-0.4-Lfio);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
%
mi_selectsegment(hpk/2+Di3/2,-bpko/2-0.4);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
%
mi_selectsegment(hpk+Di3/2-4-Hfio,-bpko/2-0.4);%Bobina
mi_copyrotate2(0,0,np,Pol,1);
mi_clearselected;
%
mi_addnode(D1/2,b01/2);%Criando nó ranhura
96
mi_addnode(D1/2,-b01/2);%Criando nó ranhura
mi_addnode(D1/2+h31,b01/2);%Criando nó ranhura
mi_addnode(D1/2+h31,-b01/2);%Criando nó ranhura
%
mi_selectnode(D1/2,-b01/2);%Polo
mi_copyrotate2(0, 0,nran, N1, 0 );
mi_clearselected;
%
mi_addsegment(D1/2,b01/2,D1/2+h31,b01/2);
mi_addsegment(D1/2,-b01/2,D1/2+h31,-b01/2);
mi_addsegment(D1/2+h31,b01/2,D1/2+h31,-b01/2);
mi_addsegment(D1/2,b01/2,D1/2,2*b01);
%
mi_selectsegment(D1/2,-b01/2);
mi_copyrotate2(0, 0,nran, N1, 1 );
mi_clearselected;
%
mi_selectsegment(D1/2,b01/2);
mi_copyrotate2(0, 0,nran, N1, 1 );
mi_clearselected;
%
mi_selectsegment(D1/2,b01);
mi_copyrotate2(0, 0,nran, N1, 1 );
mi_clearselected;
%
mi_selectsegment(D1/2+h31,0);
mi_copyrotate2(0, 0,nran, N1, 1 );
mi_clearselected;
%
np1 = np*pi/180;
d=0;
while Pol > d
97
a= -Di3/2*cos(d*np1);
b= -Di3/2*sin(d*np1);
mi_selectarcsegment(a,b);
mi_deleteselected;
mi_clearselected;
d=d+1;
end
%
mi_selectnode(Di3/2,0);
mi_deleteselected;
mi_clearselected;
mi_selectnode(-Di3/2,0);
mi_deleteselected;
mi_clearselected;
%
iformula = (Iexct/Ncpp);
%
mi_addmaterial('Air', 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0);
mi_addmaterial('Coil', 1, 1, 0, 0, 58*0.65, 0, 0, 1, 0, 0, Wcue);
% mi_addmaterial('Coils', 1, 1, 0, 0, 58*0.65, 0, 0, 1, 0, 0, Wcu);
mi_addmaterial('LinearIron', 2100, 2100, 0, 0, 0, 0, 0, 0, 0, 0, 0);
mi_addcircprop('icoil',iformula, 1);
% mi_addcircprop('icoils',I/Ncp, 1);
mi_addmaterial('1020', 2100, 2100, 0, 0, 0, 0, 0, 0, 0, 0, 0);
bhcurve = [ 0,0.25,0.93,1.25,1.39,1.53,1.71,1.87,1.96,2.02,2.11,2.23,2.43;
0, 239, 796, 1592, 2387, 3979, 7958, 15916, 23873, 39789, 79578, 159155, 318310];
mi_addbhpoints('1020', bhcurve);
muo = pi*4.e-7;
%
%Adicionando materiais aos blocos
%
98
mi_copyrotate2(0, 0, np,1, 2 );
mi_clearselected;
%
mi_addblocklabel(Di3/2+hpk/3,-bpko/2-dif/3);
mi_selectlabel(Di3/2+hpk/3,-bpko/2-dif/3);
mi_setblockprop('Coil', 0, 10, 'icoil', 0, 0,-Z2);
mi_copyrotate2(0, 0, np, 1, 2 );
mi_clearselected;
%
mi_selectlabel(Di3/2+hpk/3,-bpko/2-dif/3);
mi_selectlabel(Di3/2+hpk/3,bpko/2+dif/3);
mi_deleteselected;
mi_clearselected;
%
mi_selectlabel(Di3/2+hpk/3,bpko+dif);
mi_copyrotate2(0, 0, np*2, Pol, 2 );
mi_clearselected;
%
mi_selectlabel(cos(1.2*np*pi/180)*Di3/2+hpk/3,sin(1.2*np*pi/180)*Di3/2+hpk/3);
mi_copyrotate2(0, 0, np*2, Pol, 2 );
mi_clearselected;
%
mi_addblocklabel(Di3/2,0);
mi_selectlabel(Di3/2,0);
mi_setblockprop('LinearIron', 0, 25, '<None>', 0, 0, 0);
mi_clearselected;
%
mi_addblocklabel(De1/2.05,0);
mi_selectlabel(De1/2.05,0);
mi_setblockprop('LinearIron', 0, 25, '<None>', 0, 0, 0);
mi_clearselected;
100
mi_savedxf('projeto.dxf')
mi_saveas('projeto.fem');
%
mi_analyze ;
mi_loadsolution;