Universidad Nacional de Ingeniería
Facultad de Ingeniería Mecánica
PRIMERA PRÁCTICA
Cálculo por Elementos Finitos (MC516)
Tema:
Tracción Simple
Estudiante:
SALLO AMAU DENNIS FRANK 20180229F
Sección: E
Profesor:
Cueva Pacheco, Ronald
2021-II
Lima, 23 de septiembre.
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
ÍNDICE
ENUNCIADO DEL PROBLEMA ........................................................ 3
MODELADO DEL CUERPO REAL ................................................... 3
DIAGRAMA DE FLUJO ..................................................................... 6
CÓDIGO MATLAB ............................................................................. 7
CONCLUSIONES .............................................................................. 9
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
ENUNCIADO DEL PROBLEMA
Dada una placa triangular (Fig. 1.) de acero con
espesor constante t=150mm y las siguientes
caracteristicas.-
𝐸 = 3.0𝑥105𝑁/𝑚𝑚2
𝛾 = 8.0 𝑔𝑟 − 𝑓/𝑐𝑚3
Para una carga P=50kN, calcular los esfuerzos
para 3 elementos finitos.
Fig. 1
MODELADO DEL CUERPO REAL
Para modelar el cuerpo con tres elementos finitos, se selecciono las siguientes
alturas para cada elemento.-
𝑙1 = 750𝑚𝑚
𝑙2 = 375𝑚𝑚
𝑙3 = 375𝑚𝑚
Y el ancho de cada elemento se calculó teniendo en cuenta que el volumen de
la placa triangular tiene que ser igual a la suma del volumen de los tres elementos
finitos.
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
Para ello, calculamos el area promedio de cada base e igualamos al area
correspondiente a cada elemento.-
(1000+500)
𝑏 ∗𝑡 = ∗ 𝑡 𝑏 = 750𝑚𝑚
1 2 1
(500+250)
𝑏 ∗𝑡 = ∗ 𝑡 𝑏 = 375𝑚𝑚
2 2 2
(250+0)
𝑏 ∗𝑡 = ∗ 𝑡 𝑏 = 125𝑚𝑚
3 2 3
Con los datos de las bases de los elementos y las alturas, el modelo del triangulo
con 3 elementos finitos está dado por.-
Así, tabla de conectividad es.-
e Nodos GDL 𝑙𝑒(mm) 𝐴𝑒(𝑚𝑚2)
(1) (2) (1) (2)
1 1 2 1 2 750 1125
2 2 3 2 3 375 5625
3 3 4 3 4 375 3125
Así mismo, de acuerdo al DCL del modelo calcularemos el vector de carga
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
Global.-
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
Las fuerzas del vector de carga global son.-
𝛾(𝐴 ∗ 𝑙)1
𝐹1 = − +𝑅1
2
𝛾(𝐴 ∗ 𝑙)1 𝛾(𝐴 ∗ 𝑙)2
𝐹2 = − − −𝑃
2 2
𝛾(𝐴 ∗ 𝑙)2 𝛾(𝐴 ∗ 𝑙)3
𝐹3 = − −
2 2
𝛾(𝐴 ∗ 𝑙)3
𝐹4 = −
2
Donde el vector de carga es.-
𝐹𝑖 = [𝐹1 𝐹2 𝐹3 𝐹4]
Y la matriz de Rigidez está dada por.-
1 1 0 0
0 0 0 0 0 0 0 0
AE 1 1 0 0 AE 0 1 1 0 AE 0 0 0 0
K
i
l 1 0 0 0 0 l 2 0 1 1 0 l 3 0 0 1 1
0 0 0 0 0 0 0 0 0 0 1 1
Con los datos obtenidos hasta ahora del modelado es mas que suficiente para
hacer los calculos de los desplazamientos y la reacción R1 con un codigo
implementado en el programa MATLAB, en el cual resolveremos la ecuación de
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
rigidez Fi Ki Q , por ello no es necesario seguir haciendo los calculos a
mano.
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
DIAGRAMA DE FLUJO
Inicio
Introducimos geometria de cada elemento, datos del material y
carga aplicada y condiciones de contorno
Calculamos la matriz de
rigidez global (Kt)
Creación y calculo del vector de
Fuerza 3x3 con F2 F3 y F4 (F3x3)
Creación de matriz de rigidez
de 3x3 (K3x3)
Calculo los desplazamientos globales Q2 Q3 Q4 con
F3x3 y K3x3
Creación el vector de desplazamientos
global completo
Calculo el vector de carga completo
Calculo de R1 usando el primer termino del vector F
Calculo de esfuerzos en cada elemento finito
Mostrar resultados:
Desplazamientos,
Cargas, esfuerzos y
reacción de apoyos
FIN
NOTA: Debido a que R1 es una variable no se puede aplicar rapidamente la
ecuación de rigidez, es por ello que se procede a calcular Q2 Q3 Y Q4 primero
con un vector de carga y matriz de rigidez 3x3.
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
CÓDIGO MATLAB
clc; clear all; close all;
%
%Resolución de la Primera Practica MC516
%Tema: Tracción Simple
%
%1. DATOS
%Introducir datos de los elementos
%Elemento 1
l1=750; %mm
A1=750*150; %mm^2
%Elemento 2
l2=375; %mm
A2=375*150; %mm^2
%Elemento 3
l3=375; %mm
A3=125*150; %mm^2
%Datos del material
la=8*9.81e-6; %N/mm^3
E=3e5; %N/mm^2
%Carga
P=50000; %N
%Condiciones de contorno
Q1=0;
%
%2. CALCULOS
%Crear matriz de Rigidez
c1=E*A1/l1;
c2=E*A2/l2;
c3=E*A3/l3;
Kt=[c1 -c1 0 0
c1 c1+c2 -c2 0
0 -c2 c2+c3 -c3
0 0 -c3 c3];
%Mitad de peso de los elementos
g1 = la*A1*l1/2;
g2 = la*A2*l2/2;
g3 = la*A3*l3/2;
%
%Calculos de Q2 Q3 Q4
%Crear Vector de Fuerzas 3x3 con F2 F3 F4
F33=[-g1-g2-P -g2-g3 -g3];
%Matriz de rigidez 3x3
K33=Kt(2:4,2:4);
%Hallamos los desplazamientos Globales Q2 Q3 Q4
V=K33^-1*F33';
%
%Calculo de R1
%Vector desplazamientos global completo
Q=[Q1
V];
%Hallamos el vector de carga
%Aplicando la ecuacion de rigidez
F=Kt*Q;
%Con lo anterior hallamos R1
R1=F(1)+g1;
%
%Calculamos esfuerzos
sigma=[];
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
sigma(1)=(E/l1)*(-Q(1)+Q(2));
sigma(2)=(E/l2)*(-Q(2)+Q(3));
sigma(3)=(E/l3)*(-Q(3)+Q(4));
%
disp(' RESULTADOS ')
disp('1. Valor de la reacción en el apoyo "R1" (en N)')
disp(R1)
disp('2. Vector de desplazamiento "Qj" (en mm)')
disp(Q')
disp('3. Vector de carga "Fi" (en N)')
disp(F')
disp('4. Vector de Esfuerzos para cada E.F "sigma_e" (en N/mm^2)')
disp(sigma)
Los resultados obtenidos gracias a este código son los siguientes.-
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
----------------------------------------------------------------------------------------------------------------------------
CONCLUSIONES
La precisión en el resultado con respecto al uso del Matlab es muy
alta, lo que significa que para el cuerpo estudiado el número de
elementos finitos (tres) es suficiente gracias a su geometría
simple.
Para otras figuras, la precisión será directamente proporcional al
número de elementos finitos en que se divida, pues entre más se
escojan, menor error en los cálculos.
Para sistemas empotrados en tracción simple, el procedimiento
será el mismo que el usado anteriormente.
Con lo anterior dicho, se observa que MATLAB nos ahorra
bastante tiempo realizando calculos, debido a que el mismo
código se puede aplicar a otro sistema de tres elementos por lo
cual es muy útil.
BIBLIOGRAFÍA
CHANDRUPATLA, T. “Introducción al Estudio de los Elementos Finitos en
Ingeniería”, Prentice Hall, 1999 ZIENKIEWCTZ, O.
“The Finite Element Method”, New Cord, Mec Graw – Hill, 1977.
ZIENKIEWCTZ, O. and MORGAN K
