ANALISIS MATRICIAL DE ESTRUCTURAS POR

EL METODO DE LA RIGIDEZ
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 1

Introduccin
Los mtodos clsicos de anlisis estructural desarrollados a fines del siglo XIX, t
ienen las
cualidades de la generalidad, simplicidad lgica y elegancia matemtica. Desgraciad
amente,
conducan a menudo a clculos muy laboriosos cuando se los aplicaba en casos prc
ticos, y
en aquella poca, esto era un gran defecto.
Por esta razn sucesivas generaciones de ingenieros se dedicaron a tratar de
reducir el
conjunto de clculos. Muchas tcnicas ingeniosas de gran valor prctico fueron apa
reciendo
(Mtodo de Cross), pero la mayora de las mismas eran aplicable slo a determinados t
ipos de
estructuras.
La principal objecin a los primeros mtodos de anlisis fue que los mismos co
nducan a
sistemas con un gran nmero de ecuaciones lineales, difciles de resolver manualmente.
Con los computadores, capaces de realizar el trabajo numrico, esta objecin no tien
e ahora
sentido, mientras que la generalidad de los mtodos permanece. Esto explica po
r qu los
mtodos matriciales deben en su tratamiento bsico de las estructuras ms al siglo XIX
que al
XX.
El empleo de la notacin matricial presenta dos ventajas en el clculo de estructuras. D
esde el
punto de vista terico , permite utilizar mtodos de clculo en forma compacta, pre
cisa y, al
mismo tiempo, completamente general. Esto facilita el tratamiento de la teora de est
ructuras
como unidad, sin que los principios fundamentales se vean oscurecidos por operaci
ones de
clculo, por un lado, o diferencias fsicas entre estructuras, por otro.
Desde el punto de vista prctico, proporciona un sistema apropiado de anlisis de estr
ucturas
y determina una base muy conveniente para el desarrollo de programas de computacin
.

En contraste con estas ventajas, debe admitirse que los mtodos matriciales se car
acterizan
por una gran cantidad de clculo sistemtico
Las virtudes del clculo con computadora radican en la eliminacin del la preocupacin
por las
operaciones rutinarias, el ingenio necesario para preparar el modelo con que se
pretende
representar la realidad y el anlisis crtico de los resultados.
Se debe ser consciente que sin un modelo adecuado o sin una interpretaci
n final, el
refinamiento en el anlisis carece de sentido.

Mtodo de la Rigidez
Hiptesis: Estructura lineal- Todos los movimientos y esfuerzos son funciones lineales
de las
cargas- Pequeas deformaciones (ecuaciones de equilibrio en la estructura no distorsion
ada).
Las barras son rectas y de seccin constante.
Para estudiar una estructura por el mtodo de la rigidez, al igual que en cu
alquier otro
problema elstico, disponemos de tres conjuntos de ecuaciones que deben cumplirse.
Ecuaciones de compatibilidad
Ecuaciones constitutivas
Ecuaciones de equilibrio
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 2

Las ecuaciones de compatibilidad relacionan las deformaciones de barra

s con los
desplazamientos nodales. Introduciendo estas relaciones en las ecuaciones con
stitutivas,
relacionamos las fuerzas en los extremos de barras con los desplazamientos nodales
Introduciendo estas ltimas relaciones en las ecuaciones de equilibrio se obtiene un c
onjunto
de ecuaciones de fuerzas nodales en funcin de desplazamientos nodales, que pu
eden ser
consideradas como Ecuaciones de Equilibrio de la estructura en funcin de desplazamie
ntos.
La resolucin de este sistema de ecuaciones nos permite obtener el valor de las i
ncgnitas
(desplazamientos nodales), a partir de los cuales se obtienen las solicitaciones de la
s barras
de la estructura, as como las reacciones.

Cuando se van a calcular las relaciones esfuerzos de extremo de barra - desplazamie

ntos, es
natural escoger un sistema de coordenadas que haga estas ecuaciones lo ms
sencillas
posible.
Tomaremos por lo tanto como eje x el que coincide con el eje geomtrico de la pieza y
los ejes
y y z coincidentes con los ejes principales de la seccin transversal.
Tal sistema pertenece a la barra, y no depende de la orientacin de la misma en la estr
uctura y
lo denominaremos sistemas de ejes locales.
Por el contrario, cuando las piezas se unen entre s para formar la estructura, es
necesario
tener un sistema de coordenadas comn para todos los movimientos y esfuerzos de
extremo
de barras para poder aplicar las condiciones de equilibrio y compatibilidad. A dicho sist
ema lo
denominaremos sistema de ejes globales.
Dichos esfuerzos de extremos de barras y desplazamientos dependern del tipo de est
ructura
que estamos resolviendo, para barras de:
a) Reticulado Plano: tendremos dos desplazamientos por nudo
b) Reticulado Espacial: tres desplazamientos por nudo.
En ambos casos slo tendremos esfuerzos normales.
c) Prtico Plano: tres desplazamientos por nudo. (una rotacin en el plano del prti
co y dos
traslaciones), como solicitaciones de extremo de barra una fuerza axial, un esfuerzo de
corte y
un momento flector.
d) Prtico Espacial: seis desplazamientos por nudo, tres traslaciones y tres rotacione
s. como
solicitaciones de extremo de barra una fuerza axial, dos esfuerzos de corte dos
momentos
flectores y un momento torsor.
e) Emparrillado de vigas: tres desplazamientos nodales (un corrimiento normal al pla
no de la
grilla) y dos rotaciones alrededor de los ejes contenidos en el plano mencio
nado). Los
esfuerzos son un cortante y dos momentos (un torsor y un flector).

