UNIVERSIDAD NACIONAL

DE INGENIERÍA
Facultad de Ingeniería Mecánica
TRACCIÓN SIMPLE

Primera Práctica Calificada

PROFESOR : Ing. Jorge Vera Ermitaño

CURSO : Cálculo por Elementos Finitos

CODIGO DEL CURSO : MC516
SECCIÓN : ” H”

ALUMNO : Sevillano Armas José Leoncio

CODIGO DEL ALUMNO : 19781269I
No DE ORDEN ( #) : 9

2010
 INTRODUCCIÓN

El MATLAB es un programa para realizar cálculos numéricos con vectores

y matrices. En este nuestro caso, las capacidades de cálculo científico y
técnico de este programa se ven utilizadas en una amplia posibilidad de
problemas que plantean los cálculos de la resistencia de materiales y el
diseño de ingeniería.

Iniciamos nuestro recorrido del Cálculo por Elementos Finitos con la

relación más elemental, el caso de la tracción pura, que, sin embargo, es
simplemente el inicio de la amplia gama de problemas que trataremos en el
curso.

CÁLCULO POR ELEMENTOS FINITOS

PRIMERA PRÁCTICA CALIFICADA

TRACCIÓN SIMPLE

1. ENUNCIADO DEL PROBLEMA

En la figura se tiene una barra en la cual se aplica una carga

P= ( 200+5#)x103 N.
Determine:
a) Desplazamientos nodales.
b) Esfuerzos en cada material.
c) Fuerzas de reacción.

(300+#/2) mm (400+#/3)mm

1 2 P 3

1 2
Datos según el numero de orden # = 9

P = ( 200 + 5x#)x103 N = ( 200 + 5x9)x103 N = 245000 N

ALUMINIO
A1 = 2400 mm2
E1 = 70x109 N/m2 = 70x103 N/mm2
L1 =300+#/2 =300+9/2 = 304.5 mm

ACERO
A2 = 600 mm2
E2 = 200x109 N/m2 = 200x103 N/mm2
L2 =400+#/3 =400+9/3 = 403 mm

2. MODELADO DEL CUERPO REAL

Se considerarán 2 elementos finitos. Para facilitar los cálculos , los

elementos finitos tendrán longitud de 200 mm, 104.5 mm y 403 mm.

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

(300+#/2) mm (400+#/3)mm

1 2 P 3

2
1
Y las áreas se calculan de la siguiente relación:

Ai = bi x t

Así, el Cuadro de Conectividad queda constituido de la siguiente manera:

Nodos Grados de
le Ae
e Libertad (GDL)
(1) (2) 1 2 (mm) (mm2)
1 1 2 1 2 304.5 240
2 2 3 2 3 403 600

3. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)

Observando el gráfico anterior, el vector desplazamiento será:

0 
Q = Q  (m ) m
 2
 0 
Donde Q1=0, Q3=0 pues la placa está empotrada y los demás
desplazamientos son incógnitas que tendrán que ser calculadas.

4. VECTOR CARGA

Analizando las fuerzas en cada elemento finito:

γ ( A x l )1
F11 = + R1 ( N ) = 1544 .48 + R N
2

γ ( A x l )1
F21 = ( N ) = 1544 .48 N
2
 γ ( A x l )2
F22 = ( N ) = 1103 .20 N
2

γ ( A x l )2
F32 = + P ( N ) = 1103 .20 + P N
2

γ ( A x l )3
F33 = ( N ) = 882 .56 N
2

γ ( A x l )3
F43 = ( N ) = 882 .56 N
2

Ahora, analizamos las fuerzas para todo el cuerpo:

F1 = F11 = 1544 .48 + R N

F2 = F21 + F22 = 2647 .68 N

F3 = F32 + F33 = 1985 .76 + P N = 11985 .76 N

γ ( A x l )3
F4 = F43 = ( N ) = 882 .56 N
2

Entonces, el vector carga se expresa de la siguiente manera:

 F1  1544 .48 + R 
F   2647 .68 
F1 =  2
=   (N )
 F3   11985 .76 
   
 F4   882.56 

5. 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 0 0
 0 0 0 0 0
 A E  −1 1 0 +  AE  1 −1 +  AE  0 0 
K =     
 l 1  0 0 0 0  l 2 0 −1 1 0  l 3 0 0 1 −1
     
0 0 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
−1 0 0 −1 0
157500 x 3 x 10 5
  1 0 + 112500 x 3 x 10 5
  1 
K =
 
 
 

 250 1 0 0 0 0  250 2 0 −1 1 0
   
0 0 0 0 0 0 0 0

0 0 0 0
0 0
 45000 x 3 x 10 5
  0 0 
+ 
 

 500 3 0 0 1 −1
 
0 0 −1 1

 1884 −1884 0 0 
−1884 2030 −1350 0 
K = 10 5 x   N / mm
 0 −1350 1620 − 270 
 
 0 0 − 270 270 

6. ECUACIÓN DE RIGIDEZ Y CONDICIONES DE CONTORNO

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

Fi = K ij Q j

Aplicando nuestros valores calculados obtenemos:

1544 .48 + R   1884 − 1884 0 0  0

 2647 .68 


 = 10 5 x − 1884 2030 − 1350 0  Q2 
 11985 .76   0 − 1350 1620 − 270  Q3 
     
 882 .56   0 0 − 270 270  Q4 

Para obtener los desplazamientos tomamos la siguiente submatriz:

 2647.68   2030 − 1350 0  Q2 

11985 .76 = 105 x − 1350 1620 − 270 Q 
     3
 882.56   0 − 270 270  Q4 
Resolviendo este sistema de ecuaciones, se obtiene:

Q2 = 22.8179 x 10-5 mm

Q3 = 32.3502 x 10-5 mm

Q4 = 35.6197 x 10-5 mm

En MATLAB, la expresión sería:

>>A=[2647.68; 11985.76; 882.76]

A=

1.0e+004 *

0.2648
1.1986
0.0883

>> B=[2030 -1350 0; -1350 1620 -270; 0 -270 270]

B=

2030 -1350 0
-1350 1620 -270
0 -270 270

>> inv(B)*A

ans =

22.8179
32.3502
35.6197

Que debemos multiplicar por el factor 10-5.

Y, para obtener la reacción en el empotramiento. Tomamos la siguiente

submatriz:
 0
Q 
[1544 .48 + R ] = 10 5 x [1884 − 1884 0 0]  2 
Q3 
 
Q4 

Resolviendo, se obtiene:

R = –44533 N

En MATLAB, esto se podría expresar de la siguiente manera:

>> A1=[1884 -1884 0 0]

A1 =

Columns 1 through 3

1884 -1884 0

Column 4

>> Q=[0; 22.8179; 32.3502; 35.6197]

Q=

0
22.8179
32.3502
35.6197

>> A1*Q -1544.48

ans =

-4.4533e+004

>>
Respuesta parcial a la que debe añadirse el factor 10–5 de los valores de Qi
que se neutralizan con el factor 105 de la ecuación matricial. Finalmente,
teniendo en cuenta el término independiente del miembro izquierdo, se
halla el valor de R.

7. ESFUERZOS

Para calcular los valores de los esfuerzos por elemento, aplicamos la

siguiente ecuación:
e
E  Qi 
σ e
=   [ − 1 1] Q 
l   i +1 

Y, se obtiene lo siguiente:

 3 x 10 5   0 
σ1 =   [ − 1 1]   x 10
−5
→ σ1 = 0.274 N / mm 2
 250 1 22.8179 

En MATLAB:

>> U=[-1 1]

U=

-1 1

>> Q1=[0; 22.8179]

Q1 =

0
22.8179

>> U*Q1

ans =

22.8179

>>

Resultado que debe aplicarse con los factores numéricos expresados.

  3 x 10 5  22.8179 
σ 2 =   [ − 1 1]   x 10
−5
→ σ 2 = 0.1144 N / mm 2
 250 2 32.3502 

En MATLAB:

>> U

U=

-1 1

>> Q2=[22.8179; 32.3502]

Q2 =

22.8179
32.3502

>> U*Q2

ans =

9.5323

>>

Resultado que debe aplicarse con los factores numéricos expresados.

 3 x 10 5  32.3502 
σ 3 =   [ − 1 1]   x 10
−5
→ σ 3 = 0.02 N / mm 2
 500 3 35.6197 

En MATLAB:

>> U

U=

-1 1

>> Q3=[32.3502; 35.6197]

Q3 =
 32.3502
35.6197

>> U*Q3

ans =

3.2695

>>

Resultado que, finalmente, también debe aplicarse con los factores

numéricos expresados.

8. RESULTADOS

Finalmente, los resultados son los siguientes:

R = –44533 N

σ 1 = 0.274 N/mm2

σ 2 = 0.1144 N/mm2

σ 3 = 0.02 N/mm2
 CONCLUSIONES

En primer lugar, el MatLab demostró ser una poderosa herramienta de

cálculo de las operaciones matriciales que fueron requeridas en varias
oportunidades durante el desarrollo de la solución del problema planteado.

Este programa demostró que puede proporcionar un ahorro enorme en el

tiempo de cálculo, sin los requisitos, en algunos casos tediosos, de los
lenguajes de programación tradicionales.

Así, la relativa simplicidad del aprendizaje del MatLab sumada a su

excelente rendimiento hace que se constituya en la herramienta de cálculo
en la ingeniería de nuestros tiempos.

La resistencia de materiales se ve potenciada en sus capacidades por la

presencia de esta herramienta de cálculo y método utilizado, el de
elementos finitos.
DIAGRAMA DE FLUJO

INICIO

SE EVALÚAN LOS DATOS DEL

PROBLEMA

MODELADO DEL CUERPO REAL

(Se establecen las dimensiones de cada
elemento finito y el Cuadro de
Conectividad)

SE ESTABLECEN LOS GRADOS DE

LIBERTAD NODALES Y EL
VECTOR DESPLAZAMIENTO

SE ESTABLECE EL VECTOR
CARGA

SE CONSTRUYE LA MATRIZ DE
RIGIDEZ

SE ESTABLECEN LAS
CONDICIONES DE CONTORNO Y
LA ECUACIÓN DE RIGIDEZ

SE CALCULAN LOS ESFUERZOS

EN LA ESTRUCTURA

FIN