E3 cap4

ANALISIS MATRICIAL DE ESTRUCTURAS POR
EL METODO DE LA RIGIDEZ

ESTABILIDAD III – CAPITULO IV: ANALISIS MATRICIAL DE ESTRUCTURAS Pág. 1
Introducción
Los métodos clásicos de análisis estructural desarrollados a fines del siglo XIX, tienen las
cualidades de la generalidad, simplicidad lógica y elegancia matemática. Desgraciadamente,
conducían a menudo a cálculos muy laboriosos cuando se los aplicaba en casos prácticos, y
en aquella época, esto era un gran defecto.
Por esta razón sucesivas generaciones de ingenieros se dedicaron a tratar de reducir el
conjunto de cálculos. Muchas técnicas ingeniosas de gran valor práctico fueron apareciendo
(Método de Cross), pero la mayoría de las mismas eran aplicable sólo a determinados tipos de
estructuras.
La principal objeción a los primeros métodos de análisis fue que los mismos conducían a
sistemas con un gran número de ecuaciones lineales, difíciles de resolver manualmente.
Con los computadores, capaces de realizar el trabajo numérico, esta objeción no tiene ahora
sentido, mientras que la generalidad de los métodos permanece. Esto explica por qué los
métodos matriciales deben en su tratamiento básico de las estructuras más al siglo XIX que al
XX.
El empleo de la notación matricial presenta dos ventajas en el cálculo de estructuras. Desde el
punto de vista teórico , permite utilizar métodos de cálculo en forma compacta, precisa y, al
mismo tiempo, completamente general. Esto facilita el tratamiento de la teoría de estructuras
como unidad, sin que los principios fundamentales se vean oscurecidos por operaciones de
cálculo, por un lado, o diferencias físicas entre estructuras, por otro.
Desde el punto de vista práctico, proporciona un sistema apropiado de análisis de estructuras
y determina una base muy conveniente para el desarrollo de programas de computación.
En contraste con estas ventajas, debe admitirse que los métodos matriciales se caracterizan
por una gran cantidad de cálculo sistemático
Las virtudes del cálculo con computadora radican en la eliminación del la preocupación por las
operaciones rutinarias, el ingenio necesario para preparar el modelo con que se pretende
representar la realidad y el análisis crítico de los resultados.
Se debe ser consciente que sin un modelo adecuado o sin una interpretación final, el
refinamiento en el análisis carece de sentido.
Método de la Rigidez
Hipótesis: Estructura lineal- Todos los movimientos y esfuerzos son funciones lineales de las
cargas- Pequeñas deformaciones (ecuaciones de equilibrio en la estructura no distorsionada).
Las barras son rectas y de sección constante.
Para estudiar una estructura por el método de la rigidez, al igual que en cualquier otro
problema elástico, disponemos de tres conjuntos de ecuaciones que deben cumplirse.
Ecuaciones de compatibilidad
Ecuaciones constitutivas
Ecuaciones de equilibrio

Las ecuaciones de compatibilidad relacionan las deformaciones de barras con los
desplazamientos nodales. Introduciendo estas relaciones en las ecuaciones constitutivas,
relacionamos las fuerzas en los extremos de barras con los desplazamientos nodales
Introduciendo estas últimas relaciones en las ecuaciones de equilibrio se obtiene un conjunto
de ecuaciones de fuerzas nodales en función de desplazamientos nodales, que pueden ser
consideradas como Ecuaciones de Equilibrio de la estructura en función de desplazamientos.
La resolución de este sistema de ecuaciones nos permite obtener el valor de las incógnitas
(desplazamientos nodales), a partir de los cuales se obtienen las solicitaciones de las barras
de la estructura, así como las reacciones.
Cuando se van a calcular las relaciones esfuerzos de extremo de barra - desplazamientos, es
natural escoger un sistema de coordenadas que haga estas ecuaciones lo más sencillas
posible.
Tomaremos por lo tanto como eje x el que coincide con el eje geométrico de la pieza y los ejes
y y z coincidentes con los ejes principales de la sección transversal.
Tal sistema pertenece a la barra, y no depende de la orientación de la misma en la estructura 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 común para todos los movimientos y esfuerzos de extremo
de barras para poder aplicar las condiciones de equilibrio y compatibilidad. A dicho sistema lo
denominaremos sistema de ejes globales.
Dichos esfuerzos de extremos de barras y desplazamientos dependerán del tipo de estructura
que estamos resolviendo, para barras de:
a) Reticulado Plano: tendremos dos desplazamientos por nudo
b) Reticulado Espacial: tres desplazamientos por nudo.
En ambos casos sólo tendremos esfuerzos normales.
c) Pórtico Plano: tres desplazamientos por nudo. (una rotación en el plano del pórtico y dos
traslaciones), como solicitaciones de extremo de barra una fuerza axial, un esfuerzo de corte y
un momento flector.
d) Pórtico Espacial: seis desplazamientos por nudo, tres traslaciones y tres rotaciones. 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 plano de la
grilla) y dos rotaciones alrededor de los ejes contenidos en el plano mencionado). Los
esfuerzos son un cortante y dos momentos (un torsor y un flector).
Método de la Rigidez utilizando una
computadora
Una de las características más importantes del método de la rigidez es la forma en que las
propiedades elásticas de las piezas, y su orientación dentro de la estructura, son introducidas
en el cálculo antes de que se efectúe ninguna consideración sobre el equilibrio o la
compatibilidad de los nudos.

