2da Practica Tracción Con Deformación Térmica

UNIVERSIDAD NACIONAL DE INGENIERA
FACULTAD DE INGENIERA MECNICA
Lab. Clculo por Elementos Finitos MC516-C

TRACCIN CON DEFORMACIN TRMICA
Alumno
Cdigo
CARDENAS LLICN, Jol Jess
Docente:
Ing. Ronald Cueva
Ciclo:
2015-II
20101029I
SEGUNDA PRCTICA
(TRACCIN CON DEFORMACIN TRMICA)
ENUNCIADO DEL PROBLEMA

De la siguiente placa triangular de espesor constante, t=150mm.
Calcular:
Los esfuerzos en cada elemento finito y la reaccin en el apoyo, con un
incremento de temperatura de 120C. Utilizar n elementos finitos.
Sabiendo que:
P =30000 N
T (espesor) = 150 mm
E = 3.0x105 N/mm2
= 8.0 gr-f/cm3 = 7,848x10-5 N/mm3
= 11x10-6 C-1
SOLUCIN:
1. MODELADO DEL CUERPO REAL
Consideramos seis elementos finitos de longitud de 125, 125, 125,
125, 250 y 250 mm desde la base hasta la punta.
El ancho de cada elemento lo calculamos tomando el punto medio
de cada elemento finito:
b1
b2
b3
b4
b5
b6
1200 1050
1125 mm
2
1050 900 975 mm
2
(900 750)
825 mm
2
750 600 675 mm
600 300
450 mm
2
(300 0)
150 mm
2
Luego:
Conectividad:
NODOS
e
le
Ae
(mm) (mm2)
GDL
(2)
(1)
Q1
Q2
125
Q2
Q3
125
Q3
Q4
125
Q4
Q5
125
5
6
5
6
6
7
Q5
Q6
Q6
Q7
250
250
Primer nodo
Segundo
nodo
16875
0
14625
0
12375
0
10125
0
67500
22500
2. GRADOS DE LIBERTAD NODALES.- (GDL)

(VECTOR DESPLAZAMIENTO)
En el siguiente grfico se muestran los vectores desplazamientos
nodales globales
El vector de desplazamiento ser:

0
Q 2
Q3
mm
Q 4
Q5
Q6
Q7
Donde Q1= 0 debido a que la placa esta empotrada y los dems

desplazamientos son incgnitas donde procederemos a calcularlos.
3. VECTOR CARGA
Analizando las fuerzas en cada elemento finito:
F11=
( Al )1
( EA T )1 + R1=66825827,72+ R1
2
F2 =
( Al )1
+ ( EA T )1=66824172,28
2
F22 =
( Al )2
( EA T )2=57915717,36
2
F3 =
( Al )2
+ ( EAT )2=57914282,64
2
F33 =
( Al )3
( EA T )3 =49005606,99
2
F 4=
( Al )3
+ ( EA T )3 =49004393,01
2
F 44=
( Al )4
( EA T )4=40095496,63
2
F5 =
( Al )4
+ ( EAT )4 P=40064503,37
2
F55 =
( Al )5
( EA T )5 =26730662,18
2
F6 =
F66 =
( Al )5
+ ( EA T )5=26729337,83
2
( Al )6
( EA T )6=8910220,725
2
6
F7 =
( Al )6
+ ( EA T )6 =8909779,275
2
Ahora analizamos las fuerzas para todo el cuerpo:

1
F1=F 1=66825827,7+ R1
1
F2 =F 2 + F 2=8908454,93
F3 =F3 + F 3=8908675,65
F 4=F 34 + F 44=8908896,38
F5 =F45 + F 55=13333841,2
F6 =F56 + F 66=17819117,1
F7 =F67 =8909779,28
Entonces, el vector carga se expresara de la siguiente manera
[ ][ ]
F 1 66825827,7+ R 1
F2
8908454,93
F3
8908675,65
F= F 4 =
8908896,38
13333841,2
F5
17819117,1
F6
8909779,28
F7
1. MATRIZ DE RIGIDEZ
A continuacin pasamos a calcular la matriz de Rigidez Global, que
est determinada por la siguiente ecuacin:
K IJ =
( AEL )
( AEL )
0
0
0
0
0
0
0
][ ]
][ ]
1 1 0 0 0 0 0
1 1 0 0 0 0 0
0
0 0 0 0 0 0
AE
0
0 0 0 0 0 0+
L
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
0 0
0 0 0 0
0 1 1 0 0 0
AE
0 1 1 0 0 0 +
L
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
0 0 0 0 0
0 1 1 0 0 0 0
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
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
0 1 1 0 0
0 1 1 0 0
0 0 0
0 0
0 0 0
0 0
( AEL )
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
0 0
0 0
0 0
AE
0 0
0 0+ L
0 1 1 0
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
0
0
0
0
0
0
0
0 0
0 0
0
0 0
0
0 0
0
0 0
0
0 1 1
0 1 1
Reemplazando para los valores calculados y utilizando la tabla de

