• ASIGNATURA:

• ANÁLISIS ESTRUCTURAL II
• TEMA:

• El método de elementos finitos.

• Introducción. Conceptos de la teoría de la
elasticidad.
• DOCENTE:
• MCs. ING. RONALD J. SALAS BERROSPI
Desarrollo del Contenido
El método de elementos finitos.
 El método de elementos finitos.

INTRODUCCIÓN AL MÉTODO DE LOS ELEMENTOS FINITOS

SISTEMAS DISCRETOS Y SISTEMAS CONTINUOS
Al efectuar una clasificación de las estructuras, suelen dividirse en discretas o
reticulares y continuas. Las primeras son aquéllas que están formadas por un
ensamblaje de elementos claramente diferenciados unos de otros y unidos en
una serie de puntos concretos, de tal manera que el sistema total tiene forma de
malla o retícula. La característica fundamental de las estructuras discretas es que
su deformación puede definirse de manera exacta mediante un número finito de
parámetros, como por ejemplo las deformaciones de los puntos de unión de
unos elementos y otros. De esta manera el equilibrio de toda la estructura
puede representarse mediante las ecuaciones de equilibrio en las direcciones de
dichas deformaciones.

Estructura reticular discreta y estructura continua

 El método de elementos finitos.

Como contrapartida, en los sistemas continuos no es posible

separar, a priori, el sistema en un número finito de elementos
estructurales discretos. Si se toma una parte cualquiera del sistema,
el número de puntos de unión entre dicha parte y el resto de la
estructura es infinito, y es por lo tanto imposible utilizar el mismo
método que en los sistemas discretos, pues los puntos de unión
entre los distintos elementos, que allí aparecían de manera natural,
no existen ahora.
Las estructuras continuas son muy
frecuentes en ingeniería, como, por
ejemplo: bastidores de máquinas,
carrocerías de vehículos, losas de
cimentación de edificios, vasijas de
reactores, elementos de máquinas
(bielas, poleas, carcasas...), y para su
análisis es necesario disponer de un
método que tenga en cuenta su
naturaleza continua.
 El método de elementos finitos.

Hasta la llegada del Método de los Elementos Finitos (MEF),

los sistemas continuos se abordaban analíticamente, pero
por esa vía sólo es posible obtener solución para sistemas
con geometría muy sencilla, y/o con condiciones de
contorno simples. También se han utilizado técnicas de
diferencias finitas, pero éstas plantean problemas cuando
los contornos son complicados.
Como precursores del MEF debe
citarse a Argyris y Kelsey
(Stuttgart, 1955) y a Turner,
Clough, Martin y Topp (Boeing,
1956), aunque con posterioridad el
número de autores en el campo
del MEF ha sido enorme, siendo
uno de los campos de la ingeniería
a los que más esfuerzos de
investigación se han dedicado.
 El método de elementos finitos.

HIPÓTESIS DE DISCRETIZACIÓN
En una estructura discreta, su deformación viene definida por un número finito
de parámetros (deformaciones y/o giros), que juntos conforman el vector de
deformaciones Δ, y la estructura tiene tantas formas de deformarse como
términos tenga dicho vector. Un medio continuo tiene infinitas formas posibles
de deformarse, independientes unas de otras, ya que cada punto puede
desplazarse manteniendo fijos cualquier número finito de los puntos restantes,
por grande que sea este último. Por lo tanto, la configuración deformada de la
estructura no puede venir dada por un vector finito Δ como el anterior, sino que
es una función vectorial u, que indica cuáles son las deformaciones de cualquier
punto, y que tiene tres componentes escalares.
Esta función es la solución de la ecuación
diferencial que gobierna el problema, y si éste
𝑢 (𝑥, 𝑦, 𝑧) está bien planteado, cumplirá las condiciones
𝑢= 𝑣 (𝑥, 𝑦, 𝑧) … (1.1) de contorno impuestas, pero en principio no
𝑤 (𝑥, 𝑦, 𝑧) puede asegurarse que esta función u tenga
una expresión analítica manejable, ni siquiera
que pueda calcularse. Por lo tanto, la función u
no podrá conocerse en general.
 El método de elementos finitos.

El Método de los Elementos Finitos recurre a la

hipótesis de discretización, que se basa en lo siguiente:

➢ El continuo se divide por medio de líneas o superficies

imaginarias en una serie de regiones contiguas y
disjuntas entre sí, de formas geométricas sencillas y
normalizadas, llamadas elementos finitos.

➢ Los elementos finitos se unen entre sí en un número

finito de puntos, llamados nudos.

➢ Los desplazamientos de los nudos son las incógnitas

básicas del problema, y éstos determinan unívocamente
la configuración deformada de la estructura. Sólo estos
desplazamientos nodales se consideran independientes.
 El método de elementos finitos.

