DESCRIPCION GENERAL

La solución mediante métodos numéricos ha permitido desarrollar diferentes formas de

solucionar un problema de manera optima y eficiente, de esta forma ha permitido dar
desarrollo a problemas complejos con soluciones aproximadas. “Un método numérico es un
procedimiento el cual se obtiene, casi siempre de una manera aproximada, la solución de
ciertos problemas realizando cálculo puramente aritméticos y lógicos” (Seminario,
Ricardo,2017)
La solución de problemas mediante el uso de métodos numéricos ha permitido formar
herramientas y técnicas para la solución de problemas, como lo son herramientas iterativas,
estadísticas y la técnica del método de elementos finitos (FEM).
En este trabajo nos enfocaremos en la técnica del método de elementos finitos, dando a
conocer en que consiste, que tipo de ecuaciones son utilizadas en las herramientas
aplicativas, y como se aplica en los análisis de modelos bidimensionales y tridimensionales,
al igual de como se realiza el estudio de análisis modal y de fatiga, basándonos en estudios
realizados, libros, artículos y programas que permiten realizar el estudio de esta técnica.

ANALISIS DE MODELOS BIDIMENSIONALES Y TRIDIMENSIONALES

La amplia aplicación del método de elementos finitos permite que analicemos estructuras,
sólidos en revolución y flexión en placas y laminas, entre muchas otras aplicaciones. En la
ingeniería se encontrarán problemas tanto bidimensionales como tridimensionales, en los
cuales hay que tener en cuenta los parámetros de estudio para determinar la forma en la que
se va a estudiar el problema.
MODELOS BIDIMENSIONALES
En los modelos bidimensionales se analiza casos en un medio continuo plano, este se
considera que está situado en el plano XY. El espesor t, en dirección transversal, este es
considerado despreciable, en cuanto al plano XY.
La posición de un punto esta definida por dos coordenadas ( x , y ) , y su deformación tiene
dos componentes u ( x , y ) , v ( x , y ), en las direcciones x, y respectivamente. El campo de
deformaciones es por lo tanto un vector:

Figura . Vector de campo de deformaciones

Dentro de la elasticidad en dos dimensiones existes dos problemas diferentes:
 TENSIÓN PLANA: Cuando la tensión σ z en sentido perpendicular al plano XY es
cero, ya que el solido puede dilatarse libremente en el sentido de su espesor. Por lo
tanto, existe una deformación unitaria ϵ z no nula en dicha dirección.
(Celigüeta,2011)

Para que una pieza se encuentre bajo un estado de tensión plana, se tienen que
cumplir las siguientes condiciones:
1. La geometría de la pieza debe tener un plano de simetría XY.
2. El espesor debe ser mucho menor que las dimensiones de la pieza en el
plano XY. No es necesario que el espesor sea uniforme en toda la pieza.
3. Las cargas y condiciones del contorno deben estar aplicadas en el plano de
simetría XY. (Abasolo, Coria, Heras,2019)

 DEFORMACION PLANA: Cuando en el sentido del espesor del sólido no hay

posibilidad de deformación, es decir,
ϵ z =0 por lo que se genera una TENSIÓN en dicha dirección σz no nula.
(Celigüeta,2011)

Para que una pieza se encuentre en un estado de deformación plana, se tienen que
cumplir las siguientes condiciones:
1. La pieza debe tener una geometría con simetría respecto de un eje.
2. Las cargas y condiciones de contorno también deben tener una simetría
respecto al eje. (Abasolo, Coria, Heras,2019)

DEFORMACIONES UNITARIAS
Las deformaciones unitarias en un punto del elemento finito son:

Figura . Vector deformación unitaria

Donde ε es la deformación unitaria, δ son las deformaciones, u es componente del vector
de campo de deformaciones, v es componente del vector de campo de deformaciones.
Al sustituir las deformaciones u en función de las deformaciones nodales a través de
funciones de interpolación, se obtiene:

Figura . Ecuación deformación unitaria.

Donde δ e son las deformaciones nodales y B representa el campo de deformaciones.
Las deformaciones unitarias en el interior del elemento se pueden poner en función de las
deformaciones nodales en la forma:

Figura . Ecuaciones de las componentes de las deformaciones unitarias nodales.

Donde ε x es la deformación unitaria en el eje de la x, ε y es la deformación unitaria en el eje
de la y, U y V son las funciones de interpolación. (Celigüeta,2011)

RELACION TENSIÓN-DEFORMACIÓN
La relacion entre TENSIÓNes y deformaciones se deduce de la ecuación constitutiva de la
elasticidad tridimensional.

