q= 16 KN/m

3 5

Solución
Para este ejercicio, se plantean las ecuaciones de pendiente deflexión para cada tramo:

por las condiciones del problema, las condiciones de contorno dictan que

El giro en el punto A es cero debido a que en el empotramiento la viga es incapaz de rotar. Por otro lado el
momento Mba es igual pero opuesto en signo al momento Mbc. Esta oposición de signos se da para cumplir e
equilibrio de la fórmula. Finalmente, como no existen asentamientos en este ejemplo, los asentamientos
(deformaciones verticales en nudos) en cada tramo son cero.

Las ecuaciones anteriores entonces quedan así:

Además, reemplazando numéricamente el resto de los datos del problema:

Tramo I ( MEP)

si EI = 2000 KN*m2
L1 = 3m
L2 = 5m

q= 16

+ -
L1= 3

Ma = 12 Ma = -12

Tramo II ( MEP)

si EI = 2000 KN*m2
L1 = 3m
L2 = 5m

q= 16

+ -
L2= 5

Ma = 33.33 Ma = - 33.33

Reemplazando todos los datos en el sistema de ecuaciones:

Mab = 2 2000 (qb) + 12

3
Mba = 2 2000 (2*qb) - 12
3

Mbc = 2 2000 (2*qb+qc) + 33.33

5

Mcb = 2 2000 (qb+2*qc) - 33.33

5

Mba + Mbc = 0

Y si escribimos el sistema de cinco ecuaciones con cinco incógnitas de forma matricial, queda:

Mab Mba Mbc Mcb (qb) (qc) =

1 0 0 0 -1333.33333 0 = 12
0 1 0 0 -2666.66667 0 = -12
0 0 1 0 -1600 -800 = 33.33
0 0 0 1 -800 -1600 = - 33.33
0 1 1 0 0 0 = 0

matriz
1 0 0 -1333.33333 0
0 1 0 -2666.66667 0
0 0 1 -1600 -800
0 0 0 -800 -1600
0 1 1 0 0

matriz inversa

1 -0.34482759 -0.34482759 0.17241379 0.34482759 12

0 0.31034483 -0.68965517 0.34482759 0.68965517 -12
0 -0.31034483 0.68965517 -0.34482759 0.31034483 33.33
0 -0.00025862 -0.00025862 0.00012931 0.00025862 - 33.33
0 0.00012931 0.00012931 -0.00068966 -0.00012931 0

Mab #NAME? 0
Mba #NAME? 3
Mbc #NAME? 3
(qb) #NAME?
(qc) #NAME?
Mc = 0 5

38.2068966
 0

-1.10344828
de rotar. Por otro lado el
signos se da para cumplir el
plo, los asentamientos
cial, queda:

Mab
Mba
Mbc
(qb)
(qc)
 PRIMERA DISTRIBUCION:

2 T/m 1 2 T/m 1 2 T/m

CV= 2.00 T/ml. CV= 2.00 T/ml. CV= 2.00 T/ml.
W1= 1.00 T/ml. W1= 1.00 T/ml.

I I I I I
5.0 m 5.0 m 7.50 m 5.0 m 5.0 m
1 2 3 4 5 6

5 5 7.5 5 5
FORMULA GENERAL

NUDO 2 W1= 2.00 W2= 1.00

M1= 0 M2= M2 M3= M3 L1 = 5.00 L2 = 5.00
M12 6.25 M23 3.125
Ar = 20.8333 Ar = 10.417
FACTORES DE CARGA
0 + 20 M2 + 5 M3 = -93.75 -62.50 -31.25

NUDO 3 W1= 1.00 W2= 2.00

M1= M2 M2= M3 M3= M4 L1 = 5.00 L2 = 7.50
M12 3.125 M23 14.063
Ar = 10.4167 Ar = 70.313
FACTORES DE CARGA
5 M2 + 25 M3 + 7.5 M4 = -242.19 -31.25 -210.938

NUDO 4 W1= 2.00 W2= 1.00

M1= M3 M2= M4 M3= M5 L1 = 7.50 L2 = 5.00
M12 14.0625 M23 3.125
Ar = 70.3125 Ar = 10.417
FACTORES DE CARGA
7.5 M3 + 25 M4 + 5 M5 = -242.19 -210.94 -31.25

