Algoritmo para determinação do amortecimento de materiais cerâmicos pela técnica das freqüências naturais de vibração via excitação por impulso

Algorithm to determine the damping of ceramic materials by the impulse excitation technique

Algoritmo para determinação do amortecimento de materiais cerâmicos pela técnica das freqüências naturais de vibração via excitação por impulso

Algorithm to determine the damping of ceramic materials by the impulse excitation technique

A. H. A. Pereira^I; B. Musolino^II; C. D. Maciel^II; J. A. Rodrigues^I

^IGrupo de Engenharia de Microestrutura de Materiais - GEMM, Departamento de Engenharia de Materiais, Universidade Federal de S. Carlos, Rod. Washington Luiz, km 235, S. Carlos, SP 13565-905. [e-mail], [e-mail]

^IIDepartamento de Engenharia Elétrica, Escola de Engenharia de S. Carlos, Universidade de S. Paulo, S. Carlos, SP 13566-590

RESUMO

O amortecimento ou atrito interno é uma das propriedades mais sensíveis de materiais e estruturas tanto em escala macro como microscópica e um bom indicador para a caracterização de danos estruturais, sendo particularmente sensível à presença de trincas e micro-trincas. O objetivo deste trabalho foi desenvolver um algoritmo que permitisse incorporar a caracterização do amortecimento à técnica das freqüências naturais de vibração via excitação por impulso, já bem estabelecida pela norma ASTM E1876 e amplamente empregada na determinação não-destrutiva dos módulos elásticos de materiais cerâmicos. Foi desenvolvido um algoritmo para a determinação do amortecimento fundamentado na análise de sinal no domínio de tempo-freqüência, no método do decremento logarítmico e no modelo viscoelástico equivalente. O algoritmo foi validado empregando um gerador de ondas de formato arbitrário e aplicado na caracterização de dois materiais, sendo um concreto refratário de alta alumina, com e sem dano por choque térmico, e um aço SAE 1020. O algoritmo foi eficiente ao determinar a frequência e o amortecimento de sinais arbitrários e dos corpos de prova, e em limitar a análise ao modo de vibração de interesse, minimizando assim a interferência de outros modos. O algoritmo permitiu integrar com eficiência o cálculo do amortecimento à técnica de excitação por impulso.

Palavras chave: algoritmo, amortecimento, materiais cerâmicos, freqüências naturais, ASTM E1876.

ABSTRACT

Damping is one of the most sensitive properties of materials and structures, in macro and microscopic scale, which makes it particularly sensitive to the presence of cracks and micro-cracks and a good parameter to the characterization of structural damages. The aim of this work was to develop an algorithm to allow the incorporation of the damping characterization to the natural frequencies of vibration technique (by impulse excitation), already well established for the non-destructive determination of elastic moduli by ASTM E1876 and correlated standards for ceramic materials. It was developed an algorithm for determining the damping based on the time-frequency domain signal analysis, logarithmic decrement method and on the equivalent viscoelastic model. The algorithm was validated using an arbitrary wave form signal generator and applied to characterize two materials, a high alumina refractory castable with and without thermal shock damage, and a SAE-1020 steel. The algorithm was effective to recovery the frequencies and respective damping of the arbitrary simulated signals and to characterize the specimens, limiting the analysis to the vibration mode of interest, thereby minimizing the modes interference. This algorithm allowed the successful integration of the damping characterization to the impulse excitation technique.

Keywords: algorithm, damping, ceramic materials, natural frequencies, ASTM E1876.

INTRODUÇÃO