Mtodo de la Rigidez utilizando una

computadora

Una de las caractersticas ms importantes del mtodo de la rigidez es la forma e

n que las
propiedades elsticas de las piezas, y su orientacin dentro de la estructura, son intro
ducidas
en el clculo antes de que se efecte ninguna consideracin sobre el eq
uilibrio o la
compatibilidad de los nudos.
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 3

Esto nos permite establecer relaciones entre las fuerzas de extremo de b

arras y los
desplazamientos de nudo. Estas relaciones expresadas en forma matricial se de
nomina o
conforma la matriz de rigidez de barra.
Al considerar la interrelacin de cada barra con las dems se obtiene un sistema
global de
ecuaciones que define el comportamiento de toda la estructura y nos conduce a la solu
cin del
problema.
Podemos considerar seis etapas fundamentales en la solucin de un problema:
1) Identificacin estructural
2) Clculo de la matriz de rigidez de barra y del vector de cargas nodales equivalentes
3) Clculo de la matriz de rigidez global y del vector de cargas global de la estructura.
4) Introduccin de las condiciones de borde
5) Solucin del sistema de ecuaciones
6) Clculo de solicitaciones en los extremos de barras y reacciones nodales.

Identificacin estructural.
Esta etapa consiste en definir a travs de nmeros y datos las barras de la estructura.
a) Definir un sistema de ejes globales para la estructura. Las coordenadas de los
nudos se
refieren a dicho sistema.
b) Conectividad de los elementos, identificando para cada barra el nudo inicial y
el final. A
cada barra est asociado un sistema de ejes locales al cual se refieren todas las dime
nsiones
y caractersticas de la barra. El mismo queda definido automticamente por
el orden
establecido para la numeracin de los nudos de la barra.
El eje x local coincide con el eje geomtrico de la barra, siendo el sentido positivo el qu
e va del
nudo inicial (nudo de menor numeracin) al final (nudo de mayor numeracin). Los otr
os ejes

locales debern coincidir con los ejes principales de Inercia de la seccin transver
sal de la
barra formando un triedro directo.
c) Propiedades de la seccin transversal de cada barra. Dependiendo del tipo de
estructura
(reticulado, prtico plano, prtico espacial, emparrillado) se debe dar el rea de l
a seccin
transversal, los momentos de inercia en relacin a los ejes principales y la inercia a la t
orsin.
d) Propiedades del material. Se debe indicar, para cada barra, el mdulo de
elasticidad
longitudinal y/o el mdulo de elasticidad transversal.
e) Especificacin de los vnculos: se debe indicar el nombre del nudo que tiene u
na o ms
restricciones y cuales son las mismas.
f) Descripcin de la carga: se da el nombre del nudo y los componentes de globa
les de las
cargas externas y las reacciones de empotramiento perfecto en relacin a los ejes loc
ales de
la barra, si hay cargas en el tramo.

Matriz de Rigidez y Vector de Cargas Nodales Equiv.

a) Barra de reticulado plano
Consideremos una barra de reticulado plano, supongamos que la misma est arbitrari
amente
orientada con relacin a un sistema de ejes globales X e Y.
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 4

Supondremos que la barra es recta, de seccin transversal constante y que

el material
responde a la ley de Hooke.

Fig. n 1-Sistema local de ejes

DXJ inicial es el j y el final es el k , quedan

En Fla
i de 0la figura 0el nudo
X barra

J
L
F
0 x DYJ L
0
do definida la
F

XK
orientacin
de los ejes
0 locales x e y.
F
YK
Considerando
que no existen0 deformaciones
iniciales y que la deformacin es
DYK

elstica el
alargamiento de la barra i estar dado por:
L

(1)

Donde DxkL y DxjL son los desplazamientos del nudo k y j respectivamente en la direcci
n local
l
x.
Para una barra de reticulado existe una sola solicitacin posible que es el esfue
rzo axil o
normal.
Suponiendo un material elstico lineal sometido a esfuerzo de traccin tendremos
para los
nudos j y k respectivamente:

Li DXk EA
DX j
Li
FXJ
L

FXJ EA
L

FXK EA
L

FXk

EA L
i
L

(2)

LL

(3)

LL

(4)

(DXK DXJ )

Donde:

E= Mdulo de elasticidad

