UNIVERSIDADE FEDERAL DE PERNAMBUCO

CENTRO DE TECNOLOGIA E GEOCINCIAS

PROGRAMA DE PS-GRADUAO EM ENGENHARIA CIVIL

THIAGO JOS SOUZA CORDEIRO

LISTA 09 Anlise de viga sobre base elstica com carregamento

distribudo Anlises com MATLAB e Ftool

Recife - PE
2016

THIAGO JOS SOUZA CORDEIRO

LISTA 09 Anlise de viga sobre base elstica com carregamento

distribudo Anlises com MATLAB e Ftool

Lista 09 da disciplina Mtodos

matemticos para engenharia, do Programa
de Ps-Graduao em Engenharia Civil da
UFPE, lecionada pelo Prof. Dr. Paulo
Marcelo V. Ribeiro.

Recife - PE
2016

Lista 09
QUESTO: Pede-se para elaborar um cdigo utilizando o MATLAB (utilizando o
mtodo das diferenas finitas) para determinao dos valores de deflexes de vigas
sobre base elstica nas situaes a) biapoiada e b) com as extremidades livres. Alm
disso, pede-se para comparar os resultados com os valores obtidos no Ftool. Por fim,
tambm foi pedido para verificar o que ocorre quando K igual a zero (caso de viga
biapoiada).
a) Viga biapoiada sobre base elstica:
Figura 1 Viga biapoiada sobre base elstica

Fonte: Autor

b) Viga com extremidades livres sobre base elstica:

Figura 2 Viga com extremidades livres sobre base elstica

Fonte: Autor

Viga com biapoiada sobre base elstica

a.1) Processo de obteno do cdigo

Para a elaborao de um cdigo do MATLAB utilizando o Mtodo das Diferenas
Finitas, para o clculo das deflexes em uma viga biapoiada sobre base elstica, procedeu-se
da maneira que ser descrita nesta seo.
A malha de diferenas finitas a ser considerada ter a forma vista na figura 3:
Figura 3 Malha de diferenas finitas

Fonte: Autor

importante mencionar que o MATLAB trabalha com elementos em forma de

matrizes. Portanto, a demonstrao que ser mostrada tender a ser nessa forma. Com isso,
sabe-se que os termos das matrizes so expressos em funo de seus sub-ndices, os quais
indicam a posio do elemento na matriz. Os sub-ndices, por sua vez, so nmeros naturais
maiores ou iguais a um. Portanto, a deduo ser feita levando em conta esse fato.
Para a resoluo de vigas apoiadas sobre base elstica atravs do Mtodo das
Diferenas Finitas, sero necessrias as seguintes equaes (j adaptadas forma como elas
sero utilizadas no MATLAB):

Equao diferencial para o estudo das vigas sobre base elstica, considerando a
rigidez flexo e o carregamento aplicado constantes:
. + . =

(1)

Relaes do Mtodo das Diferenas Finitas (todas com erro da ordem h):
()

() =
() =
() =
() =

2 ()

3 ()
3
4 ()
4

=
=

=
=

(+1 1 )
2
(+1 2 +1 )
2

(+2 2+1 +21 22 )

23
(+2 4+1 +6 41 +2 )
4

(2)
(3)
(4)
(5)

Com as equaes apresentadas anteriormente, torna-se possvel elaborar o cdigo

desejado. Em sala de aula, foi apresentada o seguinte caso, visto na Figura 4:
Figura 4 Exemplo descrito em sala de aula

Fonte: Autor

Nota-se, ento, que para aplicar a equao (1) no exemplo da figura 4, torna-se
necessria a introduo de pontos virtuais/pontos-fantasmas. Tal artifcio mostrado na
figura 5.
Figura 5 Introduo de pontos virtuais para auxiliar a resoluo

Fonte: Autor

