“Año de la universalización de la salud”
UNIVERSIDAD NACIONAL DE INGENIERIA
FACULTAD INGENIERIA MECÁNICA
1era Práctica Calificada
CURSO: CÁLCULO POR ELEMENTOS FINITOS
PROFESOR: CUEVA PACHECO Ronald
ALUMNO: AVILA TARRILLO Marlon
CÓDIGO: 20161188F SECCIÓN: “E”
LIMA
Junio 2020
Índice
Enunciado del Problema............................................................................. 3
Solución...................................................................................................... 4
Grados de Libertad Nodales....................................................................... 5
Vector Carga............................................................................................... 6
Matriz de Rigidez........................................................................................ 7
Ecuación de Rigidez y Condición de Contorno........................................... 8
Esfuerzos y Resultados.............................................................................. 9
Diagrama de Flujo....................................................................................... 10
Uso de Matlab............................................................................................. 11
Conclusiones……………………………………………………………………. 14
2
PRIMERA PRÁCTICA CALIFICADA
ENUNCIADO DEL PROBLEMA
En la figura se muestra una placa triangular simétrica; con espesor constante igual a 120
mm; está sometida a su peso propio y a una fuerza concentrada, tal como se indica.
HALLAR:
La distribución de esfuerzos a lo largo de la barra; modelándola con cuatro (04)
elementos finitos (unidimensionales).
La fuerza de reacción en el apoyo.
DATOS:
3
SOLUCIÓN:
1. MODELADO DEL CUERPO REAL
Se considerarán cuatro elementos finitos. Para facilitar los cálculos los elementos
finitos tendrán longitud de 750, 375, 250 y 125 mm.
Y los espesores lo calculamos tomando el punto medio de cada elemento finito:
( 1000+500 )
b1 = =750 mm
2
( 500+250 )
b2 = =375 mm
2
(250+83. 333 )
b3 = =166 . 666 mm
2
83 . 333
b4 = =41 .666 mm
2
Entonces, el modelado del cuerpo sería el siguiente:
Y las áreas se calculan de la siguiente relación:
A 1 = b1 x t
4
Cuadro de conectividad:
NODOS GDL le Ae
e (mm) (mm2)
(1) (2) 1 2
1 1 2 1 2 750 90000
2 2 3 2 3 375 45000
3 3 4 3 4 250 20000
4 4 5 4 5 125 5000
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á:
Q [0¿][Q2¿][Q3¿][Q4¿]¿ ¿
∫¿=¿ ¿ ¿
Donde Q1= 0 pues la placa esta empotrada y los demás desplazamientos son
incógnitas que tendrán que ser calculadas.
3. VECTOR CARGA
5
Analizando las fuerzas en cada elemento finito:
y ( Axl )1
F11 = +R1 = 2647 . 7+ R1 N
2
y ( Axl )1
F12 = +P A = 52647 .7 N
2
y ( Axl )2
F22 = = 661. 9 N
2
y ( Axl )2
F23 = = 661. 9 N
2
y ( Axl )3
F33 = = 196 .125 N
2
y ( Axl )3
F34 = = 196 . 125 N
2
y ( Axl )4
F 44 = = 24 . 51 N
2
y ( Axl )4
F54 = = 24 . 51 N
2
Ahora analizamos las fuerzas para todo el cuerpo:
6
F1 = F 11 = 2647 . 7 + R1 N
F2 = F 12 + F22 = 53309 . 6 N
F3 = F 23 + F33 = 858 .025 N
F 4 =F34 + F 44 = 220 .635 N
F5 =F 45 = 24 .51 N
Entonces, el vector carga se expresaría de la siguiente manera
F =¿ [ F 1¿ ] [ F 2¿ ] [ F 3¿ ] [ F 4¿ ] ¿ ¿
1
¿
4. MATRIZ DE RIGIDEZ
A continuación pasamos a calcular la matriz de Rigidez Global, que está determinada
por la siguiente ecuación:
K ¿ 0 0 0 0 0
[ ] [ ]
1 −1 0 0 0
0 1 −1 0 0
−1 1 0 0 0 AE
i ¿=
∫
AE
l 1( )
0 0 0 0 0
0 0 0 0 0
+ ( )
l 2
0 −1 1 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
( )
l 3
[ 0
0
0
0
0 0 0 0
0 1 −1 0
0 −1 1 0
0 0 0 0
] [ +
AE
( )
l 4
0
0
0
0
0
0
0
0
0 0 0
0 0 0
0 1 −1
0 −1 1
]
7
Reemplazando para los valores calculados y utilizando la tabla de conectividad
obtenemos finalmente:
[K 2.40 -2.40 0 0 0¿][ -2.40 4.80 -2.40 0 0¿][ 0 -2.40 4.0 0 -1.60 0¿][ 0 0 -1.60 2.40 -0.80 ¿] ¿
i∫¿=105 x ¿ ¿
5. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO
La ecuación de rigidez está determinada por la siguiente ecuación:
Fi = K i∫ Q∫ ¿ ¿
Lo que con nuestros valores calculados tenemos:
[¿2647. +R1¿][53 09.6¿][85 .025¿][2 0.635¿] ¿ =107x¿[ 2.40 -2.40 0 0 0¿][ -2.40 4.80 -2.40 0 0¿][ 0 -2.40 4.0 0 -1.60 0¿][ 0 0 -1.60 2.40 -0.80 ¿] ¿ ¿[0¿][Q2¿][Q3¿][Q4¿] ¿
¿ ¿ ¿
Para obtener los desplazamientos tomamos la siguiente submatriz:
[¿ 53 09.6¿][858.025¿][2 0.635¿] ¿ =107 x¿[ 4.80 0 -2.40 0 0 0¿][ -2.40 0 4.0 0 -1.60 0 0¿][ 0 -1.60 0 2.40 0 -0.80 0¿] ¿ ¿[Q2¿][Q3¿][Q4¿] ¿
¿ ¿ ¿
Resolviendo este sistema de ecuaciones obtenemos:
Q2 = 22.671 x10−4 mm
Q3 =23 .131 x10−4 mm
Q4 =23 . 284 x10−4 mm
Q5 =23 .315 x10−4 mm
Resolviendo obtenemos:
R1 = −57060. 5 N
8
6. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la siguiente
ecuación:
e
E
e
σ = [−1 1 ] ¿ [ Q i ¿ ] ¿ ¿
()
l ¿
Y obtenemos lo siguiente:
2 x 105
(
σ1 =
750 ) [−1 1 ] ¿ [ 0¿ ¿ ] ¿ ¿
1
5
2 x 10
σ =(
2 ) [−1 1 ] ¿ [ 22 . 671 ¿ ] ¿ ¿
375 2 ¿
5
2 x 10
σ =(
3 ) [ −1 1 ] ¿ [ 23. 131 ¿ ] ¿ ¿
250 3 ¿
5
2 x 10
σ =(
4 ) [−1 1 ] ¿ [ 23. 284 ¿]¿¿
125 4 ¿
7. RESULTADOS
Finalmente, los resultados son mostrados en la siguiente tabla:
R1 = −57060 .5 N
N
σ 1 = 0 . 6046
mm2
N
σ 2 = 0 . 0245
mm 2
N
σ 3 = 0. 0123
mm2
N
σ 4 = 0 . 0049
mm2
9
8. DIAGRAMA DE FLUJO
INICIO
INGRESO DE DATOS
CONSTANTES : E, f, t
VECTORES : L, A, P
CALCULO DE VECTORES
1
AL γ EA 1 EA 1
+ R1
[ ] [
0 0 0
]
2 −
L1 L1
AL2 γ AL1 γ EA 1 EA 2
EA 1 EA 2
+ +PA − + − 0 0
2 2 L1 L 2
L1 L2
3 2
AL γ AL γ EA 2 EA 2 EA 3 EA 3
+ 0 − + − 0
2 2 L2 L2 L3 L3
AL4 γ AL3 γ EA 3 EA 3
EA 4 EA 4
+ 0 0 − 3 + −
2 2 L L 3
L4 L4
AL4 γ EA 4 EA 4
0 0 0 −
2 L4 L4
F= ; K=
TRAFORMACION DE ECUACION MATRICIAL
1
AL γ EA
1
[ ][
0 0 0
][]
2 −1 −
L1
AL 2 γ AL1 γ EA 2
EA 1 EA 2
+ +PA 0 + − 0 0
2 2 L2 L1 L2 R1
AL 3 γ AL2 γ EA 2 EA 2
EA 3 EA 3
2
+
2
0 − + − 0 Q2
L2 L 2
L3 L3
AL γ AL3 γ
4
+ 0 0 −
EA 3 EA 3 EA 4
+ −
EA 4 Q3
2 2 L3 L3 L4 L4
Q4
AL4 γ EA 4 EA 4
0 0 0 −
2 L4 L4 Q5
=
IMPRESIÓN DE RESULTADOS
R 1 , Q 2 , Q 3 , Q 4 , E1 , E 2 , E 3
10
FIN
9. USO DEL PROGRAMA DE MATLAB
SCRIPT
clc
clear all
R1=sym('R1');
%datos de entrada
b0=1000 %input('Ingrese base superior(mm):')
bn=0 %input('Ingrese base inferior(mm):')
t=120 %input('Ingrese espesor(mm):')
h=1500 %input('Ingrese altura(mm):')
n=4 %input('Ingrese numero de elementos finitos:')
E=200000 %input('Ingrese modulo de elasticidad(N/mm2):')
y=0.00007845 %input('Ingrese densidad(N/mm3):')
Pa=50000 %input('Ingrese carga(N):')
%calculo de bases y áreas de elementos
le=zeros(n,1); ho=zeros(n,1); bo=zeros(n,1); b=zeros(n,1); a=zeros(n,1);
Fe=zeros(n+1,1);
bo(1)=b0; ho(1)=h;
for i=1:n
if n>i
le(i)=input('Ingrese longitud del elemento finito(mm):');
b(i)=(bo(i)+bn+(bo(i)-bn)*(ho(i)-le(i))/ho(i))/2;
a(i)=b(i)*t;
ho(i+1)=ho(i)-le(i);
bo(i+1)=2*b(i)-bo(i);
else
le(i)=ho(i);
b(i)=(bn+bo(i))/2;
a(i)=b(i)*t;
end
end
disp('Bases(mm):')
disp(b')
disp('Longitudes(mm):')
disp(le')
disp('Areas(mm^2):')
disp(a')
%calculo de las fuerzas
for i=1:n
Fe(i)=y*a(i)*le(i)/2;
end
for i=1:n+1
if i==1
F(i)=Fe(i);
elseif i==n+1
F(i)=Fe(i-1);
else
F(i)=Fe(i-1)+Fe(i);
end
11
end
F(2)=F(2)+Pa;
disp('El vector de fuerzas(N):')
disp(F')
%calculo de la matriz rigidez
k=zeros(n+1);
for i=1:n
x=zeros(n+1);
x(i,i)=1;x(i+1,i)=-1;x(i,i+1)=-1;x(i+1,i+1)=1;
k=k+(a(i)*E/(le(i)))*x;
end
disp('La matriz de rigidez es(N/mm):')
disp(k)
%calculo de desplazamientos
inv(k(2:n+1,2:n+1));
((F(2:n+1))');
Q=inv(k(2:n+1,2:n+1))*((F(2:n+1))');
Q=[0;Q];
disp('Los desplazamientos de los nodos son(mm):')
disp(Q)
%calculo de la reaccion
k(1,:)*Q;
R1=k(1,:)*Q-F(1);
disp('La reaccion en el extremo es:')
disp(R1)
%calculo de esfuerzos
for i=1:n
e(i)=(E/(le(i)))*[-1 1]*[Q(i); Q(i+1)];
end
disp('Los valores de los esfuerzos son(N/mm^2):')
disp(e');
12
VISTA EN EL COMMAND WINDOW DE MATLAB
13
CONCLUSIONES
Podemos apreciar, al utilizar más nodos, que las respuestas no varían
enormemente, solo aumentan la precisión con la cual se presentan.
Se recomienda utilizar un número moderado de nodos, ya que las operaciones con
matrices se vuelven demasiado engorrosas al ser de orden n x n donde n es el
número de nodos.
Se puede apreciar que las deformaciones son realmente pequeñas (décimas de
micras), además todas son hacia abajo que es el sentido positivo asumido como
referencia.
En el ejemplo no desarrollamos todo estrictamente en el SI, nos referimos
específicamente al uso de los metros, debido a las magnitudes de las elongaciones
y esfuerzos; es por ello que se utilizó en el desarrollo milímetros en vez de metros.
Los esfuerzos son positivos, lo que indica esfuerzos de compresión para nuestro
sistema de referencia.
14