Esto nos permite establecer relaciones entre las fuerzas de extremo de barras y los
desplazamientos de nudo. Estas relaciones expresadas en forma matricial se denomina o
conforma la matriz de rigidez de barra.
Al considerar la interrelación de cada barra con las demás se obtiene un sistema global de
ecuaciones que define el comportamiento de toda la estructura y nos conduce a la solución del
problema.
Podemos considerar seis etapas fundamentales en la solución de un problema:
1) Identificación estructural
2) Cálculo de la matriz de rigidez de barra y del vector de cargas nodales equivalentes
3) Cálculo de la matriz de rigidez global y del vector de cargas global de la estructura.
4) Introducción de las condiciones de borde
5) Solución del sistema de ecuaciones
6) Cálculo de solicitaciones en los extremos de barras y reacciones nodales.
Identificación estructural.
Esta etapa consiste en definir a través de números 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 dimensiones
y características de la barra. El mismo queda definido automáticamente por el orden
establecido para la numeración de los nudos de la barra.
El eje x local coincide con el eje geométrico de la barra, siendo el sentido positivo el que va del
nudo inicial (nudo de menor numeración) al final (nudo de mayor numeración). Los otros ejes
locales deberán coincidir con los ejes principales de Inercia de la sección transversal de la
barra formando un triedro directo.
c) Propiedades de la sección transversal de cada barra. Dependiendo del tipo de estructura
(reticulado, pórtico plano, pórtico espacial, emparrillado) se debe dar el área de la sección
transversal, los momentos de inercia en relación a los ejes principales y la inercia a la torsión.
d) Propiedades del material. Se debe indicar, para cada barra, el módulo de elasticidad
longitudinal y/o el módulo de elasticidad transversal.
e) Especificación de los vínculos: se debe indicar el nombre del nudo que tiene una o más
restricciones y cuales son las mismas.
f) Descripción de la carga: se da el nombre del nudo y los componentes de globales de las
cargas externas y las reacciones de empotramiento perfecto en relación a los ejes locales 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é arbitrariamente
orientada con relación a un sistema de ejes globales X e Y.

Supondremos que la barra es recta, de sección transversal constante y que el material
responde a la ley de Hooke.
Fig. nº 1-Sistema local de ejes
En la barra i de la figura el nudo inicial es el j y el final es el k , quedando definida la
orientación de los ejes locales x e y.
Considerando que no existen deformaciones iniciales y que la deformación es elástica el
alargamiento de la barra i estará dado por:
Lj
X
L
ΔLi = DXk − D (1)
L y Dxj
Donde Dxk
L son los desplazamientos del nudo k y j respectivamente en la dirección local
xl.
Para una barra de reticulado existe una sola solicitación posible que es el esfuerzo axil o
normal.
Suponiendo un material elástico lineal sometido a esfuerzo de tracción tendremos para los
nudos j y k respectivamente:
F EA J k
L F EA
X i X Li
= − Δ = Δ (2)
L
L
F EA L
X
LJ(D D )
L
= − − (3)
XJ XK
EA F LJ
(D D )
L
= − (4)
X
L
XK XK
Donde:
E= Módulo de elasticidad
L= Longitud de la barra
A= Area de la sección transversal de la barra.
Como en la dirección yl para barras de reticulado no existen solicitaciones podemos expresar
las ecuaciones anteriores en forma matricial:


    


D
x
    


     


     

EA
−
0 −
EA
=

    


    

J
L
Y
L
J
K
Y
L
D
X
L
X
F
X
F
Y
F
X
J
J
K
Y
K
K
D
D
0
L
L
0 0 0 0
0
0 EA
L
EA
L
0 0 0 0
F
(5)
La expresión (5) corresponde a la ecuación matricial de la barra i en coordenadas locales y
expresa las Fuerzas de extremo de barra L
FI en función de los desplazamientos de nudos
~
L
~
DI
A la Matriz que relaciona L
FI y L
~
DI se la denomina matriz de Rigidez de barra de reticulado
~
en coordenadas locales L
SI . Expresando en forma compacta o simbólica:
~
L
~
L
~
L
~
FI = SI .DI (6)
La ecuación (6) define las fuerzas de extremo Fj y Fk para cualquier pareja de corrimientos Dj ,
Dk dados. Estas ecuaciones son simétricas, como podíamos esperar a partir del teorema de
reciprocidad. No es posible, sin embargo, ”resolverlas” y obtener los desplazamientos (D) en
términos de las fuerzas (F), puesto que la matriz S es singular. Esto refleja el hecho de que la
pieza puede sufrir un movimiento arbitrario de conjunto, sin afectar las fuerzas que actúan en
sus extremos.
Interpretación física de la Matriz de Rigidez de barra
Si en la ecuación (5), hacemos nulos todos los desplazamientos excepto Dxj
L
que es igual a la
unidad, entonces los esfuerzos en los extremos de la barra serán los indicados en la figura
Fig.nº2 - Corrimiento unitario Dxj
y F EA
F EA XJ XK = = − (7)
L
L

que corresponden a la primera columna de S.
De la misma forma podemos hacer Dyj
L
=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 doblemente
articulada y pequeños desplazamientos. Por esta razón los cuatro valores de la segunda
columna son nulos.
Si en cambio hacemos DxK
L
=1 y el resto de los desplazamientos nulos, los esfuerzos serán:
F XJ = − EA y F XK = EA
(8)
L
L
F 0 y F 0 YJ YK = = (9)
que corresponden a la tercera columna de la matriz de rigidez de la barra i
Fig. nº3 - Corrimiento unitario Dxk
En forma análoga se puede analizar la cuarta columna aplicando un desplazamiento Dyk=1
Asociando los desplazamientos y reacciones de nudos en las direcciones indicadas en la
figura, podemos deducir el significado físico de la matriz de rigidez de la barra.
Con lo cual podemos observar que los elementos Sij de la matriz de rigidez, representan las
fuerzas que se generan en i al aplicar un desplazamiento unitario en j, permaneciendo fijos los
restantes.

