Latest revision as of 12:10, 22 July 2016

Resumo

Neste artigo é proposto um modelo de elementos finitos simplificado para representação de juntas rebitadas em análise dinâmica de estruturas. O rebite é modelado com elementos que combinam mola e amortecedor. Vários modelos numéricos são estudados com diferentes quantidades de rebites (um, 2, 3 e 5) e diferentes quantidades de elementos mola‐amortecedor (4, 6, 8, 12, 16 e 20) por cada rebite. Amostras constituídas por 2 placas de alumínio ligadas entre si por diferentes quantidades de rebites (um, 2, 3 e 5) são construídas e submetidas a análise modal experimental, a fim de serem conhecidas as suas características modais – frequências naturais e formas naturais de vibração. É usada uma metodologia de melhoramento de modelos de elementos finitos, também conhecido na linguagem anglo‐saxónica por updating , com o fim de obter a melhor representação numérica da junta rebitada relativamente às amostras experimentais. Uma avaliação da melhor representação numérica da junta rebitada é efetuada com base na comparação entre as frequências naturais experimentais e as frequências naturais numéricas, para cada modo de vibração, após aplicação do processo de melhoramento de elementos finitos. Mostra‐se que o modelo de junta composto por 8 e 12 elementos mola‐amortecedor por rebite consegue uma melhor aproximação ao comportamento dinâmico observado experimentalmente. Um valor numérico da constante de rigidez é obtido para as juntas em estudo.

Abstract

In this paper, a simplified finite element model to represent a riveted lap joint in structural dynamic analysis field is proposed. The rivet is modeled by spring‐damper combination elements. Several numerical models are studied with different quantities of rivets (1, 2, 3 and 5) and spring‐damper combination elements (4, 6, 8, 12, 16 and 20) per rivet. Experimental samples of two aluminum material plates connected by different quantities of rivets (1, 2, 3 and 5) are built and tested in order to be known its modal characteristics ‐ natural frequencies and mode shapes. A finite element model updating methodology is used in order to get the best numerical riveted lap joint representation relatively to the experimental samples. An evaluation of the best numerical riveted lap joint is carried out based on the comparison between the experimental and numerical resonance frequencies after updating. It is shown that the riveted lap joints composed by eight and twelve spring‐damper combination elements per rivet have the best representation. A stiffness constant value is obtained for the riveted lap joints in study

Palavras‐chave

Junta rebitada ; Melhoramento de modelos de elementos finitos ; Análise modal ; Otimização

Keywords

Revited Lap Joint ; Finite Element Model Updating ; Modal Analysis ; Optimization

1. Introdução

As juntas rebitadas possuem ampla aplicação em diversos setores industriais, principalmente na indústria automóvel e aeroespacial. As estruturas envolvidas nesses setores de atividade possuem normalmente um largo número de ligações rebitadas. Para os projetistas desse tipo de estruturas não é prático modelar detalhadamente o rebite em elementos finitos, devido à necessidade de usar um grande número de graus de liberdade, aumentando o custo computacional da análise.

Alguns trabalhos que descrevem o comportamento estático e o comportamento à fadiga de juntas rebitadas, usam uma modelação simplificada [1] , [2] , [3]  and [4] . Todos eles utilizam elementos de mola‐amortecedor ou elementos de viga para efetuar a ligação entre os nós das placas a unir e assim simular o rebite. Contudo, a ligação rebitada não fica totalmente definida se o valor da constante de rigidez também não for definido. O valor da constante de rigidez da junta nunca é claramente explícito. Para aumentar a precisão do modelo numérico de uma junta rebitada, o valor da constante de rigidez deve ser definido, uma vez que este parâmetro modal tem grande influência nas frequências naturais e nas formas naturais de vibração do sistema.

Um trabalho semelhante ao proposto pelos autores, mas não envolvendo juntas rebitadas, descreve uma avaliação de modelos simplificados de elementos finitos para juntas soldadas por pontos [5] . Nesse trabalho, 4 tipos de modelos simplificados para simular aquele tipo de soldadura são avaliados pelos autores: o modelo de barra rígida múltipla ou multiple rigide‐bar (MRB), o modelo de casca e barra‐rígida ou rigid bar‐rigid shell (RB‐RSH), o modelo sólido tridimensional ou solid nugget (SN) e o modelo de barra rígida ou rigid bar (RB).

