UNIVERSIDADE FEDERAL DO PARÁ INSTITUTO DE TECNOLOGIA FACULDADE DE ENGENHARIA MECÂNICA ROBERTA TAMARA DA COSTA NERY / 03021003801 DESENVOLVIMENTO DE UMA FERRAMENTA VIRTUAL PARA ANÁLISE DE VIBRAÇÃO EM MÁQUINA ROTATIVA: APLICAÇÃO EM UMA BANCADA DIDÁTICA BELÉM 2008 ROBERTA TAMARA DA COSTA NERY / 03021003801 DESENVOLVIMENTO DE UMA FERRAMENTA VIRTUAL PARA ANÁLISE DE VIBRAÇÃO EM MÁQUINA ROTATIVA: APLICAÇÃO EM UMA BANCADA DIDÁTICA Trabalho de Conclusão de Curso apresentado à Faculdade de Engenharia Mecânica do Instituto de Tecnologia da Universidade Federal do Pará para obtenção do grau de Engenheiro Mecânico. Orientador: Prof. Dr. Newton Sure Soeiro BELÉM 2008 ROBERTA TAMARA DA COSTA NERY / 03021003801 DESENVOLVIMENTO DE UMA FERRAMENTA VIRTUAL PARA ANÁLISE DE VIBRAÇÃO EM MÁQUINA ROTATIVA: APLICAÇÃO EM UMA BANCADA DIDÁTICA Trabalho de Conclusão de Curso apresentado para obtenção do grau de Engenheiro Mecânico pela Universidade Federal do Pará. Submetido à banca examinadora do Colegiado constituída pelos MEMBROS: ____________________________________________ Prof. Dr. Newton Sure Soeiro (Orientador) ____________________________________________ Prof. Dr. Alexandre Luiz Amarante Mesquita (Membro) ____________________________________________ Prof. Walter dos Santos Sousa (Membro) ____________________________________________ Engª. Keliene Maria Sousa de Jesus (Convidada Externa) Julgado em: ____ de ______________ de 2008 Conceito: __________________ BELÉM 2008 Aos meus queridos pais, Jane da Costa Nery e Álvaro Celso de Oliveira Nery, que sempre me apoiaram. AGRADECIMENTOS Primeiramente, a Deus por estar sempre presente no meu coração, mostrar-me o melhor caminho e me confortar nos momentos difíceis. Aos meus queridos pais, Jane da Costa Nery e Álvaro Celso de Oliveira Nery, por todo amor, dedicação e esforço para que eu tivesse a oportunidade de estudar e realizar este sonho. Ao meu irmão Rodrigo Fabrício da Costa Nery pelo incentivo. Ao meu namorado e colega de turma Otávio Nascimento de Oliveira Neto pelo amor e companheirismo, sempre presente e disposto a me ajudar. Ao professor Newton Sure Soeiro, pela orientação neste trabalho e, principalmente, por ter me dado a oportunidade de fazer parte do Grupo de Vibrações e Acústica (GVA), onde pude complementar minha formação acadêmica. Aos meus amigos do GVA, professor Gustavo de Melo, Keliene Sousa, Diana Moraes, Juliana Vera, Alan do Vale, Carlos Cordeiro, Walter Sousa, Alexandre Sá, Márcio Mafra, Erlison Alves, Adriano da Silva, Fábio Setúbal e Helder Cardoso. Aos meus amigos, engenheiros eletricistas, Anderson Sena e Roger da Silva, sempre dispostos a me ensinar e esclarecer dúvidas sobre o Programa LabVIEW. Ao mecânico Odilon Oliveira Silva, do Laboratório de Engenharia Mecânica, fundamental na execução da parte prática deste trabalho. À Eletronorte pelo aprendizado e incentivo financeiro durante o período de projetos de pesquisa e por ter disponibilizado a bancada didática para este trabalho. Agradeço também ao professor Alexandre Luiz Amarante Mesquita, coordenador do projeto do qual esta bancada é produto. À banca examinadora, pelo convite aceito para avaliação deste trabalho. Não saber é ruim, Não querer saber é pior. PROVÉRBIO CHINÊS RESUMO A análise de vibração é um dos mais eficientes métodos de identificação de defeitos em máquinas, visto que ela se baseia no conhecimento das condições do equipamento, através de medições periódicas, de forma a evitar paradas inesperadas e substituição desnecessária de peças. Atualmente, as empresas estão utilizando cada vez mais instrumentação virtual, em seus departamentos de Engenharia de Manutenção. Uma ferramenta virtual é uma rotina desenvolvida em uma linguagem de programação, que simula um equipamento real no computador, pois apresenta botões, controles, LED’s, gráficos e demais funções. Em relação a máquinas rotativas, pode-se desenvolver várias metodologias de diagnóstico de defeitos, uma vez que as falhas mais comuns são desbalancemanto de massa, desalinhamento de eixos, falhas em rolamentos e em engrenagens. Em virtude disso, este trabalho tem como objetivo o desenvolvimento de uma ferramenta virtual que auxilie no diagnóstico de defeitos, que tem como base a plataforma LabVIEW. Esta ferramenta é composta por diversas rotinas responsáveis pela aquisição de sinais de vibração, arquivamento de dados e visualização de sinais. Além disso, foi elaborado um aplicativo que calcula as freqüências características de defeitos em rolamentos, considerando que estes são componentes mecânicos utilizados na maioria das máquinas rotativas e constituem as fontes mais comuns de falhas. Foram aplicadas também técnicas de análise de sinais, tais como: Nível Global, Fator de Crista, Curtose, Análise Espectral, Técnica do Envelope e Cepstrum. Palavras-chave: LabVIEW, Máquinas Rotativas, Manutenção Preditiva, Rolamentos, Análise de Sinais. LISTA DE FIGURAS Figura 3.1 Discrepância entre ER e EPI...................................................................... Figura 3.2 Efeito da força centrífuga no rotor............................................................. Figura 3.3 Desbalanceamento estático........................................................................ Figura 3.4 a) Binário de desbalanceamento; b) Binário de desbalanceamento causado por componente externo................................................................................ Figura 3.5 Desbalanceamento quase-estático............................................................. Figura 3.6 Desbalanceamento dinâmico..................................................................... Figura 3.7 Forças de reação nos mancais.................................................................... Figura 3.8 Sistema generalizado para alinhamento de eixos...................................... Figura 3.9 Sistema alinhado........................................................................................ Figura 3.10 Desalinhamento paralelo puro................................................................. Figura 3.11 Desalinhamento angular puro.................................................................. Figura 3.12 Desalinhamento misto............................................................................. Figura 3.13 Desalinhamento axial............................................................................... Figura 3.14 Acoplamento danificado por desalinhamento......................................... Figura 3.15 Tipos de defeitos mais comuns em rolamentos....................................... Figura 3.16 Pista interna danificada por desalinhamento........................................... Figura 3.17 Pista interna danificada por fadiga.......................................................... Figura 3.18 Rolamento na forma esquemática........................................................... Figura 4.1 Filtro ideal.................................................................................................. Figura 4.2 Filtro real................................................................................................... Figura 4.3 Filtro passa-baixa....................................................................................... Figura 4.4 Filtro passa-alta.......................................................................................... Figura 4.5 Filtro passa-banda...................................................................................... Figura 4.6 Filtro rejeita-banda..................................................................................... Figura 4.7 Amostragem de um sinal senoidal............................................................. Figura 4.8 Sinal que apresenta número inteiro de ciclos............................................ Figura 4.9 FFT do sinal que apresenta número inteiro de ciclos................................ Figura 4.10 Sinal que não apresenta número inteiro de ciclos.................................... Figura 4.11 FFT do sinal que não apresenta número inteiro de ciclos....................... Figura 4.12 Funções janela......................................................................................... Figura 4.13 FFT do sinal que não apresenta número inteiro de ciclos, com a janela Hanning....................................................................................................................... Figura 5.1 Curva de tendências................................................................................... Figura 5.2 Princípio do fator de crista......................................................................... Figura 5.3 Variação típica do fator de crista durante o processo de degradação de um rolamento.............................................................................................................. Figura 5.4 Esquema do redutor................................................................................... Figura 5.5 Simulação do sinal de vibração em um redutor......................................... Figura 5.6 Espectro do sinal....................................................................................... Figura 5.7 Modulação de amplitude............................................................................ Figura 5.8 Espectro de um sinal modulado................................................................. Figura 5.9 Esquema da técnica do envelope............................................................... Figura 5.10 Simulação de defeito em rolamento........................................................ Figura 5.11 Espectro do sinal..................................................................................... Figura 5.12 Sinal filtrado............................................................................................ Figura 5.13 Espectro do sinal filtrado........................................................................ 22 22 24 25 26 27 28 29 29 30 30 31 31 32 34 34 34 37 39 39 39 40 40 40 42 43 43 44 44 45 46 47 48 48 51 51 52 54 54 56 56 57 57 58 Figura 5.14 Envelope do sinal.................................................................................... Figura 5.15 Espectro do envelope.............................................................................. Figura 5.16 Sinal modulado........................................................................................ Figura 5.17 Espectro de potência do sinal.................................................................. Figura 5.18 Cepstrum do sinal.................................................................................... Figura 6.1 Exemplo de programação através de linhas de comando.......................... Figura 6.2 Exemplo de programação através de ícones............................................. Figura 6.3 Painel frontal.............................................................................................. Figura 6.4 Diagrama de bloco..................................................................................... Figura 6.5 Ícone.......................................................................................................... Figura 6.6 Conector.................................................................................................... Figura 6.7 Tela inicial do aplicativo........................................................................... Figura 6.8 Botões da tela inicial................................................................................. Figura 6.9 Tela de aquisição de sinais........................................................................ Figura 6.10 Controles para configuração de amostragem do sinal............................. Figura 6.11 Informação do diretório onde o arquivo deve ser salvo.......................... Figura 6.12 Exemplo de arquivo (.lvm) que contém os dados do sinal de vibração.. Figura 6.13 Configuração de filtro............................................................................. Figura 6.14 Configuração do canal de aquisição........................................................ Figura 6.15 Diagrama de bloco do módulo de aquisição de sinais............................ Figura 6.16 Tela de visualização de sinais em aceleração no domínio do tempo e da freqüência............................................................................................................... Figura 6.17 Tela de visualização de sinais em velocidade e deslocamento no domínio da freqüência................................................................................................. Figura 6.18 Configuração do sensor........................................................................... Figura 6.19 Controle que abre o arquivo.................................................................... Figura 6.20 Indicadores de valores de Pico, RMS, Fator de Crista e Curtose do sinal............................................................................................................................. Figura 6.21 Diagrama de bloco do módulo de visualização de sinais........................ Figura 6.22 Diagrama de bloco do módulo de visualização de sinais........................ Figura 6.23 Tela de cálculo de freqüências características de defeitos em rolamentos................................................................................................................... Figura 6.24 Diagrama de bloco do módulo de cálculo de freqüências características de defeitos em rolamentos................................................................... Figura 6.25 Tela da técnica do envelope..................................................................... Figura 6.26 Configuração de filtro............................................................................. Figura 6.27 Diagrama de bloco do módulo da técnica do envelope........................... Figura 6.28 Tela do Cepstrum.................................................................................... Figura 6.29 Diagrama de bloco do módulo do cepstrum........................................... Figura 6.30 Tela do nível global................................................................................. Figura 6.31 Tela do fator de crista.............................................................................. Figura 6.32 Tela da curtose......................................................................................... Figura 6.33 Diagrama de bloco que executa o nível global........................................ Figura 6.34 Diagrama de bloco que executa o fator de crista..................................... Figura 6.35 Diagrama de bloco que executa a curtose................................................ Figura 7.1 Esquema do sistema de aquisição de sinais de vibração........................... Figura 7.2 Bancada experimental................................................................................ Figura 7.3 Esquema do acelerômetro piezoelétrico.................................................... Figura 7.4 Acelerômetro piezoelétrico, modelo 784A................................................ Figura 7.5 Bancada experimental com sensores de vibração..................................... 58 58 60 60 60 61 62 63 63 63 63 64 65 66 66 67 67 68 68 68 69 69 70 70 71 71 72 73 73 74 74 75 79 80 81 81 82 82 83 83 84 85 87 87 88 Figura 7.6 Condicionador de sinais............................................................................ Figura 7.7 Placa de aquisição de sinais...................................................................... Figura 7.8 Massa desbalanceadora............................................................................. Figura 7.9 Disco rotativo com as massas desbalanceadoras...................................... Figura 7.10 Acoplamento desalinhado na bancada.................................................... Figura 7.11 Acoplamento desalinhado....................................................................... Figura 7.12 Acoplamento sem desalinhamento.......................................................... Figura 7.13 Rolamento utilizado para indução de defeitos........................................ Figura 8.1 Calibração do sistema............................................................................... Figura 8.2 Espectro em deslocamento / sinal sem desbalanceamento / posição axial............................................................................................................................. Figura 8.3 Espectro em deslocamento / sinal com desbalanceamento / posição axial............................................................................................................................. Figura 8.4 Espectro em deslocamento / sinal sem desbalanceamento / posição radial horizontal.......................................................................................................... Figura 8.5 Espectro em deslocamento / sinal com desbalanceamento / posição radial horizontal.......................................................................................................... Figura 8.6 Espectro em deslocamento / sinal sem desbalanceamento / posição radial vertical............................................................................................................... Figura 8.7 Espectro em deslocamento / sinal com desbalanceamento / posição radial vertical............................................................................................................... Figura 8.8 Espectro em deslocamento / sinal sem desalinhamento / posição axial.. Figura 8.9 Espectro em deslocamento / sinal com desalinhamento / posição axial.. Figura 8.10 Espectro em deslocamento / sinal sem desalinhamento / posição radial horizontal..................................................................................................................... Figura 8.11 Espectro em deslocamento / sinal com desalinhamento / posição radial horizontal..................................................................................................................... Figura 8.12 Espectro em deslocamento / sinal sem desalinhamento / posição radial vertical......................................................................................................................... Figura 8.13 Espectro em deslocamento / sinal com desalinhamento / posição vertical......................................................................................................................... Figura 8.14 Tela do aplicativo que calcula as freqüências características dos defeitos em rolamento................................................................................................. Figura 8.15 Domínio do tempo / posição axial / defeito na pista interna................... Figura 8.16 Domínio da freqüência / posição axial / defeito na pista interna............. Figura 8.17 Envelope / Domínio do tempo / posição axial / defeito na pista interna. Figura 8.18 Envelope / Domínio da freqüência / posição axial / defeito na pista interna.......................................................................................................................... Figura 8.19 Domínio do tempo / posição axial / defeito na pista externa................... Figura 8.20 Domínio da freqüência / posição axial / defeito na pista externa............ Figura 8.21 Envelope / Domínio do tempo / posição axial / defeito na pista externa Figura 8.22 Envelope / Domínio da freqüência / posição axial / defeito na pista externa......................................................................................................................... Figura A1 Domínio do tempo / sinal sem desbalanceamento / posição axial............. Figura A2 Domínio do tempo / sinal com desbalanceamento / posição axial............ Figura A3 Domínio do tempo / sinal sem desbalanceamento / posição radial horizontal..................................................................................................................... Figura A4 Domínio do tempo / sinal com desbalanceamento / posição radial horizontal..................................................................................................................... 89 89 90 90 91 81 92 92 94 95 95 96 96 96 97 98 98 98 99 99 99 100 101 102 102 102 103 104 104 104 109 109 110 110 Figura A5 Domínio do tempo / sinal sem desbalanceamento / posição radial vertical......................................................................................................................... Figura A6 Domínio do tempo / sinal com desbalanceamento / posição radial vertical......................................................................................................................... Figura A7 Espectro em aceleração / sinal sem desbalanceamento / posição axial..... Figura A8 Espectro em aceleração / sinal com desbalanceamento / posição axial..... Figura A9 Espectro em aceleração / sinal sem desbalanceamento / posição radial horizontal..................................................................................................................... Figura A10 Espectro em aceleração / sinal com desbalanceamento / posição radial horizontal..................................................................................................................... Figura A11 Espectro em aceleração / sinal sem desbalanceamento / posição radial vertical......................................................................................................................... Figura A12 Espectro em aceleração / sinal com desbalanceamento / posição radial vertical......................................................................................................................... Figura A13 Espectro em velocidade / sinal sem desbalanceamento / posição axial... Figura A14 Espectro em velocidade / sinal com desbalanceamento / posição axial.. Figura A15 Espectro em velocidade / sinal sem desbalanceamento / posição radial horizontal..................................................................................................................... Figura A16 Espectro em velocidade / sinal com desbalanceamento / posição radial horizontal..................................................................................................................... Figura A17 Espectro em velocidade / sinal sem desbalanceamento / posição radial vertical......................................................................................................................... Figura A18 Espectro em velocidade / sinal com desbalanceamento / posição radial vertical......................................................................................................................... Figura A19 Domínio do tempo / sinal sem desalinhamento / posição axial............... Figura A20 Domínio do tempo / sinal com desalinhamento / posição axial.............. Figura A21 Domínio do tempo / sinal sem desalinhamento / posição radial horizontal..................................................................................................................... Figura A22 Domínio do tempo / sinal com desalinhamento / posição radial horizontal..................................................................................................................... Figura A23 Domínio do tempo/ sinal sem desalinhamento / posição radial vertical.. Figura A24 Domínio do tempo / sinal com desalinhamento / posição vertical.......... Figura A25 Espectro em aceleração / sinal sem desalinhamento / posição axial....... Figura A26 Espectro em aceleração / sinal com desalinhamento / posição axial....... Figura A27 Espectro em aceleração / sinal sem desalinhamento / posição radial horizontal..................................................................................................................... Figura A28 Espectro em aceleração / sinal com desalinhamento / posição radial horizontal..................................................................................................................... Figura A29 Espectro em aceleração / sinal sem desalinhamento / posição radial vertical......................................................................................................................... Figura A30 Espectro em aceleração / sinal com desalinhamento / posição vertical... Figura A31 Espectro em velocidade / sinal sem desalinhamento / posição axial....... Figura A32 Espectro em velocidade / sinal com desalinhamento / posição axial...... Figura A33 Espectro em velocidade / sinal sem desalinhamento / posição radial horizontal..................................................................................................................... Figura A34 Espectro em velocidade / sinal com desalinhamento / posição radial horizontal..................................................................................................................... Figura A35 Espectro em velocidade / sinal sem desalinhamento / posição radial vertical......................................................................................................................... Figura A36 Espectro em velocidade / sinal com desalinhamento / posição vertical.. 111 111 111 112 112 112 113 113 114 114 114 115 115 115 116 116 117 117 117 118 118 118 119 119 120 120 120 121 121 121 122 122 SUMÁRIO 1 INTRODUÇÃO..................................................................................................... 1.1 Apresentação....................................................................................................... 1.2 Justificativa......................................................................................................... 1.3 Objetivos.............................................................................................................. 1.3.1 Objetivo geral.................................................................................................... 1.3.2 Objetivos específicos......................................................................................... 1.4 Metodologia utilizada......................................................................................... 1.5 Estrutura do trabalho........................................................................................ 14 14 15 15 15 16 16 17 2 REVISÃO BIBLIOGRÁFICA............................................................................. 2.1 Introdução........................................................................................................... 2.2 Defeitos em máquinas rotativas......................................................................... 2.3 Técnicas de análise de sinais.............................................................................. 2.4 Instrumentação virtual e LabVIEW.................................................................. 2.5 Considerações finais........................................................................................... 18 18 18 19 20 20 3 FUNDAMENTAÇÃO TEÓRICA SOBRE DEFEITOS MECÂNICOS MAIS COMUNS EM MÁQUINAS ROTATIVAS................................................ 3.1 Introdução........................................................................................................... 3.2 Desbalanceamento de massa.............................................................................. 3.2.1 Tipos de desbalanceamento............................................................................... 3.2.1.1 Desbalanceamento estático (Static Unbalance)............................................. 3.2.1.2 Binário de desbalanceamento (Couple Unbalance)....................................... 3.2.1.3 Desbalanceamento quase-estático (Quasi-Static Unbalance)...................... 3.2.1.4 Desbalanceamento dinâmico (Dynamic Unbalance)..................................... 3.2.2 Identificação do defeito de desbalanceamento através do sinal de vibração..... 3.3 Desalinhamento de eixos.................................................................................... 3.3.1 Alinhamento de eixos........................................................................................ 3.3.2 Tipos de desalinhamentos.................................................................................. 3.3.2.1 Desalinhamento paralelo puro....................................................................... 3.3.2.2 Desalinhamento angular puro........................................................................ 3.3.2.3 Desalinhamento misto.................................................................................... 3.3.3 Importância do alinhamento.............................................................................. 3.3.4 Identificação do defeito de desalinhamento através do sinal de vibração......... 3.4 Defeitos em rolamentos...................................................................................... 3.4.1 Identificação de defeitos em rolamentos através do sinal de vibração.............. 22 4 FUNDAMENTAÇÃO TEÓRICA SOBRE AQUISIÇÃO E ANÁLISE DE SINAIS....................................................................................................................... 4.1 Introdução........................................................................................................... 4.2 Filtros................................................................................................................... 4.2.1 Filtro ideal e filtro real....................................................................................... 4.2.2 Tipos de filtros segundo a banda passante......................................................... 4.3 Digitalização de um sinal analógico.................................................................. 4.3.1 Filtragem anti-aliasing....................................................................................... 4.4 Função Janela...................................................................................................... 22 22 24 24 25 26 27 28 29 29 30 31 31 32 33 33 34 36 39 39 39 39 40 42 42 43 5 FUNDAMENTAÇÃO TEÓRICA SOBRE AS TÉCNICAS DE ANÁLISE DE SINAIS PARA IDENTIFICAÇÃO DE DEFEITOS...................................... 5.1 Introdução........................................................................................................... 5.2 Técnicas de análise de sinais.............................................................................. 5.2.1 Nível global....................................................................................................... 5.2.2 Fator de crista (Crest-Factor)............................................................................ 5.2.3 Curtose............................................................................................................... 5.2.4 Análise espectral................................................................................................ 5.2.4.1 Análise espectral em aceleração, velocidade e deslocamento....................... 5.2.5 Técnica do envelope.......................................................................................... 5.2.5.1 Efeitos de modulação de amplitude em rolamentos....................................... 5.2.5.2 Demodulação de sinais para a detecção de defeitos em rolamentos (Envelope)................................................................................................................... 5.2.6 Cepstrum............................................................................................................ 6 DESCRIÇÃO DO APLICATIVO DESENVOLVIDO PARA AQUISIÇÃO E ANÁLISE DE SINAIS DE VIBRAÇÃO............................................................ 6.1 Introdução........................................................................................................... 6.1.1 Instrumentos virtuais........................................................................................ 6.2 Descrição do aplicativo desenvolvido para aquisição e análise de vibração..................................................................................................................... 6.2.1 Módulo de aquisição de sinais........................................................................... 6.2.2 Módulo de visualização de sinais...................................................................... 6.2.3 Módulo de cálculo de freqüências características de defeitos em rolamentos................................................................................................................... 6.2.4 Módulo da técnica do envelope......................................................................... 6.2.5 Módulo do cepstrum......................................................................................... 6.2.6 Módulo do nível global, fator de crista e curtose.............................................. 48 48 48 48 49 51 52 54 55 56 57 61 65 65 66 68 69 72 76 77 79 80 7 DESCRIÇÃO DA BANCADA EXPERIMENTAL............................................ 7.1 Introdução........................................................................................................... 7.2 Bancada experimental........................................................................................ 7.3 Cadeia de medição.............................................................................................. 7.3.1 Sensor de vibração............................................................................................. 7.3.1.1 Características do acelerômetro utilizado na bancada.................................. 7.3.1.2 Posicionamento dos acelerômetros na bancada............................................ 7.3.2 Condicionador de sinais..................................................................................... 7.3.3 Placa de aquisição de sinais............................................................................... 7.4 Indução de defeitos na bancada........................................................................ 7.4.1 Desbalanceamento de massa.............................................................................. 7.4.2 Desalinhamento de eixo..................................................................................... 7.4.3 Defeitos nos rolamentos.................................................................................... 84 84 85 86 86 87 88 88 89 90 90 91 92 8 MEDIÇÕES E ANÁLISE DOS SINAIS DE VIBRAÇÃO................................ 8.1 Introdução........................................................................................................... 8.2 Calibração do sistema de aquisição de sinais................................................... 8.3 Análise dos sinais de vibração........................................................................... 8.3.1 Análise do defeito de desbalanceamento........................................................... 8.3.2 Análise do defeito de desalinhamento............................................................... 8.3.3 Análise do defeito nos rolamentos..................................................................... 93 93 93 94 95 97 100 8.3.3.1 Análise do defeito na pista interna................................................................. 101 8.3.3.2 Análise do defeito na pista externa................................................................. 103 9 CONCLUSÕES E RECOMENDAÇÕES........................................................... 9.1 Conclusões........................................................................................................... 9.2 Sugestões para trabalhos futuros...................................................................... 106 106 106 REFERÊNCIAS........................................................................................................ 107 BIBLIOGRAFIA COMPLEMENTAR.................................................................. 108 APÊNDICE A............................................................................................................. 109 14 1 INTRODUÇÃO 1.1 Apresentação A vibração é um parâmetro muito utilizado em programas de manutenção preditiva de empresas de Engenharia, visto que ela indica a condição real do equipamento, evitando, portanto, paradas inesperadas e substituição desnecessária de peças. A análise do sinal de vibração pode ser realizada através de diversas técnicas, tais como: Nível Global, Fator de Crista, Curtose, Análise Espectral, Técnica do Envelope e Cepstrum, de acordo com o que se deseja diagnosticar. Por exemplo, se a necessidade é apenas monitorar o nível de vibração, pode-se utilizar o Nível Global. No entanto, se a necessidade é identificar o componente mecânico que apresenta defeito, pode-se realizar uma Análise Espectral. Defeitos como desbalanceamento de massa e desalinhamento de eixos podem ser identificados através da Análise Espectral. Porém, falhas em rolamentos requerem uma análise mais refinada, a qual pode ser realizada através da Técnica do Envelope. Este trabalho descreve um aplicativo computacional desenvolvido para análise de vibrações em máquinas rotativas, tendo por base a plataforma LabVIEW. Foram desenvolvidas diversas rotinas, que desempenham funções de aquisição e de visualização de sinais, bem como o cálculo de freqüências características de defeitos em rolamentos e, além disso, possibilitam o uso das técnicas de análise de sinais acima citadas. A validação do aplicativo foi realizada através de medições em uma bancada experimental, composta basicamente de um motor AC de ½ HP, eixo e discos rotativos. Nesta bancada foram induzidos os seguintes defeitos: desalinhamento do eixo, desbalanceamento de massa nos discos rotativos e defeitos nos diferentes elementos dos mancais de rolamento. 15 1.2 Justificativa O cenário organizacional sofreu muitas mudanças nos últimos anos. A concorrência entre as empresas está cada vez mais acirrada e o mercado tem exigido qualidade e agilidade nos processos. É necessário, portanto, que as máquinas estejam em condições operacionais o maior tempo possível. Desta forma, as empresas precisam adequar sua forma de atuação a fim de se manterem competitivas em seus mercados. Para isso, é interessante a implantação de um programa de manutenção preditiva eficaz. Departamentos de manutenção de empresas de Engenharia têm utilizado muito a instrumentação virtual, uma ferramenta que apresenta um custo inferior quando comparada a equipamentos convencionais de medição e mostra-se prática, pois com apenas um computador portátil com uma placa de aquisição de sinais e um aplicativo computacional pode ser transportada para qualquer lugar. Considerando que a análise de vibrações em máquinas já tem apresentado grandes contribuições para o entendimento do comportamento operacional de máquinas e a identificação de possíveis falhas, este trabalho tem como objetivo o desenvolvimento de uma ferramenta virtual para aquisição e análise de sinais de vibração, tendo por base a plataforma LabVIEW. 1.3 Objetivos 1.3.1 Objetivo geral Este trabalho tem o objetivo de desenvolver um aplicativo computacional para análise de vibrações em máquinas rotativas, tendo por base a plataforma LabVIEW. 16 1.3.2 Objetivos específicos Os objetivos específicos são: • Estudar as falhas mecânicas mais comuns em máquinas rotativas e seus efeitos de vibração; • Estudar as técnicas de análise de sinais utilizadas para o tratamento de dados e diagnóstico de defeitos; • Desenvolver um sistema aplicativo, tendo por base a plataforma LabVIEW, com módulo de aquisição de sinais, arquivamento de dados, cálculo de freqüências características de defeitos em rolamentos, visualização e tratamento de sinais; • Validar o aplicativo através de medições realizadas em uma bancada experimental composta, basicamente, por um motor AC de ½ HP, eixo e discos rotativos. 1.4 Metodologia utilizada A metodologia utilizada para a elaboração deste trabalho foi, primeiramente, estudar os defeitos mecânicos que ocorrem com mais freqüência em máquinas rotativas e alguns conceitos e técnicas de análise de sinais utilizadas para identificação destes defeitos. Depois foi realizado um estudo relativo ao software LabVIEW e suas ferramentas, utilizando os seus manuais, o “Help” do programa, apostilas, trabalhos de conclusão de curso e artigos. Em seguida, iniciou-se a programação em LabVIEW, na qual foram geradas rotinas para aquisição e visualização de sinais, cálculo de freqüências características de defeitos em rolamentos e aplicação de técnicas de análise de sinais. A atividade seguinte correspondeu à indução de defeitos na bancada de testes, tais como desbalanceamento de massa, desalinhamento de eixo e defeitos nos elementos dos mancais de rolamento. A medição dos sinais de vibração foi realizada sempre após a aplicação de cada defeito. A última etapa foi a análise dos sinais de vibração, sendo que foram aplicadas as técnicas de análise de sinais e realizadas as devidas conclusões. 17 1.5 Estrutura do trabalho O capítulo 2 apresenta a revisão bibliográfica que serviu de base teórica para elaboração deste trabalho. O capítulo 3 contém um estudo relativo à teoria de defeitos como desbalancemanto de massa, desalinhamento de eixos e defeitos em rolamentos. No capítulo 4 são apresentadas algumas informações importantes sobre a teoria de aquisição e análise de sinais. O capítulo 5 refere-se às técnicas de análise de sinais para identificação de defeitos. No capítulo 6 é apresentado o aplicativo computacional desenvolvido para aquisição e análise de sinais de vibração. O capítulo 7 corresponde à descrição da bancada de testes, na qual foram induzidos defeitos mecânicos e realizadas as medições de vibração. No capítulo 8 são mostrados os resultados das medições e as suas respectivas apreciações. Finalmente, o capítulo 9 apresenta as conclusões e recomendações para trabalhos futuros. 18 2 REVISÃO BIBLIOGRÁGICA 2.1 Introdução As máquinas rotativas constituem uma classe de equipamentos muito utilizada na indústria. Portanto, é de extrema importância um programa de manutenção preditiva eficiente, que busque eliminar as paradas desnecessárias destes equipamentos em um sistema. Para uma manutenção confiável, deve-se monitorar a condição de uma máquina para encontrar possíveis defeitos e diagnosticar o estado dos componentes críticos. Esta questão tem levado ao desenvolvimento de novas tecnologias nas áreas de sensores de medição, instrumentação virtual, softwares de processamento digital, simuladores de sistemas dinâmicos e análise de dados. Além disso, muitos estudos nestas áreas têm sido desenvolvidos com o objetivo de facilitar o monitoramento das máquinas e diagnosticar defeitos em tempo hábil para manutenção. 2.2 Defeitos em máquinas rotativas PACHOLOK (2004) em sua dissertação de mestrado mostra o desalinhamento de eixos como um dos defeitos mais comuns em máquinas rotativas. Além disso, apresenta a análise de vibração como uma das técnicas mais utilizadas pela manutenção preditiva e utiliza a termografia como uma ferramenta auxiliar. MARÇAL e SUSIN (2005) apresentam um método para detectar falha em máquinas rotativas baseado no padrão de vibração do sistema e diagnostica a condição de operação por Lógica Fuzzy. Eles consideram que quaisquer modificações ocorridas, devido ao desbalanceamento, ao desalinhamento, à falta de rigidez, ao acoplamento defeituoso, às correias frouxas ou gastas, aos eixos deformados e aos desajustes, alteram a amplitude da vibração na freqüência de rotação do sistema rotativo em questão. 19 2.3 Técnicas de análise de sinais NUNES (1989), em seu trabalho utiliza a Técnica do Envelope para diagnosticar defeitos em mancais de rolamentos, através da análise de vibrações. Ele mostra que esta técnica é eficiente, sensível e segura na detecção dos picos de freqüência relacionados aos defeitos nos rolamentos. ANTONIOLLI (1999) realizou uma análise comparativa entre as técnicas de medição e análise de sinais de vibração, aplicadas na manutenção preditiva: Fator de Crista, Análise Espectral e Envelope. Para isso, simulou defeitos em rolamentos em uma bancada de testes. MESQUITA et al. (2002) desenvolveram um trabalho no qual detecta falhas em rolamentos usando uma técnica de transformada tempo-freqüência, comparando com a Técnica do Envelope. BEZERRA (2004) desenvolveu um estudo para detectar falhas em rolamentos por análise de vibração, através da aplicação de diversas técnicas, tanto as que ocorrem no domínio do tempo (Média absoluta, Nível global RMS, Fator de crista e Curtose), quanto as que ocorrem no domínio da freqüência (Densidade espectral de potência, Envelope e Cepstrum). As técnicas foram aplicadas a rolamentos com falhas induzidas nas pistas e esfera com diferentes tamanhos de falhas e submetidos a diferentes velocidades, para uma melhor compreensão das técnicas. SANTOS (2006) utilizou a FRF como auxílio na Técnica do Envelope para identificar defeitos nos rolamentos. O trabalho foi realizado tendo por base um aplicativo desenvolvido no software MATLAB, o qual foi validado com a utilização de um moto-esmeril. PONCI e CUNHA explicam as principais causas de defeitos em rolamentos, quais as freqüências características e como elas se comportam no sinal de vibração, tais como os seus efeitos de modulação. Para diagnóstico dos defeitos, aplicam a Análise Espectral e Técnica do Envelope. Além disso, mostram a variação do Fator de Crista no processo de degradação de um rolamento. 20 2.4 Instrumentação virtual e LabVIEW NATIONAL INSTRUMENTS (2000) apresenta em seus manuais do LabVIEW como criar os Instrumentos Virtuais (VI’s), os quais correspondem às rotinas desenvolvidas na Linguagem G (Linguagem de Programação Gráfica). Estes VI’s podem ser simples programas, como calcular operações básicas, quanto trabalhar com a aquisição e tratamento de sinais. AQUINO (2004) explana detalhadamente uma técnica de monitoramento baseada na função densidade de probabilidade Beta que analisa o engrenamento, dente a dente, sendo que a detecção da anomalia é feita através de um parâmetro estatístico que é extraído do sinal de vibração, visando facilitar o diagnóstico de defeitos de montagem em redutores de engrenagens. SOUSA (2005) desenvolveu um sistema aplicativo em LabVIEW para o monitoramento de máquinas rotativas com um módulo de balanceamento de rotores. Em seu trabalho, ele apresentou a teoria relativa ao Desbalanceamento Rotativo e os seus efeitos de vibração, bem como Instrumentação Virtual. SOUZA (2007) desenvolveu um projeto para implantação de uma bancada de simulação de carga composta de uma máquina primária (motor CC) e uma máquina de teste (gerador síncrono). Através deste projeto, é enfocado o estudo de vários elementos que são implantados na bancada para uma completa monitoração e diagnóstico de todas as suas variáveis, tais como: sensores, placa de aquisição de dados, acionamento eletrônico, circuito condicionadores de sinais e as máquinas elétricas girantes. 2.5 Considerações finais Para o desenvolvimento de um aplicativo computacional com funções de aquisição e análise de sinais de vibração, com o objetivo de se identificar defeitos em máquinas rotativas é necessário: • Conhecer os defeitos mais freqüentes neste tipo de equipamento e de que forma eles se apresentam; 21 • Ter conhecimento das técnicas utilizadas para análise de sinais de vibração que possibilitam o acompanhamento de uma anomalia e a identificação de defeitos antes que a máquina venha a parar; • Escolher uma linguagem de programação que facilite o desenvolvimento do aplicativo e que este seja de simples utilização por parte dos usuários. Dentre os trabalhos coletados, estes foram os que mais contribuíram para o desenvolvimento deste trabalho. No entanto, vale ressaltar que há uma vasta literatura para cada tópico abordado. 22 3 FUNDAMENTAÇÃO TEÓRICA SOBRE DEFEITOS MECÂNICOS MAIS COMUNS EM MÁQUINAS ROTATIVAS 3.1 Introdução Uma fonte comum de vibrações em máquinas rotativas é o desbalanceamento de massa, gerado por assimetrias, tolerâncias dimensionais, desvios de forma, imperfeições de matéria-prima e de montagem. Qualquer uma destas causas ou uma combinação delas acabará com a condição de perfeita distribuição de massa em torno do eixo de rotação do rotor. Outra origem comum de vibrações é o desalinhamento de eixos. Acoplamentos e mancais são os primeiros componentes a sentir os efeitos do desalinhamento, pois as cargas geradas causam desgaste prematuro nestes componentes. As vibrações e empenos nos eixos causados pelo desalinhamento afetam diretamente selagens axiais e radiais. Em engrenagens, o desalinhamento pode causar ruídos por mau engrenamento e desgastes prematuros nos dentes. Por este motivo, demonstra-se a importância dada ao alinhamento no caso de acoplamento de engrenagem. Em relação aos eixos, dependendo da robustez, podem sofrer empenos, atritos com peças estacionárias ou até mesmo fraturar por fadiga quando ficam submetidos às cargas elevadas provenientes do desalinhamento. Os rolamentos, que são componentes presentes na maioria dos equipamentos rotativos, apresentam um grande índice de falhas e, muitas vezes, prematuras. As causas mais comuns dos defeitos em rolamentos são: seleção incorreta; sobrecarga; defeito de fabricação; desalinhamento; montagem incorreta; estocagem inadequada; lubrificação inadequada, excessiva ou insuficiente; falha de vedação e descargas elétricas através dos mancais. 3.2 Desbalanceamento de massa O desbalanceamento constitui uma das fontes mais comuns nos problemas de vibração em máquinas. É um defeito gerado pela distribuição irregular de massa em torno do eixo de rotação de um rotor, ou seja, o eixo principal de inércia (EPI) não coincide com o eixo de rotação (ER), como pode ser visto na figura 3.1. 23 Figura 3.1 Discrepância entre ER e EPI Fonte: ALMEIDA e GÓZ O excesso de massa em um lado do rotor faz com que a força centrífuga atuante sobre este lado mais pesado supere a força centrífuga atuante sobre o lado oposto, forçando, então, o lado mais leve na direção do lado mais pesado, sendo a força resultante a causadora da vibração. A magnitude desta força depende da velocidade de rotação da máquina e da excentricidade do centro de gravidade do rotor. A força centrífuga é dada pela equação (3.1): Fcent = m.e.ω 2 (3.1) em que m é a massa desbalanceadora, e é a excentricidade ou a distância do centro de gravidade da massa ao eixo de giro do rotor e ω (rad/s) a rotação da máquina. Figura 3.2 Efeito da força centrífuga no rotor 24 3.2.1 Tipos de desbalanceamento Existem quatro tipos de desbalanceamento. Para facilitar a definição de cada um dos tipos de desbalanceamento, será considerado um rotor perfeitamente balanceado, onde se supõe a adição de massas de desbalanceamento. 3.2.1.1 Desbalanceamento estático (Static Unbalance) O desbalanceamento estático surge quando o eixo principal de inércia de um rotor encontra-se deslocado, no entanto paralelo em relação ao eixo de rotação do mesmo, como mostrado na figura 3.3. A distância entre o centro de gravidade (CG) e o eixo de rotação (a excentricidade e) origina a força centrífuga quando da rotação do rotor. No caso do desbalanceamento estático a massa de desbalanceamento e o CG encontram-se no mesmo plano normal ao eixo de rotação, conforme figura 3.3a. Um rotor com duas massas de desbalanceamento iguais e eqüidistantes com relação ao CG (figura 3.3b) caracteriza também um desbalanceamento estático, uma vez que o efeito de ambas é equivalente ao efeito de uma massa localizada no plano do CG. A eliminação do desbalanceamento estático é obtida movendo-se o centro de gravidade em direção ao eixo de rotação do rotor. Isto é feito através da adição ou retirada de massa do rotor, de forma que a força radial causada por esta adição ou retirada seja igual em módulo à força causada pela excentricidade do CG, mas com direção oposta. Rotores desbalanceados estaticamente podem ser submetidos a um balanceamento em um plano. 25 Figura 3.3 Desbalanceamento estático Fonte: SOUSA, 2005 3.2.1.2 Binário de desbalanceamento (Couple Unbalance) O binário de desbalanceamento ocorre quando o eixo principal de inércia corta o eixo de rotação, sendo o CG o ponto comum entre eles. Neste caso a presença das massas de desbalanceamento não desloca o CG com relação ao eixo de rotação do rotor, no entanto causa a inclinação do eixo principal de inércia. As massas de desbalanceamento geram duas forças centrífugas (Fcent1 e Fcent2), mostradas na figura 3.4. Estas forças centrífugas formam um binário, portanto se anulam, por serem paralelas e possuírem sentidos opostos. No entanto, geram um momento, já que não estão contidas no mesmo eixo. O exemplo mais simples de binário de desbalanceamento é o de massas de desbalanceamento iguais posicionadas nos extremos opostos de um rotor e deslocadas de 180º uma da outra, como mostra a figura 3.4a. Em um binário de desbalanceamento, o CG pode estar inclusive fora do espaço compreendido entre tais planos, como mostra o exemplo da figura 3.4b, onde o binário de desbalanceamento é causado por um componente externo ao rotor. Quanto à rotação de um rotor com desbalanceamento binário, cada uma das extremidades vibra em direções opostas. Este tipo de desbalanceamento não pode ser 26 corrigido com apenas uma massa. Neste caso são necessárias duas massas, cada uma devendo anular o efeito de uma das massas de desbalanceamento. Faz-se, portanto, necessário um outro binário com massas opostas de 180º, possibilitando o “giro” do eixo principal de inércia sobre o CG, fazendo com que aquele coincida com o eixo de rotação. O balanceamento deste tipo de rotor é chamado de balanceamento em dois planos ou balanceamento dinâmico. Figura 3.4 a) Binário de desbalanceamento; b) Binário de desbalanceamento causado por componente externo Fonte: SOUSA, 2005 3.2.1.3 Desbalanceamento quase-estático (Quasi-Static Unbalance) O desbalanceamento quase-estático ocorre quando o eixo principal de inércia corta o eixo de rotação do rotor em um ponto, que não passa pelo CG, como pode ser visto na figura 3.5a. Este tipo de desbalanceamento também pode ser originado pela combinação dos desbalanceamentos estático e binário, conforme mostrado na figura 3.5b. A figura 3.5c mostra o desbalanceamento quase-estático causado por um componente externo acoplado. O balanceamento de um rotor deste tipo só é possível se realizado em dois planos. 27 Figura 3.5 Desbalanceamento quase-estático Fonte: SOUSA, 2005 3.2.1.4 Desbalanceamento dinâmico (Dynamic Unbalance) O desbalanceamento dinâmico ocorre quando o eixo principal de inércia não encontra o eixo de rotação. Este tipo de desbalanceamento é o mais freqüente. As massas de desbalanceamento, neste caso, encontram-se deslocadas de um ângulo diferente de 180º, como apresentado nas figuras 3.6a e 3.6b. A figura 3.6b mostra um binário de desbalanceamento combinado com um desbalanceamento estático deslocados de um ângulo diferente de 180º, o que resulta em um desbalanceamento dinâmico. Como nos dois tipos 28 anteriores de desbalanceamentos, o desbalanceamento dinâmico também só pode ser corrigido com massas de correção em pelo menos dois planos perpendiculares ao eixo de rotação. Figura 3.6 Desbalanceamento dinâmico Fonte: SOUSA, 2005 3.2.2 Identificação do defeito de desbalanceamento através do sinal de vibração O defeito de desbalanceamento pode ser identificado através da análise espectral e caracteriza-se por ocorrer na freqüência de rotação da máquina. Deve-se medir a vibração da máquina com o sensor nas posições radial e axial. Se o sinal referente à posição radial apresentar maior amplitude, na freqüência de rotação da máquina, o defeito é desbalanceamento. Isto se explica devido ao fato da força centrífuga, gerada pelo desbalanceamento, ocorrer na posição radial. 29 3.3 Desalinhamento de eixos O desalinhamento é a condição em que os eixos da máquina motriz e da máquina motora não estão localizados na mesma linha de centro. Segundo PACHOLOK (2004), o desalinhamento de eixos de máquinas dinâmicas provoca forças de reação nos mancais, conforme ilustrado na figura 3.7. A alta carga nos mancais tem como conseqüência o aumento do consumo de energia fornecida pelo acionador. Figura 3.7 Forças de reação nos mancais Fonte: PACHOLOK, 2004 3.3.1 Alinhamento de eixos O alinhamento é o processo pelo qual dois eixos são posicionados de forma que suas linhas de centro fiquem colineares quando em operação. Para a conceituação de alinhamento, considera-se a posição relativa dos eixos em dois planos imaginários, horizontal e vertical, conforme a figura 3.8. Denomina-se um dos eixos como “acionador” e o outro como “acionado”. Considera-se que um conjunto de máquinas acionador-acionado está alinhado quando o desvio de concentricidade dos eixos, em todo o seu comprimento e com uma revolução completa, está com valores de desvio dentro dos limites de tolerância. A figura 3.9 ilustra um sistema perfeitamente alinhado, o qual deve ser considerado nos planos horizontal e vertical. 30 Figura 3.8 Sistema generalizado para alinhamento de eixos Fonte: PACHOLOK, 2004 Figura 3.9 Sistema alinhado Fonte: PACHOLOK, 2004 3.3.2 Tipos de desalinhamentos O desalinhamento pode ocorrer de três formas: 31 3.3.2.1 Desalinhamento paralelo puro O desalinhamento paralelo puro ocorre quando suas linhas de centro estão paralelas entre si, porém não coincidentes, conforme a figura 3.10. Figura 3.10 Desalinhamento paralelo puro Fonte: PACHOLOK, 2004 3.3.2.2 Desalinhamento angular puro O desalinhamento angular puro ocorre quando as linhas de centro dos eixos formam um ângulo entre si, conforme a figura 3.11. Figura 3.11 Desalinhamento angular puro Fonte: PACHOLOK, 2004 32 3.3.2.3 Desalinhamento misto O desalinhamento misto, também chamado de desalinhamento combinado, ocorre quando existe a associação dos dois desalinhamentos anteriores, ou seja, as linhas de centro dos eixos não estão coplanares e formam um ângulo entre si, conforme a figura 3.12. É o desalinhamento mais encontrado na prática. Figura 3.12 Desalinhamento misto Fonte: PACHOLOK, 2004 Também é importante considerar o desalinhamento axial dos eixos, conforme a figura 3.13. Esse tipo de desalinhamento corresponde ao desvio de posição axial dos eixos em relação a uma posição de referência. Por exemplo, um dos eixos está deslocado para trás. Figura 3.13 Desalinhamento axial Fonte: PACHOLOK, 2004 33 3.3.3 Importância do alinhamento O alinhamento de eixos de equipamentos dinâmicos é muito importante para o aumento do tempo médio entre uma falha e outra, ou seja, quebra menos freqüente dos equipamentos, principalmente de mancais e acoplamentos. E, conseqüentemente, baixos custos de manutenção e aumento da confiabilidade dos equipamentos, garantindo, desta forma, a qualidade dos processos. A figura 3.14 mostra um acoplamento danificado por desalinhamento. Figura 3.14 Acoplamento danificado por desalinhamento Fonte: PACHOLOK, 2004 3.3.4 Identificação do defeito de desalinhamento através do sinal de vibração O defeito de desalinhamento, assim como o desbalanceamento, pode ser identificado através da análise espectral, apresentando as seguintes características: • A freqüência da vibração é normalmente igual a 1x rotação. No entanto, se o desalinhamento for severo surgem também em 2x rotação e 3x rotação; • A amplitude é proporcional à quantidade de desalinhamento; • A amplitude de vibração pode ser alta na direção axial bem como na radial. O desalinhamento, mesmo com acoplamentos flexíveis, produz forças axiais e radiais 34 que, por sua vez produzem vibrações radiais e axiais. Sempre que a amplitude da vibração axial for maior que a metade da maior amplitude radial, deve-se suspeitar de desalinhamento ou empenamento. • O desalinhamento angular submete os eixos à vibração axial na freqüência 1x rotação. O desalinhamento paralelo produz uma vibração radial em uma freqüência de 2x rotação. E o desalinhamento combinado, além da vibração predominante acontecer na direção axial em 1x rotação, ocorre uma vibração significativa em 2x rotação nesta direção. 3.4 Defeitos em rolamentos Os defeitos em rolamentos podem ocorrer nos seus diferentes componentes. Geralmente, estes defeitos evoluem com certa lentidão e emitem sinais com bastante antecedência da falha final, que pode ocorrer por travamento ou ruptura dos componentes. Defeitos típicos que evoluem dessa forma são: riscos nas pistas, roletes ou esferas, trincas, corrosão, erosão e contaminação. O processo de degradação de um rolamento pode se iniciar na pista externa ou interna, num dos elementos rolantes (rolos ou esferas) ou na gaiola, alastrando-se depois para os demais componentes. A figura 3.15 apresenta os tipos de defeitos mais comuns em rolamentos. A figura 3.16 mostra uma pista interna danificada por desalinhamento. E a figura 3.17 refere-se a uma pista interna danificada por fadiga. 35 Figura 3.15 Tipos de defeitos mais comuns em rolamentos Fonte: NUNES, 1989 Figura 3.16 Pista interna danificada por desalinhamento Fonte: PACHOLOK, 2004 Figura 3.17 Pista interna danificada por fadiga Fonte: ANTONIOLLI, 1999 36 3.4.1 Identificação de defeitos em rolamentos através do sinal de vibração Os defeitos em rolamentos podem ser identificados pela presença de freqüências características. Ao contrário da maioria das freqüências de vibração geradas por componentes mecânicos, essas freqüências são verdadeiramente freqüências de defeito. Isto é, elas só estarão presentes nos espectros de vibração quando os rolamentos estiverem realmente defeituosos ou, pelo menos, quando seus componentes estiverem sujeitos a tensões e deformações excessivas que poderão induzir uma falha. Na maioria dos casos a pista externa é fixa e a pista interna gira junto com o eixo de rotação da máquina. Quando isto acontece, as freqüências podem ser calculadas através das seguintes equações: • fe = • fi = • Defeito na pista externa: n D cos β f r 1 − 2 dm (3.2) Defeito na pista interna: n D f r 1 + cos β 2 dm (3.3) Defeito nos elementos rolantes: 2 dm D fb = f r 1 − cos β D d m (3.4) 37 • fg = Defeito na gaiola: 1 D f r 1 − cos β 2 dm sendo: f r : Freqüência de rotação em Hz; n : Número de elementos rolantes; D : Diâmetro dos elementos rolantes (mm); d m : Diâmetro primitivo (mm); β : Ângulo de contato; d o : Diâmetro da pista externa; d i : Diâmetro da pista interna; Do : Diâmetro externo; Di : Diâmetro interno. A figura 3.18 ilustra um rolamento na forma esquemática: (3.5) 38 Figura 3.18 Rolamento na forma esquemática Fonte: ANTONIOLLI, 1999 39 4 FUNDAMENTAÇÃO TEÓRICA SOBRE AQUISIÇÃO E ANÁLISE DE SINAIS 4.1 Introdução Este capítulo apresenta algumas informações importantes sobre a teoria de aquisição e análise de sinais: Filtros, Digitalização de um sinal analógico e Função Janela. 4.2 Filtros Segundo MESQUITA (2008), um filtro tem a função de remover partes não desejadas do sinal, como o ruído, ou extrair partes úteis do sinal, como determinadas componentes de freqüência que estão dentro do gama de freqüência. 4.2.1 Filtro ideal e filtro real De acordo com MESQUITA (2008), um filtro ideal é aquele sistema, cuja resposta em freqüência é unitária dentro de certa banda de freqüência e exatamente zero para outras bandas, sem haver atenuação, como mostra a figura 4.1. E em um filtro real, na sua resposta em freqüência, há uma atenuação em certas freqüências e também há uma oscilação na banda passante chamada “ripple”. Geralmente a freqüência de corte é definida após um decaimento de 3 dB na resposta em freqüência, conforme a figura 4.2. 40 Figura 4.1 Filtro ideal Fonte: B&K, 1998 Figura 4.2 Filtro real Fonte: B&K, 1998 4.2.2 Tipos de filtros segundo a banda passante • Filtro passa-baixa: permite a passagem de freqüências abaixo de sua freqüência de corte. Figura 4.3 Filtro passa-baixa Fonte: MESQUITA, 2008 41 • Filtro passa-alta: permite a passagem de freqüências acima de sua freqüência de corte. Figura 4.4 Filtro passa-alta Fonte: MESQUITA, 2008 • Filtro passa-banda: permite a passagem apenas de uma faixa de freqüências. Figura 4.5 Filtro passa-banda Fonte: MESQUITA, 2008 • Filtro rejeita-banda: permite a passagem de todas as freqüências fora de uma certa faixa. Figura 4.6 Filtro rejeita-banda Fonte: MESQUITA, 2008 42 4.3 Digitalização de um sinal analógico Para digitalizar um sinal analógico são necessárias, no mínimo, quatro etapas: • Filtragem anti-aliasing (Antialias filtering); • Amostragem (Sampling); • Quantificação (Quantization); • Codificação (Encoding); 4.3.1 Filtragem anti-aliasing De acordo com o Teorema de Nyquist, a taxa ou freqüência de amostragem (número de amostras por unidade de tempo de um sinal) deve ser maior que o dobro da maior freqüência contida no sinal a ser amostrado, para que possa ser reproduzido integralmente sem erro de aliasing. Esta relação é dada pela equação (4.1): f0 ≥ 2 fm (4.1) onde: f 0 : Freqüência de amostragem f m : Maior freqüência contida no sinal A metade da freqüência de amostragem é chamada freqüência de Nyquist e corresponde ao limite máximo de freqüência do sinal que pode ser reproduzido. Como não é possível garantir que o sinal não contenha sinais acima deste limite (distorções, interferências, ruídos e etc), é necessário filtrar o sinal com um filtro passa-baixa, com freqüência de corte igual (ou menor) à freqüência de Nyquist, ou filtro anti-aliasing. A figura 4.7 mostra um sinal senoidal sendo amostrado com taxas próximas ao limite. A figura 4.7a corresponde à amostragem com freqüência maior que duas vezes a do sinal: há amostras suficientes para que o sinal possa ser reproduzido sem erro de aliasing. Na figura 4.7b, a taxa de amostragem é igual a duas vezes a freqüência do sinal: não é possível a sua reprodução. Na figura 4.7c, a freqüência de amostragem é menor que o dobro da freqüência 43 do sinal: a quantidade de amostras é insuficiente e o sinal reproduzido está errado. Este erro é causado pelo fenômeno de aliasing. Figura 4.7 Amostragem de um sinal senoidal 4.4 Função Janela Durante a conversão do sinal do domínio do tempo para o domínio da freqüência, usando a FFT, ocorre o efeito de vazamento espectral (leakage). O leakage corresponde a um fenômeno que tende a espalhar a energia contida numa freqüência em linhas de freqüências adjacentes no espectro, distorcendo-o, e fazendo com que a amplitude seja subestimada. Segundo MESQUITA (2008), o processo da FFT requer que o sinal amostrado consista de uma representação completa do sinal original no domínio do tempo ou contenha uma repetição periódica do sinal original medido. Em sinais transientes, cuja amplitude decai a zero, totalmente contidos no registro temporal, isto não implica em problemas. Porém, em sinais estacionários, como, por exemplo, sinais periódicos, pode haver problemas. O sinal periódico pode não ter um número inteiro de ciclos capturado pelo registro temporal, fazendo 44 com que a repetição do registro temporal não reconstrua o sinal original exatamente, e implicando em distorção no seu espectro. A figura 4.8 exibe um sinal senoidal que apresenta um número inteiro de ciclos. A figura 4.9 representa a sua FFT. Figura 4.8 Sinal que apresenta número inteiro de ciclos Figura 4.9 FFT do sinal que apresenta número inteiro de ciclos A figura 4.10 exibe um sinal senoidal que não apresenta um número inteiro de ciclos. A figura 4.11 representa a sua FFT, com efeito de leakage. 45 Figura 4.10 Sinal que não apresenta número inteiro de ciclos Figura 4.11 FFT do sinal que não apresenta número inteiro de ciclos Para minimizar este efeito, os sinais não periódicos e os periódicos, que não possuam número inteiro de ciclos dentro do registro temporal, devem ser multiplicados por uma função matemática que faz com que sinal tenha valor igual a zero no início e no final do registro temporal. Esta função matemática é chamada de função janela ("window function"). Janelas comumente usadas pelos analisadores digitais são: a janela uniforme (retangular), janela Hanning, Kaiser-Besel e Flat Top, mostradas na figura 4.10. 46 Figura 4.12 Funções janela no domínio do tempo Fonte: MCCONNELL, 1995 Há uma recomendação diferente de Janelas para tipo de sinal, de acordo com o quadro 4.1: Tipo de sinal Hanning Kaiser-Bessel Flat Top Periódico Em geral, é a melhor. Para precisão em freqüência. Para precisão em amplitude. É a melhor. Para precisão em freqüência. Não é recomendada. Transiente Aleatório Quadro 4.1 Uso recomendável de Janelas Retangular Apenas se houver um número inteiro de ciclos dentro de T, o que é bastante difícil de alcançar na prática. É a mais adequada (exceção: Janela exponencial em testes modais com o martelo de impacto). Não é recomendada. 47 A figura 4.13 mostra a FFT do sinal mostrado na figura 4.10, com a janela Hanning. Observa-se a minimização do efeito de leakage. Figura 4.13 FFT do sinal que não apresenta número inteiro de ciclos, com a janela Hanning 48 5 FUNDAMENTAÇÃO TEÓRICA SOBRE AS TÉCNICAS DE ANÁLISE DE SINAIS PARA IDENTIFICAÇÃO DE DEFEITOS 5.1 Introdução A medição e análise de sinais de vibração possibilitam a identificação de falhas nos estágios iniciais, antes da quebra do equipamento. Este procedimento pode ser refinado através da utilização de técnicas de análise de sinais, que indicam o nível de vibração da máquina e o componente que apresenta defeito, a partir de freqüências previamente conhecidas. As técnicas de análise de sinais comumente utilizadas para identificação de defeitos em máquinas são: Nível Global, Fator de Crista, Curtose, Análise Espectral, Técnica do Envelope e Cepstrum. 5.2 Técnicas de análise de sinais 5.2.1 Nível global O Nível Global de vibrações consiste em medir o valor eficaz ou valor RMS (raiz média quadrática) do sinal, proveniente de um transdutor. A análise é realizada através de medições periódicas, cujos valores podem ser mostrados graficamente. Desta forma, acompanha-se a evolução da anomalia em uma máquina ao longo do tempo. Esta técnica indica apenas a condição do equipamento, isto é, ela não informa o tipo de defeito ou qual o componente mecânico apresenta falha. Portanto, esta técnica constitui-se em um “primeiro passo” para o monitoramento de máquinas através da análise de vibrações. O gráfico traçado com os valores das medições é chamado de Curva de Tendência ou Curva de Banheira. Nele observam-se, três níveis: Normal, onde não há ocorrência de problemas; Alerta, onde ocorre um processo de desgaste não perigoso; e Reparo, onde é necessária a manutenção do equipamento, conforme a figura 5.1. 49 Figura 5.1 Curva de tendências Fonte: ANTONIOLLI, 1999 5.2.2 Fator de crista (Crest-Factor) O Fator de Crista (FC) é a relação entre o valor de Pico ( x Pico ) e o valor RMS ( x RMS ) do sinal de vibração, conforme a equação (5.1). FC = x Pico x RMS (5.1) O Fator de Crista deve ser obtido através de medições periódicas, tal que se possa construir um gráfico de FC ao longo do tempo, como mostrado na Figura 5.2a. De uma forma geral, rolamentos devem ser substituídos na transição do 2° para o 3° estágio, ou seja, quando uma redução do Fator de Crista se tornar evidente. Isto acontece porque quando a máquina não possui defeito o nível de vibração é baixo, ou seja, tanto o valor RMS quanto o valor de pico são baixos. O Fator de Crista aumenta a partir do momento em que aparecem pequenos defeitos localizados, que fazem com que o valor de pico aumente, porém o valor RMS mantém-se estável, aumentando a relação entre ambos os valores. Quando os defeitos passam a ficar mais intensos, o nível de vibração se eleva, aumentando também os valores de pico e RMS, o que faz com que a relação diminua. 50 Figura 5.2 Princípio do fator de crista Fonte: ANTONIOLLI, 1999 Figura 5.3 Variação típica do fator de crista durante o processo de degradação de um rolamento Fonte: PONCI e CUNHA 51 Segundo BEZERRA (2004), essa redução do Fator de crista demonstra que esta técnica não é um bom indicador de falhas em um estágio de falha severa. Uma forma encontrada para eliminar este problema foi a criação do Fator K ( Fk ), o qual é o produto do x Pico e x RMS , como mostra a equação (5.2): Fk = x Pico × x RMS (5.2) Na Figura 5.4, estão representados o Fator K e o seu comportamento em relação à condição do rolamento. Figura 5.4 Fator K Fonte: BEZERRA, 2004 5.2.3 Curtose A curtose é definida como o quarto momento estatístico central normalizado pelo desvio padrão, na quarta potência, do sinal de vibração e é representada pela equação (5.4). ∞ m = 4 4 ∫ [x(t ) − µ ] p(x )dx −∞ (5.3) 52 4 1 m k= = 4 σ σ ∞ 4 ∫ [x(t ) − µ ] p(x )dx (5.4) −∞ em que: x(t ) representa o sinal de vibração; µ é a média do sinal; p (x ) é a densidade de probabilidade de x(t ) ; σ é o desvio padrão. Segundo ANTONIOLLI (1999), em um rolamento sem defeito, as vibrações geradas pelo rolamento das esferas sobre as pistas são de natureza aleatória. Além disso, a distribuição de probabilidade de um sinal aleatório tende à gaussiana quando a sua duração é suficientemente grande. A distribuição gaussiana p (x ) é definida por: 1 x(t ) − µ 2 1 exp − p(x ) = 2πσ 2 σ (5.5) A idéia do uso da curtose (k) consiste em avaliar a ocorrência de um defeito através da função densidade de probabilidade. De acordo com MESQUITA (2002), o Fator de Curtose obtido para mancais em bom estado é igual a 3. Porém, com o surgimento e desenvolvimento da falha, os valores deste fator tendem a aumentar. É importante ressaltar que a curtose não consegue localizar um defeito, recomenda-se, portanto, o uso de técnicas complementares, uma vez que corresponde a uma técnica no domínio do tempo. 5.2.4 Análise espectral A análise espectral é uma técnica que consiste em analisar o sinal de vibração a partir do espectro de freqüência, convertendo-se o sinal do domínio do tempo, através de um 53 procedimento matemático chamado de “Transformada Rápida de Fourier” ou, simplesmente, FFT. Os parâmetros fundamentais da análise espectral são: a amplitude e a freqüência, visto que a amplitude do sinal indica a gravidade do problema e a freqüência indica o componente ou o tipo de defeito. Por exemplo, defeitos em engrenagens ocorrem na freqüência de engrenamento, que é o produto entre o número de dentes da engrenagem e a sua rotação. Baseando-se nisto, foi simulado um sinal no LabVIEW, referente a um defeito em um redutor, nas engrenagens 1, 3 e 5. O esquema do redutor é apresentado na figura 5.5. A rotação do eixo do motor é igual a 1094 RPM (18,23 Hz), o que permite determinar as freqüências de engrenamento das engrenagens 1, 3 e 5 como sendo: • Eixo 1: 18,23Hz × 24 dentes = 437,60 Hz • Eixo 2: 3,12 Hz × 20 dentes = 62,51Hz • Eixo 3: 0,65Hz × 35 dentes = 23,03Hz Figura 5.5 Esquema do redutor 54 Figura 5.6 Simulação do sinal de vibração em um redutor No espectro da figura 5.7 aparecem vários picos nas freqüências de engrenamento, bem como várias modulações. O defeito caracterizado por esse tipo de espectro é o erro no módulo e excentricidade na engrenagem. 62,5 Hz Eng. 3 437 Hz Eng. 1 23 Hz Eng. 5 Figura 5.7 Espectro do sinal 5.2.4.1 Análise espectral em aceleração, velocidade e deslocamento Ao analisar um sinal de vibração através do espectro de freqüências é importante saber que em baixas freqüências (até 4 a 5 vezes a rotação do eixo) ocorrem defeitos como desbalanceamento e desalinhamento e devem ser analisados no espectro de freqüência em 55 deslocamento. Por outro lado, defeitos que ocorrem em altas freqüências (acima de 5 vezes a rotação do eixo) devem ser analisados no espectro de freqüência em aceleração. No entanto, quando não se sabe o tipo de defeito existente na máquina, ou seja, se este é de baixa ou alta freqüência, deve-se analisar o sinal no espectro de freqüência em velocidade. A equação (5.6) corresponde ao sinal em deslocamento. Ao derivá-la, obtém-se o sinal em velocidade, equação (5.7). E derivando a equação (5.7), obtém-se o sinal em aceleração, equação (5.8). Observa-se que conforme o sinal é derivado, aumenta-se a sua amplitude. Visto que no sinal em velocidade multiplica-se pela freqüência e no sinal em aceleração multiplicase pela freqüência ao quadrado. Portanto, quanto maior esta freqüência, maior será a amplitude do sinal. Isto explica o fato do sinal em deslocamento acentuar as baixas freqüências e o sinal em aceleração acentuar as altas freqüências. Enquanto que o sinal em velocidade mantém uma média. x(t ) = X sin (ωt ) (5.6) x(t ) = ωX cos(ωt ) • (5.7) •• (5.8) x (t ) = −ω 2 X sin (ωt ) 5.2.5 Técnica do envelope A técnica do envelope é considerada uma das técnicas mais poderosas para detectar e diagnosticar falhas em rolamentos. De acordo com ANTONIOLLI (1999), o envelope corresponde à demodulação das freqüências de defeito do rolamento na alta freqüência de ressonância do mancal, ela reduz os problemas das altas freqüências para as baixas freqüências. Esta técnica mostra claramente onde o defeito está localizado e qual o nível de vibração que este defeito está causando ao mancal. A inspeção ou medição no mancal deve ser periódica para poder identificar as variações do nível de aceleração em relação ao nível inicial de quando o mancal era novo. 56 5.2.5.1 Efeitos de modulação de amplitude em rolamentos Segundo PONCI e CUNHA, Modulação de Amplitude corresponde à variação com o tempo da amplitude de uma vibração (figura 5.8). Neste caso, a componente básica é denominada “Portadora” e o sinal que descreve a variação da amplitude da portadora ao longo do tempo é denominado “Sinal Modulante” ou “Componente Moduladora” e a sua freqüência é denominada “Freqüência de Modulação”. Figura 5.8 Modulação de amplitude Fonte: PONCI e CUNHA PONCI e CUNHA afirmam ainda que o espectro de um sinal modulado (figura 5.9) consta de um pico central na freqüência da portadora (fp), ladeado por picos denominados “Bandas Laterais”, que se situam acima e abaixo da portadora e são espaçados entre si e do pico central por um valor igual à Freqüência de Modulação (fm). Se o sinal modulante for senoidal, existem apenas duas bandas laterais, uma de cada lado da portadora. Se o sinal modulante for complexo, existe uma série de bandas laterais de cada lado da portadora. 57 Figura 5.9 Espectro de um sinal modulado Fonte: PONCI e CUNHA De acordo com PONCI e CUNHA, os defeitos em rolamentos geram impactos repetitivos que excitam freqüências naturais elevadas (maior que 500 Hz), as quais se propagam pela caixa dos mancais e estrutura da máquina na forma de ondas de tensão. Assim, essas vibrações geradas pelos defeitos em rolamentos são moduladas em amplitude pela seqüência de impactos repetitivos e pelo efeito do amortecimento. As portadoras são as freqüências naturais dos componentes do rolamento (pistas e elementos rolantes) e da caixa. As moduladoras são as freqüências básicas de defeito. 5.2.5.2 Demodulação de sinais para a detecção de defeitos em rolamentos (Envelope) PONCI e CUNHA afirmam que as vibrações de alta freqüência geradas pelos defeitos de rolamentos possuem intensidade bastante reduzida em comparação com as amplitudes das componentes de baixa freqüência, relacionadas com defeitos de alta energia como: desbalanceamento, desalinhamento, e folgas. A grande vantagem do processo de demodulação é eliminar essas componentes de alta energia, permitindo detectar com maior precisão e antecedência só os defeitos de rolamentos. Além disso, a análise dos espectros de envelope permite determinar as taxas de repetição dos impactos que geram as ondas de tensão, identificando a sua origem, ou seja, os componentes defeituosos do rolamento. 58 O procedimento de aplicação da Técnica do Envelope é mostrado na figura 5.10, onde: I) Sinal de vibração no domínio do tempo medido na superfície do mancal; II) Espectro de freqüência correspondente ao sinal (I), obtido através da Transformada de Fourier; III) Sinal (I) filtrado utilizando um filtro passa-banda. Contém somente as freqüências de ressonância que são moduladas pela freqüência dos impactos; IV) Espectro de freqüência correspondente ao sinal (III); V) Sinal (III) demodulado após a aplicação da Técnica do Envelope; VI) Espectro de freqüência correspondente ao sinal (V). O envelope do sinal pode ser obtido através de um aparelho analógico, denominado “detector de envelope”, ou de forma digital, através da Transformada de Hilbert. Figura 5.10 Esquema da técnica do envelope Fonte: NUNES 1989 59 Para entender melhor esta técnica, simulou-se no LabVIEW um sinal característico de um defeito em rolamento. Este defeito corresponde ao impacto de uma falha interna excitando somente uma ressonância. O sinal é mostrado na figura 5.11 e possui duas freqüências, uma portadora com valor igual a 960 Hz que corresponderia à freqüência natural do mancal e uma moduladora com valor igual a 160 Hz que seria a freqüência de defeito no rolamento. Figura 5.11 Simulação de defeito em rolamento A figura 5.12 corresponde ao espectro do sinal, obtido através da Transformada de Fourier. Como pode ser observado, o sinal apresenta modulação, com um pico central na freqüência de 960 Hz e bandas laterais espaçadas em 160 Hz. 960 Hz Figura 5.12 Espectro do sinal O sinal foi filtrado na freqüência de 960 Hz, utilizando um filtro passa-banda com freqüências de corte de 800 Hz e 1120 Hz (figura 5.13). 60 Figura 5.13 Sinal filtrado A figura 5.14 é referente ao espectro do sinal filtrado. 960 Hz Figura 5.14 Espectro do sinal filtrado A figura 5.15 é relativa ao envelope do sinal, obtido através da aplicação da Transformada de Hilbert. Observa-se que a demodulação do sinal é feita extraindo-se somente a envoltória do sinal. Outra característica interessante desta técnica é a conversão do sinal para valores positivos, por isso também recebe o nome de “Módulo da Transformada de Hilbert”. 61 Figura 5.15 Envelope do sinal O espectro de freqüência do envelope, mostrado na figura 5.16, apresenta um pico na freqüência de 160 Hz (freqüência de defeito no rolamento) e em 320 Hz (múltiplo de 160 Hz). 160 Hz 320 Hz, múltiplo de 160 Hz Figura 5.16 Espectro do envelope 5.2.6 Cepstrum O Cepstrum é uma técnica que tem sido utilizada principalmente no diagnóstico de defeitos em engrenagens e em rolamentos. Segundo ALMEIDA e GÓZ (2000), o Cepstrum identifica e mostra famílias de picos, com espaçamento igual, presentes no espectro de freqüência de um sinal. As principais causas da geração destas famílias são as diversas formas de modulação, principalmente a de freqüência. Portanto, a análise cepstral é uma ferramenta poderosa e de uso simples na descoberta e identificação de processos de modulação. 62 O Cepstrum ou mais precisamente, o Cepstrum de Potência foi originalmente definido por B. P. Borget, como “o espectro de potência do logaritmo do espectro de potência”, ou matematicamente: C p (τ ) = ℑ{log Gxx ( f )} (5.9) 2 onde: Gxx ( f ) = ℑ{xk (t )} 2 (5.10) é o espectro de potência do sinal temporal xk (t) e τ { } representa a Transformada de Fourier da função entre chaves. Mais tarde adotou-se uma outra definição para o Cepstrum, “A transformada inversa do logaritmo do espectro de potência”, ou seja: C p (τ ) = ℑ−1{log Gxx ( xk (t ))} (5.11) ALMEIDA e GÓZ (2000) afirmam que é possível demonstrar que as duas definições geram o mesmo resultado, isto é, o mesmo Cepstrum. A demonstração baseia-se no fato de que para uma função real e par, tal como o Espectro de Potência, as transformadas direta e inversa de Fourier dão o mesmo resultado, exceto por um possível fator de escala. Ainda segundo ALMEIDA e GÓZ (2000), a palavra “cepstrum” é derivada da inversão das primeiras letras da palavra “spectrum”. Pelo mesmo processo a palavra “freqüência” gerou o nome da variável independente do cepstrum, a “quefrência”, cuja unidade é τ (tau) com dimensão de segundo (tempo). A unidade da quefrência apesar de ser segundo, não é igual ao tempo normal (do relógio), pois não foi obtida pela volta ao domínio da freqüência pela transformada inversa de Fourier, mas sim de uma nova transformada aplicada no domínio da freqüência. A idéia do tempo no cepstrum representa a idéia de repetição, proximidade de componentes de freqüência no espectro de freqüência normal. Por esta concepção de tempo do cepstrum, origina-se a habilidade da análise cepstral de mostrar as famílias de picos que mantêm entre si o mesmo espaçamento de freqüência. Assim, altos valores de quefrência representam menores espaçamentos de freqüência, e baixa quefrência representa maiores espaçamentos. E a amplitude da componente cepstral é 63 proporcional ao tamanho da família identificada, ou seja, o número de componentes igualmente espaçadas. A figura 5.17 apresenta um sinal modulado por duas freqüências de 10 Hz e 100 Hz. Figura 5.17 Sinal modulado A figura 5.18 corresponde ao espectro de potência do sinal modulado, mostrando dois picos nas freqüências de 90 Hz e 110 Hz, as quais equivalem a maior freqüência (100 Hz) menos a menor (10 Hz) e a maior (100 Hz) mais a menor (10 Hz). 90 Hz 110 Hz Figura 5.18 Espectro de potência do sinal Na Figura 5.19, tem-se o cepstrum do sinal. Os picos obtidos têm um espaçamento de 0,1 s cujo valor inverso 1/0,1 = 10 Hz, que é a freqüência de modulação do sinal. 64 0,1 s Figura 5.19 Cepstrum do sinal 0,2 s 65 6 DESCRIÇÃO DO APLICATIVO DESENVOLVIDO PARA AQUISIÇÃO E ANÁLISE DE SINAIS DE VIBRAÇÃO 6.1 Introdução O aplicativo desenvolvido para aquisição e análise de sinais de vibração utiliza como base a plataforma LabVIEW. O LabVIEW é um software aplicativo baseado na linguagem G ou linguagem de programação gráfica, que emprega ícones ao invés de textos para criar aplicações. Esta linguagem é diferente das linguagens de programação baseadas em textos (linhas de comando), onde as instruções determinam a execução do programa. A figura 6.1 apresenta um exemplo de programação através de linhas de comando e a figura 6.2 refere-se à programação através de ícones. A programação em LabVIEW está baseada no fluxo de dados, onde os dados determinam a execução, trazendo algumas vantagens para aplicações científicas e de Engenharia, principalmente em aplicações de aquisição e manipulação de dados. Os aplicativos são desenvolvidos pelo usuário, utilizando um conjunto de ferramentas e objetos que possuem funções para aquisição, análise e apresentação dos dados, GPIB e controle de instrumentos seriais. Os códigos são adicionados no diagrama de bloco usando representações gráficas de funções para controlar os objetos adicionados no painel frontal. Depois de criado, o diagrama de bloco é compilado para linguagem de máquina. Figura 6.1 Exemplo de programação através de linhas de comando 66 Figura 6.2 Exemplo de programação através de ícones 6.1.1 Instrumentos virtuais Os programas do LabVIEW são chamados de Instrumentos Virtuais, ou VIs, porque sua aparência e operação imitam instrumentos físicos, tais como osciloscópios e multímetros. Cada VI tem funções que manipulam a entrada pela interface do usuário ou de outras fontes, através da aquisição de sinais, e indicam essa informação ou a movem para outros arquivos ou outros computadores. Um VI contém os três seguintes componentes: • Painel Frontal – Serve como interface do usuário (figura 6.3). • Diagrama de Bloco – Contém o código fonte gráfico que define o funcionamento do VI (figura 6.4). • Ícone e Conector – Identificam o VI, de modo que este possa ser usado em um outro VI (figuras 6.5 e 6.6). Um VI dentro de um outro VI é chamado de subVI. Um subVI corresponde a uma sub-rotina em linguagens de programação baseadas em texto. 67 Figura 6.3 Painel frontal Figura 6.4 Diagrama de bloco Figura 6.5 Ícone Figura 6.6 Conector 68 6.2 Descrição do aplicativo desenvolvido para aquisição e análise de vibração O aplicativo desenvolvido foi intitulado como “Programa de Detecção de Defeitos em Máquinas Rotativas”. Ele possui diversos módulos, responsáveis pela aquisição e visualização de sinais, cálculo de freqüências características de defeitos em rolamentos e aplicação das técnicas comumente utilizadas para análise de sinais de vibração. Cada módulo desses corresponde a um subVI. A tela inicial é composta por sete botões, como pode ser visto na figura 6.7: Figura 6.7 Tela inicial do aplicativo 69 A figura 6.8 indica a função de cada botão. ⇒ Abre o módulo de aquisição de sinais. ⇒ Abre o módulo de visualização de sinais. ⇒ Abre o módulo de cálculo de freqüências características de defeitos em rolamentos. ⇒ Abre o módulo da técnica do envelope. ⇒ Abre o módulo do cepstrum. ⇒ Abre o módulo do nível global, fator de crista e curtose. ⇒ Pára a execução do programa. Está presente em todos os módulos. Figura 6.8 Botões da tela inicial 6.2.1 Módulo de aquisição de sinais Este módulo possui a configuração responsável pela aquisição do sinal. O painel frontal, mostrado na figura 6.9, apresenta dois gráficos, um no domínio do tempo e outro no domínio da freqüência; controles para configuração da amostragem do sinal, isto é, o número de amostras e taxa de aquisição; diretório onde deve ser salvo o arquivo correspondente ao sinal de vibração aquisitado; configuração de filtro; canal de aquisição; e o botão “PARAR”. 70 Figura 6.9 Tela de aquisição de sinais A figura 6.10 apresenta os controles para configuração da amostragem, cuja informação é fundamental para a correta aquisição do sinal. Figura 6.10 Controles para configuração de amostragem do sinal A figura 6.11 apresenta o controle utilizado para escolha do diretório onde será salvo o arquivo com os dados do sinal adquirido. O botão “Desabilitar” deve ser clicado para ativar o arquivamento do sinal, ou seja, quando o programa é executado, a aquisição inicia automaticamente, porém o sinal só será salvo a partir do momento em que for habilitado. 71 Botão para escolha do diretório Figura 6.11 Informação do diretório onde o arquivo deve ser salvo O arquivo gerado possui a extensão (.lvm) e pode ser lido no programa “Bloco de notas”. Além de todos os pontos que compõem o sinal, o arquivo apresenta um cabeçalho com informações básicas, tais como a data e a hora da aquisição, o canal, o número de amostras, e a unidade do eixo Y, que neste caso é Volts. A figura 6.12 mostra um exemplo de um arquivo. A figura 6.13 refere-se à configuração do Filtro, do tipo “passa-banda”, isto é, as freqüências de corte inferior e superior. Na figura 6.14 é mostrado o controle onde é realizada a seleção do canal que fará a aquisição do sinal. A figura 6.15 apresenta o diagrama de bloco que contém o código fonte para aquisição do sinal. Figura 6.12 Exemplo de arquivo (.lvm) que contém os dados do sinal de vibração 72 Figura 6.13 Configuração de filtro Figura 6.14 Configuração do canal de aquisição Figura 6.15 Diagrama de bloco do módulo de aquisição de sinais 6.2.2 Módulo de visualização de sinais Este módulo é responsável pela visualização dos sinais, após a aquisição e arquivamento. Ele apresenta quatro gráficos: um no domínio do tempo, um espectro de freqüência em aceleração, em velocidade e outro em deslocamento. Possui também controles de configuração referentes ao sensor utilizado na aquisição; controles para abrir o arquivo salvo; indicadores que informam os valores de Pico, RMS, Fator de Crista e Curtose do sinal. As figuras 6.16 e 6.17 mostram o painel frontal do módulo de visualização de sinais. 73 Figura 6.16 Tela de visualização de sinais em aceleração no domínio do tempo e da freqüência Figura 6.17 Tela de visualização de sinais em velocidade e deslocamento no domínio da freqüência 74 A figura 6.18 exibe os controles utilizados para configuração do sensor responsável pela aquisição do sinal. Antes de realizar a leitura do sinal, deve-se informar dados como a sensibilidade do sensor, em mV/EU, e a unidade de engenharia (EU). Estes controles estão presentes também nos módulos da técnica do envelope e do cepstrum, que serão explicados posteriormente. Figura 6.18 Configuração do sensor A figura 6.19 indica o botão que abre o arquivo gerado durante a aquisição do sinal. O sinal só será mostrado a partir do momento que o botão “Desabilitar” for clicado, ou seja, ele estará habilitando a leitura do sinal. Assim como, os controles de configuração do sensor, este controle está presente nos módulos da técnica do envelope e do cepstrum. Botão que abre o arquivo Figura 6.19 Controle que abre o arquivo Na figura 6.20 são apresentados os indicadores que mostram os valores de Pico, RMS, Fator de Crista e Curtose do sinal. 75 Figura 6.20 Indicadores de valores de Pico, RMS, Fator de Crista e Curtose do sinal A figura 6.21 apresenta o diagrama de bloco que contém o código fonte para visualização do sinal em aceleração no domínio do tempo e no domínio da freqüência. E a figura 6.22 apresenta o diagrama de bloco que contém o código fonte para visualização do sinal em velocidade e em deslocamento no domínio da freqüência. A função janela é definida pela função FFT, a qual é mostrada em destaque nas figuras 6.21 e 6.22. Figura 6.21 Diagrama de bloco do módulo de visualização de sinais 76 Figura 6.22 Diagrama de bloco do módulo de visualização de sinais 6.2.3 Módulo de cálculo de freqüências características de defeitos em rolamentos Este módulo realiza o cálculo das freqüências características de defeitos em rolamento, referentes aos diferentes componentes do mancal. A figura 6.23 mostra o painel frontal e a figura 6.24 ilustra o diagrama de bloco. Dados como número de elementos rolantes, freqüência de rotação da máquina (Hz), diâmetro dos elementos rolantes (mm), diâmetro primitivo (mm) e ângulo de contato devem ser inseridos. O programa, então calcula e informa os valores das freqüências relativas a defeitos na pista externa, pista interna, gaiola e elemento rolante. 77 Figura 6.23 Tela de cálculo de freqüências características de defeitos em rolamentos Figura 6.24 Diagrama de bloco do módulo de cálculo de freqüências características de defeitos em rolamentos 6.2.4 Módulo da técnica do envelope Neste módulo foi aplicada a técnica do envelope. Ele é composto por dois gráficos: um que mostra o envelope do sinal e outro que mostra o espectro de freqüência do envelope. Além disso, possui controles para configuração do sensor, controles para abrir o arquivo salvo, e controles para configuração do filtro, conforme a figura 6.25. 78 Figura 6.25 Tela da técnica do envelope A figura 6.26 apresenta os controles de configuração do filtro, onde devem ser indicadas as freqüências de corte inferior e superior. Figura 6.26 Configuração de filtro A figura 6.27 mostra o código fonte da rotina desenvolvida para a técnica do envelope. 79 Figura 6.27 Diagrama de bloco do módulo da técnica do envelope 6.2.5 Módulo do cepstrum Este módulo é referente à técnica do cepstrum. Ele é composto por dois gráficos: um que mostra o espectro de potência do sinal e outro que mostra o cepstrum. Além disso, possui controles para configuração do sensor, e controles para abrir o arquivo salvo, como mostra a figura 6.28. Figura 6.28 Tela do Cepstrum 80 A Figura 6.29 apresenta o diagrama de bloco do módulo do Cepstrum. O espectro de potência é definido pela função FFT, a qual é mostrada em destaque na figura 6.29. Figura 6.29 Diagrama de bloco do módulo do cepstrum 6.2.6 Módulo do nível global, fator de crista e curtose Este módulo é referente às três técnicas, separadas em guias: Nível Global, Fator de Crista e Curtose. A sua composição é bem simples: possui três gráficos e três tabelas, um para cada técnica. Ao ler o sinal no módulo de visualização, o usuário possui a informação destes valores. Através de medições periódicas, ele pode inserir estes valores nas tabelas e gerar os gráficos. A figura 6.30 mostra a tela do nível global, a figura 6.31 apresenta a tela do fator de crista e a figura 6.32 é referente à tela da curtose. As figuras 6.33, 6.34 e 6.35 mostram o diagrama de bloco deste módulo. Cada técnica está contida em um ‘caso’ do ‘Case Structure’, o qual é submetido a uma condição amarrada ao ‘Case Selector’. O ‘Selector Label’ mostra o primeiro caso, correspondente à técnica do nível global, o segundo corresponde ao fator de crista e o terceiro é a curtose. 81 Figura 6.30 Tela do nível global 6.31 Tela do fator de crista 82 6.32 Tela da curtose Selector Label Case Selector Figura 6.33 Diagrama de bloco que executa o nível global 83 Figura 6.34 Diagrama de bloco que executa o fator de crista Figura 6.35 Diagrama de bloco que executa a curtose 84 7 DESCRIÇÃO DA BANCADA EXPERIMENTAL 7.1 Introdução Este capítulo é referente à descrição da bancada experimental, concebida para simulação e estudo de defeitos em máquinas rotativas e validação do aplicativo desenvolvido. Além disso, serão mostrados os sensores utilizados e de que forma os defeitos foram induzidos na bancada de testes. Antes de iniciar a descrição da bancada, é necessário entender como é realizada a aquisição do sinal de vibração. A figura 7.1 apresenta um esquema do sistema de aquisição de sinais, o qual é composto por um sensor de vibração, que neste caso é o acelerômetro; um condicionador de sinais; uma placa de aquisição de sinais; e um computador com o software de aquisição de sinais. Figura 7.1 Esquema do sistema de aquisição de sinais de vibração • Condicionador de sinais: Alimenta o sensor, amplifica e filtra o sinal antes de enviar para a placa de aquisição e ajusta a impedância elétrica. • Placa de Aquisição: Converte o sinal analógico em digital. • Computador com Software de aquisição de sinais: Interface onde são visualizados e analisados os sinais. 85 7.2 Bancada experimental A bancada experimental utilizada neste trabalho é apresentada na figura. 7.2. Esta bancada é comercializada pela Spectra Quest Inc. e foi especialmente projetada para estudar o comportamento dinâmico de rotor apoiado em mancais de deslizamento e/ou de rolamento, permitindo a simulação de problemas de desbalanceamento, desalinhamento e outros específicos ligados a um determinado componente, como, por exemplo, defeitos em rolamentos. Os seus principais elementos são: • Controlador AC de velocidade variável com mostrador (display) • Motor AC de ½ HP • Acoplamento • Eixo • Tacômetro óptico • Disco rotativo com furos rosqueados para introdução do desbalanceamento • Dispositivo de carga do mancal • Anel de fixação para retenção do rotor ao eixo • Tanque integral da bomba de óleo Figura 7.2 Bancada experimental 86 7.3 Cadeia de medição A cadeia de medição é composta pelos seguintes itens: • 3 Sensores de vibração (acelerômetros) • Condicionador de sinais • Placa de aquisição de sinais • Computador com software de aquisição de sinais. 7.3.1 Sensor de vibração O sensor de vibração utilizado na bancada é o acelerômetro, do tipo piezoelétrico. Segundo SOUZA (2007), o nome piezoelétrico vem do grego peizin, que significa esmagar, achatar. É que os elementos piezoelétricos têm a propriedade de produzir uma tensão elétrica em determinado eixo quando comprimido em outro. E vice-versa: se for aplicado um sinal elétrico no eixo elétrico, ele vai produzir um alongamento mecânico no eixo mecânico. O cristal comum de quartzo (dióxido de silício) apresenta essa propriedade e por isso foi muito usado como cristais geradores de freqüências-padrão em radiofonia, para microfones, fones e acelerômetros. Hoje os cristais de quartzo são criados em fornos com a temperatura e pressão controlados, onde se obtém uma pureza sem a qual seria impossível obter a precisão dos atuais transdutores. Algumas cerâmicas, como o titanato de bário, também apresentam propriedades piezoelétricas, desde que submetidas a um tratamento termoelétrico durante o processo de fusão. Estes últimos cristais são mais utilizados como geradores de ultra-som devido à alta freqüência de ressonância mecânica aliada à alta eficiência mecânica. A figura 7.3 apresenta um corte num acelerômetro, no qual é mostrado o cristal, de um lado apoiado no corpo do transdutor e de outro fixado a massa, também conhecida como massa sísmica. Qualquer aceleração sofrida pelo transdutor obrigará a massa a comprimir ou tracionar o cristal, que produzirá um sinal elétrico proporcional à força gerada na massa. Apenas a componente da aceleração paralela ao eixo do transdutor é que produzirá o efeito piezoelétrico; as outras componentes no plano perpendicular ao corpo do transdutor serão 87 ignoradas. Quando se deseja medir acelerações nos três eixos coordenados, deve-se utilizar um transdutor triaxial que contém três cristais independentes perpendicularmente entre si e fornece três tensões independentes de saída. Figura 7.3 Esquema do acelerômetro piezoelétrico Fonte: SOUZA, 2007 7.3.1.1 Características do acelerômetro utilizado na bancada O acelerômetro utilizado na bancada é mostrado na Figura 7.4 e apresenta as seguintes características: • Fabricante: Wilcoxon • Modelo: 784A • Tipo: Piezoelétrico ICP-PZT CERÂMICO • Sensibilidade axial: - 100mV/g +/- 20% • Sensibilidade transversal: - 3% da axial • Sensibilidade magnética: 50 µg/Gauss • Faixa de temperatura de operação: -50° C a 120° C • Alimentação: 12 VDC Figura 7.4 Acelerômetro piezoelétrico, modelo 784A. Fonte: WILCOXON A sensibilidade do sensor é um parâmetro importantíssimo na leitura do sinal, visto que um analisador de sinais por ser um equipamento eletrônico, lê o sinal adquirido na forma 88 de tensão (Volts). No entanto, um dos parâmetros de medida de vibração é a aceleração. Portanto, é necessário converter a unidade de volts para m/s². A relação é dada pela sensibilidade do sensor. Neste caso, a sensibilidade do sensor é igual a 100 mV/g, ou seja, para cada 0,1 V tem-se 9,81 m/s². 7.3.1.2 Posicionamento dos acelerômetros na bancada Foram utilizados três acelerômetros para aquisitar o sinal de vibração, nas posições axial, radial vertical e radial horizontal, conectados no mancal de rolamento, como mostrado na figura 7.5. Sensores de vibração Figura 7.5 Bancada experimental com sensores de vibração 7.3.2 Condicionador de sinais O condicionador de sinais utilizado na bancada possui características, como: 16 canais e um filtro analógico passa-baixa, com freqüência de corte igual a 2 kHz. A figura 7.6 mostra o condicionador de sinais. 89 Figura 7.6 Condicionador de sinais 7.3.3 Placa de aquisição de sinais A placa de aquisição de sinais é mostrada na figura 7.7 e apresenta as seguintes características: • Fabricante: National Instruments; • Modelo: NI PCI-6024E; • 16 Entradas analógicas; • 200 kS/s, 12-Bit; • 2 saídas analógicas 12-bit; 8 linhas digitais I/O; dois contadores 24-bit. Figura 7.7 Placa de aquisição de sinais Fonte: NATIONAL INSTRUMENTS 90 7.4 Indução de defeitos na bancada 7.4.1 Desbalanceamento de massa O desbalanceamento foi provocado através da inserção de massas desbalanceadoras, que correspondem a dois parafusos de 5,55 g e 5,56 g, aproximadamente, fixados no disco rotativo do eixo, com uma excentricidade de 4 cm e 6 cm. A figura 7.8 mostra o parafuso e a figura 7.9 apresenta o disco com os parafusos fixados. Figura 7.8 Massa desbalanceadora Massas desbalanceadoras Figura 7.9 Disco rotativo com as massas desbalanceadoras 91 7.4.2 Desalinhamento de eixo Para provocar o desalinhamento utilizou-se um acoplamento desalinhado, o qual é mostrado nas figuras 7.10 e 7.11. Acoplamento desalinhado Figura 7.10 Acoplamento desalinhado na bancada Figura 7.11 Acoplamento desalinhado A figura 7.12 apresenta o acoplamento sem desalinhamento utilizado para a realização das medições e posterior comparação de sinais. 92 Figura 7.12 Acoplamento sem desalinhamento 7.4.3 Defeitos nos rolamentos Foram utilizados dois rolamentos, um com defeito na pista interna e outro com defeito na pista externa, mostrados na figura 7.13. Figura 7.13 Rolamento utilizado para indução de defeitos Os rolamentos apresentam as seguintes características: • Fabricante: NYLA-K/MB Manufacturing • Modelo: M-BFK-1 • Número de elementos rolantes (n): 8 • Diâmetro dos elementos rolantes (D): 15,85 mm • Diâmetro primitivo (dm): 31,5 mm • Ângulo de contato (β): 0° 93 8 MEDIÇÕES E ANÁLISE DOS SINAIS DE VIBRAÇÃO 8.1 Introdução Neste capítulo serão apresentados os seguintes itens: a calibração do sistema de aquisição de sinais, as medições realizadas na bancada de teste e as análises dos sinais de vibração, através do Nível Global, Fator de Crista, Curtose, Análise Espectral e Técnica do Envelope. 8.2 Calibração do sistema de aquisição de sinais A calibração foi realizada, no Laboratório de Calibração da Eletronorte (figura 8.1), para cada acelerômetro de acordo com cada canal do sistema de aquisição de sinais. Utilizou-se para este procedimento um gerador de sinal e um “shaker”. A calibração consistiu em gerar um sinal, com uma determinada freqüência e amplitude, o qual era enviado ao “shaker”, fazendo com que este vibrasse a uma freqüência e amplitude conhecidas. O sinal de vibração aquisitado era lido no computador, e então estabelecidas as constantes de calibração, as quais são mostradas na quadro 8.1. Canal Constante 0 326,20 1 325,20 2 328,85 Quadro 8.1 Constantes de calibração 94 Figura 8.1 Calibração do sistema 8.3 Análise dos sinais de vibração Todos os defeitos foram medidos nas posições axial, radial horizontal e radial vertical, pois como já foi dito no capítulo 3, o desbalanceamento é identificado através da análise do sinal na posição radial, o qual deve ser comparado com o sinal na posição axial. Enquanto que para o desalinhamento, estas diferentes posições do sensor são importantes para se identificar o tipo de desalinhamento. A freqüência de rotação da máquina é igual a 20 Hz. A função Janela utilizada para análise dos sinais foi a Hanning. Como o condicionador de sinais possui um filtro analógico passa-baixa, com uma freqüência de corte igual a 2 kHz, foi utilizada uma taxa de aquisição de 4 kHz e um número de amostras igual a 8000. 95 8.3.1 Análise do defeito de desbalanceamento O desbalanceamento foi identificado através da técnica análise espectral em deslocamento. Além disso, são indicados os valores de nível global com e sem desbalanceamento. Os gráficos no domínio do tempo e no domínio da freqüência, em aceleração e em velocidade, podem ser vistos no APÊNDICE A. A figura 8.2 corresponde ao sinal da máquina sem desbalanceamento e a figura 8.3 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição axial. Figura 8.2 Espectro em deslocamento / sinal sem desbalanceamento / posição axial Figura 8.3 Espectro em deslocamento / sinal com desbalanceamento / posição axial A figura 8.4 corresponde ao sinal da máquina sem desbalanceamento e a figura 8.5 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial horizontal. 96 Figura 8.4 Espectro em deslocamento / sinal sem desbalanceamento / posição radial horizontal 20 Hz, freqüência de rotação da máquina, indicando o desbalanceamento Figura 8.5 Espectro em deslocamento / sinal com desbalanceamento / posição radial horizontal A figura 8.6 corresponde ao sinal da máquina sem desbalanceamento e a figura 8.7 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial vertical. Figura 8.6 Espectro em deslocamento / sinal sem desbalanceamento / posição radial vertical 97 20 Hz, freqüência de rotação da máquina, indicando o desbalanceamento Figura 8.7 Espectro em deslocamento / sinal com desbalanceamento / posição radial vertical Nas posições radiais (figuras 8.5 e 8.7), observa-se um aumento da amplitude do sinal, indicando o desbalanceamento de massa. Na verdade, nos sinais das figuras 8.4 e 8.6, há desbalanceamento, mas é um desbalancamento residual, que não compromete o bom funcionamento do equipamento. Os valores de nível global da máquina com e sem desbalanceamento são mostrados no quadro 8.2. MÁQUINA SEM DESBALANCEAMENTO Nível global (m/s²) Axial Radial Horizontal Radial Vertical 0,85 0,38 0,44 MÁQUINA COM DESBALANCEAMENTO Nível global (m/s²) Axial Radial Horizontal Radial Vertical 0,67 0,56 0,67 Quadro 8.2 Valores de nível global 8.3.2 Análise do defeito de desalinhamento O defeito de desalinhamento foi analisado também através da análise espectral em deslocamento. Além disso, são indicados os valores de nível global com e sem desalinhamento. Os gráficos no domínio do tempo e no domínio da freqüência, em aceleração e em velocidade, podem ser vistos no APÊNDICE A. A figura 8.8 corresponde ao sinal da máquina sem desalinhamento e a figura 8.9 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição axial. 98 Figura 8.8 Espectro em deslocamento / sinal sem desalinhamento / posição axial Figura 8.9 Espectro em deslocamento / sinal com desalinhamento / posição axial A figura 8.10 corresponde ao sinal da máquina sem desalinhamento e a figura 8.11 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial horizontal. Figura 8.10 Espectro em deslocamento / sinal sem desalinhamento / posição radial horizontal 99 20 Hz, 40 Hz e 60 Hz. Harmônicos da freqüência de rotação da máquina, indicando o desalinhamento. Figura 8.11 Espectro em deslocamento / sinal com desalinhamento / posição radial horizontal A figura 8.12 corresponde ao sinal da máquina sem desalinhamento e a figura 8.13 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial vertical. Figura 8.12 Espectro em deslocamento / sinal sem desalinhamento / posição radial vertical Figura 8.13 Espectro em deslocamento / sinal com desalinhamento / posição vertical O desalinhamento pôde ser identificado através dos três harmônicos da freqüência de rotação da máquina (figura 8.11), e foram significativos na posição radial horizontal, indicando um desalinhamento paralelo, no plano horizontal. 100 Os valores de nível global da máquina com e sem desalinhamento são mostrados no quadro 8.3 MÁQUINA SEM DESALINHAMENTO Nível global (m/s²) Axial Radial Horizontal Radial Vertical 0,85 0,38 0,44 MÁQUINA COM DESALINHAMENTO Nível global (m/s²) Axial Radial Horizontal Radial Vertical 2,06 5,83 2 Quadro 8.3 Valores de nível global 8.3.3 Análise dos defeitos nos rolamentos As freqüências características dos defeitos em rolamento foram calculadas utilizando o aplicativo desenvolvido, como mostra a figura 8.14. Os valores são apresentados no quadro 8.4. Figura 8.14 Tela do aplicativo que calcula as freqüências características dos defeitos em rolamento 101 Componente Freqüência (Hz) Pista Interna 120,25 Pista Externa 39,75 Gaiola 4,97 Esfera 29,68 Quadro 8.4 Freqüências calculadas Os defeitos nos rolamentos foram analisados por análise espectral e técnica do envelope. São mostrados também os valores de nível global, fator de crista e curtose. 8.3.3.1 Análise do defeito na pista interna Primeiramente analisou-se pelo espectro de freqüência em aceleração (figura 8.16). O sinal foi filtrado utilizando um filtro passa-banda, com freqüências de corte iguais a 350 e 700 Hz, considerando que esta fosse a faixa de freqüências de ressonância, pois apresenta maior amplitude. A figura 8.17 mostra o envelope do sinal e a figura 8.18 mostra o espectro do envelope, onde pode ser identificada a freqüência de 120,5 Hz, relativa ao defeito provocado na pista interna do rolamento, a qual foi previamente calculada como 120,25 Hz. Figura 8.15 Domínio do tempo / posição axial / defeito na pista interna 102 Figura 8.16 Domínio da freqüência / posição axial / defeito na pista interna Figura 8.17 Envelope / Domínio do tempo / posição axial / defeito na pista interna 120,5 Hz Figura 8.18 Envelope / Domínio da freqüência / posição axial / defeito na pista interna 103 ROLAMENTO SEM DEFEITO Axial Nível global (m/s²) 0,85 Fator de crista 5,54 ROLAMENTO COM DEFEITO NA PISTA INTERNA Axial Nível global (m/s²) 5,94 Fator de crista 3,54 Curtose 3,24 Quadro 8.5 Valores de nível global, fator de crista e curtose 8.3.3.2 Análise do defeito na pista externa Inicialmente analisou-se pelo espectro de freqüência em aceleração (figura 8.20). O sinal foi filtrado utilizando um filtro passa-banda, com freqüências de corte iguais a 500 e 900 Hz, considerando que esta fosse a faixa de freqüências de ressonância do mancal, apesar de não apresentar maior amplitude. A região que apresenta maior amplitude está entre 200 e 400 Hz, no entanto, as freqüências de ressonância tendem a ser mais altas que esses valores. A figura 8.21 mostra o envelope do sinal e a figura 8.22 mostra o espectro do envelope, onde pode ser identificada a freqüência de 39,5 Hz, relativa ao defeito provocado na pista externa do rolamento, a qual foi previamente calculada como 39,75 Hz. Figura 8.19 Domínio do tempo / posição axial / defeito na pista externa 104 Figura 8.20 Domínio da freqüência / posição axial / defeito na pista externa Figura 8.21 Envelope / Domínio do tempo / posição axial / defeito na pista externa 39,5 Hz Figura 8.22 Envelope / Domínio da freqüência / posição axial / defeito na pista externa 105 ROLAMENTO SEM DEFEITO Axial Nível global (m/s²) 0,85 Fator de crista 5,54 ROLAMENTO COM DEFEITO NA PISTA EXTERNA Axial Nível global (m/s²) 1,57 Fator de crista 4,02 Curtose 3,38 Quadro 8.6 Valores de nível global, fator de crista e curtose 106 9 CONCLUSÕES E RECOMENDAÇÕES 9.1 Conclusões Através do estudo das falhas mecânicas mais comuns em máquinas rotativas e seus efeitos de vibração e das técnicas de análise de sinais utilizadas para o tratamento de dados e diagnóstico de defeitos, foi possível desenvolver o sistema aplicativo proposto para aquisição e análise de sinais de vibração, tendo por base a plataforma LabVIEW. Em virtude das medições realizadas na bancada experimental, pôde-se validar o sistema aplicativo, verificando os diversos módulos desenvolvidos, tais como o de aquisição de sinais, de visualização de sinais e das técnicas de análise de sinais. De acordo com os defeitos induzidos na bancada, pôde-se avaliar as técnicas de análise de sinais utilizadas no aplicativo. Com o desbalanceamento de massa e desalinhamento de eixo avaliou-se o Nível global e Análise espectral; com os defeitos nos rolamentos verificou-se o Nível global, Fator de crista, Curtose, Análise espectral e Técnica do envelope. A única técnica que não pôde ser validada foi a do Cepstrum, visto que esta é mais indicada para análise de defeitos em engrenagens. Este aplicativo mostrou-se prático e viável, uma vez que instalado em um computador portátil, que possua uma placa de aquisição de sinais, pode ser transportado com facilidade para qualquer lugar. Além disso, requer pouco espaço no disco rígido do computador. E apresenta um custo inferior em relação a equipamentos convencionais de medição. 9.2 Sugestões para trabalhos futuros Considerando que a única técnica que não foi avaliada foi a do Cepstrum, recomendase: • Utilizar uma bancada que contenha um Redutor, para indução de defeitos nas engrenagens e aplicação da técnica do Cepstrum; • Utilizar banco de dados para armazenamento dos sinais de vibração; • Aplicar a ferramenta virtual em um equipamento de maior porte. 107 REFERÊNCIAS ALMEIDA, Márcio Tadeu; GÓZ, Ricardo Damião Sales. Análise de Vibrações I – Medidas e Diagnósticos. Apostila do Curso de Análise de Vibração da FUPAI, Itajubá – MG. 2000. ALMEIDA, Márcio Tadeu; GÓZ, Ricardo Damião Sales. Análise de Vibrações II – Rolamentos e Engrenagens. Apostila do Curso de Análise de Vibração da FUPAI, Itajubá – MG. 2000. ANTONIOLLI, Edilar Bento. Análise Comparativa das Técnicas de Medição e Análise das Vibrações, aplicadas na Manutenção Preditiva. 1999. 124 f. Dissertação (Mestrado em Engenharia Mecânica) – Universidade Federal de Santa Catarina, Santa Catarina. AQUINO, Marcelo Braz. Aplicação da Função Densidade de Probabilidade Beta no Diagnóstico de Defeitos de Montagem em Redutores de Engrenagens. 2004. 100 f. Dissertação (Mestrado em Engenharia Mecânica) – Universidade Estadual Paulista, São Paulo. BEZERRA, Roberto de Araújo. Detecção de Falhas em Rolamentos por Análise de Vibração, 2004. 152 f. Tese (Doutorado em Engenharia Mecânica) – Universidade Estadual de Campinas, São Paulo. DIGITALIZAÇÃO de um sinal analógico. Disponível em: <http://www.ceset.unicamp.br/~leobravo/ST%20050/Aula%20Digitaliza%E7%E3o%20de%2 0Sinais%20Anal%F3gicos%20.doc>. Acesso em: 07 Abr 2008. MARÇAL, Rui F. M.; SUSIN, Altamiro A. Detectando falhas incipientes em máquinas rotativas. Revista Gestão Industrial. v. 01, n. 021, p. 87-96, 2005. Disponível em: <http://www.pg.cefetpr.br/ppgep/revista/revista2005/PDF2/Art07Vol1Nr2.pdf>. Acesso em: 07 Abr 2008. MECÂNICA das vibrações, p 114-170. MESQUITA, Alexandre L. A et al. Detecção de falhas em rolamentos usando transformadas tempo-freqüência – Comparação com análise de envelope. Mecânica Computacional. v. 21, p. 1938-1954, 2002. Disponível em: <http://www.cimec.org.ar/ojs/index.php/mc/article/viewFile/1066/1012>. Acesso em: 07 Abr 2008. MESQUITA, Alexandre Luiz. Amarante. Sensores e Filtros. Apostila da disciplina Análise de Vibrações em Máquinas – Universidade Federal do Pará, Belém – PA, 2008. MESQUITA, Alexandre Luiz. Amarante. Processamento digital de sinais. Apostila da disciplina Análise de Vibrações em Máquinas – Universidade Federal do Pará, Belém – PA, 2008. NATIONAL INSTRUMENTS. LabVIEW – User Manual, 2000. 108 NUNES, Marcus Antonius da Costa. Diagnóstico de Defeitos em Rolamentos pela Técnica do Envelope. 1989. 110 f. Dissertação (Mestrado em Engenharia Mecânica, área de concentração Vibrações e Ruído) – Universidade Federal de Santa Catarina, Santa Catarina. PACHOLOK, Mariano. Uso da termografia para avaliação do desalinhamento de eixos de máquinas rotativas: Uma ferramenta auxiliar à análise de vibrações. 2004. 83 f. Dissertação (Mestre em Engenharia Mecânica) – Pontifícia Universidade Católica do Paraná, Paraná. PONCI, Luciano Paiva; CUNHA, Paulo Mario Rodrigues. Previsão de Falha de Rolamentos por Análise Espectral e de Envelope. Disponível em: <http://www.teknikao.com/SDAV%20e%20Rolamentos.pdf>. Acesso em: 10 Dez 2007. SOUSA, Walter dos Santos. Desenvolvimento de um Sistema aplicativo em LabVIEW para o Monitoramento de Máquinas Rotativas com um Módulo de Balanceamento de Rotores. 2005. 115 fp. Trabalho de Conclusão de Curso – Universidade Federal do Pará, Pará. SOUZA, Vitor Hugo Pereira. Projeto e Implementação de uma Bancada de Máquinas Elétricas Rotativas para Monitoração da Geração e Diagnóstico. 2007. 247 f. Trabalho de Conclusão de Curso – Universidade Federal do Pará, Pará. BIBLIOGRAFIA COMPLEMENTAR ALMEIDA, Márcio Tadeu; GÓZ, Ricardo Damião Sales. Curso de Balanceamento de Rotores. Apostila do Curso de Análise de Vibração da FUPAI, Itajubá – MG. BRÜEL & KJAER. Vibration measurement and analysis. 1998 CASTRO, Marcelo Miranda. Bombas Centrífugas – Noções de Operação e Manutenção. Apostila de Gestão da Performance Total – Alunorte - Alumina do Norte do Brasil S.A., Barcarena – PA. MCCONNELL, Kenneth G. Vibration Testing – Theory and Practice. New York: WileyInterscience, 1995. SOEIRO, Newton Sure. Balanceamento de Rotores Rígidos. Apostila da disciplina Vibrações Mecânicas – Universidade Federal do Pará, Belém – PA. 109 APÊNDICE A – Gráficos no domínio do tempo e da freqüência, em aceleração e em velocidade, dos defeitos de desbalanceamento e desalinhamento 1 Desbalanceamento 1.1 Domínio do tempo A figura A1 corresponde ao sinal da máquina sem desbalanceamento e a figura A2 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição axial. Figura A1 Domínio do tempo / sinal sem desbalanceamento / posição axial Figura A2 Domínio do tempo / sinal com desbalanceamento / posição axial 110 A figura A3 corresponde ao sinal da máquina sem desbalanceamento e a figura A4 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial horizontal. Figura A3 Domínio do tempo / sinal sem desbalanceamento / posição radial horizontal Figura A4 Domínio do tempo / sinal com desbalanceamento / posição radial horizontal A figura A5 corresponde ao sinal da máquina sem desbalanceamento e a figura A6 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial vertical. 111 Figura A5 Domínio do tempo / sinal sem desbalanceamento / posição radial vertical Figura A6 Domínio do tempo / sinal com desbalanceamento / posição radial vertical 1.2 Espectro em aceleração A figura A7 corresponde ao sinal da máquina sem desbalanceamento e a figura A8 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição axial. Figura A7 Espectro em aceleração / sinal sem desbalanceamento / posição axial 112 Figura A8 Espectro em aceleração / sinal com desbalanceamento / posição axial A figura A9 corresponde ao sinal da máquina sem desbalanceamento e a figura A10 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial horizontal. Figura A9 Espectro em aceleração / sinal sem desbalanceamento / posição radial horizontal Figura A10 Espectro em aceleração / sinal com desbalanceamento / posição radial horizontal 113 A figura A11 corresponde ao sinal da máquina sem desbalanceamento e a figura A12 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial vertical. Figura A11 Espectro em aceleração / sinal sem desbalanceamento / posição radial vertical Figura A12 Espectro em aceleração / sinal com desbalanceamento / posição radial vertical 1.3 Espectro em velocidade A figura A13 corresponde ao sinal da máquina sem desbalanceamento e a figura A14 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição axial. 114 Figura A13 Espectro em velocidade / sinal sem desbalanceamento / posição axial Figura A14 Espectro em velocidade / sinal com desbalanceamento / posição axial A figura A15 corresponde ao sinal da máquina sem desbalanceamento e a figura A16 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial horizontal. Figura A15 Espectro em velocidade / sinal sem desbalanceamento / posição radial horizontal 115 Figura A16 Espectro em velocidade / sinal com desbalanceamento / posição radial horizontal A figura A17 corresponde ao sinal da máquina sem desbalanceamento e a figura A18 corresponde ao sinal da máquina com desbalanceamento, com o sensor na posição radial vertical. Figura A17 Espectro em velocidade / sinal sem desbalanceamento / posição radial vertical Figura A18 Espectro em velocidade / sinal com desbalanceamento / posição radial vertical 116 2 Desalinhamento 2.1 Domínio do tempo A figura A19 corresponde ao sinal da máquina sem desalinhamento e a figura A20 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição axial. Figura A19 Domínio do tempo / sinal sem desalinhamento / posição axial Figura A20 Domínio do tempo / sinal com desalinhamento / posição axial A figura A21 corresponde ao sinal da máquina sem desalinhamento e a figura A22 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial horizontal. 117 Figura A21 Domínio do tempo / sinal sem desalinhamento / posição radial horizontal Figura A22 Domínio do tempo / sinal com desalinhamento / posição radial horizontal A figura A23 corresponde ao sinal da máquina sem desalinhamento e a figura A24 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial vertical. Figura A23 Domínio do tempo/ sinal sem desalinhamento / posição radial vertical 118 Figura A24 Domínio do tempo / sinal com desalinhamento / posição vertical 2.2 Espectro em aceleração A figura A25 corresponde ao sinal da máquina sem desalinhamento e a figura A26 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição axial. Figura A25 Espectro em aceleração / sinal sem desalinhamento / posição axial Figura A26 Espectro em aceleração / sinal com desalinhamento / posição axial 119 A figura A27 corresponde ao sinal da máquina sem desalinhamento e a figura A28 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial horizontal. Figura A27 Espectro em aceleração / sinal sem desalinhamento / posição radial horizontal Figura A28 Espectro em aceleração / sinal com desalinhamento / posição radial horizontal A figura A29 corresponde ao sinal da máquina sem desalinhamento e a figura A30 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial vertical. 120 Figura A29 Espectro em aceleração / sinal sem desalinhamento / posição radial vertical Figura A30 Espectro em aceleração / sinal com desalinhamento / posição vertical 2.3 Espectro em velocidade A figura A31 corresponde ao sinal da máquina sem desalinhamento e a figura A32 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição axial. Figura A31 Espectro em velocidade / sinal sem desalinhamento / posição axial 121 Figura A32 Espectro em velocidade / sinal com desalinhamento / posição axial A figura A33 corresponde ao sinal da máquina sem desalinhamento e a figura A34 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial horizontal. Figura A33 Espectro em velocidade / sinal sem desalinhamento / posição radial horizontal Figura A34 Espectro em velocidade / sinal com desalinhamento / posição radial horizontal A figura A35 corresponde ao sinal da máquina sem desalinhamento e a figura A36 corresponde ao sinal da máquina com desalinhamento, com o sensor na posição radial vertical. 122 Figura A35 Espectro em velocidade / sinal sem desalinhamento / posição radial vertical Figura A36 Espectro em velocidade / sinal com desalinhamento / posição vertical