ANÁLISIS ESTRUCTURAL

AVANZADO

INTRODUCCIÓN A LOS ELEMENTOS FINITOS

Método de Rayleigh-Ritz
Ing. Jesús Daniel Villalba Morales, [Link].
Ing. Alejandro Suárez Reyes, [Link].

2022-1
El principio de la mínima energía potencial [1,2]
• La energía potencial total de un cuerpo elástico viene dada por la contribución de la
energía de deformación (U) y la energía potencial de las fuerzas externas (V), así:

 = U +V
• Un cuerpo alcanza su configuración de equilibrio estable, ante las cargas aplicadas,
cuando su energía potencial total es mínima, lo cuál se conoce como “Principio de la
Mínima Energía Potencial”.

• El principio se establece matemáticamente como:

 = 0

Análisis Estructural Avanzado 2

El principio de la mínima energía potencial (II)
Algunas observaciones sobre el principio:
• Se puede aplicar en problemas donde el material tenga comportamientos elásticos.

• Produce los mismos resultados que el método basado en principio de los trabajos
virtuales (residuos ponderados), el cuál es válido para cualquier material.

• El sistema debe ser conservativo. Es decir, que el trabajo realizado por las fuerzas externas
no dependa de su recorrido.

• Matemáticamente, la expresión para la energía potencial total corresponde a un

funcional.

Análisis Estructural Avanzado 3

Energía de deformación
• La energía de deformación es aquella que un cuerpo absorbe cuando se deforma sobre la
acción de un conjunto de cargas y proviene del trabajo que realizan las fuerzas internas
sobre las respectivas deformaciones.

Análisis Estructural Avanzado 4

Energía de deformación (II)
Para materiales elásticos lineales , la energía de deformación se expresa como:

1
U =     dV
T

2V
• Donde  y  corresponde a las tensiones internas y las deformaciones, respectivamente.
La integral se debe realizar sobre el volumen del cuerpo V.

Expandiendo la integral:

U =  ( x x +  y y +  z z +  xy xy +  yz yz +  xz xz )dV

1
2V

Análisis Estructural Avanzado 5

Energía de deformación (III)
Para una estructura formada de elementos reticulares como el mostrado en la figura

L b

𝑏, ℎ ≪ 𝐿

Estructura con
Área transversal
elementos reticulares
[1]
constante

Análisis Estructural Avanzado 6

Energía de deformación (III)
Realizando la integración en la sección transversal, tenemos:

1  M 2 cQ 2 T 2 N 2 
L
U =  + + +   dx
2 0  EI GA GJ EA 

M, Q, T y N son las funciones para el momento, el cortante, la torsión y la fuerza axial a lo

largo del elemento, respectivamente.

E y G son el módulo de elasticidad y el de cortante del material, respectivamente.

A, I y J son el área, el momento de inercia y el momento polar de inercia de la sección. c es

una constante que depende de la forma del área.

Análisis Estructural Avanzado 7

Energía de deformación (III): barras sometidas
únicamente a fuerzas axiales
• La energía de deformación para una barra sometida solamente a fuerzas axiales viene
dada por
1 N2
U=  dx
2 EA
• Existe una relación entre los desplazamientos en la barra (u) y la fuerza axial (N)

N = E  A  u´
• Por lo tanto, U se puede calcular como

1
U =  E  A  (u´) 2 dx
2

Análisis Estructural Avanzado 8

Energía de deformación (IV): vigas
• La energía de deformación para una viga considerando únicamente el efecto del
momento es
1 M2
U=  dx
2 EI
• Existe una relación entre la deflexión de la viga (v) y el momento en la sección (M)

𝑀 = 𝐸𝐼𝑣´´
• Por lo tanto, U se puede calcular como

1
U =  EI (v´´) 2 dx
2

Análisis Estructural Avanzado 9

Energía de deformación (V): apoyos elásticos
• Resorte traslacional
m
1
=  kdi ui2
m
1
Ut _ v =  kdi vi2 Ut _ h
i =1 2 i =1 2

[3]
• Resorte rotacional
[2]

l
1
U r =  krj v´2j
j =1 2

• Base Elástica [4]

Lbe
1
U BE = 
0
2
k BE v 2 dx
Análisis Estructural Avanzado 10
Energía potencial de las fuerzas externas
• La energía potencial de las fuerzas externas equivale al negativo del trabajo que realizan:

V = −  u B dV −  u T  dS − u Pi

T T T

V S i

Donde la primera integral se refiere al trabajo de las fuerzas de cuerpo B. La segunda integral
corresponde al trabajo realizado por fuerzas de superficie T. Finalmente, la sumatoria se refiere al
trabajo realizado por fuerzas puntuales P. u representa los desplazamientos en el cuerpo.

