346 Curso de análisis estructural

además se emplean métodos numéricos adecuados para la resolución del sistema de

ecuaciones.

8.21 EJERCICIOS RESUELTOS

8.21.1 Calcular los esfuerzos en la estructura de la figura 8.33. Todas las barras
son del mismo material, con módulo de elasticidad E=200000 MPa y tienen
igual área transversal A=10 cm2.

8 kN 5 kN 2y 1y
2 B
B
1 2x
1x

E
E
A
C A C
10 m

D D
3y
4y
3 4
10 m 3x 4x

Figura 8.33

Los grados de libertad de la estructura son:

T T
1x 1y 2x 2y 3x 3y 4x 4y I F

Todos los elementos son barras articuladas en los dos extremos, cuya matriz de rigidez
viene dada por la ecuación (7.39).
Elemento A: Nudo inicial: 3. Nudo final: 2. Ángulo: 90º. Longitud: 10 m. Área: 0.001 m2.
cos =0. sin =1.
0 0 0 0 3x

EA 0 1 0 1 3y
K GA
10 0 0 0 0 2x
0 1 0 1 2y

Elemento B: Nudo inicial: 2. Nudo final: 2. Ángulo: 0º. Longitud: 10 m. Área: 0.001 m2.
cos =1. sin =0.
1 0 1 0 2x

EA 0 0 0 0 2y
K GB
10 1 0 1 0 1x
0 0 0 0 1y

Elemento C: Nudo inicial: 4. Nudo final: 1. Ángulo: 90º. Longitud: 10 m. Área: 0.001 m2.
cos =0. sin =1.
Método de rigidez 347

0 0 0 0 4x

EA 0 1 0 1 4y
K CG
10 0 0 0 0 1x
0 1 0 1 1y

Elemento D: Nudo inicial: 3. Nudo final: 1. Ángulo: 45º. Longitud: 10 2 m.

Área: 0.001 m2. cos =1/ 2 . sin =1/ 2 .
1/ 2 1/ 2 1/2 1/ 2 3x

EA 1/ 2 1/ 2 1/2 1/ 2 3y
K GD
10 2 1/ 2 1/ 2 1/2 1/ 2 1x
1/ 2 1/ 2 1/2 1/ 2 1y

Elemento E: Nudo inicial: 2. Nudo final: 4. Ángulo: -45º. Longitud: 10 2 m.

Área: 0.001 m2. cos =1/ 2 , sin =-1/ 2 .
1/ 2 1/ 2 1/2 1/ 2 2x

EA 1/ 2 1/ 2 1/2 1/ 2 2y
K GE
10 2 1/ 2 1/ 2 1/2 1/ 2 4x
1/ 2 1/ 2 1/2 1/ 2 4y

Matrices de ordenación de las barras:

0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0
0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0
BA BB
0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0

0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0
0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0
BC BD
1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0

0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
BE
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1

La contribución de cada elemento a la matriz de rigidez de la estructura es:

T
K (Ge) Be K eC
G B
e
e A, B, C, D, E
A continuación se muestran las distintas matrices:
348 Curso de análisis estructural

0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0
0 0 0 0 0 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 0
EA 0 0 0 1 0 1 0 0 EA 0 0 0 0 0 0 0 0
K (GA ) K (GB)
10 0 0 0 0 0 0 0 0 10 0 0 0 0 0 0 0 0
0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 1/ 2 1/ 2 0 0 1/ 2 1/ 2 0 0
0 1 0 0 0 0 0 1 1/ 2 1/ 2 0 0 1/ 2 1/ 2 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
EA 0 0 0 0 0 0 0 0 EA 0 0 0 0 0 0 0 0
K (GC ) K (GD)
10 0 0 0 0 0 0 0 0 10 2 1/ 2 1/ 2 0 0 1/2 1/ 2 0 0
0 0 0 0 0 0 0 0 1/ 2 1/ 2 0 0 1/2 1/ 2 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 1/2 1/ 2 0 0 1/ 2 1/ 2
EA 0 0 1/2 1/2 0 0 1/ 2 1/ 2
K (GE )
10 2 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 1/2 1/2 0 0 1/ 2 1/ 2
0 0 1/2 1/ 2 0 0 1/ 2 1/ 2

La matriz de rigidez de la estructura se obtiene sumando la colaboración de los distintos

elementos.

0.1354 0.0354 0.1 0 0.0354 0.0354 0 0

0.0354 0.1354 0 0 0.0354 0.0354 0 0.1
0.1 0 0.1354 0.0354 0 0 0.0354 0.0354
0 0 0.0354 0.1354 0 0.1 0.0354 0.0354
K EA
0.0354 0.0354 0 0 0.0354 0.0354 0 0
0.0354 0.0354 0 0.1 0.0354 0.1354 0 0
0 0 0.0354 0.0354 0 0 0 .0354 0.0354
0 0.1 0.0354 0.0354 0 0 0.0354 0.1354

La matriz KII se obtiene a partir de la matriz de rigidez completa, extrayendo los

coeficientes correspondientes a los grados de libertad no nulos, que son los cuatro primeros:
Método de rigidez 349

2.7071 0.7071 2.0 0

0.7071 2.7071 0 0
K II 10 7
2.0 0 2.7071 0.7071
0 0 0.7071 2.7071

El vector de fuerzas nodales se obtiene de forma directa a partir de las fuerzas exteriores:
F1x 0
F1y 5000
F2 x 8000
FI F2 y 0
F
FF F3 x 0
F3 y 0
F4 x 0
F4 y 0

El sistema de ecuaciones de equilibrio, después de aplicar las condiciones de ligadura de

desplazamiento conocido, es:
2.7071 0.7071 2.0 0 1x 0
0.7071 2.7071 0 0 1y 5000
K II I FI 10 7
2.0 0 2.7071 0.7071 2x 8000
0 0 0.7071 2.7071 2y 0

