UNIVERSIDAD

NACIONAL DE
INGENIERIA

Clculo por Elementos

Finitos
|

FACULTAD DE INGENIERIA
MECNICA
CLCULO POR ELEMENTOS FINITOS
INFORME:

Laboratorio N 02 Segunda Prctica

Calificada

CDIGO:

SECCION:

DOCENTE:

MC 516

Ing. Ronald Cueva Pacheco

ALUMNO:
Cruz Saravia James J.

20112093E

UNI 2015 - I

NDICE
Pg.

ndice
. 1
Enunciado
del
Problema
.. 2
Solucin

... 3
Modelado
del
cuerpo
real...
3
Grados
de
libertad
nodales
.. 4
Vector
carga
..... 5
Matriz
de
Rigidez
. 7
Ecuaciones
de
rigidez
y
condiciones
contorno....8

de

Esfuerzos
..

. 9
Diagrama
de
flujo
.. 10
1

Uso
de
MATLAB
....11
Conclusiones
...13

SEGUNDA PRCTICA CALIFICADA

(TRACCION CON VARIACION DE TEMPERATURA)
ENUNCIADO DEL PROBLEMA
Dado la siguiente placa triangular, cuyo espesor es constante t=150mm se
encuentra afectado por la carga mostrada y a un aumento de temperatura
de 150 C. Calcular los esfuerzos en cada elemento finito y la reaccin en el
apoyo. Utilizar los elementos finitos necesarios.

Considerar:
Pa = 50KN
t (espesor) = 150 mm
E = 3.0x105 N/mm2
=11 106 C1
Y = 8.0gr-f/cm3 = 78,48x10-6 N/mm3

T =150 C

SOLUCIN
1. MODELADO DEL CUERPO REAL:
Segn el problema, tomaremos 4 elementos finitos para conseguir errores bajos.
Por ende, lo dividiremos segn las alturas

le

de 375 mm, 375mm, 375mm y

375mm respectivamente para cada elemento y as se facilitaran los clculos.

Luego de la divisin, hallaremos sus respectivas dimensiones de anchura b de

cada elemento finito, y por consiguiente tendremos el modelado del cuerpo como
se aprecia en la figura:

b 1=

1000+75 0
=875 mm
2

b 2=

750+500
=625 mm
2

b 3=

500+25 0
=375 mm
2

b 4=

25 0
=12 5 mm
2

El rea de cada elemento finito e lo calculamos geomtricamente de la forma:

NODOS

Elemento
e

