UNIVERSIDAD NACIONAL DE

INGENIERÍA
FACULTAD DE INGENIERÍA MECÁNICA
Especialidad de Ingeniería Mecánica-Eléctrica

CALCULO POR ELEMENTOS FINITOS

(MC516 - F)
Práctica N°5: Vigas
PRESENTADO POR:

Apellidos y nombres Código

Mayta Aliaga Diego José 20184029A

PROFESOR
Cueva Pacheco, Ronald
LIMA, 2021

Página | 0
INDICE

1. DESARROLLO DEL PROBLEMA ........................................................................2

1.1 Enunciado del problema ..........................................................................................2

1.2 Modelado del cuerpo................................................................................................3

1.3 Grados de libertad nodales .......................................................................................4

1.4 Vector carga .............................................................................................................4

1.5 Matrices de Rigidez Locales ....................................................................................5

1.6 Matriz de Rigidez del Modelo .................................................................................7

1.7 Ecuación de Rigidez y condiciones de contorno .....................................................8

1.8 Gráfica del momento flector ....................................................................................8

1.9 Gráfica de la fuerza cortante ....................................................................................9

1.10 Deformada de la viga .............................................................................................9

1.11 Cálculo del esfuerzo máximo.................................................................................9

1.12 Resultados ............................................................................................................10

2. DIAGRAMA DE FLUJO DEL PROGRAMA ......................................................12

3. DIGITACIÓN Y EJECUCIÓN (MATLAB) .........................................................15

3.1 Código de Matlab ...................................................................................................15

3.2 Vista en el “Command Window”...........................................................................20

4. CONCLUSIONES ....................................................................................................22

Página | 1
1. DESARROLLO DEL PROBLEMA

1.1 Enunciado del problema

Utilizando el CÁLCULO POR ELEMENTOS FINITOS, y eligiendo el número de
elementos finitos para el modelo de la viga mostrada.
DATOS:
𝑁
𝐸 = 3 ∙ 105
𝑚𝑚2

HALLAR:
• Las fuerzas de reacción en los apoyos
• Los gráficos de momento flector y fuerza cortante
• Los esfuerzos máximos en cada elemento finito (𝜎𝑥 )
• Las pendientes de la deformada en (x=0) y en (x=L 1)
• La deformada de la viga

FIGURA:
𝑴𝑨 = 𝟐𝟎𝟎 𝟎𝟎𝟎 𝑵 − 𝒎𝒎

𝟕𝟎𝟎 𝐦𝐦 𝟑𝟎𝟎 𝐦𝐦

𝟒𝟎 𝐦𝐦
𝟐𝟎 𝐦𝐦

𝟖𝟎 𝐦𝐦

𝟐𝟎 𝐦𝐦

𝟏𝟎𝟎 𝐦𝐦

Página | 2
1.2 Modelado del cuerpo
Para este análisis se consideró un sistema de 4 elementos finitos, 5 nodos y 10 grados de
libertad.

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

1 2 3 4

𝟐𝟓𝟎 𝐦𝐦 𝟐𝟓𝟎 𝐦𝐦 𝟐𝟎𝟎 𝐦𝐦 𝟑𝟎𝟎 𝐦𝐦

Tabla de conectividad del modelo

A continuación, se muestra la tabla de conectividad de la viga.

Nodos GDL
(e) le (mm) Ae (mm2) Ee (N/mm2)
(1) (2) 1 2 3 4

1 1 2 Q1 Q2 Q3 Q4 250 3600 3 ∙ 105

2 2 3 Q3 Q4 Q5 Q6 250 3600 3 ∙ 105

3 3 4 Q5 Q6 Q7 Q8 200 3600 3 ∙ 105

4 4 5 Q7 Q8 Q9 Q10 300 3600 3 ∙ 105

Página | 3
 1.3 Grados de libertad nodales

(Q1) (Q3) (Q5) (Q7) (Q9)

(Q2) (Q4) (Q6) (Q8) (Q10)

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

El sistema presenta 10 grados de libertad. Considerando las condiciones de contorno

(apoyos en 1 y 5), el vector Qj quedaría de la siguiente manera:

𝑸𝒋 = [𝟎 𝑸𝟐 𝑸𝟑 𝑸𝟒 𝑸𝟓 𝑸𝟔 𝑸𝟕 𝑸𝟖 𝟎 𝑸𝟏𝟎 ]

1.4 Vector Carga

(F1) (F3) (F5) (F7) (F9)

