UNIVERSIDAD ALAS PERUANAS

FACULTAD DE INGENIERÍA Y ARQUITECTURA

ESCUELA PROFESIONAL DE INGENIERÍA CIVIL

TRABAJO:

METODO DE ELEMENTOS FINITOS

DOCENTE
INGENIERO EDWIN RICARDO RODRIGUEZ PLASENCIA

ALUMNO:
CABANILLAS QUESQUEN JOSE CARLOS

CAJAMARCA – PERÚ

2019
INTRODUCCIÓN AL MÉTODO DE LOS ELEMENTOS
FINITOS

1.1 INTRODUCCIÓN

El método de elementos finitos es una esencia para un método número para las
resoluciones de ecuaciones diferenciales como tal vez ya sabremos las ecuaciones
diferenciales rigen el comportamiento de un gran número de fenómenos estudiados por
ejemplo el comportamiento de flexión de viga o el comportamiento de tensión de
deformación que está regido por la ecuación de Poison tanto como en 2 o 3 dimensiones
todos estas ecuaciones tienen una forma particular de sintetizarse y viendo el dominio
por ejemplo que podría ser una estructura monolítica en un numero finito de partes de
ahí viene el nombre del método en contraste con el análisis diferencial donde se
sintetiza en el número infinita de partes con elementos infinitos decimal pequeños el
método de elementos finitos divide el dominio en los elementos más grandes podrás
intuir entonces que más elementos utilicen para dividir la estructura más exacta de todas
maneras la solución exacta es útil para propósitos prácticos, por lo tanto se deberá
escoger una división de malla adecuada a la precisión que se busque en cada caso
particular.
Dependiendo del propósito del análisis una estructura podrá modelarse en base a líneas
simulando vigas o columnas o como placas simulando losas o como elementos
volumétricos que nos arrojen resultados de esfuerzos internos en el análisis estructural
por ejemplo una vez dividido la estructura cada elemento tiene propiedades geométricas
y física mecánicas propias estas propiedades pueden condensarse en matrices de rigidez
elementales que se obtienen a través de aproximaciones de la ecuación diferencial junto
con funciones de formas y las leyes competitivas del material, cada una de estas
matrices contienen todas las propiedades físicas mecánicas y geométricas de un
elemento estas propiedades físicas mecánicas de cada elemento actúan entre sí a través
de los nodos es entonces que de todo los elementos de la estructura van conectándose de
unos a otros para formar la estructura completa para ser analizadas, y de la misma
manera como se conectan los elementos entre sí matemáticamente pasa lo mismo con
las matrices de rigidez individuales de cada elemento conectándose una con otra para
formar la matriz de rigidez global.
En el caso estático lineal la matriz de rigidez global esta multiplicada por los
desplazamientos de los nudos que son la incógnita de la estructura y todo es igualado a
las fuerzas externas aplicadas, todo este proceso genera un sistema lineal de ecuaciones
por lo general muy grande para hacerlo a mano es por eso por lo que este método se
populariza por los ordenadores.
Una vez resuelto el sistema de ecuaciones se tiene los desplazamientos de todos los
nudos de toda la estructura y mediante las ecuaciones de compatibilidad que son las
ecuaciones que relacionan desplazamientos de los nudos con las deformaciones de los
elementos se tiene las deformaciones internas de cada uno de los elementos finalmente
mediante relaciones de mecánica de medio continuo puede en base a las deformaciones
obtenerse los esfuerzos de la estructura y con esto resolver el problema. Cabe mencionar
que en Análisis Estructural lo que se busca son los esfuerzos internos de la estructura.

BREVE HISTORIA DEL MÉTODO DE LOS ELEMENTOS

FINITOS

Aunque el nombre del MEF se ha establecido recientemente, el concepto se ha usado