Análisis Estructural Avanzado 11

Energía potencial de las fuerzas externas (II)
• En el cálculo de la energía potencial de las fuerzas externas, el signo negativo se debe a la reducción
de la capacidad de las fuerzas de hacer trabajo que se origina en el desplazamiento que tuvo que
realizar la fuerza.

• Es importante anotar que el sistema se considera conservativo, es decir, el trabajo no depende de la

trayectoria de las fuerzas.

Análisis Estructural Avanzado 12

Fundamentos del método de Rayleigh-Ritz
• El método de Rayleigh-Ritz se utiliza para encontrar una solución aproximada a
problemas cuyo modelo matemático corresponde a un funcional. Tenemos 3 aspectos
importantes en el desarrollo del método:

Análisis Estructural Avanzado 13

Fundamentos del método de Rayleigh-Ritz (II)
ASPECTO 1: Definición del funcional
• Un funcional se expresa como

I = f ( x, y, y ', y '', y ''',...)

• En el caso de sólidos deformables, se escoge el funcional correspondiente

a la energía potencial total:

 = U +V
con el fin de determinar los desplazamientos inducidos en el sólido por un
conjunto de cargas.
Análisis Estructural Avanzado 14
Fundamentos del método de Rayleigh-Ritz (III)
ASPECTO 2: Definición de una función de aproximación
• En el método de Rayleigh-Ritz se pretende encontrar cuál es la función que minimiza el
funcional correspondiente al problema analizado. La figura presenta varias soluciones
candidatas, pero una sola de ellas minimizará el funcional.

u ( x) u1 ( x )

u2 ( x )

u3 ( x )

a x
b
Análisis Estructural Avanzado 15
Fundamentos del método de Rayleigh-Ritz (IV)
ASPECTO 2: Definición de una función de aproximación

• La función seleccionada describe el comportamiento de la variable de interés en el

dominio del problema.

• En el caso de sólidos, la función corresponde a como se desplaza cada punto material del
cuerpo cuando se somete a un determinado conjunto de cargas.

Análisis Estructural Avanzado 16

Fundamentos del método de Rayleigh-Ritz (V)
ASPECTO 2: Definición de una función de aproximación
• El siguiente paso consiste en definir una función de aproximación, es decir, la solución
encontrada no es necesariamente la exacta. Para un problema en una dimensión.

d ( x)  d ( x)
• En un problema de mecánica de sólidos d indica el campo de deformaciones del sólido
para cualquier punto x en el dominio

¿Cómo debe ser la forma de una función de

aproximación?

Análisis Estructural Avanzado 17

Fundamentos del método de Rayleigh-Ritz (V)
ASPECTO 2: Definición de una función de aproximación
• El siguiente paso consiste en definir una función de aproximación “genérica” para la
solución, como, por ejemplo

d ( x ) = A  x3 + B  x 2 + C  x + D

• Al aplicar el cálculo de variaciones se busca determinar la combinación de valores de los

coeficientes desconocidos (A, B, C y D) que dan origen a la función que minimiza el
funcional.

• La utilización de polinomios facilita las operaciones matemáticas involucradas.

¿Y el grado del polinomio?

Análisis Estructural Avanzado 18

Fundamentos del método de Rayleigh-Ritz (V)
• De forma general, la función de aproximación corresponde a una combinación lineal de
funciones  (llamadas funciones de forma), así:
n
d ( x ) =  aii ( x )
i =1
donde n es la cantidad de funciones de forma involucradas y ai son coeficientes que se
denominan parámetros generalizados.

• El cálculo de los valores de los parámetros generalizados lleva a la forma final de la

función de aproximación para la variable de interés (y por ende la solución al problema).

Análisis Estructural Avanzado 19

Fundamentos del método de Rayleigh-Ritz (VI)
• En el caso de que se utilicen 3 parámetros generalizados tendríamos:

d ( x ) = a11 ( x ) + a22 ( x ) + a33 ( x )

• La función de aproximación se escoge de tal forma que en la medida en que se utilice una
mayor cantidad de términos se obtenga una mejor aproximación.

• Por ejemplo:

v( x) = a2 x + a3 x
2 3
 1 = x 2 2 = x3

Análisis Estructural Avanzado 20

Fundamentos del método de Rayleigh-Ritz (VII)
• Las funciones de forma deben ser linealmente independientes y cada una debe satisfacer
las condiciones de contorno esenciales. Por ejemplo,

i ( x1 ) = i ( x2 ) = 0 i = 1, 2,..., n