Además para un desplazamiento del nudo k obtenemos una reacción en j que es la misma
que la obtenida en k para un desplazamiento en j, lo cual nos es expresado por la simetría de
la matriz de rigidez.

     
SIL (10)


     

EA
−
0 −
EA
0
L
L
0 0 0 0
0
0 EA
L
EA
L
0 0 0 0
1
=
~
2
3
4
También podemos ver que una columna j está formada por las reacciones debidas a un
desplazamiento unitario impuesto en la dirección j, y una fila i no es más que las reacciones en
i debido a corrimientos unitarios impuestos en las distintas direcciones.
b) Barra de Pórtico Plano
En base al significado físico de los elementos de la matriz de rigidez, deduciremos la Matriz de
Rigidez para una barra de Pórtico Plano en coordenadas locales.
Para este tipo de elemento corresponden tres desplazamientos por nudo (2 traslaciones y una
rotación 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.

Fig. nº 5 - Desplazamientos unitarios p/ pórtico plano
Las reacciones mostradas en la figura nº5 constituyen las respectivas columnas de la matriz
de rigidez de la barra de Pórtico Plano de la ecuación (11)

       


       


             


             

0 0 EA
0 −
12EI
6EI
6EI
EA
0 12EI
3 2 3 2
0 6EI
4EI
0 6EI
2 2
0 0 EA
0 12EI
6EI
EA
0 12EI
2EI
− − −
6EI
3 2 3 2
0 −
6EI
−
−
−
=

       


       

L
Xj
L
Yj
L
Zj
D
D
D
L
Xk
L
yk
L
Zk
2EI
0 6EI
2 2
L
Xj
F
L
Yj
L
Zj
F
F
L
kj
F
L
yk
F
L
Zk
D
D
D
.
4EI
L
L
L
L
L
L
L
L
0 0
L
L
L
L
L
L
L
L
L
L
0 0
L
L
F
(11)
Escribiendo en forma compacta:
L
~
L
~
L
~
FI = SI .DI (12)
Esta matriz relaciona las fuerzas de extremo de barra con los desplazamientos nodales en
ejes locales.

Fig.nº 6 - Solicitaciones para barra de Pórtico Plano
Cargas nodales equivalentes
Hasta ahora hemos supuesto que las cargas estaban aplicadas en los nudos, y por lo tanto
existe una correspondencia biunívoca entre los puntos de aplicación de las cargas y los
desplazamientos que están siendo calculados. Si esto no ocurriera, por ejemplo tuviéramos
cargas en el tramo de las barras, en forma distribuida o concentrada, debemos sustituir las
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 superposición, que es válido por haber supuesto que el sistema es
lineal, podemos descomponer las cargas tal como se indica en la figura:
Fig.nº 7 - Barra de pórtico con cargas en el tramo
Como podemos observar las cargas, reacciones y deformaciones de la estructura a) serán
equivalentes a la suma de los dos estados b) y c).
Como las deformaciones de nodos en b) son nulas, serán iguales las deformaciones de los
casos c) y a). O sea que las cargas de c) producen la misma respuesta estructural en lo
referente a desplazamientos de nudos que las cargas originales. Estas serán entonces las
cargas equivalentes en los nudos, que no son más que las reacciones de empotramiento
perfecto cambiadas de signo.

Los esfuerzos en los extremos de barra se obtienen por la suma de los casos (b) y (c).
F = F + F c
c
~
b
~
a
~
~
a
~
D = D
Por lo tanto a la ecuación (12) habrá que adicionarle las fuerzas de empotramiento perfecto
del caso b).
L L L L
FI = + AI
(13)
~
.DI
~
SI
~
~
L
~
AI representa el vector de fuerzas de empotramiento perfecto de la barra en coordenadas
locales.
Rotación de ejes en el plano
Hasta el momento expresamos la matriz de rigidez del elemento barra, tanto de reticulado
plano como de pórtico plano, según un sistema de ejes locales, estando los desplazamientos y
esfuerzos de extremo de barra referidos a los mismos.
No obstante, según ya lo mencionamos, antes de poder aplicar las condiciones de equilibrio
en los nudos y de compatibilidad de desplazamientos es necesario transformar las fuerzas y
corrimientos a un sistema de ejes globales.
Fig. nº 8 - Rotación de un vector
Supongamos el vector V de la figura referido al sistema de ejes X e Y.
Las componentes del mismo serán
Vx=V.cos α Vy=V.sen α (14)
En el sistema de ejes xL e yL, las componentes serán
Vx
L=V.cos(α-θ) Vy
L=V.sen(α-θ) (15)

Luego:
Vx
L=V.cos α.cos θ + V.sen α.sen θ (16)
Vy
L=V.sen α.cos θ - V.cos α.sen θ (17)
Teniendo en cuenta la ecuación (14) y escribiendo en forma matricial:



V
.



cos sen

− θ θ

θ θ

= 


X
Y
L
X
V
L
Y
V
sen cos
V
(18)
o en forma compacta:
~ ~
L
~
V = R.V (19)
A la matriz de cosenos directores ~ R
la llamaremos matriz de rotación.
~ R
Si queremos el pasaje del sistema local al global, en este caso la matriz de rotación es la
transpuesta de .



. V





cos sen
θ − θ
θ θ

= 