➢ El desplazamiento de un punto cualquiera viene

unívocamente determinado por los desplazamientos de
los nudos del elemento al que pertenece el punto. Para
ello se definen para cada elemento, unas funciones de
interpolación que permiten calcular el valor de cualquier
desplazamiento interior por interpolación de los
desplazamientos nodales. Estas funciones de
interpolación serán de tal naturaleza que se garantice la
compatibilidad de deformaciones necesaria en los
contornos de unión entre los elementos.

➢ Las funciones de interpolación y los

desplazamientos nodales definen unívocamente el
estado de deformaciones unitarias en el interior del
elemento. Éstas, mediante las ecuaciones constitutivas
del material definen el estado de tensiones en el
elemento y por supuesto en sus bordes.
 El método de elementos finitos.

➢ Para cada elemento, existe un sistema de fuerzas

concentradas en los nudos, que equilibran a las tensiones
existentes en el contorno del elemento, y a las fuerzas
exteriores sobre él actuantes.

Los dos aspectos más importantes de esta hipótesis, sobre

los que hay que hacer hincapié son:

➢ La función solución del problema u es aproximada de

forma independiente en cada elemento. Para una
estructura discretizada en varios elementos, pueden
utilizarse funciones de interpolación distintas para cada
uno de ellos, a juicio del analista, aunque deben
cumplirse ciertas condiciones de compatibilidad en las
fronteras entre los elementos.
 El método de elementos finitos.

➢ La función solución es aproximada dentro de cada

elemento, apoyándose en un número finito (y pequeño)
de parámetros, que son los valores de dicha función en
los nudos que configuran el elemento y a veces sus
derivadas.

El proceso de discretización descrito tiene una justificación

intuitiva, pero lo que de hecho se sugiere es la minimización
de la energía potencial total del sistema, para un campo de
deformaciones definido por el tipo de elementos utilizado
en la discretización.

Con independencia de que más adelante se estudien en

detalle, se representan a continuación algunos de los
elementos más importantes.
 El método de elementos finitos.

Elasticidad Unidimensional

Elementos para elasticidad unidimensional

Elasticidad Bidimensional

Elementos para elasticidad bidimensional

 El método de elementos finitos.

Elasticidad Tridimensional

Elementos para elasticidad tridimensional

Elasticidad con Simetría de Revolución

Elemento axisimétrico
 El método de elementos finitos.

Vigas

Elemento Viga
Flexión de placas planas

Elementos Placa Plana

 El método de elementos finitos.

Cáscaras Laminares Curvas

Elemento cáscara curva

FUNCIONES DE INTERPOLACIÓN
Consideremos un elemento finito cualquiera, definido por un
número de nudos n. Para facilitar la exposición se supondrá un
problema de elasticidad plana. Un punto cualquiera del
elemento tiene un desplazamiento definido por un vector u, que
en este caso tiene dos componentes:

𝑢 (𝑥, 𝑦) … (1.2)
𝑢=
𝑣 (𝑥, 𝑦)
 El método de elementos finitos.

Consideremos un elemento finito cualquiera, definido por

un número de nudos n. Para facilitar la exposición se
supondrá un problema de elasticidad plana. Un punto
cualquiera del elemento tiene un desplazamiento definido
por un vector u, que en este caso tiene dos componentes:

𝜹𝒆 = 𝑼𝟏 𝑽𝟏 𝑼𝟐 𝑽𝟐 … 𝑼𝒏 𝑽𝒏 T … (1.3)

En este ejemplo se supone que como deformaciones de los

nudos se emplean sólo los desplazamientos, pero no los
giros, lo cual es suficiente para elasticidad plana, como se
verá más adelante. En otros elementos (p.e. vigas o
cáscaras) se emplean además los giros.
 El método de elementos finitos.

𝜹𝒆 = 𝑼𝟏 𝑽𝟏 𝑼𝟐 𝑽𝟐 … 𝑼𝒏 𝑽𝒏 T … (1.3)

El campo de deformaciones en el interior del elemento se

aproxima haciendo uso de la hipótesis de interpolación de
deformaciones:
 El método de elementos finitos.

𝒖 = ෍ 𝑵𝒊 𝑼𝒊 𝒗 = ෍ 𝑵𝒊 𝑽𝒊 … (1.4)

Donde Ni son las funciones de interpolación del

elemento, que son en general funciones de las
coordenadas x,y. Nótese que se emplean las mismas
funciones para interpolar los desplazamientos u y v, y que
ambos desplazamientos se interpolan por separado, el
campo u mediante las Ui y el campo v mediante las Vi. Es
decir que la misma Ni define la influencia del
desplazamiento del nudo i en el desplazamiento total del
punto P, para las dos direcciones x e y.

La interpolación de deformaciones (1.4) puede ponerse en

la forma matricial general:

𝒖 = 𝑵𝜹𝒆 … (1.5)
 El método de elementos finitos.

La matriz de funciones de interpolación N tiene tantas filas

como desplazamientos tenga el punto P y tantas columnas
como grados de libertad haya entre todos los nudos del
elemento.

Las funciones de interpolación son habitualmente

polinomios, que deben poderse definir empleando las
deformaciones nodales del elemento. Por lo tanto, se
podrán usar polinomios con tantos términos como grados
de libertad tenga el elemento.
Para problemas de elasticidad la estructura de esta matriz
es normalmente del tipo:

𝑵𝟏 𝟎 𝑵𝟐 𝟎 … 𝟎 𝑵𝒏 𝟎
𝑵= … (1.6)
𝟎 𝑵𝟏 𝟎 𝑵𝟐 𝟎 … 𝟎 𝑵𝒏
 El método de elementos finitos.

Sin embargo, el aspecto de esta matriz puede ser distinto

para otros elementos, como las vigas o las placas a flexión.

Las funciones de interpolación están definidas únicamente

para el elemento, y son nulas en el exterior de dicho
elemento. Estas funciones tienen que cumplir
determinadas condiciones y aunque éstas se verán en
detalle más adelante, con la expresión anterior se puede
deducir que la función de interpolación Ni debe valer 1 en
el nudo i y 0 en los restantes nudos. Esta condición resulta
evidente si se tiene encuentra que los términos del vector
δe son grados de libertad y por lo tanto son
independientes, y deben poder adoptar cualquier valor
La matriz de funciones de interpolación N tiene tantas filas
como desplazamientos tenga el punto P y tantas columnas
como grados de libertad haya entre todos los nudos del
elemento.
 El método de elementos finitos.

Las funciones de interpolación son habitualmente

polinomios, que deben poderse definir empleando las
deformaciones nodales del elemento. Por lo tanto se
podrán usar polinomios con tantos términos como grados
de libertad tenga el elemento.
Para problemas de elasticidad la estructura de esta matriz
es normalmente del tipo:

Función de interpolación
 El método de elementos finitos.

CRITERIOS DE CONVERGENCIA
Antes de estudiar los criterios para garantizar la convergencia en
el MEF es necesario definir dicho concepto, en el ámbito del
MEF. Se dice que un análisis por el MEF es convergente si al
disminuir el tamaño de los elementos, y por lo tanto aumentar el
número de nudos y de elementos, la solución obtenida
tiende hacia la solución exacta.

Hay que indicar que, en el análisis por el MEF, se introducen,

además de la hipótesis de discretización, otras aproximaciones,
que son fuentes de error en la solución: integración numérica,
errores de redondeo por aritmética finita... El concepto de
convergencia aquí analizado se refiere solamente a la hipótesis
de discretización, prescindiendo de los otros errores, que deben
ser estudiados aparte, y cuyo valor debe en todo caso acotarse.
 El método de elementos finitos.

Las funciones de interpolación elegidas para representar el

estado de deformación de un medio continuo deben satisfacer
una serie de condiciones, a fin de que la solución obtenida por el
MEF, converja hacia la solución real.

Criterio 1
Las funciones de interpolación deben ser tales que cuando los
desplazamientos de los nudos del elemento correspondan a un
movimiento de sólido rígido, no aparezcan tensiones en el
elemento.

Este criterio se puede enunciar también de forma más sencilla:

las funciones de interpolación deben ser capaces de representar
los desplazamientos como sólido rígido, sin producir tensiones
en el elemento.
 El método de elementos finitos.

Por ejemplo, en la figura, los elementos del extremo se

desplazan como un sólido rígido, al no existir tensiones más allá
de la fuerza aplicada.

Deformación de sólido rígido

Empleando la formulación desarrollada más adelante, si se aplican unas
deformaciones en los nudos de valor R que representan un movimiento de
sólido rígido, las deformaciones unitarias en el interior del elemento son:

∈𝑹 = 𝑩𝜹𝒆 … (1.7)
 El método de elementos finitos.

Según este criterio, las tensiones correspondientes deben ser

nulas en todo punto del elemento:

𝝈 = 𝑫 ∈𝑹 = 𝑫𝑩𝜹𝒏 = 𝟎 … (1.8)

Criterio 2
Las funciones de interpolación deben ser tales que cuando los
desplazamientos de los nudos correspondan a un estado de
tensión constante, este estado tensional se alcance en realidad en
el elemento.

Claramente, a medida que los elementos se hacen más pequeños,

el estado de tensiones que hay en ellos se acerca al estado
uniforme de tensiones. Este criterio lo que exige es que los
elementos sean capaces de representar dicho estado de tensión
constante.
 El método de elementos finitos.

Se observa que este criterio de hecho es un caso particular del

