Índice

Enunciado del Problema....................................................................3

Solución.............................................................................................4

Grados de Libertad Nodales..............................................................5

Vector Carga......................................................................................6

Matriz de Rigidez................................................................................8

Ecuación de Rigidez y Condición de Contorno..................................9

Esfuerzos y Resultados....................................................................10

Diagrama de Flujo.............................................................................11

Uso de Matlab...................................................................................12

Conclusiones................................................................................... 14

1
 PRIMERA PRÁCTICA CALIFICADA
(TRACCION SIMPLE)

ENUNCIADO DEL PROBLEMA

Dado la siguiente placa triangular, cuyo espesor es constante,
t=150mm, calcular los esfuerzos en cada elemento finito y la
reacción en el apoyo. Utilizar tres elementos finitos.

2
 Considerar:
PA = 10KN
t (espesor) = 150 mm
E = 3.0x105 N/mm2
Y = 8.0gr-f/cm3 = 78,45x10-6 N/mm3

SOLUCION:
1. MODELADO DEL CUERPO REAL
Se consideraran tres elementos finitos. Para facilitar los cálculos los
elementos finitos tendrán longitud de 500, 250 y 250mm.

Y los espesores lo calculamos tomando el punto medio de cada

elemento finito:

b1 
1200  600 900 mm
2
b2 
 600  300 450 mm
2
300
b1   150 mm
2

Entonces, el modelado del cuerpo sería el siguiente:

3
Y las áreas se calculan de la siguiente relación:
A1  b1 x t

Cuadro de conectividad:

NODOS GDL le Ae
e (1) (2) 1 2 (mm) (mm2)
1 1 2 1 2 500 135000
2 2 3 2 3 250 67500
3 3 4 3 4 250 22500

2. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)

4
A través del grafico se muestran los grados de libertad nodales
globales:

Luego el vector de desplazamiento será:

0 
Q2
Q   mm
Q3

Q4

5
Donde Q1= 0 pues la placa esta empotrada y los demás
desplazamientos son incógnitas que tendrán que ser calculadas.

3. VECTOR CARGA

Analizando las fuerzas en cada elemento finito:

6
 y  Axl  1
F11  R1  2647.68  R1 N
2
y  Axl  1
F21   PA  12647.68 N
2
y  Axl  2
F22   661.92 N
2
y  Axl  2
F32   661.92 N
2
y  Axl  3
F33   220.64 N
2

y  Axl  3
F43   220.64 N
2

Ahora analizamos las fuerzas para todo el cuerpo:

F1  F11  2647.68  R1 N
F2  F21  F22  13309.6 N
F3  F32  F33  882.56 N
F4  F43  220.64 N

Entonces, el vector carga se expresaría de la siguiente manera

7
 F1 2647.68R1
F2  13 09.6 
F1    N
F3 8 2.56 
  
F4 2 0.64 

4. MATRIZ DE RIGIDEZ

A continuación pasamos a calcular la matriz de Rigidez Global, que

esta determinada por la siguiente ecuación:

1 1 0 0  0 0 0 0   0 0 0 0 
 1 1 0 0     
 AE     AE  0 1 1 0   AE   0 0 0 0 
Ki       
 l 1  0 0 0 0   l  2  0  1 1 0   l  3  0 0 1  1 
   0 0 0 0   0 0 1 1 
0 0 0 0    

8
Reemplazando para los valores calculados y utilizando la tabla de
conectividad obtenemos:

1 1 0 0  0 0 0 0 
5  
135000x3x10  1 1 0 0   67500x3x105  0 1 1 0 
Ki       
 500 1  0 0 0 0   250 2  0 1 1 0 
 
0 0 0 0 0 0 0 0 

0 0 0 0
 
 22500x3x105   0 0 0 0 
  
 250 3  0 0 1 1 
 
 0 0 1 1 
Finalmente:

9
  810  810 0 0 
 810 1620  810 0 
5   N
Ki 10 x
 0  810 1080  270 mm
 
 0 0  270 270 

5. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO

La ecuación de rigidez esta determinada por la siguiente ecuación:

Fi  K i  Q

Lo que con nuestros valores calculados tenemos:

 810 810 0 0 0 

2647.68  R1 
 13309.6   810 1620  810 0 Q2
  
   10 5 x   
  0  810 1080  270 Q3 
 882.56 

 220.64 
  
 0 0  270 270 Q4
10
Para obtener los desplazamientos tomamos la siguiente submatriz:

 1620 810 0 Q2

 882.56 5 
