UNIVERSIDAD NACIONAL DE
INGENIERA
FACULTAD DE INGENIERA
MECNICA
TRACCIN CON DEFORMACIN TRMICA
CURSO: CLCULO POR ELEMENTOS FINITOS
SECCIN: B
INTEGRANTE:
Tesen Atencia, Walter Junior
20090091E
FECHA: 25/09/2013
PROFESOR: Ing. Edwin Abreg
2013-2
MC 516-B
Pgina 1
NDICE
Enunciado del Problema....................................................................3
Preproceso.........................................................................................4
Proceso..............................................................4
Solucin.......................................................................................5
Ecuaciones de Rigidez y Contorno.....................................................6
Conclusiones................................................................................... 7
Bibliografa....................................................................................... 7
MC 516-B
Pgina 2
CLCULO POR ELEMENTOS FINITOS
Prob: Considere una barra con una seccin transversal que soporta una carga P, como se
muestra en la figura. La barra se fija en un extremo y lleva la carga P en el otro. Designamos el
ancho de la barra en la parte superior como w1, y la parte inferior como w2, el espesor como t
y la longitud por L. El modlo de elasticidad de la barra se denota por E. Considere que hay
un incremento de 110 C. Se pide determinar cunto se desviar la barra en varios puntos de
su longitud cuando se somete a la carga P.
Datos:
*w1=2 in=50.8 mm
*w2=1 in=25.4 mm
*t=0.125 in=3.175 mm
*L=10 in=254 mm
*E=10.4x106 lb/in2=7311.2 N/mm2
*P=1000 lb= 4445 N
*T=110 C
MC 516-B
Pgina 3
PREPROCESO
-
Discretizando el dominio de solucin en elementos finitos.
T=110 C
PROCESO
Con el fin de obtener los valores numricos de los desplazamientos nodales, se utilizan los
valores al comienzo mencionados, completndose la siguiente tabla:
MC 516-B
Pgina 4
La variacin del rea de la seccin transversal de la barra en la direccin y puede ser
expresado por:
25.4 50.8
) ) = (50.8 +
) (3.175)
254
() = 161.29 0.3175
() = ( + (
Con esto, podemos hallar la seccin transversal en cada nodo, que seran las siguientes:
A1=161.29 mm2
A2=141.129 mm2
A3=120.968 mm2
A4=100.806 mm2
A5=80.645 mm2
Luego, hallamos los k para cada elemento con la ecuacin:
= 17409.809 N/mm
= 15088.532 N/mm
= 12767.197 N/mm
= 10445.862 N/mm
+ )
2
Y la matriz elemental sera:
[K](1)=[K1 -K1;-K1 K1]
y as sucesivamente.
Juntando todas las matrices elementales se obtiene la matriz global de rigidez:
Analizando en cada elemento para hallar el vector de carga:
( )
( ) + = 28622.1 +
2
( )
=
( ) = 28622.105
2
( )
=
( ) = 24805.822
2
MC 516-B
Pgina 5
( )
( ) = 24805.822
2
( )
=
( ) = 20989.54
2
( )
=
( ) = 20989.54
2
( )
=
( ) = 17173.267
2
( )
=
( ) = 17173.267 + 4445
2
=
Ahora analizamos las fuerzas para todo el cuerpo:
= = 28622.1 + 1
= + = 3816.273
= + = 3816.292
= + = 2816.273
= = 21618.267
ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO
Se determina por:
Obteniendose los valores de los desplazamientos:
Q1 = 0 mm
Q2 = 0.19377 mm
Q3 = 0.39181 mm
Q4 = 0.59568 mm
Q5 = 0.80801 mm
Resolviendo el sistema de ecuaciones se halla el vector
hallamos de la siguiente submatriz.
MC 516-B
y para obtener el valor de R1 lo
Pgina 6
[28622.1 + 1] = [17409.809
17409.809 0
0 0]
[ ]
= 25248.60
Para hallar los esfuerzos resolvemos la siguiente ecuacin:
= ( ) [1 1] [
] ( )
Resolviendo para cada nodo se halla los siguientes valores:
= 34.32
= 40.90
= 50.46
= 29.50
MC 516-B
Pgina 7
CONCLUSIONES
Como se puede observar los cambios en la deformacin trmica de la barra es
considerable comparado con la deformacin que hay sin influencia de la temperatura.
MC 516-B
Pgina 8
BIBLIOGRAFA
-
INTRODUCCIN AL ESTUDIO DEL ELEMENTO FINITO EN INGENIERA, Tirupathi R.
Chandrupatla, Ashok D. Belegundu, Pearson Educacin, 1999
FINITE ELEMENT ANALYSIS THEORY AND APPLICATION WITH ANSYS, Saeed
Moaveni, 1999 by Prentice-Hall, Inc.
MC 516-B
Pgina 9
PROGRAMA EN MATLAB
clc
%constantes
format long;
E = 71722.73752;
f=7850*9.81*(10^-9);
t=3.175;
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) :');
dt=input('Ingrese la variacion de temperatura:');
%calculo de numero de nodos
nnodos=[1:n+m+1];
%calculo del vector L
L=[(127/n).*ones(1,n),(127/m).*ones(1,m)]
%calculo del vector H
HN1=[1000-250/n:-500/n:500+250/n];
HN2=[500-250/m:-500/m:0+250/m];
h=[HN1,HN2]
%calculo del vector P
P=zeros(1,n+m+1);
P(n+1)=4449.73752;
%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 fuerzapor 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
MC 516-B
Pgina 10
disp('La matriz de Esfuerzos en N/mm2
esf=(E.*diff(Q))./L'-E*x*dt
MC 516-B
es
: ' )
Pgina 11