Cálculo por elementos finitos MC-516
UNIVERSIDAD NACIONAL DEINGENIERIA
FACULTAD DE INGENIERIA MECANICA
CALCULO POR ELEMENTOS FINITOS
SECCION: B
INTEGRANTES: FLORES DONAYRE RONNY KINO
GOMEZ HUARCAYA DANIEL ENRIQUE
LIZARBE FABIAN LEONARDO ANDRES
PROFESOR: ING ABREGÚ LEANDRO
Páá giná 1
Cálculo por elementos finitos MC-516
Índice
Enunciado del Problema....................................................................3
Solución.............................................................................................4
Grados de Libertad Nodales..............................................................5
Vector Carga......................................................................................6
Matriz de Rigidez................................................................................8
Ecuación de Rigidez y Condición de Contorno..................................9
Esfuerzos y Resultados....................................................................10
Diagrama de Flujo.............................................................................11
Uso de Matlab...................................................................................13
Conclusiones................................................................................... 15
Bibliografía….................................................................................... 16
Páá giná 2
Cálculo por elementos finitos MC-516
PRIMERA PRÁCTICA CALIFICADA
(TRACCION SIMPLE)
PROBLEMA
En siguiente placa trapezoidal sometida a una carga en su punto
medio y otra en su base menor, calcular los esfuerzos en cada
elemento finito y la reacción en el apoyo.
Consideraciones:
P =100 lb
ρ = 0.2836 lb/pulg3
t = 1”
E = 30x106 psi
Páá giná 3
Cálculo por elementos finitos MC-516
SOLUCION:
1. MODELADO DEL CUERPO REAL
Se consideraran 4 elementos finitos. El ancho de cada elemento será
el punto medio del mismo:
b
6 5.25
5 .62 5 "
1
2
b2
5.25 4 .5
4 .875
2
b3
4 .5 3.75
4 .12 5
2
b4
3.75 3
3. 37 5 "
2
Luego obtenemos el siguiente modelado: Dividimos la figura por la
cantidad de nodos indicados para el problema.
Cuadro de conectividad:
Ai bi x t
t = 1”
NODOS GDL be
l e
Ae
Páá giná 4
Cálculo por elementos finitos MC-516
e 1 2 1 2 (”) (”) (”2)
1 1 2 1 2 5.625 1 5.625
2 2 3 2 3 4.875 1 4.875
3 3 4 3 4 4.125 1 4.125
4 4 5 4 5 3.375 1 3.375
2. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)
A través del grafico se muestran los grados de libertad nodales
globales:
Luego el vector de desplazamiento será:
Páá giná 5
Cálculo por elementos finitos MC-516
Q1
Q2
Q Q3 "
Q4
Q5
Donde Q1= 0 pues la placa esta empotrada y los demás
desplazamientos son incógnitas que habrán de ser calculadas.
3. VECTOR CARGA
Páá giná 6
Cálculo por elementos finitos MC-516
Analizando las fuerzas en cada elemento finito:
Formaremos la matriz para las fuerzas respectivas
(A1l1) R1
2
(A1l1) (A2l2 ) 0
2 2
(A2l2 ) (A3l3)
F P lb
2 2
(A l ) (A l )
3 3 4 4 0
2 2
(A4l4 )
0
2
Páá giná 7
Cálculo por elementos finitos MC-516
4.78 R1
4.78 4.14 0
F 4.143.5 10 lb
3.5 2.87 0
2.87 0
Entonces, el vector carga se expresaría de la siguiente manera
F1 4.78R1
F2 8.92
F1 F3107.64 lb
F4 6.37
F5 2.87
Páá giná 8
Cálculo por elementos finitos MC-516
4. MATRIZ DE RIGIDEZ
A continuación pasamos a calcular la matriz de Rigidez Global, que
está determinada por la siguiente ecuación:
1 1 0 0 0 0 0 0 0 0
1 1 0 0
0 0 1 1 0 0
AE AE
K iJ 0 0 0 0 0 0 1 1 0 0
l 1 l 2
0 0 0 0 0 0 0 0 0 0
0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0
0 0 0 0 0 0
AE AE
0 0 1 1 0 0 0 0 0 0
l 3 l 4
0 0 1 1 0 0 0 0 1 1
0 0 0 0
0 0 0 0 1 1
Reemplazando para los valores calculados y utilizando la tabla de
conectividad obtenemos:
1 1 0 0 0 0 0 0 0 0
1 1 0 0 0 0 1 1 0 0
5.625 xE 4.875 xE
K iJ 0 0 0 0 0 0 1 1 0 0
6 1 6 2
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
Páá giná 9
Cálculo por elementos finitos MC-516
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
4.125 xE 3.375 xE
0 0 1 1 0 0 0 0 0 0
6 3 6 4
0 0 1 1 0 0 0 0 1 1
0 0 0 0 0 0 0 0 1 1
Luego:
0.9375 0.9375 0 0 0
0.9375 1.75 0.8125 0 0
lb
K iJ E x 0 0.8125 1.5 0.6875 0
in
0 0 0.6875 1.25 0.5625
0 0 0 0.5625 0.5625
Donde: E 30 x10 6 psi
5. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO
La ecuación de rigidez está determinada por la siguiente ecuación:
Fi K iJ Q J …………………(1)
Sabemos que:
Páá giná 10
Cálculo por elementos finitos MC-516
F1 4.78R1 F: fuerzas
Q: cargas
F2 8.92
K: rigidez
F1 F3107.64 lb
F4 6.37
F5 2.87
0.9375 0.9375 0 0 0
0.9375 1.75 0.8125 0 0
lb
K iJ E x 0 0.8125 1.5 0.6875 0
pu lg
0 0 0.6875 1.25 0.5625
0 0 0 0.5625 0.5625
Reemplazando:
Páá giná 11
Cálculo por elementos finitos MC-516
4.78 R1
8.92
107.64
6 .37
2.87
09375 0.9375 0 0 0 0
0.9375 1.75 0.8125 0 0 Q
2
(30 x10 6 ) 0 0.8125 1.5 0.6875 0 Q3
0 0 0.6875 1.25 0.5625 Q 4
0 0 0 0.5625 0.5625 Q5
Resolviendo este sistema de ecuaciones obtenemos:
R1 130.57975 lb
6
Q2 4.47288 x10 "
6
Q3 9.26796 x10 "
6
Q4 9.71596 x10 "
6
Q5 9.88603 x10 "
6. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la
siguiente ecuación:
e
E Qi
e 1 1 ……………….. (2)
l Qi 1
Y obtenemos lo siguiente:
Páá giná 12
Cálculo por elementos finitos MC-516
30x10 0 6 lb
6
1 11 x10 1 2.2 364 2
6 1 4. 728 pulg
30x10 4. 728 6 lb
6
2 11 x10 2 3.2 9754 2
6 2 9.26796 pulg
30x10 9.26796 6 lb
6
3 11 x10 3 2. 40 2
6 3 9.71596 pulg
30x10 9.72596 6 lb
6
4 11 x10 4 0.80 35 2
6 4 9.8 603 pulg
7. RESULTADOS
Finalmente, los resultados son mostrados en la siguiente tabla:
Páá giná 13
Cálculo por elementos finitos MC-516
R1 130.57975 lb
lb
1 22.3644
pu lg 2
lb
2 23.9754
pu lg 2
lb
3 2.2400
pu lg 2
lb
4 0.80035
pu lg 2
8. DIAGRAMA DE FLUJO
INICIO
INGRESO DE DATOS
CONSTANTES : E, ρ, t
VECTORES : L, A, P
CÁLCULO DE VECTORES
( A1l1 )
2 1 R EA1 EA1
L1 1 0 0 0
( A l
11 2 2 ) ( A l ) L
EA1 EA EA1
2
EA2
2 2 2 1 2 0 0
L1 L L L
( A l ) ( A l ) 2
EA EA2
3
EA3
F P1 ; K 0
EA
2 2
3 3
2 3 2 0
2 2 L L L L3
(A l ) (A l ) 0 0
EA 3
EA4 EA3
4 3
EA4
33 44 L3 L L L4
2 2
EA4 EA5
0
( A4l4 )
0 0
L4 L5
2
TRAFORMACION DE ECUACION MATRICIAL
Páá giná 14
Cálculo por elementos finitos MC-516
( A1l1 )
2 R1 EA1 EA1
L1
L1
0 0 0
( A l )
1 1 2 2 EA1( A l )
EA 2 EA1 EA 2
2 2 L1 2 1
L2
0 0
0
L L Q
( A l ) ( A
2 2 3 3 P1 = 0l ) EA 2
EA3 EA 2 EA3 2
2 3 2 0 Q3
2 2 L L L L3
Q
4
(A l ) (A l ) 0 0
EA 3
EA 4 EA3
4 3
EA 4 Q5
33 44 L3 L L L4
2 2 EA 4
EA5
( A4l4 ) 0 0 0
L4 L5
2
IMPRESIÓN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , 1 , 2 , 3
FIN
Código fuente en MATLAB:
%PRIMERA PRACTICA CALIFICADA
%CALCULO DE ELEMENTOS FINITOS
%MC-516 "A"
%NUMERO DE NODOS: "5 nodos"
clc
disp(' ***************');
disp(' TRACCION SIMPLE');
disp(' ***************');
disp(' ');
H= input('Ingrese la altura de la placa [in]= ');
B1=input('ingrese la base mayor [in]= ');
B2=input('ingrese la base menor [in]= ');
t= input('Ingrese el espesor de la placa[in]= ');
disp(' ');
p1=input('Ingrese carga P [lb]= ');
disp(' ');
y=input('Ingrese el peso especifico [lb/in3]= ');
E=input('Ingrese el modulo de Young [ psi]= ');
f=y;
h=H/4;
s=0;
w=zeros(5);K55=zeros(5);
Páá giná 15
Cálculo por elementos finitos MC-516
c=(B1-B2)/4;
for i=1:4
b(i)=B1-(2*i-1)*c/2;
A(i)=t*b(i);
end
b
A
F=[];
F(1)=f*A(1)*h/2
F(2)=f*A(1)*h/2+f*A(2)*h/2
F(3)=f*A(2)*h/2+f*A(3)*h/2+p1
F(4)=f*A(3)*h/2+f*A(4)*h/2
F(5)=f*A(4)*h/2
F=F';
for i=1:4
w(i,i)=1;
w(i,i+1)=-1;
w(i+1,i)=-1;
w(i+1,i+1)=1;
K55=K55+A(i)*E/h*w;
w=zeros(5);
end
k55=K55(2:5,2:5);
Q=k55^-1*F(2:5,1);
Q=[0;Q];
R=K55(1,1:5)*Q-F(1);
es=[];
for i=1:4
es(i,1)=E/h*[-1 1]*Q(i:i+1,1);
end
%MOSTRANDO LOS RESULTADOS
disp(' ');
disp(' **********');
disp(' RESULTADOS');
disp(' **********');
disp(' ');
disp('...........................................................');
disp('VECTOR FUERZA [F]');
disp(F);
disp('...........................................................');
disp('VECTOR DESPLAZAMIENTO [Q]');
disp(Q);
disp('...........................................................');
disp('MATRIZ DE RIGIDEZ');
disp(K55);
disp('...........................................................');
disp('REACCION EN EL EMPOTRADO [R]');
disp(R);
disp('...........................................................');
disp('EL VECTOR DE ESFUERZOS(N/mm2) ');
disp(' e1 e2 e3 e4');
disp(es');
Páá giná 16
Cálculo por elementos finitos MC-516
9. USO DEL PROGRAMA DE MATLAB
Entrada de datos:
Ingrese el espesor [in]= 1
**********************
TRACCION SIMPLE Ingrese carga P [lb]= 100
**********************
Ingrese peso específico [lb/in3]= 0.2836
Ingrese la altura [in]= 24
ingrese la base mayor [in]= 6 Ingrese el M de Young [ psi]= 30’000,000
ingrese la base menor [in]= 3
Salida de datos:
**********
RESULTADOS
**********
...........................................................
VECTOR FUERZA [F]
4.7858
8.9334
107.6572
6.3810
2.8715
...........................................................
VECTOR DESPLAZAMIENTO [Q]
1.0e-05 *
0
0.4474
0.9271
0.9719
0.9889
Páá giná 17
Cálculo por elementos finitos MC-516
...........................................................
MATRIZ DE RIGIDEZ
28125000 -28125000 0 0 0
-28125000 52500000 -24375000 0 0
0 -24375000 45000000 -20625000 0
0 0 -20625000 37500000 -16875000
0 0 0 -16875000 16875000
...........................................................
REACCION EN EL EMPOTRADO [R]
-130.6288
...........................................................
EL VECTOR DE ESFUERZOS(N/mm2)
e1 e2 e3 e4
22.3721 23.9815 2.2430 0.8508
Páá giná 18
Cálculo por elementos finitos MC-516
CONCLUSIONES
Las deformaciones obtenidas obtienen valores por el orden de las
micras (muy bajos), además se obtiene todos los sentidos hacia
abajo como era de esperarse.
Los esfuerzos son positivos, lo que indica esfuerzos de tracción
para nuestro sistema de referencia.
Se puede apreciar que los resultados de MATLAB arrojan un error
de 1.846x10-14, valor que en la práctica equivale a cero, por lo
tanto se considera como exacto.
Para geometrías simples como esta figura, la solución se obtiene
de manera sencilla, pero, a medida que se va analizando cuerpos
de geometría compleja, el análisis se hace también más complejo;
es allí donde se demuestra el poder del análisis por elementos
finitos.
Páá giná 19
Cálculo por elementos finitos MC-516
BIBLIOGRAFIA
CHANDRUPATLA, Introducción al estudio de elementos finitos en
ingeniería, Tirupathi R. Chandrupatla
Apuntes en clase
Páá giná 20