Los desplazamientos de los nudos son:

1x 0.8167
1y 0.3980 3
I 10 m
2x 0.9647
2y 0.2520

Las reacciones en los apoyos se obtienen de la ecuación (8.94):

RF K FI I FF
R1x 0.7071 0.7071 0 0 0.8167 0 2960
R1y 0.7071 0.7071 0 2.0 0.3980 0 8000
10 7 10 3
R2 x 0 0 0.7071 0.7071 0.9647 0 5040
R2 y 0 2.0 0.7071 0.7071 0.2520 0 13000

Esfuerzos en el elemento A. Las deformaciones de los nudos de este elemento son:

3x 0.0
A 3y 0.0
BA 10 3

2x 0.9647
2y 0.2520
350 Curso de análisis estructural

3x 0 1 0 0 0. 0.
A 3y A 1 0 0 0 0. 3 0. 3
T2 10 10
2x 0 0 0 1 0.9647 0.2520
2y 0 0 1 0 0.2520 0.9647

con lo que los esfuerzos en sus extremos, calculados en el sistema local, son:
P3 x 1 0 1 0 0. 5039
P3 y EA 0 0 0 0 0. 0
PA K LA A
10 3
P2 x 10 0 0 1 0 0.2520 5039
P2 y 0 0 0 0 0.9647 0

que corresponden a un esfuerzo axial de tracción de 5039 N, como se muestra en la figura

8.34. También pueden calcularse los esfuerzos en el sistema general, empleando las
deformaciones y la matriz de rigidez referidas a dicho sistema:
F3 x 0 0 0 0 0. 0
F3 y EA 0 1 0 1 0. 5039
FA K GA A
10 3
F2 x 10 0 0 0 0 0.9647 0
F2 y 0 1 0 1 0.2520 5039

2y =0.2520 10-3
5039
2x =0.9647 10-3

2y=0.2520 10-3 1y=-0.398 10-3

B
A A
2x=0.9647 10 -3
1x=0.8167 10-3
YL

-5039 -2960 B +2960

YL

Figura 8.34

Esfuerzos en el elemento B. Al estar orientado según el eje X, su sistema local coincide con
el general de la estructura, y las deformaciones de los nudos son (figura 8.34):

2x 0.9647
B 2y 0.2520
BB B
10 3
1x 0.8167
1y 0.3980

con lo que los esfuerzos en sus extremos son:

Método de rigidez 351

P2 x 1 0 1 0 0.9647 2960
P2 y EA 0 0 0 0 0.2520 0
PB K BL B
10 3
P1x 10 0 0 1 0 0.8167 2960
P1y 0 0 0 0 0.3980 0

que corresponden a un esfuerzo axial de compresión de 2960 N.

Esfuerzos en el elemento D (figura 8.35). Las deformaciones de los nudos de este elemento,
en los sistemas general y local, son:

3x 0.0
D 3y 0.0
BD 10 3
1x 0.8167
1y 0.3980

3x 1 1 0 0 0. 0.
D 3y D 1 1 1 0 0 0. 3 0. 3
T2 10 10
1x 2 0 0 1 1 0.8167 0.2960
1y 0 0 1 1 0.3980 0.8589

con lo que los esfuerzos en sus extremos, calculados en el sistema local, son:
P3 x 1 0 1 0 0. 4186
P3 y EA 0 0 0 0 0. 0
PA K LA A
10 3
P2 x 10 2 0 0 1 0 0.2960 4186
P2 y 0 0 0 0 0.8589 0

que corresponden a un esfuerzo axial de tracción de 4186 N. También pueden calcularse los
esfuerzos en el sistema general, empleando las deformaciones y la matriz de rigidez
referidas a dicho sistema:
F3 x 0.5 0.5 0.5 0.5 0. 2960
F3 y EA 0.5 0.5 0.5 0.5 0. 2960
FD K GD D
10 3
F1x 10 2 0.5 0.5 0.5 0.5 0.8167 2960
F1y 0.5 0.5 0.5 0.5 0.3980 2960

1y=-0.398 10-3 1y=-0.859 10-3

1x=0.8167 10-3 1x=0.296 10 -3

D D
YL

3y=0 45º 3y =0

3x=0 3x=0
Figura 8.35
352 Curso de análisis estructural

2960
+4186
2960

D D

-2960 YL
-4186
45º

-2960
Figura 8.35 (cont)
Los esfuerzos en las restantes barras se calculan siguiendo el mismo procedimiento. Los
valores finales se indican en la figura 8.36.

5000
8000 -2960

+4186

+5039 -7960
-7127

2960
5040
8000 13000
Figura 8.36

8.21.2 Calcular los esfuerzos en el pórtico de la figura 8.37. Todas las barras son
del mismo material, con módulo de elasticidad E=200000 MPa. Sus
propiedades resistentes se indican en la figura.

3000 N/m
5 kN
2
1
A= 80 cm2
I=15000 cm 4
A= 100 cm 2
A= 100 cm 2 I=10000 cm4
4m I=10000 cm 4

3 4
4m 2m
Figura 8.37
Para facilitar el desarrollo se emplean únicamente los grados de libertad no nulos, que son
los siguientes:
T
I 1x 1y 1 2x 2y 2

Las unidades empleadas en la resolución son N y m.

Método de rigidez 353

1y 2y
1x B 2x

1
2

C
A

3 4
Figura 8.38
La matriz de rigidez correspondiente a estos grados de libertad es:

K GA11 K GB11 K GB12

K II
K GB 21 K GB 22 K CG 22

Todos los elementos son barras empotradas en los dos extremos, cuya matriz de rigidez
viene dada por la ecuación (7.58).
Elemento A. Nudo inicial: 3. Nudo final: 1. Ángulo: 90º. Longitud: 4 m. Área: 0.01 m2.
Inercia: 10-4 m4. cos =0, sin =1.
0.0375 0 0.075 1x
8
K GA11 10 0 5.0 0 1Y
0.075 0 0.2 1

