Capítulo 3: Método

de aná lisis de
desplazamiento o
rigidez

135
 Método de Análisis de Desplazamiento

CAPÍTULO 3
Método de análisis de desplazamiento o
de rigidez

3.1 Introducción
El método de desplazamiento es un método potente para el análisis estructuras
sometidas a diferentes solicitaciones externas, siendo especialmente eficaz para el
análisis de estructuras con alto grado de indeterminación estática.
El método de desplazamiento se introdujo en 1826 por C. L. Navier (1785-1836). A.
Clebsch (1833-1872) presentó el método de desplazamiento para el análisis de cerchas
(1862). Este método fue desarrollado por H. Manderla (1880), E. Winkler (1892), y O.
Mohr (1892). La aplicación del método de desplazamiento al análisis de pórticos fue
presentado por A. Bendixen (1914) y fue llamado el método de la pendiente-deflexión.
Más tarde este método fue desarrollado y presentado en su forma clásica canónica por
A. Ostenfeld (1926). Una contribución importante en el desarrollo del método fue
añadido por S. Leites (1936).
Los fundamentos de método de desplazamiento ofrecen un algoritmo conveniente,
unificado y riguroso para el análisis de estructuras entre ellos se encuentran las vigas,
pórticos, armaduras, emparrillados, etc. El método de desplazamiento se puede aplicar
de manera efectiva para el análisis de estructuras sometidas a cargas diferentes,
cambios de temperatura, asentamientos de apoyos y errores de montaje. La aplicación
del método de desplazamiento es eficaz para para la construcción de líneas de
influencia y conduce a un procedimiento simple y clara. Por último, el método de
desplazamiento es una herramienta muy eficaz para los tópicos especiales de análisis
estructural como son el estudio de la estabilidad y la dinámica estructural.

3.2 Descripción del método

1. Determinar el grado de indeterminación cinemática. Luego establecer un sistema de
coordenadas para identificar la ubicación y dirección de los desplazamientos de los

136
 Método de Análisis de Desplazamiento

nudos. En las coordenadas se introducen fuerzas restringentes en igual número que el

grado de indeterminación cinemática para impedir el desplazamiento de los nudos.
Debemos hacer notar que, a diferencia del método de la fuerza, el procedimiento
anterior no exige que se haga una selección respecto a las fuerzas restringentes. Este
hecho favorece el uso del método de rigidez en programas generales de computadora
para el análisis de una estructura.
2. En la estructura restringida se aplican las solicitaciones externas: cargas, variación
de temperatura, asentamiento de apoyos, errores de montaje, etc. (estos efectos
pueden considerarse por separado o juntos); y se determinan las fuerzas restringentes
como una suma de las fuerzas en los extremos fijos que se juntan en el nudo. Para la
mayoría de los casos prácticos, las fuerzas en los extremos fijos se pueden calcular
con ayuda de tablas.
Debemos recordar que las fuerzas restringentes son aquellas que se necesitan
para impedir el desplazamiento en las coordenadas debido a todos los efectos.
Si se va a efectuar el análisis para el efecto del movimiento de los apoyos, por
ejemplo un asentamiento, se incluyen en las fuerzas restringentes las fuerzas en las
coordenadas que se necesitan para mantener el apoyo en la posición desplazada.
También se determinan las fuerzas internas de las barras en los lugares necesarios
con los nudos en la posición restringida.
3. Se supone ahora que la estructura está deformada de tal modo que un
desplazamiento en una de las coordenadas es igual a la unidad y todos los demás
desplazamientos tienen valor cero, y se determinan las fuerzas necesarias para
mantener la estructura en esta configuración. Estas fuerzas se aplican en las
coordenadas que representan los grados de libertad. Se determinan las fuerzas
internas en los lugares necesarios correspondientes a esta configuración. Se repite el
procedimiento para un valor unitario de desplazamiento en cada una de las
coordenadas por separado.
4. Se determinan los valores de los desplazamientos necesarios para eliminar las
fuerzas restringentes introducidas en (2). Esto demanda el uso de ecuaciones de
superposición en que se suman los efectos de los desplazamientos separados sobre
las fuerzas restringentes.

137
 Método de Análisis de Desplazamiento

5. Finalmente, se obtienen las fuerzas sobre la estructura original sumando las fuerzas
aplicadas sobre la estructura restringida a las fuerzas producidas por los
desplazamientos de nudos determinados en (4)
Nomenclatura.
i : Coordenada, número asociado a un grado de libertad
Ri : Fuerza restringente en la coordenada i
Di : Desplazamiento en la coordenada i (incógnitas)
K ij : Fuerza de restricción que surge en la coordenada i producida al dar un
desplazamiento solo a la coordenada j. Si Dj = 1 entonces Kij se conoce como
coeficiente de rigidez.
El uso de del procedimiento anterior se explica mejor con el siguiente ejemplo

Ejemplo 3.1
Construir los diagramas de fuerza cortante y momento flector para la viga continúa
mostrada en la figura 3.1a. Tómese EI constante

Figura 3.1 Viga continúa analizada en el ejemplo 3.1. (a) Viga estáticamente indeterminada (b) sistema
de coordenadas

Considerando que la rigidez axial de la viga es infinita el grado de indeterminación

cinemática es dos, puesto que hay dos desplazamientos posibles en los nudos, como
se muestra en la figura, que ilustra también el sistema de coordenadas seleccionado.

138
 Método de Análisis de Desplazamiento

Las fuerzas restringentes son iguales a la suma de las fuerzas de los extremos en los
nudos, se calculan con ayuda de tablas y se consideran positivas cuando su sentido
coincide con el de las coordenadas. Para ilustrar la relación entre las fuerzas de los
extremos y las fuerzas restringentes consideremos el diagrama de cuerpo libre de las
barras y los nudos (véase figura 3.2). Sobre los nudos A y B actúan las fuerzas
restringentes R1 y R2 y fuerzas iguales y opuestas a la de los extremos de las barras.
R1 y R2 se calculan por equilibrio en la dirección de cada coordenada.

Figura 3.2 (a) Estructura restringida sometida a las cargas externas. (b) Diagramas de cuerpo libre

2 2
∑ M A =0 : R1 + q12L =0 R 1= −q12L
2 2
∑ M B=0 : R2 + q12L −
qL
12
=0 R 2=0

Se impone ahora una rotación unitaria ( D1=1) en la dirección de la coordenada 1

manteniendo nulo el desplazamiento de la coordenada 2 (D2=0) (véase figura 3.3).
Determinamos las fuerzas necesarias (coeficientes de rigidez K 11 y K 21) necesarias
para mantener la estructura en esta configuración.

139
 Método de Análisis de Desplazamiento

Se aplica ahora un giro unitario (D2=1) en la dirección de la coordenada 2

manteniendo nulo el desplazamiento de la coordenada 1 (D1=0) (véase figura 3.4).
Hallamos las fuerzas necesarias (K12 y K22) para mantener la viga en esta configuración

k11 k21

A B C

k11 k21
4EI/L 2EI/L

6EI/L2 6EI/L2

Figura 3.3

4 EI 4 EI
∑ M A =0 :k 11− L
=0 k 11=
L

2 EI 2 EI
∑ M B=0 :k 21− L
=0 k 21=
L

k12 k22

A B C

k12 k22
2EI/L 4EI/L 4EI/L 2EI/L

6EI/L2 6EI/L2 6EI/L2 6EI/L2

Figura 3.4

2 EI 2 EI
∑ M A =0 :k 12− L
=0 k 12=
L
4 EI 4 EI 8 EI
∑ M B=0 :k 22− L
−
L
=0 k 22=
L

140
 Método de Análisis de Desplazamiento

En la estructura real los nudos A y B sufren las rotaciones D1 y D2 y no hay fuerzas

restringentes. Por lo tanto, la suma de las fuerzas restringentes y las fuerzas debido a
los desplazamientos verdaderos a lo largo de cada coordenada debe der igual a cero.
Así obtenemos las ecuaciones de equilibrio que expresan el hecho de que las fuerzas
restringentes desaparecen cuando tienen lugar los desplazamientos D1 y D2

R1 +k 11 D 1+ k 12 D2=0
R2 +k 21 D 1+ k 22 D2=0 } (3.1)

Matriz de rigidez
Las ecuaciones de equilibrio se pueden expresar en forma de matriz

{ R } + [ K ] { D }={0 }
o bien
[ K ] { D }={−R } (3.2)

El vector columna { R } depende de la carga aplicada a la estructura.

Los elementos de la matriz [ K ] son fuerzas que corresponden a valores unitarios de
desplazamientos. Por lo tanto, [ K ] depende de las propiedades de la estructura y
representa su rigidez. Por esta razón, [ K ] se denomina matriz de rigidez y sus
elementos se llaman coeficientes de rigidez
Los elementos del vector {D} son los desplazamientos desconocidos y se calculan
resolviendo las ecuaciones 3-2, es decir,

−1
{ D }= [ K ] {−R } (3.3)

En un caso general, si el número de restricciones que se introducen en la estructura

es n, el orden de las matrices { D } , [ K ] y { R } es nx 1, nxn y nx 1, respectivamente.
Para el caso del ejemplo tenemos:

{ }
2
qL
{ R } = 12 [ K ] =
0
EI 4 2
L 2 8 [ ]
141
 Método de Análisis de Desplazamiento

[ ]
2 −1
−1L 7 14
[K] =
EI −1 1
14 7
Entonces de acuerdo con la ecuación 3-3

[ ]{ }
2 −1 2
−q L
−1
{ D }= [ K ] {−R }=
L 7 14
EI −1 1
12 =
0
q L3 1/142
EI −1/168 { }
14 7

Las fuerzas finales en cualquier miembro se pueden determinar por superposición

de la condición restringida y del efecto de los desplazamientos de los nudos

Ai=A ri +( Aui 1 D1 + Aui 2 D 2+ …..+ A uin D n) (3.4)

Donde
Ai: Cualquier acción i, que es reacción en un apoyo, fuerza cortante, fuerza axial,
momento, desplazamiento en una sección de la estructura real
Ari : La misma acción que Ai: pero en la condición restringida sometida a las cargas
externas
Aui 1, Aui 2, …. Auin : Fuerzas en los miembros correspondientes a desplazamientos
unitarios. Específicamente Auij es la fuerza correspondiente a un desplazamiento
unitario D j=1 mientras que los valores de todos los otros desplazamientos son cero.
La ecuación de superposición para todas las barras en forma matricial es

{ A } ={ A r }+ [ A u ] { D} (3.5)
Donde
{ A } : Vector cuyos elementos son las acciones finales en las barras
{ A r } : Vector cuyos elementos son las acciones en las barras en la condición restringida
[ Au ] : son las facciones en las barras correspondientes a desplazamientos unitarios.

142
 Método de Análisis de Desplazamiento

Específicamente, los elementos de la columna j de [ Au ] son las acciones en las barras

correspondientes a un desplazamiento D j=1 mientras que todos los otros
desplazamientos son nulos.
En el ejemplo, para dibujar el diagrama de fuerzas internas se necesitan los valores
de las fuerzas cortantes y momentos flectores en los extremos de todos los miembros.
Suponemos que, en el extremo, la fuerza cortante es positiva si está dirigida hacia
arriba y el momento es positivo si actúa en el sentido horario. Por lo tanto, usando la
ecuación 3-5 tenemos
Fuerzas cortantes
{ V }= {V r } + [ V u ] {D }

Los elementos de { V r }, se determinan fácilmente de las figuras 3.2b y los elementos

de [ V u ] de las figuras 3.3b y 3.4b

{}
1
2

{} [ ] { }
V AB 1

{ }
−6−6 11/28
3 1
V BA 2 EI 6 6 q L 17/28
=qL + 2 42 =qL
V BC 1 L 0−6 EI 15/28
06 −1 /168 13/28
V CB 2
1
2

Momentos flectores
{ M }= { M r }+ [ M u ] {D }

Igualmente, Los elementos de { M r } , se determinan fácilmente de las figuras 3.2b y los

elementos de [ M u ] de las figuras 3.3b y 3.4b

143
 Método de Análisis de Desplazamiento

{}
−1

{ }
12

{} []
M AB 0
1

{ }
42 3 1 3/28
M BA EI 2 4 q L
=q L2 12 + 42 =q L2 −3/28
M BC −1 L 0 4 EI
−1/168 1
M CB 12 02
14
1
12

Los diagramas de fuerza cortante y momento flector se muestran en las figuras 3.5.
Las ordenadas del momento flector se trazan al lado de las fibras en tensión.

15qL/28
11qL/28

3.04 2.32

A 1.96 B 2.68 C

13qL/28
17qL/28

3qL2/28

qL2/14
qL2/8
qL2/8
A B C

L/2 L/2

Figura 3.5 Diagramas de fuerza cortante y momento flector

3.3 Grado de indeterminación cinemática y sistema de coordenadas

En el primer paso del método de desplazamiento se determina el número n de
desplazamientos independientes de los nudos de la estructura (grado de
indeterminación cinemática, sección 1.5) y se define un sistema de coordenadas. Una

144
 Método de Análisis de Desplazamiento

coordenada se representa con una flecha que muestra la ubicación y la dirección

positiva de un desplazamiento D o una fuerza F . El sistema de coordenadas se indica
en una figura que muestra la estructura real con n flechas en los nudos.
Los pasos restantes del método de desplazamiento implican la generación y el uso
de matrices: {R}, {D} y {k}. Los elementos de estas matrices son fuerzas o
desplazamiento en las coordenadas. Por lo tanto, es imposible seguir o comprobar los
cálculos, particularmente los signos de los parámetros, cuando el sistema de
coordenadas no está definido. Por esta razón, se recomienda que el sistema de
coordenadas se muestre en una figura que represente sólo las n flechas numeradas;
Las fuerzas aplicadas sobre la estructura no deben mostrarse en esta figura.

3.4 Generalización del método de desplazamiento

Una vez determinado el grado de indeterminación cinemática y definido el sistema
de coordenadas la respuesta de la estructura se obtiene por superposición: Sumando
el problema primario más el problema complementario.
Problema primario. Consiste en analizar la estructura con los grados de libertad
totalmente restringidos y actuando sobre ella todas las solicitaciones externas tales
como cargas, variación de temperatura, errores de montaje, asentamiento de apoyos,
pre fatiga, etc. (véase figura 3.6 b). Para obtener esta configuración se generan las

fuerzas restringentes Ri en la dirección y ubicación de cada coordenada.

Problema complementario. Consiste en analizar la estructura solicitada con

fuerzas iguales y opuestas a las fuerzas restringentes obtenidas en el problema
primario (ver figura 3.6c). Así, en el problema complementario la estructura sólo
esta solicitada por cargas aplicadas en los nudos, actuando en las coordenadas.
El problema complementario, a su vez, se puede plantear como la superposición
de los estados mostrados en la figura 3.6d, e y f
Sumando fuerzas en la dirección de cada coordenada obtenemos las ecuaciones
de equilibrio

145
 Método de Análisis de Desplazamiento

R1 +k 11 D1+ k 12 D2 + k 13 D 3=0
R2 +k 21 D1+ k 22 D2 +k 23 D 3=0
R3 +k 31 D1+ k 32 D2 +k 33 D 3=0

0 en forma matricial
{ R } + [ K ] { D }={0 }

Resolviendo esta ecuación hallamos los desplazamientos { D } de la estructura real y

usando la ecuación 3.5 se determinan las acciones finales en la estructura.

q R1 q R2
(a) D3 (b) R3
D1 D2

T° T°
= +
P P

(conocido) (conocido)

-R1 -R2 k11 k21

(c) -R3 (d) k31
D3
D1 D2

k12 k22 k13 k23

(e) + k32 (f) = 1 1 k33 D1+

Figura D2+ D3
3.6 (a)

Estructura estáticamente indeterminada. (b) Problema primario. (c) Problema complementario. (d)
D1=1 , D2=0 , D3=0 . (e) D2=1 , D1=0 , D3=0 . (f) D3=1 , D1=0 , D2=0

146
 Método de Análisis de Desplazamiento

Ejemplo 3.2. Para el pórtico plano de la figura 3.7a construir los diagramas de fuerza
cortante y momento flector. Las rigideces relativas se indican para cada miembro y se
considera que estos son inextensibles (EA=∞)

