Capítulo B-7

Método de las
Matrices de Rigidez.
Ecuaciones de Equilibrio.

ft)
ra
B-7.1 Introducción.
Como ya se ha estudiado anteriormente, el método del equilibrio se basa en exigir las condi-
D
ciones de equilibrio en ciertos puntos de la estructura, quedando como incógnitas los desplaza-
mientos generalizados (traslaciones, giros, etc.). Por este motivo también recibe el nombre de
método de los desplazamientos (ver Figura B7.1).
r(
n n
Pi   K i j  j   Pi j (B-7.1)
j 1 j 1
lle

Siendo K i j la rigidez en el punto i cuando se produce un desplazamiento en el punto

j . Esta rigidez se suele definir como la fuerza necesaria Pi j que hay que hacer en i (para ga-
O

rantizar su condición de compatibilidad i  0 ), cuando se aplica un desplazamiento unitario en

el punto j (  j  1 ).

P1
S.

Pi δ1

i P2
δ2

δj

Pj
Figura B7.1 – Sólido continuo en equilibrio – Acciones y movimientos.
 S. Oller Capítulo B-7 B7-3

Como puede verse en la ecuación (B-7.1), es una expresión que contiene una suma fuer-
zas en función de los desplazamientos.
Presentada la idea fundamental en la que se basa el método del equilibrio, a continua-
ción se escribirán las ecuaciones de equilibrio para una barra estructural en función de los
movimientos de sus extremos. Es importante tener en cuenta que esta formulación se basa
en:
1. Ecuaciones de la elástica de Euler-Bernoulli, válidas para secciones transversales que
mantienen la planaridad y ortogonalidad al eje de la barra, antes y después de la
deformación.
2. Los desplazamientos trasnacionales y angulares son pequeños en relación al tama-
ño de la barra.
3. Los materiales son elásticos e isótropos.

ft)
B-7.2 Relación entre las ecuaciones de compa-

ra
tibilidad y equilibrio para una barra simple.
Ecuaciones de la elástica de una barra en
D
función de sus movimientos.
r(
En este apartado se deducen las ecuaciones de la elástica de una barra y se muestra la duali-
dad que hay entre la formulación del método de la compatibilidad y el método del equilibrio. Para
ello se aísla una barra perteneciente a una estructura hipotética, sometida a un estado de
cargas exteriores, más una acciones que actúan en sus extremos y que provienen de la inter-
lle

acción con el resto de la estructura.

Primeramente se aplica el principio de superposición y el concepto de compatibilidad
de desplazamientos en los extremos de la barra. De aquí se obtiene la formulación de las
O

ecuaciones elásticas que compatibilizan los giros y desplazamientos en los extre-

mos de barra.
En un segundo paso, se aplica el principio de superposición y se introduce el concepto
S.

de equilibrio de esfuerzos de extremo de barra y se deducen la ecuaciones elásticas que

equilibran los esfuerzos de momento, cortante y axil en los extremos de la barra.
A continuación se presenta una barra perteneciente a una estructura cualquiera, a partir
de la cual se desarrollarán las ecuaciones elásticas de compatibilidad y equilibrio.

Figura B7.2 – Movimientos y esfuerzos de extremo en una barra perteneciente a una estructura.
B7-4 S. Oller Método de las Matrices de Rigidez

En la Figura B7.3 se muestra la descomposición del comportamiento estructural de la ba-

rra de la Figura B7.2 con el objetivo de obtener: A) una formulación compatible cuyas
incógnitas sean los esfuerzos de extremo de la barra y b) una formulación del equili-
brio cuyas incógnitas sean los movimientos de extremo de la barra.

A) Ecuación de la elástica compatible con los B) Ecuación de la elástica equilibrada con los esfuerzos
desplazamientos de extremo de barra de extremo de barra

