Análisis de estructuras, método de la rigidez

Capítulo 2
CERCHAS

2.1 Ejercicio 1. Cercha sencilla con tres elementos

Para la cercha en acero mostrada en la figura 1.1-a. Determine el

desplazamiento horizontal y vertical en el punto C debido a la acción de la
carga que allí actúa, considere el módulo de elasticidad del acero (Es) igual a
200 000 MPa.

Figura 2.1- a

Resolución del ejercicio:

 Propiedades de la sección

Área de la sección: = 4 𝜋ɸ
Área= 1.0 cm ≈ 0.0001 m2
2

21
 Análisis de estructuras, método de la rigidez

Discretización de la estructura

Se enumera los elementos de la cercha y luego sus grados de libertad

empezando por aquellos que tienen restricción cinemática (que tendrán lugar a
las reacciones) cuyo desplazamiento será nulo, para dar facilidad a las
operaciones matriciales posteriores que permitirán calcular los desplazamientos
y reacciones de esta cercha.

Figura 2.1-b.

Longitud y ángulos de rotación de los elementos

Figura 2.1-c.

22
 Análisis de estructuras, método de la rigidez

Elemento No 1: (ver figura 2.1-c)

L= 1.0 m
Angulo de rotación (α1):
α1= 0° (no hay rotación permanece en su posición local)
α1= 0 rad

Elemento No 2: (ver figura 2.1-c)

L= 1.0 m
Angulo de rotación (α2):
α2= 90°
α2= 1.570 rad

Elemento No 3: (ver figura 2.1-c)

L= √𝟏𝟐 + 𝟏𝟐
L= 1.414 m
Angulo de rotación (α3):
.0
α3 = 𝑇𝑎𝑛− + 90° (Respecto al eje global X positivo)
.0
α3= 135°
α3= 2.356 rad

 Resumen de las propiedades geométricas de los elementos

ELEMENTO ÁREA (m2) LONGITUD (m) ÁNGULO

Elemento 1 0.0001 1.0 0°
Elemento 2 0.0001 1.0 90°
Elemento 3 0.0001 1.414 135°

23
 Análisis de estructuras, método de la rigidez

Matriz de rigidez local y global de los elementos

La matriz de rigidez local de un elemento cercha expresando sus grados de

libertad numéricamente como se expresó en capítulo 1, está dada por

1 2 3 4

0 0 1
-
0 0 0 0 2
[k]=
- 0 0 3

0 0 0 0 4

1 Figura
2 2.1-d.
3 4
Donde

A: es el área de la sección transversal del elemento

E: módulo de elasticidad del elemento
L: longitud del elemento

Remplazando los valores de área, longitud y módulo de elasticidad de los

elementos se obtiene la matriz de rigidez local de los elementos.

24
 Análisis de estructuras, método de la rigidez

 Elemento 1

E= 200000000 kPa

L= 1,00 m

A= 1,00 cm2

A= 0,00010 m2

Ѳ= 0,00 °

Ѳ= 0,00 rad

𝑬𝑨 𝟐𝑿𝟏𝟎𝟖 ∗ 𝟎. 𝟎𝟎𝟎𝟏
= = 𝟐𝟎 𝟎𝟎𝟎 𝒌𝑵/𝒎
𝑳 𝟏. 𝟎

Asociando el valor de la rigidez del paso anterior a la matriz local del elemento
tipo cercha se obtiene la matriz de rigidez local del elemento No 1 en kN/m.

1 2 3 4
20000 0 -20000 0 1

0 0 0 0 2
[ k1 ] =
-20000 0 20000 0 3

0 0 0 0 4

La matriz de rigidez del elemento se encuentra en coordenadas locales como se

aprecia en la figura 2.1-d. para pasar la matriz a coordenadas globales sería
necesario el uso de la matriz de rotación ó transformación de coordenadas en
caso que el elemento se encontrara inclinado respecto del eje global X positivo.

25
 Análisis de estructuras, método de la rigidez

Como el elemento 1 está en una posición horizontal, no hay necesidad de

aplicar la matriz de rotación ya que sus coordenadas locales coincide con las
Globales, aun así se realizara el ejercicio con fines académicos para vislumbrar
el proceso algebraico de operaciones con matrices.

La matriz de rotación del sistema para un elemento tipo cercha está dada por

cosƟ senƟ 0 0

-senƟ cosƟ 0 0
[T]=
0 0 cosƟ senƟ

0 0 -senƟ cosƟ

Reemplazando el valor del ángulo del elemento No 1 (0°) se obtiene la matriz

de rotación, resulta entonces

