ESTRUCTURAS HIPERESTÁTICAS ANALIZADAS POR EL MÉTODO DE
LAS FUERZAS
1. Introducción
El método de fuerza presenta un método poderoso para analizar estructuras
hiperestáticas, con un comportamiento elástico lineal. El método es muy atractivo
porque tiene significado físico claro, que se basa en un procedimiento eficiente y
bien ordenado de cálculo de desplazamientos de estructuras deformables.
El método de las fuerzas se utiliza para resolver estructuras hiperestáticas. Consiste
en “liberar” una cantidad de vínculos (pueden ser externos o internos de la
estructura) hasta llegar a una estructura isostática. El número de vínculos
eliminados será igual al grado de hiperestaticidad de la estructura .
Para reproducir el problema original en esta “nueva” estructura aparecen como
fuerzas externas, además de las cargas que están aplicadas sobre la estructura, las
reacciones o fuerzas internas provocadas por los vínculos liberados.
Estas fuerzas que sustituyen a los vínculos liberados constituyen incógnitas que es
necesario determinar.
A la nueva estructura se le imponen las condiciones en desplazamientos que son
compatibles con los vínculos liberados. De esa manera se reproduce el problema
original y se obtienen ecuaciones que permiten calcular las incógnitas.
2. Energía de deformación por flexión y directa
Las tensiones normales en un estado de flexión y directa puede expresarse de la
forma:
N ( x) M z ( x) M ( x)
x = + y+ y z (1)
A Iz Iy
Por ahora no consideraremos la energía de deformación generada por torsión (que
no aparece en problemas planos y solo aparece en problemas espaciales) ni por
cortante (que es generalmente despreciable en relación a la generada por flexión).
Sabemos que la energía de deformación se puede expresar como:
1 1 1
U=
2V x x dV = x x dV =
2V E
2E V
x2 dV
(2)
1 N ( x) M z ( x)
2
L
M y ( x)
2 E 0 A A
U= + y + z dA dx
Iz Iy
Método de las Fuerzas - Página 1 de 42
En el caso que estamos considerando, donde el origen de los ejes se encuentra en
el baricentro de la sección y los ejes z e y son ejes principales de la sección,
tenderemos que:
1 N 2 ( x) M z 2 ( x) 2 M y ( x) 2
L 2
2 E 0 A A2
U= + y + z dA dx
Iz2 I y2
1 N 2 ( x) M z 2 ( x) M y ( x)
L 2
2 E 0 A
U= + + dx (3)
Iz I y
L L L 2
1 N 2 ( x) 1 M z 2 ( x) 1 M y ( x)
2 0 EA 2 0 EI z 2 0 EI y
U= dx + dx + dx
En definitiva obtenemos (cuando tomamos ejes principales) que la energía de
deformación de la acción combinada de directa y flexión en ambas direcciones, es la
suma de la energía de deformación producida por cada una de las tres
solicitaciones.
3. Energía de deformación por cortante.
La energía de deformación producida por una fuerza cortante V(x) en una barra de
longitud L sección constante de área A , tiene la expresión:
1 1 2
2 V 2 V G
Uc = dV = dV (4)
Aplicando la fórmula de Jouravsky, la tensión rasante tiene la expresión:
V ( x)
=− (5)
I b
Por lo cual:
1 V ( x) 2 2 1 V ( x) 2 2 L V ( x) 2
L
2G V I 2 b 2 2G 0 I 2 A b 2 2G 0 A
Uc = dV = dA dx = dx
(6)
A 2
donde = 2 2 dA
I A b
En la fórmula anterior, α se llama factor de forma de corte y toma diferentes valores
según la forma de la sección transversal, lo calcularemos a continuación para una
sección rectangular de altura h y ancho b.
Método de las Fuerzas - Página 2 de 42
Recordando que el momento de primer orden μ, es igual al área rayada en la figura
anterior por la distancia del eje z a su baricentro, tendremos que:
h 1h b h2
= A yG = b − y y + − y = − y 2 (7)
2 2 2 2 4
Sustituyendo en la fórmula del factor de forma de corte se tiene:
2 2
A 2 h2 h2 2
h /2 h /2
bh b 2 1
= 2 2 dA = 144 2 6 2 b − y 2 dy = 36 5 − y dy
I A b b h 4b − h /2 4 h − h /2
4
1 h5 6 (8)
= 36 5 =
h 30 5
4. Teorema de Castigliano (segundo teorema)
En 1879, Alberto Castigliano, ingeniero italiano de ferrocarriles, publicó un libro en el
que exponía un método para determinar la deflexión o la pendiente en un punto en
una estructura, en una armadura, una viga o un marco.
Este método, conocido como el segundo teorema de Castigliano, solo es aplicable a
las estructuras que tienen una temperatura constante, soportes que no ceden y
respuestas material elástica lineal.
Si debe determinarse el desplazamiento de un punto, el teorema establece que este
es igual a la primera derivada parcial de la energía de deformación en la estructura
con respecto a una fuerza que actúa en el punto y en la dirección del
desplazamiento. De una manera parecida, la pendiente en un punto de una
estructura es igual a la primera derivada parcial de la energía de deformación en
Método de las Fuerzas - Página 3 de 42
estructura con respecto a un momento que actúa en el punto y con la dirección de
rotación.
Considerermos un sólido en equilibrio bajo la acción de un sistema de cargas
externas Fi, como muestra la figura. Sean ui, los desplazamientos de los puntos del
sólido.
El principio de mínima energía potencial complementaria total establece que dada
una estructura con ciertas condiciones de apoyos y un conjunto de fuerzas externas
aplicadas, la distribución de fuerzas internas en equilibrio compatible con los
vínculos de los desplazamientos es aquella que minimiza la energía potencial
complementaria total de la estructura.
Para un material elástico lineal, la energía potencial complementaria total es igual a
la energía potencial total.
n m
( Fi , M i ) = U ( Fi , M i ) − Fu
i i − M j j
i =1 j =1 (9)
U ( Fi , M i ) = U ( Fi , M i ) ( Fi , M i ) = ( Fi , M i )
El segundo Teorema de Castigliano consiste básicamente en las condiciones de
optimalidad del problema de mínimo de la energía potencial total, las cuales
consisten en plantear gradiente nulo:
n
m
Fk
F u + F M
i =1
i i j j = uk =
Fk
U ( Fi , M i )
k j =1
(10)
n
m
M k
i =1
Fi ui +
M k
j =1
M j j = k =
M k
U ( Fi , M i )
Demostramos el segundo Teorema de Castigliano, que se puede enunciar como:
La derivada parcial de la Energía de Deformación respecto de una de las fuerzas, es
igual al desplazamiento de su punto de aplicación medido en la dirección de la
fuerza.
Método de las Fuerzas - Página 4 de 42
5. Método de la carga unitaria.
El proceso para determinar los desplazamientos mediante la aplicación directa del
segundo teorema de Castigliano es engorroso si actúan más de dos cargas en la
estructura. Ello se debe a que la evaluación de energía de deformación y su
cuadrado es muy extensa. Se logra una simplificación considerable en la
determinación de desplazamientos si aplicamos el segundo teorema de Castigliano
antes de integrar la expresión para el momento flector. Para ejemplificar este punto,
consideremos un pórtico sobre el cual solamente son importantes las deformaciones
por flexión; entonces la energía de deformación está dada por la siguiente ecuación:
l
1 M2
2 0 EI
U= dx (11)
Para obtener la deflexión δi correspondiente a la carga Pi, debemos calcular la
derivada parcial de la energía de deformación U con respecto a P i. Por lo cual:
U M2 M M
l l
i = =
Pi Pi 0 2 EI
dx =
0
EI Pi
dx (12)
Observando que el diagrama de momentos es lineal en P i, es decir que si M̅ es el
diagrama de momentos producido por la fuerza P j, y M1 el producido cuando Pj es
unitaria se cumple la siguiente relación:
M ( Pi , q)
M ( Pi , q) = M ( Pi ) + M (q) = PM 1 (1) + M (q ) = M 1 (1) (13)
Pi
1
Se obtiene la ecuación de la carga unitaria:
U
l
MM 1
i = = dx (14)
Pi 0 EI
La deflexión δi correspondiente a la carga Pi, se calcula como la integral del
diagrama del momento flector M debido a las cargas reales sobre la estructura (P i y
q) por el diagrama del momento flector M1 debido a una carga unitaria
correspondiente al desplazamiento deseado.
Se efectúan deducciones similares para los efectos de deformaciones axiales, por
cortante y torsionales, generalizando la ecuación de la carga unitaria a la fórmula:
NN TT
l
MM VV
i = 1 + 1
+ 1 + 1 dx (15)
0
EA EI GA GJ
Observemos que si aplicamos el principio de los trabajos virtuales, considerando
como carga virtual a la carga Pi considerada como unitaria, y como desplazamientos
virtuales a los desplazamientos reales debidos a las cargas del problema real,
obtenemos la ecuación (15), este resultado es importante para más adelante poder
calcular desplazamientos en estructuras hiperestáticas.
Método de las Fuerzas - Página 5 de 42
6. Método de Vereschaguin.
El defecto principal de la determinación de los desplazamientos por el Teorema de
Castigliano consiste en la necesidad de plantear las expresiones analíticas de las
funciones integrando. Sin embargo, cuando la barra consta de tramos rectos de
rigidez constante en cada uno de ellos, la operación de integración se puede
simplificar. Esta simplificación se basa en el hecho de que los diagramas de las
derivadas de los momentos flectores suelen ser linales en los tramos rectos, y se
puede aplicar el siguiente Método de Vereschaguin.
Supongamos que en un tramo de longitud L, se necesita calcular la integral del
producto de dos funciones f1(x) y f2(x), con la condición de que por lo menos una de
la funciones es lineal. Supongamos que f2(x)=mx+n.
l l l
I = f1 ( x) f 2 ( x)dx = n f1 ( x)dx + m x f1 ( x )dx (16)
0 0 0
La primera de estas dos integrales constituye el área limitada por la curva f 1(x), que
llamaremos Ω1, mientras que la segunda constituye el momento estático de esta
área respecto del eje z, que vale Ω1.yg1 , siendo yg1 la coordenada del centro de
gravedad del primer diagrama. Entonces
l l
I = n f1 ( x)dx + m x f1 ( x)dx = n1 + m1 y g1 = 1 ( my g1 + n ) = 1 f 2 ( y g1 ) (17)
0 0
Así pues, en el Método de Vereschaguin, la integración se sustituye por calcular el
producto del área del primer diagrama por la ordenada del segundo (lineal) en el
centro de gravedad del primer diagrama.
A continuación se presenta una tabla con los valores de áreas y posiciones del
baricentro de los diagramas que más aparecen en la práctica.
Método de las Fuerzas - Página 6 de 42
En forma alternativa, puede utilizarse la tabla presentada en el anejo 2.
Método de las Fuerzas - Página 7 de 42
7. Aplicación del Método de la Carga Unitaria.
El pórtico plano ABC representado en la figura, está empotrado en A y soporta una
carga vertical P en el extremo libre C. Las barras AB y AC están rígidamente unidos
en B, y ambas tienen rigidez flexional constante EI, y área A. Determinar la deflexión
vertical δv, a) despreciando la deformación axial y la deformación por corte, b)
discutir la contribución de la deformación axial para el caso de una sección
rectangular de ancho b y altura h, c) discutir la contribución de la deformación por
corte.
a) Los momentos flectores M0 producidos por la carga P se muestran en la siguiente
figura, en la cual los diagramas de momento flector se trazan sobre los lados
traccionados de las barras.
Método de las Fuerzas - Página 8 de 42
La carga unitaria correspondiente a la deflexión horizontal así como su
correspondiente diagrama de momento flector M1 se muestran en la siguiente figura.
Aplicando la fórmula de la carga unitaria junto con la integración de Vereschauguin
se calcula la deflexión en el extremo libre como:
Método de las Fuerzas - Página 9 de 42
P L2 2L L PL ( L + 3H )
L H 2
M 0 M1 M M
vf = dx + 0 1 dx = . + LH . = (18)
0
EI 0
EI EI 2 3 2 3EI
b) Para calcular la deflexión vertical provocada por las fuerzas de directa
construimos el diagrama correspondiente a la carga P
Y el diagrama correspondiente a la carga unitaria
Luego, la deflexión vertical debida a las deformaciones axiales es:
Método de las Fuerzas - Página 10 de 42
H
N 0 N1 P PH
vd = dx = ( H .1) = (19)
0
EA EA EA
Esta cantidad debe sumarse al resultado obtenido en la parte anterior para evaluar
la deflexión vertical total del extremo libre C.
PL2 ( L + 3H ) PH
v = vf + vd = + (20)
3EI EA
Sustituyendo la inercia y el área por sus valores en una sección rectangular se tiene
que:
bh3 4 PL2 ( L + 3H )
I= vf =
12 Ebh3
(21)
PH
A = bh vd =
Ebh
Si hacemos el cociente entre ambas deflexiones se tiene:
vd h
2
H
=
vf L 4 ( L + 3H )
2
(22)
h
h L 1 d f
L
Como la altura de la barra t estructuralmente debe ser mucho menor que su luz L
(t<L/10), el efecto de la deformación axial es extremadamente pequeño comparado
con el de flexión. Por esta razón es práctica común considerar únicamente los
efectos de deformaciones por flexión cuando se analizan pórticos planos y omitir
completamente la contribución de las deformaciones axiales.
c) Para calcular la deflexión vertical provocada por las fuerzas de corte construimos
el diagrama de cortantes correspondiente a la carga P
Método de las Fuerzas - Página 11 de 42
Y el diagrama correspondiente a la carga unitaria
Luego, la deflexión vertical debida a las deformaciones por corte es:
Método de las Fuerzas - Página 12 de 42
L
V0V1 P PL
vc = 2 dx = 2 ( L.1) = 2 (23)
0
GA GA GA
Comparando la flecha debida a flexión con la debida al cortante suponiendo que la
sección es rectangular, se tiene:
36 2 PL
vc = (1 + )
25 Ebh
PL2 12
vf = ( L + 3H ) (24)
3E bh3
vf L2 H 25 1
= 2 1 + 3
vc h L 18 (1 + )
Como en vigas esbeltas L>10h, entonces L2/h2>100 lo que implica que la flecha por
cortante es despreciable frente a la flecha por flexión.
8. Método de la carga ficticia.
Con frecuencia es necesario determinar el desplazamiento en un punto de una
estructura en el que no está aplicada ninguna carga. Un procedimiento que resuelve
este problema es el método de la carga ficticia.
Para aplicar dicho procedimiento, una carga o momento ficticio se aplica a la
estructura en el punto y dirección del desplazamiento que se desea calcular. La
energía de deformación, es determinada, incluyendo el efecto de la carga o las
cargas ficticias al escribir las ecuaciones de equilibrio para determinar las fuerzas
internas, después de lo cual se aplica el segundo teorema de Castigliano de la
manera acostumbrada.
Cuando se aplica esté método para calcular desplazamientos de estructuras deben
considerarse dos sistemas de cargas. El primero consiste en la estructura sometida
a las cargas reales que producen el desplazamiento que se desea calcular. El
segundo sistema (al aplicar el teorema de Castigliano en su variante como Método
de la carga unitaria) consiste en una carga unitaria que actúa sola sobre la
estructura. Dicha carga unitaria, es una carga ficticia que se considera solamente
con el propósito de calcular el desplazamiento δ de la estructura debida a las cargas
reales.
9. Aplicación del Método de la Carga Ficticia.
Podemos utilizar el método de la carga ficticia para determinar el desplazamiento
horizontal que se produce en el extremo libre C del pórtico del ejemplo anterior.
Para ello consideremos una fuerza ficticia unitaria horizontal aplicada en el extremo
libre como se muestra en la siguiente figura.
Método de las Fuerzas - Página 13 de 42
Recordemos que despreciamos la deformaciones axiales, por lo cual solo nos
interesa determinar su diagrama de momentos que llamaremos M2 y es el siguiente:
Aplicando la fórmula de la carga unitaria junto con la integración de Vereschauguin
se calcula el desplazamiento horizontal del extremo libre C como:
H
M 0M 2 PHL2
h = dx = (25)
0
EI 2 EI
10. Teorema de Menabrea
El segundo teorema de Castigliano se puede generalizar para el análisis de
estructuras hiperestáticas. Las modificaciones necesarias consisten en expresar la
energía de deformación no solo como una función de las en fuerzas aplicadas (y / o
momentos) P1, P2,…, Pi, ..... ,Pn , sino también de las p fuerzas redundantes
estáticamente indeterminadas (y / o momentos) X1, X2 ...,Xp.
Método de las Fuerzas - Página 14 de 42
Un posible sistema de tales fuerzas se muestra en la siguiente figura.
El número de soportes seleccionados para mantener el equilibrio estático de la
sección primaria están excluidos de la enumeración de las fuerzas redundantes.
Sobre esta base, la energía de deformación U para un sistema indeterminado se
puede definir como:
U = U ( P1 , P2 ,..., Pn , X1 , X 2 ,..., X j ,... X p ) (26)
Usando esta función y el segundo teorema de Castigliano, se pueden encontrar los
p desplazamientos en los puntos de aplicación de fuerzas redundantes Xj en la
dirección de estas fuerzas.
Imponiendo que estos desplazamientos sean nulos, se tiene:
U ( P1 , P2 ,..., Pn , X 1 , X 2 ,..., X j ,... X p )
= 0, j = 1, 2,..., p (27)
X j
El número p de estas ecuaciones es igual al grado de hiperestaticidad de la
estructura. Al resolver estas ecuaciones simultáneamente, se obtienen las
magnitudes de las fuerzas redundantes.
Probamos el Teorema de Menabrea que se puede enunciar de la siguiente manera:
En una estructura hiperestática con vínculos rígidos sometida a fuerzas externas,
las reacciones hiperestáticas toman valores que hacen mínima la Energía de
deformación.
11. Aplicación del teorema de Menabrea.
Se considera la viga biempotrada de la figura, se quiere determinar sus momentos
de empotramiento, aplicando el teorema de Menabrea.
Método de las Fuerzas - Página 15 de 42
Para eso vamos a calcular la energía de deformación en función de la reacción
superabundante X .
2
U
L L
1 1
U=
2 EI (M
0
o + XM 1 ) dx =
X EI (M
0
o + XM 1 )M 1dx =
1
L
X
L M o M 1dx (28)
=
EI 0
M o M 1dx +
EI 0
M 1M 1dx = 0 X =− 0
L
M M dx
0
1 1
Los diagramas de solicitaciones son para Mo:
Y para M1:
Método de las Fuerzas - Página 16 de 42
Aplicando el Método de Vereschaguin, o la tabla de integrales de distribución de
momentos dada al final del capítulo se tiene que:
L L
4 qL2 L
0 o 1
M M dx = − .
3 8 2
. ; M M dx = L
0
1 1
(29)
qL2
X=
12
El diagrama de momentos será:
12. Método de las Fuerzas
Para aplicar el método de las fuerzas, se selecciona un conjunto de incógnitas
redundantes Xi, a los que llamaremos factores de fuerza desconocidos , cuya
eliminación de la estructura hiperestática la convierte en isostática. Este modelo
isostático se llama el sistema principal. A modo de ilustración, la siguiente figura
muestra algunos de los principales sistemas de una viga hiperestática.
Método de las Fuerzas - Página 17 de 42
Seleccionado un sistema principal, el método de las fuerzas consiste en escribir las
ecuaciones de compatibilidad de desplazamiento en las direcciones de cada fuerza
redundante aplicando el Teorema de Menabrea, y de acuerdo con el procedimiento
de superposición.
Adoptando la notación δ para representar desplazamientos, la próxima figura ilustra
esta superposición en el caso de la viga de la figura anterior, cuando las reacciones
de los apoyos centrales se eligen como redundantes.
En este caso, como los desplazamientos de estos apoyos son nulos, se puede
escribir las ecuaciones de compatibilidad de desplazamientos de la forma:
1 = 10 + 11 X1 + 12 X 2 = 0
(30)
2 = 20 + 21 X1 + 22 X 2 = 0
Entendiendo que los desplazamientos δij son positivos cuando sus direcciones
coinciden con las direcciones positivas elegidas para las fuerzas redundantes X i.
δi0 es el desplazamiento del punto de aplicación de la fuerza redundante X i en su
propia dirección cuando se aplica al sistema principal el cargamento original, en lo
que se llama "estado" E0 , nos referiremos a los esfuerzos y desplazamientos que
ocurren en ese sistema con esa carga.
Continuando en la misma figura, los desplazamientos δ ij con j distinto de cero,
representa el desplazamiento de punto de aplicación de la fuerza redundante X i en
su propia dirección cuando se aplica una fuerza unitaria al sistema principal en el
Método de las Fuerzas - Página 18 de 42
punto de aplicación y dirección de la fuerza redundante Xj, en lo que se llama estado
Ej. En notación análoga, E representa el estado de la estructura original.
La ecuación anterior (7) llamada ecuación canónica, se escribe en forma matricial
como:
11 12 X1
= − 10 ΔX = -δ0 (31)
21 22 X 2 20
Donde Δ es la matriz de flexibilidad, por el teorema de Betti, se tiene que dicha
matriz es simétrica. δ0 son los coeficientes de carga.
Como las fuerzas redundantes son independientes entre sí y el sistema principal
está correctamente vinculado (por ser isostático y por lo tanto no teniendo
desplazamientos de cuerpo rígido), la matriz de flexibilidad no es singular, con una
solución única garantizada.
Una vez que se han determinado estas incógnitas redundantes, los esfuerzos y
desplazamientos en la estructura original se pueden obtener mediante la
combinación lineal:
n
E = E0 + X i Ei (32)
i =1
donde i varía de 1 al número total de fuerzas redundantes.
Alternativamente, conociendo los valores de las fuerzas redundantes, la estructura
originalmente hiperestática se convierte en isostática, permitiendo la determinación
directa de las otras reacciones de apoyo y de cualquier solicitación utilizando las
ecuaciones de la estática.
Los coeficientes δi0 y δij se pueden obtener por el método de la carga unitaria. Por lo
cual los coeficientes δi0 y δij valen:
Método de las Fuerzas - Página 19 de 42
l
Ni N j M i M j VV TT
ij = + + i j + i j dx (33)
0
EA EI GA GJ
con i variando de 1 al número total de fuerzas redundantes y j variando de 0 al
número total de fuerzas redundantes. Ni, Mi, Vi y Ti representan los esfuerzos
solicitantes en el estado Ei, y Nj, Mj, Vj y T representan los esfuerzos solicitantes en
el estado Ej.
En resumen, el método de fuerza tiene la siguiente sistematica:
I. Elección de un sistema estructural isostático, conocido como el sistema principal.
Mediante la eliminación de un conjunto de fuerzas redundantes de la estructura
hiperestática para transformarla en isostática. Estas fuerzas redundantes son las
principales incógnitas a determinar.
II. Cálculo de la matriz de flexibilidad y los coeficientes de carga.
III. Montaje y resolución del sistema de ecuaciones de compatibilidad de
desplazamientos para obtener dichas fuerzas redundantes.
IV. Obtención de esfuerzos finales.
13. Aplicación del Método de las Fuerzas.
Volviendo al ejemplo que estábamos desarrollando tenemos:
Estado E0, en la estructura isostática dibujamos el diagrama de momentos M0.
Estado E1, en la estructura isostática dibujamos el diagrama de momentos M1.
Método de las Fuerzas - Página 20 de 42
Estado E2, en la estructura isostática dibujamos el diagrama de momentos M2.
Calculo de δ10
En este caso debemos integrar el producto de los siguientes diagramas de
momentos:
Método de las Fuerzas - Página 21 de 42
Se tiene:
3l
M 1M 0 ql 4
10 = dx = − (34)
0
EI 4 EI
Calculo de δ20
En este caso debemos integrar el producto de los siguientes diagramas de
momentos:
Se tiene:
3l
M 2M 0 5ql 4
20 = dx = − (35)
0
EI 24 EI
Método de las Fuerzas - Página 22 de 42
Calculo de δ11
En este caso debemos integrar el producto de los siguientes diagramas de
momentos:
Se tiene:
3l
M 1M 1 4l 3
11 = dx = (36)
0
EI 9 EI
Calculo de δ12 y de δ21
En este caso debemos integrar el producto de los siguientes diagramas de
momentos:
Método de las Fuerzas - Página 23 de 42
Se tiene:
3l
M 1M 2 7l 3
12 = 21 = dx = (37)
0
EI 18EI
Calculo de δ22
En este caso debemos integrar el producto de los siguientes diagramas de
momentos:
Se tiene:
3l
M 2M 2 4l 3
22 = dx = (38)
0
EI 9 EI
Solución del sistema de ecuaciones:
Debemos resolver el sistema:
11 12 X 1 10 l 3 8 7 X 1 ql 4 6
= − =
21 22 X 2 20 18EI 7 8 X 2 24 EI 5
13ql
X1 = (39)
8 7 X 1 3ql 6 20
=
7 8 X 2 4 5 X2 = −
ql
10
Aplicando ahora que:
M = M 0 + X 1M 1 + X 2 M 2 (40)
Método de las Fuerzas - Página 24 de 42
Se obtiene el siguiente diagrama de momentos flectores
14. Sistema principal obtenido al introducir rótulas.
En esta otra elección para el sistema principal, son eliminados vínculos internos de
continuidad de rotación de la curva elástica (configuración deformada) de la viga. En
este caso, se introducen dos rótulas en las secciones de los dos apoyos internos.
Las incógnitas hiperestáticas X1 y X2 son momentos flectores asociados con la
continuidad de rotación de la viga en estas secciones, como se muestra en la
siguiente figura.
Siguiendo la metodología del Método de Fuerzas, la solución al problema radica en
determinar los valores que deben tener los momentos flectores X1 y X2 para que,
juntos con la carga actuante, se restablece la continuidad de rotación en la elástica
de la viga. Los mismos pasos mostrados para la solución considerando la opción
sistema principal anterior se realiza para esta opción. Esto se muestra a
continuación:
Estado E0, en la estructura isostática dibujamos el diagrama de momentos M0.
Método de las Fuerzas - Página 25 de 42
δ10 → rotación relativa entre las secciones adyacentes a la rótula asociada con X1
debida a la carga externa en el estado 0.
δ20 → rotación relativa entre las secciones adyacentes a la rótula asociada con X2
debida a la carga externa en el estado 0.
Estado E1, en la estructura isostática dibujamos el diagrama de momentos M1.
δ11 → rotación relativa entre las secciones adyacentes a la rótula asociada con X1
debida a X1 = 1 en el estado 1.
δ21 → rotación relativa entre las secciones adyacentes a la rótula asociada con X2
debida a X1 = 1 en el estado 1.
Estado E2, en la estructura isostática dibujamos el diagrama de momentos M2.
Método de las Fuerzas - Página 26 de 42
δ12 → rotación relativa entre las secciones adyacentes a la rótula asociada con X1
debida a X2 = 1 en el estado 2
δ22 → rotación relativa entre las secciones adyacentes a la rótula asociada con X2
debida a X2 = 1 en el estado 2.
Calculo de δ10
3l
M 1M 0 ql
10 = dx = − (41)
0
EI 24 EI
Calculo de δ20
3l
M 2M 0
20 = dx = 0 (42)
0
EI
Calculo de δ11
3l
M 1M 1 2l
11 = dx = (43)
0
EI 3EI
Calculo de δ12 y δ21
3l
M 1M 2 l
12 = 21 = dx = (44)
0
EI 6 EI
Calculo de δ22
3l
M 2M 2 2l
22 = dx = (45)
0
EI 3EI
Solución del sistema de ecuaciones:
Debemos resolver el sistema:
11 12 X 1 10 l 4 1 X1 ql 3 1
= − 1 4 X = 24 EI 0
21 22 X 2 20 6 EI 2
2
ql
X1 = (46)
4 1 X 1 ql 2 1 15
=
1 4 X 2 4 0 ql 2
X2 = −
60
Se observa que los valores de X1 y X2 corresponden exactamente a los valores de
los momentos flectores en las secciones de los apoyos internos de la viga continua,
como se indica en la sección anterior, por lo tanto, esta opción del sistema principal,
Método de las Fuerzas - Página 27 de 42
como no podía ser de otra manera, da la misma solución que el sistema principal
adoptado en la sección anterior.
Una ventaja de esta segunda opción de sistema principal es la facilidad para dibujar
el diagrama de los momentos flectores finales. En las secciones donde se
introdujeron las rótulas, el valor del momento flector final es el valor de la incógnita
hiperestática correspondiente a cada rótula. La gráfica del diagrama a lo largo de las
barras se obtienen mediante una simple superposición de los diagramas de casos
básicos.
En el primer tramo hay una superposición de un triángulo con una parábola, el
segundo es una superposición de dos triángulos y el tercero es solo un triángulo.
15. Sistema principal obtenido al introducir rótulas en
pórticos.
Esta sección vamos a estudiar un marco externamente isostático, como el
representado en la siguiente figura, de tal manera que, para la obtención del sistema
principal, es necesario eliminar vínculos de continuidad interna.
La figura adjunta muestra una opción para el sistema principal de la estructura en la
que, como hay un contorno cerrado, el grado de hiperestaticidad sería tres sino
estuviera la rótula, al tener la rótula el grado de hiperestaticidad se reduce a dos.
Por lo tanto, este sistema principal se obtiene agregando dos rótulas al marco de la
estructura. Los momentos flectores en las secciones donde se agregaron las rótulas
son las incógnitas hiperestáticas de esta resolución.
Método de las Fuerzas - Página 28 de 42
Cabe señalar que las rótulas se pueden colocar en otros dos puntos diferentes de
los anteriores, siempre que no se encuentren alineados en la misma barra, lo que
produciría inestabilidad al transformarse en un mecanismo. La figura (a) muestra
otra sistema principal válido obtenido al introducir dos rótulas en la estructura. La
figura (b) indica un sistema principal inválido ya que las tres articulaciones están
alineadas en la barra pórtico superior.
Otra observación importante con respecto a la solución que se obtiene usando un
sistema principal que se genera por la introducción de rótulas es que, en general, en
la solución de los casos básicos, es necesaria la descomposición del marco
isostático compuesto en marcos isostáticos simples. En el caso general, esta
descomposición daría lugar a marcos biapoyados, triarticulados o empotrados con
ménsulas. Para el sistema principal adoptado, una posible descomposición sería en
un marco biapoyado y en otro triarticulado, como se muestra en los estados E0 y E1
representados a continuación.
Las interpretaciones físicas de coeficientes de carga y de flexibilidad para esta
opción de sistema principal, se puede hacer genéricamente de la siguiente manera:
Estado E0: actúa solamente la carga externa en el sistema principal, se deben
calcular los coeficientes de carga δi0.
δi0 → rotación relativa entre las secciones adyacentes a la rótula asociada con la
incógnita hiperestática Xi provocada por la solicitación externa en el estado E0
Método de las Fuerzas - Página 29 de 42
La figura adjunta indica la solución del estado E0 de la presente opción para el
sistema principal. Se observa que para resolver este problema isostático es
conveniente descomponer el marco compuesto de la figura anterior en un marco
triarticulado que está soportado por un marco biapoyado con una barra en ménsula
a la izquierda. El marco compuesto es separado en dos partes por las rótulas
asociadas con las incógnitas hiperestáticas X1 y X2. Los apoyos del marco
triarticulado son ficticios, pero sirven para indicar que hay dos fuerzas de conexión
y el orden en que se cargan los marcos simples, en las secciones de conexión de
las rótulas separadas, la parte que contiene el apoyo ficticio es la porción apoyada.
Para resolver el problema, se deben determinar las "reacciones" de apoyo en el
marco triarticulado y aplicar estas reacciones como si fueran cargas actuando en el
marco apoyado.
De hecho, cada par reacción-carga en un apoyo ficticio de la descomposición
representa un esfuerzo interno de conexión en una rótula. En el estado E0 de este
ejemplo solo hay esfuerzos de conexión verticales, como se muestra en la siguiente
figura.
Estado E1 – Actúa solamente la incógnita hiperestatica X1 en el sistema principal.
La solución del estado E1 para esta opción de sistema principal es similar a la
solución de estado E0. La descomposición del marco compuesto en el estado E1 se
muestra en la siguiente figura.
Se deben calcular los coeficientes de flexibilidad δi1
δi1 → rotación relativa entre secciones adyacentes a la rótula asociada con la
incógnita hiperestática Xi provocada por X1 = 1 en el estado E1.
Método de las Fuerzas - Página 30 de 42
Estado E2 – Actúa solamente la incógnita hiperestática X2 en el sistema principal.
La solución del estado E2 para esta opción de sistema principal es similar a la
solución del estado E0. La descomposición del marco compuesto en el estado 2 se
muestra en la siguiente figura.
Se deben calcular los coeficientes de flexibilidad δi2
δi2 → rotación relativa entre secciones adyacentes a la rótula asociada con la
incógnita hiperestática Xi provocada por X2 = 1 en el estado E2.
Luego se resuelve el sistema de ecuaciones de compatibilidad de desplazamientos
para obtener las incógnitas hiperestáticas, que en éste caso serán los valores del
momento flector final en las secciones donde se introdujeron las rótulas.
Hay al menos dos ventajas en obtener el sistema principal agregando rótulas. La
primera es que, en general, la introducción de una rótula resulta en un cálculo más
simple de los términos de carga y coeficientes de flexibilidad. La segunda ventaja es
que en la determinación de los momentos finales de flexión, se obtienen
superponiendo los diagramas de los casos básicos, lo cual es muy sencillo. En los
puntos donde se introdujeron las rótulas, el valor del diagrama de momentos
flectores finales es el mismo de la incógnita hiperestática correspondiente a esa
rótula.
Método de las Fuerzas - Página 31 de 42
16. Energía de deformación por torsión.
La energía de deformación producida por un momento torsor Mt(x) en una barra de
longitud L sección constante de momento de inercia torsional J , tiene la expresión:
G 2
d d
2
2
L L
Ut =
1
(
xy xy + xz xz )
dAdx = − z + y + dAdx (47)
20 o
2 dy
dz
G 2
L
Ut = Jdx
0
2
donde
G2 G2 2
2
M J
L L L
M t = G J Ut = Jdx = Jdx = t dx
0
2 0
2G 0 J 2G
L
(48)
1 M2
U t = t dx
2 0 GJ
Comparación Flexión-Torsión
Compararemos una viga sometida a un momento flector constante Mf y a un
momento torsor constante Mt. En este caso se tiene que:
1 M t2 1 M2
Ut = L, Uf = L (49)
2 GJ 2 EI
Si para comparar hacemos los momentos torsor y flector iguales tendremos que:
Ut EI U t 2(1 + ) I
= = (50)
U f GJ Uf J
Si observamos que I y J son del mismo orden tendremos que las dos energías de
deformación también lo son.
Conclusión
Podemos concluir que para vigas de dimensiones habituales (relación entre la
longitud y la altura del orden de 10 o mayor), sometidas a solicitaciones habituales,
las energías de deformación producidas por la flexión y la torsión son de orden
equivalente. Ambas son a su vez de un orden mayor que la producida por directa y
cortante, que normalmente pueden ser despreciadas.
Método de las Fuerzas - Página 32 de 42
17. Aprovechamiento de las propiedades de simetría en los
cálculos de estructuras hiperestáticas.
Consideremos un pórtico geométricamente simétrico, como el representado en la
figura:
Al calcular estas estructuras, es posible simplificar la solución del problema,
reduciendo el número de factores de fuerza desconocidos X1, X2, …,Xn.
Veamos el caso de solución de un pórtico con cargas simétrica o antisimétrica. Se
entiende por carga simétrica aquella, para la cual todas las cargas exteriores
aplicadas en la parte derecha del pórtico, coinciden con la imagen de las fuerzas
aplicadas a la parte izquierda. Por carga antisimétrica se entiende aquella para la
cual las fuerzas aplicadas a la parte derecha del pórtico son también la imagen de
las fuerzas aplicadas a la parte izquierda, pero de signo contrario.
De manera análoga se clasifican los factores de fuerza interiores, como se estudió
en cursos anteriores, los momentos flectores la fuerza normal son simétricos
mientras que el momento torsor y los cortantes son antisimétricos.
En un pórtico simétrico, en el plano de simetría, los factores de fuerza
antisimétricos correspondientes a una carga exterior simétrica son iguales a cero,
mientras que en el caso de una carga exterior antisimétrica serán iguales a cero los
factores de fuerza antisimétricos.
Método de las Fuerzas - Página 33 de 42
Esto resulta claro, si tenemos en cuenta que en una estructura dada, el diagrama de
momentos flectores correspondientes a factores antisimétricos será también
antisimétrico, mientras que el diagrama correspondiente a factores simétricos,
también será simétrico. Al multiplicar estos diagramas, obtendremos como resultado
cero; mientras que al multiplicar dos diagramas simétricos o antisimétricos,
obtendremos un resultado diferente de cero.
18. Aplicación a un emparrillado simétrico.
Consideremos un ejemplo sencillo de emparrillado. Se trata de un balcón,
compuesto de material constante, secciones de inercias flexionales I1 e I2 y de
inercias torsionales Jt1 y Jt2, donde actúan una carga uniformemente distribuida de
magnitud q perpendiculares al plano. La estructura se indica en la figura.
Nuestro problema, tiene grado de hiperestaticidad tres, por lo cual en principio,
tendríamos que determinar, por ejemplo, las tres incógnitas hiperestáticas que se
muestran en la siguiente figura.
Método de las Fuerzas - Página 34 de 42
Sin embargo, en este caso se observa que el plano que pasa por el punto medio C
perpendicular a la barra BD se convierte en un plano de simetría de la estructura y
de las cargas. Esto va a permitir simplificar el problema.
En primer lugar por tratarse de un emparrillado en cualquier punto de la estructura
solo pueden existir solicitaciones no planas. De esa manera en el punto C podría
haber un cortante perpendicular al plano de la estructura, un momento torsor y un
momento flector que deforma la estructura fuera del plano. Los desplazamientos
posibles son los correspondientes a esa solicitaciones, o sea, un desplazamiento
perpendicular al plano de la estructura, un giro que torsiona la barra y un giro
producido por la flexión fuera del plano.
Si cortamos un pequeño tramo de la barra BD a un lado y otro del punto C
tendremos que las fuerzas de cortante a un lado y otro del punto C debieran ser
ambas en el mismo sentido por razones de simetría del problema. Pero por otro lado
por razones de equilibrio debieran tener sentido contrario. Luego resulta que el
cortante en el punto C deberá ser nulo.
El mismo razonamiento se puede aplicar para el momento torsor. Por simetría debe
tener el mismo sentido pero por equilibrio debe tener sentido contrario y en
consecuencia es nulo también. El momento flector por razones de simetría tendrá
sentido contrario y por razones de equilibrio también. Luego puede ser distinto de
cero.
Resumiendo la anterior se puede decir que el momento torsor y el cortante en C son
nulos y la única solicitación que existe en C es el momento flector.
Luego si tomamos la mitad de la estructura y tomamos como incógnita el momento
flector en C tendremos un problema con una sola incógnita. Si no hubiéramos
realizado consideraciones de simetría el problema hubiera tenido tres incógnitas.
Método de las Fuerzas - Página 35 de 42
A continuación podemos construir los diagramas de momento flector y momento
torsor producidos por las cargas aplicadas (Estado de Cargas 0 produce en la barra
k Mfl0k y Mt0k) y los producidos por un momento flector en C unitario (Estado 1
produce en la barra k Mfl1k y Mt1k).
Estado E0, en la estructura isostática dibujamos los diagrama de momentos
flectores y torsores.
Estado E1 en la estructura isostática dibujamos los diagrama de momentos flectores
y torsores.
Método de las Fuerzas - Página 36 de 42
Calculo de δ10
En este caso debemos integrar el producto de los siguientes diagramas de
momentos
L1 L2
M fl 1 M fl 0 M t1M t 0 1 qL12 1 L1 1 qL12
10 = dx + dx = − . . . .1.2 − . .L2 .1.2 =
0
EI1 0
GJ t 2 EI1 8 3 2 GJ t 2 8
qL13 EI1 L2 (51)
=− 1 + 6 .
24 EI1 GJ t 2 L1
Calculo de δ11
En este caso debemos integrar el producto de los siguientes diagramas de
momentos
L EI L
L1 L2
M fl 1 M fl 1 M t1M t1 L 2L
11 = dx + dx = 1 + 2 = 1 1 + 2 1 . 2 (52)
0
EI1 0
GJ t 2 EI1 GJ t 2 EI1 GJ t 2 L1
Solución de la ecuación:
Debemos resolver el sistema:
10
11 X 1 + 10 = 0 X 1 = −
11
EI1 L2
2
1+ 6
qL GJ t 2 L1 qL2
X1 = 1
= 1 donde
24 1 + 2 EI1 L2 24
GJ t 2 L1
(53)
EI1 L2
1+ 6
GJ t 2 L1
=
EI L2
1+ 2 1
GJ t 2 L1
Aplicando ahora que:
M fl = M fl 0 + X 1 M fl 1
M t = M t 0 + X 1M t1 (54)
Q = Q0 + X 1Q1
Se obtiene los siguientes diagramas de momentos flectores, momentos torsores y
cortantes.
Método de las Fuerzas - Página 37 de 42
19. Cálculo de desplazamientos en estructuras
hiperestáticas.
Supongamos que queremos calcular en una estructura hiperestática, como la de la
figura, el desplazamiento de un punto i en la dirección s.
Si aplicamos el método de la carga ficticia,
Método de las Fuerzas - Página 38 de 42
Debemos aplicar una fuerza unitaria como la representada en la figura anterior, y
obtendríamos la ecuación (55)
NN TT
l
MM VV
i = 1 + 1
+ 1 + 1 dx (55)
0
EA EI GA GJ
Que como ya observamos es la aplicación del principio de los trabajos virtuales,
considerando como carga virtual a la carga Pi (unitaria), y como desplazamientos
virtuales a los desplazamientos reales debidos a las cargas del problema real.
El problema es que el procedimiento es un poco laborjoso, porque deberíamos
calcular los diagramas de solicitaciones generados por la carga unitaria sobre la
estructura hiperestática. Sin embargo como aplicamos el principio de los trabajos
virtuales, alcanza con considerar, como campo de tensiones virtuales, cualquier
campo de tensiones que equilibre la carga unitaria. Es evidente que las
solicitaciones generadas por la carga untaria lo cumplen, pero además lo cumplen
todas las solicitaciones generadas por la carga unitaria sobre cualquier estructura
isostática obtenida a partir de la estructura hiperestática original.
20. Aplicación del cálculo de desplazamientos.
Calcular la flecha máxima e la viga de la figura:
Método de las Fuerzas - Página 39 de 42
Aplicando el método de la carga ficticia, debemos determinar el diagrama de
momentos M1 producido por una carga ficticia, unitaría aplicada en el centro de una
viga isostática como se representa en la siguiente figura:
El diagrama de momentos Mo, fue calculado en la aplicación 11 y es:
Aplicando la fórmula (15), y utilizando las tablas de integrales de distribución de
momentos, se tiene:
l
MM 1 qL4
i = dx = (56)
0
EI 384 EI
21. Cálculo de barras curvilíneas hiperestáticas.
En el cálculo analítico de estructuras hiperestáticas que contienen barras curvilíneas
se recomienda valerse de las ecuaciones canónicas del método de las fuerzas. Sin
embargo, en {este caso los desplazamientos que forman parte de las ecuaciones
canónicas, no pueden calcularse por el método de Vereschaguin, ni por medio de
las tablas que figuran en anejo 2.
Para determinar los desplazamientos de una barra, cuya configuración coincide con
un arco de circunferencia, resulta necesario calcular las integrales de diversas
combinaciones de funciones trigonométricas más simples. Como éstas
combinaciones son bastante típicas, conviene tabular las integrales más frecuentes
en los problemas de éste tipo, las que se presentan en el anejo A.
Método de las Fuerzas - Página 40 de 42
22. Aplicación del cálculo en barras curvilíneas.
Determinar el desplazamiento horizontal del punto A de la estructura de la figura si
la rigidez EI de todos los tramos es constante.
Para calcular el desplazamiento horizontal del punto A, podemos aplicar el método
de la carga ficticia. Para ello, aplicamos una carga unitaria en el punto A como seb
muestra en la siguiente figura:
Entonces el desplazamiento en A tendrá el valor:
l
MM 1
A = dz (57)
0
EI
El momento flector originado por la fuerza P en el tramo AB, será igual a cero. En el
tramo BC, será M=Pz, y en el tramo CD, será M=PR(1+senφ).
El momento originado por la fuerza unitaria en el tramo AC es nulo, mientras que en
el tramo CD, será: M1=-1.R.(1-cosφ). El signo negativo indica que el momento
flector unitario está orientado en dirección opuesta al de M.
El Producto M.M1 en el tramo Ac es igual a cero. Por lo tanto, la integral se lleva a
cabo solamente en el tramo CD. Sustituyendo dz por R.dφ, y utilizando la tabla del
anejo 1 se obtiene δA, el valor del desplazamiento horizontal del punto A:
/2 /2
MM 1 PR3 − 1 PR 3
A = 0
EI
Rd =
EI (1 + sen )(1 − cos ) d = −
0
2 EI
(58)
Método de las Fuerzas - Página 41 de 42
Anejo A.
Método de las Fuerzas - Página 42 de 42