Figura . Relacion tensión- deformación

Donde σ es la tensión, D es la matriz de constantes elásticas y ε es la deformación unitaria.
Del teorema de Maxwell-Betti1 se deduce que D es siempre simétrica y d 12=d21. Para la
elasticidad isótropa se tiene:
 TENSIÓN PLANA

Figura . Elasticidad Isótropa para la tensión plana

Donde E es el módulo de elasticidad del material y v el coeficiente de Poisson.
DEFORMACIÓN PLANA

Figura . Elasticidad Isótropa para la deformación plana.

Donde E es el módulo de elasticidad del material y v el coeficiente de Poisson. (Solidos
Bidimensionales)

1.Teorema Maxwell-Betti: Describe la relación entre los trabajos realizados por dos estados de
carga actuantes sobre una misma estructura para que la energía de deformación sea independiente
del orden de aplicación de estos estados. (Análisis estructural)

ELEMENTO TRIANGULAR
Este elemento tiene seis desplazamientos en los nudos, que forma un vector:

Figura . Vector deformaciones nodales

Donde δ ees el vector de deformaciones nodales y U, V son las componentes del vector.
Los desplazamientos de un punto cualquiera dentro del elemento se pueden representar en
función de estos seis valores nodales, mediante una expresión polinómica. Dado que hay
seis deformaciones nodales, el polinomio solo podrá tener seis términos:

Figura . Expresiones polinómicas

Donde u y vson las expresiones del vector de deformaciones nodales, α son las longitudes
de la pieza.
Al ser tres funciones de interpolación planos el campo de desplazamientos en el interior del
elemento es también un plano que pasa por lo tres valores nodales del desplazamiento. En
consecuencia, si se emplea este elemento, el campo de deformaciones en una dirección
cualquiera u o v se aproxima mediante una superficie poliédrica de facetas triangulares.

Figura . Funciones de interpolación del elemento triangular

Donde N1, N2, N3 son funciones de interpolación correspondiente a los tres nudos del
triángulo.
El estado de tensiones en un punto cualquiera del elemento es:

Figura . Ecuación de estado de tensiones en un punto cualquiera.

Donde σ es la tensión, D es la matriz de constantes elásticas y δ e son las deformaciones
nodales. (Celigüeta,2011)
ELEMENTO RECTANGULAR
El elemento rectangular al igual que el elemento triangular tiene un análisis de nodos, la
diferencia es que en el elemento rectangular se analizarán cuatro nodos.

Figura . Elemento placa rectangular.

Se supone que el elemento rectangular tiene lados 2a y 2b. (Formulación de elementos
finitos)
APLICACIONES
Las aplicaciones de los elementos finitos bidimensionales son muy amplias, a
continuación, presentaremos algunas aplicaciones por medio de artículos científicos:
-Optimization of finite element bidimensional models: an approach based on genetic
algorithms (Annicchiarico,W, Cerrolaza,M, 1998)
Se observa como la aplicación del método de elementos finitos no solo se usa para
determinar tensiones o esfuerzos en una pieza, también es empleado en la biología
genética, aplicando el método de elementos finitos para optimizar el cruce, la mutación ya
la reproducción en una especie.
- A stability study of some mixed finite elements for large deformation elasticity problems.
(Auricchio, F, Beirao de Veiga, L, Lovadina, C, Reali, A, 2004)
Se denota que el estudio de elementos finitos permite estudiar materiales incompresibles,
realizando un análisis de manera bidimensional, además estudiando elementos finitos
mixtos buscando soluciones discretas.
MODELOS TRIDIMENSIONALES
El cálculo de tensiones y deformaciones en un solido tridimensional es un problema que no
tiene mayor complejidad conceptual que el caso bidimensional, por lo que el MEF se aplico
desde un principio a este tipo de problemas.

Las ecuaciones diferenciales que controlan el problema son similares a las del
problema bidimensional, pero incluyendo una tercera coordenada z, y un tercer
desplazamiento en dicha dirección w. Estas ecuaciones son de orden m=2. La energía
potencial del sistema contiene a las derivadas primeras de las deformaciones, es decir
n=1, con lo que es suficiente con emplear funciones de interpolación de tipo C0 para
asegurar la convergencia en este problema. (Celigüeta,2011)

Un análisis tridimensional tiene varias desventajas comparado con un análisis

bidimensional:
 El modelo tridimensional tiene mas nodos y mas grados de libertad, con lo que el
coste asociado es mayor.
 La visualización de los resultados es mas compleja, ya que los resultados de los
puntos interiores solo se pueden visualizar realizando cortes en la pieza.
 Se complica el proceso de mallado especialmente con refinados locales porque