Se liberan los grados de libertad de los extremos Se restringen todos los grados de libertad de los extremos
de barra → Estructura estáticamente determi- de barra → Estructura cinemáticamente determinada.
nada. Incógnitas los desplazamientos de ex- Incógnitas los esfuerzos de extremo.
tremo.

ft)
ra
D
r(
lle
O
S.

6 EI
12 EI M iij  
L Q  2 
i
ij



Figura B7.3 – Descomposición del comportamiento de una barra: A) Sometida acciones y mo-
mentos externos (Método de la compatibilidad), B) Sometida a movimientos de extremo (Método
del equilibrio).

A partir de compatibilizar los movimientos que se muestran en la Figura B7.3 (A) y de

equilibrar los esfuerzos que se muestran en la Figura B7.3 (B), resultan las siguientes ecua-
ciones de la elástica en dos campos (ver Figura B7.4),
 S. Oller Capítulo B-7 B7-5

A) Ecuaciones de compatibilidad. Incógnita esfuerzos de B) Ecuaciones de equilibrio. Incógnita desplazamien-

extremo de barra → Matriz de Flexibilidad. tos de extremo de barra → Matriz de Rigidez.

Compatibilidad de movimientos en extremos: Equilibrio de cortantes en extremo:

"i ": i         "i ": Qi  Qi  Qi  Qi  Qi
q i j q i j ij
i i i
 
"j ":  j   j   j   j   "j ": Q j  Q j  Q j  Q j  Q j
q i j q i j ij

Sustituyendo por sus expresiones:

Sustituyendo por sus expresiones,
 6 EI 6 EI 12 EI
Qi  Qi   2 i   2  j  3 (v j  vi )
q

 M i  M j   v j  vi 
 i  i 
q
   
 3 EI 6 EI    Q  Q q  6 EI   6 EI   12 EI (v  v )
  j j
2
i
2
j
3
j i
   q  M i   M j    v j  vi 
 j j  
 6 EI 3 EI    Equilibrio de momentos en extremo:

ft)
"i ": M i  M i q  M i i  M i j  M i ij

"j ": M j  M j  M j  M j  M j
q i j ij
Y en forma matricial resulta la siguiente forma para la ecua-

ra
ción de compatibilidad, Sustituyendo por sus expresiones:

 q  v j  vi  4 EI 2 EI 6 EI
     M i  M i   i    j   2 (v j  vi )
q

 i     
i       3 EI 6 EI  M  
 
 j   q  v j  vi

  
D 
  M
i


M  M q  2 EI   4 EI   6 EI (v  v )
j
 j j i j
2
j i
 j       6 EI  
   3 EI 

r(
Y en forma matricial resulta la siguiente forma para la
ecuación de compatibilidad,
Ec. de compatibilidad
(B-7.2)  12 EI 6 EI 12 EI 6 EI 
  
U  U q  F  F  3 2 
lle

2 3
Qi  Qi   6 EI 
q
4 EI 6 EI 2 EI  vi 
   q    
M i  M i    2  2   i 
    q    v 
Q j  Q j    12 EI 6 EI 12 EI 6 EI   j 
Siendo U los desplazamientos de extremo conocidos, M j  M q    3
   j   2 3  2   j 
O

6 EI 2 EI 6 EI 4 EI 
U q aquellos debido a las cargas exteriores más desplaza-  2 
   2  
mientos de extremo (ambos conocidos), F la matriz de
flexibilidad y F los esfuerzos incógnitas. Ec. de equilibrio
(B-7.3)
S.

Solución: F  Fq  K  U
1 q 
FF
( U  U ) q
Siendo F los esfuerzos de extremo conocidos, F los
K esfuerzos de extremo conocidos de fijación perfecta
debido a las cargas exteriores, U los desplazamientos
de extremo incógnitas, K la matriz de rigidez.
1
 4 EI 2 EI     
 3 EI  Solución:
    F 6 EI 
K 
1
  UK 1
 (F  F )
q

 2 EI 4 EI     
      6 EI 3 EI 
F


Figura B7.4 – Formulación de las ecuaciones del método de: A) compatibilidad, B) equilibrio.
B7-6 S. Oller Método de las Matrices de Rigidez

