Análisis Estructural Por El Método de Rigidez en 2D

ANALISIS EST RUCTURAL 2
PAG. 1
Anlisis Estructural por el mtodo de Rigidez en 2D

.
2.0 MATRIZ DE RIGIDEZ DE UN ELEMENTO ESTRUCTURAL

Una estructura en el plano puede subdividirse en elementos estructurales es decir, podemos desensamblarla y
aislarla y por equilibrio esttico podemos decir: si una estructura se encuentra en equilibrio entonces cualquier
parte de ella tambin lo est. El mtodo de rigidez est basado en este principio bsico, por lo tanto, un sistema
estructural es el resultado de un proceso de ensamble de todos sus elementos estructurales y sus cargas
correspondientes, cuyos desplazamientos de todos los elementos con nudos conectados s on compatibles. La
matriz de rigidez de un miembro estructural no vara y es constante solo si la estructura tiene un co mportamiento
elstico. En base a cada matriz de rigidez de los elementos (miembros) estructurales en sistemas de coordenadas
globales permite ensamblar las ecuaciones equilibrio de nudo de la estructura total y del mismo modo se
ensambla la matriz de cargas equivalentes de extremo fijo y las cargas nodales todos en base a equilibrio esttico.
Para ello se encuentran las matrices de rigidez en coordenadas locales del mismo las cargas las mimas que luego se
transforman a un sistema de coordenadas globales .
2.1
FUERZAS INTERNAS DE EXTREMO EN UNA ESTRUCTURA
Una estructura cualquiera puede ser dividida en sub-estructuras mnimas, a las cuales denominaremos miembros
estructurales. Por otro lado una estructura reticular puede estar compuesta por diversos tipos de miembros
estructurales como son: vigas, columnas, muros de cortante, armaduras, con diversas condiciones de extremo, etc.
Estos miembros estructurales estn sometidos a fuerzas internas (esfuerzos) y solo nos interesan para este mtodo
los esfuerzos o fuerzas internas en sus extremos. Estos esfuerzos en sus extremos son desconocidos, sin embargo
estos esfuerzos se pueden expresar en funcin a: las matrices de rigideces y de los desplazamientos que se
producen en los nudos o juntas de la estructura de acuerdo a los Grados de Libertad ( GDL).
Cualquier estructura por el mtodo de rigideces se puede resolver sin limitaciones para ello lo importante es
discretizar en forma adecuada la estructura en elemento estructurales. Asimismo la respuesta de la estructura
ante cargas estticas depende de la rigidez de sus elementos estructurales, los cuales tienen dos nudos en sus
extremos por lo tanto tenemos 3 esfuerzos por nudo y seis por cada elemento estructura l que son incgnitas por
determinar. Para resolver el problema se requieren ecuaciones de equilibrio y se pueden obtener tres ecuaciones
de equilibrio por nudo, cuyas direcciones son dos desplazamientos lineales y un giro. Del mismo modo que las
direcciones consideradas de las fuerzas en los nudos se deben los desplazamientos segn los Grados De libertad y
para un mismo nudo con diversos miembros estructurales que llegan a un mismo nudo rgido, entonces de dichos
desplazamientos son compatibles. Gracias a este concepto el tamao de incgnitas se logra reducir al mnimo. Por
tanto es ms oportuno expresar los esfuerzos en funcin de sus desplazamientos y sus matrices de rigideces
correspondientes. Lo que nos permite reducir el nmero de incgnitas en lugar de considerar a los esfuerzos de
extremo como incgnitas se expresan a estos en funcin de los desplazamientos que resulta un nmero de
incgnitas de esfuerzos.
Para resolver problemas por este mtodo se va desarrollar un procedimiento para obtener las matrices de
rigideces para cada miembro estructural segn diversa condiciones de borde, en principio se hallan dichas
matrices en coordenadas locales y se encuentran sus matrices de transformacin de coordenadas para expresar las
matrices de rigideces, cargas y desplazamientos en un sistema de coordenadas globales, esto es importante para
poder hacer la integracin de toda la estructura llamada ms comnmente ensamblaje. Como se indic
anteriormente sabemos que los desplazamientos en todos los nudos rgidos son compatibles para todos los
elementos que convergen en un mismo nudo. El Ensamble de la estructura se realiza en coordenadas globales
miembro por miembro incluyendo los nudos encontramos la ecuacin matricial de equilibrio y cuya matriz de
rigidez de nudo de toda la estructura se extrae de los coeficientes del vector de desplazamientos y la matriz de
PAG. 2
cargas de extremo fijo y las cargas nodales se ensamblan rpidamente ya que las cargas son datos de entrada. En
base a estas matrices se establece la ecuacin de equilibrio total de la estructura y se ensambla por equilibrio de
nudos, para posteriormente resolver el sistema de ecuaciones cuyas incgnitas son los desplazamientos, una vez
conocidos podemos encontrar para cada miembro estructural los esfuerzos de extremo (ver Figs. 1 y 2).
Fig. 1
Estructura ensamblada (completa) y codificada segn sus miembros estructurales.
Fig. 2 Estructura desensamblada del prtico de la fig. 1.
2.2
FUERZAS RESTAURADORAS DE ELEMENTOS ESTRUCTURALES
Estas fuerzas restauradoras o elsticas de un elemento estructural son dependientes de la rigidez y de los
desplazamientos elemento estructural del material, la geometra y las caractersticas geomtricas de la seccin
transversal de dichos elementos estructurales se originan debido a desplazamientos en sus extremos, es decir, los
desplazamientos (DxJ, DyJ, DJ), que se producen en los nudos inducen fuerzas restauradoras que dependen del
material y la geometra de los miembros estructurales. Se consideran tres grados de libertad por nudo, puesto que
cualquier miembro estructural I simo tiene dos nudos extremos : J y K, por lo tanto se tienen seis grados de libertad
por miembro, ya que existen casos en que varios miembros estructurales tienen nudos comunes , es decir, los
desplazamientos son compatibles en esto casos se simplifican variables por esta razn . Las fuerzas elsticas o
restauradoras en un miembro estructural se infieren debido a los desplazamientos que s ufren sus extremos. Dichas
fuerzas dependen de la rigidez y de los desplazamientos y stas se calculan aplicando el principio de la
superposicin. Las fuerzas restauradoras de extremo se pueden calcular para cada grado de libertad, aplicando
desplazamientos exclusivos para cada grado de libertad, tal como se ilustra en la fig. 3 y posteriormente se calcula el
efecto (o fuerzas) que producen en la direccin de los grados de libertad considerados. Y finalmente el efecto total
se calcula por la superposicin de los seis desplazamientos.
PAG. 3
Fig. 3
Desplazamientos sucesivos para los seis grados de libertad de un miembro estructural en coordenadas
locales.
En un miembro estructural cualquiera se inducen fuerzas internas en sus nudos extremos debido a los
desplazamientos en dichos nudos . Estos esfuerzos en los extremos se denominan vector de fuerzas restauradoras o
elsticas de un miembro estructural con nudos rgidos, cuya estructura matemticas:
fxj

fyj

fj

fxk
f
yk
f
k
AE
L
0
12
EI
3
EI
EI
L
L
AE
L
0
EI
3
L
0
EI
2
EI
12
L
0
AE
L
EI
3
L
2
EI
L
EI
3
0
12
EI
3
L
0
EI
6
d
2 xj
L
d
yj
EI
2
L d
j

d
0 xk
d
EI yk
6
3 d k
L
EI
4
2
L
0
12
L
4
EI
L
0
AE
EI
3
Esta relacin matricial puede particionarse en sub-matrices, segn sus nudos extremos: J y K. Donde los
smbolos y nomenclatura usada esta con letra minscula indica que el sistema de coordenadas es local y si el sistema
de coordenadas es global el texto esta expresado con letra mayscula. Superpuesto sobre los vectores de:
desplazamientos, fuerzas elsticas y la matriz de rigidez, nos indica que estn en un sistema de coordenadas es
locales. Y el primer sub-ndice de las fuerzas o desplazamientos representa la direccin, que puede ser: "x", "y", "".
Y el segundo sub-ndice se refiere al cdigo del nudo. El vector de fuerzas restauradoras expresado en sub -matrices
es:
(2)
Donde:
{IfJ}, {IfK} = fuerzas internas en los extremos de las barras nudas: J o K en coordenadas locales.
{Id J}, {Id K} = desplazamientos en los nudos J o K en coordenadas locales.
[Ib ]
= sub-matriz de 3x3 de rigidez del miembro I-simo en coordenadas locales.
PAG. 4
La ecuacin (2) en forma compacta, puede escribirse as:

