02/11/2017
Métodos Numéricos para Engenharia
e
Método dos Elementos Finitos
Prof° MSc Edmilton Romão da Silva
Determinando os valores para a montagem da
Matriz de Rigidez Global
Prof. Edmilton Romão da Silva
Identificando as equações Algébricas lineares
das Forças de Reação e Forças Aplicadas
Prof. Edmilton Romão da Silva
1
02/11/2017
Determinando os deslocamentos
Resolvendo o sistema pelo método de Gauss (escalonamento)
Definindo as condições de contorno de deslocamentos e restrições
6300 *1,36 = 8568
6300 * 0,48 = 3024
Prof. Edmilton Romão da Silva
Generalização do Procedimento para Qualquer
Tipo de Estrutura
Como {FII}, {∆ I}, {KI,I} e {KII, II} são conhecidos, resolve-se primeiro a
equação (2), determina ∆ II, substitui ∆ II na equação (1) e encontra FI
Prof. Edmilton Romão da Silva
Determinando os deslocamentos da Estrutura
inteira pelo Método de Gauss
Método de Gauss ou Método das Eliminações Sucessivas:
Elimina a 1ª incógnita das equações, com exceção da 1ª equação
Elimina a 2ª incógnita das equações, com exceção da 2ª equação
Elimina a 3ª incógnita das equações, com exceção da 3ª equação
Processo Forward Reduction
Quando ficar só uma incógnita, terá se encontrado o valor desta, aí então, basta
voltar substituindo as incógnitas no processo inverso (Backsubstitution)
Prof. Edmilton Romão da Silva
2
02/11/2017
Determinando os deslocamentos da Estrutura
inteira pelo Método de Gauss
Multiplicado pela primeira linha e adicionado na segunda
Multiplicado pela primeira linha e adicionado na terceira
Prof. Edmilton Romão da Silva
Determinando os deslocamentos da Estrutura
inteira
Multiplicado pela segunda linha e adicionado na terceira
Multiplicado pela terceira linha e adicionado na quarta
Prof. Edmilton Romão da Silva
Determinando os deslocamentos e das Forças
de Reações da Estrutura
Sabendo que ∆ 1 = ∆ 2= ∆ 3 = ∆ 4 =0, pode se calcular as Forças
de reações da Estrutura, substituindo os valores de deslocamentos
encontrados nas respectivas equações.
F1 = −615,41 kgf ( sentido oposto ao eixo X )
F2 = −1333,36 kgf ( sentido oposto ao eixo Y)
F3 = −384,62 kgf ( sentido oposto ao eixo X )
F4 = 833,31 kgf (mesmo sentido do eixo Y)
Prof. Edmilton Romão da Silva
3
02/11/2017
Sistema de Equilíbrio da Estrutura
Prof. Edmilton Romão da Silva
Determinando as Forças Internas de cada
Elemento
AE Força interna em um elemento de barra qualquer – sistema local
Fe = (u 2 − u1 )
L
Para calcular a força interna na Barra, os deslocamentos U1, V1, U2 e V2
globais, devem ser transformados para o sistema local
AE
Fe = (u 2 − u1 )
L
Prof. Edmilton Romão da Silva
Determinando as Forças Internas de cada
Elemento
Fazendo
Portanto:
Prof. Edmilton Romão da Silva
4
02/11/2017
Determinando as Forças Internas de cada
Elemento
Prof. Edmilton Romão da Silva
Determinando as Forças Internas de cada
Elemento
Prof. Edmilton Romão da Silva
Determinando as Forças Internas de cada
Elemento
Pegando como exemplo o elemento b, tem-se:
Lembrando que :
A = 300mm 2
L = 1000mm
E = 21000kgf / mm 2
O sinal negativo indica que existe uma tensão de compressão nesse elemento
Prof. Edmilton Romão da Silva
5
02/11/2017
Determinando as Forças Internas de cada
Elemento
Prof. Edmilton Romão da Silva
Aspectos Gerais na Formulação de Qualquer
Elemento Finito
Prof. Edmilton Romão da Silva