DXJ )
L= Longitud (D
deXK
la barra
A= Area de la seccin transversal de la barra.
l

Como en la direccin y para barras de reticulado no existen solicitaciones podemos e

xpresar
las ecuaciones anteriores en forma matricial:
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 5

EA

EA L 0

L
EA

00

EA
L
0

00

DXK
L

(5)

La expresin (5) corresponde a la ecuacin matricial de la barra i en coordenadas

locales y
L
expresa las Fuerzas de extremo de b en funcin de los desplazamientos de nud
arra FI

DI

os

L
L
A la Matriz que relacion y D se la denomina matriz de Rigidez de barra de reticula
en coordenadas
locales SI . Expresando
do en forma compacta o simblica:
a FI
I
~

FI

SI
LL

.DI

~
L

(6)

~
~

La ecuacin (6) define las fuerzas de extremo Fj y Fk para cualquier pareja de corrimie
ntos Dj ,
Dk dados. Estas ecuaciones son simtricas, como podamos esperar a partir del teo
rema de
reciprocidad. No es posible, sin embargo, resolverlas y obtener los desplazamiento
s (D) en
trminos de las fuerzas (F), puesto que la matriz S es singular. Esto refleja el hecho d
e que la
pieza puede sufrir un movimiento arbitrario de conjunto, sin afectar las fuerzas que act
an en
sus extremos.

Interpretacin fsica de la Matriz de Rigidez de barra

Si en la ecuacin (5), hacemos nulos todos los desplazamientos excepto DxjL que es ig
ual a la
unidad, entonces los esfuerzos en los extremos de la barra sern los indicados en la fig
ura

Fig.n2 - Corrimiento unitario Dxj

EA
L

EA
L

ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 6

que corresponden a la primera columna de S.

De la misma forma podemos hacer DyjL =1 y el resto de los corrimientos nulos, siendo
en este
caso nulos los esfuerzos en los extremos de barra, ya que se considera una barra dobl
emente
articulada y pequeos desplazamientos. Por esta razn los cuatro valores de la
segunda
columna son nulos.
Si en cambio hacemos DxKL =1 y el resto de los desplazamientos nulos, los esfuerzos s
ern:

FXJ EA
L

FYJ 0

FXK

EA
L

FYK 0

(8)

(9)

que corresponden a la tercera columna de la matriz de rigidez de la barra i

Fig. n3 - Corrimiento unitario Dxk

FXJ

FXK
(7)
y
En forma anloga se puede analizar la cuarta columna aplicando un desplazamiento Dyk
=1

Asociando los desplazamientos y reacciones de nudos en las direcciones indicad

as en la
figura, podemos deducir el significado fsico de la matriz de rigidez de la barra.

1 EA 0
L
2
0 0 0
3
EA
EA
0
4

0
0=

Con lo cual podemos observar que los elementos Sij de la matriz de rigidez, repres
entan las
fuerzas que se generan en i al aplicar un desplazamiento unitario en j, permaneciendo
fijos los
restantes.
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 7

Adems para un desplazamiento del nudo k obtenemos una reaccin en j que es

la misma
que la obtenida en k para un desplazamiento en j, lo cual nos es expresado por la si
metra de
la matriz de rigidez.

SI
~

EA
L

(10)

L
0

L
0

Tambin podemos ver que una columna j est formada por las reacciones deb
idas a un
desplazamiento unitario impuesto en la direccin j, y una fila i no es ms que las reacci
ones en
i debido a corrimientos unitarios impuestos en las distintas direcciones.

b) Barra de Prtico Plano

En base al significado fsico de los elementos de la matriz de rigidez, deduciremos la M
atriz de
Rigidez para una barra de Prtico Plano en coordenadas locales.

Para este tipo de elemento corresponden tres desplazamientos por nudo (2 traslaciones
y una
rotacin en el plano).

Fig. n 4
La matriz de rigidez se obtiene dando desplazamientos unitarios de a uno por
vez en las
direcciones de la figura mientras los otros permanecen nulos.
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

FXj
L

FZk

2EI DYj
L . DZj

Pg. 8

6EI
L2

Lyk
D
Zk

Fig. n 5 - Desplazamientos unitarios p/ prtico plano

Las reacciones mostradas en la figura n5 constituyen las respectivas columnas de l
a matriz
de rigidez de la barra de Prtico Plano de la ecuacin (11)

F
Y

yk

EA

FL
Fkj

EA
L

12EI

6EI

6EI
2
L

4EI
L

12EI

6EI

6EI
2
L

2EI
L

EA
L
0
0
EA
L
0
0

0
6EI

12EI

DXk
DL

4EI
L
6EI D
L
Xj
2
2
L
L

L
0
12EI
3

L
6EI
2
L

(11)

Escribiendo en forma compacta:

LL

(12)

~
~
~

Esta matriz relaciona las fuerzas de extremo de barra con los desplazamientos n
odales en
ejes locales.
L
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

FI

SI

Pg. 9

.DI

Fig.n 6 - Solicitaciones para barra de Prtico Plano

Cargas nodales equivalentes

Hasta ahora hemos supuesto que las cargas estaban aplicadas en los nudos, y po
r lo tanto
existe una correspondencia biunvoca entre los puntos de aplicacin de las car
gas y los

desplazamientos que estn siendo calculados. Si esto no ocurriera, por ejemplo tu

viramos
cargas en el tramo de las barras, en forma distribuida o concentrada, debemos su
stituir las
F F F D D
cargas en las mismas por un sistema de cargas equivalentes aplicadas en los
nodos que
produzca en la estructura el mismo efecto que las cargas originales.
Aplicando el principio de superposicin, que es vlido por haber supuesto que el si
stema es
lineal, podemos descomponer las cargas tal como se indica en la figura:

Fig.n 7 - Barra de prtico con cargas en el tramo

Como podemos observar las cargas, reacciones y deformaciones de la estructura
a) sern
equivalentes a la suma de los dos estados b) y c).
Como las deformaciones de nodos en b) son nulas, sern iguales las deformacion
es de los
casos c) y a). O sea que las cargas de c) producen la misma respuesta estru
ctural en lo
referente a desplazamientos de nudos que las cargas originales. Estas sern ento
nces las
cargas equivalentes en los nudos, que no son ms que las reacciones de emp
otramiento
perfecto cambiadas de signo.
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 10

Los esfuerzos en los extremos de barra se obtienen por la suma de los casos (b) y (c).
abca
~
~
~
~

Por lo tanto a la ecuacin (12) habr que adicionarle las fuerzas de empotramiento
perfecto
del caso b).
LLL

FI ~ SI

. DI AI

AI

(13)

representa
el vector de fuerzas de empotramiento perfecto de la barra en coordenadas
~

~
L

locales.

Rotacin de ejes en el plano

Hasta el momento expresamos la matriz de rigidez del elemento barra, tanto de
reticulado
plano como de prtico plano, segn un sistema de ejes locales, estando los desplazami
entos y
esfuerzos de extremo de barra referidos a los mismos.
No obstante, segn ya lo mencionamos, antes de poder aplicar las condiciones de
equilibrio
en los nudos y de compatibilidad de desplazamientos es necesario transformar las f
uerzas y
corrimientos a un sistema de ejes globales.