O amortecimento ou atrito interno é uma das propriedades mais sensíveis de materiais e estruturas, tanto em escala macro quanto microscópica [1], sendo particularmente sensível à presença de trincas e micro-trincas [2]. Além da aplicação clássica no estudo de metais e nas áreas de engenharia mecânica e engenharia civil, a caracterização do amortecimento também vem sendo empregada no estudo de materiais cerâmicos refratários para a avaliação do dano por choque térmico [3, 4]. A tensão mecânica induzida pelo gradiente de temperatura do choque térmico provoca a nucleação e propagação de micro-trincas e trincas que degradam as propriedades mecânicas do material determinando em grande parte a sua vida útil [5, 6]. A nucleação e evolução destas micro-trincas e trincas podem ser monitoradas com a caracterização do amortecimento [3, 4, 7]. A caracterização do amortecimento também pode ser empregada para o controle de qualidade em linhas de produção. O amortecimento é o fenômeno pelo qual a energia mecânica de um sistema dinâmico é dissipada [8]. O amortecimento de um sistema ou material pode ser classificado de três formas principais: interno, estrutural e fluídico. O interno está associado aos defeitos na microestrutura, efeitos anelásticos e a efeitos termoelásticos causados por gradientes locais de temperatura. Já o estrutural está associado a perdas de energia por atrito em juntas, parafusos e articulações semi-rígidas. Por último, o fluídico ocorre por resistência ao arraste em meio fluídico, por exemplo, a conversão de energia cinética de um pêndulo em energia térmica para o ar. Neste trabalho o interesse é pela avaliação do amortecimento interno ou atrito interno [8].

Existem diversos métodos para a determinação do amortecimento, sendo os mais utilizados em análise modal os do decremento logarítmico e da largura de banda de meia potência [9, 10]. A escolha do método depende da faixa do amortecimento, da frequência de vibração, da praticidade requerida e do hardware disponível [8, 11]. O método do decremento logarítmico é mais adequado para baixo e médio amortecimento e consiste no cálculo a partir da atenuação da vibração do material ou estrutura em função do tempo após uma excitação por impulso. O método da largura de banda de meia potência é mais adequado para médio e alto amortecimento. O método do decremento logarítmico é o mais prático, contudo demanda de uma capacidade superior de processamento de sinal. Ambos os métodos consideram um modelo para os cálculos, normalmente o modelo de amortecimento viscoelástico equivalente, em que o decremento da amplitude de vibração em função do tempo é tido como exponencial e o amortecimento não depende da freqüência de vibração [8, 11, 12]. A técnica das freqüências naturais de vibração via excitação por impulso consiste essencialmente em excitar, através de um leve e breve impacto mecânico, o modo de ressonância flexional e/ou torcional fundamental de uma amostra em formato de barra, e em calcular os módulos elásticos dinâmicos [13] empregando as equações de Pickett [14] a partir das frequências encontradas na resposta acústica, das dimensões e da massa da amostra. A pancada deve ser curta o suficiente para que apresente um espectro plano de frequência na faixa de interesse. O corpo de prova reage à pancada como um filtro passa-faixas revelando as suas freqüências naturais de vibração através da emissão de uma resposta acústica transitória, que também contém a informação para a determinação do amortecimento pelo método do decremento logarítmico. A excitação do modo de vibração desejado é realizada impondo-se condições de contorno mecânicas e de excitação adequadas, o que consiste em apoiar o corpo de prova nas linhas nodais de vibração e excitá-lo onde a amplitude de vibração será máxima [15]. Contudo, dificilmente se consegue excitar um único modo de vibração, mesmo porque quando é desejado se obter a razão de Poisson é necessário excitar modos flexionais e torcionais. Em função desta dificuldade, a aplicação do método do decremento logarítmico padrão (no domínio do tempo) possui um componente de incerteza significativa associada à contribuição para a resposta acústica de modos de vibração menos favorecidos pelas condições de contorno, visto que estes modos apresentam maior interação com o suporte de corpo de prova. Estes modos menos favorecidos também aumentam a incerteza dos módulos elásticos quando o cálculo da frequência realizado pelo sistema de medição é feito através do período do sinal captado, como acontece em métodos fundamentados em zero-crossing [16].