criterio 1, ya que un movimiento como sólido rígido (con tensión
nula) es un caso particular de un estado de tensión constante. En
muchas ocasiones ambos criterios se formulan como un sólo
criterio.

A los elementos que satisfacen los criterios 1 y 2 se les llama

elementos completos.
 El método de elementos finitos.

Criterio 3
Las funciones de interpolación deben ser tales que las
deformaciones unitarias que se produzcan en las uniones entre
elementos deben ser finitas. Esto es lo mismo que decir que debe
existir continuidad de desplazamientos en la unión entre
elementos, aunque puede haber discontinuidad en las
deformaciones unitarias (y por lo tanto en las tensiones, que son
proporcionales a ellas).

La siguiente figura ilustra las posibles situaciones, para un caso

unidimensional donde la única incógnita es el desplazamiento u
en la dirección x. En la situación de la izquierda existe una
discontinuidad en el desplazamiento u, que da lugar a una
deformación unitaria infinita: esta situación no está permitida por
el criterio 3. En la situación mostrada a la derecha la deformación
es continua, aunque la deformación unitaria no lo es: esta
situación está permitida por el criterio 3.
El método de elementos finitos.

Criterio de convergencia 3 en una dimensión

 El método de elementos finitos.

Criterio 3
Las funciones de interpolación deben ser tales que las
deformaciones unitarias que se produzcan en las uniones entre
elementos deben ser finitas. Esto es lo mismo que decir que debe
existir continuidad de desplazamientos en la unión entre
elementos, aunque puede haber discontinuidad en las
deformaciones unitarias (y por lo tanto en las tensiones, que son
proporcionales a ellas).

La siguiente figura ilustra las posibles situaciones, para un caso

unidimensional donde la única incógnita es el desplazamiento u
en la dirección x. En la situación de la izquierda existe una
discontinuidad en el desplazamiento u, que da lugar a una
deformación unitaria infinita: esta situación no está permitida por
el criterio 3. En la situación mostrada a la derecha la deformación
es continua, aunque la deformación unitaria no lo es: esta
situación está permitida por el criterio 3.
 El método de elementos finitos.

Este criterio debe cumplirse para poder calcular la energía

elástica U almacenada en toda la estructura, como suma de la
energía de todos los elementos.
𝟏 𝟏
𝑼 = න 𝝈 𝜺𝒅𝒗 = ෍ න 𝝈𝑻 𝜺𝒅𝒗 + 𝑼𝒄𝒐𝒏𝒕 … (1.9)
𝑻
𝟐 𝒗 𝟐
𝒆 𝒗𝒆
Donde el sumando Ucont representa la energía elástica
acumulada en el contorno entre los elementos, que debe ser
nula.

Si este requerimiento no se cumple, las deformaciones no son

continuas y existen deformaciones unitarias infinitas en el
borde entre elementos. Si la deformación unitaria en el
contorno es infinita, la energía que se acumula en él es:
𝟏
𝑼𝒄𝒐𝒏𝒕 = න 𝝈𝑻 ∞ 𝒅𝒗 = 𝒊𝒏𝒅𝒆𝒕𝒆𝒓𝒎𝒊𝒏𝒂𝒅𝒐 … (1.10)
𝟐 𝒗=𝟎
 El método de elementos finitos.

Ya que, aunque el volumen de integración (volumen del

contorno) es nulo, su integral puede ser distinta de cero, con lo
que se almacena energía en los bordes entre elementos, lo cual
no es correcto.

Sin embargo, si la deformación unitaria en el contorno es finita

(aunque no sea continua entre los elementos unidos), la energía
que se acumula es:
𝟏
𝑼𝒄𝒐𝒏𝒕 = න 𝝈𝑻 𝜺𝒄𝒐𝒏𝒕 𝒅𝒗 = 𝟎 … (1.11)
𝟐 𝒗=𝟎

En el caso plano o espacial este requerimiento aplica a la

componente del desplazamiento perpendicular al borde entre
elementos, ya que ésta es la única que acumula energía.
 El método de elementos finitos.

Para elementos de tipo celosía o viga, este requerimiento es fácil de cumplir

pues la unión entre elementos se hace en puntos discretos, y se usan los
mismos desplazamientos y giros para todos los elementos que se unen en un
nudo.

Para elasticidad plana la ecuación diferencial es de orden m=2, con lo que

energía potencial es de orden n=1. En efecto esta última se expresa en
términos de 𝜀, que son las derivadas primeras de las deformaciones.
Luego las funciones de interpolación deben ser continuas en los contornos
de tipo C0, es decir no se exige continuidad a la derivada de la función de
interpolación en el contorno del elemento, sino sólo a la propia función.

Para problemas de flexión de vigas y de placas delgadas, la ecuación diferencial