1,00 0,00 0,00 0,00

0,00 1,00 0,00 0,00

[T]=
0,00 0,00 1,00 0,00

0,00 0,00 0,00 1,00

Nota: como el ángulo es cero la matriz de rotación resulta ser la matriz

identidad.

La matriz de rigidez en coordenadas globales de un elemento está dada por:

[ K global ] = [ T’ ]*[ K local]*[ T ]

26
 Análisis de estructuras, método de la rigidez

Donde [T’] es la traspuesta de la matriz de rotación del sistema, para pasar de

[T] a [T’] se reordena de modo que las columnas de la matriz original se
convierten en las filas correspondientes de la matriz.

Se esta manera se obtiene que la matriz traspuesta de [ T ] será:

Para efectuar la multiplicación de matrices y para facilidad en la compresión del

método se plantean dos pasos:

Primero, la multiplicación de la matriz de rigidez local por la matriz de rotación

del sistema: [ K local]*[ T ]

1 2 3 4
20000 0 -20000 0 1 1,00 0,00 0,00 0,00

0 0 0 0 2 0,00 1,00 0,00 0,00

[ k 1 ]*[ T ] = X
-20000 0 20000 0 3 0,00 0,00 1,00 0,00

0 0 0 0 4 0,00 0,00 0,00 1,00

Para poder realizar la operación matricial de la multiplicación, el índice de

columnas de [k1] debe ser igual al número de filas de [ T ], recordando que
dentro del algebra matricial [A]*[B] ≠ [B]*[A].

27
 Análisis de estructuras, método de la rigidez

Multiplicación de [ k1 ]*[ T ]
Fila No 1

[k]*[ T ]1,1 = k1,1* T1,1 + k1,2* T2,1 + k1,3* T3,1 + k1,4* T4,1
[k]*[ T ]1,1 = (20000)*(1.0)+ (0.0)* (0.0) + (-20000)*(0.0) + (0.0)*(0.0)
[k]*[ T ]1,1 = 20 000 kN/m

[k]*[ T ]1,2 = k1,1* T1,2 + k1,2* T2,2 + k1,3* T3,2 + k1,4* T4,2
[k]*[ T ]1,2 = (20000)*(0.0) + (0.0)* (1.0) + (-20000)*(0.0) + (0.0)*(0.0)
[k]*[ T ]1,2 = 0.0

[k]*[ T ]1,3 = k1,1* T1,3 + k1,2* T2,3 + k1,3* T3,3 + k1,4* T4,3
[k]*[ T ]1,3 = (20000)*(0.0) + (0.0)* (0.0) + (-20000)*(1.0) + (0.0)*(0.0)
[k]*[ T ]1,3 = -20 000 kN/m

[k]*[ T ]1,4 = k1,1* T1,4 + k1,2* T2,4 + k1,3* T3,4 + k1,4* T4,4
[k]*[ T ]1,4 = (20000)*(0.0) + (0.0)* (0.0) + (-20000)*(0.0) + (0.0)*(1.0)
[k]*[ T ]1,4 = 0.0

Fila No 2

[k]*[ T ]2,1 = k2,1* T1,1 + k2,2* T2,1 + k2,3* T3,1 + k2,4* T4,1
[k]*[ T ]1,1 = (0.0)*(1.0) + (0.0)* (0.0) + (0.0)*(0.0) + (0.0)*(0.0)
[k]*[ T ]1,1 = 0.0

[k]*[ T ]2,2 = k2,1* T1,2 + k2,2* T2,2 + k2,3* T3,2 + k2,4* T4,2
[k]*[ T ]1,2 = (0.0)*(0.0) + (0.0)* (1.0) + (0.0)*(0.0) + (0.0)*(0.0)
[k]*[ T ]1,2 = 0.0

[k]*[ T ]2,3 = k2,1* T1,3 + k2,2* T2,3 + k2,3* T3,3 + k2,4* T4,3
[k]*[ T ]2,3 = (0.0)*(0.0) + (0.0)* (0.0) + (0.0)*(1.0) + (0.0)*(0.0)
[k]*[ T ]2,3 = 0.0

[k]*[ T ]2,4 = k2,1* T1,4 + k2,2* T2,4 + k2,3* T3,4 + k2,4* T4,4

28
 Análisis de estructuras, método de la rigidez

[k]*[ T ]2,4 = (0.0)*(0.0) + (0.0)* (0.0) + (0.0)*(0.0) + (0.0)*(1.0)

[k]*[ T ]2,4 = 0.0

Fila No 3