{If} = [Ik] {Id}
(3)
Donde:
{If} = vector de fuerzas elsticas en coordenadas locales del miembro i.
{Id} = vector de desplazamientos en coordenadas locales del miembro i.
[Ik] = matriz de rigidez en coordenadas locales del miembro i.
Fig. 4 Codificacin de fuerzas y desplazamientos en coordenadas locales de un miembro estructural I-simo.
La matriz de rigidez de un miembro estructural I-simo en coordenadas locales expresada en sub-matrices

en funcin de los nudos conectivos J y K, es:
[Ik] =
[Ib JJ] [Ib JK]

--------- ---------[Ib KJ]
2.3
(4)
[Ib KK]
MATRICES DE RIGIDEZ EN COORDENADAS LOCALES
La matriz de rigidez de un miembro en coordenadas locales, se obtiene, del coeficiente del vector de
desplazamientos, al producir desplazamientos unitarios para cada grado de libertad, tal como se ilustra en la fig. 2.
Para prticos planos en el caso ms general se consideran tres grados de libertad por nudo, a saber: dos
desplazamientos lineales en la direccin de los ejes x, y; y un desplazamiento angular en la direccin .
PAG. 5
2.3.1 Matriz de Rigidez de un miembro estructural, con los nudos extremos J y K

rgidos: (Tipo de Miembro: MT=0)
AE
L
0
EI
12
EI
EI
L
0
EI
3
L
0
EI
12
EI
L
AE
EI
3
EI
L
EI
3
L
0
12
L
2
L
0
EI
3
L
0
EI
6
2
L
EI
2
L
EI
6
3
L
EI
4
2
L
0
12
L
AE
EI
L
0
AE
EI
3
Esta matriz se ensambla para cada grado de libertad aplicando el principio de superposicin se encuentra el
efecto total. Para una mejor ilustracin se indican los grados de libertad en el encabezamiento de la matriz de rigidez
(ver Fig. 5).
Fig. 5
Desplazamientos unitarios para cada grado de libertad por separado y el efecto que produce en sus extremos, para un
miembro con los nudos J y K rgidos.
2.3.2 Matriz de Rigidez para un miembro con articulacin en el nudo J y rgido en el

nudo K. La matriz de rigidez se obtiene por el principio de superposicin cuya estructura
matemtica se muestra a continuacin (Ver Fig. 6): (MT=1)
AE
L
0
EI
AE
L
L
0
AE
L
0
EI
3
AE
L
L
0
EI
2
EI
EI
3
3
3
2
L
L
0
0
0
0
EI
EI
3
3
3
2
L
L
EI
EI
3
3
2
L
L
PAG. 6
Fig. 6
Desplazamientos unitarios para cada grado de libertad y el efecto que producen sus los extremos, para un miembro
con el nudo J articulado y el nudo K rgido.
2.3.3 Matriz de Rigidez para un miembro estructural, con el nudo J rgido y con
articulacin el nudo K articulado (Ver Fig. 7), (MT=2)
AE
L
0
EI
EI
L
0
EI
2
L
0
0
3
EI
EI
L
EI
2
L
0
EI
EI
L
AE
L
0
L
AE
AE
EI
3
L
0
Fig. 7
Desplazamientos unitarios para cada grado de libertad por separado y el efecto que produce en los extremos, para un
miembro con el nudo J rgido y con articulacin en el nudo K.
PAG. 7
2.3.4 Matriz de Rigidez de un miembro tipo armadura, (articulado en los nudos: J y K) aplicando el
principio de superposicin se obtiene la matriz de rigidez de miembro. Alternativamente tambin puede
obtenerse haciendo EI=0, a la matriz de la ecuacin (5), (MT=3):
AE
L
0 0
0 0
0 0
AE
0 0
0 0
0 0
AE
0 0
AE
L
0 0
0 0
0 0
0 0
0 0
Fig. 8 Desplazamientos unitarios en un miembro tipo armadura y el efecto que produce con los nudos: J y K articulados.
PAG. 8
2.3.5 Matriz de Rigidez considerando la deformacin por esfuerzo cortante. - esto es,
para miembros tipo: muros de cortante. Esta matriz es usada para prticos que incluyen muros de cortante (ShearWall); y tambin viga de corte su estructura matemtica en coordenadas locales, es: (Ver Fig. 9, siendo el nmero de
cdigo MT=4)).
Fig. 9
Miembro estructural considerando la deformacin por esfuerzos: axial, flexin y corte, (viga de Timoshenko).
s1
s 1
0
2
3
0
s
s
3
5
0
2
s
s
s
3
s
4
s
3
s
5
0
2
3
2
3
Cuyos parmetros son:
1 =
4 =
(2 )
(1 + )
12
2
2 =
12
(1 + )3
5 =
(4 + )
(1 + )
3 =
6
(1 + )2
Donde:
E
G
v
Ff
=
=
=
=
Mdulo de elasticidad del material.

E/(2(1+v) (Mdulo de rigidez).
Relacin de Poisson (v=1/6 para el concreto y 0.3 para el acero).
Factor de forma por corte que depende de la seccin transversal.
Arw = rea Reducida por deformacin por corte.
PAG. 9
2.3.6 Matriz de Rigidez con zonas rgidas considerando la deformacin por

esfuerzos: cortante, axial y flexin (este caso se presenta en los prticos con muros de cortante y/o en vigas de
corte, ver Fig. 9) Cdigo: Member Type: MT=5
cL
aL
bL
Fig. 10 Miembro estructural que contempla a un elemento con extremos rgidos y zonas flexibles.
1
0
0
Tf
0
0
0 0
0 0
d L 1 0
0 1
0 0
0 0 0 b L 1
1
0
0
Td
0
0
0 0
1 d L 0 0
0
0 0
1 0
0 1
0 0
0
0
b L
1
0
AE
AE
0
0
0
0
c L
c L
12 E I
6 E I
12 E I
6 E I
0
0
3
2
3
2
( c L)
E I ( c L)
E I ( c L)
E I ( c L)
6 E I
( 4 ) E I
6 E I
2
0
0
2
E I c L
2
E I c L
E I ( c L)
E I ( c L)
1
KLw
1 AE
AE
0
0
0
0
c L
c L
12
I
6
I
12
I
6
I
0
3
2
3
2
E I ( c L)
E I ( c L)
E I ( c L)
E I ( c L)
6 E I
2
6 E I
( 4 ) E I
0
0
2
E I c L
2
E I c L
E I ( c L)
E I ( c L)
AE
AE
0
0
0
0
L c ( 1)
L c ( 1)
12
I
6
12
6
0
0
3 3
2 2
3 3
2 2
L c ( 1)
L c ( 1)
L c ( 1)
L c ( 1)
6
12 E I d
6 d
4
6
12 d
6 d
2
2 2
2 3
2
L c ( 1)
2 2
2 3
2
E I L c ( 1)
L c ( 1)
L c ( 1)
L c ( 1)
L c ( 1)
L c ( 1) L c ( 1)
Tf KLw
AE
AE
0
0
0
0
L c ( 1)
L c ( 1)
12
6
12
6
3
3
2
2
3
3
2
2
L c ( 1)
L c ( 1)
L c ( 1)
L c ( 1)
6
12 b
6 b
2
6
12 b
6 b
4
2 2
2 3
2
E I L c ( 1)
2 2
2 3
2
L c ( 1)
L
1
)
L
1
)
L
1
)
L
1
)
L
1
)
L
1
)
3.
PAG. 10
CARGAS EQUIVALENTES DE EXTREMO FIJO
3.1 ESFUERZOS DE EXTREMO, CARGAS EQUIVALENTES DE EXTREMO Y