Um dos 4 modelos apresentados tem particular interesse para aplicação neste trabalho, o modelo MRB. No modelo MRB, 2 chapas são conectadas por vários elementos de barra de modo a simular a soldadura por pontos. Este modelo tem a vantagem de apresentar vários elementos de barra, o que permite obter mais establidade na junta de ligação devido à quantidade de graus de liberdade que são restringidos.

O modelo SN é composto por elementos sólidos e, portanto, mais complexo de construir e consequentemente de análise computacional mais demorada.

Para determinar o valor da constante rigidez da mola que simula o rebite, uma técnica de melhoramento de modelos de elementos finitos, também conhecida por updating , é utilizada. Essa técnica está assente num processo de otimização baseado em sucessivas iterações, sendo por isso denominado por método iterativo ou método indireto [6] . Os métodos indiretos podem ser divididos em 3 grupos principais: o método da sensibilidade [7] , o método de função de resposta [8] e o método no domínio do tempo [9] . Neste trabalho é utilizada uma metodologia de melhoramento de modelos de elementos finitos desenvolvida com base no método da sensibilidade [10] .

A metodologia, desenvolvida em Matlab^® , trabalha com um otimizador que controla os parâmetros a otimizar no modelo numérico. Por sua vez, o otimizador está dotado com um solver de procura global estocástico (GlobalSearch ) e usa o algoritmo Sequential Quadratic Programmning (SQP) para encontrar o mínimo global da função objetivo [11] . Neste trabalho o parâmetro a otimizar é a constante de rigidez da mola do elemento mola‐amortecedor. O processo de melhoramento de modelos de elementos finitos é descrito na secção 4.

É efetuada a análise modal experimental, como descrito na secção 5, de modo a conhecer o comportamento dinâmico dos modelos físicos com junta rebitada. Os modelos são constituídos por 2 chapas de alumínio unidas por diferentes quantidades de rebites (um, 2, 3 e 5).

Os modelos de elementos finitos são construídos em ANSYS^®Parametric Design Language (APDL), conforme descrito na secção 6. Um valor inicial para a constante de rigidez é atribuído aos elementos mola‐amortecedor. O amortecimento em materiais metálicos é considerado insignificante, pelo que pode ser desprezado, não sendo assim considerado no elemento mola‐amortecedor. De modo a justificar essa opção foram determinados os fatores de amortecimento ξ para todos os modos de vibração obtidos experimentalmente de um dos casos de estudo apresentado. Foi usado o método do ponto de meia potência para obtenção dos fatores de amortecimento. Para cada modo de vibração é também apresentado o respetivo valor do fator de qualidade Q . Considera‐se que para valores de ξ  < 0.1 o amortecimento do sistema é muito pequeno e para valores de ${\textstyle Q>{\frac {1}{2}}}$ o sistema é do tipo subamortecido [12] .

Com o programa de melhoramento de modelos de elementos finitos é estimado o valor da constante rigidez da mola, sendo os resultados apresentados na secção 7.

O modelo de elementos finitos da junta rebitada que melhor representa os modelos físicos é assim obtido, como se conclui na secção 8.

2. Problema teórico

Para um sistema dinâmico, com vários graus de liberdade, a equação de movimento que governa a resposta dinâmica do sistema, quando sujeito a excitações externas, é dada sob a forma da matriz generalizada [13]

onde, M , C e ${\textstyle K\in \mathbb {R} ^{n\times m}}$ , sendo n   o número de graus de liberdade, denotam, respetivamente, a matriz de massa, amortecimento e rigidez do sistema, ${\textstyle e{\mbox{ }}{\overset {..}{x}}\left(t\right),{\mbox{ }}{\overset {\mbox{.}}{x}}\left(t\right),{\mbox{ }}x\left(t\right){\mbox{ }}e{\mbox{ }}f\left(t\right)\in \mathbb {R} ^{n\times 1}}$ , denotam, respetivamente, os vetores aceleração, velocidade, deslocamento e excitação externa em função do tempo.

Para um sistema não amortecido e sem perturbação externa, a Eq. (1) é reduzida a

Uma solução do sistema linear homogéneo de segunda ordem (2) pode ser escrita como

onde ψ_i é um ângulo de mudança de fase, ϕ_i representa um vetor próprio e ω_i está relacionado com o valor próprio do problema de valores e vetores próprios como se explica a seguir.

Calculando a segunda derivada de (3) em ordem a t obtém‐se

que em conjunto com (3) pode ser inserido na expressão (2), obtendo‐se