13309.6
  
  10 x  810 1080  270 Q3
 220.64    
 0 270 270 Q4
Resolviendo este sistema de ecuaciones obtenemos:
Q2  12.3511 x10 5 mm
Q3  24.6882 x10 5 mm
Q4  24.6990 x10 5 mm

Y para obtener la reacción en el empotramiento tómanos la siguiente

submatriz:

0 
Q 
2647.68  R1 105 x810  810 0 0 2 
Q3
 
Q4 
Resolviendo obtenemos:
R1   10007.06 N

11
6. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la
siguiente ecuación:

e
 E   Qi 
    1 1  
e

 l  Qi  1
Y obtenemos lo siguiente:

12
  3 x 10   0   5 N
5
1   11  x10  1 0. 7406 2
 50 1  2.1 351  m
3x10   2.1 351  5 N
5
2    11  x10  2 0.1481 2
 250 2  4.2 68 2 m
3x10   4.2 68 2 5 N
5
31   11  x10  3 0. 0 0 98 2
 250 3  4.2 69 0 m

7. RESULTADOS
Finalmente, los resultados son mostrados en la siguiente tabla:

R1   10007.06 N
N
 1  0.074
mm 2
N
 2  0.1481
mm 2

13
 N
 3  0.00000981
mm 2

8. DIAGRAMA DE FLUJO

INICIO

INGRESO DE DATOS
CONSTANTES : E, f, t
VECTORES : L, A, P

CALCULO DE VECTORES

 AL1   EA1 EA1 

  R1   1  1 0 0 
 2   L 1 L 
 AL   AL  
2 1
 EA EA 2 EA1 EA 2
 1  0 
F=
 2 2  ; K=
 L1 L2 L L2 
 AL3 AL2    EA 2
EA EA 2
3
EA 3
   PA   0  2  2  3 
 2 2   L L3 L L 
 AL 
3
  0 EA 3
EA 3 
  0  3
2  L L3 

TRAFORMACION DE ECUACION MATRICIAL

14
   1 
1
 AL1 EA
    0 0 
2  L1 
  2 1
EA 2  R1 
 AL   AL    0 EA  EA
2 1
 2 0  Q 
 2 2 = L2 L1 L   2
 AL3 AL2    EA 2 EA3 EA 2 EA 3
Q3 
   PA  0  2  2  3   
 2 2   L L3 L L  Q4 
 AL3
   EA 3
EA3 
   0 0  3
2  L L3 

IMPRESIÓN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , E1 , E 2 , E 3

FIN
Luego escribimos la siguiente función en MATLAB:

H=input('Ingrese la altura de la placa= ');

B=input('ingrese la base de la placa= ');
pa=input('Ingrese la carga PA= ');
pb=input('Ingrese la carga PB= ');
t=input('Ingrese el espesor de la placa= ');
j=input('Ingrese la densidad del material= ');
E=input('Ingrese el modulo de elasticidad= ');
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;
%MOSTRANDO LOS RESULTADOS
disp('..............................');
disp(' RESULTADOS');

15
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');

9. USO DEL PROGRAMA DE MATLAB

=====1200========
\ PRIMERA /
\ PARTE /
\____600_____/
\SEGUNDA /
\ PARTE /
\__300_/
\ / TERCERA
\ / PARTE
\/

Luego ejecutamos el programa y obtenemos:

>>[E,F,Q,S]=traccion(3*10^11,78400,0.1125,0.05625,0.01875,
0.4,0.2,0.2,0.06,15000)

E =

1.8463e-014

F =

1.0e+004 *

-1.7940
1.7205
0.0588
0.0147

16
Q =

1.0e-006 *

0
0.2126
0.2213
0.2266

S =

1.0e+005 *

1.5947
0.1307
0.0784

Donde S es la matriz de esfuerzos, F es la matriz de fuerzas en los

nodos y Q es la matriz de deformaciones y E es el error al calcular la
reacción en el nodo 1.

CONCLUSIONES

Se puede apreciar que las deformaciones son realmente pequeñas

(décimas de micras), además todas son hacia abajo que es el sentido
positivo asumido como referencia.

Los esfuerzos son positivos, lo que indica esfuerzos de compresión para

nuestro sistema de referencia.

Se puede apreciar que los resultados de MATLAB arrojan un error de

1.846x10-14, prácticamente cero al momento de calcular la reacción en
el nodo 1; mientras que haciendo los cálculos manualmente se obtiene
un error de cero por ciento, con lo cual se podría afirmar que la
aproximación a tres elementos finitos es totalmente exacta.

17