Portanto, tais pontos virtuais sero, por meio das condies de contorno, relacionados
aos pontos internos, de tal modo que o clculo seja possvel e a matriz seja formada. Ento,
os pontos 1 e 5, por serem apoios, podem ser expressos pelas seguintes condio de contorno:

Ponto 1:
1 = 0
1 = 0

Ponto 5:
5 = 0
5 = 0

Ento, aplicando a equao (3), chega-se s seguintes expresses:

1 =
5 =

2 (1)
2
2 (5)
2

(2 21 +0 )
2
(6 25 +4 )
2

= 0 0 = 2

(6)

= 0 6 = 4

(7)

Com as expresses obtidas em (6) e (7), pode-se dar prosseguimento ao

desenvolvimento da questo. Da figura 5, sero utilizados, inicialmente, os pontos internos
em que h as molas, ou seja, os pontos 2, 3 e 4, como os pivs. Obtm-se, ento:

Piv no ponto 2:
(10 41 + 62 43 + 14 )
+ 2 =
4

Piv no ponto 3:
(11 42 + 63 44 + 15 )
+ 3 =
4

Piv no ponto 4:
(12 43 + 64 45 + 16 )
+ 4 =
4
possvel perceber, ento, que as expresses seguem um determinado padro. esse

padro que ser aproveitado no cdigo a ser criado. Das expresses dos pivs mostrados,
notrio que os pontos 2 e 4, que so os pontos imediatamente vizinhos dos pontos dos apoios,
dependem dos valores dos pontos virtuais. Isso algo que tambm se repete como um padro.
Portanto, introduzindo nas equaes dos pivs 2 e 4 as expresses obtidas em (6) e (7), temse:

Piv no ponto 2:
(2 41 + 62 43 + 14 )
+ 2 =
4
(41 + 52 43 + 14 )
+ 2 =
4

Piv no ponto 4:
(12 43 + 64 45 14 )
+ 4 =
4
(12 43 + 54 45 )
+ 4 =
4
Os pontos extremos podem ser considerados por meio das condies de apoio. O

procedimento adotado foi o de colocar os pivs nos pontos 1 e 5 e considerar os pontos

virtuais -1 e 7, como pode ser visto na figura 5.
De forma semelhante apresentada anteriormente, possvel chegar s seguintes
expresses:
1 = 3

(8)

7 = 3

(9)

Aplicando-se a expresso (1) aos pontos 1 e 5, possvel chegar a expresses

semelhantes s dos outros pontos. Tem-se, aps a realizao dos clculos:

Piv no ponto 1:
(61 )
+ 1 =
4

Piv no ponto 5:
(65 )
+ 5 =
4
As relaes encontradas podem ser reorganizadas, colocando-se em evidncias os

termos em comum e, por fim, escritas em forma matricial, conforme pode ser visto na figura
6:
Figura 6 Matriz que representa a viga biapoiada sobre base elstica

Fonte: Autor

Dessa forma, a partir dessa matriz e do desenvolvimento utilizado at ento, possvel

a elaborao do cdigo, o qual poder considerar mais pontos. Para isso, a figura 6 ser
modificada, mostrando alguns padres a serem considerados no cdigo em questo (figura
7). Notar que n no representa o nmero de pontos, mas sim o nmero de segmentos em
que a viga ser discretizada (por exemplo, na viga em questo, n = 4).
Figura 7 Observaes importantes para elaborao do cdigo

Fonte: Autor
Como j foi dito anteriormente, uma das condies de contorno para cada apoio que
a deflexo nesses pontos nula. Portanto, sero substitudos na matriz da figura 7 os seguintes
valores:
1 = 0
5 = 0
Optou-se por deixar os termos 1 e 5 sem serem colocados diretamente na matriz
pois a mesma poder ser facilmente editada no caso da considerao de recalque em algum
dos apoios, por exemplo.
O cdigo criado , ento, mostrado no item seguinte.
Observao: Nas frmulas aqui deduzidas, foi utilizado, como distncia entre os
pontos, a letra h. No entanto, para no confundir com a altura da seo, a qual ser pedida
no cdigo, tal termo ser chamado de dx.

a.2) Cdigo criado

%--------------- VIGA BIAPOIADA SOBRE BASE ELSTICA - MDF ---------------% ------------------ DADOS DE ENTRADA DO PROBLEMA ------------------------% PARMETROS GEOMTRICOS:
b=input('Entre com o valor da largura da viga, em cm:');
h=input('Entre com o valor da altura da viga, em cm:');
L=input('Entre com o valor do comprimento da viga, em m:');
I=(b/100)*((h/100)^3)/12;
% MDULO DE ELASTICIDADE DO MATERIAL:
E1=input('Entre com o valor do mdulo de elasticidade do material, em GPa:');
E = E1*10^6;
% PARMETRO GEOTCNICO:
kv=input('Entre com o valor da constante de mola do solo, em kN/m:');
k=kv*(b/100);
% CARREGAMENTO ATUANTE:
q=input('Entre com o valor do carregamento linear distribudo, em kN/m:');
% NMERO DE ELEMENTOS PARA A DISCRETIZAO DA VIGA:
n=input('Entre com o nmero de elementos para discretizar a viga:');
% VALOR DO INCREMENTO - dx
dx = L/n;
% MONTAGEM DA EQUAO MATRICIAL PARA OBTENO DA LINHA
ELSTICA:
% V = inv[M1]*M2
% A notao matricial adotada foi a mesma utilizada em exerccio
% anterior. A matriz M1 ser quadrada e de ordem n+1, uma vez que
% esto sendo considerados os apoios.
M1 = zeros(n+1,n+1)
M2 = zeros(n+1,1)

% Atravs do comando "for", criou-se um lao para calcular a regio

% interna da matriz M1. Notar que o intervalo representa a regio
% compreendida entre os pontos vizinhos aos apoios, conforme a figura
% j mostrada no texto.
for i = 3:n-1
M1(i,i) = 6*E*I+(k*dx^4);
M1(i,i-1) = -4*E*I;
M1(i,i+1) = -4*E*I;
M1(i,i-2) = E*I;
M1(i,i+2) = E*I;
M2(i,1) = q*dx^4;
end
% Alm disso, necessrio completar a matriz com os pontos referentes
% aos apoios e os pontos imediatamente vizinhos a estes. Para os pon% tos referentes ao apoio da esquerda, tem-se:
M1(1,1) = 6*E*I+(k*dx^4)
M1(2,2) = 5*E*I+(k*dx^4)
M1(2,3) = -4*E*I
M1(2,4) = E*I
M2(1,1) = q*dx^4;
% Para o apoio direita, procedeu-se de forma anloga ao caso anterior:
M1(n,n) = 5*E*I+(k*dx^4)
M1(n,n+1) = -4*E*I
M1(n,n-1) = -4*E*I
M1(n,n-2) = E*I
M2(n+1,1) = q*dx^4;
M1(n+1,n+1) = 6*E*I+(k*dx^4)
% A matriz com os deslocamentos para os pontos da viga dada por:
V = inv(M1)*M2
% As condies de contorno so, ento, colocadas:
V(1,1) = 0
V(n+1,1) = 0
V
% Por fim, a curva da linha elstica pode ser obtida:
X = 0:dx:L
plot (X,-V)

xlabel('Comprimento da viga - m')

ylabel('Deflexo da viga (linha elstica) - m')
title ('Linha elstica de viga biapoiada em base elstica')
% Sendo a amplitude mxima do deslocamento, em mm, igual a:
Deflexao_max_mm = max(abs(V))*1000
% COMENTRIO: Por ser um mtodo numrico que depende de condies
% de contorno, o mtodo das diferenas finitas acaba por ser limi% tado nesse aspecto. Para que seus resultados dem maior preciso,
% necessrio que o valor de "n" aqui pedido seja considervel,
% de tal modo que o valor do incremento "dx" seja suficientemente
% pequeno. Com isso, torna-se possvel obter os valores de deflexo
% nas vigas com uma excelente aproximao.

a.3) Resultados obtidos com o cdigo para o exemplo dado

Para os valores fornecidos pelo professor, os quais podem ser vistos na figura 4, foram
colocadas as unidades a serem consideradas. Portanto, os dados iniciais do problema foram:

b = 20 cm;

h = 50 cm;

L = 5 m;

E = 20 GPa;

Kv = 10000 kN/m;

q = 100 kN/m;
Com tais termos, foram gerados o diagrama da linha elstica da viga sobre base

elstica (figura 8), bem como calculado o seu valor mximo (figura 9).
Figura 8 Linha elstica da viga sobre base elstica

Fonte: Autor
Figura 9 Deflexo mxima obtida, em mm

Fonte: Autor

a.4) Modelagem do problema no Ftool

O problema foi, ento, modelado no Ftool com 3 valores de dx diferentes: dx = 0,25
m, dx = 0,50 m e dx = 1,00 m, cujas linhas elsticas e deflexes mximas so mostradas,
respectivamente, da figura 10 at a 15 (foram feitas as devidas alteraes nos coeficientes de
molas para considerar a variao do dx, para cada caso).
Figura 10 Linha elstica dx = 0,25 m

Fonte: Autor
Figura 11 Deflexo mxima dx = 0,25 m

Fonte: Autor

Figura 12 Linha elstica dx = 0,50 m

Fonte: Autor

Figura 13 Deflexo mxima dx = 0,50 m

Fonte: Autor

Figura 14 Linha elstica dx = 1,00 m

Fonte: Autor

Figura 15 Deflexo mxima dx = 1,00 m

Fonte: Autor

a.5) Comparao dos resultados obtidos

A tabela 1 foi, ento, preparada para resumir todos os valores obtidos at ento,
comparando-os.
Tabela 1 Resumos dos valores de deflexo obtidos

Fonte: Autor

Percebe-se, com isso, que os valores obtidos com o cdigo do MATLAB foram muito
semelhantes aos obtidos com o Ftool, sendo a diferena menor que 1% quando foram
utilizados os intervalos de dx= 0,25 m e dx = 0,50 m entre as molas.

a.6) Anlise extra: Caso em que K = 0 (ou seja, em que as molas no so

consideradas)
Utilizando-se o cdigo mostrado anteriormente para o caso em que K = 0 (declara-se
no cdigo, no caso, kv = 0, uma vez que K = kv.b), foi possvel obter a linha elstica (figura
16) e o valor da deflexo mxima (figura 17) encontrados:

Figura 16 Linha elstica obtida K = 0

Fonte: Autor
Figura 17 Deflexo mxima obtida, em mm, para K=0

Fonte: Autor

A deflexo mxima obtida pode ser comparada ao valor obtido atravs da frmula
analtica (facilmente encontrada em livros de Resistncia dos Materiais) e vista em (10).
=

5..4
384.

(10)

Com os valores j mostrados no incio da seo a.3, tem-se:

=

510054
38420106

1
480

= 0,01953125 19,5312

Portanto, a diferena entre os resultados obtidos com o cdigo e a soluo analtica

fica em torno de:
19,5312
= 1,00087629 0,09%
19,5141

a.7) Concluses parciais

A partir das diferenas percentuais analisadas, possvel constatar que o cdigo
criado apresenta excelente preciso para o clculo de deflexes em vigas biapoiadas sobre
base elstica.

Viga com extremidades livres sobre base elstica

b.1) Processo de obteno do cdigo

O processo de obteno do cdigo para a viga com extremidades livres sobre base
elstica bastante semelhante ao j descrito anteriormente para vigas biapoiadas. Por isso,
ser analisado um caso com apenas 5 pontos, e, em seguida, o raciocnio expandido para o
caso de n+1 pontos, em que n o nmero de segmentos de discretizao das vigas. Para tal,
foram preparadas as figuras 18 e 19.
Figura 18 Viga com bordos livres apoiada sobre base elstica

Fonte: Autor
Figura 19 Discretizao da viga de bordos livres sobre base elstica

Fonte: Autor

Como se pode ver, so considerados 4 pontos virtuais para esse caso. As equaes
que sero aqui utilizadas j foram mostradas (Equaes (1) a (5)). As condies de contorno
desse problema, por sua vez, so diferentes do caso anteriormente analisado. So elas:

Ponto 1:
1 = 0
1 = 0

Ponto 5:
5 = 0
5 = 0
Ento, aplicando-se as equaes do Mtodo das Diferenas Finitas, possvel chegar

s seguintes expresses:

Ponto 1:
1 = 0
1 =

0 21 + 2
=0

0 = 21 2

(11)

1 = 0

1 =

11 + 20 + 01 22 + 13
=0
2()

Utilizando-se a expresso obtida em (11), chega-se seguinte expresso:

1 = 41 42 + 13

(12)

Ponto 5:
5 = 0
5 =

4 25 + 6
=0

6 = 25 4

(13)

5 = 0
5 =

13 + 24 + 05 26 + 17
=0
2()

Substituindo-se a expresso obtida em (13), chega-se a:

7 = 45 44 + 13

(14)

Aps a obteno das expresses vistas em (11), (12), (13) e (14), pode-se calcular,
agora, os pivs dos pontos entre os apoios, ou seja, os pontos 2, 3 e 4:

Piv no ponto 2:
(10 41 + 62 43 + 14 )
+ 2 =
4

Piv no ponto 3:
(11 42 + 63 44 + 15 )
+ 3 =
4

Piv no ponto 4:
(12 43 + 64 45 + 16 )
+ 4 =
4
Como se v, as expresses, num primeiro momento, so absolutamente iguais s

obtidas no caso das vigas biapoiadas. No entanto, para os pontos vizinhos aos extremos, ou
seja, os pontos 2 e 4 do exemplo discutido, haver diferenas para o caso biapoiado porque
os valores de 0 e 6 no sero mais os mesmos. No entanto, o valor do piv 3 permanece
inalterado e esse um dos padres para que o cdigo das vigas biapoiadas seja reaproveitado.
Portanto, para os pivs 2 e 4, teremos as seguintes equaes (ao substituir a equao
(11) no piv 2 e a equao (13) no piv 4):

Piv no ponto 2:
(10 41 + 62 43 + 14 )
+ 2 =
4
(21 2 41 + 62 43 + 14 )
+ 2 =
4
(21 + 52 43 + 14 )
+ 2 =
4

Piv no ponto 4:
(12 43 + 64 45 + 16 )
+ 4 =
4

 (12 43 + 64 45 + 25 4 )
+ 4 =
4
(12 43 + 54 25 )
+ 4 =
4
Portanto, j se tm 3 equaes (pivs 2, 3 e 4). O raciocnio utilizado ser, agora, feito
para os pontos dos extremos. Ao substituir as expresses (11) e (12) no piv 1 e as expresses
(13) e (14) no piv 5, obtm-se as seguintes expresses:

Piv no ponto 1:
(21 42 + 23 )
+ 1 =
4

Piv no ponto 5:
(23 45 + 25 )
+ 5 =
4
As relaes encontradas podem ser reorganizadas, colocando-se em evidncias os

termos em comum e, por fim, escritas em forma matricial, conforme pode ser visto na figura
20:
Figura 20 Matriz encontrada para o caso de vigas com bordas livres sobre base elstica

Fonte: Autor

De forma semelhante apresentada no item anterior, a figura 20 ser modificada de

tal modo a identificar os padres presentes na mesma e que serviro de base para o cdigo a
ser criado (figura 21).

