ANALISIS MATRICIAL DE ESTRUCTURAS

CAPITULO VI

CLASE N° 11:

11.- APOYOS CON DESPLAZAMIENTO INICIAL:

11.1. APOYOS CON DESPLAZAMIENTO INICIAL:. Introducción.

11.2 Problema 11.2.1

11.1.- Αpoyos elásticos.

Un asentamiento uniforme de toda la estructura no produce deformaciones ni daño en la

estructura; sin embargo, el hecho de que existan diferencias entre estos, provoca
irremediablemente esfuerzos adicionales a la estructura y, por ende, daño estructural si no se
realiza las medidas adecuadas correspondientes.

Quien desarrolla el análisis y cálculo de la estructura deberá hacer un análisis posterior

tomando en cuenta el desplazamiento (asentamiento) como un dato ya conocido o medido.

Este proceso se realiza haciendo un traslado de submatrices en la ecuación matricial

fundamental del análisis matricial, considerándose que inicialmente fue una incógnita
(desplazamiento) como un valor ya conocido en la EQUACION FUNDAMENTAL DE
METODO MATRICIAL DE RIGIDEZ.

Es decir, se parte de la ecuación ya conocida: la {F} = {K} * {U} y se descompone la MATRIZ

DE RIGIDEZ. Por consiguiente, las otras matrices en la partición de submatrices es como
sigue:

𝑭𝒏 𝑲𝒏𝒏 ‫׀‬ 𝑲𝒏𝒂 𝑼𝒏

{―} = { ― ― ― }∗{―}
𝑭𝒂 𝑲𝒂𝒏 ‫׀‬ 𝑲𝒂𝒂 𝑼𝒂

Docente : Mg. Ing. Evangelos Manouris clase 11 - 1 | 13

Donde:
• {Fn} = : Vector de cargas aplicadas (conocido).
• {Fa} = : Reacciones en los apoyos (desconocidas).
• {Un} = : Desplazamientos de los nodos libres (desconocidos).
• {Ua} = : Desplazamientos de los apoyos (conocidos y casi siempre cero; en nuestro caso,
con desplazamientos iniciales, existirán valores iniciales diferentes de cero).
Al expandir la expresión anterior se obtiene lo siguiente:
{Fn} = {Knn} * {Un} + {Kna} * {Ua} Ec. 11.1
{Fa} = {Kan} * {Un} + {Kaa} * {Ua} Ec. 11.2
Con las dos ecuaciones determinadas, combinadas y reemplazadas en forma adecuada, se
Obtienen los resultados deseados. Esto es a partir de la Matriz de Rigidez en coordenadas
Global de toda la estructura.

12.2.- Problema 11.2.1.

Resolver la estructura del reticulado plano mostrado a continuación. Dicha estructura
presenta en uno de sus apoyos un desplazamiento de 0.4 cm. Considerar: E = 196.20 GPa.
Las barras tienen 6, 9 y 10 cm2, respectivamente.

Docente : Mg. Ing. Evangelos Manouris clase 11 - 2 | 13

Solución-11.2.1
1. Primero colocamos el sistema de ejes global.
2. Numeramos los nodos y miembros de la estructura.
3. Calculamos las coordenadas de cada nodo.

Docente : Mg. Ing. Evangelos Manouris clase 11 - 3 | 13

 Nudo Inicial Nudo Final
Coordenadas Coordenadas (A*E )/ L
Num. Num. Num.
ο
Barra Nudo X (m) Y (m) Nudo X (m) Y (m) L (cm) A (cm²) E(KN/m²) l m l*m l² m² (KN / cm) θ ( )
1 1 -400,00 0,00 4 0,00 400,00 565,69 6,00 19620,00 0,7071 0,7071 0,5000 0,5000 0,5000 208,10 45,00
2 2 -231,00 0,00 4 0,00 400,00 461,91 9,00 19620,00 0,5001 0,8660 0,4331 0,2501 0,7499 382,28 59,99
3 3 0,00 0,00 4 0,00 400,00 400,00 10,00 19620,00 0,0000 1,0000 0,0000 0,0000 1,0000 490,50 90,00

En seguida, a fin de estructurar la MATRIZ DE RIGIDEZ de cada elemento, se apoyará en la

tabla siguiente donde se ha calculado los valores de los componentes que posibilitarán la
construcción de la MATRIZ DE REGIDEZ de cada barra.