conectividad obtenemos:
K IJ =
1687503105
125
1237503105
125
67500310 5
250
[
[
] [ ]
] [ ]
] [ ]
1 1 0 0 0 0 0
1 1 0 0 0 0 0
0
0 0 0 0 0 0
146250310 5
0
0 0 0 0 0 0+
125
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
0
0 0
0 0 0 0
0 0
0 0 0 0
0 1 1 0 0 0
1012503105
0 1 1 0 0 0 +
125
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
0
0
0
0
0
0
0
0
0
0
0
0
0 0
0 0
0 0
0 0
0 0
0 0
225003105
+
0 0
0 0
250
0 1 1 0
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
0
0
0
0
0
0
0 0
0 0 0 0 0
0 1 1 0 0 0 0
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
0 0 0 0 0
0 0
0 0 0
0 0
0 0 0
0 0
0 0 0
0 1 1 0 0
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
0 0
0 0
0
0 0
0
0 0
0
0 0
0
0 1 1
0 1 1
Finalmente:
1350 1350
0
0
0
0 0
1350 2520 1170
0
0
0
0
0
1170 2160 990
0
0
0
5
K IJ =310
0
0
990 1800 810
0
0
810 1080 270
0
0
0
0
0
270
360
90
0
0
0
0
0
90
90
0
0
0
2. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO

La ecuacin de rigidez est determinada por la siguiente ecuacin:
Fi K i Q
Lo que con nuestros valores calculados tenemos:
[ ][
1350 1350
0
0
0
0 0
66825827,7+R 1
1350
2520
1170
0
0
0
0
8908454,93
0
1170
2160
990
0
0
0
8908675,65
5
=310
0
0
990 1800 810
0
0
8908896,38
810 1080 270
0
13333841,2
0
0
0
17819117,1
0
270 360 90
0
0
0
8909779,28
0
0
90
90
0
0
0
10
Para obtener los desplazamientos tomamos la siguiente submatriz:
[ ]
8908454.93
8908675,65
8908896.38
13333841,2
17819117,1
8909779,28
2520
1170
1170
2160
3 x 10 5 x
0
0
990
0
Q2
990
0
0
0 Q3
1800 810
0
0 Q4

810
1080 270
0 Q5
Q6
0
270 360
90

0
0 90 90
Q7
0
Resolviendo este sistema de ecuaciones obtenemos:

Q2 16491.053 x10^ (-5) mm
Q3 32981.1697 x10^ (-5) mm
Q4 49469.9353x10^ (-5)mm
Q5 65956.6586x10^ (-5) mm
Q6 98955.2961 x10^ (-5) mm
Q7 131954.4786 x10^ (-5) mm
Y para obtener la reaccin en el empotramiento tmanos la siguiente

submatriz:
11
0
Q
2
66825827,7 R1 3x10 5 x1350
Q3
1350
0 0 Q 4
Q5
Q6
Q
7
Resolviendo obtenemos:
R1 37063.05 N
3. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos
la siguiente ecuacin:
E

l
e
Qi
e
( E ) T
Q
1
i
1 1
Donde:
( E ) T
e
= (3*105*11*10-6)*150=396
N
mm2
Donde obtenemos lo siguiente:
3 x 10 5
125
1 1
1
396
0.1649105
12
1 0.2148
N
mm 2
3 x 10 5
2
125
3 x 10 5
125
3 x 10 5
125
3 x 10 5
250
3 x 10 5
6
250
0.1649105
396
0.3298116
2 - 0.23736
N
mm 2
1 1
0.3298116
396
0.4946993
3 -0,26952
N
mm 2
1 1
1 1
0.4946993
N
396 4 0.31872
mm 2
0.6595665
1 1
0.6595665
N
396 5 -0.01632
mm 2
0.9895529
1 1
3
0.9895529
N
396 6 -0.00984
mm 2
1.3195447
4. RESULTADOS
Finalmente, los resultados son mostrados en la siguiente tabla:
R1 37063.5 N
1 0.2148
N
mm 2
2 - 0.23736
N
mm 2
3 - 0.26952
N
mm 2
4 0.31872
N
mm 2
5 - 0.01632
N
mm 2
13
N
mm 2
6 - 0.00984
5. DIAGRAMA DE FLUJO
INICIO
INGRESO DE DATOS
CONSTANTES : E, f, t, T
VECTORES: L, A, P
CALCULO DE VECTORES
AL
R1
2
1
AL AL
2
2
AL3 AL2
PA
2
2
AL3
F=
EA1
L1
EA1
1
L
K=
EA1
L1
2
EA
EA1
L2
L1
2
EA
2
L
AL
2
1
AL
AL
2
2
AL3
AL2
PA
2
2
AL3
0
=
EA1
L1
2
EA
EA1
L2
L1
2
EA
2
L
EA
L2
3
EA
EA 2
L3
L2
3
EA
3
L
14
0
EA 2
2
L
EA3 EA 2
2
L3
L
EA3
3
L
EA3
3
L
EA3
L3
R1
Q
2

Q3
EA3
L3
3
EA
L3
TRAFORMACION DE ECUACION MATRICIAL

1
Q
4
Q5
Q
6
IMPRESIN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , Q5 , Q6 , Esf 1 , Esf 2 , Esf 3 , Esf 4 , Esf 5
FIN
6. DIGITACIN EN MATLAB
clc
%constantes
format long;
E = 300000;
f=8*9.81*(10^-6);
t=150;
dt=120;
x=11*(10^-6);
%condicion de contorno
contorno=[1];
n=input ('Ingrese el N de elementos finitos de la primera parte (n) :');
m=input ('Ingrese el N de elementos finitos de la segunda parte (m) :');
%calculo de numero de nodos
nnodos=[1:n+m+1];
%calculo del vector L
L=[(500/n).*ones(1,n),(500/m).*ones(1,m)]
%calculo del vector H
HN1=[1200-300/n:-600/n:600+300/n];
HN2=[600-300/m:-600/m:0+300/m];
h=[HN1,HN2]
%calculo del vector P
P=zeros(1,n+m+1);
P(n+1)=-30000;
%calculo del vector h
A=t.*h;
%Calculo del vector fuerza masica
pesoe=-f*0.5.*A.*L;
pesoglobal=[pesoe,0]+[0,pesoe];
%calculo del vector fuerza por deformacion termica
15
dtele=(x*dt*E).*A;
dtglobal=-[dtele,0]+[0,dtele];
%calculo de vector fuerza total
F=pesoglobal+P+dtglobal;
%Calculo de la matriz de rigidez
kele=E.*A./L;
disp('La matriz de Rigidez es : ')
kglobal=diag([kele,0]+[0,kele])-diag(kele,1)-diag(kele,-1)
%resolucion de la ecuacion matricial [F]=[K].[Q]
%introduccion de las condiciones de contorno
ncon=setdiff(nnodos,contorno);
KM=kglobal(ncon,ncon);
FM=F(ncon);
QM=inv(KM)*(FM)';
%calculo de los dezplazamientos
Q=zeros(size(nnodos,2),1);
Q(ncon)=QM;
R=kglobal*Q-F'; %Matriz de Reacciones
disp('La matriz de desplazamiento en mm es :' );
disp(Q)
disp('La fuerza (N) punto de apoyo es : ')
disp(R(1))
%calculo de esfuerzos
disp('La matriz de Esfuerzos en N/mm2 es : ' )
esf=(E.*diff(Q))./L'-E*x*dt
Resultados obtenidos en Matlab:

Ingrese el N de elementos finitos de la primera parte (n) :4
Ingrese el N de elementos finitos de la segunda parte (m) :2
L=
125 125 125 125 250 250
h=
Columns 1 through 5
16
1125
975
825
675
450
Column 6
150
La matriz de Rigidez es :
kglobal =
Columns 1 through 5
405000000 -405000000
0
0
0
-405000000 756000000 -351000000
0
0
0 -351000000 648000000 -297000000
0
0
0 -297000000 540000000 -243000000
0
0
0 -243000000 324000000
0
0
0
0 -81000000
0
0
0
0
0
Columns 6 through 7
0
0
0
0
-81000000
108000000
-27000000
0
0
0
0
0
-27000000
27000000
La matriz de desplazamiento en mm es :
0
0.164910529675926
0.329811696609687
0.494699352985949
0.659566585779159
0.989552960779159
1.319544785779159
La fuerza (N) punto de apoyo es :
3.706319999998808e+004
La matriz de Esfuerzos en N/mm2 es :
esf =
17
-0.214728777777680
-0.237199358974294
-0.269624696969800
-0.318641296295937
-0.016349999999932
-0.009810000000016
18
7. CONCLUSIONES
A medida que se toma mayor nmero de elementos finitos mejor

ser el resultado encontrado el error va a ser mnimo conforme
tomamos ms elementos finitos.
El aumento de temperatura hace cambiar notablemente el valor

deformacin por ende de los esfuerzos en 1000 veces ms, por
tanto es muy necesario hacer un estudio con anlisis con
elementos finitos a una pieza mecnica si este va a trabajar a
constantes cambios de temperatura.
Al comparar la reaccin hallada analticamente con que la que el

Matlab se encuentra un error, aunque es pequeo, podra influir en
luego en errores futuros, este error es muy probablemente
producido por error en los clculos y el uso de aproximaciones
decimales.
Cuando solo haba traccin en el cuerpo, se poda notar fcilmente

que los desplazamientos eran negativos, ya que la traccin estaba
en contra del sistema de referencia, pero ahora al haber un
incremento de temperatura, los desplazamientos son positivos,
esto
significa
que
el
efecto
de
temperatura
gana
en
desplazamiento al efecto de traccin en el cuerpo y en cada

elemento.
Aun as, los esfuerzos son negativos, debido a que de todos modos
el efecto de traccin produce un esfuerzo en contra del eje de
referencia.
Se nota que con o sin variacin de temperatura obtuvimos un

valor de reaccin igual a 37063,2 N.
19
20
8. BIBLIOGRAFA
CHANDRUPATLA, T. Introduccin al Estudio de los Elementos
Finitos en Ingeniera, Prentice Hall, 1999
ZIENKIEWCTZ, O. The Finite Element Method, New Cord, Mec
Graw Hill, 1977.
ZIENKIEWCTZ, O. and MORGAN K. Finite Elements and
Approximation, New Cork, Wiley, 1982.
LIVESLEY, R. Finite Element: An Introduction for Engineers,
Cambridge, Great Britain, Cambridge University Press, 1983.
21

2da Practica Tracción Con Deformación Térmica

Cargado por

Copyright:

Formatos disponibles

2da Practica Tracción Con Deformación Térmica

Cargado por

Información del documento

Descripción 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

2da Practica Tracción Con Deformación Térmica

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE INGENIERA

FACULTAD DE INGENIERA MECNICA

Lab. Clculo por Elementos Finitos MC516-C

CARDENAS LLICN, Jol Jess

Ing. Ronald Cueva

ENUNCIADO DEL PROBLEMA

2. GRADOS DE LIBERTAD NODALES.- (GDL)

El vector de desplazamiento ser:

Donde Q1= 0 debido a que la placa esta empotrada y los dems

Analizando las fuerzas en cada elemento finito:

Ahora analizamos las fuerzas para todo el cuerpo:

Reemplazando para los valores calculados y utilizando la tabla de

2. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO

Lo que con nuestros valores calculados tenemos:

Para obtener los desplazamientos tomamos la siguiente submatriz:

Resolviendo este sistema de ecuaciones obtenemos:

Y para obtener la reaccin en el empotramiento tmanos la siguiente

66825827,7 R1 3x10 5 x1350

Donde obtenemos lo siguiente:

TRAFORMACION DE ECUACION MATRICIAL

R1 , Q2 , Q3 , Q4 , Q5 , Q6 , Esf 1 , Esf 2 , Esf 3 , Esf 4 , Esf 5

Resultados obtenidos en Matlab:

A medida que se toma mayor nmero de elementos finitos mejor

El aumento de temperatura hace cambiar notablemente el valor

Al comparar la reaccin hallada analticamente con que la que el

Cuando solo haba traccin en el cuerpo, se poda notar fcilmente

desplazamiento al efecto de traccin en el cuerpo y en cada

Se nota que con o sin variacin de temperatura obtuvimos un

También podría gustarte