es de orden m=4, luego la energía potencial es de orden n=2. Por lo tanto, las
funciones de interpolación elegidas deben ser continuas C1 en el contorno del
elemento: tanto la función como su derivada primera deben ser continuas.

A los elementos que cumplen este tercer criterio se les llama compatibles.
 El método de elementos finitos.

Resumen de los tres criterios

Los criterios 1 y 2 pueden agruparse de manera más matemática

diciendo que las funciones de interpolación deben permitir
representar cualquier valor constante de su derivada n‐sima en el
interior del elemento (siendo n el orden de la derivada necesaria
para pasar de las deformaciones a las deformaciones unitarias, que
es el orden de derivación de las deformaciones en el potencial).

Esto puede comprobarse mediante el siguiente sencillo

razonamiento: los criterios 1 y 2 obligan a representar cualquier
valor constante (incluido el valor nulo) de la tensión , lo cual
equivale a representar cualquier valor constante (incluido el valor
nulo) de la deformación unitaria . Pero la deformación unitaria es la
derivada n‐sima de la deformación, luego en consecuencia es
necesario poder representar cualquier valor constante de dicha
derivada. Esto se satisface siempre, si se adoptan como funciones
de interpolación, polinomios completos de orden n como mínimo.
 El método de elementos finitos.

Por ejemplo, si se adopta una función lineal N=Ax+B, sólo es

válida para problemas de elasticidad (n=1), ya que se
representa cualquier valor constante de la derivada n‐ sima
dN/dx=A. Sin embargo, no vale para problemas de flexión de
vigas ni de placas delgadas (n=2), pues siempre es d2N/dx2=0,
es decir que no se puede representar cualquier valor constante
de la derivada segunda.

El criterio 3 exige que las deformaciones unitarias sean finitas

en el contorno entre los elementos. Como estas
deformaciones unitarias son las derivadas n‐simas de las
deformaciones, lo que se exige es que haya continuidad de las
deformaciones y sus derivadas hasta orden n-1 en el
contorno del elemento. Esto es equivalente a imponer la
compatibilidad de desplazamientos en el contorno.
 El método de elementos finitos.

Como resumen de los tres

criterios, para problemas de
elasticidad (n=1) es necesario
emplear polinomios completos
de orden 1, con continuidad C0
entre ellos para garantizar la
convergencia. Es suficiente con
usar funciones del tipo lineal,
que aproximan la solución
mediante una línea quebrada,
aunque se produzcan
discontinuidades en las tensiones
entre los elementos, como se Compatibilidad en
indica en la figura elasticidad unidimensional
 El método de elementos finitos.

Para problemas de flexión de

vigas y placas, (n=2) es
necesario emplear como
mínimo polinomios de grado 2,
con continuidad C1 entre ellos,
es decir que hay que garantizar
la continuidad de la flecha y el
giro entre los elementos. En la
práctica, para la flexión de vigas
planas se usan 4 parámetros
para ajustar la solución (flecha y Compatibilidad en flexión de vigas
giro en cada extremo) por lo que
el tipo de funciones empleadas
son polinomios de grado 3.
 El método de elementos finitos.

Con respecto a la velocidad de convergencia se pueden

resumir las siguientes conclusiones de los análisis efectuados.
Si se utiliza una discretización uniforme con elementos de
tamaño nominal h, y se usa para interpolar los
desplazamientos un polinomio completo de grado c (que
representa exactamente variaciones del desplazamiento de
dicho grado), el error local en los desplazamientos se estima
que es del orden 0(hc+1).

Respecto a las tensiones, son las derivadas n‐simas de los

desplazamientos, luego el error en ellas es de orden 0(hc+1-m)
Introducción. Conceptos de la
teoría de la elasticidad.
 Introducción. Conceptos de la teoría de la
elasticidad.

INTRODUCCIÓN ELASTICIDAD
El dominio continuo que se analiza se extiende según una
única dimensión x, teniendo el material un área variable
con dicha coordenada A(x) (ver figura) . Como fuerzas
exteriores pueden actuar:

• Fuerzas por unidad de volumen q, en sentido axial.

• Fuerzas puntuales aplicadas Fci.

El campo de deformaciones es una función u(x) de una sola

variable, que define la deformación axial de un punto
cualquiera del dominio. La deformación unitaria tiene
sólo la componente axial Ɛ=du/dx

Problema de elasticidad unidimensional

Introducción. Conceptos de la teoría de la
elasticidad.

El problema de elasticidad unidimensional no es de gran

aplicación práctica, pero su estudio tiene interés pues permite
conocer las peculiaridades del MEF con ejemplos muy sencillos.

La ecuación diferencial que gobierna el problema se obtiene