• Las funciones de aproximación deberán ser continuas hasta el grado m-1, siendo m el
mayor orden de derivadas en el funcional.

• Las funciones que satisfacen las condiciones anteriores son llamadas de funciones
admisibles.

Análisis Estructural Avanzado 21

Fundamentos del método de Rayleigh-Ritz (VIII)
• Las siguientes figuras muestran ejemplos de funciones que se pueden o no utilizar como
funciones de aproximación, respectivamente.

Real Aproximada

ub ub ub
ua ua ua

a b a b a b

Análisis Estructural Avanzado 22

Fundamentos del método de Rayleigh-Ritz (VIII)
ASPECTO 3: Minimización del funcional
• Como se ha comentado, en la versión presentada del método es necesario determinar el
valor de los parámetros generalizados.

• El método utiliza el cálculo de variaciones para determinar una función que haga
estacionario el funcional que describe el problema analizado.

• Se debe entonces sustituir la expresión definida para la función de aproximación en el

funcional de la energía potencial total, tal que:

 ( a1 , a2 ,..., an )

Análisis Estructural Avanzado 23

Fundamentos del método de Rayleigh-Ritz (IX)
• Se obtiene entonces la primera variación del funcional con respecto a los parámetros
generalizados:

  
 =  a1 +  a2 + ... +  an = 0
a1 a2 an

• Dado que las variaciones de los parámetros generalizados podrían tomar cualquier valor,
obtendremos un sistema de ecuaciones homogéneas de la forma:


= 0 ; i = 1, 2..., n
ai

Análisis Estructural Avanzado 24

Fundamentos del método de Rayleigh-Ritz (X)
• El número de ecuaciones que se genera es igual al número de parámetros generalizados
utilizados en la definición de la función.

• Lo anterior permite que se genere un sistema lineal de ecuaciones cuya solución lleva a
determinar los valores de los parámetros generalizados desconocidos.

• Conociendo los parámetros, la función de aproximación para el campo estudiado queda

totalmente definido. Otros campos se pueden derivar a partir de la función anterior

Análisis Estructural Avanzado 25

Fundamentos del método de Rayleigh-Ritz (X)
• Otro aspecto muy importante del método es la división del dominio del problema en sub-
dominios. A partir de esta división se formula finalmente el método de los elementos
finitos.

• Sin embargo, en los ejemplos siguientes, se aplica una versión del método donde existe
un único dominio.

• En caso de que se requiera dividir el dominio como consecuencia de la naturaleza del

problema, se seguirá utilizando la versión mostrada, pero definiendo una función de
aproximación para cada subdominio.

Análisis Estructural Avanzado 26

Esquema de aplicación del Método de Rayleigh-Ritz
para una barra uniaxial
1. Defina la función de aproximación para el desplazamiento en la barra y verifique que
las condiciones de contorno son satisfechas.

u ( x) = a0 + a1 x + a2 x 2 + ... + an x n

u ( x) = a11 ( x ) + a22 ( x ) + ... + ann ( x ) ; i = 1, 2..., n

2. Determine la energía de deformación U de la barra en función de los términos ai.

3. Encuentre el potencial de las fuerzas externas V

Análisis Estructural Avanzado 27

Esquema de aplicación del Método de Rayleigh-Ritz
para una barra uniaxial
4. Calcule la energía potencial total

 = U +V
5. Aplique el principio de la mínima energía potencial


= 0 ; i = 1, 2..., n
ai
6. Resuelva el sistema de ecuaciones generado en término de los valores de los parámetros
generalizados desconocidos.

7. Calcule la expresión para la fuerza axial

N ( x ) = AEu '( x)
Análisis Estructural Avanzado 28
Ejemplo: Barra sometida a una fuerza axial distribuida
uniforme y carga puntual en el extremo
• Material con módulo de Elasticidad E.

• La sección de la barra es constante y tiene

área A.

• El movimiento de la barra es posible a lo

largo del eje de la barra.

Representación:
• q es una carga constante a lo largo de la
barra.

Análisis Estructural Avanzado 29

Ejemplo: Barra sometida a una fuerza axial distribuida
uniforme y carga puntual en el extremo

1. Función de aproximación escogida y verificación de la condición de contorno

u( x) = ao + a1 x + a2 x2
Para u(0) tenemos que
u (0) = ao + a1 (0) + a2 (0) 2 = 0  ao = 0 1 = x
u ( x) = a1 x + a2 x 2 → u´( x) = a1 + 2a2 x 2 = x 2
Análisis Estructural Avanzado 30
Ejemplo: Barra sometida a una fuerza axial distribuida
uniforme y carga puntual en el extremo
2. Cálculo de la energía de deformación

 2 4 2 2