[k]*[ T ]3,1 = k3,1* T1,1 + k3,2* T2,1 + k3,3* T3,1 + k3,4* T4,1
[k]*[ T ]3,1 = (-20000)*(1.0) + (0.0)* (0.0) + (20000)*(0.0) + (0.0)*(0.0)
[k]*[ T ]3,1 = -20 000

[k]*[ T ]3,2 = k3,1* T1,2 + k3,2* T2,2 + k3,3* T3,2 + k3,4* T4,2
[k]*[ T ]3,2 = (-20000)*(0.0) + (0.0)* (1.0) + (20000)*(0.0) + (0.0)*(0.0)
[k]*[ T ]3,2 = 0.0

[k]*[ T ]3,3 = k3,1* T1,3 + k3,2* T2,3 + k3,3* T3,3 + k3,4* T4,3
[k]*[ T ]3,3 = (-20000)*(0.0) + (0.0)* (0.0) + (20000)*(1.0) + (0.0)*(0.0)
[k]*[ T ]3,3 = 20 000

[k]*[ T ]3,4 = k3,1* T1,4 + k3,2* T2,4 + k3,3* T3,4 + k3,4* T4,4
[k]*[ T ]3,4 = (-20000)*(0.0) + (0.0)* (0.0) + (20000)*(0.0) + (0.0)*(1.0)
[k]*[ T ]3,4 = 0.0

Fila No 4

[k]*[ T ]4,1 = k4,1* T1,1 + k4,2* T2,1 + k4,3* T3,1 + k4,4* T4,1
[k]*[ T ]4,1 = (0.0)*(1.0) + (0.0)* (0.0) + (0.0)*(0.0) + (0.0)*(0.0)
[k]*[ T ]4,1 = 0.0

[k]*[ T ]4,2 = k4,1* T1,2 + k4,2* T2,2 + k4,3* T3,2 + k4,4* T4,2
[k]*[ T ]4,2 = (0.0)*(0.0) + (0.0)* (1.0) + (0.0)*(0.0) + (0.0)*(0.0)
[k]*[ T ]4,2 = 0.0

[k]*[ T ]4,3 = k4,1* T1,3 + k4,2* T2,3 + k4,3* T3,3 + k4,4* T4,3
[k]*[ T ]4,3 = (0.0)*(0.0) + (0.0)* (0.0) + (0.0)*(1.0) + (0.0)*(0.0)
[k]*[ T ]4,3 = 0.0

29
 Análisis de estructuras, método de la rigidez

[k]*[ T ]4,4 = k4,1* T1,4 + k4,2* T2,4 + k4,3* T3,4 + k4,4* T4,4
[k]*[ T ]4,4 = (0.0)*(0.0) + (0.0)* (0.0) + (0.0)*(0.0) + (0.0)*(1.0)
[k]*[ T ]4,4 = 0.0

Reemplazando los valores obtenidos de los cálculos anteriores se ensambla la

matriz [k1]*[T]
1 2 3 4

1
20000 0 -20000 0
2
0 0 0 0
[K][ T ] =
3
-20000 0 20000 0
4
0 0 0 0

Finalmente para obtener la matriz de rigidez en coordenadas globales del

elemento se realiza la operación faltante, [K global]= [ T’]*[K local]*[T ] resulta

1 2 3 4 1 2 3 4
1 0 0 0 20000,0 0,0 -20000,0 0,0 1

0 1 0 0 0,0 0,0 0,0 0,0 2