L
X
L
Y
X
V
Y
V
sen cos
V
(20)
L
~
T
V = R V (21)
~ ~
Premultiplicando la ecuación (19) por R~
−1
L
~
1
V R V = − (22)
~ ~
Lo que implica que R~
es una matriz ortogonal (su inversa es igual a su traspuesta):
T
~
1
R− = R (23)
~
Ecuación fundamental de la barra en el sistema global
Tanto las solicitaciones como los desplazamientos pueden expresarse como vectores en el
plano, podemos entonces aplicar la transformación lineal antes vista para llevar los esfuerzos y
desplazamientos de extremos de barra del sistema local al global.
R.D
~ ~ K
~ D =
L
K
~ ~ J
L
J
=
D . R ~ D
(24)

Donde L
J ~ D
y L
K ~ D
representan respectivamente los desplazamientos de los nudos j y k en
coordenadas locales y
J ~ D
.,
K ~ D
representan los mismos desplazamientos en coordenadas
globales.
Escribiendo en forma compacta:
DIL = (25)
.DI
~
~
RT
~
Siendo:

 


=
 

0 ~ R
~ R
0
~
~
R~T (26)
La matriz de rotación de la barra por incluir ambos nudos de la misma.
L
~
Di
DI contiene los desplazamientos de los dos nudos de la barra en coordenadas locales y ~
los desplazamientos de los mismos en coordenadas globales.
Para las solicitaciones tendremos:
FIL = (27)
.FI
~
~
RT
~
L
~
FI contiene las solicitaciones de los dos extremos de la barra en coordenadas locales y ~
FI
en coordenadas globales.
La inversa de la matriz de rotación de la barra será:

 


= −
 

−
−
1
1
1
0 ~ R
~ R
0
~
~
R~T (28)
Pero como 1
~ R
− = T
~ R
será 1 T
RT
~
RT− =
~
L L L L
FI = + AI
(29)
~
.DI
~
SI
~
~
T L
~
FI = RT
.FI
(30)
~
~
T L L T L
FI = + .AI
(31)
~
RT
~
.DI
~
.SI
~
RT
~
~
T L T L
FI = + .AI
(32)
~
RT
~
.DI
~
~
.RT
.SI
~
RT
~
~
F~I = SI + (33)
.D~I A~I
~

SI = T L .RT
(34)
~
.SI
~
RT
~
~
T L
A~I R~T .A~I = (35)
Matriz de rotación para barra de reticulado plano

   
R~T (36)


   

cos sen 0 0
sen cos 0 0
0 0 cos sen
θ θ
− θ θ
θ θ
− θ θ
=
0 0 sen cos

    
Si .EA/L (37)


    

2 2
c c.s c c.s
2 2
c.s s c.s s
2 2
c c.s c c.s
− −
− −
− −
− −
=
2 2
c.s s c.s s
~
Matriz de rotación para barra de pórtico plano

       
R~T (38)


       

cos sen 0 0 0 0
sen cos 0 0 0 0
0 0 1 0 0 0
0 0 0 cos θ sen θ
0
0 0 0 − sen θ cos θ
0
θ θ
− θ θ
=
0 0 0 0 0 1

Matriz de Rigidez de barra de Pórtico Plano
2 2 3 3 2 2 2 3 3 2
(EAc / L + 12EIs / L ) (EAcs / L − 12EIcs / L ) − 6EIs / L ( − EAc / L − 12EIs / L ) ( − EAcs / L + 12EIcs / L ) −
6EIs / L
3 2 2 3 2 3 2 2 3 2
(EAcs / L − 12EIcs / L ) (EAs / L − 12EIc / L ) 6EIc / L ( − EAcs / L + 12EIcs / L ) (EAs / L −
12EIc / L ) 6EIc / L
2 2 2 2 2 2 2
6EIs / L 6EIc / L 4EI / L 6EIs / L 6EIc / L 2EI / L
−
2 2 3 3 2 2 2 3 3 2
( − EAc / L − 12EIs / L ) ( − EAcs / L + 12EIcs / L ) 6EIs / L (EAc / L + 12EIs / L ) (EAcs / L −
12EIcs / L ) 6EIs / L
3 2 2 3 2 3 2 2 3 2
( − EAcs / L + 12EIcs / L ) ( − EAs / L − 12EIc / L ) − 6EIc / L (EAcs / L − 12EIcs / L ) (EAs / L + 12EIc / L ) −
6EIc / L
6EIs / L 6EIc / L 2EI / L 6EIs / L 6EIc / L 4EI / L
SI
2 2 2 2
~
− −
=
Tabla nº 1

17
Matriz de Rigidez y Vector de Cargas Global de la Estruct.
Fig.º 9 - Reticulado Plano
Las ecuaciones de equilibrio exigen que las cargas externas aplicadas en los nudos deben
ser iguales a la suma de las solicitaciones de extremo de las barras que concurren al nudo.
Siendo el vector de cargas externas aplicadas en j :

   
A j (39)


= −
   

−
−
Xj
A
Yj
A
~
La ecuaciones matriciales de las barras a) y b) son:
a a a a
AI
.DI
~
= +
b b b b
AI
~
.DI
~
~
SI
FI
FI
~
~
~
SI
~
= +
(40)
Que pueden representarse en función de los nudos de la siguiente forma:

  

a
i
A~I
AI

+
  


+

  

a
~ Di
.
Dj

Dj
.
  


  


=
Sii
Sji

=
  



  



  


  

  


  


  


  

  

  

  

b
j
b
k
b
b
a a
Sij
Sjj
b b
Sjk
b b
a
b
b
a
j
a
a a
a
AI
~
AI
~
~
Dk
~
Skk
~
Sjj
Skj
~
~
~
~ Fi
Fj
Fj
~
Fk
~
~
~
~
~
~
~
~
(41)
dónde:
a =
~
Fi vector solicitación en el extremo i de la barra (a)
a =
~
Di vector deformación en el extremo i de la barra (a)
a =
A~Ii vector fuerza de empotramiento perfecto en el extremo i de la barra (a)

