Í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
 PRIMERA PRÁCTICA CALIFICADA
(TRACCION SIMPLE)

ENUNCIADO DEL PROBLEMA

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

Considerar:
PA = 1000 lb-f
t (espesor) = 0.125 in
E = 19x106 psi
Y = 0.28901 lb-f/in3
SOLUCIÓN:
1. MODELADO DEL CUERPO REAL
Se considerarán cuatro elementos finitos. Para facilitar los cálculos
los elementos finitos tendrán longitud de 2.5 in.
Y los espesores lo calculamos tomando el punto medio de cada
elemento finito:
(2 + 1.75)
𝑏1 = = 1.875 𝑖𝑛
2
(1.75 + 1.5)
𝑏2 = = 1.625 𝑖𝑛
2
(1.5 + 1.25)
𝑏3 = = 1.375 𝑖𝑛
2
(1.25 + 1)
𝑏4 = = 1.125 𝑖𝑛
2

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

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 (in) (in2)
1 1 2 1 2 2.5 0.234375
2 2 3 2 3 2.5 0.203125
3 3 4 3 4 2.5 0.171875
4 4 5 4 5 2.5 0.140625
2. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)

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

globales:

Luego el vector de desplazamiento será:

0 
Q 2
 
Q Q3  in 
 
Q 4
Q5 

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:

 y  Axl1
F11  R1  8.467110 2  R1 lb
2
y  Axl1
F21   8.467110 2 lb
2
y  Axl2
F22   7.338110 2 lb
2
y  Axl2
F32   7.338110 2 lb
2
y  Axl3
F33   6.2092 10 2 lb
2
F43 
y  Axl 3  6.2092 10 2 lb
2
F44 
y  Axl 4  5.0803 10 2 lb
2
F54 
y  Axl 4  P  1000.0508 lb
A
2

Ahora analizamos las fuerzas para todo el cuerpo:

F1  F11  8.467110 2  R1 lb
F2  F21  F22 15.8052 10 2 lb
F3  F32  F33 13.5473 10 2 lb
F4  F43  F44 11.2895 10 2 lb
F5  F54  1000.0508 lb

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

 2

 F1  8.467110  R1 
 F 2 15.8052 10 2 
   
F1   F 3   13.5473 10  lb 
2

   2 
 F 4 11.2895 10 
 F 5  1000.0508 
 
4. MATRIZ DE RIGIDEZ

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

está determinada por la siguiente ecuación:

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

Reemplazando para los valores calculados y utilizando la tabla de

conectividad obtenemos:

 1 1 0 0 0 
 1 1 0 0 0 
 0.234375 x19 x10 6   
K i      0 0 0 0 0 
 2.5 1  0 0 0 0 0

 
 0 0 0 0 0 

0 0 0 0 0
0 1 1 0 0 
 0.203125 x19 x10 6  
   0 1 1 0 0 
 2.5 2  
0 0 0 0 0
 0 0 0 0 0 

0 0 0 0 0
0 0 0 0 0 
 0.171875 x19 x10  
6
    0 0 1 1 0 
 2.5 3  0 0 1 1 0


 0 0 0 0 0 

0 0 0 0 0
0 0 0 0 0 
 0.140625 x19 x10  
6
    0 0 0 0 0
 2.5 4  0 0 0 1  1


 0 0 0  1 1 
Finalmente:
 178.25  178.125 0 0 0 
 178.125 332.5  154.375 0 0 
  lb
K i  10 4 x  0  154.375 285  130.625 0 
  in
 0 0  130.625 237.5  106.875
 0 0 0  106.875 106.875 

5. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO

La ecuación de rigidez está determinada por la siguiente ecuación:

Fi  K i  Q

Lo que con nuestros valores calculados tenemos:

8.467110 2  R1 
 2

15.8052 10 
 2 
 13. 5473 10 
11.2895 10 2 
 
1000.0508 
 
 178.125  178.125 0 0 0  0 
 178.125 332 .5  154 .375 0 0  Q 
  2
 10 x 
4
0  154.375 285  130.625 0  Q3 
  
 0 0  130.625 237.5  106.875 Q4 
 0 0 0  106.875 106.875  Q5 

Para obtener los desplazamientos tomamos la siguiente submatriz:

15.8052 10 2   332.5  154.375 0 0  Q2 
    Q 
13.5473 10  10 4 x   154.375  130.625
2
285 0  3
11.2895 10  2   0  130.625 237.5  106.875 Q4 
    
1000.0508   0 0  106.875 106.875  Q5 

Resolviendo este sistema de ecuaciones obtenemos:

 Q2  5.62 x10 4 in
Q3  12.10 x10  4 in
Q4  19.7 x10  4 in
Q5  29.11 x10  4 in
Y para obtener la reacción en el empotramiento tómanos la siguiente
submatriz:
0 
Q 
 2
8.467110 2

 R1 10 x178.25  178.25 0 0 0Q3 
4

 
Q4 
Q5 
 
Resolviendo obtenemos:
R1  1001.84 lb

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

 Q 
e
E
 e     1 1  i

l  Qi  1
Y obtenemos lo siguiente:

 19 x10 6   0
 1     1 1    10 5   1  427.12 psi
 2.5 1 5.62 

 19 x10 6  5.62 
 2     1 1   10 5   2  492.48 psi
 2.5  2 12.10

 19 x 10 6   12.10
 3     1 1    10   3  577.6 psi
5

 2 .5 3 19 .7 

 19 x10 6   19.7 
 4     1 1    10 5   4  715.16 psi
 2.5  4 29.11
7. RESULTADOS
Finalmente, los resultados son mostrados en la siguiente tabla:

R1   1000.6 lb
 1  427.12 psi
 2  492.48 psi
 3  577.6 psi
 4  715.16 psi
 8. DIAGRAMA DE FLUJO

INICIO

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

CALCULO DE VECTORES
 AL 1
  EA1 EA1 
  R 1   1  0 0 0 
 22 1   L 1 L1 
 AL   AL  
2
 EA EA EA1 EA 2
 1  2 0 0 
 2 2   L1 L2 L L 
 AL3 AL2    EA 2 EA3 EA 2 EA3 
F=    ; K=  0  2  2  0 
 24 2   L L3 L L3 
 AL  AL 3
  0 EA 3
EA 4
EA3 EA 4 
 0  3   4
 2 2   L L4 L3 L 
 4  
 AL   P  EA 4 EA 4 
 0 0 0  4 
 2 A
  L L4 

TRAFORMACION DE ECUACION MATRICIAL

 AL 1
  EA1 
    1  0 0 0 
 2 2 1  
L1 
 AL   AL    0 EA  EA
2 1
EA 2  R1 
 2 0 0 
 2 2   L2 L1 L  Q 
 AL3 AL2    EA 2
EA3 EA 2 EA3   2
  = 0  2  2  0  Q3 
 24 2   L L3 L L3  Q 
 AL  AL 3
   EA 3
EA 4
EA3 EA 4   4
 0 0  3   4
 2 2   L L4 L3 L  Q5 
 4  
 AL   P  0 EA 4
EA 4 
0 0  4 
 2 A
  L L4 

IMPRESIÓN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , Q5 ,  1 ,  2 ,  3 ,  4

FIN
Luego escribimos la siguiente función en MATLAB:

clc
clear all;
format long
B=input('longitud de la base mayor =');
b=input('longitud de la base menor =');
espesor=input('espesor de la placa triangular =');
L1=input('longitud elemento 1(in) =');
L2=input('longitud elemento 2(in) =');
L3=input('longitud elemento 3(in) =');
L4=input('longitud elemento 4(in) =');
E=input('modulo de elasticidad(psi) =');
Pa=input('carga puntual(lb-f) =');
y=input('peso especifico(lb-f/in3) =');
B1=(b*L1+B*(L2+L3+L4))/(L1+L2+L3+L4);
B2=(b*(L1+L2)+B*(L3+L4))/(L1+L2+L3+L4);
B3=(b*(L1+L2+L3)+B*(L4))/(L1+L2+L3+L4);
A1=(B+B1)*espesor/2;
A2=(B1+B2)*espesor/2;
A3=(B2+B3)*espesor/2;
A4=(B3+b)*espesor/2;
%matriz de rigidez para el elemento 1
w=zeros(5);
w(1,1)=1;w(1,2)=-1;w(2,1)=-1;w(2,2)=1;
m1=w;
K1=(E*A1/L1)*m1;
%matriz de rigidez para el elemento 2
t=zeros(5);
t(2,2)=1;t(2,3)=-1;t(3,2)=-1;t(3,3)=1;
m2=t;
K2=(E*A2/L2)*m2;
%matriz de rigidez para el elemento 3
x=zeros(5);
x(3,3)=1;x(3,4)=-1;x(4,3)=-1;x(4,4)=1;
m3=x;
K3=(E*A3/L3)*m3;
%matriz de rigidez para el elemento 4
z=zeros(5);
z(4,4)=1;z(4,5)=-1;z(5,4)=-1;z(5,5)=1;
m4=z;
K4=(E*A4/L4)*m4;
%matriz de rigidez global
K=K1+K2+K3+K4;
%fuerzas de cuerpo
f1=(A1*L1*y)/2;
f2=(A2*L2*y)/2;
f3=(A3*L3*y)/2;
f4=(A4*L4*y)/2;
f=[f1;f1+f2;f2+f3;f3+f4;f4+Pa];
ff=[f1+f2;f2+f3;f3+f4;f4+Pa];
KK=zeros(4);
KK(1,1)=K(2,2);KK(1,2)=K(2,3);KK(2,1)=K(3,2);
KK(2,2)=K(3,3);KK(2,3)=K(3,4);KK(3,2)=K(4,3);
KK(3,4)=K(4,5);KK(4,3)=K(5,4);
KK(3,3)=K(4,4);KK(4,4)=K(5,5);
Q=inv(KK)*ff;
Q1=0;
Q2=Q(1,1);
Q3=Q(2,1);
Q4=Q(3,1);
Q5=Q(4,1);
T1=(E/L1)*([-1 1]*[Q1;Q2]);
T2=(E/L2)*([-1 1]*[Q2;Q3]);
T3=(E/L3)*([-1 1]*[Q3;Q4]);
T4=(E/L4)*([-1 1]*[Q4;Q5]);
R1=([K(1,1) K(1,2) K(1,3) K(1,4) K(1,5)]*[Q1;Q2;Q3;Q4;Q5])-f1;
%MOSTRANDO LOS RESULTADOS
disp('..............................');
disp(' RESULTADOS');
disp('============');
disp('*EL VECTOR DESPLAZAMIENTO (mm)');
disp('Q1:');
disp(Q1);
disp('Q2:');
disp(Q2);
disp('Q3:');
disp(Q3);
disp('Q4:');
disp(Q4);
disp('Q5:');
disp(Q5);
disp('*Matriz de rigidez " Kij " (psi)');
disp(K);
disp('*Esfuerzo en 1(psi)');
disp(T1);
disp('*Esfuerzo en 2(psi)');
disp(T2);
disp('*Esfuerzo en 3(psi)');
disp(T3);
disp('*Esfuerzo en 4(psi)');
disp(T4);
disp('*Reaccion 1 (lb)');
disp(R1);

9. USO DEL PROGRAMA DE MATLAB

Luego ejecutamos el programa y obtenemos:

longitud de la base mayor =2

longitud de la base menor =1
espesor de la placa triangular =0.125
longitud elemento 1(in) =2.5
longitud elemento 2(in) =2.5
longitud elemento 3(in) =2.5
longitud elemento 4(in) =2.5
modulo de elasticidad(psi) =19e6
carga puntual(lb-f) =1000
peso especifico(lb-f/in3) =0.28901
..............................
RESULTADOS
============
*EL VECTOR DESPLAZAMIENTO (mm)
Q1:
0
Q2:
5.616601952850878e-04

Q3:
0.001209627269298

Q4:
0.001975302826864

Q5:
0.002911022876023

*Matriz de rigidez " Kij " (psi)

1781250 -1781250 0 0 0
-1781250 3325000 -1543750 0 0
0 -1543750 2850000 -1306250 0
0 0 -1306250 2375000 -1068750
0 0 0 -1068750 1068750

*Esfuerzo en 1(psi)
4.268617484166667e+03

*Esfuerzo en 2(psi)
4.924549762500000e+03

*Esfuerzo en 3(psi)
5.819134237500000e+03

*Esfuerzo en 4(psi)
7.111472373611112e+03

*Reaccion 1 (lb)
-1.000541893750000e+03

>>
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 tracció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.