UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA MECÁNICA

CÁLCULO POR ELEMENTOS FINITOS

MC516

“TRACCIÓN SIMPLE y DEFORMACIÓN

TERMICA”

Docente: (Informe 2)
Abregú Leandro, Edwin Asención

Alumnos:

Castro Camarena Diego 20144528G

Inga Vega, Morphi Jeisson 20152073E

Moran Salvatierra, Victor Angel 20141268D

Tello Vargas Alvaro André 20130183B

SECCIÓN: F

2018 - Il
 INDICE

ENUNCIADO DEL PROBLEMA.................................................................................................... 3

SOLUCIÓN DEL PROBLEMA....................................................................................................... 4
1. Modelado del cuerpo real .................................................................................................... 4
2. Formación de la matriz rigidez global .............................................................................. 5
3. Vector de carga global.......................................................................................................... 8
4. Ecuación de equilibro ........................................................................................................... 8
5. Calculo de esfuerzos con variación de temperatura .................................................... 9
6. Resultados ................................................................................Error! Bookmark not defined.
SOLUCIÓN POR MATLAB .......................................................................................................... 11
a) Diagrama de flujo (operación del proceso) .................................................................. 11
b) Código .................................................................................................................................... 11
SOLUCION USANDO ANSYS ........................................................Error! Bookmark not defined.
CONCLUSIONES........................................................................................................................... 12
ENUNCIADO DEL PROBLEMA

Una placa de acero está sujeta a una carga axial, como se muestra en la figura
[Link] las deflexiones y esfuerzos a lo largo de la placa, así mismo calculas el
factor de seguridad (F.S). La placa tiene un espesor de 1/16 pulgadas y un módulo
de elasticidad 𝐄 = 𝟐𝟗 𝐱𝟏𝟎𝟔 𝐥𝐛/𝐢𝐧𝟐 .

Moldear este problema usando 4 nodos y cuatro elementos, como se muestra en la

figura 1.

Figura 1. Placa de Acero

SOLUCIÓN DEL PROBLEMA

1. Modelado del cuerpo real

La placa de acero será dividida usando cuatro nodos y cuatro elementos, como se
muestra en la Figura 2 la cual le asumiremos su comportamiento como un sistema
de resortes.

Figura 1. Subdivisión de la placa de acero

 Figura 2. Comportamiento de Resortes

2. Formación de la matriz rigidez global

-La placa El coeficiente de rigidez está definido por la siguiente ecuación:

(𝑨𝒆 )𝑬
𝒌𝒆 =
𝒍𝒆
De donde:

𝑨𝒆 : Área del elemento “e”.

𝑬: Módulo de elasticidad del acero.

𝒍𝒆 : Longitud del elemento “e”.

-Hallamos el área transversal de cada elemento, para ello empleamos la siguiente

ecuación:

𝑨=𝒃∗𝒕

Donde:

𝒃: Base de cada elemento.

𝒕: Espesor de la placa.
-Construimos un cuadro conectividad relacionando todos los datos:
CUADRO DE CONECTIVIDAD

ELEMENTO 𝒍𝒆 𝒃 𝑨 (𝑨𝒆 )𝑬
𝒌𝒆 =
𝒍𝒆
1 1 5 0.3125 𝟗𝟎𝟔𝟐𝟓𝟎𝟎 𝐥𝐛/𝐩𝐮𝐥𝐠

2 4 2 0.125 𝟗𝟎𝟔𝟐𝟓𝟎𝟎 𝐥𝐛/𝐩𝐮𝐥𝐠

3 4 2 0.125 𝟗𝟎𝟔𝟐𝟓𝟎𝟎 𝐥𝐛/𝐩𝐮𝐥𝐠

4 2 5 0.125 𝟒𝟓𝟑𝟏𝟐𝟓𝟎 𝐥𝐛/𝐩𝐮𝐥𝐠

-Con los datos obtenidos construimos la matriz de rigidez para cada elemento:

𝟏 −𝟏 𝟗𝟎𝟔𝟐𝟓𝟎𝟎 −𝟗𝟎𝟔𝟐𝟓𝟎𝟎
𝑲𝟏 = 𝒌𝟏 [ ]=[ ]
−𝟏 𝟏 −𝟗𝟎𝟔𝟐𝟓𝟎𝟎 𝟗𝟎𝟔𝟐𝟓𝟎𝟎

𝟏 −𝟏 𝟗𝟎𝟔𝟐𝟓𝟎 −𝟗𝟎𝟔𝟐𝟓𝟎
𝑲𝟐 = 𝒌𝟐 [ ]=[ ]
−𝟏 𝟏 −𝟗𝟎𝟔𝟐𝟓𝟎 𝟗𝟎𝟔𝟐𝟓𝟎

𝟏 −𝟏 𝟗𝟎𝟔𝟐𝟓𝟎 −𝟗𝟎𝟔𝟐𝟓𝟎
𝑲𝟑 = 𝒌𝟑 [ ]=[ ]
−𝟏 𝟏 −𝟗𝟎𝟔𝟐𝟓𝟎 𝟗𝟎𝟔𝟐𝟓𝟎

𝟏 −𝟏 𝟒𝟓𝟑𝟏𝟐𝟓𝟎 −𝟒𝟓𝟑𝟏𝟐𝟓𝟎
𝑲𝟒 = 𝒌𝟒 [ ]=[ ]
−𝟏 𝟏 −𝟒𝟓𝟑𝟏𝟐𝟓𝟎 𝟒𝟓𝟑𝟏𝟐𝟓𝟎
-Con los datos obtenidos construimos la matriz de rigidez global:

𝒌𝟏 −𝒌𝟏 𝟎 𝟎 𝟎 𝟎 𝟎 𝟎
−𝒌 𝒌𝟏 𝟎 𝟎] 𝟎 𝒌𝟐 −𝒌𝟐 𝟎
𝑲𝟏𝑮 = [ 𝟏 𝑲𝟐𝑮 = [ ]
𝟎 𝟎 𝟎 𝟎 𝟎 −𝒌𝟐 𝒌𝟐 𝟎
𝟎 𝟎 𝟎 𝟎 𝟎 𝟎 𝟎 𝟎

𝟎 𝟎 𝟎 𝟎 𝟎 𝟎 𝟎 𝟎
𝟎 𝒌𝟑 −𝒌𝟑 𝟎 𝟎 𝟎 𝟎 𝟎
𝑲𝟑𝑮 = [ ] 𝑲𝟒𝑮 = [ 𝒌𝟒 −𝒌𝟒 ]
𝟎 −𝒌𝟑 𝒌𝟑 𝟎 𝟎 𝟎
𝟎 𝟎 𝟎 𝟎 𝟎 𝟎 −𝒌𝟒 𝒌𝟒

-Finalmente formamos la matriz de rigidez global

𝟒
𝒌𝟏 −𝒌𝟏 𝟎 𝟎
−𝒌𝟏 𝒌𝟏 + 𝒌𝟐 + 𝒌𝟑 −𝒌𝟐 − 𝒌𝟑 𝟎
𝑲𝑮𝒍𝒐𝒃𝒂𝒍 = ∑ 𝑲𝒆 = [ ]
𝟎 −𝒌𝟐 − 𝒌𝟑 𝒌𝟐 + 𝒌𝟑 + 𝒌𝟒 −𝒌𝟒
𝟏
𝟎 𝟎 −𝒌𝟒 𝒌𝟒

𝟗, 𝟎𝟔𝟐, 𝟓𝟎𝟎 −𝟗, 𝟎𝟔𝟐, 𝟓𝟎𝟎 𝟎 𝟎

−𝟗, 𝟎𝟔𝟐, 𝟓𝟎𝟎 𝟏𝟎, 𝟖𝟕𝟓, 𝟎𝟎𝟎 −𝟏, 𝟖𝟏𝟐, 𝟓𝟎𝟎 𝟎
[ ]
𝟎 −𝟏, 𝟖𝟏𝟐, 𝟓𝟎𝟎 𝟔, 𝟑𝟒𝟑, 𝟕𝟓𝟎 −𝟒, 𝟓𝟑𝟏, 𝟐𝟓𝟎
𝟎 𝟎 −𝟒, 𝟓𝟑𝟏, 𝟐𝟓𝟎 𝟒, 𝟓𝟑𝟏, 𝟐𝟓𝟎
3. Vector de carga global

𝑭𝟏 𝑹
𝑭𝟐
𝑭=[ ]=[ 𝟎 ]
𝑭𝟑 𝟎
𝑭𝟒 𝟖𝟎𝟎

4. Ecuación de equilibro

[𝑭] = [𝑲𝑮𝒍𝒐𝒃𝒂𝒍 ] × [𝑸]

Reemplazamos:

𝟗, 𝟎𝟔𝟐, 𝟓𝟎𝟎 −𝟗, 𝟎𝟔𝟐, 𝟓𝟎𝟎 𝟎 𝟎 𝑸𝟏 𝑹

−𝟗, 𝟎𝟔𝟐, 𝟓𝟎𝟎 𝟏𝟎, 𝟖𝟕𝟓, 𝟎𝟎𝟎 −𝟏, 𝟖𝟏𝟐, 𝟓𝟎𝟎 𝟎 𝑸𝟐
[ ]×[ ]= [ 𝟎 ]
𝟎 −𝟏, 𝟖𝟏𝟐, 𝟓𝟎𝟎 𝟔, 𝟑𝟒𝟑, 𝟕𝟓𝟎 −𝟒, 𝟓𝟑𝟏, 𝟐𝟓𝟎 𝑸𝟑 𝟎
𝟎 𝟎 −𝟒, 𝟓𝟑𝟏, 𝟐𝟓𝟎 𝟒, 𝟓𝟑𝟏, 𝟐𝟓𝟎 𝑸𝟒 𝟖𝟎𝟎

Aplicando condición de frontera 𝑄1 = 0, y finalmente:

𝑄2 = 0.105 𝑚𝑚

𝑄3 = 0.2098 𝑚𝑚

𝑄4 = 0.3146 𝑚𝑚

𝑄5 = 0.4194 𝑚𝑚

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

−398521.96 + 𝑅1 = −38.03625𝑥105 𝑥0.105

Resolviendo obtenemos:
 𝑅1 = −819.29 𝑁

5. Calculo de esfuerzos

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

𝐸 𝑒 𝑄
𝜎 𝑒 = ( 𝑙 ) [−1 1] [ 𝑖 ] − (𝐸𝛼)𝑒 ∆𝑇…………(5)
𝑄𝑖+1

Donde:

σ: Esfuerzo uniaxial

E: Módulo de Young(N/mm2)

L: espesor del elemento(mm)

α: Coeficiente de dilatación lineal (°C-1)

Δt: Variación de la temperatura(°C)

(𝐸𝛼)𝑒 ∆𝑇 es una constante pues los elementos finitos son del mismo material.

𝑁
(𝐸𝛼)𝑒 ∆𝑇 = 2.07x105 𝑥11x10−6 𝑥150 = 341.55
𝑚𝑚2

Usando la ecuación número 5 para calcular los esfuerzos

2.07𝑥105 0 𝑁
𝜎1 = ( ) [−1 1] [ ] − 341.55 = 0.73346
63.5 1
0.105 𝑚𝑚2

2.07𝑥105 0.105 𝑁
𝜎2 = ( ) [−1 1] [ ] − 341.55 = 0.0815
63.5 2
0.2098 𝑚𝑚2

2.07𝑥105 0.2098 𝑁
𝜎3 = ( ) [−1 1] [ ] − 341.55 = 0.0815
63.5 1
0.3146 𝑚𝑚2
 2.07𝑥105 0.3146 𝑁
𝜎4 = ( ) [−1 1] [ ] − 341.55 = 0.08145
63.5 2
0.4194 𝑚𝑚2

2.07𝑥105 0.3146 𝑁
𝜎4 = ( ) [−1 1] [ ] − 341.55 = 0.08145
63.5 2
0.4194 𝑚𝑚2

7. Resultados

En esta tabla se resume todo lo importante hallado en el proceso de los

cálculos

Deformaciones Esfuerzos en cada

Fuerzas nodales (N)
nodales (mm) elemento (N/mm2)

F1 =-399343.6 y R1 = - 819.29 Q1 = 0 𝜎1 = 0.73346

F2 = 52932.77 Q2 = 0.105 𝜎2 = 0.0815

F3 = 53795.23 Q3 = 0.2098 𝜎3 = 0.0815

F4 = 53140.45 Q4 = 0.3146 𝜎4 = 0.08145

F5 = 239116.60 Q5 = 0.4194
SOLUCIÓN POR MATLAB
a) Diagrama de flujo (operación del proceso)
INICIO

