XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Minimização do Fluxo de Potência Utilizando

Algoritmos Genéticos

Autores: Andrey R. da Silva UFSC / LVA

Arcanjo Lenzi UFSC / LVA
Edmar Baars EMBRACO
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Introdução
OBJETIVO : Implementar um algoritmo de otimização evolutiva (Algoritmo
Genético) dentro do ambiente ANSYS usando APDL

FOCO : Minimização do fluxo de potência

vibratória através do tubo de
descarga de um compressor
hermético
VANTAGEM : Poder executar todo o processo de
otimização dentro do ANSYS sem a
utilização de interface com outros
programas

Esquema de corte transversal de um compressor hermético

XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Metodologia
- Construção e validação experimental dos modelos em
EF do tubo e da carcaça do compressor.
- Implementação do A.G. no software ANSYS usando
APDL
- Otimização
- Análise dos resultados
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

O Algoritmo Genético
A idéia do Algoritmo Genético surgiu através do modelo darwiniano de
evolução das espécies.
O processo de otimização baseia-se em quatro etapas principais
(operadores básicos):

Geração
Aleatória
Escolha
(Fitness)
Cruzamento
(Cross-Over)
Mutação

Nova
Geração

Processo de otimização a aprtir do Algoritmo genético

XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

A Aplicação do AG na
Minimização do Fluxo de Potência

A Função objetivo
A função objetivo é média do fluxo de potência total dentro de cada faixa de frequência
analisada.
O fluxo de potência é determinado a partir dos resultados nodais (V e F) da análise harmônica.
O nó escolhido para a determinação do fluxo é aquele que se localiza na junta entre o tubo e a
carcaça.

BANDA N. de pontos

ANÁLISE 1 250 Hz 50

ANÁLISE 2 500 Hz 100

ANÁLISE 3 1000 Hz 200

ANÁLISE 4 2000 Hz 400

XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

A Aplicação do AG na
Minimização do Fluxo de Potência

As variáveis de Otimização
Como variáveis de otimização são Deslocamento
utilizados os keypoints que compõe o Prescrito (F)
modelo numérico do tubo. Estes
Keypoints podem variar no espaço
dentro de um limite de +- 2mm, com
passo de 1mm, nas três direções (x,y e
z).
O modelo do tubo possui ao todo 30
keypoints, sendo que os keypoints das
extremidades do tubo não variam.
Há, portanto, 84 variáveis no sistema. CARCAÇA
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

A Aplicação do AG na
Minimização do Fluxo de Potência
O Espaço de Otimização
Espaço de Otimização para 2 Variáveis

A figura ao lado ilustra o espaço

de otimização. Neste caso de 8

otimização consideram-se apenas 6

2 variáveis. 4

Função Otimização
2

Ou seja, variação de -2mm a 0

2mm de apenas 1 keypoint nas -2

-4
direções x e y, com passo de -6

1mm. 5

4 5
4.5
4
3
Neste caso teríamos um espaço 2
2
2.5
3
3.5

1.5
com 25 possibilidades de Espaço em Y [mm]
1 1
Espaço em X [mm]

combinação.
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

A Aplicação do AG na
Minimização do Fluxo de Potência
O Espaço de Otimização

Como, no caso do problema do tubo de descarga, existem 84

variáveis, cada qual podendo variar num espaço de 5 possibilidades,
há ao todo um espaço de otimização com 5E58 !
Caso fosse necessário varrer todo este espaço de otimização a
procura da melhor combinação (mínimo global), supondo que, para
cada combinação o computador demore 1 segundo para calcular o
fluxo de potência (suposição otimista), levaríamos 1,63E51 anos!!
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

6 – Fluxograma da Otimização
Geração :
gera 10 geometrias aleatórias

Análise :
Analise os F.P. em X,Y eZ
Cria um Vetor Média – Vetor Probabilidade
Médias e
Fitness : Mínimos
Escolhe 10 indivíduos
(10 sorteios na roleta)

Cross-over :
Binariza as geometrias dos escolhidos
Cruza genes aleatóriamente
Mutação :
Faz a inversão de um bit de um indivíduo
Aleatóriamente. Prob = 0,01
Transforma para hexa
Geometrias