Figura 21 Observaes interessantes sobre o comportamento da matriz

Fonte: Autor

Com as informaes mostradas anteriormente, pode-se expandir o modelo

apresentado para n+1 pontos e fazer o cdigo no MATLAB. O cdigo final ser mostrado no
item seguinte.

b.2) Cdigo criado

%------- VIGA DE BORDOS INFINITOS SOBRE BASE ELSTICA - MDF ------% ------------------ DADOS DE ENTRADA DO PROBLEMA ------------------------% PARMETROS GEOMTRICOS:
b=input('Entre com o valor da largura da viga, em cm:');
h=input('Entre com o valor da altura da viga, em cm:');
L=input('Entre com o valor do comprimento da viga, em m:');
I=(b/100)*((h/100)^3)/12;
% MDULO DE ELASTICIDADE DO MATERIAL:
E1=input('Entre com o valor do mdulo de elasticidade do material, em GPa:');
E = E1*10^6;
% PARMETRO GEOTCNICO:
kv=input('Entre com o valor da constante de mola do solo, em kN/m:');
k=kv*(b/100);
% CARREGAMENTO ATUANTE:
q=input('Entre com o valor do carregamento linear distribudo, em kN/m:');
% NMERO DE ELEMENTOS PARA A DISCRETIZAO DA VIGA:
n=input('Entre com o nmero de elementos para discretizar a viga:');
% VALOR DO INCREMENTO - dx
dx = L/n;
% MONTAGEM DA EQUAO MATRICIAL PARA OBTENO DA LINHA
ELSTICA:
% V = inv[M1]*M2
% A notao matricial adotada foi a mesma utilizada em exerccio
% anterior. A matriz M1 ser quadrada e de ordem n+1, uma vez que
% esto sendo considerados os apoios.
M1 = zeros(n+1,n+1)
M2 = zeros(n+1,1)

% Atravs do comando "for", criou-se um lao para calcular a regio

% interna da matriz M1. Notar que o intervalo representa a regio
% compreendida entre os pontos vizinhos aos apoios, conforme a figura
% j mostrada no texto.
for i = 3:n-1
M1(i,i) = 6*E*I+(k*dx^4);
M1(i,i-1) = -4*E*I;
M1(i,i+1) = -4*E*I;
M1(i,i-2) = E*I;
M1(i,i+2) = E*I;
M2(i,1) = q*dx^4;
end
% Alm disso, necessrio completar a matriz com os pontos referentes
% aos apoios e os pontos imediatamente vizinhos a estes. Para os pon% tos referentes borda livre esquerda, tem-se:
M1(1,1) = 2*E*I+(k*dx^4)
M1(1,2) = -4*E*I
M1(1,3) = 2*E*I
M1(2,1) = -2*E*I
M1(2,2) = 5*E*I+(k*dx^4)
M1(2,3) = -4*E*I
M1(2,4) = E*I
M2(1,1) = q*dx^4;
% Para a borda livre direita, procedeu-se de forma anloga ao caso anterior:
M1(n,n) = 5*E*I+(k*dx^4)
M1(n,n+1) = -2*E*I
M1(n,n-1) = -4*E*I
M1(n,n-2) = E*I
M1(n+1,n+1) = 2*E*I+(k*dx^4)
M1(n+1,n) = -4*E*I
M2(n+1,1) = q*dx^4;
% A matriz com os deslocamentos para os pontos da viga dada por:
V = inv(M1)*M2
V
% Por fim, a curva da linha elstica pode ser obtida:
X = 0:dx:L
plot (X,-V)
xlabel('Comprimento da viga - m')

ylabel('Deflexo da viga (linha elstica) - m')

