Método de Rayleigh-Ritz

Cálculo de deformaciones
Solución de estructuras

Análisis Estructural I
Mg. Ing. Civil D. Quintanilla H.
FIC – UNFV
 Esquema
• Potencial total o funcional del problema
• Condiciones de borde
• Configuración admisible
• Principio de energía potencial estacionaria
• Método de Rayleigh-Ritz
 Potencial total o Funcional del problema
Expresión integral que implícitamente contiene ecuaciones diferenciales que describen el
problema. En mecánica estructural el funcional comúnmente usado es la expresión de la
energía potencial (П𝑝 ) [2, p. 69].
La energía potencial total incluye la energía de deformación de distorsión elástica y el
potencial poseído por las cargas aplicadas (por su capacidad de realizar trabajo si se
desplaza una distancia [2, p. 71].
La solución de las ecuaciones diferenciales (forma fuerte) de un problema satisfacen las
condiciones del problema en cada punto. En cambio la solución del funcional (*forma débil)
de un problema satisface las condiciones en una manera promedio [2, p. 69].
El funcional puede ser derivado de la ecuación diferencial del problema para establecer una
ecuación integral para lo cual será necesario usar integración por partes (teorema de la
divergencia). Una revisión final sería que la condición estacionaria en el funcional П𝑝
derivado produzca en efecto la ecuación diferencial que gobierna el problema [1, p. 115].
*Porque (i) reduce o debilita el orden de la continuidad requerida para la función de
desplazamientos y (ii) el cumplimiento de las condiciones de borde de fuerza es hecho de
manera aproximada en lugar de exacta [4, p.32])
 Potencial total o Funcional del problema
П𝑝 = 𝑈 − 𝑊
Donde:
П𝑝 : Potencial total, energía potencial total o funcional del
problema
𝑈: Energía de deformación del sistema
𝑊: Potencial total de las cargas (trabajo externo)
 Energía de deformación: Axial

1 1 2
En una sección: 𝑑𝑈 = 𝜎 𝑑𝜖 𝑈 = 𝜎𝜖 𝑈 = 𝐸𝜖
2 2
𝐿
1 2 1
En toda la barra: 𝑈 = 𝐸𝜖 𝑑𝑣 𝑈= 𝐸𝜖𝑥 2 𝐴𝑑𝑥
𝑉 2 0 2
1 2
*Para el caso general: 𝑈 = 𝐸𝜖 − 𝜖𝐸𝜖𝑜 + 𝜖𝜎𝑜 𝑑𝑣
𝑉 2
 Potencial total: Axial
П𝑝 = 𝑈 − 𝑊

𝐿 𝐿
1
П𝑝 = 𝐸𝜖𝑥 2 𝐴𝑑𝑥 − 𝑞 𝑢 𝑑𝑥 − 𝑃𝑢
0 2 0
Energía de deformación: Flexión
𝐿
1
𝑈= 𝐸𝐼𝑤,𝑥𝑥 2 𝑑𝑥
0 2

Potencial total: Flexión

П𝑝 = 𝑈 − 𝑊
𝐿
1
П𝑝 = 𝐸𝐼𝑤,𝑥𝑥 2 𝑑𝑥 − 𝑞 𝑤 𝑑𝑥 − 𝑃𝑤 − 𝑀𝜃
0 2
Condiciones de borde [2, p. 69]
Son de dos tipos: esenciales (geométricas o cinemáticas) y
naturales (fuerzas)

Esenciales:
En 𝑥 = 0; 𝑤 = 0 y 𝑤,𝑥 = 0
Naturales:
En 𝑥 = 𝐿; 𝑤,𝑥𝑥 = 0 y 𝑤,𝑥𝑥𝑥 = 0
Ya que el momento EI𝑤,𝑥𝑥 y el cortante EI𝑤,𝑥𝑥𝑥 son cero
en 𝑥 = 𝐿
Configuración admisible
Configuración que satisface la compatibilidad interna y las condiciones de borde
esenciales. La curva superior (inadmisible) en la figura tiene cuatro faltas. Las dos
inferiores son admisibles aunque la inferior tiene más sentido físico. Las
configuraciones admisibles no necesitan cumplir condiciones de borde naturales.
Entonces en 𝑥 = 𝐿, ninguna de las dos curvas inferiores necesitan cumplir que
𝑤,𝑥𝑥 = 0 ni que 𝑤,𝑥𝑥𝑥 = 0 [2, p. 70].