(F2) (F4) (F6) (F8) (F10)

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

Considerando las condiciones de contorno (apoyos en 1 y 5 y el momento concentrado

en el nodo 4), el vector Fi quedaría de la siguiente manera:

𝑭𝒊 = [𝑹𝟏 𝟎 𝟎 𝟎 𝟎 𝟎 𝟎 𝟐𝟎𝟎𝟎𝟎𝟎 𝑹𝟗 𝟎]

Página | 4
1.5 Matrices de Rigidez Locales

Las matrices de rigidez locales para cada elemento finito, se expresan como:

Para el sistema estudiado, ya que está conformado por 4 elementos finitos, se tendrían 4
matrices de rigidez locales.

Cálculo del momento de inercia de la sección con respecto al eje neutro

En primer lugar, hallamos la ubicación del eje que contiene al centroide y al eje neutro:

𝟒𝟎 𝐦𝐦 Se divide el área y se aplica la siguiente

𝟐𝟎 𝐦𝐦 relación:

10 ∙ (100 ∙ 20) + 60 ∙ (80 ∙ 20)

𝑌𝑐 =
(100 ∙ 20 + 80 ∙ 20)
𝟖𝟎 𝐦𝐦
𝟐𝟗𝟎
𝒀𝒄 = = 𝟑𝟐. 𝟐𝟐𝟐𝟐 𝒎𝒎
𝟗

𝟐𝟎 𝐦𝐦 Por lo tanto, el centroide se encuentra

aproximadamente a 32.22 mm del borde
inferior.
𝟏𝟎𝟎 𝐦𝐦

A continuación, hallamos el momento de inercia con respecto al eje neutro:

𝟒𝟎 𝐦𝐦 Se divide el área y se aplica la siguiente

𝟐𝟎 𝐦𝐦 relación:

1
𝐼𝑧′ = (2 ∙ 4 ∙ 23 + 2 ∙ 103 ) ∙ 104 = 688 ∙ 104 𝑚𝑚4
3

𝟖𝟎 𝐦𝐦 Luego, aplicando el teorema de Steiner:

𝒛
𝐼𝑧 = (688 − (8 ∙ 2 + 2 ∙ 10) ∙ 3.2222 ) ∙ 104

𝟐𝟎 𝐦𝐦
𝒛′ 𝑰𝒛 = 𝟑𝟏𝟒. 𝟐𝟐𝟐 ∙ 𝟏𝟎𝟒 𝒎𝒎𝟒

𝟏𝟎𝟎 𝐦𝐦

Página | 5
Cálculo de la matriz de rigidez de cada elemento finito

Para el elemento finito 1

𝟏𝟐 𝟏𝟓𝟎𝟎 −𝟏𝟐 𝟏𝟓𝟎𝟎

𝟏𝟓𝟎𝟎 𝟐𝟓𝟎𝟎𝟎𝟎 −𝟏𝟓𝟎𝟎 𝟏𝟐𝟓𝟎𝟎𝟎
𝒌𝟏 = 𝟔𝟎𝟑𝟑𝟎. 𝟔𝟔𝟔𝟔𝟕 ∙ [ ]
−𝟏𝟐 −𝟏𝟓𝟎𝟎 𝟏𝟐 −𝟏𝟓𝟎𝟎
𝟏𝟓𝟎𝟎 𝟏𝟐𝟓𝟎𝟎𝟎 −𝟏𝟓𝟎𝟎 𝟐𝟓𝟎𝟎𝟎𝟎

Para el elemento finito 2

𝟏𝟐 𝟏𝟓𝟎𝟎 −𝟏𝟐 𝟏𝟓𝟎𝟎

𝟏𝟓𝟎𝟎 𝟐𝟓𝟎𝟎𝟎𝟎 −𝟏𝟓𝟎𝟎 𝟏𝟐𝟓𝟎𝟎𝟎
𝒌𝟏 = 𝟔𝟎𝟑𝟑𝟎. 𝟔𝟔𝟔𝟔𝟕 ∙ [ ]
−𝟏𝟐 −𝟏𝟓𝟎𝟎 𝟏𝟐 −𝟏𝟓𝟎𝟎
𝟏𝟓𝟎𝟎 𝟏𝟐𝟓𝟎𝟎𝟎 −𝟏𝟓𝟎𝟎 𝟐𝟓𝟎𝟎𝟎𝟎

Para el elemento finito 3

𝟏𝟐 𝟏𝟐𝟎𝟎 −𝟏𝟐 𝟏𝟐𝟎𝟎