O objetivo deste trabalho foi desenvolver um algoritmo integrável à técnica de excitação por impulso que contornasse os problemas mencionados. O desenvolvimento do algoritmo foi fundamentado na análise de sinal no domínio tempo-freqüência, através da qual o decremento logarítmico da amplitude é obtido e avaliado exclusivamente para a frequência correspondente ao modo de vibração mais favorecido pelas condições de contorno, garantindo, assim, o cálculo do amortecimento com maior precisão. Este recurso é particularmente bem vindo no caso de materiais cerâmicos com alto dano microestrutural, por exemplo, concretos refratários submetidos a choques térmicos e concretos estruturais submetidos a cargas elevadas, nos quais o dano aumenta o amortecimento de forma expressiva o que degrada a razão sinal/ruído e aumenta o acoplamento entre os modos de vibração. A conveniência de unir o método do decremento logarítmico no domínio tempo-freqüência à técnica das frequências naturais via excitação por impulso reside no fato da resposta acústica conter a informação necessária para a determinação do amortecimento, e no fato da geometria do corpo de prova e as condições de contorno do ensaio serem padronizadas [13]. Adicionalmente, a incorporação da medição do amortecimento não demanda de procedimentos experimentais adicionais ou altera a duração do ensaio, mantendo assim a praticidade e rapidez da técnica.

MATERIAIS E MÉTODOS

Modelo viscoelástico equivalente

O modelo viscoelástico equivalente [8] é o alicerce do algoritmo proposto. Neste modelo a força de amortecimento é proporcional à velocidade. Um exemplo clássico de aplicação do modelo é o sistema massa-mola-amortecedor com um grau de liberdade descrito pela equação de segunda ordem [17]:

em que:

sendo x a posição em função do tempo t, ζ a taxa de amorte-cimento ou apenas amortecimento, ω₀ a frequência natural de vibração, m a massa e κ a constante elástica da mola. Assumindo a solução x (t) = e^jyt para e equação (A) (j= número complexo unitário), chega-se a γ como dado por:

O comportamento do sistema massa-mola-amortecedor depende da solução de γ, que pode ter duas soluções reais ( ζ > 1), uma solução real (ζ = 1) ou duas soluções complexas (0 < ζ < 1). Os casos sobreamortecido (ζ > 1) e criticamente amortecido (ζ = 1) são não-oscilatórios e irrelevantes para este trabalho. Quando 0 < ζ < 1 o sistema é sub-amortecido e a solução da equação (F) pode ser escrita como:

em que A₀ é a amplitude inicial de vibração, φ é a fase inicial da vibração e ω é a frequência natural amortecida descrita por:

Esse modelo também é conhecido como sistema linear amortecido com um grau de liberdade [8, 17]. Generalizando para um sistema com K graus de liberdade, a amplitude de vibração de um determinado ponto da amostra, quando captada por um sensor, é dada por:

em que A_i, ζ _i, ω_0i, ω_i, φ_i são, respectivamente, amplitude inicial, amortecimento, frequência natural de vibração, frequência natural amortecida e fase inicial do i-ésimo grau de liberdade ou modo de vibração naquele ponto. O termo R_wn é um ruído descorrelacionado do sinal captado ou gerado pelo sensor, assumido como sendo ruído branco.

O modelo viscoelástico equivalente consiste em em-pregar as equações do modelo viscoelástico para sistemas que apresentam decremento exponencial da amplitude de vibração, que é um comportamento típico de sistemas viscoelásticos. Contudo estes sistemas também podem apresentar mecanismos de amortecimento não-viscosos, por exemplo, mecanismos anelásticos dependentes da frequência ou da amplitude de vibração. Esta simplificação por um lado desconsidera a possibilidade de ocorrência de mecanismos não-viscosos, mas por outro viabiliza caracterizações experimentais comparativas para corpos de prova com faixa de frequência de ressonância e dimensões similares [18].

Algoritmo desenvolvido

Conforme ilustrado no fluxograma da Fig. 1, o algoritmo desenvolvido foi estruturado em 05 passos principais visando aproveitar a padronização das condições de contorno mecânicas do ensaio para simplificar o processamento de sinal. Este aproveitamento se deu com a determinação do decremento logarítmico baseando-se no acompanhamento da amplitude das freqüências encontradas no primeiro segmento do sinal, em vez de se aplicar filtros digitais integradores ao redor destas freqüências ou de utilizar transformadas especiais como é usualmente realizado em análise modal [11, 12].