3t/m
(a)
D E F
1 8I 23 8I 4 1
5 (b)
I
4.00

8t
B C6 2
7 8I 8 9 10t.m
4.00

I
6.00

A 10
8.00 8.00

R1 3t/m
(c)

8t
R2
10t.m D1=0, D2=0

Figura (d) (e) 3.7 Pórtico plano analizado en el

ejemplo K11 K12 3.2. (a) Dimensiones y
cargas. (b)Sistema de
K21 D1+ K22 D2
D1=1, D2=0 D1=0, D2=1
coordenadas. (c),
(d) y (e) Superposición:
problema primario

más problema complementario

El grado de indeterminación cinemática del pórtico es dos: las rotaciones de los nudos
C y E . Las direcciones positivas para los componentes de desplazamiento, así como
para las fuerzas restringentes se eligen arbitrariamente como se indica en la figura 3.7b
La solución del problema se obtiene por superposición del problema primario (figura
3.7c l el problema complementario (figura 3.7c), Sumando fuerzas a lo largo de cada
coordenada obtenemos las siguientes ecuaciones de equilibrio

147
 Método de Análisis de Desplazamiento

R1 +k 11 D1+ k 12 D2=0
R2 +k 21 D1+ k 22 D2=0

En forma matricial
{R }+[ K ]{D }=0

Solución del problema primario

Las fuerzas restringentes {R } (véase figura 3.7c) son iguales a la suma de las
fuerzas de los extremos de los nudos y se calculan con ayuda del apéndice B. Se
consideran positivas cuando su dirección coincide con el de las coordenadas.
Considerando los diagramas de cuerpo libre de cada miembro y nudo de la estructura
(véase figura 3.8a), en los nudos además de las fuerzas restringentes actúan fuerzas
iguales y opuestas a la de los extremos de las barras. Las fuerzas restringentes se
calculan considerando el equilibrio de cada nudo como se muestra a continuación.
Del diagrama de cuerpo libre del nudo E

∑ M =0 : R1 +16=0 → R 1=−16

y del diagrama de cuerpo libre del nudo C

∑ M =0 : R2 +10−8=0 → R2=−2

Por lo tanto

{ }
{ R } = −16 tm
−2

Es evidente que para obtener las fuerzas restringentes es suficiente sumar las
fuerzas de los extremos en cada nudo y no es necesario obtenerlas a partir de los
diagramas de cuerpo libre. Si existe una fuerza externa aplicada en el nudo, necesita
una fuerza restringente igual y opuesta.

Solución del problema complementario

148
 Método de Análisis de Desplazamiento

Los elementos de la matriz de rigidez son las fuerzas necesarias en la ubicación y

dirección de las coordenadas para mantener la estructura en las configuraciones
deformadas de las figuras 3.7d y e. Estas fuerzas se obtienen sumando las fuerzas en
los extremos, que se toman del apéndice C. Los coeficientes de rigidez se calculan
considerando el equilibrio de cada nudo, como se muestra a continuación.

Para D1=1 , D2=0 (véase figura 3.8b)

Equilibrio del nudo E

∑ M =0 : K 11−4 EI −EI −4 EI =0 → K 11 =9 EI

Equilibrio del nudo C

EI EI
∑ M =0 :− + K 21=0 → K 21=
2 2

Para D 2=1 , D1=0 (véase figura 3.8c)

Equilibrio del nudo E

EI EI
∑ M =0 : K 12− =0→ K 12=
2 2

Equilibrio del nudo C

2 EI 17 EI
∑ M =0 : K 22−EI −4 EI − =0→ K 22=
3 3
Co estos resultados escribimos la matriz de rigidez

[ ]
1
9
2
[ K ] =EI
1 17
2 3

R1 16t.m 3t/m 16t.m

(a) D
E F
12 tn 12 tn

8t
8t.m 8t.m R2 D1=D2=0
B C 10t.m
4 tn 4 tn
149

A
 Método de Análisis de Desplazamiento

4EI k11 4EI 2EI

(b) 2EI
D E F
3 3
4EI 4EI 3 EI 34EI 3
4EI
8EI

3
1
2EI 8EI
D1=1, D2=0
B C k21

k12
(c)
D E 3
8EI
F
3 1
8EI 2EI

EI 3
8EI
2EI 4EI
B C k22
2 D1=0, D2=1
3EI
3 3 EI
4EI 4EI 6

EI EI
3 6
A

Figura 3.8. Diagramas de cuerpo libre

La inversa de la matriz de rigidez[ K ]−1 y el vector {−R } son

150
 Método de Análisis de Desplazamiento

[ ]
68 −2
−1
[K] =
1 609
EI −2 36 { }
203 , {−R }= 16
2
203 203

En esta etapa del problema se han generado todas las matrices necesarias para la
solución del problema. Lo que sigue requiere solo aplicar el álgebra matricial

Cálculo del vector desplazamiento {D }

De acuerdo con la ecuación 3.3
−1
{ D }= [ K ] {−R }

Sustituyendo los valores de[ K ]−1 y {−R } en la ecuación anterior se obtiene

{ }
{ D }= 1 1.766831 radianes
EI 0.197044

Las fuerzas finales en los extremos de los miembros se obtienen por superposición
usando la ecuación 3.5.
Las fuerzas cortantes finales en los extremos de los miembros se obtienen
sumando las fuerzas cortantes en las mismas ubicaciones del problema primario y
complementario:
{ V }= {Vr } + [ V u ] { D } (a).
Los elementos de las matrices { Vr } que corresponden al problema primario se
obtienen fácilmente con ayuda de la figura 3.8a. Así para los extremos 1,2,…,10 (véase
figura 3.7a) y asumiendo que fuerzas cortantes son positivas cuando están dirigidas
hacia arriba y hacia la derecha se tiene

{ Vr }={ 0 0 12 12 0 0 4 4 0 0 } tn

Para escribir la matriz de las fuerzas cortantes de extremo debidos a los

desplazamientos unitarios[ V u ] , ponemos los valores en los extremos 1, 2,…,10 en la
primera y segunda columnas respectivamente las fuerzas cortantes que se muestran
en las figuras 3.8c y d.

151
 []
Método de Análisis de Desplazamiento

−3
0
4
+3
0
4
−3
0
4
+3
0
4
−3 −3
[V u ]=EI 8 8
+3 +3
8 8
−3
0
4
+3
0
4
−1
0
6
+1
0
6

Reemplazando { Vr }, [V u ] y { D } en la ecuación (a) obtenemos:

152
 []
Método de Análisis de Desplazamiento

−3
0
4
+3
0
4

{} {}
−3
0
0 4 −1.33
0 +3 1.33
0
12 4 10.67
12 −3 −3 13.33
{ V }= 0
0
+ EI 8 8
{
1 1.766831
+ 3 +3 EI 0.197044
= −0.74
0.74
ton
}
4 8 8 3.85
4 −3 4.15
0 0 −0.03
4
0 +3 0.03
0
4
−1
0
6
+1
0
6

Así obtenemos similarmente los momentos finales en los extremos de las barras se
obtienen sumando los momentos del problema primario y complementario

{ M }= { Mr }+ [ M u ] { D } (b)

Un momento de extremo es positivo si actúa en sentido horario. Los valores de los

momentos en los extremos 1, 2, …, 10, para la condición restringida (figura 3.8a) son

{ Mr }={ 0 0 −16 16 0 0 −8 8 0 0 } tn . m

La matriz de los momentos de extremo debidos a los desplazamientos unitarios se

genera con los momentos en los extremos 1, 2,. . .,10, de las figuras 3.8c y d:

153
 Método de Análisis de Desplazamiento

[]
2 0
4 0
4 0
2 0
1
1
2
[ Mu ]=EI 1 1
2
0 2
0 4
2
0
3
1
0
3

Remplazando { Mr }, [ Mu ] y { D } en la ecuación (b) obtenemos

Los diagramas de fuerza cortante y momento flector se muestran en la figura 3.9

154
Método de Análisis de Desplazamiento

155
 Método de Análisis de Desplazamiento

Figura 3.9 Diagramas de fuerzas cortantes y momentos flectores

3.5 Análisis para diferentes estados de carga.

El método de rigidez se usa eficazmente cuando se requiere analizar una estructura
bajo diferentes estados de carga, pues siendo la matriz de rigidez (y su inversa)
independiente de la carga aplicada solo será necesario resolver el problema primario
para cada estado de carga y emplear la ecuación matricial 3-2 Si el número de estados
de carga es p, la solución se puede combinar en la siguiente ecuación matricial:

[ K ] n ×n [ D ]n × p=[−R ]n∗ p

De la cual despejando [D], se tiene

−1
[ D ] n × p= [ K ]n × n [−R ]n∗ p (3.6)

En la que cada columna de [D] y [-R] corresponde a un estado de carga

3.6 Análisis para efectos ambientales.

El método de rigidez se puede usar para analizar los efectos separados o
combinados de cambio de temperatura, errores de montaje, asentamiento de apoyos,
pre esfuerzo, etc. En estos casos sigue siendo válida la ecuación 3-3, y las fuerzas
restringentes del problema primario son las fuerzas necesarias para impedir los
desplazamientos de los nudos debidos a los efectos dados.
Cuando se efectúa el análisis para el caso del movimiento de un apoyo, también se
puede aplicar la ecuación 3-3, siempre que el desplazamiento no se presente en la
dirección de una de las coordenadas. Cuando el desplazamiento ocurre en la dirección
de una coordenada, se necesita modificar la ecuación 3-3 como se explicará en el
ejemplo 3.

3.7 Cinco pasos del método de desplazamiento

El análisis por el método de desplazamiento consta de cinco pasos que se resumen
como sigue:

156
 Método de Análisis de Desplazamiento

Paso 1 Definir el sistema de coordenadas que representan los desplazamientos

desconocidos de los nudos. Definir también las acciones requeridas [ A ] m × p , así como su
convención de signos (si es necesario).
Paso 2 Resolver el problema primario, es decir determinar [ R ] n × p, y [ Ar ]m × p.
Paso 3 Introducir uno por uno los desplazamientos unitarios en las coordenadas y
generar [ K ] n ×n y [ Au ]m × n.
Paso 4 Resolver las ecuaciones de equilibrio:

[ K ] n ×n [ D ]n × p=−[ R ] n × p (3.7)

Paso 5 Calcular las acciones requeridas por superposición;

[ A ] m × p =[ A r ]m × p+ [ A u ]m ×n [ D ] n× p
(3.8)

De manera semejante al método de fuerza (sección 2.6), al finalizar el paso 3, se

han generado todas las matrices necesarias para el análisis. Los pasos 4 y 5 solo
requieren simple algebra matricial.
Para una referencia rápida, los símbolos usados en esta sección se definen
nuevamente como sigue:
n , p , m=¿ número de desplazamientos desconocidos, número de casos de
carga, y número de acciones requeridas.
[ A ] =¿ acciones requeridas (respuesta del problema)
[ A s ]=¿ valores de las acciones en la estructura restringida debido a las
solicitaciones externas.
[ Au ] =¿ Valores de las acciones debido a desplazamientos unitarios
actuando por separado en cada coordenada.
[ R ] =¿ Fuerzas en las coordenadas necesarias para impedir los
desplazamientos de los nudos cuando actúan las solicitaciones
externas.
[ k ]=¿ matriz de rigidez.

157
 Método de Análisis de Desplazamiento

Ejemplo 3.3 Obtener el diagrama de momento flector para el pórtico plano de la figura
3.10 debido al efecto por separado de las siguientes solicitaciones: (1) las cargas
mostradas; (2) Un asentamiento δ D = 10 mm. Considere que la rigidez a la flexión EI =
constante y no considere el cambio de longitud de los miembros
80kN
1 2
(a) (b)
B 2.50 C 3

2 .0 0
20kN

2 .0 0
A 3.00 5.00 D

80kN R2 R2
R1 R1
(c) R3 (d) R3
B C B C 10mm

20kN

D1=0, D2=1, D3=0 D1=0, D2=1, D3=0

A D A D mm

(e) (f)
k11 k21 k22
k31 k12 k32
1
1 B 1 C B C
1

D1+ D2
D1=1 D2=1

A D A D

(g)
k23
1 1 k33
B C
k13

Figura 3.10 (a) + D3 Pórtico

analizado en el D3=1 ejemplo 5.3. (a)
Dimensiones y carga del
pórtico. (b) A D
Sistema de
coordenadas. (c) y (d) Problema

158
 Método de Análisis de Desplazamiento

primario caso 1 y 2 respectivamente: D 1=D 2 =D 3=0. (e), (f) y (g): Problema complementario cuando
(d) D 1=1 , D 2=1 y D 3 =1, respectivasmente
Paso 1 La estructura tiene tres grados de libertad: dos rotaciones y una traslación
como se muestra en la figura 3.10b que también muestra el sistema de coordenadas.
Para dibujar los diagramas de momentos flectores se requieren los valores momentos
en los extremos de todos los miembros; estos se consideran positivos si actúan en el
sentido horario.

[{ } { } ]
M AB M AB
M BA M BA
M BC M BC
{ A}= (a)
M CB M CB
M CD M CD
M DC caso 1 M DC caso 2

Paso 2 Las fuerzas en los extremos en la estructura restringida para los casos 1 y 2 se
han calculado con ayuda del apéndice B y C se muestran en las figuras 3.11a y b. Las
fuerzas de restricción R1, R2 y R3 necesarias para anular los desplazamientos de los
nudos (D1=D2=D3=0). (Véase figuras 3.10c y d) se determinan por considerando
diagramas de cuerpo libre de todos los miembros y nudos de la estructura restringida.
Caso 1. De los diagramas de cuerpo libre mostrados en la figura 3.11a
Equilibrio del nudo B

∑ Fy=0 :−40−10 cosθ+ N AB senθ=−40−10 ( 0.6 ) + N AB ( 0.8 )=0 → N AB=57.5 kN

∑ Fx=0 : N AB cosθ+10 senθ−N BC =N AB ( 0.6 )+10 ( 0.8 )−N BC =0→ N BC =42.5 kN
∑ M =0 : R1 +50−8.333=0→ R1=−41.667 kN . m

Equilibrio nudo C

∑ Fx=0 : N BC +10+ R3 =0 → R3 =−52.5 t

∑ M =0 : R2 −50−10=0 → R2 =60 t−m

159
 Método de Análisis de Desplazamiento

Caso 2. En la estructura restringida el asentamiento δ D (figura 3.10d) solo induce

momentos en el miembro BC que se calculan con ayuda de apéndice C.

−6 EI δ D −6 ×100000 ×0.04
M B C =M CB = 2
= 2
=−240 kN . m
l 5

NBC 80kN NBC NBC

R1 50kN.m 50kN.m R2 R1 240kN.m 240kN.m NBC R2
(a) (b)
B C R3 B C R3
40kN 40kN 96kN 96kN
NAB 10kN.m NAB
10kN

20kN

Sen

Cos 

A 10kN 10kN.m A
D D
NAB

NBC
K11 0.8EI 0.4EI K 21 K12 0.4EI 0.8EI NBC K 22
(c) (d)
B C K 31 B C K32
NBC 6 6
NBC
6 6
25EI 25EI 25EI 25EI
0.8EI NAB NAB
3 EI
6
25 EI 8EI

0.4EI A 3
6 A 8EI 0.5EI
25EI D D

NBC 240kN.m 240kN.m NBC

(b) R1 R2 K11 0.8EI 0.4EI NBC K21
B C R3 (c)
B C K31
96kN 96kN NBC 6 6
NAB 25EI 25EI
0.8EI NAB
6
25 EI

A 0.4EI
D A 6
25 EI
D

K12 0.4EI 0.8EI NBC K22

(d) K 32
B C
NBC 6 6
25EI 25EI
NAB EI
3
8EI

A 3
8EI 0.5EI
D

K13 0.18EI 0.18EI NBC K 23

(e) K 33
B C
NBC 9 9
0.3EI 125EI 125EI
NAB 3
3 0.375EI 16EI
25EI

0.3EI A 0.375EI 3
D 16EI
3
25EI

160
 Método de Análisis de Desplazamiento

Figura 3.11 Diagramas de cuerpo libre correspondientes a ; (a) Problema primario, caso 1. (b) Pronlema
primario caso 2. (c) D 1=1 , D 2=D 3=0. (d) D 2=1 , D 1=D 3=0. (e) D 3=1 , D 1=D 2=0.

De los diagramas de cuerpo libre que se muestran en la figura 3.11b, se tiene

Equilibrio del nudo B
∑ Fy=0 : N AB senθ−96=N AB ( 0.8 )−96=0→ N AB=120 kN
∑ Fx=0 : N AB cosθ−N BC =120 ( 0.6 )−N BC =0 → N BC =72 kN
∑ M =0 : R1 +240=0 → R1 =−240 kN . m

Equilibrio nudo C
∑ Fx=0 : N BC + R3 =0 → R3 =−72 kN
∑ M =0 : R2 + 48=0 → R2=−240 kN . m

Por lo tanto, la matriz [ R ] queda definida como

[ ]
−41.667 −240
[ R ]= 60 −240 (b)
−52.5 −72

Los momentos de extremo para la estructura restringida se generan con los valores
de los momentos en los extremos de las figuras 3.11a y b; así obtenemos

[ ]
−8.333 0
8.333 0
[ M r ]= −50
50
−240 kN . m
−240
(c)
10 0
−10 0

161
 Método de Análisis de Desplazamiento

Paso 3 Las fuerzas necesarias para inducir por separado desplazamientos unitarios en
las coordenadas (véase figura 3.10e, f y g) se obtienen considerando el equilibrio de los
nudos en los diagramas de cuerpo libre de la figura 3.11c, d y e

Para D 1=1 , D2=D3=0(véase figura 3.11c)

Equilibrio del nudo B

6 6 6 6 3 EI
∑ Fy=0 : EI − EIcosθ−N AB senθ= EI − EI (0.6)−N AB (0.8)=0→ N AB=
25 25 25 25 25
6 −3 EI 6 3 EI
∑ Fx=0 :−N AB cosθ+ EIsenθ−N BC = ( 0.6 )+ EI ( 0.8 )−N BC =0 → N BC =
25 25 25 25
4 EI 4 EI 8 EI
∑ M =0 : K 11− − =0 → K 11= =1.6 EI
5 5 5

Equilibrio nudo C
−3 EI
∑ Fx=0 : N BC + K 31=0 → K 31 = =−0.12 EI
25
2 EI 2 EI
∑ M =0 : K 21− =0 → K 21 = =0.4 EI
5 5

Para D 2=1 , D1=D3=0(véase figura 3.11d)

Equilibrio del nudo B

6 EI 6 EI 3 EI
∑ Fy=0 : −N AB senθ= −N AB (0.8)=0 → N AB =
25 25 10
−3 EI 9 EI
∑ Fx=0 :−N AB cosθ+ N BC = ( 0.6 )+ N BC =0 → N BC =
10 50
2 EI 2 EI
∑ M =0 : K 12− =0 → K 12= =0.4 EI
5 5

Equilibrio nudo C
9.375 EI −9 EI 3 EI −39 EI
∑ Fx=0 :−N BC + + K 32= + + K 32 =0 → K 32= =−0.195 EI
L
2
50 8 200
4 EI 9 EI
∑ M =0 : K 22− −EI =0 → K 22= =1.8 EI
5 5

162
 Método de Análisis de Desplazamiento

Para D 2=1 , D1=D3=0(véase figura 3.11e)

Equilibrio del nudo B

9 EI 3 EI 9 EI 3 EI 9 EI
∑ Fy=0 : −N AB senθ + cosθ= −N AB (0.8)+ (0.6)=0 → N AB=
125 25 125 25 50
3 EI −9 EI 3EI 51 EI
∑ Fx=0 :−N AB cosθ− senθ + N BC = ( 0.6 ) − ( 0.8 ) + N BC =0 → N BC =
25 50 25 250
3 EI 9 EI −3 EI
∑ M =0 : K 13 + − =0 → K 13= =−0.12 EI
10 50 25

Equilibrio nudo C
3 −51 EI 3 783 EI
∑ Fx=0 :−N BC − EI + K 33= − EI=0→ K 33= =0.3915 EI
16 250 16 2000
9 EI 3 EI −39 EI
∑ M =0 : K 32− + =0 → K 32= =−0.195 EI
50 8 200

Por lo tanto, la matriz de rigidez y su inversa resultan.

[ ] [ ]
1.6 0.4 −0.12 −1 1 0.668788 −0.133622 0.138437
[ K ] =EI 0.4 1.8 −0.195 [ K ] = −0.133622 0.613940 0.264837
EI
−0.12 −0.195 0.3915 0.138437 0.264837 2.728622

Los elementos de la matriz [ M u] se determinan fácilmente con los valores de los

momentos de extremo de la figura 3.11c, d y e. Por lo tanto

[ ]
0.4 0 −0.3
0.8 0 −0.3
[ Mu ]=EI 0.8 0.4 0.18
0.4 0.8 0.18
0 1 −0.375
0 0.5 −0.375

Paso 4 Sustituyendo [ K ]−1 y [ R ] en la ecuación 3.6 y resolviendo se obtiene [ D ]

−1
[ D ] =[ K ] [ R ]=
EI[−0.133622 0.613940 0.264837
0.138437
60
0.264837 2.728622 −52.5 ][
1 0.668788 −0.133622 0.138437 −41.667 −240
−240
−72 ]
163
 Método de Análisis de Desplazamiento

[ D ]=
EI [
1 43.151484 138.407367
−28.50006 134.344529
133.130693 293.246659 ]
Paso 5 Los momentos finales en los extremos de las barras se obtienen Sustituyendo
[ M r ] , [ Mu ] y [ D ] en la ecuación [ M ] =[ M r ] + [ M u ] [ D ] . Así se tiene

[ ] [ ][
−8.333 0 0.4 0 −0.3

]
8.333 0 0.8 0 −0.3
43.151484 138.407367
[ M ]= −50 −240 + EI 0.8 0.4 0.18 1 −28.50006 134.344529
50 −240 0.4 0.8 0.18 EI
133.130693 293.246659
10 0 0 1 −0.375
−10 0 0 0.5 −0.375

[ ]
−31.01 −32.61
2.92 22.75
[ M ] = −2.92 −22.75 kN . m
68.42 −24.38
−68.42 24.38
−74.17 −42.80

Donde las columnas de la matriz [ M ] son los momentos finales en los extremos de
los miembros para los casos 1 y 2 respectivamente. Los diagramas de momento flector
se muestran en la figura 3.12.

(a) Caso 1

68.42
P L /4 = 1 0 0

2.92 68.42

64.35 PL/4=15

31.01
74.17

DMF en kN.m

164
 Método de Análisis de Desplazamiento

(b) caso 2

22.75

24.38

24.38

32.61

DMF en kN.m 42.80

Figura 3.12 Diagrama de los momentos flectores

En los pasos 2 y 3 se calcularon la fuerza de restricción R3 y los coeficientes de

rigidez k 31 , k 32 y k 33 por consideraciones de equilibrio. En lo que sigue se muestra el
cálculo de estas fuerzas usando el principio de los desplazamientos virtuales.
M AB=−8.333 kN .m
M BA =8.333 k N .m
M BC =−50 kN .m

165
 Método de Análisis de Desplazamiento

M CB =50 kN . m
M CD =10 kN . m
M BC =−10 kN .m

(C)

M BC =M BC =−240 kN . m
(d)

M AB=0.4 EI ; M BA=0.8 EI

M BC =0.8 EI ;M CB =0.4 EI

(e)

M BC =0.4 EI ; M CB=0.8 EI
M CD =EI ; M CD =0.5 EI

(f)

M AB=0.3 EI ; M BA =0.3 EI
M BC =M CB =0.18 EI

M CD =M DC =−0.375 EI

Figura 3.13, Calculo de la fuerza restringente R 3, y los coeficientes de rigidez k 31 , k 32 y k 33 , usando trabajo
virtual. (a) Fuerzas actuando sobre el mecanismo para el caso 1. (b) Desplazamiento virtual. (c) fuerzas

166
 Método de Análisis de Desplazamiento

externas sobre el mecanismo para el caso 2. (d), (e) y (f) Fuerzas sobre el mecanismo cuando
D1=1 , D2=1 y D 3=1 respectivamente
Para hallar la fuerza R3 se aplica la ecuación del trabajo al mecanismo mostrado en la
figura 3.15a en el que se han introducido articulaciones en los nudos. Los miembros de
este mecanismo están sometidos a las mismas fuerzas externas la estructura primaria
(véase figura 3.10c), en donde los momentos de extremo están representados como
cargas externas. Las fuerzas en la figura 3.15a que incluye la fuerza R3 y las
reacciones en A y D (no mostradas) constituye un sistema en equilibrio. Ahora se
impone un desplazamiento unitario virtual en la coordenada 3 que origina que los
miembros roten y se trasladen como cuerpos rígidos (sin deformación) como se
muestra en la figura 3.15b. Si el mecanismo está en equilibrio el trabajo que realizan
las fuerzas (cargas o pares) en la figura 5.15a cuando se mueven a través de los
desplazamientos virtuales en la figura 4.15b debe ser igual a cero; por lo tanto,
podemos escribir
θ AB ( M AB + M BA ) +θ BC ( M BC + M CB ) +θCD ( M CD + M DC ) + R3 × 1+ ( 20× 0.625+80 × 0.375+ 20× 0.5 )=0

Reemplazando los valores de los momentos en los extremos de los miembros en la

ecuación anterior se tiene
0.25 (−8.333+8.333 )−0.15 (−50+50 ) +0.25 ( 10−10 )+ R3 × 1+ ( 20× 0.625+80 × 0.375+ 20× 0.5 )=0
De donde obtenemos
R3=−52.5 kn

De manera similar multiplicando las fuerzas de la figura 3.15c por los desplazamientos
virtuales (figura 3.15b) obtenemos la fuerza de restricción R3 para el caso 2
0.25 ( 0+0 )−0.15 (−240−240 ) +0.25 ( 0+0 ) + R3 ×1=0

De donde obtenemos
R3=−72kN
Para obtener los coeficientes de rigidez k 31 , k 32 y k 33 , multiplicamos las fuerzas de las
figuras 3.15d, e y f que representan las configuraciones deformadas para D1=1 , D2=1 y
D3=1 con los desplazamientos virtuales de la figura 3.15b. Así se tiene

167
 Método de Análisis de Desplazamiento

0.25 ( 0.4 EI + 0.8 EI )−0.15 ( 0.8 EI +0.4 EI )+ k 31 ×1=0

k 31=−0.12 EI

−0.15 ( 0.4 EI + 0.8 EI ) +0.25 ( EI + 0.5 EI ) × k 32 × 1=0

k 32=−0.195 EI

0.25 (−0.3 EI −0.3 EI ) −0.15 ( 0.18 EI +0.18 EI ) +0.25 (−0.375 EI −0.375 EI ) +¿

k 33 ×1=0
k 33=0.3915 EI
Que son los mismos resultados obtenidos considerando el equilibrio de los nudos.

168
 Método de Análisis de Desplazamiento

Ejemplo 3.4 Pórtico plano solicitado por cargas y movimiento de los apoyos
Dibujar los diagramas de momento flector para el pórtico plano mostrado en la figura
3.14a, sometido a las siguientes solicitaciones: (a) La carga indicada en la figura 3.37a;
(b) el apoyo A se asienta 20 mm, se mueve hacia la derecha 10 mm y gira 0.01
radianes en el sentido horario. Tómese EI =¿

(a) C (b) C

P=8kN
B 1 2
B
D D
2kN/m

A A

Figura 3.14. (a) Estructura estáticamente indeterminada. (b) Sistema de coordenadas

Paso 1. El pórtico tiene dos grados de libertad, que se indican con las coordenadas 1 y
2 (véase figura 3.14b). Para dibujar el diagrama de momento flector se requieren los
valores de los momentos en los extremos de todos los miembros.

[{ } { } ]
M AB M AB
[ M ] = M BA M BA
M BC M BC
M BD a M BD b

Además M CB =M DB =0. Se supone que un momento de extremo es positivo cuando

actúa en el sentido horario.

Paso 2 Las fuerzas en los extremos de los miembros para la condición restringida
cuando la estructura está sometida a las cargas y a los movimientos de apoyo (véase
la figura 3.15a y b) se encuentran con ayuda de los apéndices B y C y las fuerzas
restringentes se calculan como la suma de las fuerzas en los extremos en cada nudo

169
 Método de Análisis de Desplazamiento

como se indica en las figuras 3.16a y b

(a) C (b) C

P=8kN
B R1 R2 B R1 R2

D D
m r
apim
aio
r

m
le
b

N/
ro

k
P

A Caso a) Cargas A Caso b) Movimiento