aplicando el equilibrio de fuerzas a un elemento diferencial (ver
figura):
𝐅 = 𝐅 + 𝐝𝐅 + 𝐪𝐀𝐝𝐱 … (1.1)
𝒅𝑭
+ 𝒒𝑨 = 𝟎 … (1.2)
𝒅𝒙
Sustituyendo el valor de la fuerza en función de la tensión F = σA, y ésta
en función de la deformación unitaria σ = Eε se obtiene:

𝒅 𝒅𝒖
𝑬𝑨 + 𝒒𝑨 = 𝟎 … (1.3)
𝒅𝒙 𝒅𝒙
Introducción. Conceptos de la teoría de la
elasticidad.

Equilibrio de fuerzas en elasticidad unidimensional

Las condiciones de contorno pueden ser de dos tipos:

• Desplazamiento conocido en algunos puntos u=uc
• Fuerza aplicada conocida en algunos puntos F = Fci →
du
EA = Fci
dx
La energía elástica acumulada en el material es:
𝟐
𝟏 𝟏 𝟐
𝟏 𝒅𝒖
𝑼= න 𝝈𝜺𝒅𝒗 = න 𝑬𝜺 𝒅𝒗 = න 𝑬𝑨 𝒅𝒙 … (1.4)
𝟐 𝟐 𝟐 𝒅𝒙
Introducción. Conceptos de la teoría de la
elasticidad.

El potencial total es:

𝟐
𝟏 𝒅𝒖
𝑰𝑰 = 𝑼 + 𝑽 = න 𝑬𝑨 𝒅𝒙 − න 𝒒𝒖𝑨𝒅𝒙 − ෍ 𝑭𝒄𝒊 𝒖𝒄𝒊 … (1.5)
𝟐 𝒅𝒙

El orden de derivación del desplazamiento en el potencial es 1, por

lo que se requieren funciones de interpolación de tipo C0: sólo se
exige continuidad a la función para garantizar la convergencia del
método.
Introducción. Conceptos de la teoría de la
elasticidad.

Elemento de Dos Nudos

El elemento más simple para este problema tiene dos nudos,
con desplazamientos U1 y U2 en ellos. La interpolación del
desplazamiento dentro del elemento es:

𝒖 = 𝑵𝟏 𝑼𝟏 + 𝑵𝟐 𝑼𝟐 … (1.6)

𝑼𝟏
𝒖 = 𝑵𝟏 𝑵𝟐 = 𝑵𝜹𝝐 … (1.7)
𝑼𝟐

Elemento de dos nudos

Introducción. Conceptos de la teoría de la
elasticidad.

Con dos parámetros nodales, se puede interpolar una línea

recta: u = a x b
Particularizando en los nudos 1 y 2 se obtiene:
𝑼𝟏 = 𝒂𝒙𝟏 + 𝒃 … (1.8)

𝑼𝟐 = 𝒂𝒙𝟐 + 𝒃 … (1.9)

De estas expresiones se obtienen a y b,

𝑼𝟐 − 𝑼𝟏 … (1.10) 𝒙𝟐 𝑼𝟏 − 𝒙𝟏 𝑼𝟐
𝒂= 𝒂= … (1.11)
𝑳 𝑳

Sustituyendo en la expresión inicial y reordenando, se obtiene

la ley de interpolación:
𝒙𝟐 − 𝒙 𝒙𝟐 − 𝒙
𝒂= 𝑼 + 𝑼 … (1.12)
𝒙 𝟐 − 𝒙𝟏 𝟏 𝒙 𝟐 − 𝒙 𝟏 𝟐
Introducción. Conceptos de la teoría de la
elasticidad.

Las funciones de interpolación son líneas rectas que valen 1

en el nudo i y 0 en el otro nudo:

𝒙𝟐 − 𝒙 𝒙𝟐 − 𝒙 … (1.14)
𝑵𝟏 = … (1.13) 𝑵𝟐 =
𝒙 𝟐 − 𝒙𝟏 𝒙𝟐 − 𝒙 𝟏

Funciones de interpolación. Elemento de dos nudos

La deformación unitaria es:

𝒅𝒖 𝒅𝑵𝟏 𝒅𝑵𝟐 𝑼𝟏 … (1.15)
𝜺= = 𝑵𝜹𝝐 =
𝒅𝒙 𝒅𝒙 𝒅𝒙 𝑼𝟐
Introducción. Conceptos de la teoría de la
elasticidad.

Se define la matriz B como:

𝒅𝑵𝟏 𝒅𝑵𝟐 −𝟏 𝟏
𝑩= = … (1.16)
𝒅𝒙 𝒅𝒙 𝑳 𝑳
Siendo L=x2 –x1 la longitud del elemento.
La matriz elástica es sencillamente D = [ E]
La matriz de rigidez resulta ser:
−𝟏
−𝟏 𝟏
𝐊 = න 𝑩𝑻 𝑫𝑩𝒅𝒗 = න 𝑳 𝑬 𝑨𝒅𝒙 … (1.17)
𝟏 𝑳 𝑳
𝑳
Suponiendo E y A constantes se obtiene:
𝑬𝑨 −𝟏 𝟏 𝑬𝑨 −𝟏 𝟏
𝐊= ‫׬‬ 𝒅𝒙 = … (1.18)
𝑳𝟐 𝑳 𝑳 𝑳 𝑳 𝑳