En forma independiente y desacoplada de las ecuaciones anteriores, se obtienen las

ecuaciones de equilibrio axial de la barra, las cuales completan la ecuación de equilibrio B-
7.3 anteriormente obtenida.

ft)
ra
Figura B7.5 – Descomposición del comportamiento de una barra sometida a desplazamientos
axiales en sus extremos (Método del equilibrio).
D
Resultando así las siguientes ecuaciones de equilibrio de axiles en los extremos,

"i ": N i  N i p  N i i  N i j
r(

"j ": N j  N j  N j  N j
p i j

Sustituyendo por sus expresiones:

lle

 EA EA
 N i  N i   ui   u j
p


 N  N p  EA u  EI u
 j
O

j i j
 
Insertando y reordenando estas expresiones en la ecuación B-7.3, resulta la siguiente
ecuación de equilibrio para la barra i  j , que en adelante llamaremos ecuación de equili-
brio del elemento i  j o simplemente ecuación de equilibrio elemental e .
S.

e
 EA EA 
  0 0  0 0 

e
 
e
N i q   0 12 EI 6 EI 12 EI 6 EI  e
N i   0   2 ui 
   q   3 2 3    
Qi  Qi   6 EI 4 EI 6 EI 2 EI  vi 
M i  M q   0  2 0 i 
 i     2  
N   q    
u
 j  N j    EA 0 0
EA
0 0   j
Q j  Q q      v j 
   j
  12 EI 6 EI 12 EI 6 EI   
M j  M q   0  3 0   j 
 j   2 3 2 

 6 EI 2 EI 6 EI 4 EI 
 0  2 0
  2  
 S. Oller Capítulo B-7 B7-7

Que en forma compacta queda representada por la siguiente ecuación,

F e  (F q ) e  K e  U e (B-7.4)

Tal que e representa el elemento eesimo de barra de la estructura.

B-7.3 Formalización del método de las matrices

de rigidez.
A continuación se presenta la formalización del método para ser utilizado en estructuras
de barras. Esta formulación permitirá su programación en forma sencilla utilizando una
hoja de cálculo. Para ello, al final de este capítulo se detallará dicha programación para un

ft)
ejemplo en concreto.

B-7.3.1 Equilibrio del elemento de barra e en ejes locales.

ra
El equilibrio de cada elemento de barra se expresará en forma compacta expresando las
ecuaciones agrupadas en sub-matrices nodales. Así, la ecuación (B-7.4) se escribe en sub-
matrices agrupadas por nodo,
D
K ije   U ie   Fi  
 q e
Fie   K iie
r(
 e  e     (B-7.5)
F j  K ji  
K ejj  U ej   F q e 
 i 
lle

Siendo en esta,
e
 uie   N ie   N iq 
O

 e  e  q
 vi   Qi  e
 Qi 
 U   e 
e
F   M e 
e
F   M e 
q
U e   ei    ei  ; F e   ie    ei  ; F 
q e
  iq    qi 
U j   u j  F j   N j  Fi   N j 
S.

 v ej   Q ej   Q qj 
 e  e  q
 j   M j   M j 

Para tener en cuenta las distintas orientaciones de la barra en el plano, se define el siguiente
cambio de base,
 B7-8 S. Oller Método de las Matrices de Rigidez

y x j  xi
j cos  
y-j ( x j  xi ) 2  ( y j  y i ) 2
y j  yi
e sen 
( x j  xi ) 2  ( y j  y i ) 2

i α
y-i
 x e   cos  sen 0  x 
 e  
x x e  Tk  x   y    sen cos  0  y  (B-7.6)
 e   0 0 1 
  
x-i x-j

ft)
B-7.3.2 Definición del equilibrio elemental en ejes globales.

