Primera Practica MC516
Primera Practica MC516
Primera Practica MC516
PRIMERA PRÁCTICA
Cálculo por Elementos Finitos (MC516)
Tema:
Tracción Simple
Estudiante:
Robles Loro, José André 20162007E
Sección:
I
Profesor:
Cueva Pacheco, Ronald
2019-I
Lima, 25 de abril.
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
2
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
------------------------------------------------------------------------------------------------------- ---------------------
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.
3
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
------------------------------------------------------------------------------------------------------- ---------------------
4
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
------------------------------------------------------------------------------------------------------- ---------------------
𝐹𝑖 = [𝐹1 𝐹2 𝐹3 𝐹4 ]
Y la matriz de Rigidez está dada por.-
1 1 0 0 0 0 0 0 0 0 0 0
0 AE 0 1 1
0 AE 0
0 0
AE 1 1 0 0
Ki
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
rigidez Fi K i Q , por ello no es necesario seguir haciendo los calculos a
mano.
5
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
------------------------------------------------------------------------------------------------------- ---------------------
DIAGRAMA DE FLUJO
Inicio
Calculamos la matriz de
rigidez global (Kt)
Mostrar resultados:
Desplazamientos,
Cargas, esfuerzos y
reacción de apoyos
FIN
6
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=[];
7
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)
8
Cálculo por Elementos Finitos UNIVERSIDAD NACIONAL DE INGENIERÍA
MC516 - I Facultad de Ingeniería Mecánica
------------------------------------------------------------------------------------------------------- ---------------------
CONCLUSIONES
Se puede observar que las deformaciones son del orden de milesimas,
por ello la placa triangular al ser expuesta a la carga P dada se mantiene
como un cuerpo rígido.
Los desplazamientos y esfuerzos calculados son negativos, esto se debe
a que la carga P y el peso de la placa son antiparalelos con la dirección
de nuestra línea de acción x, lo que significa que la placa está en
compresión.
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.