𝟏𝟐𝟎𝟎 𝟏𝟔𝟎𝟎𝟎𝟎 −𝟏𝟐𝟎𝟎 𝟖𝟎𝟎𝟎𝟎
𝒌𝟏 = 𝟏𝟏𝟕𝟖𝟑𝟑. 𝟑𝟑𝟑𝟑𝟑 ∙ [ ]
−𝟏𝟐 −𝟏𝟐𝟎𝟎 𝟏𝟐 −𝟏𝟐𝟎𝟎
𝟏𝟐𝟎𝟎 𝟖𝟎𝟎𝟎𝟎 −𝟏𝟐𝟎𝟎 𝟏𝟔𝟎𝟎𝟎𝟎

Para el elemento finito 4

𝟏𝟐 𝟏𝟖𝟎𝟎 −𝟏𝟐 𝟏𝟖𝟎𝟎

𝟏𝟖𝟎𝟎 𝟑𝟔𝟎𝟎𝟎𝟎 −𝟏𝟖𝟎𝟎 𝟏𝟖𝟎𝟎𝟎𝟎
𝒌𝟏 = 𝟑𝟒𝟗𝟏𝟑. 𝟓𝟖𝟎𝟐𝟓 ∙ [ ]
−𝟏𝟐 −𝟏𝟖𝟎𝟎 𝟏𝟐 −𝟏𝟖𝟎𝟎
𝟏𝟖𝟎𝟎 𝟏𝟖𝟎𝟎𝟎𝟎 −𝟏𝟖𝟎𝟎 𝟑𝟔𝟎𝟎𝟎𝟎

Página | 6
1.6 Matriz de Rigidez del Modelo (global)
Una vez calculadas las matrices, calculamos la matriz de rigidez global según la siguiente
relación:
∈
𝑲𝒊𝒋 = ∑ 𝒌𝒆𝒔𝒓
𝒆

Para el sistema, la matriz de rigidez global es de 10x10 y está dada por:

Columnas de la 1 a la 4:

Columnas de la 5 a la 8:

Columnas de la 9 a la 10:

Página | 7
1.7 Ecuación de Rigidez y condiciones de contorno
La ecuación de rigidez está determinada por la siguiente ecuación:

𝑭𝒊 = 𝑲𝒊𝒋 ∙ 𝑸𝒋

Al aplicar dicha relación, se obtiene el valor de los 10 grados de libertad del sistema.
Grados de libertad

𝑸𝟏 = 0 𝑚𝑚 𝑸𝟐 = −0.0000368761366 𝑟𝑎𝑑

𝑸𝟑 = −0.0084297332795 𝑚𝑚 𝑸𝟒 = −0.0000274045262 𝑟𝑎𝑑

𝑸𝟓 = −0.0121236613457 𝑚𝑚 𝑸𝟔 = 0.0000010103051 𝑟𝑎𝑑

𝑸𝟕 = −0.0084865629420 𝑚𝑚 𝑸𝟖 = 0.0000373812891 𝑟𝑎𝑑

𝑸𝟗 = 0 𝑚𝑚 𝑸𝟏𝟎 = 0.0000237421701 𝑟𝑎𝑑

Asimismo, se obtienen las reacciones en los apoyos.

𝑭𝟏 = 𝑹𝟏 = 200 𝑁
𝑭𝟗 = 𝑹𝟗 = −200 𝑁

Por último, se tiene el valor de las pendientes de la deformada en x=0 y x=1000 mm.

𝒙 = 𝟎 𝒎𝒎 → 𝑸𝟐 = −0.0000368761366 𝑟𝑎𝑑
𝒙 = 𝟏𝟎𝟎𝟎 𝒎𝒎 → 𝑸𝟏𝟎 = 0.0000237421701 𝑟𝑎𝑑

1.8 Gráfica del momento flector

Para obtener la gráfica del momento flector, se aplica la siguiente relación para cada
elemento finito (𝝃 ∈ < −𝟏, 𝟏 >):

Como resultado, se obtendrá la gráfica del momento flector dentro de cada elemento finito
en particular. Posteriormente, con ayuda del software Matlab, concatenaremos cada curva
obtenida y tendremos la gráfica del momento flector para todo el sistema. El gráfico
correspondiente se muestra en la sección de Resultados.

Página | 8
1.9 Gráfica de la fuerza cortante
Para obtener la gráfica de la fuerza cortante, se aplica la siguiente relación para cada
elemento finito (𝝃 ∈ < −𝟏, 𝟏 >):