Esta es la matriz de rigidez del elemento de celosía, ya que

este elemento finito de dos nudos coincide con dicho elemento
estructural.
Introducción. Conceptos de la teoría de la
elasticidad.

Elemento con área variable

Los grados de libertad, las funciones de interpolación y la
matriz B son las mismas que en el elemento de área
constante.

Elemento de área variable

La matriz de rigidez es:

−𝟏
−𝟏 𝟏 𝑬 −𝟏 𝟏
𝐊=න 𝑳 𝑬 𝑨𝒅𝒙 = 𝟐 න 𝑨𝒅𝒙 … (1.19)
𝟏 𝑳 𝑳 𝑳 𝑳 𝑳
𝑳
Introducción. Conceptos de la teoría de la
elasticidad.

En la integral se identifica el área media del elemento Am, con

lo que se obtiene:
𝑬𝑨𝒎 −𝟏 𝟏
𝐊= … (1.20)
𝑳 𝑳 𝑳

Se obtiene la misma expresión que para el elemento de área

constante, pero usando el área media del elemento.
 Introducción. Conceptos de la teoría de la
elasticidad.

Tensiones
La tensión en el elemento (no incluyendo el efecto de las
temperaturas) es:
𝐝𝐮 −𝟏 𝟏 𝐔𝟏
𝛔 = 𝐄𝛆 = 𝐄 = 𝐄𝐁𝛅𝛜 = 𝐄 … (1.21)
𝐝𝐱 𝐋 𝐋 𝐔𝟐
𝐄
𝛔 = 𝐔𝟐 = 𝐔𝟏 … (1.22)
𝐋
Se observa que el elemento produce un campo de tensión
constante en todo su interior.

Además, la tensión también es constante si el elemento es de

área variable, en contradicción con la estática, pues lo que
debe ser constante en este caso es la fuerza axial N, y no la
tensión. Esto es debido a la hipótesis efectuada de variación
lineal de la deformación u, que no es correcta para un
elemento de área variable.
Introducción. Conceptos de la teoría de la
elasticidad.

Ley de Hooke
Originalmente formulada para casos de estiramiento
longitudinal, establece que la deformación Ɛ de un material
elástico es directamente proporcional a la fuerza aplicada F:
∆𝐋 𝐅
𝛆= = … (1.23)
𝐋 𝐀𝐄
∆𝐋: Alargamiento longitudinal
𝐋: Longitud original
𝐄:Módulo de Young o módulo de elasticidad

La ley se aplica a materiales elásticos hasta un límite

denominado límite de elasticidad.
Introducción. Conceptos de la teoría de la
elasticidad.

Ley de Hooke Para Resortes

La forma más común de representar matemáticamente la Ley
de Hooke es mediante la ecuación del resorte, donde se
relaciona la fuerza F ejercida por el resorte con la distancia
adicional producida por alargamiento del siguiente modo:
𝑭 = −𝒌∆𝒙 … (1.24)
𝑨𝑬 … (1.25)
Siendo, 𝒌 =
𝑳
Donde k se llama constante del resorte (también constante de
rigidez) y ∆𝒙 es la separación de su extremo respecto a su
longitud natural. La energía de deformación o energía
potencial elástica Uk asociada al estiramiento del resorte viene
dada por la siguiente ecuación
𝟏 𝟐
𝑼𝒌 = 𝒌𝒙 … (1.26)
𝟐
Introducción. Conceptos de la teoría de la
elasticidad.

Para los resortes reales, esta ley anterior y la ecuación de la

energía solo son válidas por debajo de un cierto valor del
cociente de la tensión F/A <σE, tras superar ese límite el
material sufre internamente transformaciones
termodinámicas irreversibles y pierde la capacidad de
recuperar su longitud original al retirar la fuerza aplicada,
persistiendo un remanente de deformación denominado
deformación plástica.

Originalmente la ley se utilizaba solo para resortes sometidos

a tracción, pero también es válida en resortes o materiales
sometidos a compresión.
Introducción. Conceptos de la teoría de la
elasticidad.

Ley de Hooke en Sólidos Elásticos

En la mecánica de sólidos deformables elásticos la distribución
de tensiones es mucho más complicada que en un resorte o
una barra estirada sólo según su eje. La deformación en el
caso más general necesita ser descrita mediante un tensor de
deformaciones mientras que los esfuerzos internos en el
material necesitan ser representados por un tensor de
tensiones.