Xi Yi
𝑨∗𝑬 𝟐
{𝑲} = [𝒌𝒊𝒋 ] = ∗{ 𝒍 𝒍 ∗ 𝒎} Xj
𝑳 𝒍∗𝒎 𝒎𝟐 Yj
Y la MATRIZ DE RIGIDEZ completa es la siguiente:
𝒍𝟐 𝒍∗𝒎 −𝒍² −𝒍 ∗ 𝒎
[𝒌𝟏,𝟏 ]𝟏,𝟏 [𝒌𝟏,𝟐 ]𝟏,𝟐 𝑨∗𝑬 𝒎𝟐
{𝑲} = [ ]= ∗{ 𝒍∗𝒎 −𝒍 ∗ 𝒎 −𝒎² }
[𝒌𝟐,𝟏 ]𝟐,𝟏 [𝒌𝟐,𝟐 ]𝟐,𝟐 𝑳 −𝒍² −𝒍 ∗ 𝒎 𝒍² 𝒍∗𝒎
−𝒍 ∗ 𝒎 𝒎² 𝒍∗𝒎 𝒎²

Reemplazando los valores obtenidos en la tabla, se tiene para cada barra la MATRIZ DE
RIGIDEZ en sistema global. Además, en la parte superior se ha anotado los grados de libertad
a que corresponde cada fila y columna. Asimismo, se observa la división de las SUBMATRICES
DE RIGIDEZ de cada barra.

Para Barra 1
U 1x U 1y U 4x U 4y
104,05 104,05 ⁞ -104,05 -104,05 F 1x
[ K1,1 ]1 ⁞ [ K1,2 ]1 104,05 104,05 ⁞ -104,05 -104,05 F 1y
[ K ]1 = ―――― ⁞ ―――― = ―――― ―――― ⁞ ―――― ―――― KN / cm
[ K2,1 ]1 ⁞ [ K2,2 ]1 -104,05 -104,05 ⁞ 104,05 104,05 F 4x
-104,05 -104,05 ⁞ 104,05 104,05 F 4y

Docente : Mg. Ing. Evangelos Manouris clase 11 - 4 | 13

 Para Barra 2
U 2x U 2y U 4x U 4y
95,61 165,55 ⁞ -95,61 -165,55 F 2x
[ K1,1 ]2 ⁞ [ K1,2 ]2 165,55 286,67 ⁞ -165,55 -286,67 F 2y
[ K ]2 = ―――― ⁞ ―――― = ―――― ―――― ⁞ ―――― ―――― KN / cm
[ K2,1 ]2 ⁞ [ K2,2 ]2 -95,61 -165,55 ⁞ 95,61 165,55 F 4x
-165,55 -286,67 ⁞ 165,55 286,67 F 4y

Para Barra 3
U 3x U 3y U 4x U 4y
0,00 0,00 ⁞ 0,00 0,00 F 3x
[ K1,1 ]3 ⁞ [ K1,2 ]3 0,00 490,50 ⁞ 0,00 -490,50 F 3y
[ K ]3 = ―――― ⁞ ―――― = ―――― ―――― ⁞ ―――― ―――― KN / cm
[ K2,1 ]3 ⁞ [ K2,2 ]3 0,00 0,00 ⁞ 0,00 0,00 F 4x
0,00 -490,50 ⁞ 0,00 490,50 F 4y

Para ensamblar las matrices donde existen quatro nudos, se formará una MATRIZ DE
RIGIDEZ de 4 × 4 submatrices, donde cada uno de sus componentes es el siguiente:

U U1 U2 U3 U4
[𝐾1,1 ] [𝐾1,2 ] [𝐾1,3 ] [𝐾1,4 ] F1
[𝐾2,1 ] [𝐾2.2 ] [𝐾2,3 ] [𝐾2,4 ] F2
[𝐾 ] =
[𝐾3,1 ] [𝐾3,2 ] [𝐾3,3 ] [𝐾3,4 ] F3
[[𝐾4,1 ] [𝐾4,2 ] [𝐾4,3 ] [𝐾4,4 ]] F4
Cabe mencionar que la MATRIZ DE RIGIDEZ es simétrica; por lo tanto, es necesario hallar
solo los valores de la MATRIZ SUPERIOR , siendo la matriz inferior (color gris) simétrica a la
parte superior (color negro).
Para hallar los valores de las submatrices, se suman todas las submatrices en la que existe
coincidencia de fila y columna en cada MATRIZ DE RIGIDEZ. De esta manera se tiene lo
siguiente:

Docente : Mg. Ing. Evangelos Manouris clase 11 - 5 | 13

 En la fila 1: ( F1 )
(F1-U1)→ [K 1,1] = [K 1,1]1
(F1-U2)→ [K 1,2] = 0 , (F1-U3)→ [K 1,3] = 0,
(F1-U4)→ [K 1,4] = [K 1,2]1,

En la fila 2: ( F2 )
(F2-U2)→ [K 2,2] = [K 1,1]2
(F2-u3)→ [K 2,3] = 0
(F2-U4)→ [K 2,4] = [K 1,2]2,

En la fila 3: ( F3 )
(F3-u3)→ [K 3,3] = [K 1,1]3
(F3-U4)→ [K 3,4] = [K 1,2]3,

En la fila 4: ( F4 )
(F4-u4)→ [K 4,4] = [K 2,2]1 + [K 2,2]2 + [K 2,2]3
Siendo la MATRIZ ENSAMBLADA DE RIGIDEZ la siguiente: { K } = KN/cm *

U_1 U_2 U_3 U_4

U1x U1y U2x U2y U3x U3y U4x U4y
F1x 104,05 104,05 0,00 0,00 0,00 0,00 -104,05 -104,05
F_1
F1y 104,05 104,05 0,00 0,00 0,00 0,00 -104,05 -104,05

F2x 0,00 0,00 95,61 165,55 0,00 0,00 -95,61 -165,55

F_2
F2y 0,00 0,00 165,55 286,67 0,00 0,00 -165,55 -286,67

F3x 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00

F_3
F3y 0,00 0,00 0,00 0,00 0,00 490,50 0,00 -490,50

F4x -104,05 -104,05 -95,61 -165,55 0,00 0,00 199,66 269,61

F_4
F4y -104,05 -104,05 -165,55 -286,67 0,00 -490,50 269,61 881,23
Submatriz correspondiente al nudo libre

Docente : Mg. Ing. Evangelos Manouris clase 11 - 6 | 13

Aplicando EQUACION FUNDAMENTAL DE METODO MATRICIAL DE RIGIDEZ. {F} = {K} * {U}
se obtiene la siguiente expresión:

- -
F1x 104,05 104,05 0,00 0,00 0,00 0,00
104,05 104,05
0 U1x
- -
F1y 104,05 104,05 0,00 0,00 0,00 0,00
104,05 104,05
0 U1y
-
F2x 0,00 0,00 95,61 165,55 0,00 0,00 -95,61
165,55
0 U2x
- -
F2y 0,00 0,00 165,55 286,67 0,00 0,00 0 U2y
F3x
= 0,00 0,00 0,00 0,00 0,00 0,00
165,55 286,67
0,00 0,00
* 0 U3x
-
F3y 0,00 0,00 0,00 0,00 0,00 490,50 0,00
490,50
0,4 U3y
- - -
22,49 KN
104,05 104,05
-95,61
165,55
0,00 0,00 199,66 269,61 U4x U4x
- - - - -
-31.47 KN
104,05 104,05 165,55 286,67
0,00
490,50
269,61 881,23 U4y U4y

Se efectúa la traslación de la submatriz correspondiente al nodo libre hasta ubicarlo en la

primera fila y columna para su respectiva solución. Asimismo, se nombró a cada una de las
submatrices.
Docente : Mg. Ing. Evangelos Manouris clase 11 - 7 | 13
 U4x U3x U1x U1y U2x U2y U3x U3y
- -
F1x 104,05 104,05 0,00 0,00 0,00 0,00
104,05 104,05
0 U1x
- -
F1y 104,05 104,05 0,00 0,00 0,00 0,00
104,05 104,05
0 U1y
-
F2x 0,00 0,00 95,61 165,55 0,00 0,00 -95,61
165,55
0 U2x
- - U2y
F2y 0,00 0,00 165,55 286,67 0,00 0,00 0

F3x
= 0,00 0,00 0,00 0,00 0,00 0,00
165,55 286,67
0,00 0,00
* 0 U3x
-
F3y 0,00 0,00 0,00 0,00 0,00 490,50 0,00
490,50
0,4 U3y
- - -
22,49
104,05 104,05
-95,61
165,55
0,00 0,00 199,66 269,61 U4x U4x
- - - - -
-31.47
104,05 104,05 165,55 286,67
0,00
490,50
269,61 881,23 U4y U4y