Fig. n 8 - Rotacin de un vector

Supongamos el vector V de la figura referido al sistema de ejes X e Y.

Las componentes del mismo sern

Vx=V.cos

Vy=V.sen
L

(14)

En el sistema de ejes x e y , las componentes sern

VxL=V.cos(-)
VyL=V.sen(-)

cos
sen

V
X
X
ESTABILIDAD
III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS
L
.

sen cos VY

(15)
Pg. 11

Luego:

V VxLR=V.cos
. V .cos + V.sen .sen

(16)

VyL=V.sen .cos - V.cos .sen

(17)

Teniendo en cuenta la ecuacin (14) y escribiendo en forma matricial:

(18)

V
VX Y cos sen VX

V
cos VYL
o en forma
compacta:
sen
L

(19)

VR
~

V~

A la matriz de cosenos directores R la llamaremos matriz de rotacin.

~
~

Si queremos el pasaje del sistema local al global, en este caso la matriz de rota

V es
R la
cin
V de R .
transpuesta
~

(20)

(21)

~
~
~

Premultiplicando la ecuacin (19) p

or R
1
~
~

(22)

Lo que implica que R es una matriz ortogonal (su inversa es igual a su traspuesta):
Donde D 1 y DT

R
~
R

(23)

Ecuacin fundamental de la barra en el sistema global

Tanto
solicitaciones
como los desplazamientos pueden expresarse como vector
DI lasRT.
DI
es en el
plano, podemos entonces aplicar la transformacin lineal antes vista para llevar los esfu
erzos y
desplazamientos
RT ~ R de extremos de barra del sistema local al global.
L ~0
~
(24)

~J

~ J

~ K

~K

DI

contiene los desplazamientos de los dos nudos de la barra en coordenadas locales y Di

ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 12

representan respectivamente los desplazamientos de los nudos j y

k en

~J ~K

coordenadas locales y D ., D representan los mismos desplazamientos en coorde

nadas

~J

globales.

~K

Escribiendo en forma compacta:

L

(25)

~
R

0
(26)

~
Siend
o:

R.D
~

R.D

La matriz de rotacin de la barra por incluir ambos nudos de la misma.

L
~

los desplazamientos de los mismos en coordenadas globales.

Para las solicitaciones tendremos:
L

(27)

~
~
~
L

~
~
en coordenadas globales.
La inversa de la matriz de rotacin de la barra ser:

FI

FI como
RT.FI =
Pero
R

ser

R
~ ~

RT1
~

contiene las solicitaciones de los dos extremos de la barra en coordenadas locales y FI

LLL
L
(29)

~
~

RT
~

~
0
~

RT

~
FI

SI
~

(28)

R~

. DIL AI

(30)

FI RT .FI
~
FI RT .SI
~
~
~

. DI RT . AI

TLLT

(31)

SI RT .SI
T L T . RT
~

(32)

AI RT . AI
~
~

RT ~

(33)

sen cos

sen cos

~
~
~
~
~
~
~
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

TL

Pg. 13

(34)

~
~
~

(35)

~
Matriz de rotacin para barra de reticulado plano
FI RT .SI . RT. DI RT . AI
cos sen
FI SI . DI AI

0
0

0
0

c.s

c.s

2
c c.s
.EA/L
2
c.s
c.s c2

2
2
s
c.s
s

c.s

(37)

Matriz de rotacin para barra de prtico plano

cos

~
Si
~

0
0
2
c0
0

sen 0

0
0
0
0

0
(38)

1
0
0
0
0 cos
s2
0 c.ssen cos 0

0
0
0
1

ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 14

0
sen cos 0
0
0
Matriz de Rigidez de barra de Prtico Plano
RT
sen 0

(EAc / L 12 (EAcs / L 12EI 6EIs / (EAc / L 12EIs / L ) (EAcs / L 1

2
3
3
2
/ L3 ) 6EIs / L2
EIs / L )
cs / L )
L
2EIcs
2
2
2
(EAcs / L 12EI (EAs / L 12EIc 6EIc / (EAcs / L 12E (EAs / L 12EIc 6EIc / L
3
2
3
2
3
2
3
cs / L )
/L )
L
Ics / L )
/L )
2EI / L
S
2
2
2
2
2
2
2
6EIc / L 4EI /
6EIs / L
6EIc / L 6EIs / L
6EIs / L
I 2
2
(EAcs / L 12EI L
(EAc / L 12EI (EAcs / L 12EI

2
3
2
3
3
~ (EAc / L 12EI cs / L )
6EIs / s / L )
cs / L )
2
2
3
L
s /L )
2
2
3
2
(EAcs / L (EAs / L 12EIc / L ) (EAcs / L 12EI (EAs / L 12EI 6EIc /
3
2
3
2
3
2
12EIcs / L ) 6EIc / L
cs / L )
c /L )
L
2
2
2
2
6EIs / L 6EIc / L
2EI / L
6EIs / L
6EIc / L
4EI / L
Tabla n 1
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Matriz de Rigidez y Vector de Cargas Global de la Estruct.

Pg. 17

Fig. 9 - Reticulado Plano

Las ecuaciones de equilibrio exigen que las cargas externas aplicadas en los nudos debe
n
ser iguales a la suma de las solicitaciones de extremo de las barras que concurren al nud
o.
Siendo el vector de cargas externas aplicadas en j :

A j Xj
~ A
Yj

(39)

La ecuaciones matriciales de las barras a) y b) son:

a

(40)

FI~ SI~ . DI
~ AI
~

FI~ SI~ . DI
~ AI
~

Que pueden representarse en funcin de los nudos de la siguiente forma:

~ a Sij
~Fi a Sii
~ a .
Fja Sji a Sjj


~
~
b Sjk
~Fjb
~Sjj ~
b
Fk b Skj
Skk


~
~

a
a

Di
~ AI
~ i
a
Dj AI

Djb AIb
~
.
~ bj
b
AI

Dk

Sij submatriz conteniendo los coeficientes de rigidez del nudo i que provienen de
dnde:
un
a

~
a

~
a

~
17
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 18

~
desplazamiento unitario del nudo j de la barra (a)
condiciones de equilibrio para el reticulado de la figura en el nudo j resultan:
Las

A j Fj Fj

A j Fj Fj

A j Sji
AI
D

. Di Sjj . Dj AI Sjj . Dj Sjk

D D

. Dk AI

(42)

~
~
a

A~j ~Sji ~. Di (Sjj Sjj ). Dj Sjk

AI
~

. Dk

(43)

de compatibilidad
de desplazamientos:
FiPor condicin
vector solicitacin
en el extremo
i de la barra (a)

Di

a
b
ac
bc
;
;
~
j
~
j
~
j
~
i
~
i
~
i
~
k
vector deformacin en el extremo i de la kbarra~(a)

~k

(44)

Reemplazando
la (44)
la (43) y extendiendo
el razonamiento
a los
nudos(a)restantes:
AIi
vector fuerza
deen
empotramiento
perfecto en
el extremo i de
la barra
a

~
A i
~~

A j
~
~

a
c
Sii a
Sii
Sij
~
~
~
~
~
~

c Sji
Sjj
b Sjjc
~ ~~
~
~
Skj
Ski

~ j ~ j

c

c
a Di c AIa AI
Sik
i
i
~
~
~ i ~ i

.
Dj
AI
AI
b

a
b

c j
b
b
j
~
k AIc
~Skk ~ Skk~ ~ k ~ AI

Dk~ k
k
~

(45)

Si la organizamos en forma matricial:

S .D

18
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

c
a

a
c

A k

~
~ ~

Pg. 19

~~
Sjk

A irepresenta
Sij . Dj
Sii global
). Dide
Sikecuaciones
. Dk AI de
la estructura, escribiendo en fo
Que
el (Sii
sistema
AI
rma

compacta:
A k Ski . Di (Skk Skk ). Dk Skj . Dk AI A
I

A A
~ eq

~ ~

Dnde:

A vector de acciones externas aplicadas en los nudos.

~
S = matriz de rigidez global de la estructura
~

D = vector deformacin de nudos

~

(47)

A = vector de cargas nodales equivalentes (fuerzas de empotramiento perf

~ eq ecto
cambiadas de signo) debido a las cargas aplicadas en los tramos de barras.
Haciendo:
(48)

A A A
~ ~ ~ eq
~

(49)

~ ~

El vector A representa las cargas nodales equivalentes ms las cargas aplicadas en

los
nudos.
La ecuacin (49) no puede ser resuelta (la matriz S es singular) si no se aplic
an las
condiciones de contorno o de borde de la estructura.

Esquema de la formacin de la matriz de rigidez global y del vector

de cargas
Consideremos el reticulado de la figura n 9. La matriz de rigidez global del retic
ulado
tendr 6 filas y 6 columnas ya que son 3 nudos y cada uno tiene 2 grados de libertad
(un
desplazamiento segn el eje x y uno segn el eje y). La contribucin de cada nudo debe
r
ser colocada en las posiciones que se indican en la figura n10.

S . Di,j,k definen las filas y columnas de la matriz (submatrices)

LosAnudos

19
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

I
i
S=

j
k

Pg. 20

a+c

a+c

a+c

a+c

a+b

a+b

a+b

a+b

b+c

b+c

b+c

b+c

Fig. n 10: Esquema Formacin Matriz Reticulado

J
c
~ ii

~ ij

~ ik

~ ii

S=

j
k

aa

~ ji ~ jj ~ jj

~ jk

bc

~ ki

~ kj ~ kk ~ kk

Las submatrices tendrn tantas filas y columnas como grados de libertad tenga cada nud
o.
Sa + S
S
S
i
En el caso del reticulado este nmero es 2.

~ jj

b~

~ jk

2(j-1)+2

Sb

Sc

~ ki

~ kk

2(k-)+2

a
~ ij
a

2(k-)+1

a~

Sa

2(j-1)+1

b
Analicemos la posicin en figura
S n 10 de Slas +submatrices
S
S de la matriz de rigidez de la bar
ra
b
b de la figura n 9.
S
S
S +S
c

~ ik
c

(50)

Las submatrices son de orden 2x2. La posicin S 1 en la matriz S est indicada

en al

S
S
~
~
ii S
S1 a
S1 b

S1

ii

jj
~
c
b
figura n 10 a partir del siguiente esquema:

S~ ji S

S~ kj

~ kk

~ jj
b
~ jk
b

S ~ kj= indica fila j columna j

b

S ~ kk
= indica fila j columna k
Ntese que en las posiciones dnde ya existen valores provenientes de otras barras, esto
= indica fila k columna j
s
son sumados a aquellos.
S = indica fila k columna k

20
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 21

Generalizando, debemos considerar una barra de reticulado que une los nudos j y k seg
n
la figura n11.

Figura n 11
Las posiciones ocupadas por la matriz de rigidez de la barra (m) en la matriz de ri
gidez
global son las siguientes:

1
2
.
.

S
~

2(j-1)+1

2(j-1)+2

=
2(k-1)+1

2(k-1)+2

Figura n 12 - Posiciones ocupadas por una barra de Reticulado Plano en la m

atriz
de Rigidez Global

3(k-1)+3

3(k-1)+2

3(k-1)+1

3(j-1)+3

os la
siguiente configuracin:

3(j-1)+2

3(j-1)+1

Para un prtico plano o estructuras con 3 grados de libertad por nudo tendram

21
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 22

1
2
.
.

S
~

3(j-1)+1

3(j-1)+2

3(j-1)+3

=
3(k-1)+1

3(k-1)+2

3(k-1)+3

Figura n 13 - Posiciones ocupadas por una barra de Prtico Plano en la matriz

de
Rigidez Global
Del mismo modo, un vector de cargas equivalentes en los nudos para una barra de la figu
ra
n11 ocupa la posicin indicada en la figura n 14 del vector acciones globales nodales.
1
2
.

3(j-1)+3

3(j-1)+1

3(j-1)+2

A =
~

1
2
.
.

3(k-1)+1 X
3(k-1)+2 X
3(k-1)+3 X

a) Prtico plano