[T'][k1][T] = x
0 0 1 0 -20000,0 0,0 20000,0 0,0 3

0 0 0 1 0,0 0,0 0,0 0,0 4

Multiplicación de [T’]* [ k1 ]*[ T ]

Fila No 1

[ T’ ]*[k]*[ T ]1,1 = [T’]1,1*[k]*[ T ] 1,1+ [T’]’1,2* [k]*[ T ] 2,1 +[T’]1,3* [k]*[ T ]

3,1+[T’]1,4* [k]*[ T ] 4,1
[T’]*[k][ T ]1,1 = (1.0)*(20000) + (0.0)* (0.0) + (0.0)*(-20000) + (0.0)*(0.0)

[T’]*[k][ T ]1,1 = 20 000 kN/m

30
 Análisis de estructuras, método de la rigidez

[ T’ ]*[ [k]*[ T ]1,2 = [T’]1,1* [k]*[T] 1,2+[T’]1,2* [k]*[T] 2,2+ [T’]1,3* [k]*[T]
3,2+[T’]1,4* [k]*[ T ] 4,2

[T’]*[k][ T ]1,2 = (1.0)*(0.0) + (0.0)* (0.0) + (0.0)*(0.0) + (0.0)*(0.0)

[T’]*[k][ T ]1,2 = 0,0

[ T’ ]*[ [k]*[ T ]1,3 = [T’]1,1* [k]*[T] 1,3+[T’]1,2* [k]*[T] 2,3+[T’]1,3* [k]*[T]

3,3+[T’]1,4* [k]*[ T ] 4,3

[T’]*[k][ T ]1,3 = (1.0)*(-20000) + (0.0)* (0.0) + (0.0)*(20000) + (0.0)*(0.0)

[T’]*[k][ T ]1,3 = -20 000 kN/m

[ T’ ]*[ [k]*[ T ]1,4 = [T’]1,1* [k]*[T] 1,4+[T’]1,2* [k]*[T] 2,4 + [T’]1,3* [k]*[T]
3,4 +[T’]1,4* [k]*[ T ] 4,4

[T’]*[k][ T ]1,4 = (1.0)*(0.0) + (0.0)* (0.0) + (0.0)*(0.0) + (0.0)*(1.0)

[T’]*[k][ T ]1,4 = 0,0

Fila No 2

[ T’ ]*[ [k]*[ T ]2,1 = [T’]2,1* [k]*[ T ] 1,1+[T’]2,2* [k]*[ T ] 2,1 +[T’]2,3* [k]*[ T
] 3,1 +[T’]2,4* [k]*[ T ] 4,1

[T’]*[k][ T ]1,1 = (0.0)*(20000) + (1.0)* (0.0) + (0.0)*(-20000) +

(0.0)*(0.0)

[T’]*[k][ T ]2,1 = 0.0

[ T’ ]*[ [k]*[ T ]2,2 =[T’]2,1* [k]*[ T ] 1,2+[T’]2,2* [k]*[ T ] 2,2 +[T’]2,3* [k]*[ T
] 3,2+ [T’]2,4* [k]*[ T ] 4,2

[T’]*[k][ T ]1,2 = (0.0)*(0.0) + (1.0)* (0.0) + (0.0)*(0.0) + (0.0)*(0.0)

31
 Análisis de estructuras, método de la rigidez

[T’]*[k][ T ]2,2 = 0.0

[ T’ ]*[ [k]*[ T ]2,3 = [T’]2,1* [k]*[ T ] 1,3+[T’]2,2* [k]*[ T ] 2,3 +[T’]2,3* [k]*[ T
] 3,3+[T’]2,4* [k]*[ T ] 4,3

[T’]*[k][ T ]2,3 = (0.0)*(-20000) + (1.0)* (0.0) + (0.0)*(20000) +

(0.0)*(0.0)

[T’]*[k][ T ]2,3 = 0.0

[ T’ ]*[ [k]*[ T ]2,4 =T’]2,1* [k]*[ T ] 1,4+[T’]2,2* [k]*[ T ] 2,4 +[T’]2,3* [k]*[ T ]
3,4+[T’]2,4* [k]*[ T ] 4,4

[T’]*[k][ T ]2,4 = (0.0)*(0.0) + (1.0)* (0.0) + (0.0)*(0.0) + (0.0)*(0.0)

[T’]*[k][ T ]2,4 = 0.0

Fila No 3

[ T’ ]*[ [k]*[ T ]3,1 =[T’]3,1* [k]*[ T ] 1,1+[T’]3,2* [k]*[ T ] 2,1+[T’]3,3* [k]*[ T ]

3,1+[T’]3,4* [k]*[ T ] 4,1

[T’]*[k][ T ]3,1 = (0.0)*(20000) + (0.0)* (0.0) + (1.0)*(-20000) + (0.0)*(0.0)

[T’]*[k][ T ]3,1 = -20 000 kN/m

[ T’ ]*[ [k]*[ T ]3,2 =[T’]3,1* [k]*[ T ] 1,2+[T’]3,2* [k]*[ T ] 2,2 +[T’]3,3* [k]*[ T
] 3,2 +[T’]3,4* [k]*[ T ] 4,2

[T’]*[k][ T ]3,2 = (0.0)*(0.0) + (0.0)* (0.0) + (1.0)*(0.0) + (0.0)*(0.0)

[T’]*[k][ T ]3,2 = 0.0

[ T’ ]*[ [k]*[ T ]3,3 =[T’]3,1* [k]*[ T ] 1,3+[T’]3,2* [k]*[ T ] 2,3 +[T’]3,3* [k]*[ T
] 3,3 +[T’]3,4* [k]*[ T ] 4,3
[T’]*[k][ T ]3,3 = (0.0)*(-20000) + (0.0)* (0.0) + (1.0)*(20000) + (0.0)*(0.0)

32
 Análisis de estructuras, método de la rigidez

[T’]*[k][ T ]3,3 = 20 000

[ T’ ]*[ [k]*[ T ]3,4 =[T’]3,1* [k]*[ T ] 1,4+[T’]3,2* [k]*[ T ] 2,4 +[T’]3,3* [k]*[ T
] 3,4 +[T’]3,4* [k]*[ T ] 4,4

[T’]*[k][ T ]3,4 = (0.0)*(0.0) + (0.0)* (0.0) + (1.0)*(0.0) + (0.0)*(0.0)

[T’]*[k][ T ]3,4 = 0.0

Fila No 4

[ T’ ]*[ [k]*[ T ]4,1 =[T’]4,1* [k]*[ T ] 1,1+[T’]4,2* [k]*[ T ] 2,1 +[T’]4,3* [k]*[ T
] 3,1 +[T’]4,4* [k]*[ T ] 4,1

[k]*[ T ]4,1 = (0.0)*(20000) + (0.0)* (0.0) + (0.0)*(-20000) +

(1.0)*(0.0)

[k]*[ T ]4,1 = 0.0

[ T’ ]*[ [k]*[ T ]4,2 =[T’]4,1* [k]*[ T ] 1,2+[T’]4,2* [k]*[ T ] 2,2 +[T’]4,3* [k]*[ T
] 3,2 + [T’]4,4* [k]*[ T ] 4,2

[T’]*[k][ T ]4,2 = (0.0)*(0.0) + (0.0)* (0.0) + (0.0)*(0.0) + (1.0)*(0.0)

[T’]*[k][ T ]4,2 = 0.0

[ T’ ]*[ [k]*[ T ]4,3 =[T’]4,1* [k]*[ T ] 1,3+[T’]4,2* [k]*[ T ] 2,3 +[T’]4,3* [k]*[ T
] 3,3 +[T’]4,4* [k]*[ T ] 4,3

[T’]*[k][ T ]4,3 = (0.0)*(-20000) + (0.0)* (0.0) + (0.0)*(20000) +

(1.0)*(0.0)

[T’]*[k][ T ]4,3 = 0.0

33
 Análisis de estructuras, método de la rigidez

[ T’ ]*[ [k]*[ T ]4,4 =[T’]4,1* [k]*[ T ] 1,4+[T’]4,2* [k]*[ T ] 2,4 +[T’]4,3* [k]*[ T
] 3,4 +[T’]4,4* [k]*[ T ] 4,4
[T’]*[k][ T ]4,4 = (0.0)*(0.0) + (0.0)* (0.0) + (0.0)*(0.0) + (1.0)*(0.0)

[T’]*[k][ T ]4,4 = 0.0

Reemplazando los valores obtenidos de los cálculos anteriores se ensambla la

matriz de rigidez global del elemento 1 que será igual a la local porque el
elemento permanece en la misma posición, al final de asocian los grados de
libertad globales establecidos en la discretizacion de la figura 2.1-b.

Globales

locales
Locales 1 2 3 4

Globales 1 2 4 3

20000,00 0,00 -20000,00 0,00 1 1

0,00 0,00 0,00 0,00 2 2

[K 1 ] =
-20000,00 0,00 20000,00 0,00 4 3

0,00 0,00 0,00 0,00 3 4

34
 Análisis de estructuras, método de la rigidez

La matriz de rigidez global del elemento 1 resulta entonces:

1 2 4 3

20000,00 0,00 -20000,00 0,00 1

0,00 0,00 0,00 0,00 2

[K 1 ] =
-20000,00 0,00 20000,00 0,00 4

0,00 0,00 0,00 0,00 3

De esta manera se realizan las operaciones matriciales para obtener la matriz

de rigidez global para cada elemento, como se puede apreciar son operaciones
tediosas de manera analítica por tal motivo es ventajoso el uso de programas
donde se puedan realizar operaciones matriciales de manera eficiente como lo
es Excel, Matlab, scilab entre otros.

 Elemento 2

E= 200000000 kPa

L= 1,00 m

A= 1,00 cm2

A= 0,00010 m2

Ѳ= 90,00 °

Ѳ= 1,57 rad

𝑬𝑨 𝟐𝑿𝟏𝟎𝟖 ∗ 𝟎. 𝟎𝟎𝟎𝟏
= = 𝟐𝟎 𝟎𝟎𝟎 𝒌𝑵/𝒎
𝑳 𝟏. 𝟎

35
 Análisis de estructuras, método de la rigidez

Matriz de rigidez local en kN/m

1 2 3 4
20000 0 -20000 0 1

0 0 0 0 2
[ k2 ] =
-20000 0 20000 0 3

0 0 0 0 4

Matriz de rotación para 90°

0,00 1,00 0,00 0,00

-1,00 0,00 0,00 0,00

[T]=
0,00 0,00 0,00 1,00

0,00 0,00 -1,00 0,00

Traspuesta de la matriz de rotación [T]

0,00 -1,00 0,00 0,00

1,00 0,00 0,00 0,00

[ T ' ]=
0,00 0,00 0,00 -1,00

0,00 0,00 1,00 0,00

36
 Análisis de estructuras, método de la rigidez

Realizando la operación matricial [ K global] = [ T’ ]*[ k local]*[ T ], se obtiene

la matriz de rigidez global del elemento, en la matriz resultante de la operación
se asocian los grados de libertad globales como se realizó para el elemento No
1.

Matriz de rigidez global del elemento 2 en kN/m

1 2 5 6

0,00 0,00 0,00 0,00 1

0,00 20000,00 0,00 -20000,00 2

[K 2 ] =
0,00 0,00 0,00 0,00 5

0,00 -20000,00 0,00 20000,00 6

 Elemento 3

E= 200000000 kPa

L= 1,41 m

A= 1,00 cm2

A= 0,00010 m2

Ѳ= 135,00 °

Ѳ= 2,36 rad

𝑬𝑨 𝟐𝑿𝟏𝟎𝟖 ∗ 𝟎. 𝟎𝟎𝟎𝟏
= = 𝟏𝟒 𝟏𝟒𝟒, 𝟐𝟕𝟏 𝒌𝑵/𝒎
𝑳 𝟏. 𝟒𝟏𝟒

37
 Análisis de estructuras, método de la rigidez

Matriz de rigidez local en kN/m

1 2 3 4
14144,27 0 -14144,27 0 1

0 0 0 0 2
[ k3 ] =
-14144,27 0 14144,27 0 3

0 0 0 0 4

Matriz de rotación para 135°

-0,71 0,71 0,00 0,00

-0,71 -0,71 0,00 0,00

[T]=
0,00 0,00 -0,71 0,71

0,00 0,00 -0,71 -0,71

Traspuesta de [T]

-0,71 -0,71 0,00 0,00

0,71 -0,71 0,00 0,00

[ T^]=
0,00 0,00 -0,71 -0,71

0,00 0,00 0,71 -0,71

38
 Análisis de estructuras, método de la rigidez

La matriz de rigidez local está planteada para una posición horizontal del
elemento como se expresó en el capítulo 1, en el caso de los elementos 2 y 3
de la cercha se encuentran girados respecto al eje global x positivo, el cual será
siempre el eje de referencia para medir los ángulos que se reemplazaran en la
matriz de rotación.

Ejemplo: El elemento 3 de la cercha se encuentra inclinado respecto al eje

global x positivo un ángulo de 135°, lo que hace en esencia la matriz de
rotación realizando la operación [K global] = [T´]*[k local]*[T] es redistribuir la
rigidez que aporta el elemento de su posición local a global ( ver figura 2.1-d).

Figura 2.1-d. Posición local del elemento 3 y relación son su posición

global en la estructura

Al final se asocian los gdl de correspondientes entre locales a globales (ver

figura 2.1-e) ya que el elemento ha sido girado.
gdl Globales

gdl Locales

gdl Locales 1 2 3 4
4 3 5 6
gdl Globales

7072,14 -7072,14 -7072,14 7072,14 4 1

-7072,14 7072,14 7072,14 -7072,14 3 2

[K 3 ] =
-7072,14 7072,14 7072,14 -7072,14 5 3

7072,14 -7072,14 -7072,14 7072,14 6 4

Figura 2.1-e. Asociación de los gdl del elemento 3 de la cercha

39
 Análisis de estructuras, método de la rigidez

Matriz de rigidez global del elemento 3 en kN/m

4 3 5 6

7072,14 -7072,14 -7072,14 7072,14 4

-7072,14 7072,14 7072,14 -7072,14 3

[K 3 ] =
-7072,14 7072,14 7072,14 -7072,14 5

7072,14 -7072,14 -7072,14 7072,14 6

Matriz de rigidez de la cercha

Para obtener la matriz de rigidez de toda la estructura, se tendrá en cuenta que

la rigidez concentrada en un nodo es la suma de las contribuciones de la rigidez
de todos los elementos estructurales conectados a tal nodo, por lo tanto se
suma la rigidez que aporta cada elemento de su matriz de rigidez global, al final
esta será cuadrada y simétrica del tamaño de los grados de libertad
establecidos en la discretización de la estructura, es decir matriz de K6x6.

Ejemplo:

K1,2 = (K1,2) e + K1,2 e2 + K1,2 e3

K1,2 = (0,0) + (0,0) + (0,0)
K1,2 = 0,0 kN/m

K 5, 3 = K 5,3 e1 + K 5,3 e2 + K 5,3 e3

K 5, 3 = 0,00 + 0,00 + 7072,14
K 5,3 = 7072,14 kN/m

K 6, 6 = K 6, 6 e1 + K 6, 6 e2 + K 6, 6 e3
K 6, 6= 0,00 + 20 000 + 7072,14
K 6, 6= 27 072,14 kN/m

40
 Análisis de estructuras, método de la rigidez

De esta manera se suman todas las rigideces que aportan cada elemento y se
ensambla la matriz de rigidez de toda la estructura.

Matriz de rigidez global de la cercha en kN/m

1 2 3 4 5 6
20000,00 0,00 0,00 -20000,00 0,00 0,00 1

0,00 20000,00 0,00 0,00 0,00 -20000,00 2

0,00 0,00 7072,14 -7072,14 7072,14 -7072,14 3
[KC]=
-20000,00 0,00 -7072,14 27072,14 -7072,14 7072,14 4
0,00 0,00 7072,14 -7072,14 7072,14 -7072,14 5

0,00 -20000,00 -7072,14 7072,14 -7072,14 27072,14 6

Los grados de libertad del 1 al 3 están asociados a las reacciones de la cercha y

su desplazamiento será 0, para esto enumeró primeramente los gdl donde
tendrán lugar las reacciones y queden agrupadas dentro de la matriz de rigidez.

Vector de fuerzas actuantes en la cercha (F) en kN

Es la representación de las fuerzas que operan en la estructura asociando el

grado de libertad donde actúan, para el caso de las reacciones se representan
como incógnitas (Ax, Ay y By) y el índice hace referencia al grado de libertad
donde se presentaran.

gdl Fuerzas
1 Ax
Fuerzas
2 Ay Desconocidas
3 By (Reacciones)

4 0,0
5 25,0 Fuerzas
Conocidas
6 0,0

41
 Análisis de estructuras, método de la rigidez

Vector de desplazamientos

Se sabe que La rigidez (K) está dada por:

𝑲=𝐔
𝐅
Ec. 2.1-a Donde F es la carga y U el desplazamiento elástico
que produce dicha carga.

Despejando F, resulta

F = K*U Ec. 2.1-b

La matriz de rigidez global de la cercha está estructurada como se muestra en

la figura 2.1-f, conforme a la distribución de los grados de libertad establecidos
en la discretización.

Representando la ecuación 2.1-b con los esquemas matriciales se obtiene la

representación general de la ecuación (ver figura 2.1-g).

Fuerzas Rigidez Desplazamientos

F
Ktt Kt0 0
desconocidas

F conocidas K0t K00 U

Figura 2.1-f. Representación matricial de la ecuación de la rigidez

Fuerzas 1 2 3 4 5 6 [U]
1 F1 20000,00 0,00 0,00 -20000,00 0,00 0,00 1 0
2

3
F2
F3
Fd 0,00
0,00
Ktt
20000,00
0,00
0,00
7072,14
0,00
-7072,14
Kt00,00
7072,14
-20000,00
-7072,14
2
3
Uc 0
0
= x
0,0 -20000,00 0,00 -7072,14 27072,14 -7072,14 7072,14 4 U4

K00
4

15,0 Fc
5

6 0,0
0,00
0,00
K0t 0,00
-20000,00
7072,14
-7072,14
-7072,14
7072,14
7072,14
-7072,14
-7072,14
27072,14
5
6
Ud U5
U6

Figura 2.1-g. Representación general de la matriz de rigidez global de

la cercha.

42
 Análisis de estructuras, método de la rigidez

Resolviendo la matriz, se obtiene

0
Fd = [Ktt] [Uc] + [Kto] [Ud]
Fd= [Kto][Ud] Ec. 2.1-c
0
FC = [K0t] [Uc] + [K00][Ud]
FC= [K00][Ud] Ec. 2.1-d

Despejando los desplazamientos desconocidos (Ud) de la ecuación 2.1-d, resulta

[Ud] = [K00]-1[FC] Ec. 2.1-e (Desplazamientos desconocidos de la estructura)

Y las fuerzas desconocidas (Reacciones) se calculan aplicando la ecuación 2.1-c

[Fd]= [Kt0] [Ud] (Reacciones de la estructura)

Se sustrae la sub matriz de rigidez donde están asociadas las fuerzas conocidas
(K00), para calcular los desplazamientos que estas producen en la cercha
aplicando la ecuación 2.1-e.

4 5 6

27072,14 -7072,14 7072,14 4

[ K 00 ] = -7072,14 7072,14 -7072,14 5

7072,14 -7072,14 27072,14 6

Obteniendo la inversa de la matriz K00:

4 5 6

0,0000500 0,0000500 0,0000000 4

[ K 00 ]-1 = 0,0000500 0,0002414 0,0000500 5

0,0000000 0,0000500 0,0000500 6

43
 Análisis de estructuras, método de la rigidez

Los desplazamientos generados por las fuerzas externas aplicadas sobre la

cercha serán: [U]= [K00]-1 [P]

4 5 6 Fuerzas
0,0000500 0,0000500 0,0000000 4 0,0 4

[U]= 0,0000500 0,0002414 0,0000500 5 X 25,0 5

0,0000000 0,0000500 0,0000500 6 0,0 6

Resolviendo matricialmente se obtiene:

U4= 0,00125 m El desplazamiento horizontal y vertical en el Nodo

U5= 0,006035 m c será:
U6= 0,00125 m U5= 0,006035 m ≈ 6,035 mm H →
U6= 0,00125 m ≈ 1,25 mm V ↑

Figura 1.1-h. Deformada de la cercha debido a la aplicación de la

carga horizontal de 25 kN en el nodo C.

44
 Análisis de estructuras, método de la rigidez

Fuerza interna del elemento AC

Se sustraen los desplazamientos globales del elemento AC (elemento 1)

teniendo en cuenta el número correspondiente a cada grado de libertad.

U1= 0,0
U2= 0,0
U5= 0,006035 m
U6= 0,00125 m

Es necesario conocer los desplazamientos locales del elemento para determinar

su fuerza interna, así como establecer si el elemento está sometido a esfuerzos
de tracción o compresión, para lo anterior se multiplica matricialmente la matriz
de rotación del elemento por los desplazamientos globales calculados, resulta
entonces

[U Locales]= [T]*[U Globales]

Donde la matriz de rotación “T” es

cosƟ senƟ 0 0
-senƟ cosƟ 0 0
[T] =
0 0 cosƟ senƟ
0 0 -senƟ cosƟ

Para el ángulo de 90° que es la inclinación del elemento 2 respecto al eje global
x positivo.

45
 Análisis de estructuras, método de la rigidez

Se establece la operación matricial

Matriz de Rotación del elemento 2 UG

0,00 1,00 0,00 0,00 0,000000 1

-1,00 0,00 0,00 0,00 0,000000 2

[ u2 ] = X
0,00 0,00 0,00 1,00 0,006035 5

0,00 0,00 -1,00 0,00 0,001250 6

u1= 0,0
u2= 0,0 Estos son los desplazamientos locales
u3= 0,00125 m del elemento 2.
u4= - 0,006035 m

Para obtener la fuerza axial interna del elemento se parte de la hipótesis

principal del método, donde la rigidez es igual a una carga F sobre el
desplazamiento elástico que esta produce.

F
𝐾=
U
F = [K local]* [U local] (elemento 2).

Se obtiene

Matriz de Rigidez local del elemento 2 UL

20000,00 0,00 -20000,00 0,00 0,000000 1

0,00 0,00 0,00 0,00 0,000000 2

[ f2 ] = X
-20000,00 0,00 20000,00 0,00 0,001250 3

0,00 0,00 0,00 0,00 -0,006035 4

46
 Análisis de estructuras, método de la rigidez

Resolviendo matricialmente se obtiene la fuerza interna del elemento:

f1 = -25,0 kN
f2 = 0,0 kN
f3 = 25,0 kN
f4 = 0,0 kN

Teniendo en cuenta que los valores obtenidos anteriormente corresponden a la

fuerza interna del elemento en sus coordenadas locales se determina el tipo de
esfuerzo al que está sometido el elemento, en este caso son de tensión, ya que
f1 es negativo es decir actúa en dirección contraria a la supuesta inicialmente,
mientras que f3 es positiva como se observa en la figura 2.1-i, como se
esperaba las fuerzas f2 y f4 son cero porque es la funcionalidad de los
elementos de una cercha ó armadura.

Figura 2.1-i

Fuerza axial del elemento será 25,0 kN (Tensión)

47