Elemento B. Nudo inicial: 1. Nudo final: 2. Ángulo: 0º. Longitud: 4 m. Área: 0.008 m2.
Inercia: 1.5 10-4 m4. cos =1, sin =0.
4 0 0 4 0 0 1x
0 0.0562 0.1125 0 0.0562 0.1125 1y
K GB11 K GB12 0 0.1125 0.3 0 0.1125 0.15 1
K GB 10 8
K GB 21 K GB22 4 0 0 4 0 0 2x
0 0.0562 0.1125 0 0.0562 0.1125 2y
0 0.1125 0.15 0 0.1125 0.3 2

Elemento C. Nudo inicial: 2. Nudo final: 4. Ángulo: -63.43º. Longitud: 4.47 m.

Área: 0.01 m2. Inercia: 10-4 m4. sin =-0.894, cos =0.447.
0.9159 1.7781 0.0537 2x
8
K CG 22 10 1.7781 3.5831 0.0268 2y
0.0537 0.0268 0.1789 2

La matriz de rigidez de la estructura se obtiene ensamblando las matrices de los distintos

elementos, según los grados de libertad correspondientes, siguiendo el esquema indicado
antes.
354 Curso de análisis estructural

4.0375 0. 0.075 4. 0. 0.
0. 5.0562 0.1125 0. 0.0562 0.1125
0.075 0.1125 0.5 0. 0.1125 0.15
K II 10 8
4. 0. 0. 4.9159 1.7781 0.0537
0. 0.0562 0.1125 1.7781 3.6393 0.0857
0. 0.1125 0.15 0.0537 0.0857 0.4789

Fuerzas nodales debidas a la carga distribuida: las fuerzas y momentos de empotramiento

perfecto en la barra B son (figura 8.39)
qL qL2
P10y 6000 N M10 4000 Nm
2 12
qL qL2
P20y 6000 N M20 4000 Nm
2 12
YL
q=3000 N/m
qL 2
qL 2 =-4000
=4000 12
12 XL

qL /2=6000 qL /2=6000
Figura 8.39
El vector de fuerzas de fase 0 correspondiente a esta fuerza distribuida es:
P1x 0
P1y 6000
M1 4000
P0 B F0B
P2 x 0
P2 y 6000
M2 4000

Este vector se ensambla, con signo menos, en el vector de fuerzas nodales de la estructura.

Incluyendo la carga puntual de 5 kN 6000 6000

aplicada en el nudo 1 en dirección X, el 5000
vector de fuerzas nodales completo es:
4000 4000
F1x 5000
F1y 6000
M1 4000
FI
F2 x 0
F2 y 6000
M2 4000 Figura 8.40

Los desplazamientos de los nudos que se obtienen son:

Método de rigidez 355

1x 0.2621
1y 0.0104
1 0.1286
I K II1FI I 10 3

2x 0.2496
2y 0.1041
2 0.1169

Esfuerzos en el elemento A. Las deformaciones de los nudos de este elemento son:

3x 0.0
3y 0.0
A 3 0.0 3
10
1x 0.2621
1y 0.0104
1 0.1286

3x 0 1 0 0 0 0 0. 0.
3y 1 0 0 0 0 0 0. 0.
A 3 A 0 0 1 0 0 0 0. 3 0. 3
T2 10 10
1x 0 0 0 0 1 0 0.2621 0.0104
1y 0 0 0 1 0 0 0.0104 0.2621
1 0 0 0 0 0 1 0.1286 0.1286

Los esfuerzos en sus extremos, calculados en el sistema local, son:

PA K LA A

P3 x 5 0 0 5 0 0 0. 5224
P3 y 0 0.0375 0.075 0 0.0375 0.075 0. 18
M3 0 0.075 0.2 0 0.075 0.1 0. 679
10 8 10 3
P1x 5 0 0 5 0 0 0.0104 5224
P1y 0.0375 0.075 0 0.0375 0.075 0.2621 18
M1 0.075 0.1 0 0.075 0.2 0.1286 606

Esfuerzos en el elemento B: al estar orientado según el eje X, su sistema local coincide con
el general de la estructura. Las deformaciones de sus nudos son:

1x 0.2621
1y 0.0104
B B 1 0.1286 3
10
2x 0.2496
2y 0.1041
2 0.1169
356 Curso de análisis estructural

Los esfuerzos en sus extremos se obtienen sumando los de la fase 0 con los de la fase 1:
PB P 0B P 1B P0B K BL B

P1x 0 4 0 0 4 0 0 0.2621
P1y 6000 0 0.0562 0.1125 0 0.0562 0.1125 0.0104
M1 4000 0 0.1125 0.3 0 0.1125 0.15 0.1286
10 8 10 3
P2 x 0 4 0 0 4 0 0 0.2496
P2 y 6000 0 0.0562 0.1125 0 0.0562 0.1125 0.1041
M2 4000 0 0.1125 0.15 0 0.1125 0.3 0.1169

P1x 0 4981 4981

P1y 6000 776 5224
M1 4000 3394 606
PB
P2 x 0 4981 4981
P2 y 6000 776 6776
M2 4000 290 3710

Los esfuerzos internos en la viga son asimismo la suma de los debidos a las fases 0 y 1.
Como ejemplo, en la figura 8.41 se muestra el diagrama de momentos de esta viga,
representado con el criterio ingenieril.
4000 3710
4000

606

2000 1842 290

3394 3842

M0 M1 M

Figura 8.41

Esfuerzos en el elemento C. Las deformaciones de los nudos de este elemento en el sistema

general son:
2x 0.2496
2y 0.1041
C 2 0.1169 3
10
4x 0.
4y 0.
4 0.