Las condiciones de borde de fuerza son siempre incluidas como potencial en la

expresión para П. Esto quiere decir éstas están implícitamente contenidas en П. Sin
embargo, las condiciones de borde esenciales se fijan de manera separada [1,
p.115].
Principio de energía potencial estacionaria
El principio de la energía potencial estacionaria afirma que:
Entre todas las configuraciones admisibles de un sistema conservativo,
aquellas que satisfacen las ecuaciones de equilibrio hacen que la energía
potencial sea estacionaria (mínimo) con respecto a pequeñas variaciones
admisibles de desplazamientos.
El principio es aplicable incluso si la relación carga desplazamiento es
no lineal. Si la condición estacionaria es un mínimo relativo, el
equilibrio es estable. Notar que las cargas se mantienen constantes
mientras que se varían los desplazamientos [2, p. 71].
Mediante la minimización de la energía potencial, minimizamos (en
cierto sentido) la violación de los requerimientos de equilibrio interno
y la violación de las condiciones de borde de fuerzas [1, p. 120].

𝛿П𝑝 = 0
Principio de energía potencial estacionaria
Resorte lineal con carga axial (Tomado de [2, p. 71])

П𝑝 = 𝑈 − 𝑊
П𝑝 = 1 2 𝑘𝑢2 − 𝑃𝑢

𝛿П𝑝 = 0
𝑘𝑢 − 𝑃 = 0
𝑃
𝑢=
𝑘
Principio de energía potencial estacionaria
Para k = 10 tnf/m y P = 50 tnf
 Método de Rayleigh-Ritz
El comportamiento del continuo es descrito por ecuaciones
diferenciales parciales. A excepción de problemas simples hay
poca esperanza de encontrar el campo de esfuerzos o
desplazamientos que resuelve las ecuaciones diferenciales y
satisface las condiciones de borde. La necesidad de resolver las
ecuaciones diferenciales se puede evitar aplicando el método de
Rayleigh-Ritz a un funcional como П𝑝 que describe el problema.
El resultado es un problema substituto que tiene un número
finito de grados de libertad y es descrito por ecuaciones
algebraicas en lugar de ecuaciones diferenciales. Una solución de
Rayleigh-Ritz es raramente exacta pero llega a ser más precisa si
se usan más grados de libertad [2, p.79].
 Método de Rayleigh-Ritz
El método opera sobre el funcional del problema. Se sustituye
las funciones de prueba 𝜑 (funciones 𝑢, 𝑣 o 𝑤 para nuestro caso)
en el funcional П𝑝 y se genera n ecuaciones simultáneas para los
parámetros 𝑎𝑖 usando la condición estacionaria de П𝑝 , 𝛿П𝑝 = 0
[1, p.119].
𝑙 𝑚 𝑛

𝑢= 𝑎𝑖 𝑓𝑖 𝑣= 𝑎𝑖 𝑓𝑖 𝑤= 𝑎𝑖 𝑓𝑖
𝑖=1 𝑖=𝑙+1 𝑖=𝑚+1
[2, p.79]
𝜕П𝑝
= 0 𝑝𝑎𝑟𝑎 𝑖 = 1, 2, … , 𝑛
𝜕𝑎𝑖
Las funciones de prueba necesitan ser sólo m-veces
diferenciables y no necesitan satisfacer las condiciones de borde
de fuerza [1, p.123]
 Método de Rayleigh-Ritz
Cada una de las funciones 𝑓𝑖 tienen que ser admisibles, esto
quiere decir que deben satisfacer las condiciones de
compatibilidad y las condiciones de borde esenciales. No se
requiere que las funciones 𝑓𝑖 satisfagan las condiciones de borde
no esenciales (pero de hacerlo produce una aproximación más
precisa). Usualmente, pero no necesariamente, 𝑓𝑖 son
polinomios [2, p.79].
Después que las ecuaciones algebraicas han sido resueltas y se
ha encontrado los valores de 𝑎𝑖 , los campos de desplazamientos
están completamente definidos. La diferenciación de los campos
de desplazamientos genera las deformaciones y de acuerdo a las
relaciones esfuerzo deformación se puede obtener los esfuerzos
[2, p.79].
Barra con sección variable empotrada
 Campo de desplazamientos por partes [2, p.88]