Como resultado, se obtendrá la gráfica de la fuerza cortante dentro de cada elemento finito
en particular. Posteriormente, con ayuda del software Matlab, concatenaremos cada curva
obtenida y tendremos la gráfica de la fuerza cortante para todo el sistema. El gráfico
correspondiente se muestra en la sección de Resultados.

1.10 Deformada de la viga

Para obtener la deformada de la viga, se aplica la siguiente relación para cada elemento
finito (𝝃 ∈ < −𝟏, 𝟏 >):

Como resultado, se obtendrá la deformada dentro de cada elemento finito en particular.

Posteriormente, con ayuda del software Matlab, concatenaremos cada curva obtenida y
tendremos la deformada para todo el sistema. El gráfico correspondiente se muestra en la
sección de Resultados.

1.11 Cálculo del esfuerzo máximo

La expresión del esfuerzo para cada sección dentro del elemento finito (𝝃 ∈ < −𝟏, 𝟏 >)
esta dada por:

Considerando que el esfuerzo máximo en la sección se da en el extremo superior (𝑦𝑚𝑎𝑥 =

67.78 𝑚𝑚), se obtienen los esfuerzos máximos en cada elemento finito. Los valores
hallados se muestran en la sección de Resultados.

Página | 9
1.12 Resultados
a) Las fuerzas de reacción en los apoyos
Al aplicar la ecuación de Rigidez, se obtuvo el valor de las reacciones en los apoyos:

𝑭𝟏 = 𝑹𝟏 = 200 𝑁
𝑭𝟗 = 𝑹𝟗 = −200 𝑁

b) Los gráficos del momento flector y de la fuerza cortante

Con ayuda del software Matlab, se obtuvieron las siguientes gráficas:

Gráfica C-1. Gráfica de momento flector en la viga

Gráfica C-2. Gráfica de fuerza cortante en la viga

Página | 10
c) Los esfuerzos máximos en cada elemento finito
Considerando que el esfuerzo máximo en la sección se da en el extremo superior (𝑦𝑚𝑎𝑥 =
67.78 𝑚𝑚), se obtienen los esfuerzos máximos en cada elemento finito. Los resultados
fueron los siguientes:

𝑬𝒍𝒆𝒎𝒆𝒏𝒕𝒐 𝒇𝒊𝒏𝒊𝒕𝒐 𝟏 → −1.078500707213575 𝑀𝑃𝑎

𝑬𝒍𝒆𝒎𝒆𝒏𝒕𝒐 𝒇𝒊𝒏𝒊𝒕𝒐 𝟐 → −2.157001414427233 𝑀𝑃𝑎
𝑬𝒍𝒆𝒎𝒆𝒏𝒕𝒐 𝒇𝒊𝒏𝒊𝒕𝒐 𝟑 → −3.019801980197985 𝑀𝑃𝑎
𝑬𝒍𝒆𝒎𝒆𝒏𝒕𝒐 𝒇𝒊𝒏𝒊𝒕𝒐 𝟒 → −1.294200848656254 𝑀𝑃𝑎

Por lo tanto, el esfuerzo máximo en la viga se encuentra en el elemento finito 3:

𝑬𝒔𝒇𝒖𝒆𝒓𝒛𝒐 𝒎á𝒙𝒊𝒎𝒐 → −3.019801980197985 𝑀𝑃𝑎

d) Las pendientes de la deformada en (x=0) y en (x=1000 mm)

Al aplicar la ecuación de Rigidez, se obtuvo el valor de las pendientes de la deformada
en x=0 y x=1000 mm.

𝒙 = 𝟎 𝒎𝒎 → 𝑸𝟐 = −0.0000368761366 𝑟𝑎𝑑
𝒙 = 𝟏𝟎𝟎𝟎 𝒎𝒎 → 𝑸𝟏𝟎 = 0.0000237421701 𝑟𝑎𝑑

e) La deformada de la viga
Con ayuda del software Matlab, se obtiene la deformada:

Gráfica C-3. Deformada de la viga

Página | 11
2. DIAGRAMA DE FLUJO DEL PROGRAMA

Página | 12
Página | 13
Página | 14
3. DIGITACIÓN Y EJECUCIÓN (MATLAB)

3.1 Código de Matlab

clc
clear all

