Lab 2 Cálculo Por Elementos Finitos

UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería Mecánica

FACULTAD DE INGENIERÍA MECÁNICA
2DA PRÁCTICA CALIFICADA

CURSO: CÁLCULO POR ELEMENTOS FINITOS
PROFESOR: CUEVA PACHECO RONALD
AUTOR: LANDA CHÁVEZ GONZALO
CODIGO: 20161083J
SECCION: F
CICLO ACADÉMICO: 2021-1
1
Cálculo por Elementos Finitos Informe N°2
OBJETIVOS
 Calcular la fuerza de reacción que posee dicha placa triangular en el apoyo debido a
la carga aplicada y a la variación de temperatura.
 Desarrollar un código de programación con el fin de que el análisis se realice para N
elementos finitos.
 Calcular la distribución de esfuerzos y desplazamientos en las determinadas secciones
para dicha placa triangular de sección variable.
 Familiarizarse con el procedimiento del cálculo de los esfuerzos y deformaciones
mediante los elementos finitos en el cálculo de una tracción simple y el efecto de la
temperatura.
ENUNCIADO DEL PROBLEMA
Calcular la distribución de esfuerzos a los que está sometido la placa triangular y

la reacción que posee en el apoyo.
Datos:
 Placa triangular de espesor constante: t= 120 mm
 Material: E = 2 * 105 N/mm2 y 8 gr-f/cm3
 Efecto de temperatura: T= 125 °C.
2
DIGITACIÓN DEL PROBLEMA, CÁLCULOS Y RESULTADOS
La digitación de problema se realizó con el software Matlab, realizándose ahí los cálculos y
obteniendo los resultados.
Tenemos el enunciado del problema y el script realizado en Matlab:
3
Realizando el Modelado del Cuerpo Real

TABLA DE CONECTIVIDAD
nodos GDL
e longitude(mm) Áreae(mm2) E (N/mm2) mm3)
1 2 1 2
1 1 2 1 2 1500/N (1500/N) (L1) 200000 8
2 2 3 2 3 1500/N (1500/N) (L2) 200000 8
3 3 4 3 4 1500/N (1500/N) (L3) 200000 8
4 4 5 4 5 1500/N (1500/N) (L4) 200000 8
… … … … … … … 200000 8
N N N+1 N N+1 1500/N (1500/N) (LN) 200000 8
Lo siguiente es el código del script de Matlab:
clc
clear all
close all
%Datos del problema

H1=750;H2=750; %Altura total, se dividió en 2 segmentos para facilitar
cálculo
B=1000;pa=30000;t=120;j=8;E=200000;c=11*10^(-6);
dt=125; %variación de temperatura
%Ingreso número de elementos finitos

n1=input('Ingrese #elementos finitos parte superior: ') %#elementos parte
superior
n2=input('Ingrese #elementos finitos parte inferior: ') %#elementos parte
inferior
h1=ones(1,n1)*500/n1;
h2=ones(1,n2)*500/n2;
h=[h1 h2]; %matriz de la dstancia de cada elemnto finito
%Peso de la placa
j=j*9.81*10^(-6);s=0;
g=n1+n2; %numero total de elementos finitos
w=zeros(g+1);K44=zeros(g+1);
H=H1+H2; %Altura total
%Calculo de las areas de cada elemento finito y de la matriz de rigidez

for i=1:g
a(i)=(1-h(i)/(2*H)-s/H)*B*t;
s=s+h(i);
w(i,i)=1;w(i,i+1)=-1;w(i+1,i)=-1;w(i+1,i+1)=1;
K44=K44+a(i)*E/h(i)*w; %Matriz de rigidez
w=zeros(g+1);
end
a;
4
p=[];
K44 %Matriz de rigidez
%Cálculo de las fuerzas del vector carga global

for i=2:g
p(i)=-a(i-1)/2*h(i-1)*j-a(i)/2*h(i)*j+E*a(i-1)*c*dt-E*a(i)*c*dt;
if i==n1+1
p(i)=p(i)-pa;
end
end
p(g+1)=-a(g)/2*h(g)*j+E*a(g)*c*dt;p;
k44=K44(2:g+1,2:g+1);p=p(2:g+1);
%Calculo de los vectores desplazamiento

Q=inv(k44)*p';Q=[0;Q];r=K44(1,1:g+1)*Q;
%Calculo de la Reaccion en la base de la placa plana