Pré‐multiplicando a Eq. (5) por M^-1 pode‐se escrever a igualdade

ou introduzir a matriz identidade,

De modo a existir uma solução não‐trivial, ${\textstyle M^{-1}K-\omega _{i}^{2}I}$ deve ser singular. Assim,

A partir de det (M^-1K  − λ_iI ) = 0 pode‐se determinar os valores próprios ${\textstyle \lambda _{i}=\omega _{i}^{2}}$ em que cada valor próprio tem um vetor próprio ϕ_i  ≠ 0 associado, o qual pode então ser obtido a partir de solução de (7). Os vetores próprios ϕ_i correspondem às formas naturais de vibração, cujas componentes representam amplitudes, e ω_i representam as frequências naturais circulares medidas em rad/s, definidas como as raízes quadradas positivas dos valores próprios.

O problema de encontrar as frequências naturais e as formas naturais de vibração de um dado sistema dinâmico é aqui abordado, na medida em que a resposta modal dos modelos de placas rebitadas é usada para o processo de melhoramento de modelos de elementos finitos.

A partir da resposta do sistema também é possível, usando a função de resposta em frequência (FRF), determinar o fator de amortecimento ξ . O fator de amortecimentos é dado por [12]

onde Δω  = p₂  − p₁ é a largura de banda em frequência entre os pontos de média potência p₁ e p₂ , e ω é frequência natural circular do modo de vibração que está entre p₁ e p₂ .

Uma vez conhecido a amplitude máxima X_máx , os pontos de média potência são encontrados quando a amplitude X_p é determinada, em que

O fator de qualidade Q é dado por

3. Função objetivo

Nesta secção é descrito o formalismo desenvolvido para a função objetivo. A função objetivo f é dada pela soma de 3 funções específicas como descrito abaixo

onde a função f_φC representa e quantifica a diferença entre os pares de formas naturais de vibração correlacionadas do modelo de referência e do modelo numérico, a função f_φU representa e quantifica a diferença entre os pares de formas naturais de vibração não correlacionadas do modelo de referência e do modelo numérico, e a função f_λ representa e quantifica a diferença entre as frequências do modelo numérico e do modelo de referência.

A função objetivo f   depende do vetor ${\textstyle v}$ que contém os parâmetros a otimizar e quantifica a diferença do comportamento dinâmico entre o modelo numérico e o modelo de referência.

A função objetivo f_φC é dada por

onde N_C é o número de pares de formas naturais de vibração correlacionadas da diagonal da matriz MAC (Modal Assurance Criterion), ${\textstyle M{\overset {\mbox{ˆ}}{A}}C_{ii_{0}}}$ são os elementos da diagonal da matriz MAC inicialmente estimados, ${\textstyle {\overset {\mbox{ˆ}}{v}}_{0}}$ é o vetor que contém os valores iniciais dos parâmetros, ${\textstyle M{\overset {\mbox{ˆ}}{A}}C_{ii_{q}}}$ são os elementos da diagonal da matriz MAC estimados a cada iteração q   e ${\textstyle {\overset {\mbox{ˆ}}{v}}_{q}}$ é o vetor que contém os valores dos parâmetros estimados a cada iteração q , com q  ≥ 0. O valor do MAC entre 2 formas naturais de vibração ϕ_i and ${\textstyle \phi _{j_{q}}}$ é definido por [14]

onde ϕ_i é o vetor de referência que contém as formas naturais de vibração experimentais i   e ${\textstyle {\overset {\mbox{ˆ}}{\phi }}_{j_{q}}}$ é vetor que contém as formas naturais de vibração numéricas estimadas j a cada iteração q , com q  ≥ 0.

A função f_φU é dada por

onde N_U é o número de pares de formas naturais de vibração não correlacionadas fora da diagonal da matriz MAC.

A função f_λ é dada por