Estos dos tensores están relacionados por ecuaciones lineales

conocidas por ecuaciones de Hooke generalizadas o
ecuaciones constitutivas que caracterizan el comportamiento
de un sólido elástico lineal. Estas ecuaciones tienen forma
general.
𝑼𝒌 = ෍ 𝑪𝒊𝒋𝒌𝒍 𝜺𝒌𝒍 … (1.27)
𝒌,𝒍
Introducción. Conceptos de la teoría de la
elasticidad.

Funciones de Interpolación
Se trata de definir las funciones de interpolación para
elementos de tipo general, con dos o más nudos. Para ello
resulta muy práctico el utilizar una coordenada Ɛ local al
elemento, de tal forma que ésta siempre varíe de -1 en el
nudo inicial a +1 en el nudo final. En este sistema de
coordenadas local, el elemento siempre es de longitud 2.

Coordenadas locales
Introducción. Conceptos de la teoría de la
elasticidad.

Elemento de dos nudos

Las funciones de interpolación son:
𝟏−𝜺 𝟏+𝜺 … (1.28)
𝑵𝟏 = … (1.23) 𝑵𝟏 =
𝟐 𝟐

Elemento de dos nudos. Funciones de interpolación

Introducción. Conceptos de la teoría de la
elasticidad.

Elemento de tres nudos

La función N2 es una parábola que debe valer cero en las dos
esquinas y 1 en el nudo central 2, luego su ecuación debe ser
del tipo:
𝑵𝟐 = 𝑪(𝟏 − 𝜺)(𝟏 + 𝜺) … (1.29)

Pero en Ɛ=0 debe ser N2(0)=1, de donde se deduce que el coeficiente

C debe valer 1.
Por consiguiente, la función es:

𝑵𝟐 = 𝟏 − 𝜺 𝟏 + 𝜺 = 𝟏 − 𝜺𝟐 … (1.30)

Elemento de tres nudos. Funciones de interpolación

Introducción. Conceptos de la teoría de la
elasticidad.

La función N1 es una parábola que debe valer cero en los

nudos 2 y 3, y debe valer la unidad en el nudo 1 ( ver figura),
luego debe ser del tipo:
𝑵𝟏 = 𝑪𝜺(𝟏 − 𝜺) … (1.31)
Pero en Ɛ=-1 debe ser N1(-1)=1, luego el coeficiente C debe ser –1/2.
Por consiguiente, la función es:
Ɛ 𝟏 − 𝜺 … (1.32)
𝑵𝟏 =
𝟐

Elemento de tres nudos. Función de interpolación de nudo esquina

Introducción. Conceptos de la teoría de la
elasticidad.

Análogamente, para el nudo 3 se obtiene (ver figura):

Ɛ Ɛ−𝟏
𝑵𝟏 = … (1.33)
𝟐
Elemento de 3 nudos.
Interpolación de nudo esquina

Por lo tanto, este elemento tiene una ley parabólica para a la

interpolación del desplazamiento (ver figura).

Elemento de 3 nudos. Campo de deformaciones

 Introducción. Conceptos de la teoría de la
elasticidad.

Elemento de cuatro nudos

Sus funciones de interpolación se obtienen de forma análoga
al elemento de tres nudos.
𝟏 𝟗
𝑵𝟏 = − (𝟏 − 𝜺)(𝟏 − 𝟑𝜺))(𝟏 + 𝟑𝜺) … (1.34) 𝑵 𝟐 = (𝟏 + 𝜺)(𝟏 − 𝜺))(𝟏 + 𝟑𝜺) … (1.35)
𝟏𝟔 𝟏𝟔
La obtención de las funciones de interpolación siguiendo este método
requiere en algunas ocasiones cierta dosis de buena suerte, y fue bautizada
por Zienkiewicz como formulación “serendipity”, inspirándose en la obra de
Horace Walpole “Los Príncipes de Serendip” (1754).

Elemento de cuatro nudos

Elemento de 4 nudos. Elemento de 4 nudos.

Interpolación de nudo esquina Interpolación de nudo interior
Ejercicios
 Ejercicios

La barra de hormigón mostrada en la figura tiene una longitud de

2.80m y una sección transversal rectangular que varía
linealmente desde 60x80 cm a 30x50cm.

a) Hallar el desplazamiento exacto del extremo de la viga al ser

sometida a una carga puntual de P=80Ton
b) Modelar la barra con elementos de sección uniforme con el
área igual a la sección transversal de la barra real en el eje z
del punto medio
c) Calcular el % de error del resultado anterior
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
 Ejercicios

Resolver la siguiente viga de dos tramos usando el método de

elementos finitos. El módulo de elasticidad del material es E =
2x108ton/m2, la viga es de sección constante de 0.2m x 0.5m
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios
Ejercicios