Análisis Estructural Por El Método de Rigidez en 2D
Análisis Estructural Por El Método de Rigidez en 2D
Análisis Estructural Por El Método de Rigidez en 2D
PAG. 1
2.1
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
2.2
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 ]
PAG. 4
(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.
[Ik] =
2.3
(4)
[Ib KK]
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
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.
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
1 =
4 =
(2 )
(1 + )
12
2
2 =
12
(1 + )3
5 =
(4 + )
(1 + )
3 =
6
(1 + )2
Donde:
E
G
v
Ff
=
=
=
=
PAG. 9
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
Fig. 11
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
Donde:
{IfF} = Vector de cargas equivalentes de extremo fijo en coordenadas locales del miembro I-simo.
J
X
Fig. 13
4.
PAG. 12
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
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)
(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)
(17)
PAG. 13
(18)
{IFJ}
[0]
{IfJ}
=
{IFK}
(19)
[IA]T
[0]
{IfK}
(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)
(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
(24)
Despejando tenemos:
{Id} = [IR] {ID}
(25)
(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)
(27)
(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)
en coordenadas locales a un sistema de coordenadas globales, la (ec-11) se encuentra con la siguiente expresin:
(11)
(30)
(31)
Reemplazando por: {Id} = [IR] [ID] en la ecuacin anterior, haciendo un arreglo, se tiene:
(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:
[IK] = [IR]T [Ik] [IR]
(29)
(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:
[IK] = [IR]T [Ik] [IR]
[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:
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
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.
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)
(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
(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 }
=
=
=
=
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}
=
=
=
=
=
=
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)
{DU }
(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
DATOS GENERALES:
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
Ir
0
0
0
0
0
1
1
1
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
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
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
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
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
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
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
Matriz de transformacin
de coordenadas para
miembros verticales
(1, 2, 3, 4, 5)
Matriz de transformacin
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]
MIEMBROS 1, 2, 3 y 4:
[RT]
[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
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:
PAG. 23
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).
.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
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
3.75
1.875
F
F
0
3.75
1.875
:
:
:
:
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
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
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
2262.74 123.07
6400
82.04
8744.79 123.07
82.04
123.07
6400
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
123.07
0.03
123.07
738.45
2262.74
82.04
123.07
6400
2262.74
6400
82.04
123.07
123.07
123.07
123.07
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
3.75
1.875
0
7.5
0
0
3.75
1.875
0
3.75
1.875
0
3.75
1.875
ANALISIS ESTRUCTURA L 2
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.02125 -0.00115
29
10
0.00522
0.00874
11
12
0.00262 -0.00583
13
14
15
ANALISIS ESTRUCTURA L 2
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
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
J ( 6 7 4 5 8 4 4 1 1 2 )
K ( 4 5 1 2 3 7 5 5 2 3 )
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 )
10
FF
0
0
NNC 2
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 0 0 0 0 0 )
8
fe
10
30
ANALISIS ESTRUCTURA L 2
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
Nmero de nudos = 8
Nmero de barras = 10
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
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
con extremos empotrados
con extremos empotrados
con extremos empotrados
con extremos empotrados
con extremos empotrados
con extremos articulados
con extremos empotrados
con extremos articulados
con extremos empotrados
con extremos empotrados
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
ANALISIS ESTRUCTURA L 2
NN
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
ANALISIS ESTRUCTURA L 2
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
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
ANALISIS ESTRUCTURA L 2
T
i
KGE R KL R
i
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
ANALISIS ESTRUCTURA L 2
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
123.1
-123.1
738.4
T
i
FEG R fe
i
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
ANALISIS ESTRUCTURA L 2
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
DD
DD KG
( FN FE)
10.0215
0.0036
-0.007
10
0.0213 -0.0011
0.0052
0.0087
11
12
0.0026 -0.0058
13
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
f KL R D fe
i
Fidel Copa P.
36
14
15
ANALISIS ESTRUCTURA L 2
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
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
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
Fidel Copa P.
82.044 123.067
123.067 246.133
0
82.044 123.067
ANALISIS ESTRUCTURA L 2
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
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
16.745
0
0.358
KL d fe
1 1
1 16.745
0
0.359
Fidel Copa P.
38
ANALISIS ESTRUCTURA L 2
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
6.321
-1
6.
-0.379
-8.043
-0.932
74
9.407
-21.709
-3.669
16.74
-1
9.
61
Fidel Copa P.
39
0.379