NUDO 5 W1= 1.00 W2= 2.00

M1= M4 M2= M5 M3= 0 L1 = 5.00 L2 = 5.00
M12 3.125 M23 6.25
Ar = 10.4167 Ar = 20.833
FACTORES DE CARGA
5 M4 + 20 M5 + 5 (0) = -93.75 -31.25 -62.5

POR LO TANTO SE TIENE

M2 M3 M4 M5
 20 5 0 0 = -93.75 M2 = -2.938
5 25 7.5 0 = -242.19 M3 = -7.000
0 7.5 25 5 = -242.19 M4 = -7.000
0 0 5 20 = -93.75 M5 = -2.938

RESUMEN
MATRIZ INVERSA
0.05 -0.01 0.00 0.00 M1 = 0.000 tn-m R1 = 4.412
-0.01 0.05 -0.01 0.00 M2 = -2.938 tn-m R2 = 7.280
0.00 -0.01 0.05 -0.01 M3 = -7.000 tn-m R3 = 10.812
0.00 0.00 -0.01 0.05 M4 = -7.000 tn-m R4 = 10.812
M5 = -2.938 tn-m R5 = 7.275
M6 = 0.000 tn-m R6 = 4.412
A) Para el portico mostrado determinar el Diagrama de fuerza cortante, Momento Flector y desplaz
por ensamblaje.

2.00 T/m
f'c =
5T EA =
2 b2 3
Viga :
3.50 m Columna :
b1 b3

1 4

5.00 m

SOLUCIÓN

1.- Definimos el sistema Q-D y q-d

D1 D2
D3
b2

Q-D b1 q-d

3 GDL

2.- Generamos la matriz de transformacion de desplazamientos [A]

D1=1 Di =0 , i ≠1

0
1
[C1] = 1
0
0
0

D2=1 Di =0 , i ≠2
 0
0
[C2] = 0
1
1
0

D3=1 Di =0 , i ≠3

- 2/7
- 2/7
[C3] = 0
0
- 2/7
- 2/7

D1=1 D2=1 D3=1

0 0 - 2/7
[A1]
1 0 - 2/7
[A] = 1 0 0
[A2]
0 1 0
0 1 - 2/7
[A3]
0 0 - 2/7