aparecen elementos de mala relacion de aspecto (muy irregulares). (Abasolo, Coria,
Heras,2019)
 Figura . Grados de libertad de piezas tridimensionales.

DEFORMACIONES UNITARIAS
En un punto del solido tridimensional, el vector de deformaciones unitarias contiene seis
términos, que corresponden a las tres deformaciones unitarias ε, y las tres deformaciones
de cortadura ingenieriles γ .

Figura . Deformaciones unitarias tridimensionales.

Al igual que en el modelo bidimensional, la ecuación de deformación unitaria se define
como:

Figura . Ecuación deformación unitaria.

Donde δ e son las deformaciones nodales y B representa el campo de deformaciones.
(Celigüeta,2011)

Existe un tipo de elemento bidimensional que se utiliza en análisis tridimensionales,

elemento cascara. Es un triángulo o cuadrilátero, pero tiene la capacidad de deformarse
fuera del plano. Generalmente su usan para mallar superficies de pequeño espesor bajo
cargas y condiciones de contorno que no tienen que ser aplicadas en el plano del modelo.

Figura . Elemento cascara.

La matriz de rigidez del elemento cascara se obtiene a partir de una formulación

simplificada de la teoría de la elasticidad aplicada a su geometría. (Abasolo, Coria,
Heras,2019)

Existen tipos de elementos en el modelo tridimensional como lo son:

 Elementos prisma rectangular
 Elemento lineal de ocho nudos
 Elemento cuadrático de veinte nudos
 Elemento cubico de treinta dos nudos
 Elementos langraviatos
 Elementos con numero de nudos por lado variable
 Elementos tetraedro
 Elementos tetraedro de cuatro nudos
 Elementos tetraedros cuadrático de diez nudos
 Elemento tetraedro de veinte nudos
 Elementos prisma triangular

APLICACIONES
-Initial stress induced in periodontal tissue with diverse degrees of bone loss by an
orthodontic force: Tridimensional analysis by means of the finite element method (Cobo,
Sicilia, Arguelles, Suarez, Vijande, 2007)

Observamos como el método de elementos finitos también contribuye a lo estudios

odontológicos y avances tecnológicos de estos por medio de esta técnica iterativa

-Rehabilitation of the atrophic mandible with short implants in different

positions: A finite elements study (Peixoto, Camati, Faot, Sotto, Martinez, Peruzzo, 2017)

Se observa como aplican el metodo tridimensional de elementos finitos para obtener lo

posibles esfuerzos que se generan en un implante bucal.
 Figura . Esfuerzos ejercidos en un implante bucal.

ANALISIS

BIBLIOGRAFIA

1. Seminario, Ricardo, Métodos numéricos para ingeniería, 2017, pag.4, Recuperado

de: https://disi.unal.edu.co/~lctorress/MetNum/LiMetNu2.pdf
2. Celigüeta, Juan, Método de los elementos para análisis estructural, 2011, pag.39,
Editorial TECNUN.
3. Abasolo, Mikel, Coria, Ibai, Heras, Iker, Análisis de métodos bidimensionales,
2019, pag.3-6, Recuperado de:
https://ocw.ehu.eus/pluginfile.php/52302/mod_resource/content/0/TEMA
%204%20AN%C3%81LISIS%20DE%20MODELOS
%20BIDIMENSIONALES.pdf
4. Análisis estructural, Teorema de Maxwell-Betti, Recuperado de:
file:///C:/Users/Lenovo/Downloads/Teorema_de_Maxwell_Betti.pdf
5. Solidos Bidimensionales,pag.4-6, Recuperado de:
http://cv.structuralia.com/moodle/cotenidos_NO_scorm/Version4/pdf/tema4/cap4.p
df
6. Formulacion de elementos finitos, pag.68. Recuperado de:
http://bibing.us.es/proyectos/abreproy/70238/fichero/Capitulo+5.pdf
7. Annicchiarico,W, Cerrolaza,M, Optimization of finite element bidimensional
models: an approach based on genetic algorithms,1998.
8. Cobo,Juan, Sicilia,Alberto, Arguelles,Juan, Suarez,David, Vijande,Manuel, Initial
stress induced in periodontal tissue with diverse degrees of bone loss by an
orthodontic force: Tridimensional analysis by means of the finite element method,
2007.
9. Peixoto, Hugo, Camati, Paulo, Faot, Fernanda Sotto, Bruno, Martinez, Elizabeth,
Peruzzo, Daiane, Rehabilitation of the atrophic mandible with short implants in
different positions: A finite elements study, 2017