Proyectándolas sobre el sistema local se obtiene:

Método de rigidez 357

2x

2y

C 2
T2C C

4x

4y

2x 0.447 0.894 0 0 0 0 0.2496 0.0185

2y 0.894 0.447 0 0 0 0 0.1041 0.2698
2 0 0 1 0 0 0 0.1169 3 0.1169 3
10 10
4x 0 0 0 0.447 0.894 0 0. 0.
4y 0 0 0 0.894 0.447 0 0. 0.
4 0 0 0 0 0 1 0. 0.

Al no haber cargas aplicadas sobre esta viga, los esfuerzos en sus extremos, calculados en
el sistema local, son:
PC K CL C

P2 x 4.472 0 0 4.472 0 0 0.0185 8288

P2 y 0 0.0268 0.06 0 0.0268 0.06 0.2698 1425
M2 0 0.06 0.1789 0 0.06 0.089 0.1169 3710
10 8 10 3
P4 x 4.472 0 0 4.472 0 0 0. 8288
P4 y 0 0.0268 0.06 0 0.0268 0.06 0. 1425
M4 0 0.06 0.089 0 0.06 0.1789 0. 2664

La figura 8.42 muestra los esfuerzos en los extremos de todos los elementos así como el
diagrama de momentos flectores de la estructura.

-606 -5224 5224 YL 6776 1425

3710
4981 XL
-18 -4981 8288
-3710
606 YL
XL

YL
XL

5224 -1425

18 2664
679
-8288

Figura 8.42
358 Curso de análisis estructural

3710
3710
606
606

679 M (N-m) 2664

Figura 8.42 (cont.)

8.21.3 Calcular las deformaciones y esfuerzos producidos en la estructura del

ejercicio 8.21.2 por un calentamiento de +20ºC en la cara superior de la viga
B. Se supone que el canto de dicha viga es h=20 cm, y que el coeficiente de
dilatación térmica del material vale 10-5 ºC-1.
Los grados de libertad de la estructura y los coeficientes de rigidez correspondientes son los
mismos que se emplearon en el ejercicio 8.21.2, siendo necesario calcular únicamente el
vector de fuerzas nodales equivalentes al incremento de temperatura.
Las temperaturas media y gradiente en la viga son:
TS TI 20 0 TS TI 20 0
Tm 10 º C Tg 100 º C / m
2 2 h 0.20
Por lo tanto, el vector de fuerzas nodales equivalentes, según la expresión (8.44) es:
PT 1x 160000
PT 1y 0 YL
MT 1 30000 +20ºC 30000
PT0 B 160000 XL
PT 2 x 160000
-160000
PT 2 y 0 -30000
MT 2 30000
Figura 8.43
El vector de fuerzas de la estructura se obtiene ensamblando, con signo menos, las fuerzas
de fase 0 existentes en la barra B.
F1x 160000
F1y 0
M1 30000
FI
F2 x 160000
F2 y 0
M2 30000
Método de rigidez 359

La figura 8.44 muestra dicho vector de fuerzas, que es equivalente al calentamiento de la

barra B.
30000 + 20ºC 30000
160000 160000

Figura 8.44
Las deformaciones nodales que se
producen son (figura 8.45): 0.9207 10-3
0.9364 10-3
1x 0.3239
1y 0.0002
1 0.9364 3
I 10
2x 0.0906
2y 0.0515
2 0.9207
Figura 8.45

Los esfuerzos en la barra B se calculan siguiendo el proceso habitual, como suma de las
fases 0 y 1:
PB PT0 B K BL B

P1x 160000 4 0 0 4 0 0 0.3239

P1y 0 0 0.0562 0.1125 0 0.0562 0.1125 0.0002
M1 30000 0 0.1125 0.3 0 0.1125 0.15 0.9364
10 8 10 3
P2 x 160000 4 0 0 4 0 0 0.0906
P2 y 0 0.0562 0.1125 0 0.0562 0.1125 0.0515
M2 30000 0.1125 0.15 0 0.1125 0.3 0.9207

P1x 160000 165810 5810

P1y 0 110 110
M1 30000 13700 16300
PB
P2 x 160000 165810 5810
P2 y 0 110 110
M2 30000 14150 15850

La figura 8.46 muestra los esfuerzos en la barra B, su deformada y su diagrama de

momentos flectores.
360 Curso de análisis estructural

YL

-110 16300 15850

110
-5810 XL 5810
0.9364 10 -3
0.92 10-3
-16300 15850

Figura 8.46
Nótese que la barra tiene un momento flector que es prácticamente constante y con las
tracciones en su cara inferior, mientras que la deformada de la viga indica una curvatura
con la concavidad hacia abajo (negativa para los ejes empleados). Esta aparente
discordancia es debida a la presencia de la temperatura, con lo que la relación momento
flector - curvatura viene dada por la ecuación (5.23). Adoptando un valor medio para el
momento de valor 16000, dicha ecuación puede ponerse como:
M EIv EI Tg 16000 EIv 30000 EIv 14000

lo cual explica la curvatura negativa para un momento flector positivo.

Los esfuerzos en las restantes barras se calculan siguiendo el mismo procedimiento, y su
obtención se deja como ejercicio. El diagrama de momentos flectores de la estructura,
representado con el criterio ingenieril, se muestra en la figura 8.47.
16300
16300
15850

15850

6934 M (N-m) 7611

Figura 8.47

8.21.4 Calcular las deformaciones de los nudos y el esfuerzo en el muelle de la

estructura mostrada en la figura 8.48. Las vigas tienen una rigidez a flexión
de valor EI = 10 7 N m2 y la rigidez del muelle es de 105 N/m.

100 N/m 500 N

1 A 2 B 3 C 4

3 x 10 m
Figura 8.48
Se prescinde de los grados de libertad según X, que corresponden a la deformación axial y
que son todos nulos pues están desacoplados de la flexión y no hay cargas axiales.
Método de rigidez 361