L
1 1
U = AE  ( a1 + 2a2 x ) dx = AEL  a1 + 2 La1a2 + L a2 
2

2 0
2  3 
3. Cálculo de la energía potencial de las fuerzas externas

L 
V = −   qx ( x ) u ( x ) dx + ( − P ) u ( L ) 
0 
L 
= −   q ( a1 x + a2 x ) dx − P ( a1 L + 2a2 L ) 
2 2

o 
 L2   2 L3 
= a1  PL − q  + a2  PL − q 
 2  3
Análisis Estructural Avanzado 31
Ejemplo: Barra sometida a una fuerza axial distribuida
uniforme y carga puntual en el extremo
4. Cálculo de la energía potencial total

 = U +V

5. Aplicación del principio de la mínima energía potencial

 L2
= AELa1 + AEL a2 + PL − q = 0
2

a1 2
 4 L3
= AEL a1 + AEL a2 + PL − q = 0
2 3 2

a2 3 3
Resolviendo el sistema se obtiene
Lq − P
a1 =
AE
q
a2 = −
2 AE
Análisis Estructural Avanzado 32
Ejemplo: Barra sometida a una fuerza axial distribuida
uniforme y carga puntual en el extremo
6. La forma final para la función de aproximación de los desplazamientos es:

𝐿𝑞 − 𝑃 𝑞 2 𝑞𝐿2 𝑥 𝑥 2 𝑃𝑥
𝑢 𝑥 = 𝑥− 𝑥 = 2 − −
𝐴𝐸 2𝐴𝐸 2𝐴𝐸 𝐿 𝐿 𝐴𝐸
Determinación para la fuerza axial en la barra

𝐿𝑞 − 𝑃 𝑞
𝑁 𝑥 = 𝐸𝐴𝑢´ = 𝐸𝐴 − 𝑥 =𝑞 𝐿−𝑥 −𝑃
𝐴𝐸 𝐴𝐸

Análisis Estructural Avanzado 33

Ejemplo 2: Barra sometida a una carga
puntual en el centro
• Material con módulo de Elasticidad E.

• La sección de la barra es constante y tiene área

A.

• El movimiento de la barra es posible a lo largo

del eje de la barra.

Análisis Estructural Avanzado 34

Ejemplo 2: Barra sometida a una carga
puntual en el centro

Función de aproximación escogida y

verificación de la condición de
u1 ( x) = ao + a1 x u2 ( x) = b0 + b1 x
contorno
0 x L L xL
2 2
Las condiciones de contorno para este
problema son que el u1 (0) = 0 → a0 = 0
desplazamiento en los dos
extremos es cero u2 ( L ) = 0 → b0 = −b1L

Análisis Estructural Avanzado 35

Ejemplo 2: Barra sometida a una carga
puntual en el centro
Además, deben cumplirse las ecuaciones de compatibilidad en el punto B

L L
u1   = u2   → a1 = −b1
2 2
Por lo tanto, las funciones de aproximación quedan definidas por:

u1 ( x ) = a1 x → u1´( x ) = a1

u2 ( x ) = −  L − x  b1 =  L − x  a1 → u1´( x ) = −a1

Análisis Estructural Avanzado 36

Ejemplo 2: Barra sometida a una carga
puntual en el centro
2. Cálculo de la energía de deformación
L /2 L
1 1
U = EA  ( a1 ) dx + EA  ( −a1 ) dx
2 2

2 0
2 L /2
1
U = EALa12
2
3. Cálculo de la energía potencial de las fuerzas externas

PLa1
V =−
2

Análisis Estructural Avanzado 37

Ejemplo 2: Barra sometida a una carga
puntual en el centro
4. Cálculo de la energía potencial total

 = U +V

5. Aplicación del principio de la mínima energía potencial

 P
= 0 → a1 =
a1 2 EA
6. Ecuación para el desplazamiento y la fuerza axial
P P P
u1 ( x) = x N1 ( x) = EA =
2 EA 2 EA 2
P P
u2 ( x ) = ( L − x) N 2 ( x) = −
2 EA 2
Análisis Estructural Avanzado 38
Ejemplo 3: Barra sometida a un cargamento
lineal – polinomio de segundo grado
• Material con módulo de Elasticidad E.

• La sección de la barra es constante y tiene área A.

• El movimiento de la barra es posible a lo largo del

eje de la barra.

 x
q ( x) = q0 1 − 
 L • q varía linealmente a lo largo de la barra.

Análisis Estructural Avanzado 39

