UNIVERSIDAD NACIONAL DE INGENIERIA

FACULTAD DE INGENIERIA MECANICA

SEGUNDA PRACTICA (MC516-C)
LABORATORIO 2

Tema de laboratorio: TRACCIÓN CON DEFORMACIÓN TÉRMICA

Curso: Calculo por elementos Finitos

Profesor: Ing. Ronald Cueva Pacheco

Alumno: Obregón Aguilar Jahir 20172577I

Fecha de entrega: 2/ 05 / 2021

Índice

Problema...........................................................................................3

Solución.............................................................................................4

Diagrama de Flujo.............................................................................11

Uso de Matlab...................................................................................12

Conclusiones................................................................................... 13

2
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.

1000mm

750mm

1500mm
50000N

Datos:

E= 3.0x105 [N/mm2]

Y= 8.0gr-f/cm3= 78,45x10-6 [N/mm3]

3
SOLUCION:

1. MODELADO DEL CUERPO REAL

Se consideraran tres elementos finitos. Para facilitar los cálculos los elementos finitos tendrán
longitud de 100, 500 y 500mm. El espesor(t=120mm)

Cálculo de los espesores Calculo de las Áreas de canto

Para los primeros n/4 valores

1000+500
bi = = 750mm Ai =b i*t= 750*120=90 000 mm2
2

Para los segundos n/4 valores

500+250
b2 = = 375mm A2 =b 2*t= 375*120=45 000mm2
2

Para los últimos n/4 valores

250
b3 = = 125mm A3 =b 3*t= 125*120=15 000mm2
2

1000mm

1 A1
A2
2
A3
3
.
1500/n
.

.
1500mm .
. 1500/n
. 375mm
.
.
A(n-3)
n-2
A(n-2) 375mm
n-1
A(n-1)
n

4
Cuadro de conectividad:

NODOS GDL Le(mm) Ae (mm2)

e (1) (2) 1 2
1 2 1 2 1500/n 90 000
1

2 3 2 3 1500/n 90 000
2

3 4 3 4 1500/n 90 000
3

4 5 4 5 1500/n 45 000
4

…… …….. ….. ……. ……….. ……….

…

n-1 n-2 n-1 n-2 1500/n 15 000

n-1

n n-1 n n-1 1500/n 15 000

n

5
2. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)
Grados de libertad nodales globales

1000mm

1 A1 Q1
A2 Q2
2
A3 Q3
3
. 1500/n
.

.
1500mm .
50000
. 1500/n
. 375mm
.
.
A(n-3)
n-2
A(n-2) 375mm
n-1
A(n-1)
n

Qn

3. VECTOR CARGA

1000mm

1 A1 F1
A2 F2
2
A3 F3
3
. 1500/n
.

.
1500mm .
50000
. 1500/n
. 375mm
.
.
A(n-3)
n-2
A(n-2) 375mm
n-1
A(n-1)
n

Fn

6
4. DIAGRAMA DE FLUJO

INICIO

INGRESO DE DATOS
CONSTANTES: E,f,t
VECTORES:L.A.P

CALCULO DE VECTORES

AL1 γ EA 1 EA 1

[ ] [ ]
+ R1 − 0 0
2 L1 L1
AL2 γ AL 1 γ EA 1 EA 2
EA 1 EA 2
+ − + − 0
2 2 L1 L2 L1 L2
3
AL γ AL γ 2 EA 2 EA 3
EA 2 EA 3
+ +P A 0 − + −
2 2 L2 L 3
L2 L3
AL3 γ EA 3 EA 3
0 0 − 3 3
2 L L
F= ; K=

TRAFORMACION DE ECUACION MATRICIAL

AL1 γ EA 1

[ ][ ][]
−1 − 0 0
2 L1
AL2 γ AL1 γ 0
EA 2
+ 1
EA 1
−
EA 2
0
+ 2 R1
2 2 L L L2
AL3 γ AL2 γ EA 2 EA 3
EA 2 EA 3 Q2
+ +P A 0 − 2 3
+ −
2 2 L L L2 L3
3
EA 3 EA 3
Q3
AL γ 0 0 − 3
2 L L3 Q4
=

IMPRESIÓN DE RESULTADOS

R1 , Q 2 , Q 3 , Q 4 , E1 , E 2 , E 3

FIN

7
Programa en 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=300000;c=11*10^(-6);
dt=120; %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;
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=[];

8
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);

9
5. CONCLUSIONES

 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 .
 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.

10
11