18
a =
~
Sij submatriz conteniendo los coeficientes de rigidez del nudo i que provienen de un
desplazamiento unitario del nudo j de la barra (a)
Las condiciones de equilibrio para el reticulado de la figura en el nudo j resultan:
0
−
A j− a − Fj
b =
~
Fj
~
~
a b
Fj
~
−
A j = Fj
+
~
~
(42)
b
j
a b b b b
j
a a a a
AI
~
−
A j = + + + + .Dk
+
~
Sjk
~
.Dj
~
Sjj
~
AI
~
.Dj
~
Sjj
~
.Di
~
Sji
~
~
(43)
Por condición de compatibilidad de desplazamientos:
k
a
b
a
~ c
b
c
j j
j i
i
i k
D
k
~ D
~D ;
~ D
~ D
~D ;
~ D
~ D
~ D
= = = = = = (44)
Reemplazando la (44) en la (43) y extendiendo el razonamiento a los nudos restantes:
c
k
= + + + + +
AI
b
k
.Dk
Sjk
Sik
).Dj
).Di
Sjj
Sii
.D~i (Sjj
(Sii
.Dj
= + + + + +
c b c b
c
i
a
i
Sji
a a c c
b
j
a
j
a a b b
AI
~
AI
~
AI
.Dk
~
Skj
~
).Dk
~
Skk
~
(Skk
~
.Di
~
Sij
Ski
~
−
A j
−
Ai
Ak
~
~
~
~
~
~
~
~
~
~
~
.D~k A~I A~I
~
~
~
~
~
~
= + + + + +
−
(45)
Si la organizamos en forma matricial:

19

    


    

+
AI
AI
A~I +
A~I
+
+

   


~ Di
.
   


    


    

a c a c
Sik
Sij
Sii
a a b b
Sjk
+
+
+
=

     


     

−
A~i
−
A j
~
−
c
k
b
k
b
j
a
j
c
i
a
i
Sjj
Sjj
c b b c
AI
~
AI
~
~
~
Dj
~
Dk
~
Skk
~
Skk
~
Skj
~
Sji
Sii
Ski
~
~
~
~
~
~
~
~
~
A~k
(46)
Que representa el sistema global de ecuaciones de la estructura, escribiendo en forma
compacta:
~ D
.
~ S
A
~ eq
~ A
+ =
−
(47)
Dónde:
−
=
~ A
vector de acciones externas aplicadas en los nudos.
~ S
= matriz de rigidez global de la estructura
~ D
= vector deformación de nudos
=
eq ~ A
= vector de cargas nodales equivalentes (fuerzas de empotramiento perfecto
cambiadas de signo) debido a las cargas aplicadas en los tramos de barras.
Haciendo:
eq ~ A
−
~ A
~ A
= +
(48)
~ D
.
~ S
~ A
= (49)
El vector A representa las cargas nodales equivalentes más las cargas aplicadas en los
nudos.
La ecuación (49) no puede ser resuelta (la matriz S es singular) si no se aplican las
condiciones de contorno o de borde de la estructura.
Esquema de la formación 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 reticulado
tendrá 6 filas y 6 columnas ya que son 3 nudos y cada uno tiene 2 grados de libertad (un
desplazamiento según el eje x y uno según el eje y). La contribución de cada nudo deberá
ser colocada en las posiciones que se indican en la figura nº10.
Los nudos i,j,k definen las filas y columnas de la matriz (submatrices)

20
I j k
a+c a+c a a c c
i
a+c a+c a a c c
a A a+b a+b b b
S = j
a A a+b a+b b b
c C b b b+c b+c
k
c C b b b+c b+c
Fig. nº 10: Esquema Formación Matriz Reticulado
i J k
i
a
ii ~ S
+ c
ii ~ S
a
ij ~ S
c
ik ~ S
S = j
a
ji ~ S
b
jj ~ S
+ a
jj ~ S
b
jk ~ S
k
a
ki ~ S
b
kj ~ S
b
kk ~ S
+ c
kk ~ S
Las submatrices tendrán tantas filas y columnas como grados de libertad tenga cada nudo.
En el caso del reticulado este número es 2.
Analicemos la posición en figura nº 10 de las submatrices de la matriz de rigidez de la barra
b de la figura nº 9.

  


= a
  

S S
~ jj
a
~ ji
a
~ ij
a
a ~ ii
S1
~ S S

  


= b
  

S S
~ kk
b
~ kj
b
~ jk
b
b ~ jj
S1
~ S S

  
S1 (50)


= c
  

S S
~ kk
c
~ ki
c
~ ik
c
c ~ ii
~ S S
Las submatrices son de orden 2x2. La posición b
S1 en la matriz
~
~ S
está indicada en al
~ figura nº 10 a partir del siguiente esquema:
b
= indica fila j columna j
jj S
b
jk ~ S
= indica fila j columna k
b
kj ~ S
= indica fila k columna j
b
kk ~ S
= indica fila k columna k
Nótese que en las posiciones dónde ya existen valores provenientes de otras barras, estos
son sumados a aquellos.

Generalizando, debemos considerar una barra de reticulado que une los nudos j y k según
la figura nº11.
Figura nº 11
Las posiciones ocupadas por la matriz de rigidez de la barra (m) en la matriz de rigidez
global son las siguientes:
21
1
2
3
2(j-1)+1
2(j-1)+2
2(k-)+1
2(k-)+2
1
2
.
.
2(j-1)+1 X X X X
2(j-1)+2 X X X X
~ S
=
.
2(k-1)+1 X X X X
2(k-1)+2 X X X X
Figura nº 12 - Posiciones ocupadas por una barra de Reticulado Plano en la matriz
de Rigidez Global
Para un pórtico plano o estructuras con 3 grados de libertad por nudo tendríamos la
siguiente configuración:

22
1
2
3
3(j-1)+1
3(j-1)+2
3(j-1)+3
3(k-1)+1
3(k-1)+2
3(k-1)+3
1
2
.
.
3(j-1)+1 X X X X X X
3(j-1)+2 X X X X X X
3(j-1)+3 X X X X X X
~ S
=
.
3(k-1)+1 X X X X X X
3(k-1)+2 X X X X X X
3(k-1)+3 X X X X X X
Figura nº 13 - Posiciones ocupadas por una barra de Pórtico Plano en la matriz de
Rigidez Global
Del mismo modo, un vector de cargas equivalentes en los nudos para una barra de la figura
nº11 ocupa la posición indicada en la figura nº 14 del vector acciones globales nodales.
1
2
.
3(j-1)+1 X
3(j-1)+2 X
3(j-1)+3 X
~ A
= .
3(k-1)+1 X
3(k-1)+2 X
3(k-1)+3 X
a) Pórtico plano
1
2
.
.
2(j-1)+1 X
2(j-1)+2 X
~ A
= .
2(k-1)+1 X
2(k-1)+2 X
b) Reticulado plano
Figura nº 14 Posición del vector Acciones nodales