ra
Realizando el cambio de base del “equilibrio local elemental”, ecuación (B-7.5), se ob-
tiene el “equilibrio en ejes globales” D
F e  K e U e  K e U e   F q  e
 i ii i ij j i

F je  K eji U ie  K ejj U ej   F jq 
e
r(

Sustituyendo estas expresiones por las correspondientes en base global, resulta

 Te  Fi   K iie  Te  U i   K ije  Te  U j    Te  Fiq 

lle


 e
 T  F j   K ji  T  U i   K jj  T  U j    T  F j 
e e e e e q
O

xe
y M ; e e
Q ej ; v ej j j

N ej ; u ej
S.

yj j
e
y
e
M i e ; ie
αe
yi i
N ie ; uie
Qie ; vie
x
xi xj

Figura B7.6 – Representación del elemento de barra “e” en ejes locales.

 S. Oller Capítulo B-7 B7-9

Reordenado los términos, queda la ecuación de equilibrio local expresada en ejes globa-
les,

Fi   TeT  K iie  Te  U i   TeT  K ije  Te  U j  Fiq



F j   T  K ji  T  U i   T  K jj  T  U j  F j
eT e e eT e e q

Y en forma matricial, se escribe.

Fi    T  K ii  T  T  K ije  Te    U i  F q 
 eT e e eT

        iq 
F j   TeT  K eji  Te  T eT
 K jj  T   U j  Fi 
e e

o también,

ft)
T
Fi  Te 0   K iie K ije  Te 0   U i  Fiq 
        
F j   0 Te  K eji K ejj   0 Te  U j  Fiq 
     

ra
F Π eT Ke Πe U Fq

F   Π eT  K e  Π e   U  F q
D
F  K  U  Fq (B-7.7)
r(
Tal que estas matrices elementales definidas en referencia global, se definen como,
lle

 cos  e sin  e 0 0 0 0
 
  sin  cos  e
e
0 0 0 0
Te 0  0 0
O

0 1 0 0
Πe    
0 Te   0 0 0 cos  e sin  e 0
 0 0 0  sin  e cos  e 0
 
 0 0 0 0 0 1 
S.

 ui   Fxi   Fxiq 
     q
 vi   Fyi   Fyi 
    M   M iq 
U  Π eT  U e   i  ; F  Π eT  F e   i  ; F q  Π eT   F q e  q
u j   Fxj   Fxj 
v j   Fyj   Fyjq 
     q
 j   M j   M j 
B7- S. Oller Método de las Matrices de Rigidez
10

Fyj ; v j xe
y M j ;j

yj Fxj ; u j
j
ye q
Fyi ; vi
e

M i ; i αe
yi i
Fxi ; ui

x
xi xj

ft)
Figura B7.7 – Representación del elemento de barra “e” en ejes globales.

ra
B-7.3.3 Colocación y ensamblaje.

Una vez se ha definido el equilibrio de cada elemento “e” de barra en ejes globales, se
D
debe “ensamblar” la contribución de cada una de ellas dentro de la estructura global. A
continuación se explica este concepto siguiendo un modelo estructural simple (ver detalle
r(
del procedimiento numérico en el archivo adjunto “Ejemplo Estructura-1- Mat de Rigi-
dez.xls”)
Dada la estructura de la figura, obtener a través del “Método de las Matrices de Rigidez”
lle

las leyes de Momento Flector, Esfuerzo Cortante y Esfuerzo Axil y Desplazamientos en

todos los nudos.
O

y
Datos:
 E  2.10 10 5 MN / m 2
q=10KN/m  1 2 2 2
 A  5 10 m ; A  6 10 m ; A  8 10 m
2 2 2 3 2
P=50KN
S.

 I 1  8 10  4 m 4 ; I 2  1,5 10 3 m 4 ; I 3  2 10 3 m 4

2 3 
2

5,0m
1 3

1 4 x

3,0m 8,0m

Figura B7.8 – Datos mecánicos y geométricos de la estructura a estudiar.

 S. Oller Capítulo B-7 B7-11

Tabla de Relaciones Nodales

Barra N-i N-j x-i y-i x-j y-j cos(α-i) sen(α-i)

1 1 2 0,0 0,0 0,0 5,0 0,0x100 1,0x100
2 2 3 0,0 5,0 3,0 5,0 1,0x100 0,0x100
3 3 4 3,0 5,0 11,0 0,0 8,48x10-1 -5,3x10-1

Una vez definida las matrices elementales en base local (B-7.4) se introducen los corres-
pondientes cambios de base para expresar el equilibrio elemental en base global (B-7.7).
Una vez realizado esto, se procede a colocar y ensamblar (suma selectiva) cada matriz con
las restantes, permitiendo así considerar el equilibrio de cada una de ellas dentro de la es-
tructura global.

 Π F  
n n

FG  AΠTe  K e  Π e   UG  eA1 T q e

ft)
e
e 1
  
KG
F 
G
q
(B-7.8)

ra
FG  K G  UG   F q 
G

Siendo la rigidez global ensamblada,

D
r(
n
  Te 0 T  K e
n
K ije   Te 0  
K  A Π  K  Π   A  
G eT e e
e
  eii   
e 1 e 1  0
 T  K ji K ejj   0 Te  


   T  K ii  T   T  K ij  T   
n 
lle

eT e e eT e e

 A 
e 1  T  K eji  Te   TeT  K ejj  Te   
eT
  
O

1 2 3 4 NUDOS
 T1T K11
1
T1 T1T K12
1
T1 0 0  1
S.

 T 1 
T K T T1T K122T1  T2T K 222 T2 T2T K 223T2 0 2
K G   1 21 1 
 0 T2T K 32
2
T2 T2T K 33
2
T2  T3T K 333T3 T3 K 34 T3 
T 3
3
 T 
 0 0 T K T
3
3
43 3 T3T K 344 T3  4
B7- S. Oller Método de las Matrices de Rigidez
12

Siguiendo el mismo procedimiento, se obtiene la fuerza global ensamblada

 e
0   Fi  
 q e    i
 TeT  F q e

Π 
n n T n
 T 
F 
q G
 A eT
 F q e
 A      A 
e 1 e 1  0


Te   F q e 
 j 
e 1


 TeT   F jq 
e

NUDOS
  Fx1 q   0 
   0  1
  Fy1    
q

   0 
  M q   0 0 
 1 q   
 Fx 2    qL2 

ft)
   0   2
 T1T  F1q 
1

  y2 

q
  2 
 F
 T q 1 q 2
   qL22 
T1  F2   T2  F2    M 2   
T 0 
Π 
q

ra
n
12 
F 
q G
 Fq 
e
 A eT
       
T2T  F3q   T3T  F3q     Fx 3    0  0 
2 3 q
e 1

   q  qL2 
 T3  F4 
q 3    F     0  3

T
D  
y 3
  2 
 M 3    qL22 
q

 q  
r(
 Fx 4    12 
 q  0 
 Fy 4    0  4
 q  
 M 4    0  ;
lle

NUDOS
O

 Fx1   H 1 
F   V  1
 y1   1 
 M 1  M 1 
S.

   
 Fx 2   P 
 T1T F11   Fy 2   0 
   
  2
n  T 1
F  AΠ  F    T 2
G eT e
T 2 
T1 F2  T2 F2  M 2   0 
     ;
   
3 
e 1
T2 F3  T3 F3   Fx 3   0 
T
   
 T3T F43   Fy 3   0 
   
  3

M 3   0 
 F  H 
 x4   4  4
 F y 4   V4 
M   0 
 4  
 S. Oller Capítulo B-7 B7-13

Y el desplazamiento en coordenadas globales se expresa

NUDOS
 u1 
v  1
 1
 1 
 
 u2 
U1G   v2  2
 G  
 U   
U   G2    2 
G

U 3   u3 
U G4   v3  3
 
3 

ft)
u 
 4
 v4  4
 

ra
 4
D
B-7.3.4 Solución y obtención de los esfuerzos en los extremos de la ba-
r(
rra e.

La solución de la ecuación (B-7.8), conduce a obtener los desplazamientos de extremo

de barra,
lle

FG  K G  UG   F q  U G   K G    F G  F q  
G 1 G
 (B-7.9)
 
O

Como ya se ha visto, UG representa los desplazamientos generalizados incógnitas, K G

la matriz de rigidez global (de toda la estructura), F G las fuerzas nodales externas y
F 
G
S.

q
las fuerzas nodales de fijación, que resultan de restringir completamente los movi-
mientos de extremo de la barra (reacciones vinculares).
Conocidos los desplazamientos globales, se debe obtener a continuación los esfuerzos
interiores en la estructura para cada elemento de barra en ejes locales. Para esto hay que
utilizar la ecuación (B-7.5), luego de cambiar de base los desplazamientos globales a locales,

 U e  Te 0   U i 
U e = Π e  U   ei      (B-7.10)
U j   0 Te  U j 
B7- S. Oller Método de las Matrices de Rigidez
14

Siendo U el subconjunto de los desplazamientos globales U G , que le corresponden al

elemento e. Así, por ejemplo, para la barra e=2 (i=2 , j=3 )

 u1 
v  scattering
 1
 1 
 
 u2   u2 
 
U1   v2 
G
 v2 
 G  
 U 2  2   
U  G  
G
; U   2
U 3   u3   u3 
U G4   v3   v3 
   

ft)
3  3 
u 
 4
 v4 

ra
 
 4

Cambio de base
D
u2e  2   cos  e  2 sin  e  2 0 0 0 0   u2 
 e2    
r(
e2 e2
 v2    sin  cos  0 0 0 0   v2 
 U     0  2 
e2 e2
0 0 1 0 0
U e2   ie22    e2 2     
U j 3  u3   0 0 0 cos  e  2 sin  e  2 0   u3 
lle

 v3  
e  2
0 0 0  sin  e  2 cos  e  2 0   v3 
 e2    
3   0 0 0 0 0 1  3 
O

Sustituyendo estos desplazamientos en la ecuación de equilibrio local (B-7.4)

F  (F q )e  K e  U e , resultan los esfuerzos de extremo de barra en ejes locales,
e
S.

 N ie 
 e
 Qi 
F   M e 
e

F e   ie    ei 
F j   N j 
 Q ej 
 e
 M j 

Siguiendo este procedimiento se obtienen los esfuerzos de todas las barras de la estructura.
 S. Oller Capítulo B-7 B7-15

B-7.3.1 Solución Numérica de la Ecuación de Equilibrio Global

Dada la ecuación de Equilibrio Global y su Solución

FG  K G  UG   F q  U G   K G    F G  F q  
G 1 G
 (B-7.11)
 
Hay que tener en cuenta que se trata de un sistema de ecuaciones lineal particular que
necesita un tratamiento especial para su solución. Sus características son:

1. Es un sistema de ecuaciones de (m´ m) , siendo m = n ´3 el Número de Grados

de Libertad y n el Número de Nodos de la Estructura,
2. El Vector Desplazamiento Global tiene grados de libertad cuyos movimientos
son conocidos (Condición de Contorno del Problema) y otros movimientos que
son incógnitas. Por otro lado, el Vector de Fuerzas Global tiene algunas fuerzas

ft)
incógnitas (Reacciones Vinculares donde los Desplazamientos son Conocidos) y
otras fuerzas conocidas donde los movimientos son incógnitas:

ra
ü
ï
 u1   Fx1  ï
ï
ü
ï F  ï
v  ï
ï ï Incónitas
ý
 1 ý Datos  y1  ï