desde
hace varios siglos. El empleo de métodos de discretizado espacial y temporal y la
aproximación numérica para encontrar soluciones a problemas ingenieriles o físicos es
conocido desde antiguo. El concepto de ‘elementos finitos’ parte de esa idea.
Para encontrar vestigios de este tipo de cálculos podríamos remontarnos a la época de la
construcción las pirámides egipcias. Los egipcios empleaban métodos de discretizado
para determinar el volumen de las pirámides. Arquímedes (287-212 a.C.) empleaba el
mismo método para calcular el volumen de todo tipo de sólidos o la superficie de áreas.
En oriente también aparecen métodos de aproximación para realizar cálculos. Así el
matemático chino Lui Hui (300 d.C.) empleaba un polígono regular de 3072 lados para
calcular longitudes de circunferencias con lo que conseguía una aproximación al
número Pi de 3.1416.
El desarrollo de los elementos finitos tal y como se conocen hoy en día ha estado ligado
al cálculo estructural fundamentalmente en el campo aeroespacial. En los años 40
Courant1 propone la utilización de funciones polinómicas para la formulación de
problemas elásticos en subregiones triangulares, como un método especial del método
variacional de Rayleigh- Ritz para aproximar soluciones.
Fueron Turner, Clough, Martin y Topp2 quienes presentaron el MEF en la forma
aceptada
hoy en día. En su trabajo introdujeron la aplicación de elementos finitos simples (barras
y placas triangulares con cargas en su plano) al análisis de estructuras aeronáuticas,
utilizando los conceptos de discretizado y funciones de forma.
El trabajo de revisión de Oden3 presenta algunas de las contribuciones matemáticas
importantes al MEF. Los libros de Przemieniecki4 y de Zienkiewicz y Holister5
presentan el MEF en su aplicación al análisis estructural. El libro de Zienkiewicz y
Cheung6 o
Zienkiewicz y Taylor7 presenta una interpretación amplia del MEF y su aplicación a
cualquier problema de campos. En él se demuestra que las ecuaciones de los EF pueden
obtenerse utilizando un método de aproximación de pesos residuales, tal como el
método
de Galerkin o el de mínimos cuadrados. Esta visión del problema difundió un gran
interés
entre los matemáticos para la solución de ecuaciones diferenciales lineales y no lineales
mediante el MEF, que ha producido una gran cantidad de publicaciones hasta tal punto
que hoy en día el MEF está considerado como una de las herramientas más potentes y
probadas para la solución de problemas de ingeniería y ciencia aplicada.
Actualmente el método se encuentra en una fase de gran expansión: es ampliamente
utilizado en la industria y continúan apareciendo cientos de trabajos de investigación en
este campo. Los ordenadores han aportado el medio eficaz de resolver la multitud de
ecuaciones que se plantean en el MEF, cuyo desarrollo práctico ha ido caminando
parejo de las innovaciones obtenidas en el campo de la arquitectura de los ordenadores.
Entre éstas, además de permitir la descentralización de los programas de EF, ha
contribuido a favorecer su uso a través de sofisticados paquetes gráficos que facilitan el
modelado y la síntesis de resultados. Hoy en día ya se concibe la conexión inteligente
entre las técnicas de análisis estructural, las técnicas de diseño (CAD), y las técnicas de
fabricación.

CONCEPTOS GENERALES DEL MÉTODO

La idea general del método de los elementos finitos es la división de un continuo en un
conjunto de pequeños elementos interconectados por una serie de puntos llamados
nodos.
Las ecuaciones que rigen el comportamiento del continuo regirán también el del
elemento.
De esta forma se consigue pasar de un sistema continuo (infinitos grados de libertad),
que
es regido por una ecuación diferencial o un sistema de ecuaciones diferenciales, a un
sistema con un número de grados de libertad finito cuyo comportamiento se modela por
un
sistema de ecuaciones, lineales o no.
En cualquier sistema a analizar podemos distinguir entre:
Dominio. Espacio geométrico donde se va ha analizar el sistema.
Condiciones de contorno. Variables conocidas y que condicionan el cambio del
sistema: cargas, desplazamientos, temperaturas, voltaje, focos de calor, ...
Incógnitas. Variables del sistema que deseamos conocer después de que las
condiciones de contorno han actuados sobre el sistema: desplazamientos, tensiones,
temperaturas.

El método de los elementos finitos supone, para solucionar el problema, el dominio

discretizado en subdominios denominados elementos. El dominio se divide mediante
puntos (en el caso lineal), mediante líneas (en el caso bidimensional) o superficies ( en
el
tridimensional) imaginarias, de forma que el dominio total en estudio se aproxime
mediante el conjunto de porciones (elementos) en que se subdivide. Los elementos se
definen por un número discreto de puntos, llamados nodos, que conectan entre si los
elementos. Sobre estos nodos se materializan las incógnitas fundamentales del
problema.
En el caso de elementos estructurales estas incógnitas son los desplazamientos nodales,
ya
que a partir de éstos podemos calcular el resto de las incógnitas que nos interesen:
tensiones, deformaciones, ... A estas incógnitas se les denomina grados de libertad de
cada nodo del modelo. Los grados de libertad de un nodo son las variables que nos
determinan el estado y/o posición del nodo.
Por ejemplo si el sistema a estudiar es una viga en voladizo con una carga puntual en el
extremo y una distribución de temperaturas tal y como muestra la figura,