Ejemplo 3: Barra sometida a un cargamento
lineal – polinomio de segundo grado

Los pasos 1 y 2 para resolver este problema ya fueron realizados en el ejemplo 3.1por lo
tanto serán omitidos.
3. Cálculo de la energía potencial de las fuerzas externas
L
V = −  qx ( x ) u ( x ) dx
0

 x
L
= −  q0 1 −  ( a1 x + a2 x 2 ) dx
o  L
q0 L2
=−  ( 2a1 + a2 L )
Análisis Estructural Avanzado 12 40
Ejemplo 3: Barra sometida a un cargamento
lineal – polinomio de segundo grado
4. Cálculo de la energía potencial total

 = U +V

5. Aplicación del principio de la mínima energía potencial

 L2
= AE ( a1 L + a2 L ) − q0
2
=0
a1 6
  2 4  L3
= AE  a1 L + a2 L3  − q0 =0
a2  3  12

Resolviendo el sistema: Lq0

a1 = 0.41667
AE
q0
a2 = −0.25
Análisis Estructural Avanzado
AE 41
Ejemplo 3: Barra sometida a un cargamento
lineal – polinomio de segundo grado
6. Ecuación final para el desplazamiento en la barra

q0  1 2
u ( x) =  0.41667 Lx − x 
AE  4 

7. Determinación de la ecuación para la fuerza axial:

q0  1   1 
N ( x) = EAu´= EA  0.41667 L − x  = q0  0.41667 L − x 
AE  2   2 

Análisis Estructural Avanzado 42

Ecuación de aplicación del Método de
Rayleigh-Ritz para una viga
1. Defina la función de aproximación para la deflexión en la viga y verifique que las
condiciones de contorno son satisfechas.

v( x) = a11 ( x ) + a22 ( x ) + ... + ann ( x ) ; i = 1, 2..., n

2. Determine la energía de deformación U de la viga en función de los términos ai
3. Encuentre el potencial de las fuerzas externas V
4. Calcule la energía potencial total
 = U +V
5. Aplique el principio de la mínima energía potencial y determine los valores de los
parámetros generalizados.
6. Calcule la expresión para el momento

M ( x ) = EIv ''( x)
Análisis Estructural Avanzado 43
Ejemplo 4: viga en voladizo con carga puntual en el
extremo libre – polinomio de segundo grado

Análisis Estructural Avanzado 44

Ejemplo 4: viga en voladizo con carga puntual en el
extremo libre – polinomio de segundo grado

1. Función de aproximación escogida y verificación de la condición de contorno

v( x) = ao + a1 x + a2 x 2
Para v(0) tenemos que

v(0) = ao + a1 (0) + a2 (0) 2 = 0  ao = 0

v´(0) = a1 + 2a2 (0) 2 = 0  a1 = 0
v( x) = a2 x 2 → v´( x) = 2a2 x → v´´( x) = 2a2
Análisis Estructural Avanzado 45
Ejemplo 4: viga en voladizo con carga puntual en el
extremo libre – polinomio de segundo grado
2. Cálculo de la energía de deformación
L
1
U = EI  ( 2a2 ) dx = 2 EIa22 L
2

2 0

3. Cálculo de la energía potencial de las fuerzas externas

V = −v( L)  P = −a2 L2 P

4. Cálculo de la energía potencial total

 = 2EIa22 L − a2 L2 P

Análisis Estructural Avanzado 46

Ejemplo 4: viga en voladizo con carga puntual en el
extremo libre – polinomio de segundo grado
5. Aplicación del principio de la mínima energía potencial

= 4 EIa2 L − L2 P = 0
a2
Resolviendo la ecuación se obtiene
PL
a2 =
4 EI
6. Ecuación para la deflexión, el momento y el cortante en la viga
PL 2 PL
v( x) = x M ( x) = − → V ( x) = 0
4 EI 2
El valor exacto para la deflexión en el extremo libre es

PL3
v( x) =
3EI
Análisis Estructural Avanzado 47
Ejemplo 5: viga en voladizo con carga puntual en el
extremo libre – polinomio de tercer grado

E,I

Análisis Estructural Avanzado 48

Ejemplo 5: viga en voladizo con carga puntual en el
extremo libre – polinomio de tercer grado

1. Función de aproximación escogida y verificación de la condición de contorno

v( x) = ao + a1 x + a2 x 2 + a3 x3
Para v(0) tenemos que

v(0) = ao + a1 (0) + a2 (0) 2 + a3 (0)3 = 0  ao = 0

