Estudo de Molas e Barras no MEF
Estudo de Molas e Barras no MEF
1. Introdução
1. Introdução
Etapa 5: Montar as equações dos elementos para obter as equações globais ou totais e introduzir
condições de contorno
1. Introduction
Ver capítulo 1
É necessário escolher antecipadamente uma função matemática para representar a forma deformada do elemento finito sob carga.
Porque é difícil obter uma forma fechada ou uma solução exata, supomos uma forma de solução ou uma
distribuição de deslocamento dentro do elemento usando uma função matemática apropriada. As funções os
as polinômios são os mais comumente usados.
Em geral, o número total de coeficientes da função polinomial é igual ao número total de graus de liberdade associados a
o elemento.
1. Introdução
! = !+ " !
Exemplo 2: função de deslocamento de um elemento viga
Dans le cas d’une poutre, le nombre de degrés de liberté locaux est 4. Donc, la fonction de déplacement # derrubar
representar o deslocamento transversal através do comprimento do elemento é:
1. Introdução
Para pequenas deformações, as tensões devem estar ligadas às deformações pela lei de Hooke dada por
Geralmente, o desenvolvimento das matrizes de rigidez dos elementos e das equações dos elementos é baseado em 4 métodos:
•Método de coeficientes de influência da rigidez
Método de equilíbrio direto
Método energético
Métodos de resíduos ponderados
1. Introdução
1. Introdução
Étape 6 : Résoudre les degrés de liberté inconnus (ou déplacements généralisés) Ver capítulos 1 e 2
Etapa 7: Resolver as restrições e limitações do elemento
Veja a etapa 4, capítulos 1 e 2
Étape 8 : Interpréter les résultats
Montagem (Etapa 5)
Referências locais
Referências globais
( =( *
F=
As letras minúsculas como x, y e z As letras minúsculas como x, y e z sem o
symbole ^ désignent des variables de coordonnées O símbolo ^ designa variáveis de coordenadas
locais. globais.
Convenção de sinal ^y
+1 +2 . 2
y
. 1
x^
i=1 i=2
elemento barra
As direções positivas para as forças e os deslocamentos em cada nó são consideradas na direção positiva # como
indicado do nó 1 ao nó 2 na figura. O símbolo k (não é é chamado +de constante da mola ou rigidez do
ressort.
• +são
$ as forças nodais locais
. são os deslocamentos nodais locais
• $"
•k = AE/L, onde A representa a seção transversal da barra, E é o módulo de elasticidade e L é o comprimento da barra.
Agora queremos desenvolver uma relação entre as forças nodais e os deslocamentos nodais para um elemento
ressort. Esta relação será a matriz de rigidez. Portanto, queremos relacionar a matriz de força nodal à matriz
de deslocamento nodal conforme segue:
! = !+ " !
Então
Com e
Antes da deformação
As forças de tração T produzem um alongamento total (deformação)
do ressort.
Após deformação
Ou T T
Então
Nós temos:
Então, 11= = /
. !! = − = − (0 − 1)= k
21= − =− /
.#! = − .!! = = 0 − 1= −
.## = 22
L’équilibre mécanique :
.!# = 12
Estamos procurando k22, então :
.## = − .!#
. !" = 0 =0
22= - 12
. #" = 1
Nós temos:
Então, . ## = = (1 - 0)= k 22= = /
Finalmente:
Aqui, os $" são escritos acima das colunas e ao lado das linhas nas matrizes k indicam os graus de liberdade associados
a cada linha e coluna de elemento.
Yassine EL HALOUI Capítulo 3: Estudo de molas, barras e treliças 19
Curso: Método dos Elementos Finitos
Então:
Além disso, os elementos 1 e 2 devem permanecer conectados ao nó comum 3 durante todo o deslocamento. É isso que se chama
a exigência de continuidade ou de compatibilidade. A exigência de compatibilidade dá:
(!) (#)
Sabendo que !" = !" e =#" #" Então:
Este método confiável de montagem direta de matrizes de rigidez de elementos individuais para formar a matriz de rigidez de
A estrutura total e o conjunto total de equações de rigidez é chamada de método de rigidez direta. Esta é a etapa mais importante.
da metodologia dos elementos finitos.
Yassine EL HALOUI Chapitre 3 : Etude des ressorts, barres et treillis 22
Curso: Método dos Elementos Finitos
Para este exemplo simples, é fácil aumentar as matrizes de rigidez dos elementos e depois sobrepô-las para chegar a
matriz de rigidez total. No entanto, para problemas envolvendo um grande número de graus de liberdade, ele se tornará
dificilmente expandir cada matriz de rigidez de elemento à ordem da matriz de rigidez total.
Para evitar essa expansão da matriz de rigidez de cada elemento, sugerimos uma forma direta ou abreviada da
méthode de rigidité directe pour obtenir la matrice de rigidité totale. Pour l'exemple d'assemblage du ressort, les lignes et
as colunas de cada matriz de rigidez dos elementos são etiquetadas de acordo com os graus de liberdade que lhes são atribuídos
associados como segue:
A matriz de rigidez global K é então construída simplesmente adicionando diretamente os termos associados aos graus
de liberdade em (!) et (#) nos seus locais de grau de liberdade idênticos correspondentes em K da seguinte forma:
Condições de contorno
Devemos especificar as condições de contorno (ou de suporte) para este exemplo, ou K será singular; ou seja, que o
o determinante de K será nulo e seu inverso não existirá. Isso significa que o sistema estrutural é instável. Sem o nosso
especificação de restrições cinemáticas ou de condições de suporte adequadas, a estrutura estará livre para se mover como
um corpo rígido e não resistirá a nenhuma carga aplicada.
Resolução (passo 6)
Exercício 1
Para a montagem de três barras ilustrada na figura seguinte, determine (a) a matriz de rigidez global, (b) os
deslocamentos dos nós 2 e 3, e (c) as reações nos nós 1 e 4.
Uma força de 13340 N é aplicada na direção x no nó 2. O comprimento de cada elemento é
de 80 cm. Seja E1=e2=200 GPa e A1=A2= 25 cm2 para os elementos 1 e 2, e E3= 100 GPa e A3= 50 cm2 para
o elemento 3.
Na seção 2.6, estabelecemos as equações do elemento, bem como a matriz de rigidez global no caso em que os
Os referenciais local e global estão confundidos. Sabendo que esses referenciais geralmente não estão confundidos, uma transformação
deve ser feita a transição das equações locais (vetores) para as equações globais.
Com: e
Então:
É a matriz de transformação ou de rotação. Ela será utilizada a seguir para desenvolver a matriz.
de rigidez global para um elemento de barra com orientação arbitrária e para transformar os
deslocamentos nodais globais e as forças globais em relação às que são locais.
Agora vamos usar a relação de transformação para obter a matriz de rigidez global de um elemento barra.
Precisamos da matriz de rigidez global de cada elemento para montar a matriz de rigidez global total de
a estrutura.
Mostramos que para um elemento de barra no sistema de coordenadas local, a equação básica é:
Observamos a partir da equação anterior que um total de quatro componentes de força e quatro de deslocamento
ocorrem quando as coordenadas globais são usadas. No entanto, um total de duas componentes de força e duas de
Os deslocamentos aparecem na representação em coordenadas locais de uma mola ou de uma barra.
Usando as relações entre os componentes de força locais e globais e entre os componentes de deslocamento
locais e globais, poderemos obter a matriz de rigidez global. Sabemos da relação de transformação
seguinte :
Com Então:
No entanto, para escrever a relação final dando f(globais) em função de d(globais), precisaremos inverter a
matriz T. Isso não é possível porque a matriz T não é uma matriz quadrada. Portanto, é necessário expandir a matriz T.
Com:
k localize o elemento
Então:
k global do elemento
Mostramos que T-1=Tt
Agora, sabendo que a função de deslocamento foi considerada contínua, a matriz de rigidez global para cada
o elemento pode ser adicionado usando o método de rigidez direta para obter
K(en majuscule = matriz de rigidez da estrutura) agora relaciona as forças nodais globais F aos deslocamentos
nodaux globaux d pour l'ensemble de la structure par :
Exercice 2
Agora vamos considerar a determinação da contração em um elemento de barra. Para uma barra, as forças
As localidades estão ligadas aos deslocamentos locais pela seguinte equação:
A definição geral da tensão de tração axial é a força axial dividida pela seção transversal.
consequentemente, a constrição axial é
Onde, . #"é utilizado (em vez de ). !"porque é a força axial que puxa em
barra ( . #"est la force de traction), comme illustré à la figure ci-contre. Donc,
Com,
Exercício 3
Resposta:
5. TD
y
Exercício 1
A treliça plana representada na figura ao lado é composta por três
pinturas de mesma natureza e de mesma seção reta.
Seja E o módulo de Young do material e A a área das seções transversais.
O nó 1 é articulado e o nó 3 repousa em um apoio simples cujo L
a normal é horizontal. 1
O nó 2 carrega uma carga de P.
x
L L
1) Discretizar e selecionar o tipo do elemento finito.
2) Déterminer les matrices de rigidité locales de chaque élément choisi.
3) Déterminer les matrices de rigidité globales de chaque élément choisi.
4) Montar as matrizes de rigidez globais para obter a matriz de rigidez global total da estrutura.
5) Dar a relação F=Kd mostrando as condições de contorno da estrutura.
6) Calculer les déplacements aux nœuds 1, 2 et 3.
7) Calculer les forces de réactions.
8) Déterminer les contraintes dans chaque barre. Déduire la nature de ces contraintes (traction ou compression).
9) Verifique os resultados do exercício examinando o equilíbrio mecânico das forças.
5. TD
Exercício 2
A malha plana representada na figura ao lado é composta por três
pilares da mesma natureza e da mesma seção reta.
Soient E le module deYoung du matériau et A l’aire des sections droites.
O nó 1 é articulado e o nó 2 repousa sobre um apoio simples, cujo
a norma é horizontal.
O nó 3 carrega uma intensidade (P,3P,0).
5. TD
Exercício 3
A treliça plana representada na figura ao lado é composta por três
poutres de mesma natureza e de mesma seção reta.
5. TD
Exercice 3
5. TD
Exercice 3 0+ 0,5/racine
(2)+1=1,354
d1x
d1a
d2x
d2y d1x
d1 ano
d1x
d4y
d1x
5. TD
Exercice 3
5. TD
Exercício 3