CARGAS NODALES
Todas las estructuras estn conformadas por elementos o miembros o barras estructurales que en un caso
general pueden estar sometidas a cargas aplicadas en los nudos (denominadas carga nodal), y cargas que actan
sobre los miembros denominados por el mtodo de rigidez cargas equivalentes de extremo o de empotramiento
perfecto.
El mtodo de la rigidez se basa en una formulacin generalizada de la matriz de rigidez de miembro en
coordenadas globales para los desplazamientos en los nudos y en la direccin que sea posible describir su
deformacin y ello se demuestra con el procedimiento de flexibilidades o por el mtodo de energa de deformacin
que nos sirve para cualquier tipo de cargas. Por ello este mtodo trabaja con los desplazamientos que sufren los
nudos y que son consideradas como incgnitas llamados grados de libertad (GDL). Puesto que debemos establecer
el equilibrio en los nudos de la estructura, por ello las cargas que actan sobre el miembro 2 de la Fig. 11 por
ejemplo deben transformarse a cargas equivalentes de extremo.
Fig. 11
Cargas en un nudo sobre un miembro de una estructura.
Las fuerzas internas de extremo de los miembros en coordenadas locales { If}se calculan en base a su matriz
de rigidez [Ik] por el desplazamiento de los nudos extremos { Id} y agregando las cargas equivalentes de extremo
fijo todos en coordenadas locales {IfF}, se calculan en base a las cargas en coordenadas locales. Por tanto el vector
de fuerzas de extremo est dado por: (Ver Fig.12).
{If} = [Ik] {Id} + {IfF}
(11)
PAG. 11
Y
K
Y
J
X
Fig. 12
Cargas equivalentes de extremo fijo en un sistema de coordenadas locales.
Donde:
{IfF} = Vector de cargas equivalentes de extremo fijo en coordenadas locales del miembro I-simo.
J
X
Fig. 13
Fuerzas de extremo y desplazamientos de un miembro en un sistema de coordenadas locales.
4.
PAG. 12
TRANSFORMACION DE COORDENADAS de FUERZAS Y DESPLAZAMIENTOS
fy
fxk.Sen
fyk.Cos
Un miembro estructural puede estar orientado en cualquier direccin. Se sabe tambin que los vectores de
fuerzas y desplazamientos tienen los mismos cosenos directores. Por tanto basta con calcular los cosenos directores
del miembro geomtricamente y mediante una transformacin de coordenadas por rotacin de dichas coordenadas
locales obtenemos sus componentes en un sistema de coordenadas globales. (Ver Fig. 14).
xm
fyk.Sen
YK
k
fx
fxk.Cos
fy
YJ
k
f
fxj.Sen
fyj.Cos
j
fx
fxj.Cos
fyj.Sen
j
i
f
XJ
XK
Fig. 14a Transformacin de fuerzas de extremo de un sistema de coordenadas de locales a globales.
de la Fig. 14a en base a las coordenadas calculamos la longitud del miembro I-simo, por Pitgoras:
(xK - xJ)2 + (y K - y J)2
LI =
y sus cosenos directores:
(12)
cx = cos I = (xK - xJ)/ LI
(13)
cy = sen I = (y K - y J)/ LI
(14)
4.1 Transformacin de coordenadas del vector de fuerzas de un nudo cualquiera.- De la Fig. 14 para
un nudo cualquiera se obtienen las relaciones de transformacin de coordenadas de fuerzas en coordenadas locales a
globales de una barra (mediante el lgebra vectorial de esttica se logra esta transformacin):
Minscula -> local
Mayscula-> global
Fx = fxcx - fy cy + 0
Fy = fxcy + fy cx + 0
(15)
F = 0
+ 0
+ f
Escribiendo estas ecuaciones en forma matricial, tenemos:
Fx
Fy =
F
cx
cy
0
-cy
cx
0
0
0
1
fx
fy
f
(16)
y en forma compacta la fuerza en un nudo se expresa con la submatriz de 3x3:
{F} = [A]T {f}
(17)
PAG. 13
4.2 Transformacin de coordenadas del vector de fuerzas de un miembro I-simo,

que tiene los nudos conectivos: J y K, este vector se forma completando un arreglo matricial quedando la siguiente
expresin:
{IFJ}
------{IFK}
[IA]T {IfJ} + [0] {IfK}

-------------------------------[0] {IfJ} + [IA]T {IfK}
(18)
Factorizando los vectores de fuerzas locales de los nudos J y K, se obtiene:

[IA]T
{IFJ}
[0]
{IfJ}
=
{IFK}
(19)
[IA]T
[0]
{IfK}
Expresando en forma compacta:

{IF} = [IR]T {IfJ}
(20)
en donde la matriz [IR], representa a la matriz de transformacin de coordenadas, cuya matriz tiene la
propiedad de ser una matriz ortogonal, esto sucede cuando su determinante es igual a la unidad y tambin
cuando la inversa de dicha matriz es igual a su transpuesta, es decir:
[IR]T = [IR]-1
(21)
de ello, podemos encontrar la siguiente transformacin coordenadas globales a locales :

{If} = [IR] {IF}
(22)
Donde:
[IR] =
cx
-cy
0
cy
cx
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
cx
-cy
0
cy
cx
0
0
0
1
(23)
PAG. 14
xm
F yK
DyK
DxK
YK
F xK
DyJ
F yJ
ym
DxJ
YJ
F xJ
X
XJ
XK
Fig. 14b Fuerzas de extremo y desplazamientos de nudo no dependen del miembro en un sistema de coordenadas globales.
4.3
Transformacin de coordenadas del vector de desplazamientos
para un miembro I-simo, anlogamente al procedimiento descrito en la seccin 4.1.2, se obtiene:

{ID} = [IR]T {Id}
(24)
Despejando tenemos:
{Id} = [IR] {ID}
(25)
4.4 Transformacin de coordenadas locales del vector de fuerzas de empotramiento perfecto

a un sistema de coordenadas globales de un miembro I-simo:
{IFF} = [IR]T {IfF}
(26)
Donde:
{IfF} = Vector de cargas equivalentes de extremo fijo en coordenadas locales.
{IFF} = Vector de cargas equivalentes de extremo fijo en coordenadas globales.
4.5
Transformacin del sistema de coordenadas locales del vector de fuerzas elsticas del miembro I-simo a
coordenadas globales. Sabemos que el vector de fuerzas restauradoras en coordenadas locales se encuentra con la
expresin matricial (ecs-1 y 4):
{If} = [Ik] {Id}
(4)
Reemplazando el vector de fuerzas restauradoras a un sistema de coordenadas globales, se tiene

{IF} = [IR]T {If} = [IR]T [Ik] {Id}
(27)
transformando el vector de desplazamientos a un sistema de coordenadas globales, la ecuacin (27)

anterior, queda:
{IF} = [IR]T [Ik] {Id} = [IR]T [Ik] [IR] {I D}
(28)
PAG. 15
de donde la matriz coeficiente del vector { ID} es la matriz de rigidez del miembro I-simo en coordenadas
globales, esto es:
[IK] = [IR]T [Ik] [IR]
4.6
(29)
Transformacin de coordenadas del vector de fuerzas de extremo de un miembro I
en coordenadas locales a un sistema de coordenadas globales, la (ec-11) se encuentra con la siguiente expresin:
{If} = [Ik] {Id} + [IfF]
(11)
transformando el sistema de coordenadas de este vector de fuerzas a un sistema de coordenadas globales,

reemplazando en la ec-11 en la ec-20 se tiene:
{IF} = [IR]T [Ik] {Id} + {IfF}
(30)
efectuando el producto matricial, obtenemos:
{IF} = [IR]T [Ik] {Id} + [IR]T {IfF}
(31)
Reemplazando por: {Id} = [IR] [ID] en la ecuacin anterior, haciendo un arreglo, se tiene:
{IF} = ( [IR] T [Ik] [IR] ) {ID} + [IR]T {IfF}
(32)
en donde la matriz coeficiente del vector de desplazamiento { ID} representa a la matriz de rigidez de un
miembro en coordenadas globales, esto es:
(29)
y el vector de cargas de extremo equivalente en coordenadas globales es:

[IFF] = [IR]T [IfF]
de ello la ecuacin (32) puede ser expresada as:
{IF} = [IK] {ID} + [IFF]
(32)
4.7
PAG. 16
Aplicacin de la transformacin.
Ejemplo 1. Encontrar la matriz de rigidez de una columna para un sistema de coordenadas globales,
sabemos que la direccin de una columna es = 90; y la matriz de rigidez en coordenadas globales se encuentra
con la transformacin:
Y la matriz de rigidez de un miembro estructural en coordenadas locales puede expresarse, as:
[IR] =
a1
0
0
0
a2
a3
0
a3
2a4
-a1
0
0
0
-a2
-a3
0
a3
a4
-a1
0
0
0
-a2
a3
0
-a3
a4
a1
0
0
0
a2
-a3
0
-a3
2a4
donde:
a 2 = 12EI/L3
a1 = AE/L,
a3 = 6EI/L2
a4 = 2EI/L
Y la matriz de transformacin dada en la ecuacin (23), para: =90 tenemos que, cx=0, cy=1, ser,
reemplazando, se obtiene:
[IR] =
0
-1
0
1
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-1
0
1
0
0
0
0
1
Efectuando el triple producto matricial, obtenemos la matriz de rigidez de miembro en coordenadas globales:
[IK] = [IR]T [IK] [IR] =
a2
0
-a3
-a2
0
-a3
0
a1
0
0
-a1
0
-a3
0
2a4
a3
0
a4
-a2
0
a3
a2
0
a3
0
-a1
0
0
a1
0
-a3
0
a4
a3
0
2a4
PAG. 17
5.
EQUILIBRIO DE LA ESTRUCTURA
5.1
ECUACIONES DE EQUILIBRIO DE NUDOS DE LA ESTRUCTURA
Una estructura plana de NJ nudos que est en equilibrio esttico, tiene tres ecuaciones de equilibrio por
nudo, por lo tanto tenemos: 3 NJ ecuaciones de equilibrio de la estructura. Y en cada miembro estructural actan
seis fuerzas, tres por cada nudo. La magnitud de dichas fuerzas depende de las rigideces de los miembros
estructurales, de las cargas en los nudos y de las cargas equivalentes de extremo fijo. En consecuencia, puesto que
la estructura est en equilibrio, la suma de las fuerzas en un nudo debe ser igual a cero .
Con estas ecuaciones de nudo se ensambla la matriz de rigidez del nudo total y los coeficientes en sub matrices del vector de desplazamientos se suman.
Sea S un nudo S-simo cualquiera, al ser aislado de la estructura dicho nudo tambin debe estar en
equilibrio, es decir, la resultante de la sumatoria de las fuerzas totales que llegan a este nudo; transmitidas por los
extremos de los miembros en sentido contrario ms las cargas nodales debe ser igual a cero.
Fig. 15 Equilibrio de la estructura se da en todos sus nudos.
El vector de fuerzas de extremo de un miembro en coordenadas globales, expresado en sub -matrices, es:
{IFJ}
------- =
{IFJ}
[IBJJ] [IBJK ]
---------------[IBKJ] [IBKK ]
{DJ}
{IFFJ}
------- + ------{DK }
{IFFK }
(33)
Es importante notar que el desplazamiento de los nudos no depende del sub ndice I que representa el
cdigo de miembro debido a ello se ha obviado este sub ndice puesto que puede existir varios miembros con un
mismo nudo en comn. Luego efectuando el producto matricial en sub-matrices de la ecuacin (33) se obtienen los
vectores de fuerzas de extremo en los nudos J y K:
{IFJ} = [IBJJ] {DJ} + [IBJK ] {DK } + {IFFJ}
(34)
{IFK } = [IBKJ] {DJ} + [IBKK ] {DK } + {IFFK }
(35)
Con estos vectores de fuerzas de extremo por nudo se establecen las ecuaciones de equilibrio de un nudo cualquiera
ensimo, es decir: (n = J K)
PAG. 18
NB
NB
- {IFn } + {FN n } = {0}

n=1
(36)
n=1
donde:
{FN n }
= Vector fuerzas del nudo ensimo (3x1).
{IFJ}, {IFK } = Vectores de fuerzas de extremos de los nudos: J y K, del miembro I-simo. (3x1)
NB
= Nmero de miembros conectivos al nudo ensimo.
Con la ecuacin (34), se ensambla la matriz de rigidez de nudo total [K] separando los coeficientes de
rigidez, cargas y desplazamientos mediante un arreglo matricial, se obtiene la ecu acin matricial de equilibrio de
nudo:
-[K] {D} - {FF} + {FN } = {0}
(37)
donde:
[K]
{D}
{FF}
{FN }
=
=
=
=
matriz de rigidez de total de la estructura.

vector de desplazamientos en los GDL.
vector de cargas equivalentes de extremo fijo de toda la estructura.
vector de fuerzas nodal.
La ecuacin matricial de equilibrio de nudo total, puede particionarse en funcin de los grados de libertad.
Ya que los vectores de desplazamientos en los GDL son desconocidos y sabemos que los desplazamientos de los
nudos restringidos son nulos. Por ello, planteamos las ecuaciones de equilibrio en sub -matrices en funcin de los
desplazamientos segn los GDL y los grados restringidos, as:
[KUU ] [KUR]
- ----------------[KRU ] [KRR]
{DU }
------{DR}
{FFU }
------- +
{FFR}
{FN U }
------{FN R}
{0}
----{0}
(38)
donde:
{DU }
{DR}
{FFU }
{FFR}
{FN U }
{FN R}
=
=
=
=
=
=
vector de desplazamientos de los nudos libres.

vector de desplazamientos de los nudos restringidos.
vector de fuerzas de empotramientos de los nudos libres.
vector de fuerzas de empotramientos de los nudos restringidos.
vector de cargas en los nudos libres.
vector de cargas en los nudos restringidos o vector de reacciones.
Los grados restringidos con subndice u son desplazamientos conocidos , es decir, son nulos,
reemplazando {DR} = {0} efectuando el producto en sub matrices:
- [KUU ] {DU }
{FFU } +
{FN U } =
{0}
(39)
- [KRU ] {DU }
{FFR}
{FN R}
{0}
(40)
de la ecuacin (39) despejamos el vector de desplazamientos de nudos libres:
{DU }
[KUU ]-1 (-{FFU } + {FN U } )
y el vector de reacciones o de fuerzas restringidas se encuentra con la ecuacin (40):
(41)
{FN R}
PAG. 19
= [KRU ] {DU }
{FFR}
(42)
EJEMPLO 1
GEOMETIA Y CARGAS
Para el prtico que se muestra en la figura encontrar lo siguiente:
1. Matriz de rigidez global.
2. Desplazamientos y giros en cada nudo.
3. Momentos flectores, fuerzas cortantes y fuerzas axiales.
5
2
1
2.5Tn/m
14
13
3.00 m
20
21
23
19
24
17
16
5Tn
3.00 m
15
10
18
10Tn
11
12
10
2.5Tn/m
8
4
22
3.00 m
3.00 m
3.00 m
3.00 m
Fig. 16 Prtico codificado en nudos y miembros y las cargas a la derecha.
DATOS GENERALES:
Nmero de nudos = 8, Nmero de barras = 10
BLOQUE DE DATOS DE NUDOS:

Coordenadas
X
Y
0
6
3
6
6
6
0
3
3
3
0
0
3
0
6
0
NUDO
1
2
3
4
5
6
7
8
Restricciones
Irx
Iry
0
0
0
0
0
0
0
0
0
0
1
1
1
1
1
1
BLOQUE DE DATOS DE MIEMBROS:
Ir
0
0
0
0
0
1
1
1
I = 9.232E-5 A = 9.6E-3 E = 2E6 L = 3
Miembro
I
JI
KI
Tipo De Elemento
II
AI
LI
EI
Cx I
CyI
1
2
3
4
5
6
7
8
6
7
4
5
8
4
4
1
4
5
1
2
3
7
5
5
con extremos empotrados

con extremos articulados
I
I
I
I
I
I
I
I
A
A
A
A
A
A
A
A
L
L
L
L
2L
2L
L
2L
E
E
E
E
E
E
E
E
0
0
0
0
0
1/
1
1
1
1
1
1
1
-1/
0
-1/
2
2
9
10
1
2
2
3
PAG. 20

A
A
I
I
L
L
E
E
1
1
DATOS DE CARGAS:
Cargas nodales:
Nudo 1:
Fx = 10Tn
Nudo 4:
Fx = 5Tn
Cargas sobre los miembros:

Miembro 7:
y = -2.5Tn/m2
Miembro 9:
y = -2.5Tn/m2
Miembro 10:
y = -2.5Tn/m2
SOLUCION
Encontramos las matrices de rigidez locales [k], de cada miembro con las siguientes frmulas:
k E
A
L
0 0
A
L
0 0 0 0 0
0 0 0 0 0
A
L
0 0
A
L
0 0 0 0 0
0 0 0 0 0
0
0
k E
M atriz de rigidez para armaduras
A
L
12I
6I
4I
6I
2I
3200
0
0
KL
5 3200
0
6400
82.062
123.093
82.062
123.093
3200
10.258 30.773
10.258
30.773 123.093
30.773
3200
10.258 30.773
10.258
30.773 61.547
30.773
6I
0
12I
L
6I
L
Matriz de rigidez local para el miembro 5:
123.093 246.187
0
123.093 123.093
0
0
6400
0
0
82.062 123.093 0
82.062 123.093
123.093 123.093
0
123.093 246.187
0
6I
2
L
2I
6I
2
L
4I
L
M atriz de rigidez para miembros con extremos empotrados
Matriz de rigidez local para los miembros 1, 2,3 y 4:
6400
0
k= 0
KL1
6400
0
6I
12I
L
A
L
12I
6I
L
A
L
A
L
k=
30.773
61.547
0
30.773
123.093
0
Matriz de rigidez local para los miembros 6 y 8:
Matriz de rigidez local para los miembros 7, 9 y 10:
0
0
4525.483
k=
0
KL6
4525.483
PAG. 21
0 0 4525.483 0 0
0 0
0 0
0 0
0 0 4525.483
0 0
0 0
0
0
6400
0
k= 0
KL7
6400
0
6400
82.044
123.067
82.044
123.067 246.133
123.067
6400
82.044 123.067
82.044
123.067 123.067
123.067
123.067
123.067
123.067
246.133
0
Encontramos las matrices transformacin de coordenadas [R], para cada miembro con la siguiente frmula:
Cx
C y
0
IR
0
0
1
0
IR
0
0 0
Cx 0 0 0 0
0 1 0 0 0
0 0 Cx Cy 0
0 0 C y C x 0
0 0 0 0 1
Cy
0 0
Donde:
Cx, Cy: son los cosenos directores
L
Cx
0 1 0 0 0
0 0 0 0 0
1 0 0 0 0
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
IR
0 0 0 0 1
0 0 1 0 0
0 0 0 1 0
0 0 0 1 0
0 0 0 0 1
0 0 0 0 0
0.7071 0.7071 0
0
0
0
0.7071 0.7071 0
0
0
0
0
0
1
0
0
0
IR
0
0
0 0.7071 0.7071 0
0
0 0.7071 0.7071 0
0
0
0
0
0
1
0
Matriz de transformacin
de coordenadas para
miembros horizontales
(7, 9, 10)
XK XJ2 YK YJ2
XK XJ
L
Cy
YK YJ
L
0
0
de coordenadas para
miembros verticales
(1, 2, 3, 4, 5)
de coordenadas para
miembros inclinados
(6, 8)
En base a estos datos podemos encontrar las matrices de rigidez de los miembros en coordenadas locales, para
luego encontrar en la matriz de rigidez del elemento en coordenadas globales; con la siguiente frmula:
[IK] = [IR] T [Ik] [IR]
matriz de rigidez global de uno de los elementos.
MIEMBROS 1, 2, 3 y 4:
[RT]
[K] = [1,2,3,4 R] T [1,2,3,4 k] [1,2,3,4 R]