O primeiro passo consiste na aquisição do sinal e na sua segmentação. A digitalização é realizada através de uma placa de áudio ou pelo hardware onde o algoritmo estiver embarcado. A taxa de amostragem da digitalização deve ser pelo menos o dobro da frequência mais alta que se deseja analisar. A segmentação do sinal é realizada de acordo com a equação:

em que h é o índice do segmento, H é espaçamento (offset) entre o início de um novo segmento e o início do segmento anterior, e n_s é o comprimento de cada segmento em número de pontos. A razão entre H e n_s determina a sobreposição ou espaçamento entre os segmentos.

O segundo passo consiste na realização da transformada rápida de Fourier (FFT) [19] do primeiro segmento do sinal (s₁) com n_s amostras. É no primeiro segmento que as componentes de sinal apresentam a maior amplitude, visto que a resposta acústica do corpo de prova é um sinal transitório com decremento exponencial. Usando a FFT e zero-padding [19], é encontrado o espectro de magnitudes S₁ (e^jw) do segmento de sinal s₁ (k).

O terceiro passo consiste na detecção dos picos do espectro de magnitudes |S₁(e^jw)| e no cálculo dos módulos elásticos. Para isso é calculado o valor médio da derivada do segmento, através da equação:

com n sendo a indexação dos valores de frequência do espectro |S₁(e^jw)|. São consideradas possíveis frequências de ressonância todos os picos que possuírem em sua borda de subida inclinação maior que a derivada média.

O quarto passo consiste na análise tempo-frequência do sinal com a montagem de uma matriz que contém a informação de como a amplitude de cada frequência detectada varia com o tempo desde o momento da excitação do corpo de prova. Para o cálculo da FFT dos segmentos s_h, todos são multiplicados por uma janela flat-top [20] e submetidos ao zero-padding.

O quinto e último passo consiste no ajuste de uma curva sobre a evolução da amplitude de cada frequência ω_i registrada na matriz tempo-freqüência, e no cálculo do respectivo amortecimento a partir do decremento logarítmico da região onde ocorre o decremento exponencial do sinal. Para a realização do ajuste de curva é empregada a equação (N) com o método não-linear por mínimos quadrados proposto por Levemberg-Marquardt [21, 22]. Este método permite selecionar a região a ser priorizada pelo ajuste em função da amplitude através do parâmetro s que pode assumir valores entre 0,1 e 1.

Na equação (N) empregada pelo método Levemberg-Marquardt [21, 22] a_i (t) é a amplitude RMS (Root Mean Square) da frequência ω_i em função do tempto t, δ_i é o decremento logarítmico correspondente e A_wn o piso de ruído. O decremento logarítmico δ_i pode ser visualizado como o coeficiente angular da queda de amplitude em escala logarítmica em função do tempo. Determinado δ_i , o amortecimento ζ _i é calculado através da relação:

A execução destes cinco passos é muito rápida e tem a duração determinada essencialmente pela duração da resposta acústica emitida pelo corpo de prova.

Aparato experimental

Para a aquisição de sinal foi empregado neste trabalho um sensor acústico capacitivo com largura de banda de 20 Hz a 20 kHz. Para a digitalização do sinal foi empregada uma placa de aquisição e processamento de sinais com taxa de aquisição de até 96 ksp/s (ksamples/second - quilo amostras/segundo) e quantização em 24 bits.

As caracterizações foram realizadas utilizando o suporte de corpo de prova de precisão para barras retangulares modelo Sonelastic^® SP-B (ATCP Engenharia Física) que impõem as condições de contorno ótimas para o modo de vibração flexional fundamental. Neste suporte a amostra é apoiada exatamente sobre os nós de vibração flexional fundamental por fios metálicos na posição 0,224L, em que L é o comprimento da barra. O ajuste da posição é realizado automaticamente por um sistema mecânico fundamentado no comprimento total do corpo de prova (Fig. 2). É importante que o sistema de excitação e captação do sinal exerça a mínima influência possível sobre a vibração flexional fundamental. A excitação mecânica é provida pelo atuador eletromagnético Sonelastic^® IED.