v´(0) = a1 + 2a2 (0) + 3a3 (0) = 0  a1 = 0
v( x) = a2 x 2 + a3 x 3 → v´( x) = 2a2 x + 3a3 x 2 → v´´( x ) = 2a2 + 6a3 x
Análisis Estructural Avanzado 49
Ejemplo 5: viga en voladizo con carga puntual en el
extremo libre – polinomio de tercer grado
2. Cálculo de la energía de deformación
L
U = EI  ( 2a2 + 6a3 x ) dx = EI (12a32 L3 + 12a2 a3 L2 + 4a22 L )
1 2 1
2 0 2
3. Cálculo de la energía potencial de las fuerzas externas

V = − P * ( a2 L2 + a3L3 )

4. Cálculo de la energía potencial total

 = U +V

Análisis Estructural Avanzado 50

Ejemplo 5: viga en voladizo con carga puntual en el
extremo libre – polinomio de tercer grado
5. Aplicación del principio de la mínima energía potencial


= EI (12a3 L3 + 6a2 L2 ) − PL3 = 0
a3

= EI ( 6a3 L2 + 4a2 L ) − PL2 = 0
a2

Resolviendo la ecuación se obtiene

P PL
a3 = − a2 =
6 EI 2 EI

Análisis Estructural Avanzado 51

Ejemplo 5: viga en voladizo con carga puntual en el
extremo libre – polinomio de tercer grado
6. Ecuación para la deflexión, el momento y el cortante en la viga

P 3 PL 2 𝑃𝐿 𝑃
v( x) = − x + x 𝑀 𝑥 = −𝐸𝐼 − 𝑥 =𝑃 𝑥−𝐿
6 EI 2 EI 𝐸𝐼 𝐸𝐼

La deflexión en el extremo libre es

P 3 PL 2 PL3
v( x) = − L + L =
6 EI 2 EI 3EI

Análisis Estructural Avanzado 52

Ejemplo 6: viga simplemente apoyada con carga en la
mitad de la luz – función trigonométrica

E,I constante

Análisis Estructural Avanzado 53

Ejemplo6: viga simplemente apoyada con carga en la
mitad de la luz – función trigonométrica

1. Función de aproximación escogida y verificación de la condición de contorno

x   2 x   3 x 
v( x) = a1sen   + a2 sen   + a3 sen  
 L   L   L 

v(0) = 0 v( L) = 0

Análisis Estructural Avanzado 54

Ejemplo 6: viga simplemente apoyada con carga en la
mitad de la luz – función trigonométrica
2. Cálculo de la energía de deformación
2
1
L
 2   x  4
2
 2 x  9
2
 3 x  
U = EI   − 2 a1sen   − 2 a2 sen   − 2 a3 sen    dx
2 0 L  l  L  l  L  l 

U=
( a 2
1 + 16a22 + 81a32 ) EI 4
4 L3
3. Cálculo de la energía potencial de las fuerzas externas

 L/ 2  2 L / 2   3 L / 2  
V = − P   a1sen   + a2 sen   + a3 sen  
  L   L   L 
V = − P  ( a1 − a3 )

Análisis Estructural Avanzado 55

Ejemplo 6: viga simplemente apoyada con carga en la
mitad de la luz – función trigonométrica
4. Cálculo de la energía potencial total

 = U +V
5. Aplicación del principio de la mínima energía potencial

  a1 EI  4 
= −P +  =0
a1  2L 
3

  8a2 EI  4
= =0
a2 L3

  81a3 EI  4 
= P+ =0
a3  2L
3


Análisis Estructural Avanzado 56

Ejemplo 6: viga simplemente apoyada con carga en la
mitad de la luz – función trigonométrica
Resolviendo el sistema se obtiene

2 L3 P 2L3 P
a1 ( x) = 4 a2 ( x) = 0 a3 ( x) =
 EI 81 4 EI
6. Deflexión en el punto medio de la viga

2 L3 P   x  2L P  3 x 
3
v( x) = 4 sen  + sen  
 EI  L  81 4
EI  L 

 L  2L P   L  2L P  3 L  PL
3 3 3
v   = 4 sen  + sen  
 2   EI  2 L  81 EI  2 L  48EI
4

Análisis Estructural Avanzado 57

Ejemplo 7: viga con base elástica

E,I, Kb

Análisis Estructural Avanzado 58

Ejemplo 7: viga con base elástica

Función de aproximación escogida y verificación de la condición de contorno

v( x) = ao + a1 x + a2 x 2 + a3 x3 + a4 x 4
Para v(0)=0, v(L)=0,

v(0) = ao + a1 (0) + a2 (0) 2 + a3 (0)3 + a4 (0) 4 = 0 → ao = 0