Nova geração
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Resultados Obtidos
Otimização em banda estreita
Caso 1 : Banda de 250 Hz
10 individuos por geracao 6000 - 6250 Hz Faixa Otimizada : 6.0 - 6.25 kHz Df= 5
-80 -20

-30
-85
Menor fluxo de potencia medio de cada geracao dB W

-40
-90

-50

-95

Fluxo de Potencia dB W
-60

-100
-70

-105 -80

-90
-110

-100
-115
-110 tubo nao otimizado
faixa de otimizacao
-120 tubo otimizado
2 4 6 8 10 12 14 16 18 -120
Geracoes 5000 5500 6000 6500 7000 7500
Frequencia Hz
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

7 – Resultados Obtidos
Otimização em banda média
Caso 2 : Banda de 500 Hz
10 individuos por geracao 6000 - 6500 Hz
-65 Faixa Otimizada : 6.0 - 6.50 kHz Df= 5
-20

-30
-70
Menor fluxo de potencia medio de cada geracao dB W

-40
-75
-50

Fluxo de Potencia dB W
-80 -60

-70
-85

-80

-90
-90

-100
-95

-110 tubo nao otimizado

faixa de otimizacao
-100 tubo otimizado
2 4 6 8 10 12 14 16 18 20 -120
Geracoes 5000 5500 6000 6500 7000 7500
Frequencia Hz
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Resultados Obtidos
Otimização em banda larga
Caso 3 : Banda de 1000 Hz
10 individuos por geracao 6000 - 7000 Hz Faixa Otimizada : 6.0 - 7.0 kHz Df= 5
-60 -20

-65
-40
Menor fluxo de potencia medio de cada geracao dB W

-70

-60
-75

Fluxo de Potencia dB W
-80
-80

-85

-100
-90

-95 -120

-100
-140
-105 tubo nao otimizado
faixa de otimizacao
tubo otimizado
-110 -160
2 4 6 8 10 12 14 16 18 20 5000 5500 6000 6500 7000 7500 8000
Geracoes Frequencia Hz
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Resultados Obtidos
Otimização em banda larga
Caso 4 : Banda de 2000 Hz

10 individuos por geracao 6000 - 8000 Hz

-60 Faixa Otimizada : 6.0 - 8.0 kHz Df= 5
-20

-40
-70
Menor fluxo de potencia medio de cada geracao dB W

-60

-80

Fluxo de Potencia dB W
-80

-90
-100

-120
-100

-140

-110

-160
tubo nao otimizado
faixa de otimizacao
tubo otimizado
-120
5 10 15 20 25 -180
5000 5500 6000 6500 7000 7500 8000 8500 9000
Geracoes
Frequencia Hz
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Resultados Obtidos

Banda de Taxa de Redução média em

Otimização Otimização dB
250 Hz 55.75% 61.52
500 Hz 43.16% 41.51
1000 Hz 49.70% 52.27
2000 Hz 54.64% 60.25
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Geometrias Otimizadas

Banda de 250 Hz Banda de 500 Hz Banda de 1000 Hz

Banda de 2000 Hz Tubo Original

XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Conclusões
- É possível perceber que a minimização do fluxo de potência dentro da
banda de otimização é atingida graças a um aumento do fluxo de potência em
algumas regiões fora da banda otimizada.
- Acredita-se que a minimização dentro das faixas analisadas sejam atingidas
por dois mecanismos:
1) Os modos dentro da faixa de otimização são modificados de forma a
mover seus picos de ressonância para fora dessa região.
2) Alguns modos de vibração dentro da janela são alinhados de tal
maneira que se tenham interferências destrutivas.
XIII - Seminário de Elementos Finitos e Encontro de Usuários ANSYS

Bibliografia
BERANEK, Leo L. , Noise And Vibration Control Engineering. Principles And
Applications, JW&S INC. 1992.
ELLIOTT, S.J. and KEANE, A.J. Robustness of Optimal Design Solutions to
Reduce Vibration Transmission in a Lightweight 2-d Structure. JSV (2000)
229(3), 505 - 528
JONES, David L. G., Handbook of Viscoelastic Vibration Damping, JW&S
INC. 2001.
LENZI, Arcanjo, Apostila de Análise Estatistica Energética – Notas de Aula
MEAD, Denys J., Passive Vibration Control, JW&S INC. 1999.
.