A =

2(j-1)+1

2(j-1)+2

b)

2(k-1)+1

2(k-1)+2

Reticulado

plano

Figura n 14 Posicin del vector Acciones nodales

22
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Formacin de la MATRIZ de
banda-simtrica.

Pg. 23

RIGIDEZ GLOBAL en forma de matriz

La Matriz de Rigidez Global de la estructura es simtrica respecto de la diagonal

principal y
tiene normalmente caractersticas de matriz de bandas. Por ello no es necesario al
macenar
todos los elementos de la matriz toda vez que esto implica tiempo de procesamient
o para el
computador y adems mayor memoria de almacenamiento. Bastar entonces almac
enar los
elementos del tringulo superior (o el inferior) de la matriz. Tambin esta parte
triangular
superior (o el inferior) puede sufrir una racionalizacin del arreglo matricial de nmer
os. Esto
se corrige aprovechando la caracterstica de banda de la matriz y pre
viniendo su
reorganizacin para un arreglo bidimensional de dimensin NGLTxB dnde:
NGLT= es el nmero de grados de libertad total de la estructura.
B= es el semi largo de banda del sistema.
Sea el ejemplo de la figura n 15. La Matriz de Rigidez de la estructura en s
u forma de

almacenamiento convencional representada en la figura n 16-a y en la forma

de banda
simtrica en la figura n 16-b.

Figura n 15: Ejemplo de Reticulado Plano.

Figura n 16: Formas de almacenamiento de la Matriz de Rigidez

Ntese que el semi-ancho de banda est vinculado a la mayor diferencia de
numeracin entre
nudos que estn unidos por una barra. En el caso de la figura n 15, la mayor diferen
cia entre
dos nudos de una barra es 2 y el semi-ancho de banda es 6. De modo general el
semi-ancho
de banda es obtenida por B=(L+1)*N1. Dnde L es la mayor diferencia entre
nudos que
2
3
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 24

concurren a una misma barra y N1 es el nmero de grados de libertad del nudo. En el

caso del
Reticulado Plano N1=2, y entonces el semi ancho de banda es B=(2+1)x2=6.

Influencia de la Numeracin de nudos sobre el ancho de Banda

En el ejemplo que sigue se puede apreciar como varan tanto la cantidad de ele
mentos
a almacenar como el semiancho ancho de banda con la forma de numerar los n
udos de
la estructura.

Condiciones de Contorno o de Borde

Un sistema de ecuaciones : A S . D correspondiente a una estructura complet

a antes de

~
~ ~

(52) La
.
aplicar las~ condiciones
de contorno es INDETERMINADO, pues S es singular.
S
S
A

~ I

razn de~ II,I

~ II,II

esta singularidad es el resultado de no haber considerado las vinculaciones o apo

yos de la
estructura con el exterior. Al introducir las condiciones de vnculo de
saparece la
indeterminacin (o singularidad de la matriz de rigidez), siempre que el nmero de
vnculos
sea por lo menos el mnimo necesario para eliminar los movimientos de cuerpo rg
ido de la
estructura.
El conocimiento de determinados corrimientos nodales, disminuye el nmero de in
cgnitas,
tornndose innecesarias las ecuaciones correspondientes a estos corrimientos.
La eliminacin de la ecuacin de un desplazamiento implica la destruccin de la ba
nda de la
matriz, lo cual exigira un reacomodamiento de las incgnitas del problema.
Al sistema de ecuaciones

AS.D
~

(51)

Podemos expresarlo de la siguiente forma:

2
4
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 25

tambin pasar al vector A los coeficientes Sij . di con el signo cambiado.

~ II

~ I,I

~ I,II

~
~

~ II

Dnde:

A = reacciones de apoyo con desplazamiento impedido.

D =0, reacciones

que por el
~ I

~ I

momento desconocemos.

= agrupa las fuerzas externas sobre nodos con desplazamientos D

.
~ II

~ II

desconocidos