onde N_λ é o número de valores próprios λ correspondente aos pares de formas naturais de vibração correlacionadas (ver secção 1), f_i é a i^th frequência natural de referência obtida experimentalmente, ${\textstyle {\overset {\mbox{ˆ}}{f}}_{i_{0}}}$ frequência numérica inicial estimada definida por ${\textstyle {\overset {\mbox{ˆ}}{f}}_{i_{0}}={\frac {\sqrt {{\overset {\mbox{ˆ}}{\lambda }}_{i_{0}}}}{2\pi }}}$ , ${\textstyle {\overset {\mbox{ˆ}}{f}}_{i_{q}}}$ é a i^th frequência natural numérica estimada a cada iteração q   definida por ${\textstyle {\overset {\mbox{ˆ}}{f}}_{i_{q}}={\frac {\sqrt {{\overset {\mbox{ˆ}}{\lambda }}_{i_{q}}}}{2\pi }}}$ , ${\textstyle {\overset {\mbox{ˆ}}{\lambda }}_{i_{0}}}$ é o i^th valor próprio numérico inicial estimado, e ${\textstyle {\overset {\mbox{ˆ}}{\lambda }}_{i_{q}}}$ é o i^th valor próprio numérico estimado a cada iteração q , com q  ≥ 0. O termo quadrático na Eq. (16) é usado para acelerar a convergência e obter apenas diferenças positivas entre as frequências do modelo numérico e de referência. O denominador em (13), (15) e (16) é usado para obter diferenças normalizadas. Os vetores que contêm as formas naturais de vibração estimadas ${\textstyle {\overset {\mbox{ˆ}}{\phi }}_{i_{q}}}$ e os valores próprios estimados ${\textstyle {\overset {\mbox{ˆ}}{\lambda }}_{i_{q}}}$ são função dos parâmetros contidos no vetor ${\textstyle {\overset {\mbox{ˆ}}{v}}_{q}}$ . A relação entre eles pode ser escrita, respetivamente, como

e

onde p é o número de parâmetros a otimizar e n é o número de graus de liberdade.

Os parâmetros ${\textstyle {\overset {\mbox{ˆ}}{v}}_{q_{1}},...,{\overset {\mbox{ˆ}}{v}}_{q_{p}}}$ contidos no vetor ${\textstyle {\overset {\mbox{ˆ}}{v}}_{q}}$ estão sujeitos à restrição

onde a quantidade ${\textstyle v_{LB}}$ representa o limite inferior e a quantidade ${\textstyle v_{UB}}$ representa o limite superior.

O valor ótimo para a constante de rigidez do elemento mola‐amortecedor é obtido quando a função objetivo é minimizada. No entanto, o valor mínimo obtido na função objetivo é diferente para cada caso, uma vez que os modelos envolvidos em cada caso também são diferentes, não podendo ser considerado como uma referência direta para avaliar a melhor representação numérica da respetiva junta rebitada. O valor mínimo obtido na função objetivo apenas indica que o valor ótimo da constante de rigidez foi encontrado. A avaliação é então dada pela diferença percentual média entre as frequências naturais do modelo de referência e do modelo numérico, definida como

onde ${\textstyle {\overset {\mbox{ˆ}}{f}}_{i}}$ é a i^th frequência natural numérica estimada final.

4. Processo de melhoramento de modelos de elementos finitos

O algoritmo usado no processo de melhoramento de modelos de elementos finitos é apresentado nesta secção. O algoritmo é implementado em MATLAB^® , usando algumas ferramentas da sua Toolbox , e está preparado para interagir com o programa de cálculo de elementos finitos ANSYS^® .

4.1. Descrição do algoritmo