Los grados de libertad considerados son:

I 3y

3i

En el nudo 3 se ha considerado como incógnita el giro por la izquierda, pero no el giro por
la derecha, como se indica en la figura 8.49.

3y
2 3i

A B C 4
1

Figura 8.49
Para cada uno de los 4 elementos se calcula su colaboración a la rigidez de toda la
estructura. Todos ellos tienen la misma longitud L y rigidez a flexión EI.
Elemento A. Se trata de una viga empotrada en ambos extremos:
4 EI
0 0 2
L
K (IIA) 0 0 0 3y
0 0 0 3i

Elemento B. Se trata de una viga empotrada en ambos extremos:

4 EI 6 EI 2 EI 2
L L2 L
6 EI 12 EI 6 EI
K (IIB) 3y
L2 L3 L2
2 EI 6 EI 4 EI
L L2 L 3i

Elemento C. Se trata de una viga articulada en el nudo 3 y empotrada en el nudo 4. Su

colaboración a la rigidez de toda la estructura es:
0 0 0 2
3EI
K (IIC ) 0 0 3y
L3
0 0 0 3i

Muelle. Colabora a la rigidez de la deformación vertical del nudo 3:

0 0 0 2
K (IIm) 0 K 0 3y
0 0 0 3i

Sumando las matrices anteriores se obtiene la matriz de rigidez de la estructura:

362 Curso de análisis estructural

4 EI 4 EI 6 EI 2 EI 2
L L L2 L
6 EI 12 EI 3 EI 6 EI
K II K (IIA) K (IIB) K (IIC) K (IIm ) K 3y
L2 L3 L3 L2
2 EI 6 EI 4 EI
L L2 L 3i

El vector de fuerzas nodales equivalentes tiene contribución de las dos cargas aplicadas.
Carga distribuida sobre B. Las fuerzas de fase 0 debidas a esta carga son (figura 8.50):
qL qL2
F20y 500 N M20 833 Nm
2 12
qL qL2
F30y 500 N M30i 833 Nm
2 12
YL
q=100 N/m
qL2
qL2 =-833
=833 12
12 XL

qL/2=500 qL/2=500
Figura 8.50
El vector de fuerzas de fase 0 correspondiente a esta fuerza distribuida es:
F2 y 500
M2 833
F0 B
F3 y 500
M3 i 833

Carga puntual sobre C. Las fuerzas de fase 0 debidas a esta carga son las producidas por
una carga puntual sobre una viga articulada - empotrada (figura 8.51):
5P 11P 3PL
F30y 156 N M30d 0 F40y 344 N M40 937 Nm
16 16 16
P=500 N
YL
3PL
=-937
16

XL
5P/16=156 11P/16=344
Figura 8.51
El vector de fuerzas de fase 0 correspondiente a esta fuerza es:
P3 y 156
M3d 0
F 0C
P4 y 344
M4 937
Método de rigidez 363

El vector de fuerzas nodales equivalentes de toda la estructura se obtiene ensamblando los

dos anteriores, y se muestra en la figura 8.52.
M2 833 833
FI F3 y 500 156 656
M3i 833 833

656
833 833
A B C 4
1

Figura 8.52
Las deformaciones de los nudos son:

2 0.349
I K II1FI 3y 3.97 10 3

3i 0.213

La fuerza en el muelle es:

Fm K 3y 10 5 ( 3.97 10 3 ) 397 N

El giro por la derecha de 3 puede obtenerse como suma del giro en la fase 0 y en la fase 1
en la barra C. El giro en la fase 0 es el giro producido en el apoyo de una viga articulada -
empotrada por una carga puntual situada en su centro (figura 8.53). El valor de este giro
puede calcularse fácilmente aplicando el primer teorema de Mohr entre los puntos 3 y 4.
0 0 10 5 1562.5 7
4 3d 3d 1562.5 2500 1562.5 10 rad
2 EI 2 EI EI
500 N 500 N 1562
3d
937

156 344 156

2500

Figura 8.53

El giro de fase 1 se calcula mediante la expresión (7.116):

1 3 4
3 3y 4
3d 4y 3y 5.95 10 rad
2L 2 2L
Por lo tanto el giro total por la derecha en el nudo 3 es:
0 1 3
3d 3d 3d 0.439 10 rad
364 Curso de análisis estructural

8.21.5 Calcular las deformaciones y los esfuerzos en el pórtico de la figura 8.54

cuando el apoyo 4 sufre un descenso de 2 cm. Todas las barras son del
mismo material, con módulo de elasticidad 200 GPa, y tienen un área de
0.01 m2 y un momento de inercia de 2 10 -4 m4.

1 3
2

8m

4 5
2 cm
3m 3m
Figura 8.54
En el nudo 3 sólo se consideran sus dos desplazamientos como incógnitas, no incluyéndose
los giros. Con esto la estructura tiene 8 grados de libertad desconocidos y un grado de
libertad de valor conocido, como se muestra en la figura 8.55.
Estos grados de libertad se agrupan de 1y 2y 3y
la forma siguiente: 1x 2x 3x

1x B C
1 2
1y

2x A
I D E
I 2y
C
2

3x

3y

4y 4y = -2 cm

Figura 8.55

A continuación se obtienen las matrices de rigidez de los distintos elementos estructurales,

pero sólo se calculan aquellos términos que se corresponden con los grados de libertad.
Elemento A. Se trata de un elemento empotrado en ambos extremos.
Nudo inicial: 4. Nudo final: 1. Ángulo: 90º. Longitud: 8 m. cos =0, sin =1.
Método de rigidez 365

2.5 0 2.5 0 4y

8 0 0.0094 0 0.0375 1x
K GA 10
2.5 0 2.5 0 1y
0 0.0375 0 0.2 1