= Matriz de rigidez (cuadrada) que relaciona fuerzas conocidas (cargas extern

as A

~ II,II
~ II

con desplazamientos desconocidos D

~ II

Tcnica nmero 1:
Suprimir o eliminar de la Matriz de Rigidez del sistema las columnas correspo
ndientes a
desplazamientos impedidos y las filas correspondientes a reacciones exteriores.
Ventaja: reduce el tamao del sistema de ecuaciones a resolver (en un orden determin
ado por
los grados de libertad eliminados por los apoyos). Dicha reduccin es importante en el
caso de
pocas barras y muchos vnculos externos, fcilmente aplicable en forma manual, es la
tcnica
que emplearemos en el curso.
Desventaja: no podemos resolver as desplazamientos prescriptos.
Tcnica nmero 2:
Un artificio usado normalmente para introducir un desplazamiento conocido en la dir
eccin i,
consiste en hacer el elemento de la diagonal principal de la matriz S de la lnea i ig
ual a 1 y
~

anular todas las restantes posiciones pertenecientes a esa lnea y columna. Luego
debe ser

colocado el valor del corrimiento conocido d i en la posicin anteriormente ocupada

por ai, y

S11 .d1 S12 .d 2 .... S1i .di .... S1n .dn a1

a2
S 21 .d1 S 22 .d 2 .... S 2i .di .... S 2n .dn
.

.
S .d S .d .... S .d .... S .d
i2
2
ii
i
in
n
i
a i1 1
.

Sn1 .d1 Sn2 .d 2 .... Sni .di .... Snn .dn an

(53)

2
5

ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

.
.

Pg. 26

S11 .d1 S12 .d 2 .... 0 ....

S1n .dn a1 S1i .di

S 21 .d1 S 22 .d 2 .... 0 ....

S 2n .dn a 2 S 2i .di

(54)

.... di .... 0 di

Sn1 .d1 Sn2 .d 2 .... 0 ....

En forma matricial

Snn .dn an Sni .d

di
d. 1
...
S1n d. 2
.
...
S 2n
... .

(55)

S12 ...
S11 S
22
S 21 . ...
.

...

...
...
.
.
...

.
.
... . . .

1 ... 0
di
.... S. 2n.
.
....
.. .
0 ...
Snn . .

an Sni .di
di
.

. .
0 0
. .
. .
Sn1 Sn2

dn

a1 S1i .di

a 2 S 2i . di
.

Como caso particular, cuando un corrimiento di sea nulo, o sea, cuando exista un vncu
lo en la
direccin i, la ecuacin anterior toma la siguiente forma:

S11 S12 ... 0 ...

S 21 S 22 ... 0 d1
. . ... .
.
0
.

.
0
.

... . ...
... 1 ...
. . .

S1n

.
.
0
.

d2
.
.

a1
a2
.

0
.

(56)

. .
Sn1 Sn2

. . ..
... 0 ...

.
Snn an

dn
Es importante mencionar, que con esta forma de aplicacin de las condiciones de b
orde, las
cargas que estn aplicadas en la direccin i del apoyo quedaran sin considerar.
Ventaja: No destruye la organizacin simtrica y banda del sistema de ecua
ciones. Es
fcilmente implementable en computador.
Tcnica nmero 3:
Modificar los coeficientes de la diagonal principal en la matriz de rigidez del si
stema que
correspondan a desplazamientos impedidos.
2
6
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 27

Siendo una ecuacin cualquiera en el sistema:

ai

Sij .d j Sii .d

di a i
Sij .d j
Sii

(57)

20

Luego, haciendo Sii suficientemente grande (por ejemplo multiplicarlo por 10 ) obtene
mos un
di prximo a cero.
Ventajas: no altera la organizacin del sistema de ecuaciones, ni tampoco reduce el or
den del
sistema.

Desventaja: puede provocar problemas numricos.

Solucin del sistema global de ecuaciones

Una vez introducidas las condiciones de borde o contorno en el sistema de e
cuaciones,
podemos resolverlo mediante cualquier mtodo conocido (Gauss, Cholesky, Gauss
Jordan,
etc.)

Existen en la literatura mtodos eficientes de clculo teniendo en cuenta la simetra,

banda y
forma especfica de archivo. La solucin de la ecuacin 55 nos provee los corri
mientos o
desplazamientos de los nudos de la estructura.
Muchas veces existe ms de un estado de carga. Para esta situacin se pueden d
esarrollar
algoritmos de solucin que trabajando con ms de un vector de car
gas (trmino
independiente) nos dan como respuesta ms de un vector desplazamiento (uno p
ara cada
estado de cargas).

Clculo de solicitaciones en extremos de barras

La ecuacin matricial de una barra (m) con conexiones i,j se expresa en forma particio
nada de
acuerdo a la ecuacin 41:
m

mmm

~ i
m

mmm

~ j

Los desplazamientos
dos de la

D ,
~ i

D son obtenidos del vector desplazamiento de nu

~ j

estructura, solucin del sistema de ecuaciones (54). Particionada la matriz de rigid

ez de la
m
~ i,i

m
~ i,j

m
~ j,j

mm
~ i~ j

almacenadas al confeccionar la matriz de rigidez global y el vector de car

F S .Di S .D j A
gas global.
Entretanto, dependiendo de la capacidad de memoria del procesador, puede
S .Di S .D j A
F ms
ser
aconsejable recalcular la matriz de rigidez y el vector de reaccin de empotramiento
perfecto
de cada barra, para despus obtener con el trmino (58) las solicitaciones de extremo d
e barra
m

barra ~ im ~ j (S

,S

, S) y el vector de empotramiento perfecto (A

, A ) , deben
ser
Las solicitaciones de barras nos interesan referirlas al sistema de ejes de coordenadas l

ocales,
ya que en esta forma tendremos los esfuerzos normales, cortantes, etc. Para esto
debemos
(F hacer
y F una
) . rotacin de ejes sobre las solicitaciones obtenidas por (58), para lo cual:
2
7
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 28

 R .F

m, ~
m
L R . F .
~ i
~ i
m

m,
L
~ j

(59)

~
~ j

Esta no es la nica forma de obtener las solicitaciones de las barras. Podemos o

btener los
corrimientos D y D y transformarlos a coordenadas locales mediante el product
o con R .
~ i

~ j
~

Entonces tendramos las ecuaciones (58), todas en coordenadas locales para as obte
ner las
solicitaciones en las mismas coordendadas.

Clculo de las reacciones de nudos

Para un nudo sobre un apoyo, podemos obtener las acciones de todas las
barras que
concurren a un nudo y sumarlas.
El resultado de esta suma ser el vector reaccin del nudo.

F~ mi

m1

p
~
i

~
i

~ i

p
m1

(60)

m1

~ i

dnde p es el nmero de barras ligadas al nudo i y R

~ i

es la reaccin del nudo i.

Si hacemos esta misma suma para un nudo que no sea apoyo, tendremos como result
ado las
cargas aplicadas al nudo. Es importante que las solicitaciones del extremo de b
arra sean
referidas a un sistema global de coordenadas.

2
8
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 29

TOPICOS ESPECIALES EN ANALISIS

MATRICIAL
Efectos Trmicos
Lo analizamos en forma similar al caso de cargas aplicadas en el tramo.

F~ mi

Fig. n 17

Xj (T1 T2 .L
)
t
2

El giro d de un tramo de viga dx ser:

j
t

(T1 T2
.L
)
h

(61)

( T1 T
d
.dx
2 )
t
h

(62)

teniendo en cuenta q
ue:

d d y M
2
dx
E.I

2
d (T1 T2 )
t
y
h

(63)

integrando dos veces y particularizando para el nudo j.

dy j
t

( T1 T 2
.L
2 )
2.h

(64)

Las fuerzas de empotramiento necesarias para restaurar el extremo de la barra a su p

osicin
primitiva sern: (para prtico plano).

NLj
QLj
M Lj

0 .L

2
L Lh

dx
dx

E.A

0
12.E.I
3

L
6.E.I

0
6.E.I
L
4.E.I

(T T2 )
.L
t . 1
(T1 T2 )
2.h
(T1 T2 )
.
t

(65)
2
9
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

t .(T1 T2 ).
P S .
~ j

Pg. 30

(66)

~ jj ~ j

t .(T1 T2 ).

Cmo estas fuerzas son de empotramiento perfecto, y adems estn en coordenadas l

ocales,
hay que llevarlas al sistema global y restarlas (cargas nodales equivalentes) al
vector de
cargas en los nudos.

NLj
QLj
M Lj

E.A
2

(67)

E.A
h

Apoyos elsticos

Un apoyo elstico se caracteriza por el hecho

de que el desplazamiento que sufre es
proporcional a la fuerza que recibe, definida por
su constante de resorte k (fuerza necesaria
para producir un desplazamiento unitario)
Un desplazamiento genera una fuerza k. en
la misma direccin.
Siendo una ecuacin cualquiera en el sistema:

Fig. n 18

ai

Sij j Sii .i

(68)
2

Si establecemos ahora el equilibrio

en cuenta la fuerza del resort
t .
.L en el nudo i. teniendo
e:
0

k.i ai

Sij j Sii .i

ai

Sij j (Sii k).i

(69)
(70)

Es decir debemos sumar la rigidez del resorte k en la diagonal principal de la matriz d

e rigidez
del sistema en la fila y columna correspondiente a la direccin en la que acta el resort
e.

3
0
ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 31

Apoyos deslizantes no concordantes

Fig. n 19
Si el plano de desplazamiento no es paralelo a uno de los ejes globales, debemos mo
dificar la
ecuacin matricial de la estructura, transformando en el nodo en cuestin los vectores f
uerza y
desplazamiento al sistema de ejes global.

R
.
~1

~ 1
~1

F
F

R .L

~ 1~ 1

(71)

~ 1

dnde el 1 indica el nmero de nodo.

Con

cos sen 0
R sen cos 0
0
0 1

(72)

R .S
R R. L .S ~ 11. R
~ 12 ~ ~ ~ 1
~1
F
SSS0
~ 1
~ 21 ~ 22 ~ 23 ~ ~ 2

0
SSS
~ 1
~ ~ 32 ~ 33 ~ 34
F
00SS
~

0L

(74)

~ ~ ~ 43 ~ 44

F
~

F
~

S =3x3 para prtico plano

~ ij

En esta ecuacin es lo mismo premultip

licar

por

que postmultiplicar la prim

era
~ 1

columna de S por R . Por otra parte, premultiplicando la primera fila a ambos l

ados de la
~

~ 1

T
~
1

no se altera la igualdad y adems considerando que en siste

mas de

coordenadas ortogonales R
~

.R

1 ~ 1

3
1

R .

ESTABILIDAD III CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS

Pg. 32

(73)
~3

L
~ 1

TT
~ 1

~ 2

~ 11

~ 3

~ 1~

F
F

~ 4

~4

~ ~ 32 ~ 33 ~ 34 ~ 3
~ ~ ~ 43 ~ 44 ~ 4

~ 12 ~
~~1

igualdad por R
Resumen: Se modifica la ecuacin matricial de toda la estructura, premultiplicand
o la fila,
T
correspondiente al nodo con apoyo inclinado por R
y postmultiplicando l
a columna
~

correspondiente por R .
~ 1

Tcnica alternativa
Consiste en adicionar un barra ficticia que deber impedir desplazamientos en la direc
cin v
y permitir desplazamientos en la direccin u y giro. Para que esto sea posible, el r
ea de la
seccin transversal de la misma deber ser mucho mayor que las barras restantes y la
inercia
mucho menor a la de las otras barras. (Adems de pequea longitud para evitar acorta
mientos
de la barra ficticia)

Fig. n 20

32