𝑙𝑒𝑒𝑟: 𝑙1, 𝑙2, 𝑙3, 𝑙4;

𝐴1, 𝐴2, 𝐴3, 𝐴4, 𝐸, 𝐹

𝑘1 = 𝐴1 ∗ 𝐸/𝑙1

𝑘2 = 𝐴2 ∗ 𝐸/𝑙2

𝑘3 = 𝐴3 ∗ 𝐸/𝑙3

𝑘4 = 𝐴4 ∗ 𝐸/𝑙4

ൣ𝐾1𝐺 ൧; ൣ𝐾2𝐺 ൧; ൣ𝐾3𝐺 ൧; [𝐾4𝐺 ]

[𝑭] = [𝑲𝑮𝒍𝒐𝒃𝒂𝒍 ] × [𝑸]

𝑸𝟐
൥𝑸𝟑 ൩ ; 𝑹
𝑸𝟒

(𝑸𝒆+𝟏 − 𝑸𝒆 )
𝝈𝒆 = 𝑬
𝒍

IMPRESIÓN DE
RESULTADOS
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.

 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.

 El modelo desarrollado para obtener los objetivos propuestos en este trabajo es

apropiado, ya que se presentó adecuadamente los elementos y la mayor cantidad
de detalles estructurales presentes en la pieza real.
 El elemento finito puede servir para asegurar el nivel de confianza en las pruebas
de laboratorio. El método de elementos finitos es un método bastante cercano a la
realidad ya que los errores bastante pequeños.