UNIVERSIDAD NACIONAL DE

INGENIERÍA
FACULTAD DE INGENIERÍA
MECÁNICA

Curso: CÁLCULO POR ELEMENTOS FINITOS

Profesor: Ing. ABREGU LEANDRO EDWIN
Tema: TENSION UNIDIMENSIONAL
Fecha de entrega: y evaluación 4 de Abril del 2017
PRIMERA PRÁCTICA CALIFICADA
 ÍNDICE

1. INTRODUCCIÓN………………………………………………………………..1

2. ENUNCIADO DEL PROBLEMA……………………………………………….2

3. PROCEDIMIENTOS…………………………………………………………….2

4. RESULTADOS…………………………………………………………………9

5. DIAGRAMA DE FLUJO………………………………………………………..10

6. LENGUAJE EN MATLAB……………………………………………………..11

7. CONCLUSIÓN…………………………………………………………………..14

8. BIBLIOGRAFÍA………………………………………………………………....15
INTRODUCCIÓN

En el presente trabajo se moldeará una placa de forma trapezoidal por

medio de elementos finitos, se hallará sus esfuerzos, su matriz de rigidez
y reacciones, esta placa está sometida a una carga puntual en el medio de
ella.
El problema se planteará en 6 puntos o enunciados de preguntas, las
cuales se describen más abajo en el transcurso del informe.
Para este caso moldearemos la placa para 4 elementos finitos, a la vez he
realizado un programa en Matlab el cual simula correctamente la situación
del problema, pero lo realice con mayor campo no solo para 4 elementos
sino para “n” elementos.
Al finalizar realizare algunas conclusiones, espero que el informe cumpla
con sus expectativas.

1
 ENUNCIADO DEL PROBLEMA

Considere la placa delgada (acero) en la figura presente. La placa tiene un espesor

uniforme t=0.125 in.

Módulo de Young 𝐸 = 30 × 106 psi y densidad por peso 𝜌 = 0.2836 𝑙𝑏/𝑖𝑛3 .

Además de su propio peso, la placa está sometida a una carga concentrada P=100
lb en su punto medio.

Cuestiones a resolver:

(a) Modelar la placa con cuatro elementos finitos

(b) Escriba las expresiones para las matrices de rigidez de los elementos y los
vectores fuerza
(c) Sume las matrices de rigidez estructural de cada elemento y halle también el
vector fuerza de la placa
(d) Resuelva el vector desplazamiento Q
(e) Evalué los esfuerzos en cada elemento
(f) Determine la fuerza de reacción en el soporte

2
LEYENDA:

A1 = 2 in

A2 = 1 in

L = 10 in

t (espesor) = 0.125 in

E = 10.4 × 106 psi

𝛾 = 6.29 𝑙𝑏/𝑖𝑛3

SOLUCIÓN:

1. MODELADO DEL CUERPO REAL

Se considerarán 4 elementos finitos. Para facilitar los cálculos los elementos finitos tendrán
longitud de 2.5” cada uno.

Y los espesores lo calculamos tomando el punto medio de cada elemento finito:

b1 
2  1.75  1.875"
2
b2 
1.75  1.5 1.625"
2
(1.5  1.25)
b3   1.375"
2
(1.25  1)
b4   1.125"
2

3
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 (pulg) (pulg2)

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

4
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  pu lg 
 
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.

5
 3. VECTOR CARGA

Analizando las fuerzas en cada elemento finito:

𝑦(𝐴𝑥𝑙)1
𝐹11 = + 𝑅1 = 1.8427 + 𝑅1 𝑙𝑏𝑓
2

𝑦(𝐴𝑥𝑙)1
𝐹21 = = 1.8427
2

𝑦(𝐴𝑥𝑙)2
𝐹22 = = 1.7
2

6
 𝑦(𝐴𝑥𝑙)2
𝐹32 = + 𝑃 = 1001.7
2

𝑦(𝐴𝑥𝑙)3
𝐹33 = = 1.45
2

𝑦(𝐴𝑥𝑙)3
𝐹43 = = 1.45
2

𝑦(𝐴𝑥𝑙)4
𝐹44 = = 1.25
2

𝑦(𝐴𝑥𝑙)4
𝐹54 = = 1.25
2

𝑦(𝐴𝑥𝑙)5
𝐹55 = = 1.1
2

Ahora analizamos las fuerzas para todo el cuerpo:

𝐹1 = 𝐹11 = 𝑅1 + 1.8427 𝑙𝑏𝑓

𝐹1 = 𝐹12 + 𝐹22 = 1003.4𝑙𝑏𝑓

𝐹1 = 𝐹23 + 𝐹33 = 2.9 𝑙𝑏𝑓

𝐹1 = 𝐹34 + 𝐹44 = 2.5𝑙𝑏𝑓

𝐹1 = 𝐹45 = 1.1𝑙𝑏𝑓

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

𝐹1 𝑅1 + 1.84
𝐹2 1003.4
𝐹 = 𝐹3 = 2.9 𝑙𝑏𝑓
𝐹4 2.5
[𝐹5 ] [ 1.1 ]

7
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 0 0 0 0 0 
 AE    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 0 0 0 0 0 0 0  1 1 

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

 1 1 0 0 0 0 0 0 0 0
 1 1  0 1  1
0 0 0 0 0
 0.234375 x10.4 x106    0.203125 x10.4 x106  
K i      0 0 0 0 0  K i     0  1 1 0 0
 2.5 1  0 0 0 0 0
  2.5 1 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 
 0.171875 x10.4 x106    0.140625 x10.4 x106  
 K i     0 0 1 1 0  K i     0 0 0 0 0
 2.5 1 0 0 1 1

0  2.5 1 0 0 0 1  1

 
0 0 0 0 0 0 0 0  1 1 

8
Finalmente:

 9.75  9.75 0 0 0 
  9.75 18.2  8.45 0 0 
 lb
K i  10 5 x  0  8.45 15.6  7.15 0 
  pu lg
 0 0  7.15 13  5.85
 0 0 0  5.85 5.85 

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:

 9.75  9.75 0 0 0  0

𝑅1 + 1.84   9.75 18.2  8.45 0 0  Q 2
1003.4   
2.9 10 5 x  0  8.45 15.6  7.15 0   Q3 
2.5    
 0 0  7.15 13  5.85 Q 4
[ 1.1 ]
 0 0 0  5.85 5.85  Q5

Para obtener los desplazamientos tomamos la siguiente submatriz:

 18.2  8.45 02 0  Q 
1003.4   8.45 15.6  7.15  Q 
0  3
[ 2.9 ] 10 5 x 
2.5  0  7.15 13  5.85 Q4 
1.1   
 0 0  5.85 5.85  Q5 

9
Resolviendo este sistema de ecuaciones obtenemos:

𝑄1 = 0 𝑝𝑢𝑙𝑔

𝑄2 = 0.0010 𝑝𝑢𝑙𝑔

𝑄3 = 0.0010 𝑝𝑢𝑙𝑔

𝑄4 = 0.0010 𝑝𝑢𝑙𝑔

𝑄5 = 0.0011 𝑝𝑢𝑙𝑔

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

0
Q 2
 
R1 10 5 x9.75  9.75 0 0 0Q3
 
Q 4
Q5

Resolviendo obtenemos:

𝑅1 = −1011.79 𝑙𝑏𝑓

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

 1 1 
E
e
Qi 
  
e

l  Qi  1

Y obtenemos lo siguiente:

10.4𝑥106 0 𝑙𝑏
𝜎1 = ( ) [−1 1] [ ] = 4309.12
2.5
1
0.0010 𝑝𝑢𝑙𝑔2

10.4𝑥106 0.0010 𝑙𝑏
𝜎2 = ( ) [−1 1] [ ] = 32.04
2.5
2
0.0010 𝑝𝑢𝑙𝑔2

10.4𝑥106 0.0010 𝑙𝑏
𝜎3 = ( ) [−1 1] [ ] = 20.72
2.5
3
0.0010 𝑝𝑢𝑙𝑔2

10.4𝑥106 0.0010 𝑙𝑏
𝜎4 = ( ) [−1 1] [ ] = 7.86
2.5
4
0.0011 𝑝𝑢𝑙𝑔2

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

𝑅 = 1011.79 𝑙𝑏𝑓
𝑙𝑏
𝜎1 = 4309.12
𝑝𝑢𝑙𝑔2
𝑙𝑏
𝜎2 = 32.04
𝑝𝑢𝑙𝑔2
𝑙𝑏
𝜎3 = 20.72
𝑝𝑢𝑙𝑔2
𝑙𝑏
𝜎4 = 7.86
𝑝𝑢𝑙𝑔2

11
DIAGRAMA DE FLUJO
INICIO

INGRESO DE DATOS

CONSTANTES : E, f, t

VECTORES: L, A, P

CALCULO DE VECTORES

 AL1   EA1 EA1 

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

   PA   0  2   3 
 2 2   L L3 L2 L 
 AL3   0 EA 3
EA3 
0  3
 2   L L3 

TRAFORMACION DE ECUACION MATRICIAL

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

12
 IMPRESIÓN DE RESULTADOS
R1 , Q2 , Q3 , Q4 , E1 , E2 , E3

FIN

LENGUAJE EN MATLAB

clc

clear all

%datos del problema

h=input('Ingrese la altura de la placa en pulgadas= ');

b0=input('ingrese la base superior de la placa en pulgadas= ');

bn=input('ingrese la base inferior de la placa en pulgadas= ');

Pa=input('Ingrese la carga concentrada en lbf= ');

t=input('Ingrese el espesor de la placa en pulgadas= ');

y=input('Ingrese la densidad del material en lb/pulg^3= ');

E=input('Ingrese el modulo de elasticidad en lbf/pulg^2= ');

n=4;

%calculo de bases y áreas de elementos

le=zeros(n,1); ho=zeros(n,1); bo=zeros(n,1); b=zeros(n,1); a=zeros(n,1);

Fe=zeros(n+1,1);

bo(1)=b0; ho(1)=h;

for i=1:n

if n>i

le(i)=input('Ingrese longitud del elemento finito(pulg):');

b(i)=(bo(i)+bn+(bo(i)-bn)*(ho(i)-le(i))/ho(i))/2;

a(i)=b(i)*t;

ho(i+1)=ho(i)-le(i);

bo(i+1)=2*b(i)-bo(i);

else

le(i)=ho(i);

13
 b(i)=(bn+bo(i))/2;

a(i)=b(i)*t;

end

end

disp('Bases(pulg):')

disp(b')

disp('Longitudes(pulg):')

disp(le')

disp('Areas(pulg^2):')

disp(a')

%calculo de las fuerzas

for i=1:n

Fe(i)=y*a(i)*le(i)/2;

end

for i=1:n+1

if i==1

F(i)=Fe(i);

elseif i==n+1

F(i)=Fe(i-1);

else

F(i)=Fe(i-1)+Fe(i);

end

end

F(2)=F(2)+Pa;

disp('El vector de fuerzas(psi):')

disp(F')

%calculo de la matriz rigidez

k=zeros(n+1);

14
for i=1:n

x=zeros(n+1);

x(i,i)=1;x(i+1,i)=-1;x(i,i+1)=-1;x(i+1,i+1)=1;

k=k+(a(i)*E/(le(i)))*x;

end

disp('La matriz de rigidez es(lbf/pulg):')

disp(k)

%calculo de desplazamientos

inv(k(2:n+1,2:n+1));

((F(2:n+1))');

Q=inv(k(2:n+1,2:n+1))*((F(2:n+1))');

Q=[0;Q];

disp('Los desplazamientos de los nodos son(pulg):')

disp(Q)

%calculo de la reaccion

k(1,:)*Q;

R1=k(1,:)*Q-F(1);

disp('La reaccion en el extremo es:')

disp(R1)

%calculo de esfuerzos

for i=1:n

e(i)=(E/(le(i)))*[-1 1]*[Q(i); Q(i+1)];

end

disp('Los valores de los esfuerzos son(lbf/pulg^2):')

disp(e');

15
Ingreso de datos

Resultados parciales

16
Resultados finales

17
CONCLUSIÓN:

Podemos comparar los resultados realizados manualmente y con Matlab, el error

es mínimo, podemos ver la magnitud de los esfuerzos como varía en forma
descendente entre los elementos 1 al 4.

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.

18
Los esfuerzos son positivos, lo que indica esfuerzos de compresión para nuestro
sistema de referencia.

19
BIBLIOGRAFÍA
 Teoría brindada en las clases del Ing.Abregú Leandro de la sección F, los
días lunes de 8 a 10
 Chandrupatla, Tirupathi (2002). “Introduction to finite elements in
engineering”. 3era edición. 59-68.
 Moaveni, Saeed (2008). “Finite element analysis: theory and application with
ANSYS”.3era edición. 8-19.
 http://www.iit.upcomillas.es/~carnicero/Resistencia/Introduccion_al_MEF.pdf
 http://www.uhu.es/javier.pajon/apuntes/finitos.pdf

20