disp('................................................................
.............')
disp(' INGRESO DE DATOS')
disp('----------------------------------------------------------------
-------------')

E=input('Ingrese el valor del Módulo de elasticidad en MPa: ');

k=input('Ingrese el número de elementos finitos: ');
MA=input('Ingrese el valor del momento MA en N-mm: ');
nodo_MA=input('Ingrese el nodo en el que se ubica el momento MA: ');
I=input('Ingrese el momento de inercia de la viga con respecto al eje
neutro en mm4: ');
ymax=input('Ingrese el valor del ymax de la sección de la viga con
respecto al eje neutro en mm: ');
fprintf('\n')

%INGRESO DE LA LONGITUD DE CADA ELEMENTO FINITO

l=zeros(k,1);
for i=1:k
fprintf('Ingrese la longitud del elemento %.0f en mm : \n',i)
le = input('l = ');
l(i) = le;
end
l; %Vector columna que contiene las longitudes de los elementos
finitos
fprintf('\n')

%CONDICIONES DE CONTORNO

desplazamientos=zeros(2*(k+1),1);
desplazamientos(1)=0; %debido al apoyo simple 1
desplazamientos(2*k+1)=0; %debido al apoyo simple 2

%DETERMINACIÓN DEL VECTOR FUERZA

fuerzas=zeros(2*(k+1),1);
fuerzas(2*nodo_MA)=MA; %debido al apoyo simple 1

%HALLAMOS LA MATRIZ DE RIGIDEZ LOCAL DE CADA ELEMENTO

ksr=[];
for i=1:k
ks = (I*E/(l(i)^3))*[12 6*l(i) -12 6*l(i);
6*l(i) 4*l(i)*l(i) -6*l(i) 2*l(i)*l(i);
-12 -6*l(i) 12 -6*l(i);
6*l(i) 2*l(i)*l(i) -6*l(i) 4*l(i)*l(i)];
ksr=[ksr ks];
end

ksr; %Vector que guarda todas las matrices de rigidez locales

Página | 15
%HALLAMOS LA MATRIZ DE RIGIDEZ GLOBAL

krs=zeros(2*(k+1));
Kij=zeros(2*(k+1));

for i=1:k
ps1=i*2-1;
ps2=i*2;
ps3=(i+1)*2-1;
ps4=(i+1)*2;

krs(ps1,ps1)=12;
krs(ps1,ps2)=l(i)*6;
krs(ps1,ps3)=-12;
krs(ps1,ps4)=l(i)*6;

krs(ps2,ps1)=l(i)*6;
krs(ps2,ps2)=4*l(i)^2;
krs(ps2,ps3)=-l(i)*6;
krs(ps2,ps4)=2*l(i)^2;

krs(ps3,ps1)=-12;
krs(ps3,ps2)=-l(i)*6;
krs(ps3,ps3)=12;
krs(ps3,ps4)=-l(i)*6;

krs(ps4,ps1)=l(i)*6;
krs(ps4,ps2)=2*l(i)^2;
krs(ps4,ps3)=-l(i)*6;
krs(ps4,ps4)=4*l(i)^2;

Kij=Kij+(I*E/(l(i)^3))*krs;
krs=zeros(2*(k+1));
end

Kij; %Matriz de Rigidez Global

%RESOLUCIÓN DE LA ECUACIÓN DE RIGIDEZ

%Primero reescribimos el sistema para agrupar las variables en un solo
%vector

Kmod=Kij;
Kmod(1:2*(k+1),1)=0;
Kmod(1,1)=-1;
Kmod(1:2*(k+1),2*(k+1)-1)=0;
Kmod(2*(k+1)-1,2*(k+1)-1)=-1;

Kmod; %Matriz de rigidez modificada

%Se resuelve la ecuación para todas las variables

RESULTADOS=inv(Kmod)*fuerzas;

%HALLAMOS LOS GRADOS DE LIBERTAD

q=RESULTADOS;

Página | 16
q(1)=0;
q(2*(k+1)-1)=0;

q; %Vector que contiene los grados de libertad del sistema

%HALLAMOS EL VALOR DEL VECTOR FUERZA

fuerzas=RESULTADOS-q;

fuerzas; %Vector que contiene las fuerzas y momentos en la viga

%CÁLCULO DE ESFUERZOS MÁXIMOS

ESFUERZOS=zeros(k,1);
syms e;
for i=1:k
f(e)=(-E*ymax/l(i)^2)*(6*e*q(2*i-1)+(3*e-1)*l(i)*q(2*i)-
6*e*q(2*(i+1)-1)+(3*e+1)*l(i)*q(2*(i+1)));
if abs(f(1))>abs(f(-1))
ESFUERZOS(i)=f(1);
elseif abs(f(-1))>abs(f(1))
ESFUERZOS(i)=f(-1);
end
end

ESFUERZOS;

%GRÁFICA DE MOMENTO FLECTOR

x1=linspace(-1,7,1000);
inicio=0;
final=-1;
y1=0;
u=0;
for i=1:k
inicio=final;
final=final+2*l(i)/250;
a=2/(final-inicio);
b=1-final*a;
u=a.*x1+b;
y1=y1+((inicio<=x1)&(x1<final)).*(E*I/l(i)^2).*(6.*q(i*2-
1).*u+(3.*u-1).*l(i).*q(i*2)-6.*u.*q(2*(i+1)-
1)+(3.*u+1).*l(i).*q(2*(i+1)));
u=0;
end

figure(1)
plot(x1,y1)
title('Gráfica del momento flector en la viga');
grid on

%GRÁFICA DE FUERZA CORTANTE

x2=linspace(-1,7,1000);
inicio=0;

Página | 17
final=-1;
y2=0;
u=0;
for i=1:k
inicio=final;
final=final+2*l(i)/250;
a=2/(final-inicio);
b=1-final*a;
u=a.*x2+b;
y2=y2+((inicio<=x2)&(x2<final)).*(6*E*I/l(i)^3).*(2*q(2*i-
1)+l(i)*q(2*i)-2*q(2*(i+1)-1)+l(i)*q(2*(i+1)));
u=0;
end

figure(2)
plot(x2,y2)
title('Deformada de la Viga');
grid on

%GRÁFICA DE LA DEFORMADA

x3=linspace(-1,7,1000);
inicio=0;
final=-1;
y3=0;
u=0;
for i=1:k
inicio=final;
final=final+2*l(i)/250;
a=2/(final-inicio);
b=1-final*a;
u=a.*x3+b;
y3=y3+((inicio<=x3)&(x3<final)).*0.25.*(q(i*2-1)*(2-
3.*u+u.^3)+l(i).*0.5.*q(i*2).*(1-u-u.^2+u.^3)+q(2*(i+1)-1).*(2+3.*u-
u.^3)+l(i).*0.5.*q(2*(i+1))*(-1-u+u.^2+u.^3));
u=0;
end

figure(3);
plot(x3,y3)
title('Deformada de la Viga');
grid on

%RESULTADOS
disp('................................................................
.............')
disp(' RESULTADOS')
disp('----------------------------------------------------------------
-------------')

disp('Los grados de libertad del sistema en mm y rad:')

disp(q)

disp('La reacción en el apoyo simple 1 en N: ')

disp(fuerzas(1))

disp('La reacción en el apoyo simple 2 en N: ')

disp(fuerzas(2*(k+1)-1))

Página | 18
disp('La pendiente de la deformada en x=0: ')
disp(q(2))

disp('La pendiente de la deformada en x=L: ')

disp(q(2*(k+1)))

disp('El vector de esfuerzos máximos en MPa:')

disp(ESFUERZOS)

Página | 19
3.2 Vista en el “Command Window”

Página | 20
Gráficas

Página | 21
4. CONCLUSIONES

• Notamos que el valor de esfuerzo máximo se da en un extremo del elemento finito 3

y se da en la parte superior de la sección de la viga (el ala). Este esfuerzo máximo es
de compresión y tiene un valor igual a 3.0198 MPa.

• Las reacciones en los extremos de la viga resultaron ser iguales en magnitud (200 N)
pero opuestas en dirección. Estos valores cumplen con el equilibrio de traslación de
la viga.

• Los valores de reacciones y esfuerzos obtenidos con el método de los elementos

finitos son casi iguales a los valores obtenidos mediante la aplicación de las
ecuaciones de la resistencia de materiales. Esto indica que 4 elementos finitos son
suficientes para analizar esta viga y obtener resultados bastante precisos.

• De la gráfica de fuerza cortante, notamos que la fuerza cortante en la viga es constante

e igual a 200 N.

• A partir de la gráfica de momento flector obtenida, notamos que el momento flector

máximo se da en un extremo del elemento finito 3 y tiene un valor de 140000 N-mm.
Asimismo, cabe resaltar que el esfuerzo máximo se da en la misma sección en la cual
se da el momento flector máximo.

• Notamos que la forma de la deformada de la viga es cóncava hacia arriba y su

pendiente se acentúa más en la zona en la que actúa el momento externo de 200000
N-mm.

Página | 22