title ('Linha elstica de viga com bordos infinitos em base elstica')
% Sendo a deflexo, em mm, considerando um termo do centro da viga,
% igual a:
Deflexao_mm = V(n/2)*1000
% COMENTRIO: Por ser um mtodo numrico que depende de condies
% de contorno, o mtodo das diferenas finitas acaba por ser limi% tado nesse aspecto. Para que seus resultados dem maior preciso,
% necessrio que o valor de "n" aqui pedido seja considervel,
% de tal modo que o valor do incremento "dx" seja suficientemente
% pequeno. Com isso, torna-se possvel obter os valores de deflexo
% nas vigas com uma excelente aproximao.

b.3) Resultados obtidos com o cdigo para o exemplo dado

Os valores fornecidos pelo professor, os quais j so conhecidos (ver figura 18) so:

b = 20 cm;

h = 50 cm;

L = 40 m (mas a viga considerada na regio de 5 metros central, a qual ser

destacada em verde na figura 24);

E = 20 GPa;

Kv = 10000 kN/m;

q = 100 kN/m;
Com tais termos, foram gerados o diagrama da linha elstica da viga com bordos

livres sobre base elstica (figura 22), bem como calculado a sua deflexo (figura 23).
Figura 22 Formato da linha elstica Bordos livres

Fonte: Autor

Figura 23 Deflexo obtida

Fonte: Autor

Da figura 22 possvel fazer uma anlise mais aprofundada (vista na figura 24) e ver
a influncia das condies de apoio na curva de deslocamentos obtida. A viga foi considerada
com 40 metros de tal modo que o seu comprimento fosse consideravelmente maior que sua
altura (no caso, a altura utilizada foi de 50 cm). A viga de 5 metros de comprimento pode ser
considerada no trecho entre as abscissas 17,5 m e 22,5 m.
Figura 24 Interferncia dos bordos (vermelho) e viga do exemplo considerado (verde)

Fonte: Autor

b.4) Modelagem do problema com o Ftool

O problema foi, ento, modelado no Ftool com 2 valores de dx diferentes: dx = 1,00
m, e dx = 2,00 m, cujas linhas elsticas e deflexes mximas so mostradas, respectivamente,
da figura 25 at a 28 (foram feitas as devidas alteraes nos coeficientes de molas para
considerar a variao do dx, sendo que, nos extremos, eles foram considerados com metade
do valor que possuem nos pontos internos).
Figura 25 Linha elstica dx = 1,00 m

Fonte: Autor
Figura 26 Deflexo encontrada dx = 1,00 m

Fonte: Autor

Figura 27 Linha elstica dx = 2,00 m

Fonte: Autor
Figura 28 Deflexo encontrada dx = 2,00 m

Fonte: Autor

Aps a obteno dos resultados, eles so comparados no item seguinte.

b.5) Comparao dos resultados obtidos

A tabela 2 foi, ento, preparada para resumir todos os valores de deflexo obtidos nas
duas sees anteriores, relacionando-os.
Tabela 2 Resumo dos valores de deflexo obtidos

Fonte: Autor

Percebe-se, com isso, que os valores obtidos com o cdigo do MATLAB foram muito
semelhantes aos obtidos com o Ftool, sendo a diferena menor que 1% para as duas situaes
analisadas.

b.6) Concluses parciais

A partir das comparaes do item anterior, notvel que o cdigo criado apresenta
resultados satisfatrios.

CONCLUSES
Tendo em vista o que foi apresentado, possvel constatar que os cdigos criados
deram resultados de boa preciso. Para isso, utilizou-se um valor n (nmero de segmentos)
igual a 100 em todas as situaes, de tal modo que o valor do dx ficasse muito pequeno, o
que, de certa forma, melhorou a preciso dos valores encontrados (e isso foi constatado no
texto). Outra vantagem do mtodo aqui exposto a facilidade de se considerar recalques
existentes em qualquer um dos apoios/extremos.

REFERNCIA BIBLIOGRFICA:
TIMOSHENKO, S. P.; GERE, J. E. Mecnica dos Slidos, Vol. 1, Livros Tcnicos e
Cientficos, Rio de Janeiro, 1994.