Formación de la MATRIZ de RIGIDEZ GLOBAL en forma de matriz
banda-simétrica.
La Matriz de Rigidez Global de la estructura es simétrica respecto de la diagonal principal y
tiene normalmente características de matriz de bandas. Por ello no es necesario almacenar
todos los elementos de la matriz toda vez que esto implica tiempo de procesamiento para el
computador y además mayor memoria de almacenamiento. Bastará entonces almacenar los
elementos del triángulo superior (o el inferior) de la matriz. También esta parte triangular
superior (o el inferior) puede sufrir una racionalización del arreglo matricial de números. Esto
se corrige aprovechando la característica de banda de la matriz y previniendo su
reorganización para un arreglo bidimensional de dimensión NGLTxB dónde:
NGLT= es el número 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 su forma de
almacenamiento convencional representada en la figura nº 16-a y en la forma de banda
simétrica en la figura nº 16-b.
23
Figura nº 15: Ejemplo de Reticulado Plano.
Figura nº 16: Formas de almacenamiento de la Matriz de Rigidez
Nótese que el semi-ancho de banda está vinculado a la mayor diferencia de numeración entre
nudos que están unidos por una barra. En el caso de la figura nº 15, la mayor diferencia 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. Dónde L es la mayor diferencia entre nudos que

concurren a una misma barra y N1 es el número 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 Numeración de nudos sobre el ancho de Banda
En el ejemplo que sigue se puede apreciar como varían tanto la cantidad de elementos
a almacenar como el semiancho ancho de banda con la forma de numerar los nudos de
la estructura.
A = S .D correspondiente a una estructura completa antes de
24
Condiciones de Contorno o de Borde
Un sistema de ecuaciones :
~ ~ ~
aplicar las condiciones de contorno es INDETERMINADO, pues
~ S
es singular. La razón de
esta singularidad es el resultado de no haber considerado las vinculaciones o apoyos de la
estructura con el exterior. Al introducir las condiciones de vínculo desaparece la
indeterminación (o singularidad de la matriz de rigidez), siempre que el número de vínculos
sea por lo menos el mínimo necesario para eliminar los movimientos de cuerpo rígido de la
estructura.
El conocimiento de determinados corrimientos nodales, disminuye el número de incógnitas,
tornándose innecesarias las ecuaciones correspondientes a estos corrimientos.
La eliminación de la ecuación de un desplazamiento implica la destrucción de la banda de la
matriz, lo cual exigiría un reacomodamiento de las incógnitas del problema.
Al sistema de ecuaciones
A = S .D (51)
~ ~ ~
Podemos expresarlo de la siguiente forma:

25
A φ
~
~ II
S S
~ I,II
~ I,I
~
~ II,I ~ II,II
~ I
~ II
D
.
S S
A
= (52)
Dónde:
I ~ A
= reacciones de apoyo con desplazamiento impedido.
I ~ D
=0, reacciones que por el
momento desconocemos.
II ~ A
= agrupa las fuerzas externas sobre nodos con desplazamientos
II ~ D
desconocidos.
II , II ~ S
= Matriz de rigidez (cuadrada) que relaciona fuerzas conocidas (cargas externas
II ~ A
)
con desplazamientos desconocidos
II ~ D
.
Técnica número 1:
Suprimir o eliminar de la Matriz de Rigidez del sistema las columnas correspondientes a
desplazamientos impedidos y las filas correspondientes a reacciones exteriores.
Ventaja: reduce el tamaño del sistema de ecuaciones a resolver (en un orden determinado por
los grados de libertad eliminados por los apoyos). Dicha reducción es importante en el caso de
pocas barras y muchos vínculos externos, fácilmente aplicable en forma manual, es la técnica
que emplearemos en el curso.
Desventaja: no podemos resolver así desplazamientos prescriptos.
Técnica número 2:
Un artificio usado normalmente para introducir un desplazamiento conocido en la dirección i,
consiste en hacer el elemento de la diagonal principal de la matriz
~ S
de la línea i igual a 1 y
anular todas las restantes posiciones pertenecientes a esa línea y columna. Luego debe ser
colocado el valor del corrimiento conocido i d −
en la posición anteriormente ocupada por ai, y
también pasar al vector
~ A
Sij .d
los coeficientes i
~
−
con el signo cambiado.
S .d S .d .... S .d .... S .d a

   
S .d S .d .... S .d .... S .d a

   
.
.
S .d S .d .... S .d .... S .d a
S .d S .d .... S .d .... S .d a