U4x U3x U1x U1y U2x U2y U3x U3y

- -
F1x 104,05 104,05 0,00 0,00 0,00 0,00
104,05 104,05
0 U1x
- -
F1y 104,05 104,05 0,00 0,00 0,00 0,00
104,05 104,05
0 U1y
-
F2x 0,00 0,00 95,61 165,55 0,00 0,00 -95,61
165,55
0 U2x
- -
F2y 0,00 0,00 165,55 286,67 0,00 0,00 0 U2y
F3x
= 0,00 0,00 0,00 0,00 0,00 0,00
165,55 286,67
0,00 0,00
* 0 U3x
-
F3y 0,00 0,00 0,00 0,00 0,00 490,50 0,00
490,50
0,4 U3y
- - -
22.49
104,05 104,05
-95,61
165,55
0,00 0,00 199,66 269,61 U4x U4x
- - - - -
-31.47
104,05 104,05 165,55 286,67
0,00
490,50
269,61 881,23 U4y U4y

Docente : Mg. Ing. Evangelos Manouris clase 11 - 8 | 13

 K11 K12 U11
F11
U4x U3x U1x U1y U2x U2y U3x U3y
22,49 199,66 269,61 -104,05 -104,05 -95,61 -165,55 0,00 0,00 U4x U4x
31,47 269,61 881,23 -104,05 -104,05 -165,55 -286,67 0,00 -490,50 U4y U4y
F1x -104,05 -104,05 104,05 104,05 0,00 0,00 0,00 0,00 0 U1x
F1y -104,05 -104,05 104,05 104,05 0,00 0,00 0,00 0,00 0 U1y
F2x -95,61 -165,55 0,00 0,00 95,61 165,55 0,00 0,00 0 U2x
F2y
= -165,55 -286,67 0,00 0,00 165,55 286,67 0,00 0,00 * 0 U2y
F3x 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0 U3x
F3y 0,00 -490,50 0,00 0,00 0,00 0,00 0,00 490,50 -0,4 U3y

F21 K21 K22 U21

𝑭 𝑲 𝑲𝟏𝟐 𝑼
{ 𝟏𝟏 } = { 𝟏𝟏 } ∗ { 𝟏𝟏 }
𝑭𝟏𝟐 𝑲𝟐𝟏 𝑲𝟐𝟐 𝑼𝟐𝟏
Expandiendo la expresión anterior se obtiene:

{ F11 } = { K11 } * { U11 } + { K12 } * { U21 } EC. 1

{ F21 } = { K21 } * { U11 } + { K22 } * { U21 } EC. 2
Aplicando la ecuación (1) se obtiene lo siguiente:

22,49 199,66 269,61 U4x

-31,47
= 269,61 881,23
* U4y
+
0
0
-104,05 -104,05 -95,61 -165,55 0,00 0,00 0
+ -104,05 -104,05 -165,55 -286,67 0,00 -490,50
* 0
0
-0,4
Docente : Mg. Ing. Evangelos Manouris clase 11 - 9 | 13
 22,49 199,66 269,61 U4x 0,00
= * +
-31.47 269,61 881,23 U4y 196,20

22.49 0 199,66 269,61 U4x

- = *
-31.47 196,20KN 269,61 881,23 U4y

22,49 199,66 269,61 U4x

= *
-227.67 269,61 881,23 U4y

-1
U4x 199,66 269,61 22.49
= *
U4y 269,61 881,23 -227.67

U4x 0,00853 -0,00261 22,49 U4x 0.78638

= * → =
U4y -0,00261 0,00193 -227.67 U4y -0.49895

Para el cálculo de reacciones en los apoyos, se hace uso de la ecuación obtenida

anteriormente:

{ F21 } = { K21 } * { U11 } + { K22 } * { U21 } EC. 2

Docente : Mg. Ing. Evangelos Manouris clase 11 - 10 | 13

 F1x -104,05 -104,05 104,05 104,05 0,00 0,00 0,00 0,00 0
F1y -104,05 -104,05 104,05 104,05 0,00 0,00 0,00 0,00 0
F2x -95,61 -165,55 0,78638 0,00 0,00 95,61 165,55 0,00 0,00 0
F2y = -165,55 -286,67 * -0,49895 + 0,00 0,00 165,55 286,67 0,00 0,00 * 0
F3x 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0
F3y 0,00 -490,50 0,00 0,00 0,00 0,00 0,00 490,50 -0,4