[k]
[R]
6400
0
0
KL1
6400
0
PAG. 22
123.093 246.187
0
123.093 123.093
0
0
6400
0
0
82.062 123.093 0
82.062 123.093
123.093 123.093
0
123.093 246.187
0
6400
82.062
123.093
82.062
123.093
En el MathCad se utiliza otra nomenclatura como la que sigue
0
1
0
R
1 0
0
1 0 0 0 0
0 0 0 0
0 1 0 0
0 0 0 1
0 0 1 0
0 0 0 0
0
0
6400
0
0
KL
1 6400
0
82.062 123.093
0
82.062 123.093
123.093 246.187
0
123.093 123.093
0
0
6400
0
0
82.062 123.093 0
82.062 123.093
123.093 123.093
0
123.093 246.187
0
6400
En el MathCad
0
82.062
0
6400
123.093 0
T
R KL R
1
1 1 82.062
0
0
6400
123.093 0
123.093 82.062
123.093
6400
246.187 123.093 0
123.093
123.093 82.062
0
123.093
0
0
6400
0
123.093 123.093 0
246.187
0
MIEMBRO 5:
[5 K] = [5 R] T [5 k] [5 R]
MIEMBROS: 6 y 8:
[6,8 K] = [6,8 R] T [6,8 k] [6,8 R]
PAG. 23
[7,9,10 K] = [7,9,10 R] T [7,9,10 k] [7,9,10 R]
MIEMBROS 7, 9 y 10:
Ahora procederemos a encontrar las matrices de Cargas Equivalentes de Extremo Fijo para cada miembro:
If
If
{If
{If
=
}
If
If X
If Y
If
F
F
f X
F
F
f Y
f X
f
f Y
Donde el sub-ndice I denota el cdigo del miembro y los sub-ndices J, K denotan los nudos conectivos de
dicho miembro (cercano y lejano).
Momentos de empotramiento perfecto:
Reacciones de empotramiento perfecto:
.L2
12
.L
2
2.5 x 32
1.875Tn.m
12
2.5 x 3
3.75Tn
2
Para el ejemplo, slo los miembros 7, 9, y 10 poseen carga distribuida y su matriz de fuerzas de empotramiento
perfecto es la siguiente:
0 X
3.75 Y
Para este caso la matriz de fuerzas es la misma para los
1.875
F
elementos 7, 9 y 10, al tener similar geometra,

f
0 X
dimensiones y misma condicin de carga.
Y
3.75
1.875
En base a estos resultados podemos encontrar la matriz de fuerzas de empotramiento perfecto de cada miembro en
coordenadas globales, con la siguiente expresin:
PAG. 24
Donde:
[IFF] : Matriz de empotramiento en coordenadas globales
[IR] T: : Matriz de transformacin de coordenadas
[IfF] : Matriz de cargas equivalentes de extremo fijo en coord. locales
[IFF] = [IR] T [IfF]
Desarrollando el producto matricial obtenemos:
3.75
1.875
F
F
0
3.75
1.875
PROCEDIMIENTO DE ENSAMBLE DE LAS MATRICES DE RIGIDEZ, FUERZAS Y

DESPLAZAMIENTO
Ahora procederemos a resolver la siguiente ecuacin en forma matricial.
{IF}= [IK] {ID}+ {IFF}
Donde:
[IFF]
{ID}
[IK]
{IF}
:
:
:
:
Matriz de fuerzas de empotramiento en coordenadas globales.

Matriz de desplazamientos globales.
Matriz de rigidez de miembro en coordenadas globales.
Matriz de fuerzas del miembro en coordenadas globales.
Por simplicidad del problema la ecuacin anterior la trabajaremos con sub-matrices, todas las matrices se encuentran
en coordenadas globales, esta ecuacin se aplica para cada miembro.
{IFJ }
{IFK }
[IKJJ]
[IKJK]
{IDJ }
[IKKJ] [IKKK]
{IDK }
{IFJ }
+
{IFK }
Donde el sub-ndice I denota el cdigo del miembro y los sub-ndices J, K denotan los nudos conectivos de dicho
miembro.
1 paso.- Se preparan los casilleros de la matriz de rigidez de NJ x NJ, en este ejemplo la matriz de rigidez es de
8x8, la dividimos en sub-matrices de 3x3, para cada nudo:
PAG. 25
U
1
GRADOS DE LIBERTAD
R
5
1
2
U
3
4
5
6
R
7
8
NJx1
NJxNJ
NJx1
NJx1
Fig. 17 Casilleros para el ensamble de la matriz de rigidez del nudo de la estructura

Donde:
U = Grados de libertad libres.
R = Grados de libertad restringidos.
NJ = Nmero de nudos.
2 paso.- Cada miembro tiene 4 sub-matrices, correspondientes a sus nudos, estas sub-matrices se colocan en el
casillero correspondiente de la matriz de rigidez del prtico, Por ejemplo, el miembro I = 1, tiene las conectividades
J = 6 y K = 4, y est dividida en 4 sub-matrices las cuales se colocan en los casilleros segn los sub-ndices, como se
indica en la siguiente representacin de una sub-matriz: Entonces la sub-matriz [1 K66 ] se ubica en el casillero de la
fila 6 y la columna 6, la sub-matriz [1 K64 ] se ubica en la fila 6 y la columna 4, y as sucesivamente con las otras submatrices quedando tal como se muestra en la siguiente figura; el mismo procedimiento se realiza con las sub matrices de empotramiento, desplazamiento y fuerzas:
U
1
R
4
1
2
U
{ F4 }
3
4
[1K44]
[1K46]
{ D4 }
{ 1F4 }
{ F6 }
6
R
[1K64]
[1K66]
{ D6 }
F6 }
7
8
Fig. 18 Ensamble de la matriz de rigidez total de nudo, proceso de ensamblaje para el miembro I = 1.
Despus proseguimos ensamblando con el segundo miembro I = 2, J = 7 y K = 5, de la cual sus 4 sub matrices se aaden en la matriz de rigidez global y posteriormente proseguimos con el miembro I = 3, J = 4
y K = 1, donde la sub-matriz [K44 ] se coloca en casillero (4,4) y encontramos que tenemos dos sub matrices en el mismo casillero, esto significa que los miembros con el cdigo 1 y 3 tienen el mismo nudo
conectivo y por lo tanto se superponen sumndose todas las sub-matrices que se colocan en un mismo
casillero esto se produce solamente en la diagonal (ver Fig. 19). As sucesivamente se va ensamblando la
matriz de rigidez de nudo con los dems miembros hasta completar todos ellos.
PAG. 26
U
1
{F3}
R
4
8
F
[3K11]
[3K14]
{ D1 }
{ 3F1
[3K41]
[1K44]+
[3K44]
{ D4 }
{ 1F4 }+
F
{ 3F4 }
{F4}+
{F4}
=
5
{F6}
{F7}
[1K46]
+
[2K55]
[1K64]
[1K66]
[2K75]
[2K77]
{ 2F5
[2K57]
{ D6 }
{ 1F6 }
{ D7 }
{ 2F7
Fig. 19 Proceso de ensamblaje de las matrices con los miembros 1,2 y 3.
PAG. 27
{ F1 }
{ F2 }
[3K11] +
[8K11] +
[9K11]
[9K12]
[0]
[3K14]
[8K15]
[0]
[0]
[0]
[9K21]
[4K22] +
[9K22] +
[10K22]
[10K23]
[10K32]
[5K33] +
[10K33]
[0]
[0]
[1K44] +
[3K44] +
[6K44] +
[7K44]
[7K45]
[7K54]
[2K55] +
[4K55] +
[7K55] +
[8K55]
[0]
[0]
[1K66]
[0]
{ F3 }
{ F5 }
[8K51]
{ F6 }
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
{ 4F2 }+
F
{ 9F2 }+
F
{10F2 }
F
[5K38]
[1K64]
[1K46]
[6K47]
{ 5F3 }+
F
{10F3 }
[0]
{ D4 }
[2K57]
[0]
{ D5 }
[0]
{ 1F4 }+
F
{ 3F4 }+
F
{ 6F4 }+
F
{ 7F 4 }
{
{
{
{
2F5
}+
}+
F
7F5 }+
F
8F5 }
F
4F5
[0]
{ 1F6 }
{ D6 }
[0]
[0]
[6K74]
[2K75]
[0]
[2K77] +
[6K77]
[0]
{ D7 }
[0]
8
[0]
[0]
{ F8 }
[4K52]
[4K25]
{ D3 }
[0]
{ F7 }
[0]
[0]
{ D2 }
{ 3F1 }+
F
{ 8F1 }+
F
{ 9F1 }
{ D1 }
[3K41]
{ F4 }
[0]
[5K83]
[0]
[0]
[0]
[0]
[5K88]
{ 2F7 }+
F
{ 6F7 }
{ 5F8 }
{ D8 }
PAG. 28
{ F1 }
{ F2 }
[3K11] +
[8K11] +
[9K11]
[9K12]
[0]
[3K14]
[8K15]
[9K21]
[4K22] +
[9K22] +
[10K22]
[10K23]
[10K32]
[5K33] +
[10K33]
[0]
[0]
[1K44] +
[3K44] +
[6K44] +
[7K44]
[7K45]
[7K54]
[2K55] +
[4K55] +
[7K55] +
[8K55]
[0]
{ F3 }
{ D1 }
[3K41]
{ F4 }
[0]
[0]
[4K52]
{ 4F2 }+
F
{ 9F2 }+
F
{10F2 }
F
[0]
{ D3 }
[0]
{ F5 }
[4K25]
{ D2 }
[8K51]
{ 3F1 }+
F
{ 8F1 }+
F
{ 9F1 }
{ 5F3 }+
F
{10F3 }
{ 1F4 }+
F
{ 3F4 }+
F
{ 6F4 }+
F
{ 7 F4 }
{ D4 }
{
{
{
{
{ D5 }
2F5
}+
}+
F
7F5 }+
F
8F5 }
F
4F5
MATRIZ DE CARGAS, EN EQUILBRIO CON LA FUERZA ELASTICA Y CARGA EQUIVA LENTES DE EXTREMO FIJO DE LA ESTRUCTURA
10

0
0
0

0
0
0

FF 0
0

5
0
0

0
0
0

8744.79
2262.74
123.07
6400
0
0
0
0
0
82.04
0
123.07
2262.74
2262.74
0
Matriz de Cargas Nodales
2262.74 123.07
6400
82.04
8744.79 123.07
82.04
123.07
6400
123.07 2262.74 2262.74

0
123.07 123.07
2262.74 2262.74
123.07
492.27
123.07
123.07
12882.04
123.07
6400
82.04
82.04
123.07
6564.09
82.04
123.07
6400
123.07
123.07
123.07
738.4
123.07 123.07
123.07
6400
6410.26
82.04 123.07
3282.04 123.07
123.07
123.07
30.77
123.07 369.2
0.03
6400
82.04
123.07
30.77
123.07
8826.85 2262.74
6400
2262.74 15144.79 123.07
123.07
0.03
123.07
738.45
2262.74
82.04
123.07
6400
2262.74
6400
82.04
123.07 2262.74 15144.79
123.07
123.07
123.07
123.07
Matriz global de rigidez
8826.83 2262.74
0
123.07
0
123.07
0
123.07
0
0
0
123.07
123.07
0
123.07
738.4
0
D1
D2
D3
D
4
D
5
D6
D7
D8
D
9
D
10
D
11
D12
D13
D
14
D
15
Matriz de desplazamientos globales
ECUACIN MATRICIAL DE EQUILIBRIO DE LA ESTRUCTURA
3.75
1.875
0
7.5
0
0
3.75
1.875
0
3.75
1.875
0
3.75
1.875
Matriz de fuerzas de empotramiento
ANALISIS ESTRUCTURA L 2
INGENIERIA CIVIL DE LA UNSA
MATRICES COMPLETAMENTE ENSAMBLADAS PARA EL PORTICO

La matriz de rigidez se divide en funcin a los nudos libres y restringidos, sabemos que los nudos del 1 al 5
son libres y los nudos 6, 7 y 8 estn restringidos entonces podemos dividir dicha matriz en la forma siguiente:
Donde:
U = Grados de libertad.
R = Grados restringidos.
Resolviendo la anterior ecuacin matricial hallamos los valores correspondientes a la matriz de desplazamientos
en coordenadas globales:
{D} = [K]-1 ({F}- {FF})
T
DD
1
0.02146
Fidel Copa P.
0.0036 -0.00703
0.02131 -0.00465 -0.00223
0.02125 -0.00115
29
10
0.00522
0.00874
11
12
0.00262 -0.00583
13
14
15
0.01021 -0.00339 -0.00067
Procedimiento en MATHCAD V4.0

PROGRAMA MATRIZ DE RIGIDEZ DE MARCO PLANO 3 GDL/NUDO
DATOS DE M IEM BROS
E
I
TB
21016
9.2310 -5
1
0.0096
1
21026
9.2310 -5
2
0.0096
2
21
21036
9.2310 -5
3
0.0096
3
31
21046
9.2310 -5
4
0.0096
4
41
21056
9.2310 -5
5
0.0096
5
51
21066
9.2310 -5
6
0.0096
6
64
21076
9.2310 -5
7
0.0096
7
71
21086
9.2310 -5
8
0.0096
8
84
21096
9.2310 -5
9
0.0096
9
91
6
21010
-5
9.2310 10
0.0096
10
10
1
Donde
E=Mdulo de Elasticidad
I =Mto. de Inercia
NUMERO DE BARRAS
NB row s( A )
A=Area
TB=tipo de Miembro
NB 10
COORDENADAS DE LOS NUDOS

T
X ( 0 3 6 0 3 0 3 6 )
Y ( 6 6 6 3 3 0 0 0 )
NUMERO DE NUDOS
NN row s( X)
NN 8
NUDO INICIAL (J) Y NUDO FINAL (K) PARA CADA BARRA

T
J ( 6 7 4 5 8 4 4 1 1 2 )
K ( 4 5 1 2 3 7 5 5 2 3 )
RESTRICCONES EN (X, Y, PARA CADA NUDO

T
rx ( 0 0 0 0 0 1 1 1 )
ry ( 0 0 0 0 0 1 1 1 )
r ( 0 0 0 0 0 1 1 1 )
MATRIZ DE CARGAS NODALES
10
FF
0
0
NNC cols ( FF)
NNC 2
CARGAS EQUIVALENTES DE EXTREM O FIJO DE M IEMBROS

T
fe ( 0 0 0 0 0 0 )
1
fe ( 0 0 0 0 0 0 )
2
fe ( 0 0 0 0 0 0 )
3
fe ( 0 0 0 0 0 0 )
4
fe ( 0 0 0 0 0 0 )
5
Fidel Copa P.
fe ( 0 0 0 0 0 0 )
6
fe ( 0 3.75 1.875 0 3.75 1.875)

7
fe ( 0 0 0 0 0 0 )
8
fe ( 0 3.75 1.875 0 3.75 1.875)

9
fe
10
( 0 3.75 1.875 0 3.75 1.875)
30
Grafico del ejemplo prtico 2D
2
1
15
2.5Tn/m
13
5Tn
3.00 m
20
21
16
23
19
24
17
3.00 m
4
14
10
18
10Tn
11
12
10
2.5Tn/m
8
4
22
3.00 m
3.00 m
3.00 m
3.00 m
Fig. 16 Prtico del ejemplo

DATOS GENERALES:
Nmero de nudos = 8
Nmero de barras = 10
BLOQUE DE DATOS DE NUDOS:

Coordenadas
X
Y
NUDO
Irx
Restricciones
Iry
1
2
3
4
5
6
7
8
0
3
6
0
3
0
3
6
6
6
6
3
3
0
0
0
BLOQUE DE DATOS DE MIEMBROS:

I
1
2
3
4
5
6
7
8
9
10
J
6
7
4
5
8
4
4
1
1
2
Fidel Copa P.
K
4
5
1
2
3
7
5
5
2
3
0
0
0
0
0
1
1
1
0
0
0
0
0
1
1
1
I = 9.232E-5
Tipo De Elemento
INERCIA
I
I
I
I
I
I
I
I
I
I
31
0
0
0
0
0
1
1
1
A = 9.6E-3
AREA
A
A
A
A
A
A
A
A
A
A
E = 2E6
LONG
L
L
L
L
2L
2
L
2
L
L
E
E
E
E
E
E
E
E
E
E
E
L=3
Cx
0
0
0
0
0
1/
1/
2
2
1
1
1
Cy
1
1
1
1
1
-1/
-1/
0
0
0
2
2
NN
CALCULO DE LA MATRIZ N QUE ALMACENA

LOS GRADOS DE LIBERTAD DE CADA NUDO
N
j1
for i 1 NN
if rx
1 i
j j 1
M
1 i
0 otherwise
if ry
2 i
j j 1
M
2 i
0 otherwise
if ri
3 i
NGL max( N)
NGL
j j 1
M
3 i
0 otherwise
M
MATRICES DE TRANSFORMACION DE COORDENADAS PARA CADA M IEMBRO i
L
for i 1 NB
M
i
XK XJ 2 YK YJ 2
i
M
R
for i 1 NB
X X
K
J
C
i
YKi YJi
L
i
Ci Si 0 0
Si Ci 0 0
0 0 1 0
M
i
0 0 0 Ci
0 0 0 Si
0 0 0 0
0 0
0 0
0 0
S 0
i
C 0
i
0 1
Fidel Copa P.
32
MATRIZ DE RIGIDEZ EN COORD. LOCALES PARA CADA TIPO DE MIEMBRO i

i 1 NB
A
Ai
0
0
0
0
A
L
Ai
Li
i
i
0
0
0
0
L
Li
12 I
6 I
12 I
6 I
i
i
i
i
i
0
0
3I
3 I
3I
3
2
3
2
i
i
i
L
L
L
L
i
i
i
i
0 0
0
3
3
2
L
L
L
6 I
4 I
6 I
2 I
i
i
i
i
i
i
i
0
0
0
0
0 0
0
0
2
L
2
L
i
i
L
L
k2i Ei
i
i
A
A
k1 E
i
i A
i 0 0 i
A
0
0
i
i
L
Li
0
0
0
0
i
L
i
i
3 I
3I
3 I
i
i
i
12 I 6 I
12 I
6 I
0 0
i
i
i
i
0
3
3
2
0
0
L
L
L
i
i
i
3
2
3
2
L
L
L
L
i
i
i
i
0
0 0
0
0
0
6 I
2 I
6 I
4 I
i
i
i
i
0
0
A
A i
2
L
2
L
i
i
i
L
L
0
0
0
0
i
i
A i
Li
0
0
k4 E
i
i
Ai
L
i
0
0
A
0 0
0 0
0 0
A
0 0
0 0
0 0
i
i
3I
3I
i
i
0
0
3
2
Li Li
3I
3I
i
i
0
0
2
L
k3 E
i
i
i
Li
A
A i
i
0
0
L
Li
i
3 I
3 I
i
i
0
0
3
2
Li Li
0
0
0
0
0
0
i
0
3
Li
3 I
i
0
2
Li
0
0
3I
i
0
Li3
0
0
3 I
ASIGNACION DE LA M ATRIZ DE RIGIDEZ EN COORD. LOCALES POR MIEM BRO

KL
for i 1 NB
M k1 if TB
M k2 if TB
M k3 if TB
M k4 if TB
i
i
i
i
Fidel Copa P.
33
CALCULO DE LAS M ATRICES DE RIGIDECES GLOBALES PARA CADA MIEMBRO i

i 1 NB
T
i
KGE R KL R
i
EMSAMBLE DE LA MATRIZ DE RIGIDEZ EN COORDENADAS GLOBALES

KG
f0
for i 1 NGL
for j 1 NGL
M
i j
for i 1 NB
for j 1 3
for k 1 3
aa N
k Ji
bb N
M
j Ji
aa bb
aa bb
KGE
k j
if aa 0 bb 0
aa N
k Ki
bb N
M
j Ji
aa bb
aa bb
KGE
3 k j
if aa 0 bb 0
k 3 j
if aa 0 bb 0
aa N
k Ji
bb N
M
j Ki
aa bb
aa bb
KGE
aa N
k Ki
bb N
M
aa bb
j Ki
aa bb
KGE
i 3 k 3 j
Fidel Copa P.
34
if aa 0 bb 0
MATRIZ DE RIGIDEZ DE LA ESTRUCTURA EN COORDENADAS GLOBALES

KG
1
10
11
1 8744.8 -2262.7
123.1
-6400
-82
2-2262.7
8744.8
123.1
-82
123.1
-6400
3 123.1
123.1
492.3
-123.1
123.1
-123.1
123.1
4 -6400
12882
123.1
-6400
-82
123.1
12
13
14
15
123.1 -2262.7
2262.7
2262.7 -2262.7
-82
-123.1
6564.1
-82
123.1
-6400
123.1
123.1
123.1
738.4
-123.1
123.1
-123.1
123.1
-6400
6410.3
30.8
-82
-123.1
3282
-123.1
123.1
123.1
30.8
-123.1
369.2
10
-82
-123.1
8826.8 -2262.7
-6400
11
-6400
0 -2262.7 15144.8
123.1
-82
123.1
12 123.1
123.1
123.1
738.5
-123.1
123.1
13-2262.7
2262.7
-82
-123.1
-6400
142262.7 -2262.7
-6400
-82
15
123.1
123.1
123.1
8826.8 -2262.7
-123.1 -2262.7 15144.8
-123.1
123.1
-123.1
738.4
MATRICES DE CARGAS DE EQUIVALENTES DE EXTREMO FIJO POR MIEMBRO

i 1 NB
T
i
FEG R fe
i
ENSAM BLE DEL VECTOR TOTAL DE CARGAS DE EQUIVALENTES DE EXTREMO FIJO

FE
for i 1 NGL
M 0
i
for i 1 NB
for j 1 3
aa N
M
aa
j Ji
bb N
M
bb
aa
j Ki
FEG
bb
FEG
if aa 0
j 3
if bb 0
FET
1
1
0
Fidel Copa P.
2
3.75
3
1.875
4
0
5
7.5
6
0
7
0
8
9
3.75 -1.875
35
10
0
11
3.75
12
1.875
13
0
14
15
3.75 -1.875
ENSAM BLE DEL VECTOR TOTAL DE CARGAS NODALES

FN
for i 1 NGL
M 0
i
for i 1 NNC
for j 1 3
aa N
M
aa
j FF1 i
FF
if aa 0
j 1 i
M
T
FN
1
10
10 11 12 13 14 15
VECTOR DE DESPLAZAMIENTOS DE LA ESTRUCTURA EN COORD. GLOBALES

1
DD
DD KG
( FN FE)
10.0215
0.0036
-0.007
0.0213 -0.0046 -0.0022
10
0.0213 -0.0011
0.0052
0.0087
11
12
0.0026 -0.0058
13
DESPLAZAMIENTOS POR MIEM BRO DE LA ESTRUCTURA

D
for i 1 NB
for j 1 3
aa N
M
j i
j i
j Ji
DD
if aa 0
aa
0 otherwise
bb N
M
j 3 i
j 3 i
j Ki
DD
bb
if bb 0
0 otherwise
for i 1 NB
N submatrix( M 1 6 i i)
i
CALCULO DE FUERZAS INTERNAS EN COORD. LOCALES EN MIEM BROS
f KL R D fe
i
DESPLAZAMIENTOS DE MIEMBROS EN COORDENADAS GLOBALES
Fidel Copa P.
36
14
15
0.0102 -0.0034 -0.0007
0
0
0
D
1 0.0087
0.0026
0.0058
0
0
D
2 0.01021
0.00339
0.00067
0.00874
0.00262
0.00583
D
3 0.02146
0.0036
0.00704
0.01021
0.00339
0.00067
D
4 0.02131
0.00465
0.00223
0
0
D
5 0.02125
0.00115
0.00522
0.0087
0.0026
0.0058
D
6
0
0
0.00874
0.00262
0.00583
D
7 0.01021
0.00339
0.00067
0.02146
0.0036
0.00704
D
8 0.01021
0.00339
0.00067
0.02146
0.0036
0.00704
D
9 0.02131
0.00465
0.00223
0.02131
0.00465
0.00223
D
10 0.02125
0.00115
0.00522
MATRICES DE TRANSFORMACION DE COORD. DE MIEMBROS
0
1
0
R
1 0
0
1 0 0 0 0
0
1
0
R
5 0
0
1 0 0 0 0
0 0 0 0 0
0 1 0 0
0 0 0 1
0 0 1 0
0 0 0 0
0.707 0.707
0.707 0.707
0
0
R
6 0
0
0
0
0
0
0
0
0 0 0 0 0
0 1 0 0
0 0 0 1
0 0 1 0
0 0 0 0
1
0
0
R
7 0
0
0
0
0 0.707 0.707 0
0 0.707 0.707 0
0
0 0 0 0 0
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 0 1
0 0 0 1 0
MATRICES DE RIGIDEZ DE MIEMBROS EN COORDENADAS LOCALES

6400
0
0
KL
1 6400
0
0
3200
0
0
KL
5 3200
0
82.062 123.093
0
82.062 123.093
123.093 246.187
0
123.093 123.093
0
0
6400
0
0
82.062 123.093 0
82.062 123.093
123.093 123.093
0
123.093 246.187
0
6400
3200
10.256 30.767
10.256 30.767
30.767 123.067
30.767 61.533
3200
10.256 30.767
30.767 61.533
30.767 123.067
0
10.256 30.767
4525.483 0
0
0
0
0
KL
6 4525.483 0
0
0
0
0
6400
0
0
KL
7 6400
0
0 4525.483 0 0
0
0 0
0 0
0 4525.483 0 0
0
0
0
0 0
0 0
6400
37
82.044 123.067
123.067 246.133
123.067 123.067
6400
82.044 123.067
123.067 123.067
VECTORES DE FUERZAS INTERNAS EN MIEMBROS EN COORD. LOCALES
Fidel Copa P.
82.044 123.067
123.067 246.133
0
82.044 123.067
16.7448
0.0003
0.3585
f
1 16.7448
0.0003
0.3594
21.709
0.755
1.174
f
2 21.709
0.755
1.092
19.61
0
0
f
6 19.61
0
9.407
3.443
1.096
f
7 9.407
4.057
2.018
6.321
0.541
0.737
f
3 6.321
0.541
0.885
8.043
0.554
0.927
f
4 8.043
0.554
0.735
3.669
0.378
0.975
f
5 3.669
0.378
1.296
13.588
0
0
f
8 13.588
0
0.932
3.287
0.885
f
9 0.932
4.213
2.273
0.378
3.831
1.538
f
10 0.378
3.669
1.296
VECTORES DE DESPLAZAMIENTOS EN MIEMBROS EN COORD. LOCALES

[1 R] {1 D} = {1 d}
0
1
0
R
1 0
0
1 0 0 0 0
0 0 0 0
0 1 0 0
0 0 0 1
0 0 1 0
0 0 0 0
0
0
0
0
D
1 0.00874
0.00262
0.00583
d R D
1
0
0
d
1 0.00262
0.00874
0.00583
VECTORES DE FUERZAS INTERNAS EN MIEMBROS EN COORD. LOCALES

[1 k] {1 d} + {1 fF} = {1 f}
16.745
0
0.358
KL d fe
1 1
1 16.745
0
0.359
Fidel Copa P.
38
4.213
2.273
1.538
-
0.735
1.25
1.296
-
1.40
1.096
4.057
+ 2.019
3.287
3.831
+
0.359
1.296
0.885
0.885 -
3.669
0.737
0.927
+
0.554
0.541
1.092
1.27
3.443
1.175 +
0.358
0.975
0.756
DIAGRAMA DE MOMENTOS FLECTORES
6.321
-1
6.
DIAGRAMA DE FUERZAS CORTANTES
-0.379
-8.043
-0.932
74
9.407
-21.709
-3.669
16.74
-1
9.
61
DIAGRAMA DE FUERZAS AXIALES
Fidel Copa P.
39
0.379

Análisis Estructural Por El Método de Rigidez en 2D

Cargado por

Copyright:

Formatos disponibles

Análisis Estructural Por El Método de Rigidez en 2D

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Análisis Estructural Por El Método de Rigidez en 2D

Cargado por

Copyright:

Formatos disponibles

ANALISIS EST RUCTURAL 2

Anlisis Estructural por el mtodo de Rigidez en 2D

2.0 MATRIZ DE RIGIDEZ DE UN ELEMENTO ESTRUCTURAL

FUERZAS INTERNAS DE EXTREMO EN UNA ESTRUCTURA

ANALISIS EST RUCTURAL 2

Estructura ensamblada (completa) y codificada segn sus miembros estructurales.

Fig. 2 Estructura desensamblada del prtico de la fig. 1.

FUERZAS RESTAURADORAS DE ELEMENTOS ESTRUCTURALES

ANALISIS EST RUCTURAL 2

= sub-matriz de 3x3 de rigidez del miembro I-simo en coordenadas locales.

ANALISIS EST RUCTURAL 2

La ecuacin (2) en forma compacta, puede escribirse as:

Fig. 4 Codificacin de fuerzas y desplazamientos en coordenadas locales de un miembro estructural I-simo.

La matriz de rigidez de un miembro estructural I-simo en coordenadas locales expresada en sub-matrices

[Ib JJ] [Ib JK]

MATRICES DE RIGIDEZ EN COORDENADAS LOCALES

ANALISIS EST RUCTURAL 2

2.3.1 Matriz de Rigidez de un miembro estructural, con los nudos extremos J y K

2.3.2 Matriz de Rigidez para un miembro con articulacin en el nudo J y rgido en el

ANALISIS EST RUCTURAL 2

ANALISIS EST RUCTURAL 2

ANALISIS EST RUCTURAL 2

Cuyos parmetros son:

Mdulo de elasticidad del material.

Arw = rea Reducida por deformacin por corte.

ANALISIS EST RUCTURAL 2

2.3.6 Matriz de Rigidez con zonas rgidas considerando la deformacin por

ANALISIS EST RUCTURAL 2

CARGAS EQUIVALENTES DE EXTREMO FIJO

3.1 ESFUERZOS DE EXTREMO, CARGAS EQUIVALENTES DE EXTREMO Y

Cargas en un nudo sobre un miembro de una estructura.

ANALISIS EST RUCTURAL 2

Cargas equivalentes de extremo fijo en un sistema de coordenadas locales.

Fuerzas de extremo y desplazamientos de un miembro en un sistema de coordenadas locales.

ANALISIS EST RUCTURAL 2

TRANSFORMACION DE COORDENADAS de FUERZAS Y DESPLAZAMIENTOS

Fig. 14a Transformacin de fuerzas de extremo de un sistema de coordenadas de locales a globales.

cx = cos I = (xK - xJ)/ LI

y en forma compacta la fuerza en un nudo se expresa con la submatriz de 3x3:

{F} = [A]T {f}

ANALISIS EST RUCTURAL 2

4.2 Transformacin de coordenadas del vector de fuerzas de un miembro I-simo,

[IA]T {IfJ} + [0] {IfK}

Factorizando los vectores de fuerzas locales de los nudos J y K, se obtiene:

Expresando en forma compacta:

de ello, podemos encontrar la siguiente transformacin coordenadas globales a locales :

ANALISIS EST RUCTURAL 2

Transformacin de coordenadas del vector de desplazamientos

para un miembro I-simo, anlogamente al procedimiento descrito en la seccin 4.1.2, se obtiene:

4.4 Transformacin de coordenadas locales del vector de fuerzas de empotramiento perfecto

Reemplazando el vector de fuerzas restauradoras a un sistema de coordenadas globales, se tiene

transformando el vector de desplazamientos a un sistema de coordenadas globales, la ecuacin (27)

ANALISIS EST RUCTURAL 2

Transformacin de coordenadas del vector de fuerzas de extremo de un miembro I

{If} = [Ik] {Id} + [IfF]

transformando el sistema de coordenadas de este vector de fuerzas a un sistema de coordenadas globales,

efectuando el producto matricial, obtenemos:

{IF} = [IR]T [Ik] {Id} + [IR]T {IfF}

{IF} = ( [IR] T [Ik] [IR] ) {ID} + [IR]T {IfF}