TEORIAS DAS
ESTRUTURAS
Diego Guimarães
Revisão técnica:
André Luís Abitante
Mestre em Engenharia e Ciência dos Materiais,
ênfase em Controle de Processos
Engenheiro Civil
Rossana Piccoli
Mestre em Engenharia Civil
Engenheira Civil
T314 Teoria das estruturas / Douglas Andrini Edmundo
... [et al.] ; [revisão técnica: André Luís Abitante, Rossana
Piccoli]. – Porto Alegre : SAGAH, 2018.
392 p. : il. ; 22,5 cm.
ISBN 978-85-9502-354-3
1. Engenharia civil. I. Edmundo, Douglas Andrini.
CDU 624.01
Catalogação na publicação: Karin Lorien Menoncin CRB-10/2147
Conceitos básicos
de métodos dos
elementos finitos
Objetivos de aprendizagem
Ao final deste texto, você deve apresentar os seguintes aprendizados:
Identificar os métodos numéricos aplicados para solução de problemas
matriciais.
Determinar a matriz de rigidez.
Determinar a solução do sistema linear.
Introdução
Neste capítulo, você vai ver como analisar e conhecer os princípios bá-
sicos de métodos numéricos, o conceito de iteração e os critérios de
convergência utilizados como base para a formulação do método dos
elementos finitos.
Os problemas de engenharia de estruturas partem da premissa da
escolha de um modelo físico e matemático da estrutura analisada para
se obter a solução. Entretanto, para entender melhor o comportamento
das estruturas, você precisa avaliar diversas alternativas e pontos de
vista. Para isso, você deve utilizar modelos matemáticos e físicos mais
complexos, que por sua vez exigem a aplicação de métodos numéricos.
Os conhecimentos de métodos numéricos são importantes para o de-
senvolvimento do método dos elementos finitos, das não linearidades,
das análises dinâmicas de estruturas e da transferência de calor.
192 Conceitos básicos de métodos dos elementos finitos
Métodos dos elementos finitos
Os métodos dos elementos finitos são um processo de análise de estruturas com
base em procedimentos numéricos aplicado a estruturas e meios contínuos.
Esse processo tem como fundamento o conceito de discretização do problema.
A discretização de um problema consiste em transformar um problema
complexo na soma de diversos problemas simples buscando soluções locais.
Dessa maneira, as propriedades podem garantir sua convergência para os
problemas globais.
O primeiro método numérico com o objetivo de resolver equações diferen-
ciais parciais foi o método das diferenças finitas. Nele, o domínio da solução
é dividido em uma malha de pontos ou nós discretos, aplicado para cada nó
da malha, e as derivadas são substituídas por diferenças finitas divididas.
Conceitualmente, essa aproximação é de fácil compreensão, porém existem
alguns inconvenientes na aplicação desse método. A aplicação do método
das diferenças finitas torna-se difícil em sistemas estruturais que possuem
geometria irregular, composições heterogêneas ou condições de contorno
não usuais.
Já o método dos elementos finitos apresenta uma alternativa melhor para
aplicação em sistemas que possuam tais características. Nesse método, o
domínio da solução é dividido em formas simples de regiões ou de elementos.
Isso permite que uma solução aproximada pelo método das diferenças finitas
seja desenvolvida para cada um dos elementos, convergindo para a solução total
do sistema, por meio do ajuntamento das soluções parciais de cada elemento,
constituindo a solução final.
A utilização de elementos no método dos elementos finitos — em vez de
malhas retangulares no método das diferenças finitas — proporciona uma
aproximação melhor em sistemas estruturais com formatos irregulares. Assim,
permite inclusive a geração contínua de valores desconhecidos por meio do
domínio de uma solução inteira, em vez resultados parciais por pontos isolados.
Conceitos básicos de métodos dos elementos finitos 193
O método dos elementos finitos é mais recente que o método das diferenças finitas.
Além disso, ele abrange vários problemas que não podem ser resolvidos por diferenças
finitas, podendo ser aplicado na solução de estruturas com geometria complexa e
em ambientes onde existam mudanças de meio.
O método dos elementos finitos pode ser aplicado em diversas áreas, como
em análise de estruturas, análise de escoamentos de fluidos, distribuição de
temperaturas, eletromagnetismo e projetos de equipamentos eletromecânicos
como transformadores e contatores. Em alguns casos, ele é a única ferramenta
com capacidade de fornecer uma solução aceitável, ainda que do ponto de vista
da matemática as soluções sejam consideradas aproximadas.
Métodos numéricos
Os métodos numéricos comumente aplicados na programação de algoritmos
utilizados para a análise estrutural no método dos elementos finitos são em-
pregados para a resolução dos problemas a seguir:
raízes de equações;
sistemas de equações;
interpolação;
integração;
equações diferenciais.
Erros em cálculo numérico
Os erros em cálculo numérico ocorrem por diversos motivos e são inerentes ao
processo. A simplificação dos modelos matemáticos e modelos físicos provoca
aproximações na análise da estrutura e impõe incertezas e erros que refletem
no processo de cálculo e análise da estrutura. Entre os erros mais comuns, você
pode considerar os erros na determinação dos parâmetros básicos da análise
da estrutura, erros de truncamentos por meio de representação limitada no
194 Conceitos básicos de métodos dos elementos finitos
cálculo de séries finitas. A representação fixa de algarismos significativos nos
cálculos também pode aumentar devido à sequência de cálculos.
Quanto maior a precisão dos dados obtidos na análise das condições de
contorno da estrutura, ou seja, quanto mais preciso for o modelo matemático,
maior será a precisão dos resultados. Ainda que o método dos elementos
finitos realize uma análise muito mais detalhada que a solução analítica de
uma estrutura, ele apresenta erros inerentes ao processo. Esses erros não
impedem a análise ou o dimensionamento da estrutura; eles apenas sina-
lizam que o processo pode ser mais preciso caso sejam tomadas medidas
para dirimir os erros nas diversas etapas. Os erros no desenvolvimento do
modelo físico e do modelo matemático podem ser reduzidos com uma análise
mais criteriosa do problema real. Isso será possível uma vez que o projetista
conheça profundamente o comportamento da estrutura e saiba analisar suas
condições de contorno.
Os erros de precisão matemática muitas vezes estão vinculados à limitação
das máquinas e dos softwares de análise. Portanto, não resultam de uma mera
falta de experiência do projetista.
Cada método numérico, com base na forma com que foi desenvolvido,
apresenta um erro maior ou menor, dependendo das iterações que cada método
possuir. Alguns métodos são mais precisos que outros. Por exemplo, o método
de Euler é menos preciso que o método de Range Kutta de quarta ordem. Ou
seja, a escolha do método numérico também vai influenciar o resultado da
análise de uma estrutura.
Como em toda análise estrutural, aqui o ponto de partida sempre será a análise das
condições de contorno que envolvem as restrições e os graus de liberdade. A precisão
com que será feita essa análise influenciará diretamente os resultados finais da análise
do comportamento da estrutura.
A escolha do método numérico, cujo grau de precisão deve ser compatível com
o nível de resultados desejados para cada problema, fará a diferença entre maior
ou menor erro de convergência em relação aos resultados exatos de um problema.
Conceitos básicos de métodos dos elementos finitos 195
Matriz de rigidez
O método dos elementos finitos baseia-se no método da rigidez ou dos deslo-
camentos para compor sua análise matricial.
Considere uma barra de comprimento L submetida a uma força F, com área
da seção transversal A constante, cujo material possui comportamento linear
e elástico e módulo de elasticidade constante E. Pode-se determinar qual é a
intensidade da força necessária para provocar um deslocamento u no sentido
da força F. Você pode observar que esses conceitos também se aplicam aos
métodos de energia.
F1 F2
u1 u2
nó1 nó2
Figura 1. Elemento finito de uma barra.
Se você considerar que o deslocamento no nó 2 é maior que zero (u2 > 0)
e que o deslocamento no nó 1 é igual a zero (u1 = 0), para que se mantenha o
equilíbrio do sistema estrutural, será necessário que:
–EA u2
F1 =
L
EA
F2 = u2
L
196 Conceitos básicos de métodos dos elementos finitos
Admitindo F1 e F2 como diferentes de zero, você pode escrever a forma
matricial do elemento finito da barra:
[ ]
EA 1 –1 u1 F1
=
L –1 1 u2 F2
Portanto,
[Ke] |ue| = |Fe|
Onde:
[Ke] é a matriz de rigidez do elemento;
|ue| é o vetor dos deslocamentos nodais do elemento;
|Fe| é o vetor nas forças nodais do elemento;
A matriz de rigidez do elemento é simétrica, portanto representa uma
relação linear entre os esforços e os deslocamentos. Essa relação permite
que a solução obtida seja exata quando aplicada a barras de seção constante
submetidas a cargas pontuais.
Por outro lado, se a geometria for variável, os resultados serão aproximados
devido aos erros de aproximação da modelagem e da geometria.
Figura 2. Com a aplicação do método dos elementos finitos, você
pode determinar com maior precisão as distribuições de tensões
em todas as partes de uma estrutura por meio da modelagem
computacional de cada elemento que compõe a estrutura.
Fonte: Dvoevnore/[Link].
Conceitos básicos de métodos dos elementos finitos 197
No exemplo a seguir, você vai analisar um elemento estrutural submetido a duas
forças, F1 e F2.
Primeiro, você deve determinar os deslocamentos nos nós 1 e 2 da Figura 3.
F1
F2
01 02
1.50
Figura 3. Elemento estrutural submetido a duas forças.
Considere o valor da força F1 de 10 kN e da força F2 de 15 kN. O módulo de elasticidade
E é de 250 GPa e a área da seção transversal A é de 1,27 cm2.
250 GPa = 25.000 kN/cm2
DCL
F1 F2
u1 u2
01 02
1.50
Figura 4. Deslocamentos nos nós 01 e 02.
198 Conceitos básicos de métodos dos elementos finitos
Formulação matricial para elementos finitos:
[Ke] |ue| = |Fe|
Matriz de rigidez do elemento:
EA 1
L –1 [ ]
–1
1
25.000 ∙ 1,27 1
150 –1 [ –1
1 ]
211,67
[1
–1
–1
1 ]
K=
[211,67
–211,67
–211,67
211,67 ]
Vetor de deslocamentos nodais:
u1
u=
u2
Vetor das forças nodais do elemento:
10
F=
15
Aplicando a formulação matricial, você tem:
[Ke] |ue| = |Fe|
[ 211,67
–211,67
–211,67
211,67 ] u1
u2
=
10
15
O sistema linear de equações será:
{ 211,67u1 – 211,67u2 = 10
– 211,67u1 + 211,67u2 = 15
Resolvendo o sistema de equações, você encontra os valores dos deslocamentos
nodais:
u1 = 5,905 . 10 -3 cm
u2 = 5,905 . 10 -2 cm
Conceitos básicos de métodos dos elementos finitos 199
Solução de sistemas lineares
O método dos elementos finitos é uma simulação estrutural que utiliza um
modelo matemático desenvolvido a partir do estado inicial de um elemento
estrutural antes da primeira iteração, que adota métodos numéricos para a
convergência dos resultados.
O objetivo do método é transformar uma “superfície” em um sistema linear
de equações que possa ser resolvido por meio dos modelos matemáticos.
O método dos elementos finitos não possui limite de aplicação. Assim, pode ser utilizado
não somente na engenharia de estruturas, mas na área da saúde, com a análise de
próteses de membros e próteses dentárias. Ele pode ser aplicado no estudo de todo
um sistema estrutural ou apenas em parte desse sistema, a qual se deseja analisar com
maior grau de precisão e detalhe.
A análise geralmente é feita com uso de ferramentas computacionais que permitem
maior velocidade de análise com maiores recursos para estudos.
Todo sistema linear de equações nada mais é que um conjunto de equa-
ções lineares. Um sistema de equações é definido a partir do conjunto de m
equações lineares originadas da formulação matricial, com n incógnitas, que
tenha a forma:
{
a11x1 + a12x2 ... + a1nxn = b1
a21x1 + a22x2 ... + a2nxn = b2
am1x1 + am2x2 ... + amnxn = b
Esse sistema linear de equações pode ser representado pelo produto de
matrizes AX = B, onde:
a11 a12 ... a1n x1 b1
a21 a22 a2n x2 b2
A= ,X= eB=
...
...
...
...
...
...
am1 am2 ... amn xn bn
200 Conceitos básicos de métodos dos elementos finitos
A solução do sistema linear de equações será uma matriz coluna de nú-
meros reais:
S1
S2
S=
...
Sn
Isso significa que todas as equações do sistema são satisfeitas quando as
incógnitas são substituídas pelos valores da matriz solução S, em que
x1 = s1
x2 = s2
xn = sn
De tal forma que AS = B.
Confira, no link a seguir, um estudo sobre a análise não linear física da flexão de vigas
de concreto armado utilizando uma formulação variacional com base no método
dos elementos finitos:
Nesse estudo, o método dos elementos finitos é utilizado para prever a reação à
flexão de vigas de concreto armado por meio de um modelo de proposta que combina
a teoria clássica de laminação com um modelo de dano aplicado ao concreto.
[Link]
1. Analise as afirmativas a seguir, não numéricos, sendo
referentes a definição de observado apenas o
métodos dos elementos finitos, comportamento estrutural.
e assinale a alternativa correta. b) Além da análise estrutural,
a) Os métodos dos elementos os métodos dos elementos
finitos permitem a análise de finitos podem ser utilizados
estruturas por procedimentos
Conceitos básicos de métodos dos elementos finitos 201
na análise do comportamento individualmente pelo método
de meios descontínuos. dos elementos finitos.
c) Os métodos dos elementos 3. Analise as afirmativas a seguir,
finitos consistem em referentes aos erros em
procedimentos numéricos cálculo numérico, e assinale
baseados no conceito de a alternativa correta.
discretização dos problemas. a) A simplificação dos modelos
d) A discretização consiste na físicos e matemáticos visa à
transformação de vários eliminação de erros e incertezas
problemas simples em um único na análise estrutural.
problema complexo, de modo b) Uma das vantagens do método
a facilitar a análise estrutural. dos elementos finitos é a
e) Os métodos dos elementos eliminação de erros existentes
finitos permitem a obtenção nas soluções analíticas clássicas.
de soluções globais, c) Os erros em cálculo numérico
garantindo sua convergência podem ocorrer por diversos
para os problemas locais. motivos, sendo estes
2. Analise as afirmativas abaixo, exteriores ao processo.
referentes aos métodos de d) Os erros podem ser reduzidos
diferenças finitas e dos elementos em função da experiência
finitos, e assinale a alternativa correta. dos projetistas, contudo, não
a) A aplicação do método dos dependem unicamente desta.
elementos finitos é limitada à e) As condições de contorno
análise de sistemas estruturais. não possuem influência
b) O método das diferenças significativamente na precisão
finitas, desenvolvido a partir do modelo analítico adotado.
do método dos elementos 4. Analise as afirmativas a seguir,
finitos, permite a resolução de referentes à matriz de rigidez, e
equações diferenciais parciais. assinale a alternativa correta.
c) O método dos elementos a) A análise matricial do método
finitos pode ser utilizado em dos elementos finitos é baseada
ambientes onde ocorram no método da rigidez apenas.
mudanças no meio. b) A matriz de rigidez representa
d) O método das diferenças uma relação linear entre os
finitas é indicado para a análise esforços e os deslocamentos.
de sistemas estruturais com c) O método dos deslocamentos
geometria irregular, quando é o único que permite a
comparado ao método montagem adequada da
dos elementos finitos. matriz de rigidez utilizada no
e) No método das diferenças método dos elementos finitos.
finitas, os sistemas estruturais d) A relação linear entre os esforços
são divididos em regiões de e os deslocamentos permite que
formas simples, permitindo seja obtida uma solução exata
que estas sejam analisadas
202 Conceitos básicos de métodos dos elementos finitos
para barras de seção variável b) O método dos elementos finitos
submetidas a cargas pontuais. deve ser aplicado na análise do
e) A análise matricial do método sistema estrutural como um todo.
dos elementos finitos permite c) O sistema de equações lineares
a obtenção de resultados pode ser representado pelo
exatos, independente da produto de matrizes AB = X.
geometria da barra. d) A análise pelo método dos
5. Analise as afirmativas a seguir, elementos finitos é realizada
referentes a sistemas lineares, e de forma manual, em função
assinale a alternativa correta. de não haverem métodos
a) O método dos elementos computacionais adequados
finitos permite transformar para sua resolução.
uma “superfície” em um e) A solução do sistema de
sistema linear de equações. equações lineares é dada pela
matriz coluna X, resultante
do produto AX = B.
BEER, F. P. et al. Estática e mecânica dos materiais. Porto Alegre: AMGH, 2013.
BEER, F. P. et al. Mecânica dos materiais. 7. edição. Porto Alegre: AMGH, 2015.
LEET, K. M.;UANG, C.; GILBERT, A. M. Fundamentos da análise estrutural. 3. ed. Porto
Alegre: AMGH, 2010.
Leituras recomendadas
HIBBELER, R. C. Análise das estruturas. 8. ed. São Paulo: Pearson Prentice Hall, 2013.
ILKIU, A. M. Teoria das estruturas: parte I. [S.l.]: [s.n.], 1998. Notas de aula.
TIMOSHENKO, S. P.; GERE, J. E. Mecânica dos sólidos: volume I.. Rio de Janeiro: LTC, 1983.
Encerra aqui o trecho do livro disponibilizado para
esta Unidade de Aprendizagem. Na Biblioteca Virtual
da Instituição, você encontra a obra na íntegra.
Conteúdo: