“Año del Fortalecimiento de la Soberanía Nacional”

UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA MECÁNICA

ELEMENTOS FINITOS (MC516-E)
PRIMERA PRÁCTICA CALIFICADA
(TRACCIÓN SIMPLE UNIDIMENSIONAL)

Docente/Asesor a cargo : Ing. Miguel Ampuero Suárez

Integrantes Código

Rojas Povis, Eduardo Antonio 20204054F

Culli Echegaray Keneth Lyn 20190147B

Espinoza Rojas Luis Erasmo 20212071C

2023-I

14
 ÍNDICE

PRIMERA PRÁCTICA CALIFICADA........................................................................... 2

1. ENUNCIADO............................................................................................................... 2
2. RESOLUCIÓN DEL CÁLCULO ANALÍTICO - MODELADO DEL CUERPO REAL..3
3. RESOLUCIÓN EN MATLAB.......................................................................................7
4. SIMULACIÓN EN ANSYS.........................................................................................13
5. CONCLUSIONES....................................................................................................... 21
6. RECOMENDACIONES............................................................................................. 21
7. REFERENCIAS BIBLIOGRÁFICAS........................................................................ 22

1
 PRIMERA PRÁCTICA CALIFICADA

1. ENUNCIADO

Actividades a Desarrollar:
1.- Elaborar el cálculo analítico para 2 y 3 elementos finitos, realizar el cálculo matricial
y presentar como solución el campo de desplazamientos y esfuerzos, (vector Global
desplazamiento y esfuerzos), así como la reacción en el apoyo.
2.- Realizar la modificación en los códigos de MATLAB® del problema unidimensional
de la placa con sección variable estudiada en aulas.
3.- Realizar la simulación en ANSYS presentando los resultados el campo de
desplazamiento y esfuerzos (usar Criterio de Von Mises para la simulación).

2
2. RESOLUCIÓN DEL CÁLCULO ANALÍTICO - MODELADO DEL
CUERPO REAL
Se consideran tres elementos finitos. Para facilitar los cálculos se seleccionó las
siguientes alturas para cada elemento.

𝑙1 = 500𝑚𝑚
𝑙2 = 250𝑚𝑚
𝑙3 = 250𝑚𝑚

Y el ancho de cada elemento se calculó teniendo en cuenta que el volumen de la placa

triangular tiene que ser igual a la suma del volumen de los tres elementos finitos.

Para ello, calculamos el área promedio de cada base e igualamos al área

correspondiente a cada elemento.

𝑏1 = (1200 + 600)/2 −> 𝑏1 = 900𝑚𝑚

𝑏2 = (600 + 300)/2 −> 𝑏2 = 450𝑚𝑚
𝑏3 = (300 + 0)/2 −> 𝑏3 = 150𝑚𝑚

3
Con los datos de las bases de los elementos y las alturas, el modelo del triángulo
con 3 elementos finitos está dado por.-

Se puede calcular el área de las secciones transversales de los tres elementos finitos
delimitados anteriormente, haciendo uso de las siguientes fórmulas:

4
Cuadro de conectividad:

E NODOS GDL 𝑙𝑒 (𝑚𝑚) 𝐴𝑒 (𝑚𝑚 )

2

(1) (2) 1 2

1 1 2 1 2 900 135000

2 2 3 2 3 450 67500

3 3 4 3 4 150 22500

GRADOS DE LIBERTAD NODALES (VECTOR DESPLAZAMIENTO)

A través del gráfico se muestran los grados de libertad nodales globales:

Donde Q1 = 0 mm, pues la placa está empotrada y los demás desplazamientos son
incógnitas que tendrán que ser calculadas.

5
VECTOR CARGA GLOBAL

Analizando las fuerzas en cada elemento finito se dan bajo las siguientes fórmulas:

Donde el vector de carga se traduce de la siguiente manera:

La matriz de Rigidez se halla de la siguiente manera:

6
Finalmente se resuelve la Ecuación de Rigidez 𝐹𝑖 = 𝐾𝑖 . 𝑄𝑖

Con los datos obtenidos a partir de los cálculos y el desarrollo del modelado es
suficiente para hacer los respectivos cálculos y desplazamientos matriciales; la siguiente
parte del análisis se realizará mediante el software de Programación MatLab que facilita
la resolución de la ecuación de Carga y el cálculo de la fuerza; se toma esta medida
principalmente para no alargar de más el procedimiento.

3. RESOLUCIÓN EN MATLAB
Para este primer laboratorio se detallan los pasos que se utilizaron para la realización de
la implementación del código en MATLAB®.

Al inicio del código definimos las variables que se van a utilizar a la hora de realizar las
operaciones. En los comentarios se indican las unidades en las que se encuentran las
variables.

Luego, el código solicita al usuario introducir el número de elementos necesarios para

realizar el cálculo requerido. Con ese número, el código crea vectores columna que
contendrán las longitudes y la medida de las bases de todos los elementos.

7
El código hace uso de estructuras condicionales y bucles para introducir los datos
correspondientes de los elementos en los dos vectores previamente creados.