de apoyo

+
(c) C

B K11 K21

D D1

D1=1

+
(d)
io

C
m o
acm m
le
p n
e r
ta

B K12 K22
r
Pole
b

D
D2

D2=1

170
 Método de Análisis de Desplazamiento

Figura 3.15 Superposición: problema primario y complementario. (a), (b) Estructuras restringidas para el
caso de cargas y movimiento de apoyo respectivamente. (c) y (d) Configuraciones deformadas cuando
D1=1 y D 2=1 respectivamente
Caso a: cargas (véase figura 3.15a)

R1=M BA + M BC + M BD=4.167+0−15.36=−11.193 kN
R2=−V BA =−5 kN
Caso b) Desplazamientos de apoyos (véase figura 3.39b)

R1=M BA + M BC + M BD=0.0064 EI +0+ 0.0012 EI =+0.0076 EI

R2=−V BA =−0.00336 EI

q l2
M AB= =M BA =4.167 kN . m
12
M BC =0
(a) C

M
B BD
P=8kN
M BD =
Pab
l 2
b
2( )
a+ =15.36 kN . m

MBA
VBA ql
D V BA= =5 kN
2
N
k
2m
/

A VAB 6 EI ∆ H 4 EI θ A
M AB= + =0.0104 EI
MAB
l 2
l
(b) C
6 EI ∆ H 2 EI θ A
M BA = + =0.0064 EI
l2 l
B MBD
MBA
VBA
M BC =0
D

M BD =3 ¿ ¿
12 EI ∆ H 6 EI θ A
V BA= + =0.00336 EI
VAB A l3 l2
MAB

171
 Método de Análisis de Desplazamiento

Figura 3.16 Fuerzas en los extremos de las barras para la estructura restringida. (a) Para el caso de
carga. (b) Caso de movimiento de apoyo
Así la matriz [ R ] resulta ser

[
[ R ] = −11.193
−5
0.0076 EI
−0.00336 EI ]
Cada columna de [R] son las fuerzas restringentes para los casos (a) y (b)
respectivamente. Los momentos de extremo cuando los desplazamientos en las
coordenadas están impedidos (véase figuras 5.16a y b) son

[ ]
−4.167 0.0104 EI
[ M r ]= +4.167
0
0.0064 EI
0
−15.36 0.0012 EI

Cada columna de [ M r] son los momentos en los extremos de las barras

correspondientes al problema primario para los casos (a) y (b) respectivamente

Paso 3. Las fuerzas de extremo cuando se aplican los desplazamientos unitarios por
separado (véase las figuras 3.15c y d) se hallan con ayuda del Apéndice C , y los
coeficientes de rigidez se calculan como la suma de las fuerzas de extremo en cada
nudo como se indica en las figuras 3.17a y b, así se obtiene

Para D 1=1 , D2=0 (véase figura 3.17a)

K 11 =M BA + M BC + M BD=0.8 EI + EI +0.6 EI =2.4 EI

K 21 =−V BA +V BC =−0.24 EI + 0.3333 EI =0.0933 EI

Para D 2=1 , D1=0 (véase figura 3.17 b)

K 12 =−M BA + M BC =−0.24 EI +0.3333 EI =0.0933 EI

172
 Método de Análisis de Desplazamiento

K 22=V BA +V BC =0.096 EI +0.1111 EI =0.2071 EI

Por lo tanto, la matriz de rigidez y su inversa son

