Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
1
UNIVERSIDAD NACIONAL DE
INGENIERIA
Facultad de Ingeniera
Mecnica
PRIMER LABORATORIO
Curso : Clculo por Elementos Finitos MC516
Seccin : C
Profesor : Ing. Ronald Cueva Pacheco
Tema : Traccin Simple
Alumno(s): Apellidos y Nombres Cdigo
CALLE FLORES, Ren Eduardo
20080073D
2013 - 1
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
2
INDICE
PG
ndice 2
Enunciado del problema 3
Solucin Analtica 4
Diagrama de flujo 9
Solucin usando MATLAB 10
Conclusiones 12
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
3
ENUNCIADO DEL PROBLEMA
Dada la siguiente placa triangular de espesor constante (t = 150 mm), se pide modelar
el cuerpo con 3 elementos finitos y calcular todas las variables desconocidas, as como
los esfuerzos.
DATOS:
()
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
4
SOLUCION ANALTICA
a) MODELADO DEL CUERPO REAL
Modelaremos el cuerpo real con tres elementos finitos, tal como se muestra en
la figura siguiente:
Del grfico, "b1", "b2" y "b3", son las longitudes transversales de cada
elemento finito. stas se calcularn en base al punto medio de cada elemento
finito y por medio de relaciones geomtricas sencillas.
( )
( )
Adems "A1", "A2" y "A3" son las reas respectivas de cada elemento finito.
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
5
Construyendo el cuadro de conectividad
e
NODOS GDL
()
)
(1) (2) 1 2
1 1 2 1 2 600
112500
2 2 3 2 3 300
56250
3 3 4 3 4 300
18750
b) Grados de libertad nodales (Vector Desplazamiento)
Los grados de libertad nodales se muestran en la siguiente figura:
Entonces el vector desplazamiento ser
]
Se asume que la placa esta empotrada, entonces Q1 = 0.
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
6
c) Vector Carga
Analizando las fuerzas globales en el cuerpo, tenemos:
()
()
()
()
()
()
()
()
()
()
El vector fuerza global ser:
] [
]
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
7
d) Matriz de Rigidez
La matriz de rigidez global ser
(
(
(
(
|
.
|
\
|
=
}
0 0 0 0
0 0 0 0
0 0 1 1
0 0 1 1
1
l
AE
K
i
(
(
(
(
|
.
|
\
|
+
0 0 0 0
0 1 1 0
0 1 1 0
0 0 0 0
2
l
AE
(
(
(
(
|
.
|
\
|
1 1 0 0
1 1 0 0
0 0 0 0
0 0 0 0
3
l
AE
Reemplazando los datos, se obtiene:
5
562.5 562.5 0 0
562.5 1125 562.5 0
10
0 562.5 750 187.5
0 0 187.5 187.5
i
N
K x
mm
}
(
(
(
=
(
(
e) Matriz de Rigidez
La matriz de rigidez global est determinada por la siguiente ecuacin:
Si reemplazamos nuestros datos en la ecuacin obtenemos:
2648.7 1
13310.875
882.9
220.725
R + (
(
(
(
(
2 5
3
4
0 562.5 562.5 0 0
562.5 1125 562.5 0
10
0 562.5 750 187.5
0 0 187.5 187.5
Q
x
Q
Q
( (
( (
( (
=
( (
( (
Para obtener los desplazamientos, se toma la siguiente submatriz:
11310.875
882.9
220.725
(
(
(
(
2
5
3
4
1125 562.5 0
10 562.5 750 187.5
0 187.5 187.5
Q
x Q
Q
( (
( (
=
( (
( (
Resolviendo, se obtiene:
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
8
Para la reaccin en el empotramiento, se tiene:
| | | |
2 5
1
3
4
0
2648.7 10 562.5 562.5 0 0
Q
R x
Q
Q
(
(
(
+ =
(
(
Resolviendo tenemos:
f) Esfuerzos
Para el clculo de los esfuerzos se usar la siguiente expresin:
[
] [
]
Y obtenemos lo siguiente:
)[
] [
)[
] [
)[
] [
g) Resultados
Finalmente, los resultados son los siguientes:
1 1 1
5
2 2 1
2
5
3 3 2
2
5
4 4 3
2
0 14415.51 17063.21
25.6258 10 13310.875 0.128129
27.5878 10 882.9 0.196200
28.7650 10 220.725 0.011772
Q mm F N R N
N
Q x mm F N
mm
N
Q x mm F N
mm
N
Q x mm F N
mm
o
o
o
= = =
= = =
= = =
= = =
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
9
DIAGRAMA DE FLUJO
Inicio
Modelado del Problema Se elige el vector X, los nodos y las partes de la
figura
Crear Tabla de Conectividad
Calculo de las Matrices de rigidez Local
Calculo de Matriz de rigidez Global
Creacin de la Matriz de carga
Obtencin de las Matrices reducidas de carga, desplazamiento y de
rigidez
Calculo de los Desplazamientos Nodales
Calculo de las cargas nodales y de las reacciones en los
apoyos
Calculo de esfuerzos en cada elemento finito
Mostrar resultados:
Desplazamientos,
Cargas, esfuerzos y
reaccin de apoyos
FIN
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
10
SOLUCION USANDO MATLAB
clc, clear all,close all;
%---------------------------------------------------------------------
----
% RESOLUCION DEL PROBLEMA 1ra practica (CEF)
% Tema: Traccin Simple
%---------------------------------------------------------------------
----
% Nombre : RENE EDUARDO EDUARDO CALLE FLORES
% Curso : CALCULO POR ELEMENTOS FINITOS - MC516
% Seccin: C
%---------------------------------------------------------------------
% 1. DATOS
%---------------------------------------------------------------------
% 1.1. DIMENSIONES
h1 = 1000; % mm (base)
L = 1200; % mm (alura)
t = 150; % mm (espesor)
% 1.2. DEL MATERIAL
E = 3e5; % N/mm2
% 1.3. CARGAS
Pa = 10000; % (N)
gamma = 8; % (gr-f/cm^3)
gamma = gamma*(9.81e-6);
% 1.4. ELEMENTOS FINITOS
%
% Elemento 1
L1 = 600; % mm
A1 = 750*t; % mm^2
% Elemento 2
L2 = 300; % mm
A2 = 375*t; % mm^2
% Elemento 3
L3 = 300; % mm
A3 = 125*t; % mm^2
%---------------------------------------------------------
% 2. CODIGO PRINCIPAL
%---------------------------------------------------------
% 2.1. Vector desplazamiento
Qj = zeros(4,1); % Qj = [Q1 Q2 Q3 Q4]'
% Q1 = 0
% 2.2. Vector carga
Fi = zeros(4,1); % Fi = [F1 F2 F3 F4]'
% -> Valores de la mitad del peso de cada elemento finito
g1 = gamma*A1*L1/2;
g2 = gamma*A2*L2/2;
g3 = gamma*A3*L3/2;
% -> Por el momento slo calcularemos F2,F3 y F4
Fi(2) = g1 + g2 + Pa;
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
11
Fi(3) = g2 + g3;
Fi(4) = g3;
% 2.3. Matriz de Rigidez Global
k1 = E*A1/L1*[1 -1 0 0;-1 1 0 0;0 0 0 0;0 0 0 0];
k2 = E*A2/L2*[0 0 0 0;0 1 -1 0;0 -1 1 0;0 0 0 0];
k3 = E*A3/L3*[0 0 0 0;0 0 0 0;0 0 1 -1;0 0 -1 1];
Kij = k1 + k2 + k3;
disp('Matriz de Rigidez')
disp(Kij)
% 2.4. Usando la ecuacin de rigidez
% OJO: Normalmente el numero de ecuaciones lineales seran 4, pero
% se observa que la primera de ellas (la que contiene a "R1")
% es independiente de las otras.
% -> se obtiene el siguiente sistema de ecuaciones reducido
% Fi_r = Kij_r*Qj_r
Fi_r = Fi(2:4);
Kij_r = Kij(2:4,2:4);
Qj_r = pinv(Kij_r)*Fi_r; % Aqui Qj_r es la variable
% -> Obtenemos los valores de Q2, Q3 Y Q4
Qj(2:4) = Qj_r;
% -> Calculamos la carga F1 y la reaccion en el apoyo R1
Fi(1) = Kij(1,:)*Qj; % Valor de F1
R1 = Fi(1)-g1; % Valor de R1
% 2.5. Calculo de los esfuerzos
sigma1 = (E/L1)*[-1 1]*[Qj(1) Qj(2)]';
sigma2 = (E/L2)*[-1 1]*[Qj(2) Qj(3)]';
sigma3 = (E/L3)*[-1 1]*[Qj(3) Qj(4)]';
sigma = [sigma1 sigma2 sigma3]';
%---------------------------------------------------------
% 3. PLOTEANDO RESULTADOS
%---------------------------------------------------------
disp('----------------------RESULTADOS----------------------------')
disp('1. Valor de la reaccin en el apoyo "R1" (en N)')
disp(R1)
disp('2. Vector de desplazamiento "Qj" (en mm)')
disp(Qj)
disp('3. Vector de carga "Fi" (en N)')
disp(Fi)
disp('4. Vector de Esfuerzos para cada E.F "sigma_e" (en N/mm^2)')
disp(sigma)
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
12
La salida de MATLAB es la siguiente:
Matriz de Rigidez
56250000 -56250000 0 0
-56250000 112500000 -56250000 0
0 -56250000 75000000 -18750000
0 0 -18750000 18750000
----------------------RESULTADOS----------------------------
1. Valor de la reaccin en el apoyo "R1" (en N)
-1.7063e+04
2. Vector de desplazamiento "Qj" (en mm)
1.0e-03 *
0
0.2563
0.2759
0.2876
3. Vector de carga "Fi" (en N)
1.0e+04 *
-1.4415
1.3311
0.0883
0.0221
4. Vector de Esfuerzos para cada E.F "sigma_e" (en N/mm^2)
0.1281
0.0196
0.0118
Clculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERA
MC516 - C Facultad de Ingeniera Mecnica
13
CONCLUSIONES
- Se observa que las deformaciones son completamente pequeas, esto es lgico, pues
al tratarse de un cuerpo rgido.
- Los esfuerzos calculados son positivos, esto se debe a que el punto 2 de cada elemento
se mueve ms que el punto 1, esto significa que el elemento est a traccin en la lnea
de accin de nuestra referencia "x".
- Se observa que los resultados de MATLAB arrojan un error prcticamente cero, lo cual
nos indica que la aproximacin del cuerpo a tres elementos finitos es prcticamente
exacta.
- Se aprecia tambin que el eje de referencia "x" considerado debe ser asumido en
primera instancia con experiencia visualizando como se deformara el modelo y as
tener menos problemas en la utilizacin del mtodo.