ï  M1  ï
1 
 
ï
ï
ï
þ
D  
ï
ï
ï
þ
   ü   ü
ï
ï
ï  Fx 2  ï
ï
r(
 u2  ï
ï ï
  ï   ï
ï
ï
ï  Fy 2  ï
 v2  ï ï
ï
ï
U1  2 
G ï
ï F1   M 2 
G ï
ï
ï
ï  G  
lle

 G   ï ï
ï
U 2    
ï
ï
ï F2      ï
ï
ï
    ï     ï
UG  U 3G    u3  ï
ï ; F  F3G    Fx 3 
G
ï
ï
ï
ï    F  ï
   v  ï ï
ï
O

   3 ï
ï    y3  ï
ï
ï ý Datos
U Gn  3  ý Incógnitas
ï FnG   M 3  ï
ï   ï
ï
  ï ï
   ï
ï
ï
   ï
ï
ï    ï
  ï
S.

ï
ï ï
  ï   ï
ï
ï    ï
   ï
ï ï
ï
 un  ï
ï F  ï
ï ï
ï
  ï
ï  xn  ï
ï
 vn  ï
ï  Fyn  ï
ï
  ï
ï M  ï
ï ï
ï
 n ï
ï
þ
 n ï
þ

3. En este sistema de ecuaciones particular, ocurre que cuando se conocen los despla-
zamientos (condición de contorno), se desconocen las fuerzas (reacciones):
B7- S. Oller Método de las Matrices de Rigidez
16

U iG conocido  Fi G incógnita
U Gj incógnita  Fi G iconocido

4. Esta particularidad obliga a escribir el sistema de ecuaciones agrupado en dos sub-

conjuntos, los Grados Liberes “L” y los “Grados Restringidos “R”

 F q G 
FLG   K GLL K G
 G
UL   L 
 G G  G   q G 
LR
(B-7.12)
FR  K RL K  U R   F  
G
RR
 R 
Son incógnitas los desplazamientos de los nudos libres U GL y son incógnitas las
fuerzas de los nudos restringidos FRG

ft)
Así la solución se obtiene por condensación de variables, método que se explica-
rá a continuación, o también mediante el método de eliminación de filas y co-
lumnas de los grados conocidos

ra
B-7.3.1.1 Método de Condensación
Reescribiendo la ecuación (B-7.13) en forma canónica y luego formando dos grupos de
D
ecuaciones con sus respectivas incógnitas (L para los nudos libres cuyas incógnitas son
UGR y R para los nudos restringidos cuyas incógnitas son FRG ), se tiene
r(
 L F G  K G U G  K G U G   F q G
 L LL L LR R L
 (B-7.13)
 R FRG  K GRL U GL  K GRR U GR   FRq 
G
lle

Se despeja la incógnita U GL de la primera, resolviendo las incógnitas de desplazamiento, y se

O

la sustituye en la segunda ecuación, resolviendo las incógnitas de fuerza:

Fuerzas externas nodales

U GL   K GLL   F G  K G U G   F q G 
1
(B-7.14)
 L LR R L

S.

Fuerzas nodales de bloqueo de

desplazamientos de la barra


FRG  K GRL  K GLL 
1

 L LR R L
 RR
 F G  K G U G   F q  G   K G U G   F q G
R R (B-7.15)

De esta forma queda resuelto el sistema de ecuaciones

B-7.3.1.2 Método de Eliminación de Filas y Columnas

A continuación se comenta el método de eliminación, que es el que se aconseja utilizar
en los programas computacionales de resolución de ecuaciones estructurales.
 S. Oller Capítulo B-7 B7-17

El método de eliminación de filas y columnas es un caso particular del método de con-

densación y consiste en eliminar las filas y columnas de los grados de libertad cuyos despla-
zamientos son conocidos. Luego se reordena el “sistema de ecuaciones reducido” que
resulta de esta eliminación y finalmente se resuelve este sistema de ecuaciones.

G   G G
 Fx1     u1  0    Fx1  q 
F    K11   K12G   K1Gn   v  0   
G

   Fy1  
q
 y 1     1
 M1    1  0   
  M1  
q
    
       
          
        F q 
 Fx 2  u
 2   x2 
 
F    K 21G   K 22G    K 2Gn    v   Fy 2 q 
     
y2 2
  q

ft)
M2      2    M 2  

    G      

  
 