Elemento B. Se trata de un elemento empotrado en ambos nudos extremos.

Nudo inicial: 1. Nudo final: 2. Ángulo: 0º. Longitud: 3 m. cos =1, sin =0.
6.667 0 0 6.667 0 0 1x
0 0.1778 0.2666 0 0.1778 0.2666 1y
0 0.2666 0.5333 0 0.2666 0.2666 1
K GB 10 8
6.667 0 0 6.667 0 0 2x
0 0.1778 0.2666 0 0.1778 0.2666 2y
0 0.2666 0.2666 0 0.2666 0.5333 2

Elemento C. Se trata de un elemento empotrado en el nudo 2 y articulado en el nudo 3.

Nudo inicial: 2. Nudo final: 3. Ángulo: 0º. Longitud: 3 m. cos =1, sin =0.
6.667 0 0 6.667 0 0 2x
0 0.0444 0.1333 0 0.0444 0 2y

8 0 0.1333 0.4000 0 0.1333 0 2

K CG 10
6.667 0 0 6.667 0 0 3x
0 0.0444 0.1333 0 0.0444 0 3y
0 0 0 0 0 0 3

Elemento D. Se trata de un elemento articulado en ambos extremos.

Nudo inicial: 4. Nudo final: 2. Ángulo: 69.44º. Longitud: 8.54 m. cos =0.351, sin =0.93.
2.0522 0.7696 2.0522 4y
K GD 10 8 0.7696 0.2886 0.7696 2x
2.0522 0.7696 2.0522 2y

Elemento E. Se trata de un elemento empotrado en el nudo 5 y articulado en el nudo 3.

Nudo inicial: 5. Nudo final: 3. Ángulo: 90º. Longitud: 8 m. cos =0, sin =1.
De sus 5 posibles grados de libertad, sólo intervienen en la estructura los correspondientes
al nudo final, que son los de las filas y columnas 4 y 5 de su matriz de rigidez:
0.023 0 3x
K GE 10 8
0 2.5 3y

La matriz de rigidez de la estructura se obtiene ensamblando las matrices de los distintos

elementos, según los grados de libertad correspondientes.
366 Curso de análisis estructural

6.676 0. 0.0375 6.667 0. 0. 0. 0. 0.

0. 2.6778 0.2667 0. 0.1778 0.2667 0. 0. 2.5
0.0375 0.2667 0.7333 0. 0.2667 0.2667 0. 0. 0.
6.667 0. 0. 13.622 0.7696 0. 6.667 0. 0.7696
K 10 8 0. 0.1778 0.2667 0.7696 2.2745 0.1333 0. 0.0444 2.0522
0. 0.2667 0.2667 0. 0.1333 0.9333 0. 0.1333 0.
0. 0. 0. 6.667 0. 0. 6.669 0. 0.
0. 0. 0. 0. 0.0444 0.1333 0. 2.5444 0.
0. 2.5 0. 0.7696 2.0522 0. 0. 0. 4.5522

Donde se identifican las distintas matrices:

K II K IC
K
K CI K CC

Las fuerzas nodales correspondientes a la deformación se obtienen mediante la ecuación

(8.99).
F1x 0. 0.
F1y 2.5 5.0
M1 0. 0.
F2 x 0.7696 1.539 6
FC* K IC C 10 8 0.02 10
F2 y 2.0522 4.104
M2 0. 0.
F3 x 0. 0.
F3 y 0. 0.

Se observa que aparecen fuerzas nodales en todos 1.539 106

aquellos grados de libertad que están
relacionados, a través de la matriz de rigidez, con
5.0 106 4.10 106
el grado de libertad cuyo movimiento está
impuesto. La figura 8.56 muestra las distintas
componentes del vector de fuerzas equivalentes
al desplazamiento.
FC*
Este vector constituye el vector de fuerzas
nodales de la estructura, ya que en este caso no
hay ningún otro tipo de fuerzas aplicadas sobre
ella. 2 cm

Figura 8.56

El sistema de ecuaciones final es:

Método de rigidez 367

6.676 0. 0.0375 6.667 0. 0. 0. 0. 1x 0.

0. 2.6778 0.2667 0. 0.1778 0.2667 0. 0. 1y 5.
0.0375 0.2667 0.7333 0. 0.2667 0.2667 0. 0. 1 0.

8
6.667 0. 0. 13.622 0.7696 0. 6.667 0. 2x 6
1.539
10 10
0. 0.1778 0.2667 0.7696 2.2745 0.1333 0. 0.0444 2y 4.104
0. 0.2667 0.2667 0. 0.1333 0.9333 0. 0.1333 2 0.
0. 0. 0. 6.667 0. 0. 6.669 0. 3x 0.
0. 0. 0. 0. 0.0444 0.1333 0. 2.5444 3y 0.

Los desplazamientos de los nudos son:

1x 1.892
1y 2.008
1 0.247
2x 1.894
I K II1FC* I 10 2
2y 1.273
2 0.320
3x 1.893
3y 0.005

La reacción exterior en la dirección Y del nudo 3, es decir la fuerza que es necesario aplicar
en dicho punto para imponer el descenso, se calcula mediante la expresión:
RC K CI I K CC I FF

1.892
2.008
0.247
1.894
R4 y 10 8 0. 2.5 0. 0.7696 2.0522 0. 0. 0. 10 2
1.273
0.320
1.893
0.005

4.5522 10 8 0.02 13617 N

Esfuerzos en los elementos: como ejemplo se calculan los esfuerzos en la barra A. Las
deformaciones de los nudos de este elemento en los sistemas general y local son:
4x 0.
4y 2.
A 4 0. 2
10
1x 1.892
1y 2.008
1 0.247
368 Curso de análisis estructural