F1x -29,908 0,000 -29,91

F1y -29,908 0,000 -29,91
F2x 7,418 0,000 7,42
F2y = 12,845
+ 0,000
= 12,85
F3x 0,000 0,000 0,00
F3y 244,733 -196,200 48,53

Cálculo de los esfuerzos axiales en cada barra.

Para el cálculo de los esfuerzos axiales en cada barra se aplica la siguiente expresión:
Num. Num. Num. U Inic X U Inic Y U Fin X U Fin Y E
L (cm) l m -l -m Α cm²
Barra Inicial Final cm cm cm cm KN/c
1 1 4 565,69 0,7071 0,7071 -0,7071 -0,7071 0,0000 0,0000 0,7864 -0,4989 6 19620
2 2 4 461,91 0,50 0,87 -0,5001 -0,8660 0,0000 0,0000 0,7864 -0,4989 9 19620
3 3 4 400,00 0,00 1,00 0,0000 -1,0000 0,0000 -0,4000 0,7864 -0,4989 10 19620

Luego, empleando la expresión obtenida para cada barra (ec. 5.3 ):

(𝑼𝟏 )𝒙
𝟏 (𝑼𝟏 )𝒚
𝜺= ∗ {− 𝐜𝐨𝐬 𝜽 − 𝐬𝐢𝐧 𝜽 𝐜𝐨𝐬 𝜽 𝐬𝐢𝐧 𝜽} ∗
𝑳 (𝑼𝟐 )𝒙
{(𝑼𝟐 )𝒚 }
La deformación unitaria es la siguiente:

Docente : Mg. Ing. Evangelos Manouris clase 11 - 11 | 13

En la Barra 1 ( U1, U2 )
(𝑼𝟏 )𝒙
𝟏 (𝑼𝟏)𝒚
𝜺𝟏 = ∗ {−𝒍𝟏 −𝒎𝟏𝟐 𝒍𝟏𝟐 𝒎𝟏} ∗
𝑳𝟏𝟐 (𝑼𝟒 )𝒙
{(𝑼𝟒)𝒚 }

0,0000
ε1= (1/565,69cm) * -0,7071 -0,7071 0,7071 0,7071 * 0,0000 cm = 3,593E-04
0,7864
-0,4989
Esfuerzo unitario axial:
σ1 = ε1 * Ε = 3.593*10-4 * 19620 ΚΝ / cm² = 7,0494 KN/cm² = σ1 .
Ν1 = σ1 * Α1 = 7,0494 ΚΝ/cm² * 6 cm² = 42,29 KN = N1.

En la Barra 2 ( U2, U3 )
(𝑼𝟐)𝒙
𝟏 (𝑼𝟐 )𝒚
𝜺𝟐 = ∗ {−𝒍𝟐 −𝒎𝟐 𝒍𝟐 𝒎𝟐 } ∗
𝑳𝟐 (𝑼𝟒)𝒙
{(𝑼𝟒 )𝒚 }

0,0000
ε2= (1/461,91cm) * -0,5001 -0,8660 0,5001 0,8660 * 0,0000 cm = -8,401E-05
0,7864
-0,4989

Esfuerzo unitario axial:

σ2 = ε2 * Ε = -8,401E-05 * 19620 ΚΝ / cm² = -1,6482 KN/cm² = σ2.
Ν2 = σ2 * Α2 = -1,6482ΚΝ/cm² * 9 cm² = -14,83 KN = N2.

Docente : Mg. Ing. Evangelos Manouris clase 11 - 12 | 13

En la Barra 3( U2, U3 )
(𝑼𝟑)𝒙
𝟏 (𝑼𝟑 )𝒚
𝜺𝟑 = ∗ {−𝒍𝟑 −𝒎𝟑 𝒍𝟑 𝒎𝟑 } ∗
𝑳𝟑 (𝑼𝟒)𝒙
{(𝑼𝟒 )𝒚 }

0,0000
ε3= (1/400cm) * 0,0000 -1,0000 0,0000 1,0000 * -0,4000 cm = -0,000247
0,7864

-0,4989
Esfuerzo unitario axial:
σ3 = ε3 * Ε = -0,000247 * 19620 ΚΝ / cm² = -4,8533 KN/cm² = σ2.
Ν3 = σ3 * Α3 = -4,8533 ΚΝ/cm² * 10 cm² = -48,53 KN = N3.

Docente : Mg. Ing. Evangelos Manouris clase 11 - 13 | 13