v( L) = a1 L + a2 L2 + a3 L3 + a4 L4 = 0 → a1 = −a2 L − a3 L2 − a4 L3
Análisis Estructural Avanzado 59
Ejemplo 7: viga con base elástica
Por lo tanto, la función de aproximación

v( x) = ( x 4 − L3 x ) a4 + ( x3 − L2 x ) a3 + ( x 2 − Lx ) a2

y su segunda derivada es

v´´( x) = 12 x2 a4 + 6 xa3 + a2
2. Cálculo de la energía de deformación

L
U = EI  (12 x a4 + 6 xa3 + 2a2 ) dx
1 2 2

2 0

 (( x − L x ) a4 + ( x − L x ) a3 + ( x − L x ) a2 )
K BasElas L 2
+ 4 3 3 2 2 2
dx
2 0

Análisis Estructural Avanzado 60

Ejemplo 7: viga con base elástica
EI  144 5 2 
U=  L a4 + 12 L3 2
a3 + 4 La2
2
+ 36 L4
a a
4 3 + 16 L3
a a
4 2 + 12 L2
a3 2
a
2  5 
K BasElas  L9 2 8L7 2 L5 2 11L8 5L7 L6 
+  a4 + a3 + a2 + a4 a3 + a4 a2 + a3a2 
2 9 105 30 60 42 10 

3. Cálculo de la energía potencial de las fuerzas externas

( )
L
V = −  w ( x 4 − L3 x ) a4 + ( x 3 − L2 x ) a3 + ( x 2 − L2 x ) a2 dx
0

3 1 4 1 3
V = wL a4 + wL a3 + qL a2
5

10 4 6

Análisis Estructural Avanzado 61

Ejemplo 7: viga con base elástica
4. Cálculo de la energía potencial total

 = U +V

5. Aplicación del principio de la mínima energía potencial

  
=0 =0 =0
a2 a3 a4

Análisis Estructural Avanzado 62

Ejemplo 7: viga con base elástica
Se obtiene el siguiente sistema de ecuaciones

 144 L9 11L8 5 L7   3 5
 5 EIL + 9 K B 18 EIL + 120 K B 8 EIL + −
5 4 3
KB  wL
84  10 
   a4   
 11L8 8 L7 L 6
   1 4 
18 EIL + 120 K B 12 EIL + 105 K B 6 EIL + K B   a3  = − wL
4 3 2

20  4 
   a2   
 5 L7 L6 
5
L 1
 − wL 
 8 EIL + 84 K B 6 EIL + 4 EIL + K B 
3 2 3
KB
 20 30   6 

A partir de la solución del sistema se obtienen los valores de las incógnitas y

posteriormente puede ser establecida la ecuación para la deflexión en la viga.

Análisis Estructural Avanzado 63

Esquema de aplicación del Método de Rayleigh-Ritz
para un problema de la mecánica de sólidos
A continuación, se establece un esquema que permite utilizar el funcional dado por la
energía potencial total y el método de Rayleigh-Ritz para determinar una función que
aproxime el campo de desplazamientos de un sólido sometido a una combinación de
fuerzas externas.

El esquema es el siguiente:

1. Divida el dominio del problema en los subdominios necesarios para el análisis.

Análisis Estructural Avanzado 64

Esquema de aplicación del Método de Rayleigh-Ritz
para un problema de la mecánica de sólidos
2. Defina una función de aproximación para el campo de desplazamiento dentro de cada
uno de los subdominios y verifique que las condiciones de contorno son satisfechas.

u j ( x) = a1 j f1 j ( x ) + a2 j f 2 j ( x ) + ... + anj f nj ( x )
i = 1, 2..., n j = 1, 2,..., m

n indica el número de términos de la función (no necesariamente es el mismo para todas

las funciones) y m el número de subdominios. aij son los parámetros generalizados.

Análisis Estructural Avanzado 65

Esquema de aplicación del Método de Rayleigh-Ritz
para un problema de la mecánica de sólidos
3. Determine la energía de deformación U correspondiente a cada subdominio y calcule el
total para el cuerpo.

4. Encuentre el potencial de las fuerzas externas V para cada subdominio y calcule el total
para el cuerpo.

5. Defina la expresión para el funcional de la energía potencial total:

 = U +V

Análisis Estructural Avanzado 66

Esquema de aplicación del Método de Rayleigh-Ritz
para un problema de la mecánica de sólidos
6. Calcule la derivada del funcional de la energía potencial total con respecto a cada uno de
los parámetros generalizados desconocidos.

7. Ensamble y resuelva el sistema de ecuaciones que se origina a partir de igualar las

derivadas obtenidas en el paso anterior a cero.