3x 0 1 0 0 0 0 0. 2
3y 1 0 0 0 0 0 2. 0.
A 3 A 0 0 1 0 0 0 0. 2 0. 2
T2 10 10
1x 0 0 0 0 1 0 1.892 2.008
1y 0 0 0 1 0 0 2.008 1.892
1 0 0 0 0 0 1 0.247 0.247

con lo que los esfuerzos en sus extremos, calculados en el sistema local, son:
PA K LA A

P4 x 2.5 0 0 2.5 0 0 2 20766

P4 y 0 0.0094 0.0375 0 0.0094 0.0375 0. 8457
M4 0 0.0375 0.2 0 0.0375 0.1 0. 46208
10 8 10 2
P1x 2.5 0 0 2.5 0 0 2.008 20766
P1y 0 0.0094 0.0375 0 0.0094 0.0375 1.892 8457
M1 0 0.0375 0.1 0 0.0375 0.2 0.247 21445

La figura 8.57 muestra el diagrama de momentos flectores de la estructura y los esfuerzos

axiales en las barras. Nótese que por efecto de la deformación impuesta en el nudo 4 la
barra A está a compresión y la E a tracción. Ambos esfuerzos se equilibran con la reacción
exterior de valor 13617 N.

21445 -8456 4437

40852

+36721
-13617
-20766

46208 35497

13617 M (m-N) N (N)

Figura 8.57

8.21.6 En el pórtico de la figura 8.58 existe un empotramiento deslizante en el nudo

1, orientado a 45º con respecto a la horizontal. Calcular las deformaciones y
los esfuerzos en la barra A, bajo la acción de las cargas indicadas. Todas las
Método de rigidez 369

barras son del mismo material, con módulo de elasticidad 210000 MPa, área
0.01 m2 e inercia 10-4 m4.
1000 N/m 10 kN

2 B 3

C 3m
A

1 4
45º
3m
Figura 8.58
Se adopta un sistema de ejes local en el nudo 1, tal y como se indica en la figura 8.59, con
objeto de imponer con sencillez en él la condición de ligadura introducida por la
deslizadera inclinada. Los grados de libertad de la estructura en el sistema de ejes general y
en el sistema de ejes mixto, considerando el sistema local del nudo 1, son:
N
1x 1x 2y 3y
N
1y 1y
2x 3x
2x N 2x
2y 2y
3x 3x YN
3y 3y 1y

1x

XN

Figura 8.59

La relación entre los grados de libertad del Y

nudo 1 en los sistemas general y local se YN
muestra en la figura 8.60, y es
1Y
N
1Y =0
N
1x cos( 45) sin( 45) 1x
N X
1y sin( 45) cos( 45) 1y
1X

-45º
N
1 T1N 1
N
1X
XN

Figura 8.60

La matriz de rigidez de la estructura en el sistema general es:

370 Curso de análisis estructural

K GA11 K GA12 0
K K GA 21 K GA22 K GB 22 K GB23
0 K GB 32 K GB33 K CG33

donde todas las submatrices son de tamaño 2x2. Esta matriz, expresada en el sistema mixto
de grados de libertad, es:
T
T1N K GA11T1N T1N K GA12 0
T
N
K K GA 21T1N K GA22 K GB 22 K GB 23
0 K GB 32 K GB 33 K CG33

A continuación se obtienen las matrices de rigidez de los distintos elementos, pero sólo se
calculan aquellos términos que se corresponden con los grados de libertad no nulos.
Elemento A. Se trata de un elemento empotrado en el nudo 1 y articulado en el nudo 2. Se
adopta un sistema de ejes local con su origen en el nudo 1. La matriz de rigidez de la barra,
expresada en el sistema de ejes general es:
0.0233 0 0.0233 0 1x

K GA11 K GA12 0 7. 0 7. 1y
K GA 10 8
K GA 21 K GA22 0.0233 0 0.0233 0 2x
0 7. 0 7. 2y

Elemento B. Se trata de un elemento articulado en ambos nudos extremos. Se adopta un

sistema de ejes local con su origen en el nudo 2, con lo que coincide con el general.
7. 0. 7. 0. 2x

K GB 22 K GB 23 0. 0. 0. 0. 2y
K GB 108
K GB 23 K GB33 7 0. 7. 0. 3x
0. 0. 0. 0. 3y

Elemento C. Se trata de un elemento empotrado en el nudo 4 y articulado en el nudo 3. Se

adopta un sistema de ejes local con su origen en el nudo 4.
0.0233 0 3x
K CG 33 10 8
0 7. 3y

La matriz de rigidez de la estructura se obtiene ensamblando las matrices de los distintos

elementos, según los grados de libertad correspondientes a cada uno de ellos:
0.0233 0. 0.0233 0. 0. 0. 1x
0. 7. 0. 7. 0. 0. 1y
0.0233 0. 7.0233 0. 7. 0. 2x
K 10 8
0. 7 0. 7. 0. 0. 2y
0. 0. 7. 0. 7.0233 0. 3x
0. 0. 0. 0. 0. 7. 3y
Método de rigidez 371

Con objeto de transformar esta matriz al sistema de ejes local del nudo 1, se pre y post
multiplican las dos primeras filas y columnas por la matriz de rotación del nudo 1, T1N , con
lo que se obtiene:
N
3.5117 3.488 0.0165 4.9497 0. 0. 1x
N
3.488 3.5117 0.0165 4.9497 0. 0. 1y
0.0165 0.0165 7.0233 0. 7. 0. 2x
KN 10 8
4.9497 4.9497 0. 7. 0. 0. 2y
0. 0. 7. 0. 7.0233 0. 3x
0. 0. 0. 0. 0. 7. 3y

En esta matriz es fácil introducir la condición de ligadura impuesta por la deslizadera, que
es N1y 0 . Para ello basta con borrar la segunda fila y la segunda columna, con lo que se
obtiene la matriz de rigidez final, que está referida a los grados de libertad no nulos y entre
ellos a la deformación en la dirección local de la deslizadera.
N
3.5117 0.0165 4.9497 0. 0. 1x
0.0165 7.0233 0. 7 0. 2x
8
K IIN 10 4.9497 0. 7. 0. 0. 2y
0. 7. 0. 7.0233 0. 3x
0. 0. 0. 0. 7. 3y