Ae (

l e

GDL

mm

(mm)

Qi

Qj

(1)

(2)

Q1

Q2

375

131250

(2)

(3)

Q2

Q3

375

93750

(3)

(4)

Q3

Q4

375

56250

(4)

(5)

Q4

Q5

375

18750

A e =b e t
Con los resultados anteriores tenemos la siguiente tabla de conectividad:

2. GRADOS DE LIBERTAD NODALES (VECTOR

DESPLAZAMIENTO)
En el siguiente grafico se muestran los grados de libertad nodales globales:

Con ello se tiene el vector desplazamiento siguiente,

[] []

Q1
Q2
Qj= Q 3 ( mm ) Qj=
Q4
Q5

0
Q2
Q 3 ( mm )
Q4
Q5

Donde se tiene que por ser un cuerpo empotrado y fijo, Q1 ser igual a 0, luego de
ello los dems desplazamientos sern calculados.

3. VECTOR CARGA (Fi)

El vector carga se obtendr luego de hallar cada fuerza que acta en cada
elemento finito.

Analizando las fuerzas en cada elemento finitos por las relaciones de ecuaciones
que tenemos, se obtendr cada una de ellas como a continuacin:

N
2
mm
A nodo l e
e
1
f ( nodo )=
.Y + F +( E . A . . T )e .
donde
6 N
Y =78.48 10
2
1
3
mm
T =150 C

[ ]

E=3 10

f 1( 1) =

131250
.375 .78, 48. 106( 3.105 .131250 .11 .106 .150 ) + R 1=64966818. 6562+ R 1 ( N )
2

f 1( 2) =

131250
.375 .78, 48. 106 + ( 3. 105 .131250.11 . 106 .150 )=64970681.3438 ( N )
2

f 2( 2)=

93750
.375 .78 .48 .106 ( 3. 105 .93750 .11 . 106 .150 )=46404870.4688 ( N )
2

f 2( 3)=

93750
.375 .78,48 .106+ ( 3.10 5 .93750 .11 .106 .150 ) + Pa=46407629.5312+ Pa ( N )
2

f 3( 3) =

56250
.375 .78,48 .106 ( 3. 105 .56250 .11 . 106 .150 )=27842922.2812 ( N )
2

f 3( 4 ) =

56250
.375.78,48 . 106 + ( 3. 105 .56250 .11 . 106 .150 )=27844577.7188 ( N )
2

f 4( 4 ) =

18750
.375 .78,48 .106 ( 3. 105 .18750 .11 . 106 .150 )=9280974.0937 ( N )
2

f 4( 5 ) =

18750
.375.78,48 . 106 + ( 3. 105 .18750 .11 . 106 .150 )=9281525.9062 ( N )
2

Con ello tendremos la fuerza en todo el cuerpo.

F 1=f 1(1 ) =64966818.6562+ R 1 ( N )

F 2=f 1(2 )+ f 2(2 )=64970681.343846404870.4688=18565810.875 ( N )
F 3=f 2(3 )+ f 3( 3)=46457629.531227842922.2812=18614707.25 ( N )
F 4=f 3( 4) + f 4( 4)=27844577.71889280974.0937=18563603.6251 ( N )
F 5=f 4 (5)=9281525.9062 ( N )

Entonces el vector fuerza ser:

[ ][

F1
64966818.6562+ R 1
F2
18565810.875
Fi= F 3 =
(N)
18614707.25
F4
18563603.6251
F5
9281525.9062

4. MATRIZ DE RIGIDEZ
Para hallar la matriz de rigidez global usaremos de la siguiente ecuacin:

( A l E )

K ij =

1 1 0
1 1 0
0
0 0
0
0 0
0
0 0

0
0
0
0
0

] [ ] [ ] [

0
0
A E
0+
l
0
0

0
0

0
1

0 0
1
0 1 1 0
0 0 0 0
0 0 0 0

( A l E )

0 +
0
0

0 0
0 0

0 0
0 0
0 0 1 1
0 0 1 1
0 0 0 0

0
0
A E
0+
l
0
0

0
0

0
0
0

Reemplazando los valores ya anteriormente calculados que se muestran tambin

en la tabla de conectividad, nos resultar la matriz de rigidez:

131250 3. 105
K ij =
375

1 1 0
1 1 0
0
0 0
0
0 0
0
0 0

0
0
0
0
0

] [ ] [

0
0
93750 3. 105
0+
375
0
0

0
0

0
1

0 0
1
0 1 1 0
0 0 0 0
0 0 0 0

56250 3.10 5
0+
375
0
0

Resolviendo tenemos:

K ij =106 .

105 105
105 180
0
0
0

75
0
0

0 0 0
75 0 0
120 45 0
45 60 15
0
15 15

( mmN )

5. ECUACIN DE RIGIDEZ Y CONDICIN DE CONTORNO

La ecuacin de rigidez esta determinada por la siguiente ecuacin:

Fi=K ij Qj

Con los valores anteriormente obtenidos y usando la anterior ecuacin,

tenemos:

] [

64966818.6562+ R 1
18565810.875
Fi=
=
18614707.25
18563603.6251
9281525.9062

10 .

105 105
105 180
0
0
0

75
0
0

][ ]

0 0 0
0
75 0 0
Q2
120 45 0 Q3
Q4
45 60 15
Q5
0 15 15

Haciendo uso de propiedades, tomaremos una submatriz para hacer sencillo el

clculo de los desplazamientos.
1

0 0
0 0

) 00
3

0
0
0 0

] [

][ ]

180 75
0 0
18565810.875
Q2
18614707.25
12 0 45 0
Q3
6 75
=10 .

18563603.6251
Q4
0
4 5 6 0 15
9281525.9062
Q5
0
0
15 15

Resolviendo este sistema de ecuaciones, tendremos:

[ ][ ]

Q2
0.619291
Q3
1.238756
=
mm
Q4
1.857537
Q5
2.476305

Y para poder calcular la reaccin R1 usaremos la ecuacin:

[]

0
Q2
[ 64966818.6562+ R 1 ] =106 . [ 10 5 105 0 0 0 ] Q 3
Q4
Q5
Y as obtenemos que

R1 = -58736.3438 (N)

6. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la
siguiente ecuacin:

e=

[ ]

E
[ 1 1 ] Qi ( E . )e . T
le
Qi+1

( E . )e . T =495
Con los valores anteriores tenemos lo siguiente:

 1=

3 105 [
N
0
1 1 ]
495=0.4 328
2
375
0.619291
mm

3 105 [
N
0.619291
2=
1 1 ]
495=0.572
2
375
1.238756
mm
5

3=

3 10 [
N
1 1 ] 1.238756 495=0.0 248
2
375
1.857537
mm

4=

3 105 [
N
1.857537
1 1 ]
495=0.014 4
375
2.476305
mm2

Finalmente los resultados son los siguientes:

R 1=58736.3438 N
N
1=0.4328
mm 2
N
2=0.572
mm 2
N
3=0.0248
mm2
N
4 =0.0144
2
mm

7. DIAGRAMA DE FLUJO DEL PROCESO

INICIO
NUMERO DE ELEMENTOS FINITOS A UTILIZAR:
n = 4
MODELAMIENTO
CALCULO DE BASES:

bi

EFECTOS TERMICOS:

E A t
VECTOR
DESPLAZAMIENTO:

Qi

MATRIZ DE
RIGIDEZ:

Ki

VECTOR CARGA:

Fi

Fi =K i Qi

VALOR DE

Qi
REACCION

R1
ESFUERZOS

FIN

8. USO DEL MATLAB

Al hacer uso del programa MATLAB, una opcin es la siguiente para poder hallar
los resultados anteriores en n elementos finitos, para nuestro caso hallado
tomamos 4 elementos.
clear all
clc
%Ingreso de datos
disp('***Ingreso de Datos***');
n=input('Ingrese Nmero de Elementos Finitos=');
%Datos iniciales
L=1500; %en mm
w1=1000; %en mm

w2=0; %en mm
t=150; %en mm
P=50000; %P en Newton
E=3*(10^(5)); %E en Newton/mm^2
d=78.48*10^(-6);
alpha=11*10^(-6); %Valor de alpha en C^-1
dt=150; %Delta de temperatura en C
%calculo de las reas
l=L/n;
ll=(w1-w2)/n;
b=zeros(1,n);
for i=1:n
b(i)=(w1-i*ll);
end
a(1)=(w1+b(1))/2*t;
for j=2:n
a(j)=(b(j)+b(j-1))/2*t;
end
%Clculo de la Matriz de Rigidez y el Vector Fuerza
w=zeros(n+1);
K=zeros(n+1);
F=zeros(1,n+1);
FF=zeros(1,n);
ff=0;
for k=1:n
w(k,k)=1;w(k,k+1)=-1;w(k+1,k)=-1;w(k+1,k+1)=1;
K=K+a(k)*E/l*w;
F(k)=d*a(k)*l/2-E*a(k)*alpha*dt;
FF(k)=F(k)+2*E*a(k)*alpha*dt;
if(k>1 & k<n+1)
ff=FF(k-1);
F(k)=F(k)+ff;
end
w=zeros(n+1);
end
F(n+1)=d*a(n)*l/2+E*a(k)*alpha*dt;
h=10*25.4 % distancia de P en mm
disp('*****Clculos y Resultados*****');
disp('Resolviendo . . .');
%Aadiendo la carga P al Vector de Fuerza
for s=1:n
if(h==s*l)
F(s+1)=F(s+1)+P;
break;
end
end
%Grados de Libertad Nodales
Q=zeros(1,n+1);
%Eliminado el grado de libertad fijo

KK=K(2:n+1,[2:n+1]);
FF=F(1,[2:n+1]);
QQ=Q(1,[2:n+1]);
%Resolviendo las Ecuaciones
QQ=FF/KK;
Q=[Q(1),QQ];
%Clculo de los Esfuerzos
e=zeros(1,n);
for z=1:n
e(z)=E/l*[-1,1]*[Q(z);Q(z+1)]-E*alpha*dt;
end
%Clculando la Reaccin
R=K(1,:)*Q'-F(1);
%Mostrando los Resultados
disp('Matriz de Rigidez K');
disp(K);
disp('Vector de Fuerza F');
disp(F');
disp('Desplazamientos Q');
disp(Q');
disp('Esfuerzos');
disp(e');
disp('Reaccin en el grado de libertad 1');
disp(R);

9. CONCLUSIONES
El sentido asumido fue hacia abajo, con lo cual obtuvimos una
reaccin negativa lo que nos muestra que el sentido verdadero
de la reaccin es hacia arriba como que deba ser obtenido.
La reaccin que obtuvimos es parecida a la que se obtuvo
cuando no hubo una variacin de temperatura.
Con los resultados obtenidos podemos apreciar que las
deformaciones son pequeas (del orden de los micrmetros) y
nos resulta positivas como se deba de esperar ya que la carga
acta hacia abajo.

 Las deformaciones en este caso en el cual se varia la

temperatura nos resulta diferentes comparndolas con el caso
sin variacin de temperatura, lo cual resulta que la temperatura
afecta en gran parte a las deformaciones.
Los resultados de los esfuerzos nos resultaron positivos ya que la
carga tambin es positiva, lo que existe una traccin en dichos
puntos, estos a su vez resultaron en magnitud parecida al caso
sin variacin de temperatura.
La divisin del cuerpo en cuatro elementos nos resulta apreciable
y suficiente en la obtencin de errores comparando al usar el
MATLAB para varios elementos finitos.