[ K ] =EI [ 0.0933
2.4 0.0933
0.2071 ] [
[ K ]−1= 1 0.4241 0.0933
EI 0.0933 4.9145 ]
(a) C
VBC
2 EI
M AB= =0.4 EI
MBC l
B VCB
4 EI
MBA M BA = =0.8 EI
D l
VBA MBD
3 EI
D1=1 M BC = =EI
l
3 (2 EI )
M BD = =0.6 EI
l
A VAB 6 EI
V BA= 2
=0.24 EI
MAB l

(b)
C VBC

6 EI
M AB=M BA =0.24 EI
MBC VBC l
2

B VBA D
MBA 3 EI
M BC = 2
=0.333 EI
l
D2=1 M BD =0
12 EI
V BA= =0.096 EI
l3
A VAB
MAB

Figura 3.17 Fuerzas en los extremos de las barras para la estructura cuando. (a) D 1=1. (b) D 2=1

También la matriz [Mu] resulta ser

173
 Método de Análisis de Desplazamiento

[ ]
0.4 −0.24
[ M u ]=EI 0.8
1
−0.24
0.333
0.6 0 .00

Paso 4 Sustituyendo [ K ] , y [ R ] en la ecuación 3-6, se encuentra

−1

−1
[ D ] =[ K ] [−R ] =
EI 0.0933 4.9145 5[
1 0.4241 0.0933 11.193 −0.0076 EI
0.00336 EI ][ ]

[ ]
3.7915
−0.003865
EI
[ D]=
22.4330
0.017965
EI

Cada columna de [D] representa los desplazamientos en los casos (a) y (b).

(a) C
13.09
Pab/L=19.2

B
1.82
11.27
D
11.35

qL2/8=6.25 1.33

DMF en kN.m
8.03 A

(b) C

15.11
B
13.56 28.66
D

DMF en kN.m
A 61.32

174
 Método de Análisis de Desplazamiento

Figura 3.18 Diagramas de los momentos flectores (a) caso de cargas. (b) Caso de movimiento de apoyos

Paso 5 Los momentos finales de extremos se calculan sustituyendo [ M r ] , [ M u ] , y [ D ]

en la ecuación 3.8 y después de sustituir el valor de EI , se obtiene

[ ] [ ][ ]
−4.167 0.0104 EI 0.4 −0.24 3.7915 −0.003865
[ M ]= 4.167 0.0064 EI + EI 0.8 −0.24 EI
0 0 1 0.333 22.4330
0.017965
−15.36 0.0012 EI 0.6 0 .00 EI

[ ]
−8.03 61.32
[ M ] = 1.82 −13.56 kN . m
11.27 28.66
−13.09 −15.11

