Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
UNIVERSIDAD NACIONAL DE INGENIERA
Facultad de Ingeniera Mecnica
TRACCIN PURA
Curso Profesor Alumno :
: :
CALCULO POR ELEMENTOS FINITOS ING. CUEVA PACHECO RONALD
GARAY SALAZAR WILMER
20110079E
UNI, ABRIL DEL 2014
1. DESARROLLO NUMRICO DEL PROBLEMA
CLCULO POR ELEMENTOS FINITOS
Pgina 1
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
1.1. ENUNCIADO DEL PROBLEMA Dado la siguiente placa triangular, cuyo espesor es constante, con las caractersticas dadas, se pide calcular los esfuerzos en cada elemento finito y la reaccin en el apoyo. Utilizar tres elementos finitos.
1.2. MODELADO DEL CUERPO REAL Como piden tres elementos finitos, escogeremos tres elementos y para facilitar los clculos los dos primeros sern de longitud de 500mm y el tercero de 1000mm. Adems las reas de cada elemento finito se calculan de la siguiente manera:
El modelado del cuerpo quedara como se observa en la siguiente figura:
CLCULO POR ELEMENTOS FINITOS
Pgina 2
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
La tabla de conectividad seria el siguiente:
e 1 2 3
NODOS (1) 1 2 3
(2) 2 3 4
GDL 1 1 2 3
2 2 3 4
( ) 1000 500 500
( ) 45000 112500 157500
1.3.
GRADOS DE LIBERTAS NODALES (Vector Desplazamiento)
Pgina 3
CLCULO POR ELEMENTOS FINITOS
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
A travs del grafico se muestran los grados de libertad nodales globales:
1.4. VECTOR CARGA
CLCULO POR ELEMENTOS FINITOS
Pgina 4
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
Debido a que la densidad es:
Se hallara el peso con este valor asumiendo que este se distribuye de manera simtrica en cada nodo. Analizando las fuerzas en cada elemento finito:
( )
( ( ( ( )
) )
( )
CLCULO POR ELEMENTOS FINITOS
Pgina 5
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
( ( )
Ahora analizamos las fuerzas para todo el cuerpo:
Entonces el vector carga se expresara de la siguiente manera:
[ ]
1.5.MATRIZ DE RIGIDEZ A continuacin pasamos a calcular la matriz de Rigidez Global, que est determinada por la siguiente ecuacin:
) [
) [
) [
Reemplazando los valores calculados y utilizando la tabla de conectividad tenemos:
CLCULO POR ELEMENTOS FINITOS
Pgina 6
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
1.6. ECUACIN DE RIGIDEZ Y CONDICION DE CONTORNO La ecuacin de rigidez est determinada por la siguiente ecuacin:
Lo que con nuestros valores calculados tenemos:
] [
Resolviendo el sistema de ecuaciones tenemos:
CLCULO POR ELEMENTOS FINITOS
Pgina 7
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
1.7. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la siguiente ecuacin:
( ) [
Y obtenemos lo siguiente:
][
) [
][
) [
][
) [
][
CLCULO POR ELEMENTOS FINITOS
Pgina 8
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
2. DIAGRAMA DE FLUJO
INICIO
Ingreso de datos
CONSTANTES: E, P,H,B, t
Calculo de vectores
F=
K=
Para i desde 1 hasta 3
NO
SI
()
( ) ( )
() ) () ( ) ( ) () () ( )
Trasformacin de ecuacin matricial
CLCULO POR ELEMENTOS FINITOS
Pgina 9
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
Impresin de resultados
CONSTANTES: P; VECTORES: Q (4,1); e (3,1)
FIN
3. CODIFICACIN EN MATLAB
H=input('Ingrese la altura en mm = '); B=input('ingrese la base en mm = '); P=input('Ingrese la carga P en N = '); t=input('Ingrese el espesor en mm = '); den=input('Ingrese la densidad del material en gr-f/cm3 = '); E=input('Ingrese el modulo de elasticidad E(N/mm2) = '); h=[H/2 H/4 H/4]; den=den*9.81*10^(-6);s=0; w=zeros(4);K44=zeros(4);a=zeros(3,1); for i=1:3 a(i)=(s+h(i)/2)*B/H*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; w=zeros(4); end p=[]; p(1)=-a(1)/2*h(1)*den; p(2)=-a(1)/2*h(1)*den-a(2)/2*h(2)*den-P; p(3)=-a(2)/2*h(2)*den-a(3)/2*h(3)*den; k44=K44(1:3,1:3);Q=k44\p';Q=[Q;0];k=K44(4,1:4)*Q; R=k+a(3)/2*h(3)*den; esfuerzo=zeros(3,1); for i=1:3 esfuerzo(i,1)=E/h(i)*[-1 1]*Q(i:i+1,1); end clc;
CLCULO POR ELEMENTOS FINITOS
Pgina 10
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
%MOSTRANDO LOS RESULTADOS disp('..............................'); disp(' RESULTADOS DE CALCULO'); disp('============'); disp('EL VECTOR DESPLAZAMIENTO EN (mm)'); disp(Q); disp('LA REACCION EN EL APOYO(N)'); disp(R); disp('EL VECTOR DE ESFUERZOS(MPa)'); disp('esfuerzos en la placa'); disp(esfuerzo(1,1)'); disp(esfuerzo(2,1)'); disp(esfuerzo(3,1)');
VERIFICANDO LOS CLCULOS EN MATLAB
Ingrese Ingrese Ingrese Ingrese Ingrese Ingrese la la la el la el altura en mm = 2000 base en mm = 1200 carga P en N = 50000 espesor en mm = 150 densidad del material en gr-f/cm3 = 8 mdulo de elasticidad E(N/mm2) = 3*10^5
.............................. RESULTADOS DE CLCULO ======================== EL VECTOR DESPLAZAMIENTO EN (mm) -0.001602446984127 -0.001471646984127 -0.000645886243386 0 LA REACCIN EN EL APOYO(N) 6.412640000000001e+004 EL VECTOR DE ESFUERZOS (MPa) Esfuerzos en la placa 0.039240000000000 0.495456444444444 0.387531746031746
CLCULO POR ELEMENTOS FINITOS
Pgina 11
Universidad Nacional de Ingeniera Facultad de Ingeniera Mecnica
4. CONCLUSIONES
El MEF permite realizar un modelo matemtico de clculo del sistema real, ms fcil y simplificado sin embargo no es un mtodo exacto sino uno aproximado debido a las hiptesis bsicas del mtodo. Vemos que con el mtodo de elementos finitos podemos determinar distintos tipos de incgnitas interpretados como grados de libertad de un nodo son las variables que nos determinan el estado o posicin del nodo. Si deseamos una mejor aproximacin en el desarrollo de un problema es mejor dividir en una mayor cantidad de elementos, en este caso solo son tres por tanto el margen de error es considerable. Para este problema tambin se pudo elegir otra manera de fragmentar en tres elementos, es decir dos elementos sobre la mitad y otro elemento por debajo de la mitad, lo cual dara pequeos cambios en los clculos. El signo negativo de las deformaciones se debe a que el desplazamiento es hacia abajo, ya que los clculos se realizaron respetando el eje de las coordenadas. Vemos que las deformaciones de cada elemento es pequeo, lo cual es lgico de acuerdo a la experiencia y al curso de resistencia de materiales.
CLCULO POR ELEMENTOS FINITOS
Pgina 12