K Reducida


    
     

ra
       
         
 Fx n    un   q 

  Fx n  
    
    
 Fy n 
Mn 
  K n1 
G
 K nG2 
D   K nnG    vn  
 n   Fy n  
q

 
       q
     M n  
r(
 

G
Quedando así un sistema de ecuaciones reducido K Reducida :
lle

G
  Fx 2 q 
   u2   
G G
 Fx 2 
 y 2  
q
    v   F
 K   K 
O

G G

 Fy 2   22 2n   2  
   2   M 2  
q
M2 
        
        
  G 
  
K      
S.

            
   Reducida        
       
 un   Fx n  
q
 Fx n   
F    v   
 K n 2   K nn 
G G
 yn     n   F q 
 M n   n   
yn
 
   M q
 n  

U GReducida   K GReducida    FReducida  F q  

1 G G

 (B-7.16)
 Reducida 
 B7- S. Oller Método de las Matrices de Rigidez
18

Siendo este el sistema de ecuaciones que se debe resolver. Sustituyendo luego esta solución
reducida en el sistema de ecuaciones completo se obtiene el vector fe fuerzas nodales F G a
nivel global que permite verificar el equilibrio en cada nodo y las reacciones vinculares

FG  K G  UG   F q 
G

Siendo
FReducida  U Reducida 
G G

FG   G  ; U G
  G 
FReacciones  U Reacciones 
(B-7.17)
G
FReducida : Fuerzas nodales Prescritas (fuerzas concentradas en nodos)
G
U : Desplazamientos nodales

ft)
Reducida
G
FReacciones : Reacciones vinculares
U GReacciones : Desplazamientos nodales prescritos o impuestos

B-7.3.2 Solución de la estructura propuesta.

ra
D
Reacciones y Leyes de Esfuerzo (Ver hoja Excel)
r(
y
lle

q=10KN/m
P=50KN
3
O

2 1

1
S.

M 1  0.011667 MNm
1 4 x

F1x  H 1  0.00522 MN
F1 y  V1  0.004242 MN
F4 x  H 4  0.04477 MN F4 y  V 4  0.025757 MN
 S. Oller Capítulo B-7 B7-19

1.784×10-2
2.175m

0.424m
1.441  10 2

1.441×10-2

1.534×10-2
s M
s2
M ( s )  1.444  10 3  4.24  10  3 s  0.01
2

1.167  10 2

ft)
1.891×10-3
4.24×10-3

ra Q
2.576×10-2

s
D
Q( s )  4.24  10 3  0.01 s
r(

5.221×10-3
lle
O 4.47×10-2
S.

5.162×10-2

4.24×10-3
B7- S. Oller Método de las Matrices de Rigidez
20

B-7.4 Ejemplos simples.

A continuación se presentan dos ejemplos simples, que puedan resolverse a mano, del
método de las matrices de rigidez.

B-7.4.1 Ejemplo 1.

ft)
ra
D
r(
lle
O
S.
 S. Oller Capítulo B-7 B7-21

ft)
ra
D
r(
lle
O
S.
B7- S. Oller Método de las Matrices de Rigidez
22

ft)
ra
D
r(
lle
O
S.
 S. Oller Capítulo B-7 B7-
23

B-7.4.1 Ejemplo 2.

ft)
ra
D
r(
lle
O
S.
B7- S. Oller Método de las Matrices de Rigidez
24

ft)
ra
D
r(
lle
O
S.