el
discretizado
del dominio
puede ser:

Los grados de libertad de cada nodo serán:

Desplazamiento en dirección x
Desplazamiento en dirección y
Giro según z

Temperatura
El sistema, debido a las condiciones de contorno: empotramiento, fuerza puntual y
temperatura, evoluciona hasta un estado final. En este estado final, conocidos los
valores
de los grados de libertad de los nodos del sistema podemos determinar cualquier otra
incógnita deseada: tensiones, deformaciones,... También sería posible obtener la
evolución temporal de cualquiera de los grados de libertad.
Planteando la ecuación diferencial que rige el comportamiento del continuo para el
elemento, se llega a fórmulas que relacionan el comportamiento en el interior del mismo
con el valor que tomen los grados de libertad nodales. Este paso se realiza por medio de
unas funciones llamadas de interpolación, ya que éstas ‘interpolan’ el valor de la
variable
nodal dentro del elemento.
El problema se formula en forma matricial debido a la facilidad de manipulación de las
matrices mediante ordenador. Conocidas las matrices que definen el comportamiento
del
elemento (en el caso estructural serán las llamadas matrices de rigidez, amortiguamiento
y masa, aunque esta terminología ha sido aceptada en otros campos de conocimiento) se
ensamblan y se forma un conjunto de ecuaciones algebraicas, lineales o no, que
resolviéndolas nos proporcionan los valores de los grados de libertad en los nodos del
sistema.

Principios generales aplicados a un continuo elástico

A continuación se muestran algunas de las ideas básicas relacionadas con los
fundamentos
matemáticos del MEF aplicadas al caso estructural. En el siguiente capítulo se realiza un
ejemplo con objeto de aclarar las ideas que se muestran en este capítulo.

Ecuaciones de equilibrio. Principio de los Trabajos Virtuales

Muchos problemas de medios continuos vienen expresados mediante ecuaciones
diferenciales y condiciones de contorno sobre la función o funciones incógnita. Ante la
dificultad, y en muchos casos la imposibilidad, de encontrar una solución cerrada, se
opta
por realizar una aproximación, siendo necesaria la expresión integral del Principio de
los
Trabajos Virtuales (PTV).

Se considera un continuo elástico como el de la figura sometido a unas fuerzas


superficiales t tx, ty, tz  y a unas fuerzas por unidad de volumen
T
  
{X} = Xx , Xy , Xz T  (las fuerzas por unidad de superficie podrían ser presiones y el
peso propio sería una fuerza por unidad de volumen). El vector desplazamientos lo
notamos por {u} u,v,wT . Las deformaciones correspondientes a estos
desplazamientos son:

y las tensiones debidas a estas deformaciones serán

Las ecuaciones de equilibrio para un elemento diferencial de volumen pueden

escribirse de la forma
{s}{X} 0
donde el vector {X} incluye de forma general las fuerzas de inercia {X} {X}pu, es
decir, consideramos las fuerzas de inercia como fuerzas por unidad de volumen.
Multiplicando esta ecuación por una función de ponderación {u} e integrando

Utilizando formula de Green se puede escribir

Si se asocia la función de ponderación {u} con un desplazamiento virtual, el operador

actuando sobre él será una deformación virtual.
dedu
El equilibrio en el contorno exige que se cumpla la relación s{n} {t} y
sustituyendo en la expresión.
La formula de Green se obtiene a partir de la relación: (ab) ab ba y del teorema de la divergencia
( adv ands), que nos permite escribir:

Si se asocia la función de ponderación {u} con un desplazamiento virtual, el operador

actuando sobre él será una deformación virtual
dedu
El equilibrio en el contorno exige que se cumpla la relación s{n} {t} y
sustituyendo en la expresión.

La formula
de Green se obtiene a partir de la relación: (ab) ab ba y del teorema de la divergencia

( adv ands), que nos permite escribir:

En la relación anterior es posible introducir la ley de comportamiento de material

mediante
la matriz elástica con las propiedades de éste, [C], de forma que

siendo los vectorese0 T y s0 T las deformaciones y las tensiones iniciales

respectivamente. Introduciendo la expresión (1.3), suponiendo deformaciones y
tensiones
iniciales nulas, en la ecuación (1.2), obtenemos
que constituye la formulación del PTV y relaciona el sistema de cargas real y esfuerzos
con el virtual de desplazamientos.