3.- Calculamos el [k'] matriz de rigidez de barra en el sistema local

E = 2.17E+06 Tn/m2 2.173706512 [𝑘

Iviga = 0.003125 m4 ′]=[■8(4𝐸𝐼⁄𝐿&2𝐸𝐼⁄𝐿@2𝐸
𝐼⁄𝐿&4𝐸𝐼⁄𝐿)]
Icol = 0.000675 m4

barra 1:
1676.86 838.43
[k'1] =
838.43 1676.86

barra 2:
5434.27 2717.13
[k'2] =
2717.13 5434.27
 barra 3:
1676.86 838.43
[k'3] =
838.43 1676.86

4. Transformamos el [k'] al sistema global

[𝑘_𝑖 ]=[𝐴_𝑖 ]^𝑇∗[ 〖𝑘′〗 _𝑖 ]∗[𝐴_𝑖 ]

barra 1:
0 1 838.43
[k'1]*[A1] =
[A1]T = 0 0 1676.86
- 2/7 - 2/7

1676.86 0 -718.65
[k1] = 0 0 0
-718.65 0 410.66

barra 2:
1 0 5434.27
[k'2]*[A2] =
[A2]T = 0 1 2717.13
0 0

5434.27 2717.13 0
[k2] = 2717.13 5434.27 0
0 0 0

barra 3:
0 0 0
[k'3]*[A3] =
[A3]T = 1 0 0
- 2/7 - 2/7

0 0 0
[k3] = 0 1676.86 -718.65
0 -718.65 410.66

5.- Ensamblamos la matriz de rigidez [K]

[𝐾]=∑▒[𝑘_𝑖 ]
 [𝐾]=∑▒[𝑘_𝑖 ]

7111.13 2717.13 -718.65

[K] = [k1]+[k2]+[k3] = 2717.13 7111.13 -718.65
-718.65 -718.65 821.32

6.- Calculo de {R}

2.00 T/m
5T R1 R2
R3

2.00 T/m
5T 4.17 T-m 4.17 T-m 5T

R1 R2 R3

2 b2
5T 5T 5T 4.17 T-m 3 5T
5T
5T
5T

b3
b1

5T
5T

{q}p = problema primario

barra 1 barra 2
{r1}= 0 {r2}= -4.17
0 4.17

7.- Calculo de {D} de los desplazamientos en el sistema global

{𝑄}=[𝐾]∗{𝐷} {𝑄}=−{𝑅}
{𝐷}=[𝐾]^(−1)∗{𝑄}
 4.17 T-m 7111.13 2717.13
{Q} = -4.17 T-m [K] = 2717.13 7111.13
5.00 T -718.65 -718.65

1.72E-04 -5.55E-05 1.02E-04

[𝐾]^(−1)= -5.55E-05 1.72E-04 1.02E-04
1.02E-04 1.02E-04 1.40E-03

D1
{D} = D2
D3

Por lo tanto:
0.00146 rad
{D} = -0.00044 rad
0.00698 m

8.- Determinanos los esfuerzos y desplazamientos en el sistema local

Desplazaminetos

{𝑑_𝑖 }=[𝐴_𝑖 ]∗{𝐷}

barra 1
0 0 - 2/7 0.00146 -0.00199
{d1} = =
1 0 - 2/7 -0.00044 -0.00054
0.00698

barra 2
1 0 0 0.00146 0.00146
{d2} = =
0 1 0 -0.00044 -0.00044
0.00698

barra 3
0 1 - 2/7 0.00146 -0.00243
{d3} = =
0 0 - 2/7 -0.00044 -0.00199
0.00698
Esfuerzos (Momentos)

{𝑞_𝑖 }={𝑞_𝑖 }𝑝+{𝑞_𝑖 }𝑐

donde:
{𝑞_𝑖 }𝑐=[ 〖𝑘′〗 _𝑖 ]∗{𝑑_𝑖 }

Entonces:

{𝑞_𝑖 }={𝑞_𝑖 }𝑝+[ 〖𝑘′〗 _𝑖 ]∗{𝑑_𝑖 }

barra 1

0 1676.86 838.43 -0.00199 -3.79

{q1}= + =
0 838.43 1676.86 -0.00054 -2.57

barra 2

-4.17 5434.27 2717.13 0.00146 2.57

{q2}= + =
4.17 2717.13 5434.27 -0.00044 5.75

barra 3

0.00 1676.86 838.43 -0.00243 -5.75

{q3}= + =
0.00 838.43 1676.86 -0.00199 -5.38

Cortantes

{V}i = {V}isos - {(Mi + Mj)/Lij} i

{V}j = {V}isos+ {(Mi + Mj)/Lij}

0.00 1.82 T
0.00 -1.82 T
{V}isos = 5.00 {V}T = 3.34 T
5.00 6.66 T
0.00 3.18 T
0.00 -3.18 T

3.34

+
1.82 3.18
-
 1.67 m -

-6.66
+ +

1.82 3.18

V
Tn
ento Flector y desplazamiento lateral, mediante el metodo matricial

210 kg/cm2
∞
b (m) h (m)
0.30 0.50
0.30 0.30

q-d b3
 0 -718.65
0 -718.65

2717.13 0
5434.27 0

1676.86 -718.65
838.43 -718.65
 Hacemos equilibrio en los nudos para determinar el vector {R}
Nudo 2:
ΣM=0 ==> R1 = -4.17 T-m

Nudo 3:
ΣM=0 ==> R2 = 4.17 T-m

Σ Fx = 0 ==> R3 = -5 T

Por lo tanto
-4.17 T-m
{R} = 4.17 T-m
-5.00 T

barra 3
{r3}= 0
0
 1 2 3
-718.65 1
-718.65 2
821.32 3

rad
rad

rad
rad

rad
rad
T-m
T-m

T-m
T-m

T-m
T-m

-5.75

2.57 -
-5.75
+ -
 2.57 + -
5.35

- +

-3.79 5.38
M
Tn-m