Las columnas de esta matriz representan los momentos finales de extremo para los
casos (a) y (b. En la figura 3.18 se muestran los diagramas de momento flector

Ejemplo 3. 5. Para la armadura plana mostrada en la figura 3.19a hallar las fuerzas
axiales en todas las barras. Todos Los miembros tienen el mismo EA

1m 1m

C D C D
2
60°

B B

30°

A A
P 1

Figura 3.19 (a) Armadura plana analizada en el ejemplo 3.5. (b) Sistema coordenado

Paso 1 El grado de indeterminación cinemática es cuatro, pues para una carga

arbitraria se producirá dos traslaciones en A y dos en B. Debido a la simetría

175
 Método de Análisis de Desplazamiento

geométrica y de carga, los nudos A y B solo se desplazan a lo largo de la vertical,

reduciéndose el grado de indeterminación a tan solo dos como se muestra en la figura
3.19b, que también indica el sistema de coordenadas. Las acciones requeridas son las
fuerzas axiales en los miembros

{ A r }= {N 1 N 2 N 3 N 4 N 5 }
Las fuerzas en los miembros se consideran positivas si son de tensión.

Figura 3.20 (a) Problema primario: D1=0, D2=0. (b) configuración deformada para D1=1, D2=0. (c)
Configuración deformada para D1=0, D2=1. (d) y (e) Diagramas de cuerpo libre para D1=1, y D2=1
respectivamente. (f) fuerzas finales en los miembros

176
 Método de Análisis de Desplazamiento

Paso 2. Puesto que en las armaduras las cargas se aplican solo en los nudos, las
fuerzas restringentes en las coordenadas 1 y 2 son fuerzas iguales y opuestas a las
cargas externas. Para el ejemplo se debe aplicar una fuerza igual y opuesta a la carga
P en la coordenada 1 (véase figura 3.20a; por lo tanto, se tiene

{0}
{ R } = −P

También queda claro que las fuerzas son cero en todos los miembros, es decir

{ A r }= { 0 }

Paso 3 Las fuerzas necesarias para mantener la armadura en la configuración

deformada que se ilustra en las figuras 3.20b y c) se obtienen considerando el equilibrio
de los nudos en cada caso, como se muestra a continuación

Para D1=1 , D2=0 . En general las deformaciones y las fuerzas de las barras se
determinan como:

(EA ¿¿ i)
δ i=⃗ui ⃗d j N i= ¿
l

En las expresiones anteriores δ i es la deformación del miembro i ; u⃗ i es el vector

unitario en la dirección del miembro ij ; d⃗ jes el desplazamiento del nudo j (extremo de la

barra que se desplaza), N i es la fuerza axial en la barra i , EA i es la rigidez axial del
miembro i , y l i es la longitud del miembro i . De la figura 3.20b se deduce

δ 1=δ 2=0 → N 1=N 2=0

δ 3=δ 5= ( 12 ,− √23 ) .( 0 ,−1)= √23 → N =N = EA2 ( √23 )= √43 EA

3 5

( 1 )= √ EA
EA 3
δ 4 =( 0 ,−1 ) . ( 0 ,−1 )=1 → N 4=
2
√3 2
3

177
 Método de Análisis de Desplazamiento

Considerando el equilibrio de los nudo A y B (véase figura 3.20d)

Nudo A

∑ Fy=0 : N 3 cos 30 º+ N 4 + N 5 cos 30 º−K 11=0

( √43 EA )( √23 )+ √23 EA+( √43 EA )( √23 )−K 11

=0

K =(
4 )
3+ 2 √ 3
11 EA

Nudo B

∑ Fy=0 : N 1 cos 60+ N 2 cos 60−N 4−K 21=0+0−

√ 3 EA −K =0
21
2
−√ 3
K 21 = EA
2

Cuando D 2=1 , D1=0 (véase figura 3.15c)

δ 1=δ 2= ( √23 ,− 12 ) .( 0,1) = 12 → N =N = 2EA√3 ( 12 )= √43 EA

1 2

3
δ 3=δ 5=0 → N 3=N 5=0

(−1 ) = √ EA
EA − 3
δ 4 =( 0 ,−1 ) . ( 0 ,1 ) =−1→ N 4 =
2
√3 2
3

Considerando el equilibrio de los nudos A y B (véase figura 3.20e)

Nudo A

∑ Fy=0 : N 3 cos 30 º+ N 4 + N 5 cos 30 º−K 12=0+ ( −2√ 3 EA )+0−K 12

=0

178
 Método de Análisis de Desplazamiento

−√ 3
K 12 = EA
2

Nudo B

∑ Fy=0 : N 1 cos 60+ N 2 cos 60−N 4−K 22= ( √43 EA )( 12 )+( √43 EA )( 12 )−( −2√ 3 EA )−K 22 =0

3 √3
K 22= EA
4

Matriz de rigidez y su inversa

[ ]
3+2 √ 3 − √3
[ K ] =EA 4
− √3
2
3√3
−1
y [K ] = [
1 0.9628 0.6418
EA 0.6418 1.1977 ]
2 4

Las fuerzas axiales en los miembros producidos por valores unitarios separados de los
desplazamientos D1 y D 2, son

[]
0
√3
4
0
√3
4
[ Au ] =EA √3 0
4
√3− √3
2 2
√3 0
4

Paso 4 Los desplazamientos de los nudos {D} se determinan sustituyendo [ K ] y { R } en

la ecuación 3.7. Despejando {D} de obtiene

179
 Método de Análisis de Desplazamiento

−1
{ D }= [ K ] {−R }= [ ]{ } {
1 0.9628 0.6418 P
=
P 0.962765
EA 0.6418 1.1977 0 EA 0.641844 }
5. Las fuerzas finales en los miembros se calculan por superposición con la ecuación
3.8

{ N }={N r }+[ N u ]{D }

[]
0√
3
4

{} { }{ }
0
√3
0 4 0.27793 2223.41
0 0.27793 2223.41
{ N }= 0
0
+ EA √ 3 0
4 {
P 0.962765
EA 0.641844 }
=P 0.41689 = 3335.16 Kgf
0.27793 2223.41
0 √3 − √ 3 0.41689 3335.16
2 2
√3 0
4

Los valores finales de las fuerzas en las barras se muestran en la figura 3.20f

Ejemplo 3.6 Análisis de un emparrillado horizontal

Encontrar las reacciones en el apoyo A del emparrillado horizontal (figura 3.21a) debido
a una carga vertical uniforme de intensidad q=20 kN / m sobre el tramo AEC . Las
rigideces a la flexión EI y a la torsión GJ son constantes. Tómese E=200 GPa; y
considerar que GJ / EI =0.5

Paso 1 El grado de indeterminación cinemática del emparrillado es tres, pues hay tres
desplazamientos independientes posibles: dos rotaciones alrededor de los ejes x y z y
una traslación en la dirección del eje y del nudo E , los que están representados por las
tres coordenadas en la figura 3.21b. Las acciones requeridas y sus direcciones
positivas también se definen en la misma figura y son:

180
 Método de Análisis de Desplazamiento

{ A } ={ A 1 A 2 A 31 }

(a)

(b) A3
A2

A1 B

E
1 2

181
 Método de Análisis de Desplazamiento

Figura 3.21 (a) Emparrillado analizado en el ejemplo 3.6. (b) Sistema de coordenadas y direcciones
positivas de las componentes de reacción en A y B

Paso 2. Las fuerzas en las tres coordenadas para la condición restringida que se
muestran en la figura 3.22a se hallan sumando las fuerzas en los extremos del nudo E
(véase la figura 3.23), y que se han calculado con ayuda del apéndice B; de este modo
se tiene

−qL qL −ql −20 ×12

R1=−V EA −V EC = − = = =−120 kN
3 6 2 2
R2=0
2 2 2 2
ql q l q L 20× 12
R3=M EA−M EC = − = = =80 kN . m
27 108 36 36

(a) (b)

A A
B B
k11
R3 k31
E E
D1
R1 R2 k21

C C

D D
D1=1
D1=D2=D3=0

(c) (d)

A A
B B
k12 k13

k32 k33
E E
D2 + D3
k22 k23

C C

D D
D2=1 D3=1

182
 Método de Análisis de Desplazamiento

Fig. 3.22 Superposición. (a)Problema primario, D1=D2=D3=0. (b) D1=1, D2=D3=0. (c) D2=1, D1=D3=0. (d) D3
=1, D1 =D2=0

A
q l AE ql
MAE
V AE=V EA = =
VAE B
2 3
MEA

MEA
E
MCE q l 2AE q l 2
VED VEC M AE=M EA= =
C
VCE 12 27
D

q l EC ql
V EC =V CE = =
2 6
q l 2EC q l 2
M EC =M CE = =
12 108

Figura 3.23. Fuerzas en los extremos para la condición restringida: D 1=D 2 =D 3=0

Por lo tanto

{ }
−120 kN
{ R }= 0
80 kN . m
Las componentes de reacción en A cuando D1=D2 =D3=0

{ }{
−L

{ } }
−V AE 3 −80 kN
{ Ar }= 0 =q 0 = 0
2
−M AE −q L −106.67 kN .m
27

Paso 3 Las fuerzas en las tres coordenadas que se generan cuando se impone por
separado desplazamientos unitarios en cada coordenada (véase 3.22b, c y d) se

183
 Método de Análisis de Desplazamiento

determinan sumando las fuerzas en los extremos del nudo E (véanse las figuras 3.33a,
b y c), y que se han tomado del apéndice C; así se obtiene,

Para D 1=1 , D2=D3=0 (véase figura 3.24)

40.5 EI 40.5 EI 324 EI 324 EI 729 EI

K 11 =V EA +V EB+ V EC +V ED= 3
+ 3
+ 3 + 3 = 3
l l l l l
13.5 EI 54 EI −40.5 EI
K 21 =M ED−M EB = − 2 =
l2 l l2
−13.5 EI 54 EI 40.5 EI
K 31 =−M EA + M EC = 2
+ 2 = 2
l l l
Para D 2=1 , D1=D3=0 (véase figura 3.25)

−54 EI 13.5 EI −40.5 EI

K 12 =−V EB +V ED= + =
l2 l2 l2
12 EI 6 EI 3 GJ
K 22 =M EB+ M ED +T EA +T EC = + + + 3GJ
l l 2

Y teniendo en cuenta que GJ=0.5 EI

6 EI 12 EI 3 EI 3 EI 20.25 EI
K 22 = + + + =
l l 4l 2l l

K 32 =0

Para D 3=1 , D1=D2=0 (véase la figura 3.26)

−13.5 EI 54 EI 40.5 EI
K 13 =−V EA +V EC = 2
+ 2 = 2
L L L
K 23 =0
6 EI 12 EI 3 GJ
K 33 =M EA + M EC +T EB +T ED = + + +3 GJ
l l 2

Y teniendo en cuenta que GJ=0.5 EI

184
 Método de Análisis de Desplazamiento

6 EI 12 EI 3 EI 3 EI 20.25 EI
K 33 = + + + =
l l 4l 2l l

(a) MAE A MBE

VAE
B
V EA V EC
VBE

E MCE
V ED V EB
C
VDE VCE

D
MDE

Figura 3.24

TAE
(b)
A
MBE
B

VBE
TEA
MEB
MED
E TEC
V EB
V ED
C
MDE TCE

D
VDE

Figura 3.25

185
 Método de Análisis de Desplazamiento

(c) MAE
A
V AE V EA B TBE
MEA
TEB
TED E MEC

V EC MCE
C
VCE

D
TDE

Figura 3.26
Con los coeficientes de rigidez calculados arriba se escribe la matriz de rigidez

[ ]
729 −40.5 40.5
3 2 2
l l l
−40.5 20.25
[ K ] =EI 0
l2 l
40.5 20.25
0
l2 l

Después de reemplazar l=12 m obtenemos

[ ]
0.421875 −0.28125 0.28125
[ K ] =EI −0.28125 1.6875 0
0.28125 0 1.6875

[ ]
3.0476 0.5079 −0.5079
−1 1
[K] = 0.5079 0.6772 −0.0847
EI
−0.5079 −0.0847 0.6772

Los elementos de la matriz [ Au ] , correspondientes a las reacciones en A debido a los

desplazamientos unitarios, se hallan directamente de las figuras 3.24, 3.25 y 3.26. Así
se obtiene

186
 Método de Análisis de Desplazamiento

[ ]
−40.5 13.5
3
0 2

[ ]
l l
−0.0224 0 0.09375
−3
[ Au ] =EI 0
4l
0 =EI 0 −0.0625 0
−0.0938 0 0.25
−13.5 3
0
l2 l

Paso 4. Los desplazamientos {D} se calculan sustituyendo [ K ] y { R } en la ecuación 3.7 y

despejando { D } se tiene

[
1 3.0476
]{ } { }
−1
0.5079 −0.5079 120 1 406.3492
{ D }= [ K ] {−R }= 0.5079 0.6772 −0.0847 0 = 67.7249
EI EI
−0.5079 −0.0847 0.6772 −80 −115.1323

Paso 5 Las reacciones en el apoyo A se calculan con la ecuación 3.8

{ } [ ]{ }
−80 −0.0224 0 0.09375 1 406.3492
{ A}= 0 + EI 0 −0.0625 0 67.7249
EI
−106.67 −0.0938 0 0.25 −115.1323

{ }
−100.32kN
{ A } = −4.23 kN . m
−173.55 kNm

Los elementos de { A } representan las componentes de reacción vertical, el momento

torsor y el momento flector respectivamente en el apoyo A. Los signos negativos
significan que las fuerzas actúan en el sentido opuesto al indicado en el sistema de
coordenadas (véase figura 3.21b).

3.8 Análisis para efectos de desplazamientos en las coordenadas.

Considérese una estructura que tenga n grados de indeterminación cinemática. Se
trata de analizar la estructura cuando esté sometida a m desplazamientos impuestos de
magnitud Δ 1 , Δ 2 , . . . Δm en m coordenadas.

187
 Método de Análisis de Desplazamiento

Se genera la matriz de rigidez de modo que las coordenadas correspondientes a los

desplazamientos impuestos ocurren en las m primeras filas, así se tiene

[ ]
k 11 k 12 … k 1 m k 1 (m+1) … k1 n
k 21 k 22 … … … … …
… … … … … … …
[ K ] = k m 1 k m 2 … k mm k m(m +1) … k mn (3.9)
k m (m +1) … … … … … … …
… … … … … … … …
kn1 kn2 … … … … … k nn

Esta matriz se puede dividir a lo largo de las líneas punteadas que de muestran y por lo
tanto se puede escribir de la siguiente forma

[[
[ K ] = [ 11 ] [ 12]
K K
K 21 ] [ K 22] ] (3.10)

donde [ K ij ] son las submatrices. Inspeccionando la ecuación 3.10, se observa que el

orden de las submatrices es:

[ K 11 ]m ×m, [ K 12 ]m ×(n−m), [ K 21 ]( n−m)× m, [ K 22 ]( n−m)×(n−m)

Para imponer los desplazamientos Δ 1 , Δ 2 , . . . Δm se deben aplicar las fuerzas externas

¿ ¿ ¿
R1 , R 21 , … , R men las coordenadas 1 a m, sin fuerzas aplicadas en las demás
coordenadas, con la finalidad de que ocurran los desplazamientos D m +1 , Dm +2 , . . . Dn . La
suma de las fuerzas a lo largo de las coordenadas en las primeras mcoordenadas
¿
deben ser iguales a las fuerzas Ri 1 y en las restantes deben ser iguales a cero, por lo
tanto

[ ]{ } { }
[ K 11 ] [ K 12 ] { D¿ } = { R¿ }
[ K 21 ] [ K 22 ] { D } { 0 }
(3.11)

188
 Método de Análisis de Desplazamiento

donde { D¿ }={ Δ }, esto es, los desplazamientos impuestos; { D } es el vector de

desplazamientos desconocidos en las coordenadas m+1 , m+2 , … , n; y { R } es el vector
¿

de fuerzas desconocidas en las coordenadas 1 ,2 , … , m.

A partir de la segunda fila de la ecuación matricial de arriba, se encuentra

−1
{ D }=−[ K 22 ] [ K 21 ] { D¿ } (3.12)

Al calcular { D } se conocen los desplazamientos en las n coordenadas los

desplazamientos, luego las acciones en cualquier sección se pueden determinar con la
ecuación
{ A } =[ A u ] { D }
(3.13)

donde { A } es cualquier acción, y [ Au ] es la misma acción correspondiente a un

desplazamiento unitario en la coordenada dada solamente. Esta ecuación es la misma
ecuación 3.5 con { A r }= { 0 }, puesto que las acciones necesarias se deben solo al efecto
de los desplazamientos { D }
Las fuerzas { R¿ } se pueden calcular sustituyendo { D́ } encontrado con la ecuación 3. En
la primera fila de la ecuación 3; así

{ R¿ }={ [ K 11 ]−[ K 12 ][ K 22 ]−1 [ K 21 ] } { D¿ } (3.14)

La ecuación 3. Puede reescribirse en la forma

{ R¿ }=[ K ¿ ] { D¿ } (3.15)

donde [ K ¿ ] es la matriz de rigidez condensada de la estructura correspondiente a las

primeras m coordenadas, con la omisión de las coordenadas restantes, Por lo tanto, la
matriz de rigidez condensada es

189
 Método de Análisis de Desplazamiento

[ K ¿ ] =[ K 11 ] −[ K 12 ][ K 22 ]−1 [ K 21 ] (3.16)

Ejemplo 3.7: Estructura sometida a carga, temperatura y desplazamiento en una

coordenada
Encontrar el diagrama de momento flector del pórtico que se ilustra en la figura 3.27a
para: (a) una carga uniformemente distribuida de 12 kN/m sobre el tramo BD ; (b) un
incremento de temperatura de 40°C en las fibras inferiores y 20°C en las fibras
superiores del tramo BD y (c) una rotación horaria de B de 0.001 radianes. Ignore el
cambio de longitud de los miembros, suponer la rigidez a la flexión EI =108000 kN . m2
constante, el coeficiente de dilatación térmica α =1 ×1 0−5 /° C . y h=0.60 m (peralte de
los miembros)

(a) (b)

Figura 3.27 Pórtico considerado en el ejemplo 3.7. (a) Dimensiones y cargas. (b) sistema coordenado

Paso 1. El grado de indeterminación cinemática es dos (las rotaciones en B y C ). El

sistema coordenado se muestra en la figura 3.27b. Las acciones requeridas para
dibujar el diagrama de momento flector son

190
 Método de Análisis de Desplazamiento

[{ } { } { } ]
M AB M AB M AB
M BA M BA M BA
M BC M BC M BC
[ M ] = M CB M CB M CB M dc =0
M CE M CE M CE
M EC M EC M EC
( a) M CD a M CD b M CD c
MBC MCB MCD
B
C
D

Los subíndices a , b y c se refieren a los tres

A E casos de solicitación. Se considera que un
momento de extremo es positivo cuando
(b)
MBC MCB MCD
B Tsup=20ºc Tsup=20ºc
actúa en el sentido horario.
Tinf=40ºc C Tinf=40ºc
D

Paso 2 Para las solicitaciones

correspondientes a la carga y variación de
temperatura las fuerzas en los extremos
A E
empotrados se determinan con ayuda del
(c) apéndice B y las y las fuerzas restringentes
MBC MCB MCD
B
C
se calculan sumando las fuerzas de
MBA D
extremo en cada nudo (véase figura 3.28a
y b). Para la rotación de B no existe
problema primario puesto que esta
D1=1
MAD
A E solicitación corresponde a un
desplazamiento en la ubicación y dirección
(d)
MBC MCB de una coordenada.
B
C
MCD D
MCD

D2=1
MER M AB=M BA =0
A E
ql
2 191
M BC =M CB = =36 kN . m
12
M CE =M EC =0
 Método de Análisis de Desplazamiento

Figura 3.28. Momentos en los extremos de los miembros: (a) y (b) para las estructuras restringidas
sometido a carga y variación de temperatura. (c) y (d Cuando D 1=1 y D 2=1 , respectivamente
Caso ( a ): cargas (véase figura 3.28a)

R1=−M BC =−36 kN R2 =M CB −M CD =36−30.375=5.625 kN . m

Caso ( b ): Variación de temperatura (véase figura 3.26b)

M AB=M BA =0
R1=−M BC =−36 kN R2 =M CB −M CD =36−54=−18 kN . m
EIα
M BC =M CB = (T ¿ −T ¿ )=36 kN . m
h
Por lo tanto, la matriz
M CE =M EC =0
[ R ] para los casos (a) y (b) es
3 EIα
M CD = ( T ¿−T ¿ ) =54 kN . m 0
[ ]
−36 −36 2h
[ R ]= kN . m
5.625 −18

2 EI de B: { R } ={ 0 } 4 EI
Caso ( c ): rotación
M AB= =0.333 EI ; M BA = =0.667 EI
l l
4 EIde la matriz [ Ar ] se2 EI
Los elementos
M BC = =0.667 EI ; M CB = deducen=0.333fácilmente
EI de las figuras 3.28a y b,
l l

[ ]
M CE =M EC =0
0 0 0
M CD =0
0 0 0
−36 −36 0
[ M r ]= 36 −36 0 kN . m
0 0 0
0 0 0
M AB=M BA =0
−30.775 −54 0
2 EI 4 EI
M BC = =0.333 EI ; M CB = =0.667 EI
l l
Paso 3 Los4elementos
EI de la matriz
2 EI de rigidez se calculan sumando los momentos en
M CE = =0.667 EI ; M EC = =0.333 EI
los extremos lde los miembros necesarios l para producir las configuraciones deformadas
de lasMfiguras3 EI3.28c y d yEIcon ayuda del Apéndice C; así se tiene
CD = =0.6667
l
Para D 1=1 y D 2=0

192
 Método de Análisis de Desplazamiento

K 11 =M BA + M BC =0.6667 EI + 0.667 EI =1.333 EI K 21=M CB =0.333 EI

Para D2=1 y D 1=0

K 12 =M BC =0.333 EI

K 22 =M CB + M CE + M CD =0.667 EI +0.667 EI +0.667 EI =2 EI

Luego, la matriz de rigidez resulta

[ 0.333 EI
[ k ]= 1.333 EI 0.333 EI
2 EI ]
Los valores de los momentos en los extremos de los miembros debido a los
desplazamientos unitarios actuando por separado se determinan también de las figuras
3.26c y d, así se obtiene

[ ]
0.333 0
0.667 0
0.667 0.333
[ M u ]=¿ EI 0.333 0.667
0 0.667
0 0.333
0 0.667

Paso 4 Los desplazamientos [ D ] usando la ecuación 3.7 y despejando [ D ]

[0.333 EI ] [−5.625 ] = [ ]
−1
1 28.9076
[ D ] = 1.333 EI 0.333 EI 36 36 25.8261
2 EI 18 EI −7.6304 4.6957

Para el caso (c ), la ecuación de superposición está dada por la ecuación 3.11, así

[ 1.333 EI
0.333 EI 2 EI ]{ } { }
0.333 EI 0.008 = R¿1
D2 0

193
 Método de Análisis de Desplazamiento

De la segunda fila de la ecuación matricial anterior se tiene

0.333 EI × 0.001+2 EI × D 2=0, de la cual obtenemos: D2=−0.000167radianes

Y de la ecuación de la primera fila se tiene

1.333 EI ×0.001+0.333 EI × D2=R¿1, de la cual obtenemos: R¿1=0.001278 EI

¿
Y después de sustituir el valor de EI en el resultado anterior se obtiene R1=138 kN . m ,
que es el momento que hay que aplicar en la coordenada 1 para producir la rotación de
0.001 radianes en la misma coordenada. Finalmente, la matriz [ D ] resulta

[−7.6304 / EI
[ D ] = 28.9076 / EI 25.8261/EI
4.6957/ EI
0.001
−0.000167 ]
Paso 5. Los momentos de extremo finales se calculan por la ecuación 3.8.

[ ][ ]
0 0 0 0.333 0

[ ]
0 0 0 0.667 0 28.9076 25.8261
−36 −36 0 0.667 0.333 0.001
[ M ] = 36 EI EI
−36 0 + EI 0.333 0.667
−7.6304 4.6957
0 0 0 0 0.667 −0.000167
EI EI
0 0 0 0 0.333
−30.775 −54 0 0 0.667

[ ]
9.64 8.61 36
19.27 17.22 72
−19.27 −17.22 66
[ M ] = 40.55 47.74 24 kN . m
−5.09 3.13 −12
−2.54 1.57 −6
−35.46 −50.87 −12

Cada columna representa los momentos de extremo para cada caso. Los diagramas
correspondientes se muestran en la figura 3.27.

194
 Método de Análisis de Desplazamiento

50.87
40.55 47.74
35.46
(a) 19.27 (b) 17.22

ql 2/8=54
B 30.38 B
19.27 17.22
C 5.09 3.13 C
D D
12.24
24.09

A 9.64 2.54 E A 8.61 E 1.57

DMF en kN.m DMF en kN.m

(c) 24
12
72 B
C 12
D

66

A 36 6 E

DMF en kN.m

Figura 3.29 Diagramas de momento flector. (a) Caso de carga. (b) caso de variación de temperatura. (c)
Caso de desplazamiento en una coordenada

3.9 Análisis de estructuras simétricas

Se puede aprovechar la simetría para reducir el grado de indeterminación cinemática
cuando el análisis se realiza por el método de desplazamiento. Debido a la simetría, la
magnitud del desplazamiento en una coordenada es cero o son iguales en magnitud en
una (o más) coordenadas. Para un desplazamiento cero, la coordenada puede; en el
caso de que en dos coordenadas las magnitudes de desplazamiento sean iguales se
pueden representar por un solo número. La figura 3.20 muestra ejemplos de sistemas
de coordenadas las cuales pueden utilizarse en el análisis de estructuras simétricas
sometidas a carga simétrica usando el método de desplazamiento.

195
 Método de Análisis de Desplazamiento

Las figuras 3.30a-d representan pórticos planos en donde se desprecian las

deformaciones axiales. En la figura 3.30a, debido a la simetría de carga, las rotaciones
en B y C son iguales y opuestas y no hay desplazamiento lateral, por lo tanto, el pórtico
tiene un grado de libertad que está representado por la coordenada 1. En la figura
3.30b, debido a la simetría, las rotaciones en B y F son iguales y opuestas, mientras

que la rotación en D y el desplazamiento lateral son nulos; por lo tanto, el pórtico tiene
solo un desplazamiento desconocido.

196
 Método de Análisis de Desplazamiento

Figura 3.30 Ejemplos de grados de libertad para el análisis por el método de desplazamientos de
estructuras simétricas cargadas simétricamente.
No se muestra un sistema de coordenadas para las estructuras de las figuras 3.30c y d,
pues todos los nudos bajo carga simétrica no sufren ningún desplazamiento. Así, no se
necesita analizarlas y las fuerzas en los extremos de los miembros se pueden calcular
fácilmente usando el Apéndice B.
La viga sobre apoyos elásticos de la figura 3.30e solo tiene dos desplazamientos
desconocidos: las traslaciones verticales en A y C , que son iguales y la traslación
vertical en B. Debido a la simetría de carga la rotación en B es cero. Además, las
rotaciones en A y C serán iguales y opuestas y no constituyen grados de libertad si al
calcular las rigideces se tiene en cuenta que los miembros tienen el extremo articulado.
Si en el análisis se consideran las deformaciones axiales, se deben usar coordenadas
adicionales para representar los desplazamientos. En la figura 6.28a, además de las
rotaciones en B y C , ocurrirán los desplazamientos horizontales y verticales, pero
debido a la simetría los desplazamientos son de igual magnitud y cada uno se puede
representar con una sola coordenada; por lo tanto, el pórtico tiene tres desplazamientos
desconocidos. En el caso de las estructuras representadas en las figuras 3.30c y d se
tiene un solo desplazamiento desconocido: la traslación vertical del nudo B.
Los emparrillados horizontales mostrados en las figuras 3.30f y g, tienen uno o más
planos de simetría. El sistema de coordenadas mostrado puede ser usado para el
análisis cuando estas estructuras estén sometidas a una carga simétrica.
Cada una de las estructuras mostradas en las figuras 3.30a-g se pueden analizar
considerando solo la mitad (o un cuarto) de la estructura separándolas mediante el
plano de simetría. Los miembros situados simétricamente a la parte analizada deben
tener sus propiedades geométricas y mecánicas iguales ( A , I , J y E ).
Cuando el análisis de una estructura se realiza en forma manual, la generación de la
matriz de rigidez para la parte simétrica de la estructura puede que no incorpore una
ventaja apreciable, pero, en estructuras grandes con muchos miembros y nudos, el
análisis se realiza íntegramente por computadora y el aprovechamiento de la simetría
en el cálculo de la estructura representa una considerable ventaja.

197
 Método de Análisis de Desplazamiento

Ejemplo 3.8 Estructura simétrica con carga simétrica

Encontrar el diagrama de momento flector para el pórtico de la figura 3.31a. Además de
la carga el apoyo F se asienta 5 mm. La rigidez a la flexión EI =1.4 ×10 5kN.m2 es
constante y los miembros no se deforman axialmente.

Paso1 Sin deformaciones axiales de los miembros, debido a la simetría las rotaciones
en B y C son iguales y opuestas y en el nudo C no existe giro. Las rotaciones en A y D
no son desplazamientos independientes pues se trata de miembros con extremo
articulado y el pórtico tiene un desplazamiento desconocido: la rotación en B ( θ B ¿.
Luego se puede resolver solo la mitad de la estructura (véase la figura 3.31b). El
sistema coordenado se indica en la figura 3.29c. Los momentos de extremo requeridos
son
{ M }= { M BA , M BC M CB }

Los momentos en la mitad derecha son iguales y opuestos al valor correspondiente de

la mitad izquierda. Se considera que un momento en el extremo, es positivo si actúa en
el sentido horario.

Paso 2 El momento restringente R , se obtiene sumando los momentos en los extremos

empotrados, que se toman de apéndices B y C; en la ubicación y dirección de la
coordenada 1 (véase la figura 3.31d), así se tiene

R=67.5−293=−225.5 kN . m

Los momentos en los extremos { M r } se determinan fácilmente de la figura 3.31d

{ M r }={ 67.5 , −293 , −43 } kN .m

198
 Método de Análisis de Desplazamiento

(a) 60 kN /m (b) 60 kN/m

B C D B C
00

00
4
.

4
.

A F E A
m
5m

3.00 5.00
3.00 5.00 5.00 3.00

(c)
1
B C

(d) 60 kN/m

B MBC MCB C
MBA

( e)
MCB
MBA
B C
MBC

q l 2ABH
A
M BA = =67.5
8
2
−q l BC 6 EIδ
M BC = − 2 =−293 kN . m
12 l BC
2
q l BC 6 EIδ
M CB = − 2 =−43 kN .m
12 l BC

3 EI
M BA = =0.6 EI
5
4 EI
M BC = =0.8 EI
5
2 EI
M CB = =0.4 EI
5

199
 Método de Análisis de Desplazamiento

Figura 3.31. Pórtico simétrico del ejemplo 3.8. (a) Geometría y cargas. (b) Estructura equivalente. (c)
Sistema coordenado. (d) Fuerzas en los extremos de la estructura restringida sometida a las
solicitaciones externas. (e) Fuerzas en los extremos debido al desplazamiento unitario D 1=1

Paso 3 Las fuerzas en Loa extremos de los miembros inducidos por un desplazamiento
unitario en la coordenada 1 se determinan del apéndice C y se indican en la figura
3.31e, y el coeficiente de rigidez, se obtiene sumando las fuerzas de extremo en la
dirección y ubicación de la coordenada 1

k 11=0.6 EI +0.8 EI =1.4 EI

Los elementos de { M u } se escriben fácilmente de la figura 3.31e

{ }
0.6 EI
{ M u }= 0.8 EI
0.4 EI

Paso 4 Los desplazamientos se obtienen Sustituyendo los valores de R y k 11 en la

ecuación 3.7, y despejando { D } se tiene

−1 −1
{ D }= [ k ] { R }=[ 1.4 EI ] { 225.5 }= {161.07
EI }

Paso 5 Los momentos de extremo finales se calculan por la ecuación 3.8

{ } { }{ } { 21.43 }
67.5 0.6 EI 164.14
{ M }= { M r }+ { M u } { D }= −293 + 0.8 EI 161.07 = −164.14 kN . m
EI
−43 0.4 EI

200
 Método de Análisis de Desplazamiento

El diagrama de momento flector se muestra en la figura 3.32.

164.14 164.14
164.14 21.43 164.14
2
ql /8=187.5
B C D

94.72 94.72

A F E

Figura 3.32. Diagrama de momento flector (kN.m)

Ejemplo 3.9: Estructura simétrica con carga antisimétrica

Obtener el diagrama de fuerza cortante y momento flector para el pórtico mostrado en
la figura 3.33a. La rigidez a la flexión EI para cada miembro se indica en la figura. No
se consideran las deformaciones axiales

Paso1 Se tiene el caso de una estructura simétrica con carga antisimétrica, por lo tanto,
las rotaciones en B y C son iguales. Luego, el pórtico tiene un desplazamiento
desconocido (θ B=θ C ¿ y se analiza solo la mitad, puesto que las fuerzas en la mitad
derecha son iguales al valor correspondiente de la mitad izquierda. El sistema
coordenado se indica en la figura 3.33b. Las acciones requeridas son las fuerzas
cortantes y los momentos flectores de extremo.

{ V }= { V AB ,V BA , V BC , V BE , V EB } { M }= { M BA , M BC , M BE , M EB }

Se considera que el momento de extremo es positivo cuando actúa en el sentido

horario y la fuerza cortante es positiva cuando actúa hacia arriba o hacia la derecha.

201
 Método de Análisis de Desplazamiento

Paso 2 El momento restringente R , se obtiene sumando los momentos en los extremos

empotrados, que se toman de apéndices B y C; en la ubicación y dirección de la
coordenada 1 (figura 3.33c), así se tiene

R=M BA + M BC + M BE =0+0+72=72kN . m

Los elementos de { V r } y { M r } se determinan fácilmente con las fuerzas de extremo de la

figura 3.31c

{V r }= {0 0 0 −72 −72 } kN { M r }= { 0 0 72 −72 } kN . m

(a) (b)
3.00 3.00 3.00
A B C D A 1 1 D
EI EI EI B C
2 4 k N /m
2 4 k N /m

6 .0 0
3EI

3EI

E F E F

(c)
A B C D
V AB=V BA =0
MBE
VBE M BA =0
V BC =0
VEB D1=0
MEB M BC =0
E F
q l BE
V BE=V EB= =72 kN.m
2
2
q l BE
M BE =M EB= =72kN . m 202
12
 Método de Análisis de Desplazamiento

3 EI
V AB=V BA = =0.333 EI
l 2AB
(d) 3 EI
MBA M BA = =EI
A D l AB
B C
MBC
MBE
VBE −12 EI
V BC = =−1.333 EI
l 2BC
D1=1
6 EI
E VEB M BC = =2 EI
MEB
l BC
F

6 ( EI )BE
V BE=V EB= 2
=0.5 EI
l BE
4 ( EI ) BE
M BE = =2 EI
l BE
Figura 3.33. Pórtico analizado en el ejemplo 3.9. (a) Dimensiones y carga. (b) Sistema coordenado. (c)
Fuerzas de extremo debido a las cargas. (d) Fuerzas de extremo debido a D 1=1

203
 Método de Análisis de Desplazamiento

Paso 3 Las fuerzas en los extremos de los miembros inducidos por un desplazamiento
unitario en la coordenada 1 se determinan del apéndice C y se indican en la figura
3.33d, y el coeficiente de rigidez, se obtiene sumando las fuerzas de extremo en la
dirección y ubicación de la coordenada 1

k 11=M BA + M BC + M BE=EI +2 EI +2 EI=5 EI

Los elementos de { V u } y { M u } se escriben fácilmente con las fuerzas de extremo de la

figura 3.33d

{ } {}
−0.333
1
0.333
2
{ u}
V =EI −1.333 { M u } =EI
2
−0.5
1
0.5

Paso 4 El desplazamiento se obtiene sustituyendo los valores de R y k 11 en la ecuación

3.7 y despejando { D } resulta

−1 1 −14.4
{ D }= [ k ] {−R }= (−72 )=
5 EI EI

Paso 5 Las fuerzas finales de extremo se calculan usando la ecuación 3.8. Luego

{}{ } { }
0 −0.333 4.8
0 0.333 −14.4 −4.8
{ V }= { V r } + {V u } { D }= 0 + EI −1.333
−72 −0.5
{EI }
= 19.2 kN
−64.8
−72 0.5 −79.2

{ } { }{ } { }
0 1 −14.4
{ M }= { M r }+ { M u } { D }= 0 + EI 2 −14.4 = −28.8 kN . m
72 2 EI 43.2
−72 1 −86.4

204
 Método de Análisis de Desplazamiento

El diagrama de fuerza cortante y momento flector se muestra en la figura 3.34

28.80
(a) 19.20 (b) 14.40
4.80 4.80
A 43.20 43.20 D
64.80 B 64.80 C D B C
A
14.40
28.80

43.20 43.20
ql2/8=108

E 79.2 F 79.2 86.40 E 86.40 F

Figura 3.34. (a) Diagrama de fuerza cortante (kN). (b) Diagrama de momento flector (kN.m)

3.10 Consideración de las deformaciones axiales

En una estructura compuesta por elementos tipo barras (estructuras reticulares), las
fuerzas axiales cambian las longitudes de los miembros, así como sus valores de
rigidez. En todos los ejemplos de este capítulo, a excepción de las armaduras, se ha
despreciado el efecto de las fuerzas axiales. Cuando el valor de las fuerzas axiales es
alto es necesario tomarlos en consideración en el análisis.
En esta parte solo se considera el efecto de las fuerzas axiales en el cambio de
longitud de un miembro, más no su cambio de rigidez. Con esta limitación cuando se
analiza un pórtico con el método de desplazamiento, las incógnitas se incrementan por
que aparecen grados de libertad correspondientes a la traslación en cada nudo; por
ejemplo, un nudo sin apoyo de un pórtico plano tiene tres grados de libertad: una
rotación y dos traslaciones; y un nudo libre de un pórtico tridimensional tendrá seis
grados de libertad: tres rotaciones y tres traslaciones

Ejemplo 3.10: Análisis de un pórtico considerando las deformaciones axiales

205
 Método de Análisis de Desplazamiento

Obtener los diagramas de fuerza normal, fuerza cortante y momento flector para el
pórtico que se muestra en la figura 3.35 tomando en consideración el cambio de
longitud de los miembros. Tómese EA=7.2× 106kN y EI =4.86 × 105 kN.m2

P=900kN
(a) 30 kN/m (b) 2
1
B B 3
C C

15

A A
15

Figura 3.35. Análisis del pórtico del ejemplo 3.10. (a) Propiedades y carga del pórtico. (b) Sistema
coordenado

Paso1 Si se consideran las deformaciones axiales, el nudo B tiene tres

desplazamientos independientes como se muestra en la figura 3.35b que también
ilustra el sistema coordenado. Para dibujar los diagramas de fuerzas internas se
requieren las fuerzas en los extremos de todos los miembros.

{} {} { }
N AB V AB
M AB
N V BA
{ N }= BA { V }= { M }= M BA
N BC V BC
M BC
N CB V CB

Además, se observa que M CB =0.

En el extremo de un miembro, el momento en el sentido horario y la fuerza normal o
cortante dirigidas hacia arriba o hacia la derecha se consideran positivas.

Paso 2 Las fuerzas restringentes { R } se obtienen sumando las fuerzas de extremo en la

ubicación y dirección de las coordenadas en la figura 3.36 así se tiene

206
 Método de Análisis de Desplazamiento

R1=0 R 2=P+V BC =900+281.25=1181.25 kN R3 =M BC =843.75 kN . m

Por lo tanto
{ R } ={ 0 1181.25 kN 843.75 kN . m }

900kN
30 kN/m MBC
5 ql
C V BC = =281.25 kN
B
8
VCB VBC
3 ql
V CB = =168.75 kN
D1=D2=D3=0 8
2
ql
A M BC = =843.75 kN . m
8

Figura 3.36. Fuerzas en los extremos en la estructura restringida sometida a las cargas

Los elementos de { N r }, { V r } y { M r } se determinan fácilmente con las fuerzas de extremo

de las figura 3.36; así,

{} { }
0 0

{ }
0
0 0
{ N r }= 0 { V r } 281.25 kN { M r } 0 kN . m
843.75
0 168.75

Paso 3 Las fuerzas en los extremos de los miembros que corresponden a un

desplazamiento unitario separado en cada una de las coordenadas se determinan
ahora por el apéndice C y se indican en las figuras 3.37, 3.38 y 3.39. Las fuerzas en la
ubicación y dirección de las coordenadas son iguales a la suma de las fuerzas de
extremo.

207
 Método de Análisis de Desplazamiento

Para D 1=1 , D2=D3=0 (Véase figura 3.37)

k 11=V BA + N BC =1728+480000=481728
K 21 =0
K 31 =−M BA =−12960

N AB=N BA=0
12 EI
NCB NCB 1
V AB=V BA = =1728
C
B
MBA l 3AB
VBA

6 EI
M AB=M BA = 2
=12960
D1=1 l AB
EA
N BC =N CB = =480000
A VAB l BC
MAB

V BC =V CB =0

Figura 3.37. Fuerzas en los extremos de la estructura cuando D 1=1

Para D 2=1 , D1=D3=0 (Véase figura 3.38)

k 12=0
K 22 =N BA +V BC =480000+ 432=480432
K 32 =M BC =6480
EA
N AB=N BA= =480000
MBC l AB
C NBA
B V AB=V BA =¿0
VBC
VCB
M AB=M BA =0
D2=1 N BC =N CB =0
3 EI
V BC =V CB = 3
=432
A
l BC
NAB 3 EI
M BC = =6480
l 2BC

Figura 3.38. Fuerzas en los extremos de la estructura cuando D 2=1

208
 Método de Análisis de Desplazamiento

D3=1 , D 1=D 2=0 (Véase la figura 3.39)

k 13=V BA =−12960
K 23 =V BC =6480
K 33 =M BA + M BC =129600+97200=226800

N AB=N BA=0
6 EI
V AB=V BA = 2
=12960
l AB
2 EI
MBC

C
M AB= =64800
VBC B VBA
l AB
VCB MBA

4 EI
M BA = =129600
l AB
VAB A N BC =N CB =0
MAB

3 EI
V BC =V CB = =432
l 3BC

Figura 3.39. Fuerzas en los extremos de la estructura cuando D 3=1

Luego, la matriz de rigidez y su inversa son

[ ]
481728 0 −12960
[ K ]= 0 480432 6480
−12960 6480 226800

[ ]
−5 −9 −7
2.0791× 10 −1.6030× 10 1.1885 ×10
−1
[ K ] = −1.6030× 10 −9
2.0823× 10
−6
−5.9585 ×10
−8

−7 −8 −6
1.1885× 10 −5.9585× 10 4.4177× 10

Las matrices [ N u ] , [ V u ] y [ M u ], se escriben sin dificultad con las fuerzas de extremo de

las figuras 3.37, 3.38 y 3.39

209
 Método de Análisis de Desplazamiento

[ ] [ ]
0 −480000 0 −1728 0 12960
0 480000 0 1728 0 −12960
[ N u ]= 480000 0 0
[ V u ]= 0 432 6480
−480000 0 0 0 −432 −6480

[ ]
−12960 0 64800
[ u ] 12960
M = 0 129600
0 6480 97200

Paso 4. Los desplazamientos se hallan sustituyendo [ K ] y { R } en la ecuación 3.7

[ ]{ }
−5 −9 −7
2.0791× 10 −1.6030 ×10 1.1885× 10 0
{ D }= −1.6030 ×10−9 2.0823× 10−6 −5.9585× 10−8 1181.25
−7 −8 −6
1.1885× 10 −5.9585× 10 4.4177 ×10 843.75

{ }
−5
−9.8385 ×10 m
{ D }= −2.4094 ×10−5 m
−5
−3.6570 ×10 rad

Paso 5 Las fuerzas finales de extremo se calculan por la ecuación 3.8

[ ]{ }{ }
0 −480000 0 −9.8385 ×10−5 1156.51
{ N }= 0 480000 0 −1156.51
−2.4094 ×10 = kN
−5
480000 0 0 −5 −47.22
−480000 0 0 −3.6570 ×10 47.22

Fuerza cortante: { V }= { V r } + {V u } { D }

{ }[ ]{ }{ }
0 −1728 0 12960 −5 −47.22
−9.8385× 10
{ V }= 0 + 1728 0 −12960 47.22 kN
−2.4094 × 10 =
−5
281.25 0 432 6480 −5 256.51
168.75 0 −432 −6480 −3.6570× 10 193.49

Momento flector: { M }= { M r }+ { M u } { D }

210
 Método de Análisis de Desplazamiento

{ }[ ]{ }{ }
−5
0 −12960 0 64800 −9.8385 ×10 −235.70
{ M }= 0 + 12960 0 129600 −2.4094 ×10−5
=¿ −472.67 kN . m
843.75 0 6480 97200 −3.6570 ×10 −5
472.67

Los diagramas de fuerza internas se muestran en la figura 3.40

193.49
(a) (b)
C B 1156.51 B 47.22
(-)
47.22 47.22 C

256.51
(-)

A A 47.22
1156.51

Figura 3.40. Diagramas. (a) Fuerza normal en kN.

472.67
(b) Fuerza cortante en kN. (c) Momento flector
(kN/m)
ql2/8=843.75 472.67
B
C

Ejemplo 3.11 607.41

Construir los diagramas de fuerza

normal, fuerza 235.70 A cortante y momento
flector para el pórtico que se muestra en la
figura 3.41, tomando en cuenta el cambio
de longitud de los miembros. Junto a
las cargas indicadas el miembro inclinado sufre un incremento de temperatura de 10°C
en la cara interior y de 40°C en la cara exterior. Las rigideces axiales EA y a la flexión
EI son constantes. Tómese E=2× 106 Ton. m2, sección transversal de las barras:
0.50m x 0.50 m y α =10−5 ° /C

211
 Método de Análisis de Desplazamiento

Paso 1 Si se considera el cambio de longitud de los miembros, el pórtico tiene tres

grados de libertad: dos traslaciones y una rotación. El sistema de coordenadas se
muestra en la figura 3.41b. Para dibujar los diagramas de fuerzas internas se requieren
las fuerzas en los extremos de todos los miembros.

Figura 3.41 Pórtico con deformaciones axiales. (a) Propiedades y carga. (b) sistema coordenado

{}{} {}
N AB V AB M AB
N V BA M BA
{ N }= BA { V }= { M }=
N BC V BC M BC
N CB V CB M CB

En el extremo de un miembro, el momento en el sentido horario, la fuerza normal y

cortante dirigidas hacia arriba o hacia la derecha se consideran positivas

Paso 2. Las fuerzas de extremo para la estructura restringida cuando sobre ella actúan
las cargas y la variación de temperatura se toman del apéndice B. Las fuerzas
restringentes { R } se obtienen sumando las fuerzas de extremo en la ubicación y
dirección de las coordenadas de las figuras 3.42 y 3.43 así se tiene

a) Por temperatura (véase la figura 3.42)

NBC N AB=N BA=0

A B
V AB=V BA =0
MBC 40ºC
M AB=M BA =0
10ºC

MCB
C
NCB
N BC =N CB =α ( t +2 t ) AE=125 kN
e i

V BC =V CB =0 212

M BC =M CB =α ( )
t e−t i
h
EI =6.25 kN . m
 Método de Análisis de Desplazamiento

Figura 3.42. Fuerzas de extremo en la estructura restringida por efecto de la temperatura

R1=N BC cosθ=125 ( 35 )=75

senθ=−125 ( ) =−100
4
R2=−N BC
5
R3=M BC =6.25

{ }
75
{ R } ΔT = −100
6.25

b) Caso de las cargas (véase la figura 3.43)

N AB=N BA=0
MAB 2t/m NBC
VBC q l AB
A B MBC 10t V AB=V BA = =6 kN
2
VAB VBA
q l 2AB
M AB=M BA = =6 kN . m
12
VCB
A P
MCB
N BC =N CB = senθ=4 kN
NCB 2
P
V BC =V CB = cosθ=3 kN
2
Pl BCH
M BC =M CB = =3.75 kN
Figura 3.43. Fuerzas en los extremos de la estructura restringida. Por efecto de las cargas 8

R1=0
P
R2=V BA + =11
2
R3=M BA + M BC =6−3.75=2.25

{ }
0
{ R }carga = 11
2.25

213
 Método de Análisis de Desplazamiento

{ }{ }{ }
75 0 75 Ton
{ R } ={ R }ΔT + { R }carga +¿ −100 + 11 = −89Ton
6.25 2.25 8.5 Ton. m

Los elementos de las matrices { N r } , {V r } y { M r } se escriben fácilmente sumando las

fuerzas de extremo de la estructura restringida solicitada por la variación de
temperatura y cargas (figuras 3.42 y 3.43)

{ } { } {}
0 0 6
0 = 0 Ton { V r } = 6 Ton
{ N r }= 125−4 121 3
−125−4 −129 3

{ }{ }
−6 −6
6 6 Ton . m
{ M r }= 6.25−3.75 2.5
−6.25−3.75 −2.5

Paso 3 Las fuerzas de extremo en los miembros cuando se impone a la estructura

desplazamientos unitarios por separado en la ubicación y dirección de las
coordenadas, se toman del apéndice C. Las fuerzas en la ubicación y dirección de las
coordenadas, necesarias para mantener las configuraciones deformadas de las figuras
3.44, 3.45 y 3.46 se obtienen sumando las fuerzas en los extremos.

Para D 1=1 , D2=D3=0. De la figura 3.44 se deduce:

3 4
δ AB=1 Δ AB =0 δ BC = ΔBC =
5 5

NBC
NAB B' MBC
A B 
VBC

D1=1
 C MCB
VCB
NCB

214
 Método de Análisis de Desplazamiento

Figura 3.44. Fuerzas en los extremos de la estructura cuando D 1=1

k 11=N AB+ N BC cosθ+V BC senθ= ( ) ( )( 45 )= 179750EA + 192

EA 3 EA 3
6
+
25 5
+
48 EI
625
EI
3125

25 ( 5 ) 625 ( 5 )
−3 EA 4 48 EI 3 −12 EA 144 EI
k 21=−N BC senθ+ V BC cosθ= + = +
125 3125
−24 EI
k 31=−M BC =
125
Para D 2=1 , D1=D3=0. De la figura 3.45 se deduce:
4 3
δ AB=0 Δ AB=1 δ BC = Δ BC =
5 5

N AB=N BA=0
MBA
NBC B' MBC 12 EI EI
MAB
VBC V AB=V BA = Δ AB=
A


B
3
l AB 18
VAB

6 EI EI
M AB=M BA = Δ AB=
D2=1
2
l AB 6
 C MCB

()
VCB
NCB EA EA 4 4 EA
N BC =N CB = δ = =
l BC BC 5 5 25

V BC =V CB =
12 EI
3
l BC
Δ BC =
125 5()
12 EI 3 36 EI
=
625

Figura 3.45. Fuerzas en los extremos de la estructura cuando D 2=1

k 21=−N BC cosθ+V BC senθ=

25 5() ()
4 EA 3 36 EI 4 −12 EA 144 EI
+
625 5
=
125
+
3125

18 25 ( 5 ) 625 ( 5 ) 125
EI 4 EA 4 36 EI 3 16 EA
k 22=V BA + N BC senθ+V BC cosθ= + + = +0.090116 EI

EI 18 EI 17 EI
k 32=M BA−M BC = − =
6 125 750

Para D 3=1 , D1=D2=0. De la figura 3.46 se deduce:

δ AB=0 Δ AB=δ BC =0

215
 Método de Análisis de Desplazamiento

N AB=N BA=0
6 EI EI
V AB=V BA = =
2
l AB 6
MAB MBA

2 EI EI 4 EI 2 EI
M AB= = M BA = =
A B
VAB VBC MBC

D3=1
VBC
l AB 3 l AB 3
 C
VCB

N BC =N CB =0
MCB

6 EI 6 EI
V BC =V CB = =
2
l BC 25

Figura 3.46. Fuerzas en los extremos de la estructura cuando D 3=1

k 13 =−V BC senθ=
25 5 ()
−6 EI 4 −24 EI
=
125

k 23 =V BA −V BC cosθ=
EI 6 EI 3 17 EI
6
−
25 5
=
750 ()
2 EI 4 EI 22 EI
k 33=M BA + M BC = + =
3 5 15

Matriz de rigidez

[ ]
179 EA 192 EI −12 EA 144 EI −24 EI
+ +
750 3125 125 3125 125
12 EA 144 EI 16 EA 17 EI
[ K ]= − + +0.090116 EI
125 3125 125 750
−24 EI 17 EI 22 EI
125 750 15

Reemplazando EA=500000 t y EI=10416.6667tm 2 en la matriz anterior se obtiene

[ ]
119973.333 −47520 −2000
[ K ] = −47520 64938.704 236.111
−2000 236.111 15277.778

[ ]
−5 −6 −6
1.1767× 10 8.6058 ×10 1.4074 ×10
−1
[ K ] = 8.6058× 10−6 2.1694 ×10−5 7.9131× 10−7
−6 −7 −5
1.4074 × 10 7.9131× 10 6.5627 ×10

216
 Método de Análisis de Desplazamiento

Los elementos de las matrices [ N u ] , [V u ] y [ M u ] escriben fácilmente con las fuerzas de

extremo de las figuras 3.44, 3.45 y 3.46 y después de reemplazar EA y EI por sus
valores se tiene

[ ] [ ]
−83333.33 0 0 0 −578.70 −1736.11
{ N u }= 83333.33 0 0 V = 0
{ u} 578.70 1736.11
60000 −80000 0 800 600 −2500
−60000 80000 0 −800 −600 2500

[ ]
0 1736.11 3472.22
0 1736.11 6944.44
{ M u }= 2000 −1500 8333.33
−2000 −1500 4166.67

Paso 4 Los desplazamientos { D } se obtienen sustituyendo los valores de [ K ] y { R } en la

ecuación 3.7

[ ]{ }
−5 −6 −6
1.1767 ×10 8.6058 ×10 1.4074 ×10 −75
−1
{ D }= [ K ] {−R }= 8.6058 ×10 −6
2.1694 ×10
−5
7.9131 ×10
−7
89
−6 −7 −5
1.4074 ×10 7.9131 ×10 6.5627 ×10 8.5

{ }
−1.2859 ×10−4 m
{ D }= 1.2786 ×10−3 m
−4
−5.9296 ×10 rad

Paso 5 Las fuerzas finales en los extremos de los miembros se obtienen por
superposición usando la ecuación 3.8

Fuerza normal: { N }= { N r } + { N u } { D }

{ }[ ]{ }{ }
0 −83333.33 0 0 −1.2859× 10−4 10.72
{ N }= 0 + 83333.33 0 0 1.2786 ×10−3 = −10.72 Tonf
121 60000 −80000 0 −4 11
−129 −60000 80000 0 −5.9296× 10 −19

Fuerza cortante: { V }= { V r } + {V u } { D }

217
 Método de Análisis de Desplazamiento

{}[ ]{ }{ }
6 0 −578.70 −1736.11 −1.2859 ×10−4 6.29
6
{ V }= + 0 578.70 1736.11 1.2786 × 10−3
= 5.71 Tonf
3 800 600 −2500 −4 5.15
3 −800 −600 2500 −5.9296 × 10 0.85

Momento flector: { M }= { M r }+ { M u } { D }

{ }[ ]{ }{ }
−6 0 1736.11 3472.22 −1.2859 ×10−4 −5.84
{ M }= 6 + 0 1736.11 6944.44 1.2786 ×10−3 = 4.10 Tonf . m
2.5 2000 −1500 8333.33 −4 −4.10
−2.5 −2000 −1500 4166.67 −5.9296 ×10 −6.63

Los diagramas de fuerzas internas se muestran en la figura 3.47

218
 Método de Análisis de Desplazamiento

Figura 3.47. Diagramas de fuerzas internas. (a) Fuerza normal. (b) Fuerza cortante. (c) Momento flector

3.11 Uso del método de fuerza o de desplazamiento

Después de haber los dos métodos básicos de análisis: Fuerza y desplazamiento, y
obviamente se puede usar cualquiera de ellos para analizar cualquier estructura. No
obstante, es conveniente establecer las ventajas de cada método en particular.

3.11.1 Relación entre las matrices de flexibilidad y rigidez

Considere un sistema de n coordenadas definidas en una estructura en la que es
aplicable el principio de superposición (ver Sección 1.6). Se ha definido la flexibilidad
matriz, [ f ] y sus elementos, los coeficientes de flexibilidad. Además, se ha definido la
matriz de rigidez, [ k ] y sus elementos, los coeficientes de rigidez. Un coeficiente de
flexibilidad, f ij es el desplazamiento en la coordenada i debido a la fuerza F j=1 en la
coordenada j , sin fuerzas aplicadas a las demás coordenadas. Un coeficiente de
rigidez, k ij es la fuerza en la coordenada i cuando la estructura se deforma tal que D j =
1 mientras que los desplazamientos son cero en las coordenadas restantes. A
continuación, se establece la relación entre las matrices de flexibilidad y rigidez,
correspondientes al mismo sistema de coordenadas.

(3.17)
[ k ]= [ f ]−1

[ f ] =[ k ]−1 (3.18)

Las Ecuaciones 3.17 a 3.21 se aplican con cualquier número de coordenadas. Para
facilitar el estudio, suponer un sistema de dos coordenadas ( n=2 ) (véase la figura
3.48a). Los desplazamientos { D1 , D 2 } se pueden expresar como la suma de los
desplazamientos debidos a las fuerzas F 1 y F 2 actuando separadamente

D1=f 11 F 1+ f 12 F 2
D2=f 21 F 1+ f 22 F 2

Estas ecuaciones se pueden escribir en forma de matriz

219
 Método de Análisis de Desplazamiento

{ D }= [ f ] { F } (3.19)

Las fuerzas { F } se pueden expresar en términos de los desplazamientos resolviendo la

ecuación anterior

−1
{ F }=[ f ] { D } (3.20)

Se puede usar la ecuación 3.20 para determinar las fuerzas que forman los elementos
de la matriz de rigidez de la estructura. Si se deforma por las fuerzas F 11 y F 21 a lo largo
de las coordenadas de tal manera que el desplazamiento D 1=1, mientras D 2=0 (véase
la figura 348b), entonces

{ } {}
F11
F21
=[ f ]
−1 1

Figura 3.48 Estructura con dos coordenadas. (a) Sistema de coordenadas. (b) Fuerzas manteniendo la
estructura en una configuración deformada con D 1=1, mientras que D 2=0 . (C) Fuerzas manteniendo la
estructura en una configuración deformada con D 2=1 , mientras que D 1=0

Se puede usar la ecuación 3.20 para determinar las fuerzas que forman los elementos
de la matriz de rigidez de la estructura. Si se deforma por las fuerzas F 11 y F 21 a lo largo
de las coordenadas de tal manera que el desplazamiento D 1=1, mientras D 2=0 (véase
la figura 348b), entonces

220
 Método de Análisis de Desplazamiento

{ } {}
F11
F21
=[ f ]
−1 1

Asimismo, las fuerzas necesarias para mantener la estructura en la configuración

deformada con D2=1 mientras D1=0 (figura46c) son:

{ } [] {}
F12
F22
=f
−1 0
1

Los resultados anteriores se pueden combinar en una sola ecuación matricial:

[ F11 F 12
F21 F 22 ] [ ]
=[ f ]
−1 1 0
0 1

Los elementos de la matriz en el lado izquierdo de esta ecuación son, de hecho, los
coeficientes de rigidez, como se ha definido anteriormente. La matriz del lado derecho
de la ecuación es una matriz unitaria (matriz de identidad). Los elementos de la matriz
en el lado izquierdo de esta ecuación son, de hecho, los coeficientes de rigidez, como
se ha definido anteriormente. Así, la ecuación puede escribirse en la misma forma que
la ecuación 3.17. Inversión de ambos lados de la ecuación 3.17 da la ecuación 3.18. La
sustitución de la ecuación 3.17 en la ecuación 3.20 da:

{ F }=[ k ] { D } (3.21)

Las ecuaciones 3.17 y 3.18 muestran que la matriz de rigidez es la inversa de la matriz
de flexibilidad, y viceversa, siempre que se use el mismo sistema de coordenadas de
fuerzas y desplazamientos en la generación de las dos matrices.
Recordemos que en el método de análisis de fuerza, se introducen liberaciones para
obtener la estructura liberada y el sistema de coordenadas representa la ubicación y
dirección de estas fuerzas eliminadas. En el método de análisis de desplazamiento se
agregan fuerzas restringentes para impedir los desplazamientos de los nudos. En este
caso, el sistema de coordenadas representa la ubicación y dirección de estos

221
 Método de Análisis de Desplazamiento

desplazamientos. Se deduce que los dos sistemas de coordenadas no pueden ser

iguales para la misma estructura. Por consiguiente, la inversa de la matiz de flexibilidad
generada en el método de fuerza es una matriz cuyos elementos son coeficientes de
rigidez, pero no son los que se usan en el método de análisis de desplazamiento.
Además, observamos que la matriz de flexibilidad determinada en el método de fuerza
corresponde a la estructura liberada, mientras que en el método de desplazamiento se
utiliza la matriz de rigidez de la estructura real.
La ecuación 3.19 no de ser confundida con la ecuación: [ f ] { X }={−D0 } usado en el
método de análisis de fuerza en que aplicamos las fuerzas redundantes desconocidas
{ X }de tal magnitud que producen los desplazamientos {−D0 } para corregir las
incongruencias { D 0 } de la estructura liberada.
De manera similar, la ecuación 3.21 no debe confundirse con la ecuación: [ k ] { D }= {−R }
usado en el método de análisis de desplazamiento en la que buscamos un
desplazamiento desconocido { D } de tal magnitud que eliminará las fuerzas artificiales
de restricción { R } .

3.11.2 Comparación de los métodos de fuerza y desplazamiento

Los métodos de fuerza y desplazamiento, estudiados en los capítulos anteriores, son
los principales métodos analíticos en el análisis estructural. Ambos métodos son
ampliamente utilizados, no sólo para el análisis estático, sino también para la
estabilidad y el análisis dinámico. La comparación entre ambos métodos muestra que
tienen tanto semejanzas como diferencias.
Ambos métodos requieren la elección de una estructura básica, sin embargo, el
sistema básico en el método de fuerzas se obtiene liberando fuerzas en la estructura
dada, lo que reduce el grado de indeterminación estática, mientras que en el de
desplazamiento el sistema básico se obtiene introduciendo fuerzas restringentes, lo
que aumenta el grado de indeterminación estática.
Las incógnitas primarias en el método de fuerza, son las fuerzas redundantes X1, X2,
…, Xn y en el de desplazamientos, los desplazamientos en la dirección de los las
fuerzas restringentes añadidas D1, D2,…, Dn

222
 Método de Análisis de Desplazamiento

En ambos métodos se plantean ecuaciones idénticas por su forma, pero diferentes por
su contenido. En el método de fuerzas, las ecuaciones contienen las fuerzas
redundantes X1, …, Xn que destruyen los desplazamientos en la estructura liberada en
la dirección de las redundantes, por cuanto estos desplazamientos no existen en la
estructura dada. En el método de desplazamientos ellas contienen los desplazamientos
D1, …, Dn, que destruyen las reacciones de las restricciones introducidas en el sistema
básico, que no existen en la estructura dada.
En ambos métodos se emplean diagramas unitarios. Sin embargo, en el de fuerzas, los
diagramas unitarios se construyen por las fuerzas y, en el de desplazamientos, por los
desplazamientos.
En ambos métodos, para construir los diagramas definitivos de fuerzas internas se
emplea el principio de superposición.
Las diferencias fundamentales entre estos métodos se presentan en la Tabla 3.1

Tabla 3.1 Diferencias fundamentales entre el método de fuerza y desplazamiento

Criterios de
Método de fuerza Método de desplazamiento
comparación
Estructura primaria o Se obtiene eliminando fuerzas redundantes Se obtiene introduciendo fuerzas
básica en la estructura dada restringentes en la estructura dada.
Las fuerzas (fuerzas o momentos) liberadas Desplazamientos (traslaciones o
Incógnitas iniciales en el sistema básico rotaciones) que han sido impedidos por las
fuerzas restringentes
Igual al grado de indeterminación Igual al grado de indeterminación
Número de incógnitas cinemática
estática
Muchas. La única condición es que la Único. Se obtiene añadiendo las fuerzas
Alternativas para elegir estructura elegida debe ser estable y restringentes para impedir los
el sistema básico estáticamente determinada desplazamientos(traslaciones o rotaciones)
en la estructura dada,
Ecuaciones de compatibilidad Ecuaciones de equilibrio
D01 +f 11 X 1 +f 12 X 2 +…=0 D01 +f 11 X 1 +f 12 X 2 +…=0
Ecuaciones planteadas D01 +f 11 X 1 +f 12 X 2 +…=0 D01 +f 11 X 1 +f 12 X 2 +…=0
.. . . . . . . . . . . . . .. . . . . . . . . . . . .
Número de ecuaciones igual al número de Número de ecuaciones igual al número de
incógnitas incógnitas
Significado de las el desplazamiento en la dirección de las La reacción en la dirección de los
ecuaciones fuerzas liberadas (redundantes) causado desplazamientos eliminados causados por
por la acción las incógnitas primarias la acción de las incógnitas
(fuerzas o momentos) y la carga aplicadas primarias(traslaciones o rotaciones) y las

223
 Método de Análisis de Desplazamiento

es cero cargas aplicadas es cero

[ ] [ ]
f 11 f 12 … f 1n k 11 k 12 … k 1n

Matriz de los [ f ] = f 21 f 22 … f 2n
[ k ]= k 21 k 22 … k2n
coeficientes … … … … … … … …
f n 1 f n2 … f nn kn 1 kn 2 … k nn
Donde [ f ] es la matriz de flexibilidad Donde [ k ] es la matriz de rigidez
f ij es el desplazamiento de la estructura K ij es la fuerza en la estructura primaria,
primaria en la dirección de la redundante que se genera en la dirección de la fuerza
Significado de los Xi debido a una fuerza unitaria en la restringente i cuando se aplica un
desplazamiento unitario en la dirección de
coeficientes dirección de la. redundante X j. Se llaman la fuerza restringente en la dirección j
coeficientes de flexibilidad manteniéndose todos los restantes nulos.
Se llama coeficiente de rigidez

Significado de los
D 0 i representa el desplazamiento en la Ri representan la fuerza de restricción en
términos estructura liberada en la dirección de la la estructura primaria en la dirección de la
independientes redundante Xi debido a la acción de las de la coordenada i debido a la acción de
cargas las cargas
Unidades de los f ij : dimensión de desplazamiento k ij
: dimensión de fuerza (fuerza o
coeficientes (traslación o rotación) momento o rotación)

Un ingeniero experimentado puede elegir como estructura liberada en el método de

fuerza una estructura estáticamente indeterminada, siempre y cuando ésta, se pueda
analizar fácilmente. Igualmente se puede obtener la estructura básica del método de
desplazamiento, introduciendo fuerzas restringentes en menor número que el grado de
indeterminación cinemática, siempre que la estructura resultante resulte fácil de
analizar.
Es el momento de hacer la pregunta: ¿cuándo es más conveniente aplicar el método
de fuerza y cuando el método de desplazamiento? La respuesta general es la
siguiente: Se debe elegir el método de análisis que tenga el menor número de
incógnitas básicas.
Pero considerando que en el método de desplazamiento el sistema básico es único,
que la construcción de los diagramas unitarios y de carga son simples y siempre
estándar y que para una gran cantidad de incógnitas las ecuaciones son incompletas,
las operaciones matemáticas en el método de desplazamiento son considerablemente
más simples que en el de fuerzas. Así en muchos casos, el método de desplazamiento
es más cómodo que el de fuerzas. A pesar de las ventajas indicadas este método no
puede reemplazar íntegramente al de fuerzas por las razones que se exponen a
continuación.

224
 Método de Análisis de Desplazamiento

1. En algunas estructuras, el análisis es más fácil por el método de desplazamiento y

en otras por el método de fuerzas. Esto se puede ilustrar considerando diferentes
estructuras. La Figura 3.49a presenta un pórtico con apoyos empotrados. El pórtico, de
acuerdo con el método de fuerza, tiene nueve incógnitas, mientras que por el método
de desplazamiento solo una: la rotación del nudo rígido. El análisis de esta estructura
por el método de desplazamiento es un problema muy simple. En caso de que esta
estructura se modifique adicionando más elementos que se conecten al nudo rígido, el
número de incógnitas por el método de desplazamiento sigue siendo el mismo,
mientras que el número de incógnitas por el método de la fuerza se incrementa con la
adición de cada nuevo elemento.
En el pórtico mostrado en la figura 3.49b, el número de incógnitas por el método de
fuerza es uno, mientras que por el método de desplazamiento es de seis (cuatro
rotaciones y dos traslaciones). Además, es importante señalar que esta estructura tiene
miembros inclinados, que conducen a cálculos adicionales engorrosos. Por lo tanto, el
análisis de este pórtico es más adecuado por el método de fuerza.
El arco estáticamente indeterminado con extremos empotrados mostrado en la figura
3.49c, se puede analizar por el método de o de desplazamiento. El número de
incógnitas por el método de fuerza es de tres. Con el fin de analizar este tipo de arco
por el método de desplazamiento, su eje curvilíneo debe ser sustituido por un conjunto
de miembros rectos (ya que los elementos estándar del método de desplazamiento son
miembros rectos). Una versión de tal segmentación del arco se muestra en la figura
3.49c con líneas de trazos; en este caso, el número de incógnitas por el método de
desplazamiento es de seis. Obviamente, el método de la fuerza es más preferible aquí.
Si el arco es parte de cualquier estructura compleja, como se muestra en la figura
3.49d, esta estructura puede ser analizada fácilmente por el método de
desplazamiento. De hecho, el número de incógnitas por el método de fuerza es de seis,
mientras que por el método de desplazamiento es sólo uno (la rotación del nudo rígido).
Si se ignoran las deformaciones axiales, En el pórtico de la figura 3.49e, tiene 27
incógnitas si se analiza por el método de fuerza y solo 9 incógnitas con el método de
desplazamiento, por lo tanto, para su análisis debe elegirse el método de
desplazamiento

225
 Método de Análisis de Desplazamiento

(e)

Figura 3.49 Diferentes tipos de estructuras que pueden ser analizadas por el método de fuerza o de
desplazamiento.

2. Si el análisis se efectúa considerando las deformaciones axiales y cortantes de los

miembros, el método de desplazamientos se hace muy complicado, y pierde sus
ventajas, mientras que, las complicaciones en el método de fuerzas, en este caso, son
considerablemente menores. Por ejemplo, la consideración de las deformaciones
axiales y cortantes en el método de desplazamientos para la estructura representada
en la figura 3.49e, las incógnitas aumentan considerablemente: 18 incógnitas
adicionales, (las traslaciones de los nudos), mientras que, en el método de fuerzas, las
incógnitas no aumentan y los inconvenientes solo están ligados con la construcción de
los diagramas de las fuerzas internas y con el cálculo de estos.
3. En el análisis de armaduras, el número de incógnitas del método de desplazamiento
es siempre mayor que el de fuerzas.
4. Los elementos de la estructura básica del método de desplazamiento son
estáticamente indeterminados, y previamente deben analizarse usando el método de
fuerza. Empero, teniendo en cuenta que en método de fuerzas los cálculos son más
complicados que por el de desplazamiento, se recomienda hacer el análisis por el
método de desplazamiento en aquellas estructuras en donde el número de incógnitas
sea algo mayor que por el método de fuerza.

226
 Método de Análisis de Desplazamiento

3.11.3 Propiedades de la matriz de flexibilidad y rigidez

Si una fuerza F i se aplica gradualmente a una estructura (véase figura 3.50), de tal
manera que la energía cinética de su masa sea cero. Sea D i el desplazamiento
resultante en el punto de aplicación y en la dirección de la fuerza F i.

Fi

Figura 3.50 Estructura sometida a la fuerza Fi

Si la estructura se comporta elásticamente, la curva fuerza-desplazamiento sigue la

misma trayectoria durante la carga y la descarga (véase la figura 3.51a). Considérese
ahora que en alguna etapa de la aplicación de la carga se incrementa la fuerza F i en
Δ F i que causa un aumento correspondiente en el desplazamiento Δ Di , el trabajo
realizado por este incremento de carga es

ΔW =F i Δ Di

Como se ilustra en el rectángulo sombreado de la figura 3.51a. En el límite cuando los

incrementos son suficientemente pequeños, se puede escribir

dW =F i d Di

Se deduce que el trabajo externo total realizado por F i durante el desplazamiento D i es

el área bajo la curva fuerza-desplazamiento entre O y Di . Cuando la estructura tiene un
comportamiento elástico-lineal, la curva fuerza desplazamiento es lineal (véase la figura
3.51b), y el trabajo que realiza la fuerza F i es

227
 Método de Análisis de Desplazamiento

1
W = F i Di
2

Fuerza en i Fuerza en i

Desplazamiento en Desplazamiento en
la coordenada i la coordenada i

Figura 3.51. Relaciones fuerza-desplazamiento

Si la estructura se somete a un sistema de fuerzas F 1 , F2 , … , F n que se aplican

gradualmente desde cero hasta su valor final y causan desplazamientos D1 , D2 , … , Dn
en la ubicación y dirección de las fuerzas, el trabajo externo total resulta ser:

n
1 1 1 1
W = F 1 D 1+ F 2 D2 +…+ F n D n= ∑ Fi Di (3.22)
2 2 2 2 1

O en forma matricial

1 T
[ W ] 1 ×1= { F }n× 1 { D }n × 1 (3.23)
2

Donde { F }T es la transpuesta del vector columna { F } que representa las fuerzas. El

trabajo realizado es una cantidad escalar cuyas dimensiones son fuerza x longitud.
Los desplazamientos y las fuerzas se relacionan por la ecuación 3.19, sustituyendo en
la ecuación 3.23 se tiene

228
 Método de Análisis de Desplazamiento

1 T
[ W ] 1 ×1= { F }n× 1 [ f ] n× n { F }n ×1 (3.24)
2

Tomando la transpuesta de ambos miembros el lado izquierdo no cambia. El miembro

de la derecha se convierte en el producto de las transpuestas de las matrices en orden
inverso. Así se obtiene

1 T T
[ W ] 1 ×1= { F }n× 1 [ f ] n× n { F }n ×1 (3.25)
2

De las ecuaciones 3.24 y 3.25 se deduce que la matriz de flexibilidad y su transpuesta

son iguales

[ f ] =[ f ] T

Esto significa que para un coeficiente general de la matriz de flexibilidad

f ij =f ji (3.26)

Es decir, la matriz de flexibilidad es simétrica.

La ecuación 3.17 nos dice que la matriz de flexibilidad es la inversa de la matriz de
rigidez. Puesto que la inversa de una matriz simétrica es otra matriz simétrica se
concluye que la matriz de rigidez [ k ] también es simétrica. Así para un coeficiente de
rigidez cualquiera

k ij=k ji (3.27)

Otra propiedad importante de las matrices de flexibilidad y rigidez es que los elementos
de la diagonal principal f ii o k ii son positivos. El coeficiente f ii es el desplazamiento en
la coordenada i debido a una carga unitaria en i , así el desplazamiento y la fuerza

229
 Método de Análisis de Desplazamiento

actúan en el mismo sentido, por lo tanto f ii es positivo. El coeficiente k ii es la fuerza que

se requiere aplicar en la coordenada i para producir un desplazamiento unitario en i ; de
nuevo, la fuerza y el desplazamiento tiene el mismo sentido, por lo tanto k ii es positivo.
Como se ha visto en los fundamentos de cada método los coeficientes de flexibilidad y
rigidez dependen únicamente del tipo de estructura, pero no dependen de cargas
externas, asentamientos de apoyos, cambios de temperatura o errores de fabricación.
Los determinantes de la matriz de rigidez o flexibilidad deben ser positivos. Esta
condición describe una propiedad interna de estructuras basadas en una
ley fundamental de los sistemas elásticos: la energía potencial de una estructura
sometida a cualquier carga es positiva. Puesto que D ≠ 0, la solución a las ecuaciones
de compatibilidad o de equilibrio de cualquier estructura sometida a cualquier tipo de
carga, cambio de temperatura o asentamiento de apoyos, etc., es única.

En resumen, podemos enumerar las propiedades de las matrices de flexibilidad y

rigidez:

1. La matriz de flexibilidad y rigidez son matrices cuadradas simétricas ( f ij =f ji , k ij =k ji ).

2. Los términos de la diagonal principal de la matriz de flexibilidad y rigidez siempre
mayores que cero
3. La matriz de flexibilidad y rigidez son definida positiva (determinante diferente de
cero), es decir, su inversa existe.
4. La matriz de flexibilidad es la inversa de la matriz de rigidez, y viceversa, siempre
que se use el mismo sistema de coordenadas de fuerzas y desplazamientos en la
generación de las dos matrices.
5. Las matrices de flexibilidad y de rigidez no dependen de la solicitación externa
(cargas, asentamientos de apoyo, variaciones de temperatura). Las matrices de
flexibilidad y rigidez dependen de las propiedades de la estructura: geometría, material,
condiciones de apoyo y de las propiedades de los elementos

230