+ + + + + =
11 1 12 2 1i i 1n n 1
+ + + + + =
21 1 22 2 2i i 2n n 2
+ + + + + =
i1 1 i2 2 ii i in n i
+ + + + + =
n1 1 n2 2 ni i nn n n
.
(53)

26

−
S .d + S .d + .... + 0 + .... + S .d = a −
S .d
11 1 12 2 1n n 1 1i i
     
S .d + S .d + .... + 0 + .... + S .d = a −
S .d
21 1 22 2 2n n 2 2i i
.
.
0 + 0 + .... + d + .... + 0 =
d
i i

S .d S .d .... 0 .... S .d a S .d
n1 1 n2 2 nn n n ni i
.
     

+ + + + + = −
−
− −
−
(54)
En forma matricial
−
a −
S .d
1 1i i
a −
S .d
2 2i i
= (55)
i
.
.
d
.
.
n ni i
1
2
i
d
d
.
.
d
.
.
n
S S ... 0 ... S
11 12 1n
S S ... 0 ... S
21 22 2n
.
...
.
...
.
.
0 0 ... 1 ... 0
.
.
.
.
.
.
n1 n2 nn
a S .d
d
.
S S ... 0 ... S
.
.
.
.
.
.
.
...
.
...
.
.
−
−
−
−
−
Como caso particular, cuando un corrimiento di sea nulo, o sea, cuando exista un vínculo en la
dirección i, la ecuación anterior toma la siguiente forma:
1
2
a
a
.
. 0 .
.
n
1
2
i
d
d
.
.
d
.
.
n
S S ... 0 ... S
11 12 1n
S S ... 0 ... S
21 22 2n
.
...
.
...
.
.
0 0 ... 1 ... 0
.
.
.
.
.
.
n1 n2 nn
a
d
.
S S ... 0 ... S
.
.
.
.
.
.
.
...
.
...
.
.
− = (56)
Es importante mencionar, que con esta forma de aplicación de las condiciones de borde, las
cargas que están aplicadas en la dirección i del apoyo quedarían sin considerar.
Ventaja: No destruye la organización simétrica y banda del sistema de ecuaciones. Es
fácilmente implementable en computador.
Técnica número 3:
Modificar los coeficientes de la diagonal principal en la matriz de rigidez del sistema que
correspondan a desplazamientos impedidos.

27
Siendo una ecuación cualquiera en el sistema:
a S .d S .d d Σ Σ −
a S .d
i ij j
= + ⇒ = (57)
ii
i ij j ii i i S
Luego, haciendo Sii suficientemente grande (por ejemplo multiplicarlo por 1020) obtenemos un
di próximo a cero.
Ventajas: no altera la organización del sistema de ecuaciones, ni tampoco reduce el orden del
sistema.
Desventaja: puede provocar problemas numéricos.
Solución del sistema global de ecuaciones
Una vez introducidas las condiciones de borde o contorno en el sistema de ecuaciones,
podemos resolverlo mediante cualquier método conocido (Gauss, Cholesky, Gauss Jordan,
etc.)
Existen en la literatura métodos eficientes de cálculo teniendo en cuenta la simetría, banda y
forma específica de archivo. La solución de la ecuación 55 nos provee los corrimientos o
desplazamientos de los nudos de la estructura.
Muchas veces existe más de un estado de carga. Para esta situación se pueden desarrollar
algoritmos de solución que trabajando con más de un vector de cargas (término
independiente) nos dan como respuesta más de un vector desplazamiento (uno para cada
estado de cargas).
Cálculo de solicitaciones en extremos de barras
La ecuación matricial de una barra (m) con conexiones i,j se expresa en forma particionada de
acuerdo a la ecuación 41:
m
= + +
~ i
F S .D S .D A
m
~ j j
m
~ jj i
m
~ ji
m
~ j
m
~ i j
m
~ ij i
m
~ ii
F S .D S .D A
= + +
(58)
Los desplazamientos
i ~ D
,
j ~ D son obtenidos del vector desplazamiento de nudos de la
estructura, solución del sistema de ecuaciones (54). Particionada la matriz de rigidez de la
barra “m” (S ,S ,Sm )
~ j,j
m
~ i,j
m
~ i,i
y el vector de empotramiento perfecto (A ,Am )
~ j
m
~ i
, deben ser
almacenadas al confeccionar la matriz de rigidez global y el vector de cargas global.
Entretanto, dependiendo de la capacidad de memoria del procesador, puede ser más
aconsejable recalcular la matriz de rigidez y el vector de reacción de empotramiento perfecto
de cada barra, para después obtener con el término (58) las solicitaciones de extremo de barra
(F yFm )
~ j
m
~ i
.
Las solicitaciones de barras nos interesan referirlas al sistema de ejes de coordenadas locales,
ya que en esta forma tendremos los esfuerzos normales, cortantes, etc. Para esto debemos
hacer una rotación de ejes sobre las solicitaciones obtenidas por (58), para lo cual:

28
m
F m,L
R.F
~ i
F R.F
~ ~ j
m,L
~ j
m
~ ~ i
=
=
. (59)
Esta no es la única forma de obtener las solicitaciones de las barras. Podemos obtener los
corrimientos
i ~ D
y
j ~ D
y transformarlos a coordenadas locales mediante el producto con
~ R
.
Entonces tendríamos las ecuaciones (58), todas en coordenadas locales para así obtener las
solicitaciones en las mismas coordendadas.
Cálculo 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 reacción del nudo.
Σ Σ
= =
R − F = 0 → R =
F
p
m 1
m
~ i ~ i
p
m 1
m
~ i ~ i
Σ=
R =
F (60)
p
m 1
m
~ i ~ i
dónde p es el número de barras ligadas al nudo i y
i ~ R
es la reacción del nudo i.
Si hacemos esta misma suma para un nudo que no sea apoyo, tendremos como resultado las
cargas aplicadas al nudo. Es importante que las solicitaciones del extremo de barra sean
referidas a un sistema global de coordenadas.