Considere‐se ${\textstyle {\overset {\mbox{ˆ}}{\varphi }}_{0}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ e ${\textstyle {\overset {\mbox{ˆ}}{\lambda }}_{0}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ , respetivamente como os vetores próprios numéricos iniciais estimados e os valores próprios numéricos iniciais estimados, a partir dos quais podem ser calculados ${\textstyle M{\overset {\mbox{ˆ}}{A}}C_{ii_{0}}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ , ${\textstyle M{\overset {\mbox{ˆ}}{A}}C_{ij_{0}}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ e ${\textstyle {\overset {\mbox{ˆ}}{\omega }}_{i_{0}}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ , onde ${\textstyle M{\overset {\mbox{ˆ}}{A}}C_{ii_{0}}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ são os elementos da diagonal da matriz MAC inicial estimada e ${\textstyle M{\overset {\mbox{ˆ}}{A}}C_{ij_{0}}\left({\overset {\mbox{ˆ}}{v}}_{0}\right)}$ são os valores fora da diagonal da matriz MAC inicial estimada.

A função objetivo inicial f⁰ , definida em (12), dependendo das Eqs. (13) ‐ (16), é usada para comparar os valores numéricos iniciais estimados com os valores de referência/experimentais. Se f⁰  → −1⁺ , onde ${\textstyle -1^{+}={\underset {\in \rightarrow 0}{lim}}\left(-1+\in \right)}$ ou se a condição ${\textstyle \left\|f^{q+1}-f^{q}\right\|\leq 10^{-6}}$ é satisfeita como critério de paragem, então o processo de melhoramento de modelos de elementos é concluído. Se nenhuma das condições acima é satisfeita, então q  = q  + 1 e novas estimativas ${\textstyle {\overset {\mbox{ˆ}}{\varphi }}_{q+1}\left({\overset {\mbox{ˆ}}{v}}_{q+1}\right)}$ e ${\textstyle {\overset {\mbox{ˆ}}{\lambda }}_{q+1}\left({\overset {\mbox{ˆ}}{v}}_{q+1}\right)}$ são calculadas. Novamente, estes valores são comparados com os valores de referência usando a função objetivo f^{q +1} baseada na Eq. (12). Se para alguma iteração q  = h  ≥ 1, tal que f^h  → −1⁺ , ou se a condição ${\textstyle \left\|f_{d}^{q+1}-f_{d}^{q}\right\|\leq 10^{-6}}$ é satisfeita, então o processo de melhoramento de modelos de elementos finitos é concluído e os valores ótimos de ${\textstyle {\overset {\mbox{ˆ}}{\varphi }}_{h}\left({\overset {\mbox{ˆ}}{\nu }}_{h}\right)}$ e ${\textstyle {\overset {\mbox{ˆ}}{\lambda }}_{h}\left({\overset {\mbox{ˆ}}{v}}_{h}\right)}$ são encontrados. O mesmo é dizer que o valor ótimo do parâmetro que carateriza a constante de rigidez do elemento mola‐amortecedor e que melhor aproxima o modelo numérico ao modelo de referência/experimental é encontrado.

4.2. Fluxograma do processo

O fluxograma que descreve a interação entre o MATLAB^® e o ANSYS^® é apresentado na figura 1 .

Figura 1. Fluxograma de interação entre o Matlab® e o Ansys® .

Os principais passos do processo de melhoramento de modelo de elementos finitos são os seguintes.

• Executa o ANSYS^® com o modelo de elementos finitos inicial a otimizar, construído em APDL num ficheiro input.txt   , com os parâmetros a otimizar contidos no vetor inicial ${\textstyle v_{0}}$ .
• Lê o ficheiro de resultados output.txt obtido do ANSYS^® e compara os valores estimados com os valores de referência/experimentais, de modo a construir a função objetivo f .
• Conclui o processo se um valor ótimo é atingido, ou seja, se as condições de paragem são satisfeitas, ou vai para o passo 4.
• Obtém os novos parâmetros ${\textstyle {\overset {\mbox{ˆ}}{v}}_{q+1}}$ definidos pelo otimizador, tendo em conta os constrangimentos definidos pela condição (19).
• Modifica o modelo de elementos finitos no ficheiro input.txt   com os novos parâmetros ${\textstyle {\overset {\mbox{ˆ}}{v}}_{q+1}}$ .
• Inicia uma nova análise com o novo ficheiro input.txt , indo para o passo 1.

O processo de melhoramento de elementos finitos usa no seu otimizador um solver de procura global baseado nos métodos estocásticos, o GlobalSearch , fornecido pela Toolbox do MATLAB^® . O GlobalSearch é processado em 2 fases: a fase local e a fase global. Na fase local, os pontos de amostragem, obtidos aleatoriamente, são manipulados por um solver local, o fmincon , de modo a encontrar vários candidatos a mínimo local. Na fase global, o mínimo local com menor valor da função objetivo é usado como uma aproximação para o ótimo global. Os restantes pontos de partida são rejeitados.

5. Procedimento experimental

As amostras para análise modal experimental consistem em 2 placas de alumínio (PL1 e PL2) interligadas, com 2 mm de espessura, como mostrado na figura 2 a. As placas estão ligadas por meio de rebites de alumínio com 3 mm de diâmetro. Quatro amostras, com junta composta por um, 2, 3 e 5 rebites, são submetidas a análise modal experimental. Os testes são realizados em temperatura ambiente, cerca de 20 °C, utilizando o equipamento LMS SCADAS^® para análise modal experimental. As amostras são testadas na condição próxima de livre no espaço, estando suspensas em 2 pontos por meio de um fio de nylon de 350 mm de comprimento, de modo a não causar interferência no ensaio, conforme mostrado na figura 2 b. Os testes são realizados utilizando um martelo de impacto instrumentado da PCB^® , sendo a força de impacto aplicada no ponto P1, e a resposta medida com um laser de interferometria Doppler em 8 pontos, P1 a P8 como mostrado na figura 2 c.

Figura 2. (a) Representação esquemática de um modelo com junta composta por 5 rebites; (b) modelo submetido a análise modal experimental; (c) localização dos pontos medidos.

Os 8 pontos selecionados são os mínimos necessários para a obtenção de dados de referência/experimentais suficientes de modo a que o processo de melhoramento de modelos de elementos fintos (updating ) seja efetuado com sucesso.

Foi efetuado outro ensaio experimental em que a resposta foi medida em 12 pontos, de modo a ser possível representar graficamente com boa aproximação as 8 primeiras formas naturais de vibração dos modelos experimentais. Os dados são coletados no domínio do tempo (amplitude em função do tempo) e processados no software de análise modal LMS^® , de modo a serem convertidos para o domínio da frequência. As frequências naturais e amplitudes de cada ponto são obtidas a partir da curva ajustada da FRF.

O material foi adquirido todo do mesmo lote e os cortes efetuados por meio de tecnologia laser. Foram efetuados 3 ensaios em condições idênticas de um mesmo caso para avaliar a repetibilidade dos resultados, verificando‐se uma diferença entre eles desprezável, pelo que se procedeu a posteriormente a um único ensaio para cada caso de estudo.

As 2 placas usadas para construir as amostras foram também elas sujeitas individualmente a análise modal experimental, com o fim de verificar se as placas modeladas em elementos finitos representam bem as placas físicas.

6. Modelos de elementos finitos

Os modelos numéricos são construídos utilizando o APDL, com as mesmas propriedades geométricas e mecânicas dos modelos experimentais, como apresentado na tabela 1 . As propriedades mecânicas dos modelos experimentais, como o módulo de Young e o coeficiente de Poisson, foram estimadas por meio de ensaios experimentais de tração.

As placas são modeladas com elementos de casca (shell63). O rebite é construído com elementos que combinam mola e amortecedor (combin14). A dimensão do elemento de casca é de 0,003 m de largura 0,003 m de altura. O amortecimento em materiais metálicos como o alumínio é considerado insignificante, pelo que pode ser desprezado, não sendo assim considerado no elemento mola‐amortecedor. Para justificar essa opção são apresentados na secção 7 os valores dos fatores de amortecimento ξ , assim como os respetivos valores do fator de qualidade Q para os modos de vibração obtidos experimentalmente no modelo de junta rebitada com 5 rebites.

Um valor inicial de 5 × 10⁷  N/m é atribuído à constante de rigidez da mola, que é variável durante o processo de melhoramento de modelos de elementos finitos. A massa do rebite obtida por pesagem é dividida pontualmente pelos nós que interligam as 2 placas por meio dos elementos mola‐amortecedor e é modelada com elementos de massa (mass21), não sendo utilizada como variável ao longo do processo de melhoramento.). As propriedades do rebite modelado em elementos finitos são apresentadas na tabela 2 .

A figura 3 mostra as várias representações do rebite em elementos finitos com diferentes quantidades de elementos mola‐amortecedor (4, 6, 8, 12, 16 e 20) por rebite.

Figura 3. Representação do rebite em elementos finitos com: (a) 4, (b) 6, (c) 8, (d) 12, (e) 16 e (f) 20 elementos mola‐amortecedor por rebite.

A figura 4 mostra em corte transversal os vários modelos em elementos finitos de 2 placas interligadas por diferentes quantidades de rebites (um, 2, 4 e 5).

Figura 4. Modelo da junta rebitada em elementos finitos com: (a) um, (b) 2, (c) 3 e (d) 5 rebites.

O processo de melhoramento de modelos de elementos finitos é usado de modo a encontrar o valor da constante de rigidez k dos elementos mola‐amortecedor e, consequentemente, da junta rebitada numérica que melhor reproduz o modelo experimental. É esperado que esse valor possa variar de acordo com a quantidade de rebites e elementos mola‐amortecedor por cada rebite. Para estudar esse efeito são construídos 24 modelos diferentes de juntas rebitadas em elementos finitos: 4 juntas rebitadas com diferentes quantidades de rebites (um, 2, 3 e 5), combinados com 6 diferentes quantidades de elementos mola‐amortecedor (4, 6, 8, 12, 16 e 20) por cada rebite.).

7. Resultados e discussão

Os resultados obtidos são apresentados nesta secção. Primeiro apresenta‐se nas Tabela 3  and Tabela 4 a comparação do comportamento dinâmico entre as dua2s placas físicas (PL1 e PL2, respetivamente) e o modelo em elementos finitos que as representa. A diferença percentual média, obtida pela Eq. (20) , entre as frequências naturais da placa PL1 e PL2 e o respetivo modelo de elementos finitos é de 0,37 e 0,35%, respetivamente. O valor 1 para o MAC revela que as formas naturais de vibração também estão bem correlacionadas. Podemos concluir que o modelo em elementos finitos representa com elevada aproximação as placas reais.

De modo a justificar a opção por desprezar o coeficiente de amortecimento nos modelos de elementos finitos, são apresentados na tabela 5 os valores dos fatores de amortecimento ξ e os respetivos fatores de qualidade Q para cada modo de vibração do modelo de junta rebitada com 5 rebites.

Os fatores de amortecimento ξ apresentados são muito pequenos (ξ  < 0.1) para todos os modos de vibração. Os fatores de qualidade Q são muito elevados e superiores a 0,5, o que demonstra que o sistema é do tipo subamortecido. Assim, o amortecimento pode ser desprezado nos modelos de elementos finitos.

A figura 5 mostra uma representação das 8 primeiras formas naturais de vibração para o modelo de elementos finitos com junta rebitada composta por 5 rebites. A figura 6 mostra a representação das 8 primeiras formas naturais de vibração obtidas experimentalmente em 12 pontos para o mesmo modelo.

Figura 5. Representação das 8 primeiras formas naturais de vibração do modelo de elementos finitos com junta rebitada composta por 5 rebites.

Figura 6. Representação das 8 primeiras formas naturais de vibração do modelo experimental com junta rebitada composta por 5 rebites.

Como se pode observar por comparação entre as Figura 5  and Figura 6 , existe boa correlação entre as formas naturais de vibração entre os 2 modelos. O mesmo acontece para os restantes modelos com um, 2 e 3 rebites.

O gráfico da figura 7 mostra a relação entre a diferença percentual média, obtida por aplicação da Eq. (20), entre as frequências naturais dos modelos numéricos melhorados e dos respetivos modelos experimentais, em função da quantidade de elementos mola‐amortecedor. As juntas construídas com 8 elementos mola‐amortecedor por rebite permitem obter a melhor representação em elementos finitos relativamente aos respetivos modelos experimentais com um e 2 rebites como se pode observar nas figura 7 a e b. A diferença percentual média usando 8 elementos mola‐amortecedor é menor relativamente às outras quantidades.

Figura 7. Diferença percentual, após updating , entre o modelo de elementos finitos e o modelo experimental de junta rebitada com (a) um rebite, (b) 2 rebites, (c) 3 rebites e (d) 5 rebites.

O mesmo acontece para os casos de juntas rebitadas com 3 e 5 rebites. No entanto, as modelações com 12 elementos mola‐amortecedor para esses mesmos casos também apresentam a mesma diferença percentual média, como pode ser observado nas figura 7 c e d. No entanto, a sua maior complexidade de modelação relativamente ao rebite com 8 elementos mola‐amortecedor não justifica a sua implementação. Os gráficos mostram também que o aumento de rebites na junta reduz a diferença percentual média entre o modelo de elementos finitos e o modelo de referência.

As Tabela 6 , Tabela 7 , Tabela 8  and Tabela 9 mostram a evolução das frequências naturais dos modelos de elementos finitos e os valores MAC, antes (coluna 7) e após (coluna 8), o updating , para os modelos com um, 2, 3 e 5 rebites, respetivamente. Nessas tabelas também podem ser observadas as diferenças percentuais entre as 8 primeiras frequências naturais do modelo de referência (coluna 2) e do modelo de elementos finitos, antes do updating (coluna 4) e depois do updating (coluna 6). Apenas são mostrados os resultados para os modelos com 8 elementos mola‐amortecedor por cada rebite.

A diferença percentual média, obtida pela aplicação da Eq. (20) , é de 3,92% para o modelo com um rebite (figura 7 a), 0,36% para o modelo com 2 rebites (figura 7 b), 0,28% para os modelos com 3 rebites (figura 7 c) e 5 rebites (figura 7 d). Confirma‐se que essa diferença se reduz com o aumento do número de rebites na junta. Por outro lado, o valor de um para o MAC inicial e final revela que as formas naturais de vibração estão bem correlacionadas para todos os casos. Este fato indica que os valores ótimos obtidos para as constantes de rigidez k dos elementos mola‐amortecedor, apresentados na tabela 10 , são obtidos com elevada robustez.

8. Conclusões

Neste artigo, vários modelos de elementos finitos simplificados de juntas rebitadas são estudados e avaliados. A avaliação é efetuada com base na comparação entre as características modais desses modelos, após aplicado o processo de updating , e os respetivos modelos experimentais. Com base nos resultados obtidos conclui‐se que as juntas rebitadas com 8 elementos mola‐amortecedor por rebite apresentam a diferença média mais baixa entre o modelo numérico e o modelo experimental. Para as juntas com 3 e 5 rebites, o modelo com 12 elementos mola‐amortecedor por rebite apresenta o mesmo comportamento que o modelo com 8 elementos mola‐amortecedor por rebite. No entanto, a maior complexidade deste modelo não justifica a sua escolha.

Os resultados das avaliações permitem concluir que os modelos de elementos finitos simplificados propostos permitem representar os modelos de juntas rebitadas experimentais com elevada robustez. O modelo proposto, com 8 elementos mola‐amortecedor por rebite, permite também boa eficiência computacional devido à simplicidade do tipo de elementos usados.

Com este estudo pretendemos contribuir com informação útil para os designers na modelação de juntas rebitadas para análise em elementos finitos de dinâmica de estruturas.

Agradecimentos

Os autores agradecem ao Centre for Mechanical and Materials Technologies (CT2 M) da Universidade do Minho.

Referências

1. [1] D.V.T.G. Pavan Kumar, S. Sathiya Naarayan, S. Kalyana Sundaram, S. Chandra; Further numerical and experimental failure studies on single and multi‐row riveted lap joints; Eng. Fail. Anal., 20 (2012), pp. 9–24
2. [2] S. Sathiya Naarayan, D.V.T.G. Pavan Kumar, S. Chandra; Implication of unequal rivet load distribution in the failures and damage tolerant design of metal and composite civil aircraft riveted lap joints; Eng. Fail. Anal., 16 (2009), pp. 2255–2273
3. [3] W. Huang, T-J. Huang, Y. Garbatov, C. Guedes Soares; Fatigue reliability assessment of riveted lap joint of aircraft structures; Int. J. Fatigue, 43 (2012), pp. 54–61
4. [4] M.R. Urban; Analysis of the fatigue life of riveted sheet metal helicopter airframe joints; Int. J. Fatigue, 25 (2003), pp. 1013–1026
5. [5] S. Xu, X. Deng; An evaluation of simplified finite element models for spot‐welded joints; Finite Elem. Anal. Des., 40 (2004), pp. 1175–1194
6. [6] M.I. Friswell, J.E. Mottershead; Finite element model updating in structural dynamics; Kluwer Academic Publishers, Dordrecht (2000)
7. [7] J.D. Collins, C.G. Hart, T.K. Hasselman, B. Kennedy; Statistical identification of structures; AIAA Journal, 12 (1974), pp. 185–190
8. [8] C.P. Fritzen; Identification of mass, damping and stiffness matrices of mechanical systems; J. Vib. Acoust. Stress. Reliab. Des., 108 (1986), pp. 9–17
9. [9] E.M. Hernandez, D. Bernal; Iterative finite element model updating in the time domain; Mech. Syst. Signal Process., 34 (1‐2) (2013), pp. 39–46
10. [10] J. Meireles; Análise dinâmica de estruturas por modelos de elementos finitos identificados experimentalmente; Tese de Doutoramento em Engenharia Mecânica, Universidade do Minho, Guimarães (2008)
11. [11] M. Dourado, J. Meireles, A.M.A.C. Rocha; A global optimization approach applied to structural dynamic updating; Lect. Notes Comput. Sc., 8580 (2014), pp. 195–210
12. [12] C. F. Beards, Structural Vibration: Analysis and Damping, Arnold, London,;1; 1996.
13. [13] N.M.M. Maia, J.M.M. Silva; Theoretical and experimental modal analysis; Research Studies Press Ltd, Hertfordshire (1997)
14. [14] R. J. Allemang, L. D. Brown, A Correlation Coefficient for Modal Vector Analysis, in: Proc. 1 st Int. Modal Analysis Conference & Exhibit, Florida, Holiday Inn (1982), p. 110‐116.