El cálculo de las fuerzas nodales sigue el proceso habitual. La carga distribuida aplicada
sobre la barra A genera unas fuerzas de fase 0 cuyo valor es (figura 8.61):
XL
1000 N/m
5qL
P10y 1875 N 3qL/8=1125
8
qL2
M10 1125 Nm
8
3qL
P20y 1125 N
8
M20 0
YL 5qL/8=1875

qL 2
=1125
8

Figura 8.61

El vector de fuerzas de fase 0, referido a los sistemas de ejes local y general, es:
372 Curso de análisis estructural

P1x 0 F1x 1875

P1y 1875 F1y 0
0A M1 0A M1
P 1125 F 1125
P2 x 0 F2 x 1125
P2 y 1125 F2 y 0

La carga puntual en la barra B genera unas YL

10 kN
fuerzas en los extremos de la barra cuyo valor XL
es (figura 8.62):
PL 5 kN 5 kN
P20y 5000 N
2
M0
M20 0
PL
P30y 5000 N
2
M30 0 Figura 8.62

No aparecen momentos en los extremos, ya que la barra está articulada en ellos, pero sin
embargo sí que existe un diagrama de momentos triangular, debido a la carga puntual. El
vector de fuerzas de fase 0 para esta barra, en el sistema general, es:
F2 x 0.
F2 y 5000
F0 B
F3 x 0.
F3 y 5000

El vector de fuerzas nodales de la estructura se obtiene ensamblando, con signo menos, los
vectores de fase 0 producidos por las dos cargas estudiadas, según los grados de libertad
correspondientes.
F1x 1875 0. 1875
F1y 0. 0. 0.
F2 x 1125 0. 1125
F
F2 y 0. 5000 5000
F3 x 0. 0. 0.
F3 y 0. 5000 5000

Este vector se transforma al sistema de grados de libertad mixto, a base de proyectar las
fuerzas aplicadas en el nudo 1 al sistema local de dicho nudo.
F1Nx 1 1 1 F1x F1Nx 1 1 1 1875 1326
F1N T1N F1
F1Ny 2 1 1 F1y F1Ny 2 1 1 0 1326

La figura 8.63 muestra las distintas componentes del vector de fuerzas nodales
equivalentes.
Método de rigidez 373

5000 5000
1125

1326

1875

1326 Figura 8.63

El sistema de ecuaciones final es:
N
3.5117 0.0165 4.9497 0. 0. 1x 1326.
0.0165 7.0233 0. 7 0. 2x 1125.
8
K IIN N
F N
10 4.9497 0. 7. 0. 0. 2y 5000.
0. 7. 0. 7.0233 0. 3x 0.
0. 0. 0. 0. 7. 3y 5000.

Los desplazamientos de los nudos que se obtienen son:

N
1x 0.9032
2x 0.3440
N 2
2y 0.6394 10
3x 0.3429
3y 0.007

Esfuerzos en la barra A: para calcularlos es necesario transformar previamente los

desplazamientos del nudo 1 desde su sistema local al sistema general de la estructura. Esto
se hace empleando la matriz de orientación del sistema local del nudo 1.
T
1 T1N N
1

T N
1x cos( 45) sin( 45) 1x 1 1 1 0.9032 0.6386
N
1y sin( 45) cos( 45) 1y 2 1 1 0. 0.6386

Los valores de los grados de libertad del elemento A, en el sistema general son (figura
8.64):
1x 0.6386
1y 0.6386
A 2
1 0. 10
2x 0.3440
2y 0.6394

Su proyección sobre el sistema local de la barra es:

374 Curso de análisis estructural

1x 0 1 0 0 0 0.6386 0.6386
1y 1 0 0 0 0 0.6386 0.6386
A A 2 2
1 T2 0 0 1 0 0 0. 10 0. 10
2x 0 0 0 0 1 0.3440 0.6394
2y 0 0 0 1 0 0.6394 0.3440

Los esfuerzos en sus extremos, calculados en el sistema local, son:

PA P0A K LA A

P1x 0 7 0 0 7 0 0.6386 5000

P1y 1875 0 0.0233 0.07 0 0.0233 0.6386 5000
M1 8 2
1125 10 0 0.07 0.21 0 0.07 0 10 19500
P2 x 0 7 0 0 7 0 0.6394 5000
P2 y 1125 0 0.0233 0.07 0 0.0233 0.3440 8000

5000
8000
-0.6394 10 -2

0.344 10-2

XL

5000
YL 0.6386 10 -2

0.9032 10 -2 5000
0.6386 10 -2 19500

Figura 8.64 Figura 8.65

Se observa que, en el nudo 1, las dos fuerzas de reacción iguales según los ejes generales
producen una resultante en dirección perpendicular a la deslizadera, mientras que su
proyección en la dirección de la deslizadera es nula (figura 8.65).
El diagrama de momentos flectores en la barra A se obtiene como suma de los diagramas
de las fases 0 y 1, y se muestra en la figura 8.66.
La figura 8.67 muestra el diagrama de momentos de la estructura.
Método de rigidez 375

1000 N/m 1000 N/m

1125 20625 19500

Fase 0 Fase 1 Total

Figura 8.66

7500

19500 24000
Figura 8.67

8.21.7 Calcular las deformaciones de los nudos y los esfuerzos en las barras de la
celosía espacial de la figura 8.68. Todas las barras son del mismo material,
con módulo de elasticidad 200 GPa, y área 0.0005 m2.
4m

4m 50
50
20 40

50
50
40
20 20 8m

20

8m

8m
Figura 8.68