𝑢 = 𝑎1 𝑥 + 𝑎2 𝑝𝑎𝑟𝑎 0 ≤ 𝑥 ≤ 𝑥2
𝑢 = 𝑎3 𝑥 + 𝑎4 𝑝𝑎𝑟𝑎 𝑥2 ≤ 𝑥 ≤ 𝑥3

Para que el campo de desplazamientos sea admisible:

Condición de borde 𝑢 0 = 0 → 𝑎2 = 0
Continuidad 𝑎1 𝑥2 = 𝑎3 𝑥2 + 𝑎4 → 𝑎4 = 𝑥2 (𝑎1 − 𝑎3 )

Entonces:
𝑢 = 𝑎1 𝑥 𝑝𝑎𝑟𝑎 0 ≤ 𝑥 ≤ 𝑥2
𝑢 = 𝑎3 𝑥 + 𝑥2 𝑎1 − 𝑎3 𝑝𝑎𝑟𝑎 𝑥2 ≤ 𝑥 ≤ 𝑥3
 Conclusiones [2, p.74]
• Si se usan desplazamientos o rotaciones en los nudos como
variaciones (incógnitas) en los campos de desplazamientos, aplicar
el principio de energía potencial estacionaria proporciona
ecuaciones de equilibrio de fuerzas o momentos en los nudos.
Ejemplo: Si 𝐷𝑖 es el desplazamiento en el nudo (o rotación), la
𝜕П
ecuación 𝑝 𝜕𝐷𝑖 = 0 es una ecuación de equilibrio en el nudo
afirmando que las fuerzas (o momentos) aplicadas al nudo suman
cero en la dirección de 𝐷𝑖 . Incluidos en la suma se encuentran las
cargas aplicadas externamente y las cargas aplicadas internamente
debido a la deformación de los componentes estructurales.
• La indeterminación estática no altera el procedimiento o lo hace
más difícil.
• La energía potencial puede ser escrita con sus términos de manera
matricial.
 Ventajas [1, p.87]
• Las ecuaciones de equilibrio pueden ser generadas
simplemente añadiendo las contribuciones de todos los
elementos a П y recurriendo a la condición estacionaria. Esta
condición genera automáticamente los requerimientos de
interconectividad de elementos.
• Las ecuaciones de equilibrio que gobiernan el sistema pueden
ser generadas de manera sistemática o mecánica.
• Se tiene menos visión física del problema, entonces es
importante interpretar físicamente las ecuaciones de
equilibrio una vez que han sido obtenidas mediante este
enfoque.
Observaciones sobre el método variacional
[1, p.116]
1. Puede proporcionar una manera fácil de construir las
ecuaciones que gobiernan el sistema.
2. Puede conducir directamente a las ecuaciones que gobiernan
el sistema y condiciones de borde.
3. Proporciona una visión adicional del problema y sirve de
revisión independiente de la formulación del problema.
4. Una gran clase de funciones de prueba pueden ser utilizadas.
Estás no necesitan satisfacer las condiciones de borde de fuerzas
porque están implícitamente contenidas en el funcional (ver
sección 3.3.4).
 𝑚−1
Anexo: Problema variacional 𝐶
En el funcional del problema la derivada más alta de una variable
de estado es de orden m. Entonces tenemos un problema
variacional:
𝐶 𝑚−1
El orden de las derivadas de las condiciones de borde esenciales
a lo máximo serán m-1.
El orden de las derivadas de las condiciones de borde de fuerza
serán de m a 2m-1.
El orden de la derivada más alta en la ecuación diferencial que
gobierna el problema es 2m [1, p. 110].
Referencias
[1] K. Bathe, Finite Element Procedures
[2] Cook, Malkus and Plesha, Concepts and
applications of finite element analysis
[3] Zienkiewicz (Vol. 1) El método de los
elementos finitos
[4] Brebbia and Dominguez, Boundary elements
an introductory course.