Para introducir las longitudes, el programa localiza cuáles son los dos elementos que
son vecinos con el nodo en el que actúa la fuerza externa; y les asigna longitudes que
permitan que exista un nodo en la marca de 500 mm de longitud, y que esta marca no
esté en medio de un elemento.

Después, definimos los vectores columna que contienen los datos del área transversal de
cada uno de los elementos , y el peso de estos elementos. También se define un vector
columna que contiene los valores de las fuerzas que actúan en cada nodo; por lo cual el
número de elementos de este vector es una unidad mayor al número de elementos.

El código introduce los valores de las fuerzas que actúan en cada nodo haciendo uso de
un bucle, y sumando para cada nodo la mitad del valor de los pesos de los dos
elementos vecinos. Después, el código hace uso de una estructura condicional para
establecer dónde se debe ubicar el nodo sobre el cual actúa la fuerza externa.

8
Luego de eso, definimos un vector columna que contenga los datos de las constantes
elásticas de los elementos; y definimos también el vector asociado al desplazamiento de
los nodos, y la matriz de rigidez asociada al sistema de elementos finales. Hay que tener
en cuenta que el orden del vector y la matriz relacionados con los nodos es una unidad
mayor que el número de elementos finitos.

El programa calcula primero los valores de los cuatro elementos de la sub-matriz

cuadrada ubicada lo más cerca posible a la esquina superior de la matriz. Después de
eso, se ejecuta un bucle, de tal manera que en cada iteración de bucle se produce una
nueva suma de matrices con sub-matrices; mediante la cual se van definiendo los
valores de los elementos de la matriz rigidez del cuerpo evaluado en este informe.

Finalmente, el programa calcula el valor de sub-vector columna conformado por los

elementos del vector columna de desplazamientos cuyos valores no son conocidos; para
lo cual se calcula el producto entre la inversa de la submatriz de la matriz de rigidez de
orden “n” ubicada en la esquina inferior derecha de la matriz con el sub-vector
conformado por los elementos de valores conocidos del vector columna “F”.

Para calcular todos los elementos del vector columna “F”, se multiplican la matriz de
rigidez de los elementos junto con el vector columna de los desplazamientos de los
nodos; de tal manera que ahora se conoce el valor del primer elemento desconocido del
vector columna “F”, el cual es el único elemento de ese vector que es negativo.

9
Después, el programa calcula el valor del vector columna de los esfuerzos que afectan a
los elementos del cuerpo, mediante el uso de un bucle que realiza productos de vectores.
Además, el programa calcula el valor de la reacción que ejerce el apoyo del cuerpo, para
lo cual toma el valor del primer nodo del cuerpo; y resta de ese valor el valor del peso
correspondiente al primer elemento del cuerpo.

Finalmente, el programa verifica la precisión de los cálculos restando el valor de la

reacción del apoyo obtenida por los elementos finitos con el valor de la reacción que se
obtiene aplicando la primera condición de equilibrio; ya que en el cuerpo evaluado en
este informe actúan tres fuerzas externas, las cuales serían el peso del cuerpo, la fuerza
“P” y la reacción en el apoyo.

El programa expresa esta diferencia como el error absoluto del método de los elementos
finitos, y luego divide esta diferencia entre el valor de la reacción en el apoyo obtenida
por la primera condición del equilibrio; con lo cual se calcula el error relativo del
método de los elementos finitos.

Una vez terminados los cálculos, el programa permite al usuario visualizar las tres
incógnitas solicitadas por el ejercicio en la ventana de comandos de MATLAB. Para
ello, se hace uso de los comandos “display”; los cuales muestran los vectores columna
del desplazamiento de los nodos y de los valores de los esfuerzos de los elementos.

Mediante estos comandos, el programa también muestra al usuario el valor de la

reacción de apoyo calculada por el método de los elementos finitos; así como también
los errores absoluto y relativo asociados al método de los elementos finitos.

10
En la imagen se muestra el resultado de los cálculos realizados por el programa al
evaluar dos elementos, junto con las unidades del Sistema Internacional en las que se
expresan las unidades de medida correspondientes a esos valores. Es posible apreciar
que los errores obtenidos en el proceso de cálculo son relativamente pequeños, incluso
para una cantidad reducida de elementos finitos.

11
Los siguientes resultados se ofrecen al descomponer el triángulo en tres elementos
distintos. Se puede notar que los valores de los dos errores han incrementado
considerablemente sus valores con respecto a cuando se usaron dos elementos finitos,
debido que el número de elementos finitos empleado en esta ocasión era reducido e
impar.

También se puede ejecutar el programa para un número de elementos igual a 1000, por
poner un ejemplo. En ese caso, se puede apreciar que el valor de la reacción en el apoyo

12
converge a un valor numérico similar al que se obtuvo cuando se trabajó con dos
elementos finitos; pero con errores de mayor valor numérico.

Sin embargo, algo que se debe tener en cuenta es que la magnitud de los errores con un
número de elementos par es similar a la magnitud de los errores que se generan cuando
se usa un número de elementos impar; siempre y cuando la cantidad de elementos
empleados sea grande. Esto se puede evidenciar para el siguiente ejemplo, en el que se
usan 1001 elementos finitos:

4. SIMULACIÓN EN ANSYS
Para este primer laboratorio se detallan los pasos que se utilizaron para la realización de
la simulación en el software de Ansys.

Creamos un archivo nuevo y seleccionamos ‘Static Structural’. En el apartado de

Engineering Data aplicamos la siguiente conversión.
Conversión de

5 2 5
Modulo de Young : E = 3. 0𝑥10 𝑁/𝑚𝑚 <> 3. 0𝑥10 𝑀𝑃𝑎

13
 −6 3 −6 9 3
Peso específico : 78, 45𝑥10 𝑁/𝑚𝑚 <> 78, 45𝑥10 𝑥10 / 9, 81 𝑘𝑔/𝑚

En el apartado de ‘Geometría’ seleccionamos Line Bodies y 2D.

14
Empezamos el sketch. Asimismo ajustamos la unidad

Una vez creado el skectch darle el espesor de 150 mm vamos a ‘Model’.

15
En esta aparte podemos variar parámetros del Modelo

16
Para el desarrollo de la Primera Práctica Calificada

Al no dejar insertar una fuerza en el medio de la placa triangular, podemos obtener el

mismo resultado por medio de una fuerza nodal:

17
Indicamos el nodo por medio de un Mesh- Cuadrangular-Triangular

Insertamos el resultado para la deformación total:

18
Insertamos el resultado para la esfuerzo (Von Mises):

Tenemos los siguientes resultados; corroboramos

una respuesta similar si usamos un Mesh
Triangular.

19
Resultados Finales

20
5. CONCLUSIONES

● Se puede notar que las deformaciones son extremadamente pequeñas, del orden
de milésimas, lo que demuestra que la placa triangular, al ser sometida a la carga
P, se mantiene como un cuerpo rígido.

● Los cálculos de esfuerzos son positivos debido a que el punto 2 de cada

componente se desplaza más que el punto 1, indicando que el componente está
sometido a tracción en la dirección de nuestra referencia "x".

● Por otro lado, si los desplazamientos y esfuerzos calculados son negativos, se

debe a que la carga P y el peso de la placa son antiparalelos con la dirección de
nuestra línea de acción x, lo que significa que la placa está en compresión.

● Para sistemas sujetos a tracción simple, el procedimiento será idéntico al

previamente empleado.

● En este contexto, se puede observar que MATLAB ahorra significativamente

tiempo en los cálculos, ya que el mismo código puede aplicarse a otro sistema de
tres elementos, lo que lo convierte en una herramienta altamente útil.

● Los resultados obtenidos en MATLAB presentan un error prácticamente nulo, lo

que sugiere que la aproximación mediante tres elementos finitos es casi
perfectamente precisa.

● Un detalle interesante para tener en cuenta, es que se puede apreciar que los
errores obtenidos a la hora de ejecutar el código de MATLAB son mayores a
medida que mayor es el número de elementos finitos empleados en el método,
siempre y cuando se use un número de elementos par.

● Asimismo, para este método de resolución en particular, se puede apreciar que al

elegir una cantidad de elementos impar, se producen mayores errores que
cuando se usa una cantidad de elementos par; y también que la magnitud de los
errores se reduce mientras mayor es la cantidad de elementos empleados.

● También es evidente que el eje de referencia "x" inicialmente elegido debe ser
establecido con experiencia, visualizando cómo el modelo se deformaría en
primer lugar.

6. RECOMENDACIONES
Es importante tener en cuenta algunas recomendaciones clave para garantizar resultados
precisos y eficientes:

21
 ● Es fundamental comprender completamente el problema que estás tratando de
resolver y el comportamiento del sistema. Esto te ayudará a definir
correctamente el modelo de elementos finitos.

● La calidad de la malla (discretización) es esencial. Asegúrate de que la malla

esté bien refinada en áreas críticas y que no haya elementos de malla con
relaciones de aspecto extremas.

● Para comprender cómo ciertos parámetros afectan los resultados, realiza análisis
de sensibilidad variando las propiedades del material o las condiciones de carga.

● Asegúrate de mantener tus herramientas de software actualizadas con las últimas

correcciones de errores y mejoras.

● No se recomienda usar este método en MATLAB para calcular el valor de la

reacción en el apoyo, debido a que el valor de los errores que se producen a
partir de este método se incrementa a medida que se usan más elementos.

● Este método debería usarse solamente para determinar los desplazamientos de

los nodos del cuerpo a evaluar, junto con el valor de los esfuerzos en cada
elemento; con el fin de determinar el comportamiento mecánico del elemento
frente a una fuerza externa.

7. REFERENCIAS BIBLIOGRÁFICAS

● “Introducción a los elementos finitos en ingeniería”. Ashok D. Belegundu,

Tiruphati R. Chandrupatla.

● “Análisis de elementos finitos”. Saeed Moaveni.

● https://campusvirtualfim.com/pluginfile.php/200/course/section/874/CURSO
%20MC516%20CCLAB%202A%20%28P%29.ppsx
Ampuero Suárez, Miguel Ángel

22