EXEMPLO DE APLICAÇÃO

O algoritmo desenvolvido programado em MATLAB® foi empregado em conjunto com o aparato experimental para a determinação simultânea do amortecimento, do módulo de Young, do módulo de cisalhamento e da razão de Poisson de um concreto refratário de alta alumina e de um aço SAE 1020. O concreto refratário foi caracterizado com e sem dano por choque térmico. Os corpos de prova caracterizados foram preparados no formato de barras retangulares e estão mostrados na Fig. 3. Este formato é padrão para ensaios mecânicos de resistência mecânica, portanto muito comum.

Na Tabela I são apresentadas as dimensões e a massa dos corpos de prova com as respectivas incertezas. Os corpos de concreto foram conformados por vibração em moldes metálicos e tiveram a face mais grosseira retificada (a mesma face que está sendo mostrada pela foto da Fig. 3; a amostra de aço teve as quatro faces maiores retificadas.

Aquisição e segmentação do sinal

Conforme o fluxograma da Fig. 1, o 1º passo do algoritmo é a aquisição e a segmentação do sinal, que foram realizadas de acordo com os parâmetros da Tabela II.

Os parâmetros de aquisição e segmentação foram ajustados para que o número total de segmentos se mantivesse em torno de 250 para os três corpos de prova. O tempo de aquisição foi ajustado para que o sinal fosse captado até pelo menos no momento em que se torna menor que o piso de ruído. Desta forma, tem-se toda a informação possível para o posterior ajuste de curva sem se sobrecarregar o processamento computacional. Nota-se que a duração do sinal é muito maior para o aço, devido ao seu baixo amortecimento comparativamente ao do concreto.

Na Fig. 4 é mostrado um exemplo de sinal adquirido, mais especificamente o sinal da amostra de concreto sem dano por choque térmico (C1).

Detecção das freqüências e cálculo dos módulos

O segundo e o terceiro passo consistem na obtenção do espectro de freqüências do sinal contido no primeiro segmento, na detecção automática das freqüências presentes e no cálculo dos módulos elásticos. Na Fig. 5 é mostrado um exemplo de espectro, mais especificamente da amostra de concreto sem dano por choque térmico (C1), com as frequências detectadas na faixa de interesse destacadas com quadrados. A amplitude de referência para as escalas em dB é amplitude máxima aceita pela placa de aquisição (2 V pico-a-pico).

Na Tabela III são apresentadas as freqüências detectadas na faixa de interesse (até 15 kHz); acima de 15 kHz se encontram s freqüências harmônicas destes corpos e estas não são necessárias para os cálculos dos módulos elásticos.

Com o conhecimento da massa, das dimensões dos corpos de prova (Tabela I) e das suas frequências de ressonância (Tabela III são calculados os módulos elásticos a partir das equações de Pickett [14], que também permitem a identificação do modo de vibração a que pertence cada frequência a partir da razão entre elas. Estas razões dependem essencialmente das proporções geométricas do corpo de prova. No caso de barras retangulares de secção quadrada com aresta 6 vezes menor que o comprimento, a primeira frequência (1º pico) é o modo flexional fundamental e a segunda (2º pico) o modo torcional fundamental. Na Tabela IV são apresentados os resultados obtidos para os módulos elásticos. As incertezas apresentadas nos resultados da Tabela IV foram calculadas a partir da propagação da incerteza das dimensões, da massa e das equações de Pickett.

O quarto passo é a análise tempo-frequência do sinal, que consiste na obtenção do espectro de freqüências via FFT para cada um dos segmentos do sinal, de modo que seja possível o acompanhamento do decremento logarítmico de cada frequência presente, e a partir desta informação, realizar o cálculo dos respectivos amortecimentos. Na Fig. 6 é apresentado como exemplo o espectrograma da resposta acústica do corpo de prova de aço. O perfil do decremento é visto como uma reta e não como uma exponencial porque a escala de amplitude é logarítmica; o parâmetro taxa de decremento logarítmico (di) pode ser visualizado como o coeficiente angular da rampa formada pelo decaimento do sinal.

O quinto passo consiste no ajuste de curva sobre o decremento da amplitude da frequência de interesse em função do tempo, empregando o algoritmo não-linear por mínimos quadrados Levemberg-Marquardt [16, 17]. Esta frequência corresponde normalmente ao modo de vibração mais favorecido pelas condições de contorno do suporte de corpo de prova, no caso do aparato experimental empregado para este trabalho, o modo flexional fundamental. Na Fig. 7 é apresentado um exemplo de ajuste de curva; o coeficiente angular da região de queda linear do ajuste é o parâmetro taxa de decremento logarítmico.

Na Tabela V são apresentados os resultados encontrados para o amortecimento do modo flexional fundamental dos concretos e do aço empregando o algoritmo desenvolvido. Para efeito de comparação, é apresentado também o resultado realizando-se o ajuste sobre a amplitude do sinal no domínio do tempo.

Nota-se que o amortecimento obtido com o algoritmo desenvolvido é menor que o obtido do decremento do sinal no domínio do tempo; o que era esperado, visto que este elimina a influência indesejável dos modos de vibração menos favorecidos pelas condições de contorno mecânicas, mas que estão inevitavelmente presentes e são necessários para o cálculo simultâneo do módulo de cisalhamento e da razão de Poisson.

FAIXA DE MEDIÇÃO E VALIDAÇÃO

Faixa de medição

Os extremos para a medição do amortecimento pelo método do decremento logarítmico são 0 e 1. O extremo igual a 0, que é o menor valor possível de amortecimento, não ocorre dissipação de energia e o sistema vibra indefinidamente com amplitude constante. No extremo igual a 1 o sistema é criticamente amortecido e retorna ao equilíbrio após a excitação sem completar um ciclo de oscilação; a medição do amortecimento com o algoritmo neste extremo, ou muito próximo a ele, não é possível porque a frequência do sinal não é detectável. Portanto, determinar a faixa de medição consiste essencialmente em estimar o quão próximo é possível chegar do limite onde o amortecimento é igual a 1, em função da frequência e amplitude do sinal, e dos parâmetros da aquisição e do processamento. Como o algoritmo é fundamentado no ajuste de uma curva sobre a evolução da amplitude da frequência de interesse ω_i ao longo dos segmentos S_h, o amortecimento máximo mensurável ζ _{i max} para uma determinada frequência ω_i é limitado pela razão entre a duração do sinal e a duração do primeiro segmento mais o espaçamento entre este e o (s) segmento (s) subsequente (s). Esta razão deve ser maior ou igual a 1. Se esta razão for inferior a 1, o que significa que duração do sinal não foi suficiente para a obtenção de pelo menos dois pontos para o preenchimento da matriz tempo-frequência, não é possível o ajuste da curva para a determinação do decremento logarítmico, e consequentemente, o cálculo do amortecimento. Adi-cionalmente, os segmentos S_h não podem ser reduzidos indefinidamente porque a redução degrada a precisão na determinação da frequência e da amplitude, visto que a precisão em frequência é diretamente proporcional à taxa de amostragem e inversamente proporcional ao número de pontos do segmento.

Da mesma forma que o cálculo do amortecimento, a faixa mensurável de amortecimento depende da frequência ω_i (equação (O) ). Portanto é conveniente basear o equacionamento de ζ _{i max} em termos do decremento logarítmico máximo δ_{i max}:

e o decremento logarítmico máximo δ_{i max} em termos da amplitude do sinal relativa ao piso de ruído e da duração do primeiro segmento mais o (s) intervalo (s) de tempo entre o primeiro e os N-1 segmento (s) conseguinte (s):

nas equações (Q) e (P) A₀ é a amplitude inicial do sinal (em dB), A_wn o piso de ruído (em dB), f_a a taxa de aquisição ou frequência de amostragem do sinal e N o número mínimo de pontos que se deseja obter para realizar o ajuste da curva (N deve ser > 2). A inserção da constante 1/20 . log (e) na equação (Q), em que "e" é o número de Neper, deve-se à compatibilização da amplitude (em dB na base 10) com a base "e" do método do decremento logarítmico.

Com o amortecimento máximo mensurável equa-cionado, o próximo passo para estimar a faixa de medição é realizar o cálculo para as condições experimentais típicas. Na Tabela VI são apresentados valores típicos para os parâmetros determinantes de δ_{i max} nas situações em que o algoritmo desenvolvido será inicialmente aplicado e para que possa ser facilmente embarcado em uma eletrônica com processamento de sinal voltado para áudio.

Conhecendo δ_{i max} , pode-se calcular o amortecimento máximo em função da frequência empregando-se a equação (Q); quanto maior a frequência, menor o amortecimento máximo mensurável. Para a frequência de 5 kHz, com δ_{i max} = 2.250 dB/s , o amortecimento máximo é da ordem de 10^-2, que é um valor alto para materiais cerâmicos em temperatura ambiente. Para 500 Hz o máximo atinge 10^-1. Para o teste de materiais cerâmicos, por exemplo, concretos refratários e abrasivos, a faixa de frequência usual está entre 2 e 6 kHz e de amortecimento entre 10^-5 e 10^-2. Contudo para amostras pequenas (<10 cm), a frequência pode atingir valores superiores.

Portanto, para frequências de até 5 kHz, a faixa de medição do amortecimento pelo algoritmo proposto para materiais cerâmicos é de aproximadamente de 0 a 10^-2, sendo este valor inversamente proporcional à frequência de acordo com a equação (Q) e com a ressalva de que esta faixa pode variar em função dos parâmetros da aquisição e do processamento do sinal.

Validação

O algoritmo desenvolvido foi testado com a simulação de medições empregando sinais transitórios sintetizados com frequência e amortecimento bem definidos. Para sintetizar os sinais foi utilizado um gerador de forma de onda arbitrária modelo HP-33120A (Hewlett-Packard). Foram escolhidos amortecimentos próximos ao limite de medição, entre aproximadamente 10^-1 e 10^-3 p/ N=3, na faixa de frequência de 200 Hz a 20 kHz. A amplitude do sinal simulado com relação ao piso de ruído foi ajustada para 44 dB. Também foram simuladas medições com amortecimento igual a zero na mesma faixa de frequência. Os resultados mostraram que para sinais com valores de amortecimento de até 10^-2 e frequências entre 500 e 10 kHz, o erro médio e o desvio padrão da determinação do amortecimento é inferior a 1% e para a determinação da frequência inferior a 0,1%. Para condições extremas, tendendo ao limite mensurável teórico, o erro médio e o desvio padrão foram de até 4% para o amortecimento e 8% para a frequência. Estes valores de erro estão compatíveis ou inferiores ao descritos na literatura para métodos de medição similares [12].

CONCLUSÕES

Foi desenvolvido um algoritmo fundamentado no método do decremento logarítmico, no modelo viscoelástico equivalente e na análise tempo-frequência, para a determinação do amortecimento de materiais cerâmicos; e este incorporado à técnica das frequências naturais de vibração via excitação por impulso. O algoritmo, em associação com o aparato experimental empregado, se mostrou preciso na análise de sinais sintetizados com frequência e amortecimento nas faixas típicas das amostras de materiais cerâmicos (frequência de até 6 kHz e amortecimento de até 10^-2), na caracterização de concreto refratários de alta alumina e na caracterização de um aço SAE 1020. Graças à análise tempo- freqüência, o algoritmo foi eficiente na medição isolada do amortecimento de cada frequência presente nas respostas acústica e apresentou erros compatíveis ou inferiores aos descritos na literatura para sistema de análise modal que compartilham os mesmos princípios. É interessante notar que o aumento percentual do amortecimento do concreto caracterizado devido ao choque térmico foi aproximadamente uma ordem de grandeza superior à queda percentual dos módulos elásticos (~ +350% vs. -40%), o que evidencia a sensibilidade deste parâmetro para a avaliação do dano por choque térmico nesta categoria de materiais cerâmicos. O choque térmico induziu a formação de trincas e micro-trincas que enfraqueceram o material e elevaram o atrito interno. A integração do algoritmo não comprometeu a praticidade ou alterou os procedimentos da técnica das frequências naturais de vibração via excitação por impulso. Com esse avanço, a técnica passa a informar um parâmetro importante dos materiais cerâmicos a partir do mesmo procedimento experimental usado até então, sem a necessidade de empregar aparatos experimentais adicionais específicos para análise modal.

AGRADECIMENTOS

Ao CNPq (Procs. 301073/2006-6, 303061/2009-0 e 141868/2008-2) e à FINEP (Proj. 4028-6) pelo apoio financeiro, e às empresas Elfusa e Almatis pela doação da matéria prima.

REFERÊNCIAS

[1] B. J. Lazan, Damping of materials and members in structural mechanics, Pergamon Press, Oxford (1968).

[2] R. Dieterle, H. Banchmann, Int. Assoc. Bridge and Structural Eng. Report Working Com. 34(1981)69.

[3] J. A. Coppola, J. Am. Ceram. Soc. 56,4(1973)214.

[4] T. Tonnesen, R. Telle, Ceram. Forum Int. 84,9(2007)132.

[5] D. P. Hasselman, J. Am. Ceram. Soc. 46,1 (1963)535.

[6] W. D. Kingery, J. Am. Ceram. Soc. 38,3(1955)3.

[7] S. H. Chowdhury, PhD Thesis, Griffith University, Australia (1999).

[8] C. W. Silva, Vibration damping, control, and design, Taylor & Francis Group, Vancouver, Canada (2007) 634p.

[9] A. Puskar, Internal friction of materials, Cambridge Int. Sci. Publ., Cambridge, Inglaterra (2001) 342p.

[10] S. Gade, H. Herlufsen, H. Konstantin-Hansen, Application notes BO 0177-11, Brüel&Kraer (2002).

[11] Damping Measurements, Technical Rev. 2, Brüel&Kraer (1994).

[12] Digital Filter Techniques vs. FFT Techniques for Damping Measurements, Tech. Rev. 1, Brüel&Kraer (1994).

[13] Standard Test Method for Dynamic Young's Modulus, Shear Modulus, and Poisson's Ratio by Impulse Excitation of Vibration, ASTM E 1876-07, ASTM Int. (2007).

[14] G. Pickett, Proc. Am. Soc. Test. Mater. 45(1945)846.

[15] R. Morrel, Measuring elastic properties of advanced technical ceramics - A review, UK Nat. Phys. Lab. Report 42(1996).

[16] Operating Instructions for the GrindoSonic MK5 "Industrial" Instrument, J. W. Lemmens, Inc.

[17] D. Thorby, Structural Dynamics and Vibrations in Practice - An Engineering Handbook, Elsevier Ltd., Oxford, Inglaterra (2008)420p.

[18] A. H. A. Pereira, G. M. Fortes, B. Schickle, T. Tonnesen, B. Musolino, C. D. B. Maciel, J. A. Rodrigues, Cerâmica 56(2010)311.

[19] A.V. Oppenheim, R. W. Schafer, J. R. Buck, Discrete-Time Signal Processing, 2^nd Ed., Prentice-Hall, Inc., Upper Saddle River, EUA (1998) 870p.

[20] A. F. Harvey, M. Cerna, Application Note 041, Nat. Instr. Corp., EUA (1993).

[21] K. Madsen, H. B. Nielsen, O. Tingleff, Methods for Non-Linear Least Squares Problems, 2^nd Ed. Informatics and Mathematical Modeling, Dinamarca (2004).

[22] L. Finschi, Technical Report, Institute for Operations Research, ETH Zürich (1996), http://www.ifor.math.ethz.ch/staff/finschi/LevMar/LevMar.pdf

(Rec. 28/06/2011, Ac. 22/08/2011)

Datas de Publicação

• Publicação nesta coleção
  24 Jul 2012
• Data do Fascículo
  Jun 2012

Histórico

• Recebido
  28 Jun 2011
• Aceito
  22 Ago 2011