UNIVERSIDAD NACIONAL DE INGENIERA
FACULTAD DE INGENIERA MECNICA
1ERA PRCTICA CALIFICADA
CURSO:
CALCULO POR ELEMENTOS FINITOS
TEMA:
TRACCIN SIMPLE
ALUMNO:
DAVIRAN MELGAREJO MICHAEL SMITH
SECCION:
MC 516 - D
PROFESOR:
CUEVA PACHECO RONALD
2015-I
PROBLEMA 01:
Para la placa triangular issceles que se muestra en la figura. Calcular la reaccin en la
base y los esfuerzos en cada elemento. Modele la placa con tres elementos finitos
DATOS:
Espesor uniforme= 150mm
E =3.0x105 N/mm2
= 8.0gr-f/cm3
PA= 30000N en el punto indicado.
SOLUCIN:
1.- MODELADO DEL CUERPO REAL
Se va a considerar tres elementos finitos. Para facilitar los clculos los dos primeros sern
de longitud de 250mm y el tercero de 500mm.
El ancho de cada elemento se calcula tomando el punto medio de cada elemento finito.
Hallando las bases medias por proporcionalidad:
800+
b3 =
800
2
b2=300 mm
b1=100 mm
=600 mm
Las reas se calculan de la siguiente relacin:
A i=bi x t
Luego se obtiene el Cuadro de conectividad:
e
1
2
3
NODOS
(1)
1
2
3
GDL
(2)
2
3
4
1
1
2
3
2
2
3
4
l e ( mm)
250
250
500
2.- GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)
A travs del grafico se muestran los grados de libertad nodales globales
Luego el vector de desplazamiento ser:
[]
Q1
Q
Q= 2 mm
Q3
0
A e (mm2 )
15000
45000
90000
Donde
Q4
= 0 pues la placa esta empotrada y los dems desplazamientos son
incgnitas que tendrn que ser calculadas.
3.- VECTOR CARGA:
Debido a que la densidad es: = 8.0gr-f/cm3= 8.0x10^-3gr-f/mm3
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:
F1 =
( Axl )1
=15000
2
F 2 1=
( Axl )1
=15000
2
F2 =
( Axl )2
=45000
2
F3 =
( Axl )2
FA =45000137755.1
2
F33 =
( Axl)3
=180000
2
F4 =
( Axl )3
+ R4 =180000+ R 4
2
Ahora analizamos las fuerzas para todo el cuerpo:
F1=F 11=15000 grf
F2 =F 21 + F22=60000 grf
F3 =F32 + F33 =182755.1 grf
3
F 4=F 4 + F 4 =180000+ R 4
Entonces el vector carga se expresara de la siguiente manera:
[ ][ ] [ ]
F1
15000
147.15
F
60000
588.6
F= 2 =
grf =
N
182755.1
1792.8275
F3
1765.8+ R 4
F 4 180000+ R 4
4.- MATRIZ DE RIGIDEZ
A continuacin pasamos a calcular la matriz de Rigidez Global, que est determinada por
la siguiente ecuacin:
AE
K i=
l
( )
] ( )[
1 1 0 0
AE
1 1 0 0
+
0
0 0 0
l
0
0 0 0
] ( )[
0 0
0 0
AE
0 1 1 0
+
0 1 1 0
l
0 0
0 0
0
0
0
0
0 0
0
0 0
0
0 1 1
0 1 1
Reemplazando los valores calculados y utilizando la tabla de conectividad tenemos:
1 1 0 0
15000 x 3 x 10 1 1 0 0 45000 x 3 x 105
K i=
+
250
0
0 0 0
250
0
0 0 0
5
0 0
0 0
5
0 1 1 0 + 90000 x 3 x 10
0 1 1 0
500
0 0
0 0
18000000 18000000
0
0
54000000
0
18000000 72000000
K i=
N / mm
0
54000000 162000000 108000000
0
0
108000000 108000000
5.- ECUACION DE RIGIDEZ Y CONDICION DE CONTORNO
La ecuacin de rigidez est determinada por la siguiente ecuacin:
Fi =K i x Q
Lo que con nuestros valores calculados tenemos:
][ ]
[ ][
Q1
147.15
18000000 18000000
0
0
Q
588.6
0
18000000 72000000 54000000
N=
x 2
1792.8275
0
54000000 162000000 108000000 Q 3
1765.8+ R 4
0
0
108000000 108000000
0
Para poder resolver este sistema de ecuaciones tomamos la siguiente submatriz:
][
][ ]
Q1
147.15
18000000 18000000
0
588.6 = 18000000 72000000 54000000 x Q 2
1792.8275
0
54000000 162000000
Q3
Resolviendo el sistema de ecuaciones tenemos:
[ ][ ]
Q 1 0.452
Q 2 = 0.370 m
Q 3 0.234
0
0
0
0
0 0
0
0 0
0
0 1 1
0 1 1
Para obtener la reaccin en el empotramiento tomamos la siguiente submatriz:
[]
Q1
Q
[1765.8+ R 4 ]=[ 0 0 108000000 108000000 ] Q2
3
0
Reemplazando los valores de Q obtenemos:
[ ]
0.452
0.370
[1765.8+ R 4 ]=[ 0 0 108000000 108000000 ] 0.234 x 103
0
R4 =41725.8 grf
6.-ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la siguiente ecuacin:
e=
( El ) [1
1]
[ ]
Qi
Qi+1
Y obtenemos lo siguiente:
3.0 x 10
1=
250
x 10
) [1 1 ][0.452
0.370 ]
=0.0984 MPa
3.0 x 105
250
x 10
) [1 1 ] [0.370
0.234 ]
=0.1632 MPa
3.0 x 105
500
x 10
) [1 1] [0.234
0 ]
=0.1404 MPa
2=
3=
7.- RESULTADOS
Finalmente, los resultados son mostrados en la siguiente tabla:
INCOGNITAS
VALORES
UNIDADES
R4
41725.8
grf
0.0984
MPa
0.1632
MPa
0.1404
MPa
8.- DIAGRAMA DE FLUJO
INICIO
INGRESO DE DATOS
CONSTANTES: E, f, t
VECTORES: L, A, P
CALCULO DE VECTORES
AL1
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
EA
L2
3
EA
EA2
L3
L2
3
EA
3
L
AL1
2
1
AL
AL
2
2
AL3
AL2
PA
2
2
AL3
EA1
L1
2
EA
EA1
L2
L1
EA 2
2
L
EA
L3
EA 3
L3
IMPRESIN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , E1 , E 2 , E 3
FIN
9.- FUNCIN EN MATLAB:
H=input('Ingrese la altura de la placa= ');
B=input('ingrese la base de la placa= ');
pa=input('Ingrese la carga PA= ');
t=input('Ingrese el espesor de la placa= ');
j=input('Ingrese la densidad del material= ');
E=input('Ingrese el modulo de elasticidad= ');
EA
2
L
EA3 EA 2
2
L3
L
EA3
3
L
TRAFORMACION DE ECUACION MATRICIAL
R1
EA 3
3
L
EA3
L3
Q
2
Q3
Q4
h=[3*H/10 3*H/10 4*H/10];
j=j*9.81*10^(-6);s=0;
w=zeros(4);K44=zeros(4);
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)=pa-a(1)/2*h(1)*j; p(2)=-a(1)/2*h(1)*j-a(2)/2*h(2)*j;
p(3)=-a(2)/2*h(2)*j-a(3)/2*h(3)*j+pb;
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)*j;
es=[];
for i=1:3
es(i,1)=E/h(i)*[-1 1]*Q(i:i+1,1);
end
clc;
disp('..............................'); %MOSTRANDO LOS RESULTADOS
disp(' RESULTADOS');
disp('============');
disp('EL VECTOR DESPLAZAMIENTO');
disp(Q);
disp('LA REACCION EN EL APOYO(N)');
disp(R);
disp('..............................');
disp('EL VECTOR DE ESFUERZOS(MPa)');
disp('
e1
e2
e3');
disp(es');
USO DEL PROGRAMA DE MATLAB
10.- CONCLUSIONES
Podemos apreciar que la placa es de un material rgido debido a que las
deformaciones son pequeas, esto se debe al mdulo de Young, ya que mientras
mayor es el valor mdulo de Young, ms rgido el material.
De los resultados se aprecia que todas la deformaciones son hacia abajo que es el
sentido opuesto al asumido, esto se debe a la carga de P A= 30000N que acta
sobre la placa.
Se aprecia que las deformaciones son debidas en mayor proporcin a la carga PA
que al peso de la placa, ya que PA es casi 6 veces el peso de la placa.
Se concluye que los esfuerzos son de compresin debido a que los valores
obtenidos en los esfuerzos son positivos.