8. Obtenga la forma final de las funciones de aproximación del campo de desplazamiento

para cada uno de los subdominios mediante la sustitución de los valores calculados
para los parámetros generalizados.

9. Calcule las expresiones para otras variables de interés (como esfuerzos y deformaciones)
utilizando las funciones obtenidas para el campo de desplazamiento. Y las relaciones
entre variables.

Análisis Estructural Avanzado 67

Introducción a la idea del método del elemento finito
• Considere la viga del ejemplo 7 pero considerando que se tiene una función de
aproximación de tercer grado y dividida en 6 elementos

• 1. Función de aproximación y condiciones de contorno y compatibilidad

vi ( x) = aoi + a1i x + a2i x 2 + a3i x3

vi ( X i ) = v j ( X i ) vi ´( X i ) = v j ´( X i )
5 5
Ecuaciones
v1 ( 0 ) = v6 ( L ) = 0 Ecuaciones

2
Ecuaciones

Análisis Estructural Avanzado 68

Introducción a la idea del método del
elemento finito (II)
2. Cálculo de la energía de deformación

 1 xi+1
6
K BasElas
xi+1

U =   EI  ( v´´i ) dx +  (v )
2 2
dx 

i =1  2 2
i

xi xi 

3. Energía potencial de las fuerzas externas

6 xi+1
V = −  wv dx
i
i =1 xi
4. Energía potencial total

 = U +V
Análisis Estructural Avanzado 69
Introducción a la idea del método del
elemento finito (III)
5. Aplicación del principio de la mínima energía potencial


= 0 i = 1, 2,...,12
ai
6. La solución del sistema anterior permite conocer los valores de los parámetros
generalizados.

Se obtienen las ecuaciones para cada una de las funciones de aproximación en la viga.

7. Se derivan las expresiones para el momento

Análisis Estructural Avanzado 70

Introducción a la idea del método del
elemento finito (III)
El funcional para este problema se define como:

c = U + Q + q + h
2
1  dT   q = −  q*TdS
U=
2V K xx   dV
 dx  s2
1
Q = −  QTdV q = −  h (T − T ) dS
2

2 s3
donde V

q* corresponde a flujo de calor.

h(T-T∞) es la pérdida por convección.
Q es una fuente de calor interna.
Kxx es la conductividad térmica.

Análisis Estructural Avanzado 71

Introducción a la idea del método del
elemento finito (III)
El funcional de energía electroestática viene dado como

e = U + V
  V   V   '
2 2
1
U =    +  dV  q = −  pVdV '
 
2 V '   x   y   V'

Donde
U es la energía interna almacenada en un campo eléctrico en el volumen V’.
Wp es la energía potencial de la densidad de carga r.
V es voltaje.
 Constante eléctrica.

Análisis Estructural Avanzado 72

Ejercicio
Determine la ecuación para la fuerza axial en la barra de la figura utilizando una función de
aproximación de tercer grado. La barra está hecha con un material de módulo de
elasticidad E y tiene una sección transversal constante con área A. La carga sigue una
distribución parabólica

  x 2 
q( x) = q0 1 −   
 L 
 

Análisis Estructural Avanzado 73

Ejercicio
• Muestre la convergencia del desplazamiento en el extremo libre de la barra considerando
de 1 a 4 elementos. La barra tiene sección circular variable con radio dado por la
expresión abajo. El área para cada elemento puede ser calculada como el valor medio
entre el área inicial y el área final. Utilice una función de aproximación de primer grado.

r ( x) = 0.050 − 0.040 x

L = 1m
q = 10 KN / m
E = 100 MPa

Análisis Estructural Avanzado 74

Bibliografía Publicaciones impresas

[1] Paiva, J. Notas de Aulas. Introdução ao método dos elementos finitos. Universidade de Sao Paulo, 2012.
[2] Assan, A. Método dos elementos finitos. Editora Unicamp. 2010.
[3] T. R. Chandrupatla, Belegundu A. D. Introduction to finite elements in engineering, Pearson, 4Ed, 2011.
[4] Kim N.-H. & Sankar B. V. Introduction to finite element analysis and design, John Wiley& Sons, 2009.
[5] Logan, D. A first course in the finite element method. 5 Ed. Cengage Learning, 2012,
[6] Cook, R. D., Malkus, D. S., Plesha, M. E., Witt R. J. Concepts and Applications of Finite Element Analysis,
Wiley, 4Ed, 2002.

Análisis Estructural Avanzado 75

Bibliografía Imágenes

[1]
[Link]
ara_la_construccion.html
.

Análisis Estructural Avanzado 76