Escuela Superior Politécnica del Litoral

Facultad de Ingeniería en Mecánica y Ciencias de la Producción

Asignatura

Mecánica de Sólidos II–MECG1015

Informe de la Práctica 1

“ANÁLISIS ESTRUCTURAL DE UNA ARMADURA POR

EL MÉTODO DE ELEMENTOS FINITOS (MEF)”

Estudiantes
Jorge Guido Espinosa Morocho - [email protected]
Samir Echeverria – [email protected]
Edison Contreras Herrera – [email protected]
Docente
Ing. Jorge Marcial Hernández

Paralelo 2

Fecha de entrega
4 de diciembre de 2019
1. DESCRIPCIÓN DEL PROBLEMA
Esta práctica se realizó con la finalidad de diseñar una armadura para que los esfuerzos en sus
MN MN
elementos no excedan 150 2 en tensión ni los
120 2 en compresión, con el fin de evitar
m m
fallos por deformaciones plásticas ni pandeo respectivamente, para ello asumimos un módulo de
Young de 200Gpa y v=0.3.con esta práctica se logró determinar el área transversal de todos los
elementos de la armadura de manera analítica y por el método de elementos finitos, finalmente
se realizó un contraste entre los resultados obtenidos entre ambos métodos, la estructura que
será analizada es la siguiente:

Ilustración 1 Estructura Analizada

2. SOLUCIÓN ANALÍTICA DEL PROBLEMA

Diagrama de cuerpo libre

Ilustración 2 Diagrama de cuerpo libre

pág. 1
Reacciones
Primero definiremos el sistema de referencia:

Ilustración 3 sistema de referencia positivo

Además los momentos anti horarios son considerados positivos.

Para calcular la reacción en dirección x en el nodo F planteamos sumatoria de fuerzas en X:

∑ X=0
3 ( 135 )−Fx=0
Fx=405 KN
Para encontrar la reacción en el nodo G aplicamos una sumatoria de momentos en el nodo F:

∑ M F=0
4.5 Fg=1944
Fg=432 KN

Finalmente para encontrar la reacción en dirección y en el nodo F planteamos sumatoria de

fuerzas en Y:

∑ Y =0
Fy+ Fg=0
Fy=432 KN
Fuerzas axiales

Nodo F:

F FD =432 KN (t)

F FG=405 KN (t)

Ilustración 4 Diagrama de cuerpo libre nodo F

pág. 2
Nodo G:

∑ X=0
F FG + F DG cosθ=0

F DG =−F FG /cosθ

F DG =459 KN (c)

Ilustración 5 Nodo G

∑ Y =0
F EG + F DG senθ+ F G=0

2.4
F EG =−432+ 459( )
5.1

F EG =216 ( c )

Nodo D:

∑ X=0
F Y +135+ F X =0

F DE =405−135

F DE =270 KN (t )

Ilustración 6 Nodo D

pág. 3
 ∑ Y =0
F BD−F Y + F G −F¿ =0

F BD=F Y −F G + F¿

F BD=216 KN (t )

Nodo E:

∑ X=0
−F DE −F EB cosθ=0

F EB=F DE / cosθ

F EB=306 KN ( c)

Ilustración 7 Nodo E

∑ Y =0
F EC −F EB senθ+ F EG =0

2.4
F EC =( 306 ) −216
5.1
F EC =72 KN (c )

Nodo C:

∑ Y =0
F AC =F EC /senθ

F AC =72/(2.4 /5.1)

F AC =153 KN (c)

pág. 4
 Ilustración 8 Nodo C

∑ X=0
F BC =F AC cosθ

4.5
F BC =153( )
5.1
F BC =135 KN (t)

Nodo A ∑ Y =0
F AC senθ−F AB=0

2.4
F AB=(153)( )
5.1

F AB=72 KN (t )
Ilustración 9 Nodo A

Áreas de las secciones transversales

Se van a dimensionar las áreas transversales de cada barra acorde a los requerimientos de
diseño:

 Esfuerzo Máximo de tensión (σ t): 150 MPa

 Esfuerzo Máximo de compresión (σ c): 120 MPa
 Módulo de Young (E): 200 GPa

F
Para ello emplearemos la siguiente ecuación: A=
σ
Ecuación 1

Barra AB: A=9 X 10−4 m 2

72 X 103 N Barra CE:
A=
150 X 106 Pa
72 X 103 N
A=
A=4.8 X 10−4 m2 120 X 106 Pa
Barra AC: A=6 X 10−4 m 2
72 X 103 N Barra BE:
A=
150 X 106 Pa
306 X 103 N
A=
A=4.8 X 10−4 m2 120 X 106 Pa
Barra BC: A=2.55 X 10−3 m2
135 X 103 N Barra BD:
A=
150 X 106 Pa

pág. 5
 216 X 103 N Barra DF:
A=
150 X 106 Pa 432 X 103 N
A=
A=1.44 X 10−3 m 2 150 X 106 Pa

A=2.88 X 10−3 m2
Barra DE: Barra FG:

270 X 103 N 405 X 103 N

A= A=
150 X 106 Pa 150 X 106 Pa

A=1.8 X 10−3 m2 A=2.7 X 10−3 m2

Barra DG: Barra EG:

459 X 103 N 216 X 103 N

A= A=
120 X 106 Pa 120 X 106 Pa

A=3.825 X 10−4 m 2 A=1.8 X 10−4 m 2

Deformaciones axiales
Barra AB:

150 x 106 Pa∗2.4 m Barra CE:

δ AB=
200 x 109 Pa 120 x 10 6 Pa∗2.4 m
δ CE =
δ AB=1.8 x 10−3 m 200 x 109 Pa

δ CE =1.44 x 10−3 m

Barra AC:

120 x 106 Pa∗5.1 m Barra BE:

δ Ac =
200 x 10 9 Pa 120 x 106 Pa∗5.1 m
δ BE =
δ Ac =3.06 x 10−3 m 200 x 10 9 Pa

δ BE =3.06 x 10−3 m
Barra BD:
Barra BC: 150 x 10 6 Pa∗2.4 m
δ BD =
150 x 106 Pa∗4.5 m 200 x 109 Pa
δ BC =
200 x 10 9 Pa δ BD =1.8 x 10−3 m
δ BC =3.375 x 10−3 m
Barra DE:

pág. 6
 150 x 10 6 Pa∗4.5 m δ DF =1.8 x 10−3 m
δ DE =
200 x 109 Pa

δ DE =3.375 x 10−3 m Barra FG:

150 x 106 Pa∗4,5 m

δ Ac =
Barra DG: 200 x 109 Pa

120 x 10 6 Pa∗5.1m δ Ac =3.375 x 10−3 m

δ DG =
200 x 10 9 Pa

δ DG =3.06 x 10−3 m Barra EG:

120 x 10 6 Pa∗2.4 m
δ EG =
Barra DF: 200 x 109 Pa

150 x 106 Pa∗2.4 m

δ DF =
200 x 109 Pa

δ EG =1.44 x 10−3 m Esfuerzos axiales y esfuerzos axiales máximos

Ya que el diseño requería que el esfuerzo en cada elemento sea menor o igual a 150MPa en
tensión y 120 MPa en compresión para evitar el pandeo, los esfuerzos a los que están sometidos
los elementos de la estructura son 150MPa para los elementos que se encuentran en tensión y
120 MPa para los que se encuentran comprimidos, siendo estos valores los esfuerzos máximos
te tensión y compresión respectivamente.

Fuerza Magnitud Tensión

en de (t) Área Deformación Esfuerzos
elemento Fuerza -
[KN] Compresión [m¿¿ 2]¿ −3
x 10 [m] [MPa]
(c)
F AB 72 t 4.800 X 10−4 1.800 150
F AC 153 c 1.275 X 10−3 3.060 120
F BC 135 t 9.000 X 10−4 3.375 150
F CE 72 c 6.000 X 10− 4 1.440 120
F BE 306 c 2.550 X 10−3 3.060 120
F BD 216 t 1.440 X 10−3 1.800 150
F DE 270 t 1.800 X 10−3 3.375 150
F DG 459 c 3.825 X 10−4 3.060 120
F DF 432 t 2.880 X 10−3 1.800 150
F FG 405 t 2.700 X 10−3 3.375 150
Tabla 1 Datos
de
F EG 216 c 1.800 X 10−41.440 120 Fuerza,
Área y
Deformación de cada elemento de la estructura

Diagrama del elemento diferencial cuadrado con los esfuerzos aplicados correspondientes
Para plantear los esfuerzos en un elemento diferencial usaremos los esfuerzos máximos de
tensión y compresión:

pág. 7
  Esfuerzo Máximo de tensión (σ t): 150 MPa
 Esfuerzo Máximo de compresión (σ c): 120 MPa

Barra AB:

σy σy

Ilustración 10 estado de esfuerzos barra AB en tensión

Barra CE:

σy σy

Ilustración 11 Estado de esfuerzos barra CE

Círculo de Mohr: Esfuerzos principales y esfuerzo cortante máximo

Barra AB:

τ max=σ y /2

τ max=75 MPa

σ max=150 MPa

σ min=0 Mpa

Ilustración 12 Circulo de Mohr Barra AB

Barra CE:

τ max=σ y /2

τ max=60 MPa

σ max=−120 MPa

pág. 8
 σ min=0 Mpa

Ilustración 13 Circulo de Mohr Barra CE

3. SOLUCIÓN DEL PROBLEMA POR ANÁLISIS DE ELEMENTOS FINITOS

Descripción del modelo de elementos finitos

El modelo de elementos finitos es un método que se usa para resolver problemas ingenieriles de
manera rápida, sin la necesidad de establecer modelos matemáticos complejos de solución muy
difícil y posiblemente sujeta a errores de cálculo. El manejo de esta herramienta es relativamente
sencillo, sin embargo, para plasmar el resultado del modelo de elementos finitos en la realidad, el
análisis de dicho modelo debe ser planteado por personas que tengan conocimientos en el área de
ingeniería, de esta manera se puede garantizar que el resultado del modelo de elementos finitos
arrojara resultados con aproximaciones coherentes y bastante aproximados a los esfuerzos reales a
los que estará sometido el sólido que estamos analizando. Este método se basa en dividir un cuerpo
complejo continuo en un conjunto de elementos conectados entre sí, de esta manera, el análisis no se
realiza sobre todo el cuerpo complejo si no, el software se encarga de analizar uno a uno dichos
elementos, arrojando los parámetros de interés por el cual se realiza el análisis.

 Descripción del software utilizado (ANSYS APDL)

ANSYS Mechanical APDL (Advance Parametric Design Language) es una herramienta
computacional de elementos finitos utilizado para el análisis de estructuras estáticas
lineales, térmicos y demás parámetros necesarios para el diseño de estructuras.
 Descripción del tipo de elementos usado (grados de libertad)
La estructura es del tipo isostática, lo que significa que su grado de libertad es igual a cero.
También indica que la estructura se encuentra en equilibrio y las ecuaciones de la estática
son suficientes para determinar las fuerzas internas y reacciones.

 Materiales
Acero (lineal, elástico, isotrópico).

 Secciones (definir una sección diferente para cada miembro de la armadura de

acuerdo a las áreas calculadas en la solución analítica)
11 secciones

 Tamaño del elemento

180 Stn

 Número de elementos (un elemento por cada miembro de la armadura)

pág. 9
 11 elementos

 Número de nodos
7 nodos

 Restricciones de los grados de libertad en los apoyos

Apoyo de primer orden: Un rodillo que restringe un grado de libertad (movimiento
perpendicular) y genera una reacción perpendicular a la base o al movimiento permitido.
Apoyo de segundo orden: Una articulación que restringe dos (2) grados de libertad,
permitiendo únicamente el giro o rotación alrededor de sí misma, e impidiendo el
movimiento vertical y horizontal.

 Cargas
Tres (3) cargas horizontales de 135 [KN] cada uno en los nodos (D, B, A) o (2, 3, 4).

 Gráfico del modelo

pág. 10
Resultados obtenidos con el modelo de elementos finitos (gráficos y valores)

 Reacciones en los apoyos

pág. 11
 Fuerzas axiales (todos los elementos)

 Áreas de la secciones transversales

pág. 12
 Deformaciones axiales

 Esfuerzos axiales

pág. 13
  Esfuerzos axiales máximos de tensión y compresión

 Deflexión máxima de la armadura (valor absoluto)

Deflexión máxima de 0.029098

4. VALIDACIÓN DEL MODELO MEF

Comparación de los resultados analíticos vs. MEF (tabla de porcentajes de error)

pág. 14
 Reacciones en los apoyos

Nodo Reacciones en los Reacciones en los Porcentaje de

(MEF - dibujo original) apoyos apoyos Error%
MEF Método analítico
Nodo 1 – F Fx = - 405 [KN] Fx = - 405 [KN] 0%
Fy = - 432 [KN] Fy = - 432 [KN] 0%
Nodo 7 – G Fy = 432 [KN] Fy = 432 [KN] 0%

 Fuerzas Axiales

Barra Fuerzas Axiales Sección Fuerzas Axiales Porcentaje de

(MEF) MEF [KN] (dibujo Método analítico Error%
original) [KN]
L1 432 FD 432 0%
L2 216 DB 216 0%
L3 72 BA 72 0%
L4 -153 AC -153 0%
L5 135 BC 135 0%
L6 -72 CE -72 0%
L7 -216 EG -216 0%
L8 270 DE 270 0%
L9 -306 BE -306 0%
L10 405 FG 504 0%
L11 -459 DG -459 0%

 Áreas de las secciones transversales

Barra Áreas de sección Sección Fuerzas Axiales Porcentaje de

(MEF) transversal (dibujo Método analítico Error%
MEF [m2] original) [m2]
L1 0.288 X 10−2 FD 2.88 X 10−3 0%
L2 0.144 X 10−2 DB 1.44 X 10−3 0%
L3 0.480 X 10−3 BA 4.80 X 10−4 0%
L4 0.1275 X 10−2 AC 1.275 X 10−3 0%
L5 0.900 X 10−3 BC 9.00 X 10−4 0%
L6 0.600 X 10−3 CE 6.00 X 10− 4 0%
L7 0.180 X 10−2 EG 1.80 X 10−3 0%
L8 0.180 X 10−2 DE 1.80 X 10−3 0%
L9 0.255 X 10−2 BE 2.55 X 10−3 0%
L10 0.270 X 10−2 FG 2.70 X 10−3 0%
L11 0.3825 X 10−2 DG 3.825 X 10−3 0%

pág. 15
  Deformaciones axiales

Barra Deformaciones Sección Deformaciones Porcentaje de

(MEF – Dibujo axiales (dibujo axiales Error
original) MEF [m] original) Método
analítico [m]
L1 – FD 0.750 X 10−3 FD 1.800 X 10−3 58.3%
L2 – DB 0.750 X 10−3 DB 1.800 X 10−3 58.3%
L3 – BA 0.750 X 10−3 BA 1.800 X 10−3 58.3%
L4 – AC −0.600 X 10−3 AC −3.060 X 10−3 80.4%
L5 – BC 0.750 X 10−3 BC 3.375 X 10−3 77.7%
L6 – CE −0.600 X 10−3 CE −1.440 X 10−3 58.3%
L7 – EG −0.600 X 10−3 EG −1.440 X 10−3 58.3%
L8 – DE 0.75 X 10−3 DE 3.375 X 10−3 77.7%
L9 – BE −0.600 X 10−3 BE −3.060 X 10−3 80.4%
L10 – FG 0.75 X 10−3 FG 3.375 X 10−3 77.7%
L11 – DG −0.600 X 10−3 DG −3.060 X 10−3 80.4%

 Esfuerzos axiales

Barra Esfuerzos axiales Sección Esfuerzos Porcentaje de

(ANSYS – ANSYS [MPa] (dibujo axiales Error
Dibujo original) original) Método
analítico [MPa]
L1 – FD 150 FD 150 0%
L2 – DB 150 DB 150 0%
L3 – BA 150 BA 150 0%
L4 – AC -120 AC -120 0%
L5 – BC 150 BC 150 0%
L6 – CE -120 CE -120 0%
L7 – EG -120 EG -120 0%
L8 – DE 150 DE 150 0%
L9 – BE -120 BE -120 0%
L10 – FG 150 FG 150 0%
L11 – DG -120 DG -120 0%

5. ANÁLISIS DE RESULTADOS
La deformación de la armadura presenta inclinación hacia el lado derecho, esto se debe
a las tres fuerzas puntuales que actúan sobre los nodos (D, B, A).
Se visualiza una mayor deformación en el Nodo A. Esto se debe a que existe mayor
distancia entre el nodo y el apoyo, lo cual conlleva a que el momento generado por la
fuerza de 135 [KN] sea mayor.

pág. 16
 Las deformaciones axiales presentan valores muy disparejos a los obtenidos
analíticamente, por lo que se cree que el programa introdujo un análisis más exhaustivo
con respecto a las propiedades del material y eso causó el elevado porcentaje de error
entre ambos valores.
6. CONCLUSIONES
La armadura cumple con los requerimientos de que sus esfuerzos no excedan los 150
[MPa] en tensión, ni los 120 [MPa] en compresión.
La vida deberá ser reforzada en las barras donde se aprecia un mayor intento a
flexionarse (color rojizo).
Se debería incrementar el área de las barras que están próximas a superar el límite de
esfuerzo permitido en tensión, para de esta manera disminuirlo.

pág. 17