R=r+a(1)/2*h(1)*j+E*a(1)*c*dt;
%Calculo de los esfuerzos globales en cada elemento finito

es=[];
for i=1:g
es(i,1)=E/h(i)*[-1 1]*Q(i:i+1,1)-E*c*dt;
end
N=n1+n2; %numero total de elemntos finitos
%Resultado de los calculo anteriores

disp('..............................');
disp('RESULTADOS');
disp(''),
disp('NUMERO DE TOTAL DE ELEMENTOS FINITOS UTILIZADOS');
disp(N);
disp('VECTOR DESPLAZAMIENTO');
disp(Q);
disp('REACCION EN EL APOYO(N)');
disp(R);
disp('VECTOR DE ESFUERZOS(MPa)');
disp(es);
5
Mostramos unas imágenes del script y del software Matlab con los resultados:
matriz rígida
6
7
CONCLUSIONES
 El método del cálculo por elementos finitos es una herramienta muy útil en el
cálculo de los esfuerzos y deformaciones que un cuerpo con geometría variable
pueda tener, ya que al ser un proceso corto pero de dificultad debido a la
cantidad de ecuaciones que se necesita resolver (éstas dependen del número de
elementos finitos con los que desee trabajar, por lo tanto si se desea obtener una
mejor precisión debemos tener un número alto de elementos finitos), la cantidad
de elementos finitos dificulta el cálculo pero mejora la precisión de los
resultados, por eso es necesario utilizar un software que agilice dicho cálculo,
como el Matlab.
 El informe se realizó para un cálculo de N números de elementos finitos, por lo
que nuestros resultados tienden a aproximarse a los reales, pero debido al uso
de bastantes elementos finitos fue necesario el uso de un software y con ayuda
de los métodos numéricos lograr resolver el problema el cual se puede hacer un
cálculo para innumerables elementos finitos los que requerirá un buen procesos
a fin de lograr soportar dichos cálculos.
 El valor de la resistencia que se obtiene en el apoyo es mayor al caso visto en el
anterior laboratorio, esto es debido al efecto de la temperatura sobre los
metales, ya que se trata de un T >0, el material se dilata en todas las direcciones,
al intentar desplazarse hacia la parte inferior el apoyo lo impide ocurriendo un
aumento de la resistencia debido al efecto de temperatura.
 Si el T fuese menor a 0, ocurriría un efecto contrario al explicado anteriormente
reduciendo así la resistencia en el apoyo. Por lo cual es importante estudiar el
efecto que tiene el cambio de temperatura en los elementos que deseamos
utilizar y así obtener unos mejores resultados a los esperados, ya que debido a
las temperaturas que están sometidos algunos elementos éstos fallan,
produciendo así pérdidas cuantiosas de dinero y tiempo.
 En el script del Matlab no se muestran algunos resultados como las matrices
para el cálculo de esfuerzos debido a que si colocamos un valor elevado de
elementos finitos el Matlab tendría demasiada información en la pantalla lo que
dificultaría observar los resultados. En el caso anterior se realizó para 10
elementos finitos lo cual se obtuvo una matriz de rigidez pequeña para poder
mostrar dicha matriz, en el caso fuese demasiado grande es necesario
esconderla.
8

Lab 2 Cálculo Por Elementos Finitos

Cargado por

Copyright:

Formatos disponibles

Lab 2 Cálculo Por Elementos Finitos

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Lab 2 Cálculo Por Elementos Finitos

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE INGENIERÍA

Facultad de Ingeniería Mecánica

UNIVERSIDAD NACIONAL DE INGENIERÍA

2DA PRÁCTICA CALIFICADA

ENUNCIADO DEL PROBLEMA

Calcular la distribución de esfuerzos a los que está sometido la placa triangular y

DIGITACIÓN DEL PROBLEMA, CÁLCULOS Y RESULTADOS

Tenemos el enunciado del problema y el script realizado en Matlab:

Realizando el Modelado del Cuerpo Real

Lo siguiente es el código del script de Matlab:

%Datos del problema

%Ingreso número de elementos finitos

%Calculo de las areas de cada elemento finito y de la matriz de rigidez

%Cálculo de las fuerzas del vector carga global

%Calculo de los vectores desplazamiento

%Calculo de la Reaccion en la base de la placa plana

%Calculo de los esfuerzos globales en cada elemento finito

%Resultado de los calculo anteriores

También podría gustarte