29
TOPICOS ESPECIALES EN ANALISIS
MATRICIAL
Efectos Térmicos
Lo analizamos en forma similar al caso de cargas aplicadas en el tramo.
Fig. nº 17
L .
(T T ) 1 2
.L ; (T T )
h
2
j t
1 2
Xj t
−
θ = −α
+
δ = α (61)
El giro dθ de un tramo de viga dx será:
.dx
d (T1 T2 )
h
t
−
θ = −α (62)
teniendo en cuenta que:
M
E.I
d y
dx
d
dx
2
2
= = −
θ
2 −
(T T )
d y 1 2
= −α (63)
h
dx
2 t
integrando dos veces y particularizando para el nudo j.
−
dy (T T )
1 2 2
j t .L
= −α (64)
2.h
Las fuerzas de empotramiento necesarias para restaurar el extremo de la barra a su posición
primitiva serán: (para pórtico plano).
= − (65)
L .
. (T T )
. (T T )
. (T T )
h
L .
2.h
L .
2
.
0 0
6.E.I
3 2
4.E.I
L
E.A
0 12.E.I
0 6.E.I
L
L
L
L
N
Q
Lj
Lj
M
1 2
t
1 2 2
t
1 2
t
2
Lj
−
− α
−
− α
+
α
−

30
P = S .δ (66)
~ j ~ jj ~ j
Cómo estas fuerzas son de empotramiento perfecto, y además están en coordenadas locales,
hay que llevarlas al sistema global y restarlas (cargas nodales equivalentes) al vector de
cargas en los nudos.
.(T T ). E.A
α +
t 1 2
= (67)
.(T T ). E.A
h
0
2
N
Q
Lj
Lj
M
t 1 2
Lj
− α −
Apoyos elásticos
Un apoyo elástico 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 dirección.
Siendo una ecuación cualquiera en el sistema:
Fig. nº 18
ai = ΣSijδj + Sii .δi (68)
Si establecemos ahora el equilibrio en el nudo “i” teniendo en cuenta la fuerza del resorte:
− k.δi + ai = ΣSijδ j + Sii .δi (69)
ai =ΣSijδj + (Sii + k).δi (70)
Es decir debemos sumar la rigidez del resorte k en la diagonal principal de la matriz de rigidez
del sistema en la fila y columna correspondiente a la dirección en la que actúa el resorte.

RT no se altera la igualdad y además considerando que en sistemas de
31
Apoyos deslizantes no concordantes
Fig. nº 19
Si el plano de desplazamiento no es paralelo a uno de los ejes globales, debemos modificar la
ecuación matricial de la estructura, transformando en el nodo en cuestión los vectores fuerza y
desplazamiento al sistema de ejes global.
L
δ = R .δ F = R .F (71)
~ 1 ~ 1 ~ 1
L
~ 1 ~ 1 ~ 1
dónde el “1” indica el número de nodo.
Con
cos β sen β
0
R − β β
sen cos 0
= (72)
0 0 1
δ
~ 2
δ
~ 3
~ 4
L
~ 1 ~ 1
S S 0 0
~ 11 ~ 12 ~ ~
S S S 0
~ 21 ~ 22 ~ 23 ~
0 S S S
~ ~ 32 ~ 33 ~ 34
~ ~ ~ 43 ~ 44
R .F
F
~ 2
F
~ 3
~ 4
L
~ 1 ~ 1
R .
.
0 0 S S
F
δ
δ
= (73)
ij ~ S
=3x3 para pórtico plano
En esta ecuación es lo mismo premultiplicar L
1 ~ δ
por
1 ~ R
que postmultiplicar la primera
columna de
~ S
por
1 ~ R
. Por otra parte, premultiplicando la primera fila a ambos lados de la
igualdad por
~ 1
coordenadas ortogonales
RT .R = I
~ 1 ~ 1 ~

RT y postmultiplicando la columna
32
L
~ 1
δ
δ
= ~ 2
(74)
δ
~ 3
~ 4
R .S .R R .S 0 0
~ 12 ~ ~
T
~ 11 ~ 1 ~ 1
S .R S S 0
~ 21 ~ 1 ~ 22 ~ 23 ~
0 S S S
~ ~ 32 ~ 33 ~ 34
~ ~ ~ 43 ~ 44
T
~ 1
L
F
~ 1
F
~ 2
F
~ 3
~ 4
.
0 0 S S
F
δ
Resumen: Se modifica la ecuación matricial de toda la estructura, premultiplicando la fila,
correspondiente al nodo con apoyo inclinado por
~ 1
correspondiente por
1 ~ R
.
Técnica alternativa
Consiste en adicionar un barra ficticia que deberá impedir desplazamientos en la dirección “v”
y permitir desplazamientos en la dirección “u” y giro. Para que esto sea posible, el área de la
sección transversal de la misma deberá ser mucho mayor que las barras restantes y la inercia
mucho menor a la de las otras barras. (Además de pequeña longitud para evitar acortamientos
de la barra ficticia)
Fig. nº 20

E3 cap4

Más contenido relacionado

La actualidad más candente (18)

Destacado (20)

Similar a E3 cap4 (20)

Más de Ricardo Lujan Jeri (12)

E3 cap4