Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
“TRACCIÓN CON DEFORMACIÓN TÉRMICA”
CURSO: CÁLCULO DE ELEMENTOS FINITOS
AUTOR:
RUIZ QUISPE FRANKIE 20122202A
SECCIÓN: C
SETIEMBRE - 2014
1
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
OBJETIVOS
Analizar la estructura de una barra de sección variable.
Utilizar métodos matemáticos (modelamiento) para aproximar la forma
de la sección a una más sencilla de estudiar.
Calcular los esfuerzos en determinadas secciones de la barra de sección
variable producidos por una fuerza.
Someter a esa tracción una alteración debida a un cambio de
temperatura.
La reacción en los apoyos.
2
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
SEGUNDA PRÁCTICA
ENUNCIADO DEL PROBLEMA
De la siguiente placa triangular de espesor constante, t=150mm. Calcule:
Los esfuerzos en cada elemento finito y la reacción en el apoyo, con un incremento de
temperatura de 80°C. Utilizar n elementos finitos.
Sabiendo que:
P =30000 N
T (espesor) = 150 mm
E = 3x105 N/mm2
γ = 8.0 gr-f/cm3 = 7,848x10-5 N/mm3
α = 11*10-6 °C-1
3
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
SOLUCIÓN:
1. MODELADO DEL CUERPO REAL
Consideramos seis elementos finitos de longitud de 250, 250, 250, 250, 500 y
500 mm desde la base hasta la punta.
El ancho de cada elemento lo calculamos tomando el punto medio de cada
elemento finito:
b1
1200 1050 1125 mm
2
b2
1050 900 975 mm
2
(900 750)
b3 825 mm
2
b4
750 600 675 mm
2
b5
600 300 450 mm
2
(300 0)
b6 150 mm
2
4
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
Luego:
Conectividad:
NODOS GDL
le Ae
| (1) (2) (mm) (mm2)
1 2
Primer nodo Segundo nodo
1 1 2 Q1 Q2 250 168750
2 2 3 Q2 Q3 250 146250
3 3 4 Q3 Q4 250 123750
4 4 5 Q4 Q5 250 101250
5 5 6 Q5 Q6 500 67500
6 6 7 Q6 Q7 500 22500
2. GRADOS DE LIBERTAD NODALES.- (GDL)
(VECTOR DESPLAZAMIENTO)
En el siguiente gráfico se muestran los vectores desplazamientos nodales
globales:
5
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
El vector de desplazamiento será:
0
Q 2
Q3
Q Q 4 mm
Q5
Q 6
Q 7
Donde Q1=0 debido a que la placa esta empotrada y los demás
desplazamientos son incógnitas donde procederemos a calcularlos.
3. VECTOR CARGA
6
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
Analizando las fuerzas en cada elemento finito:
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
7
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
Analizando las fuerzas para todo el cuerpo:
Entonces, el vector carga se expresaría de la siguiente manera
[ ] [ ]
4. MATRIZ DE RIGIDEZ
A continuación pasamos a calcular la matriz de Rigidez Global, que está
determinada por la siguiente ecuación:
( ) ( )
[ ] [ ]
8
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
( ) ( )
[ ] [ ]
( ) ( )
[ ] [ ]
Reemplazando para los valores calculados y utilizando la tabla de conectividad
obtenemos:
( )
[ ]
( )
[ ]
( )
[ ]
( )
[ ]
9
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
( )
[ ]
( )
[ ]
Finalmente:
[ ]
5. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO
La ecuación de rigidez está determinada por la siguiente ecuación:
Fi K i Q
Con nuestros valores calculados tenemos:
[ ] [ ] [ ]
10
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
Resolviendo este sistema de ecuaciones obtenemos:
Q2 0.22022 mm
Q3 0.440455mm
Q4 0.660716mm
Q5 0.881008 mm
Q6 1.32106 mm
Q7 1.7611 mm
R1 46112.9 N
6. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la siguiente
ecuación:
Q
e
E
e 1 1 i
( E ) e T
l Qi 1
Donde:
N
( E )e T = (3*105*11*10-6)*80=264
mm2
Donde obtenemos lo siguiente:
3 x105 0
1 1
N
1 264 1 0.264
250 1 0.22022 mm2
3 x105 0.22022
1 1
N
2 264 2 0.282
250 2 0.440455 mm2
3 x105 0.440455
1 1
N
3 264 3 0.3132
250 3 0.660716 mm2
11
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
3 x105 0.660716
1 1
N
4 264 4 0.3504
250 3 0.881008 mm2
3 x10 5 0.881008
1 1
N
5 264 5 0.0312
500 3 1.32106 mm2
3 x10 5 1.32106
1 1
N
6 264 6 0.024
500 3 1.7611 mm2
7. RESULTADOS
R1 46112.9 N
N
1 0.264
mm2
N
2 0.282
mm2
N
3 0.3132
mm2
N
4 0.3504
mm2
N
5 0.0312
mm2
N
6 0.024
mm2
12
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
8. DIAGRAMA DE FLUJO
INICIO
INGRESO DE DATOS
CONSTANTES : E, f, t
VECTORES: L, A, P
CALCULO DE VECTORES
AL1 EA1 EA1
R1 0 0
1
L1
2 L 1 2
EA1 EA 2
AL AL EA EA
0
2 1
F= 2 2 ; K= L1 L2 L1 L2
AL3 AL2 EA 2
EA 3
EA 2 EA 3
PA 0 2 3
2 2 L L3 L2 L
AL3
0 EA 3
EA 3
0 3
2 L L3
13
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
TRAFORMACION DE ECUACION MATRICIAL
AL1 EA1
1
L1
0
0
R1
2
AL2 AL1 EA 2
EA1 EA 2 Q
0 2 0 2
2 2 = L2 L1 L
AL3 AL
2 EA 2
EA 3 EA 2 EA 3
Q3
PA 0 2 2 3
2 2 L L3 L L Q 4
AL3 0 EA 3 EA 3
0 3
2 L L3
CALCULO DE ESFUERZOS
1 1
E
e
Qi
e
l Qi 1
IMPRESIÓN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , Q5 , Q6 , Esf1 , Esf 2 , Esf 3 , Esf 4 , Esf 5
FIN
14
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
9. DIGITACIÓN EN MATLAB
clc
%constantes
format long;
E = 300000;
f=8*9.81*(10^-6);
t=150;
dt=80;
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=[(1000/n).*ones(1,n),(1000/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
15
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
pesoe=f*0.5.*A.*L;
pesoglobal=[pesoe,0]+[0,pesoe];
%calculo del vector fuerza por deformacion termica
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
16
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
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=
250 250 250 250 500 500
h=
Columns 1 through 5
1125 975 825 675 450
Column 6
150
La matriz de Rigidez es :
kglobal =
Columns 1 through 5
202500000 -202500000 0 0 0
-202500000 378000000 -175500000 0 0
0 -175500000 324000000 -148500000 0
0 0 -148500000 270000000 -121500000
0 0 0 -121500000 162000000
0 0 0 0 -40500000
0 0 0 0 0
Columns 6 through 7
0 0
0 0
0 0
0 0
-40500000 0
54000000 -13500000
-13500000 13500000
17
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
La matriz de desplazamiento en mm es :
0.220209733148148
0.440434125242165
0.660681479535095
0.880965634782008
1.321020134782008
1.761052834782009
La fuerza (N) punto de apoyo es :
-4.412640000000596e+004
La matriz de Esfuerzos en N/mm2 es :
esf =
0.251679777777838
0.269270512820640
0.296825151515066
0.340986296296421
0.032700000000034
0.019620000000089
18
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
10. CONCLUSIONES
Mayor cantidad de divisiones permitirá una aproximación más precisa del
análisis que se requiera dar a un elemento con geometría complicada.
Los cálculos a mano y los hechos en Matlab difieren en aproximadamente un
4.347%, un error aceptable (cálculos de piezas pequeñas) pues el estándar es
un error de 5%.
El aumento de temperatura hace cambiar notablemente el valor deformación
por ende de los esfuerzos en 1000 veces más, por tanto es muy necesario
hacer un estudio con análisis con elementos finitos a una pieza mecánica si
este va a trabajar a constantes cambios de temperatura.
Cuando solo había tracción en el cuerpo, se podía notar fácilmente que los
desplazamientos eran negativos, ya que la tracción 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 tracción en el cuerpo y en cada elemento.
Los esfuerzos son positivos, debido a que de todos modos el efecto de tracción
produce un esfuerzo a favor del eje de referencia.
Se nota que con o sin variación de temperatura obtuvimos un valor de reacción
igual a 441264.4 N.
La no inclusión de la carga PA=30000 N, provocaría que los cálculos se basen
en la variación de temperatura, produciendo esfuerzos en sentido contrario al
eje de referencia.
19
Segunda Práctica: TRACCIÓN CON DEFORMACIÓN TÉRMICA UNI-FIM
11. BIBLIOGRAFÍA
CHANDRUPATLA, T. “Introducción al Estudio de los Elementos Finitos en
Ingeniería”, 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.
20
