Coordenadas y vector de conectividad

Elemento xf xi b h E I ciy ciz cfy cfz

1 12 0 0.20 0.25 200 ⋅ 10 6 562.5 ⋅ 10 -6 6 5 4 2

6 -6
2 20 12 0.20 0.25 200 ⋅ 10 562.5 ⋅ 10 4 2 3 1

1. DEFINICION DE MATRICES Y VARIABLES

Matriz de rigidez en coordenas locales Longitud elemental Vector de reacciones en coordenadas locales

⎡ 12 ⋅ E ⋅ I 6⋅E⋅I 12 ⋅ E ⋅ I 6 ⋅ E ⋅ I ⎤ ⎡ Vi ⎤
⎢ ――― ――― -――― ―――⎥ L ⎛⎝xf , xi⎞⎠ ≔ ⎛⎝xf - xi⎞⎠ ⎢ ⎥
L 3
L 2
L3 L2 M
⎢ ⎥ S ⎛⎝Vi , Mi , Vj , Mj⎞⎠ ≔ ⎢ i ⎥
⎢ 6⋅E⋅I 4⋅E⋅I 6⋅E⋅I 2⋅E⋅I ⎥ ⎢ Vj ⎥
⎢ ――― ――― -――― ―――⎥ Vector de conexiones ⎢M ⎥
L2 L L2 L ⎣ j⎦
k' ((E , I , L)) ≔ ⎢ ⎥
⎢ 12 ⋅ E ⋅ I 6 ⋅ E ⋅ I 12 ⋅ E ⋅ I 6⋅E⋅I ⎥ T
⎢ -――― -――― ――― -――― ⎥ c ⎛⎝ciy , ciz , cfy , cfz⎞⎠ ≔ ⎡⎣ ciy ciz cfy cfz ⎤⎦
L3 L2 L3 L2
⎢ ⎥
6⋅E⋅I
⎢ ――― 2⋅E⋅I 6⋅E⋅I 4⋅E⋅I
――― -――― ―――⎥
⎢⎣ L2 L L2 L ⎥⎦
 Momentos y cortantes para una carga distribuida Momentos y cortantes para una carga puntual
w ⋅ L2 w⋅L P ⋅ b2 ⋅ a P ⋅ b2 ⎛ b⎞
( )
Mwi (w , L) ≔ -―― Vwi ((w , L)) ≔ -―― ( )
Mpi (P , a , b , L) ≔ -――― ( )
Vpi (P , b , L) ≔ -―― ⋅ ⎜3 - 2 ⋅ ―⎟
12 2 L2 L2 ⎝ L⎠

w ⋅ L2 w⋅L P ⋅ a2 ⋅ b P ⋅ a2 ⎛ a⎞
Mwj ((w , L)) ≔ ―― Vwj ((w , L)) ≔ -―― Mpj ((P , a , b , L)) ≔ ――― Vpj ((P , a , L)) ≔ -―― ⋅ ⎜3 - 2 ⋅ ―⎟
12 2 L2 L2 ⎝ L⎠

Número de nodos Grados de libertad por Grados de libertad total NGL Libres NGL Restringidos Matriz Global
nudo
num ≔ 3 n≔2 N ≔ n ⋅ num = 6 s≔2 t≔4 K ≔0
N,N

i ≔ 1,2‥2 ⋅ n j ≔ 1,2‥2 ⋅ n

2. PROBLEMA PRIMARIO

Elemento 1

Cosenos directores y longitud. Carga Reacciones

L1 ≔ L ⎛xf , xi ⎞ = 12 w1 ≔ -6 M1i ≔ Mwi ⎛⎝w1 , L1⎞⎠ = 72 M1j ≔ Mwj ⎛⎝w1 , L1⎞⎠ = -72
⎜⎝ 1 1⎟⎠
V1i ≔ Vwi ⎛⎝w1 , L1⎞⎠ = 36 V1j ≔ Vwj ⎛⎝w1 , L1⎞⎠ = 36
⎡ 36 ⎤
⎢ 72 ⎥
s1 ≔ S ⎛⎝V1i , M1i , V1j , M1j⎞⎠ = ⎢ ⎥
⎢ 36 ⎥
⎣ -72 ⎦
Elemento 2

Cosenos directores y longitud.

L2 ≔ L ⎛xf , xi ⎞ = 8
⎜⎝ 2 2⎟⎠
3. CALCULO DE FUERZAS EQUIVALENTES

Fuerzas de fijacion (R). Vector de cargas (Q). Vector de reacciones por fuerzas
elementales (J).
Nudo 1.
R1 ≔ 20
⎡ 0 ⎤
⎢s +0⎥ ⎡ 0⎤
Nudo 2. ⎡ R ⎤ ⎡ 20 ⎤ ⎢ 13 ⎥ ⎢ ⎥
36
R2 ≔ -s1 = 72 Qc ≔ ⎢ 1 ⎥ = ⎢ ⎥ J ≔ s1 ⎥ = ⎢ ⎥
⎢
4
⎣ R2 ⎦ ⎣ 72 ⎦ ⎢ 2 ⎥ ⎢ 72 ⎥
⎢ ⎥ ⎣ 36 ⎦
⎢ s11 ⎥
⎣ ⎦

⎡0⎤
⎢0⎥
Q ≔0 q ≔0 qc ≔ ⎢ ⎥
⎢0⎥
N,1 N,1

⎣0⎦

3. CALCULO DE LA RIGIDEZ ELEMENTAL Y GLOBAL (PROBLEMA SECUNDARIO).

Elemento 1

Cosenos directores y longitud Matriz de rigidez elemental en coordenadas Matriz de rigidez elemental en Vector de conectividad.
globales. coordenada globales

L1 = 12 k1 ≔ k' ⎛E , I , L1⎞ K1 ≔0 c1 ≔ c ⎛ciy , ciz , cfy , cfz ⎞

⎝ 1 1 ⎠ N,N ⎜⎝ 1 1 1 1⎟
⎠
T
c1 = [[ 6 5 4 2 ]]
⎡ 781.25 4687.5 -781.25 4687.5 ⎤
⎢ 4687.5 37500 -4687.5 18750 ⎥ K1 ≔ k1
c1 , c1 i,j
k1 = ⎢ ⎥ i j

⎢ -781.25 -4687.5 781.25 -4687.5 ⎥

⎣ 4687.5 18750 -4687.5 37500 ⎦
 ⎡0 0 0 0 0 0 ⎤
⎢ 0 37500 0 -4687.5 18750 4687.5 ⎥
⎢ ⎥
0 0 0 0 0 0
K1 = ⎢ ⎥
⎢ 0 -4687.5 0 781.25 -4687.5 -781.25 ⎥
⎢ 0 18750 0 -4687.5 37500 4687.5 ⎥
⎢⎣ 0 4687.5 0 -781.25 4687.5 781.25 ⎥⎦
Elemento 2

Cosenos directores y longitud. Matriz de rigidez elemental en Matriz de rigidez elemental en Vector de conectividad.
coordenadas globales. coordenada globales

L2 = 8 k2 ≔ k' ⎛E , I , L2⎞ K2 ≔0 c2 ≔ c ⎛ciy , ciz , cfy , cfz ⎞

⎝ 2 2 ⎠ N,N ⎜⎝ 2 2 2 2⎟
⎠
T
c2 = [[ 4 2 3 1 ]]

K2 ≔ k2
c2 , c2 i,j
i j

⎡ 2636.72 10546.88 -2636.72 10546.88 ⎤

⎢ 10546.88 56250 -10546.88 28125 ⎥ ⎡ 56250 28125 -10546.88 10546.88 0 0 ⎤
k2 = ⎢ ⎥ ⎢ 28125
⎢ -2636.72 -10546.88 2636.72 -10546.88 ⎥ 56250 -10546.88 10546.88 0 0 ⎥
⎢ ⎥
⎣ 10546.88 28125 -10546.88 56250 ⎦ -10546.88 -10546.88 2636.72 -2636.72 0 0 ⎥
K2 = ⎢
⎢ 10546.88 10546.88 -2636.72 2636.72 0 0 ⎥
⎢ 0 0 0 0 0 0⎥
⎢⎣ 0 0 0 0 0 0 ⎥⎦

Matriz de rigidez global

K ≔ K1 + K2

⎡ 56250.00 28125.00 -10546.88 10546.88 0.00 0.00 ⎤

⎢ 28125.00 93750.00 -10546.88 5859.38 18750.00 4687.50 ⎥
⎢ ⎥
-10546.88 -10546.88 2636.72 -2636.72 0.00 0.00 ⎥
K=⎢
⎢ 10546.88 5859.38 -2636.72 3417.97 -4687.50 -781.25 ⎥
⎢ 0.00 18750.00 0.00 -4687.50 37500.00 4687.50 ⎥
⎢⎣ 0.00 4687.50 0.00 -781.25 4687.50 781.25 ⎥⎦
 ⎢ 10546.88 5859.38 -2636.72 3417.97 -4687.50 -781.25 ⎥
⎢ 0.00 18750.00 0.00 -4687.50 37500.00 4687.50 ⎥
⎢⎣ 0.00 4687.50 0.00 -781.25 4687.50 781.25 ⎥⎦

Sub Matrices

K11 ≔ submatrix ((K , 1 , s , 1 , s)) K12 ≔ submatrix ((K , 1 , s , s + 1 , s + t))

⎡ 56250 28125 ⎤ ⎡ -10546.88 10546.88 0 0 ⎤

K11 = ⎢ K12 = ⎢
⎣ 28125 93750 ⎥⎦ ⎣ -10546.88 5859.38 18750 4687.5 ⎥⎦

K21 ≔ submatrix ((K , s + 1 , s + t , 1 , s)) K22 ≔ submatrix ((K , s + 1 , s + t , s + 1 , s + t))

⎡ -10546.88 -10546.88 ⎤ ⎡ 2636.72 -2636.72 0 0 ⎤
⎢ 10546.88 5859.38 ⎥ ⎢ -2636.72 3417.97 -4687.5 -781.25 ⎥
K21 = ⎢ ⎥ K22 = ⎢ ⎥
⎢ 0 18750 ⎥ ⎢ 0 -4687.5 37500 4687.5 ⎥
⎣ 0 4687.5 ⎦ ⎣ 0 -781.25 4687.5 781.25 ⎦

4. CÁLCULO DE DESPLAZAMIENTOS Y REACCIONES (COORD, GLOBALES)

Desplazamientos en los grados de libertad libres Reacciones

⎡0⎤ ⎡ -7.85 ⎤
⎡ 20 ⎤ ⎢0⎥ ⎡0 ⎤ ⎢ 40.21 ⎥
Qc = ⎢ ⎥ qc = ⎢ ⎥ qnc ≔ K11 -1 ⋅ ⎛⎝Qc - K12 ⋅ qc⎞⎠ = ⎢ Qnc ≔ K21 ⋅ qnc + K22 ⋅ qc + J = ⎢ ⎥
⎣ 72 ⎦ ⎣ 0.0008 ⎥⎦
⎢0⎥ ⎢ 86.59 ⎥
⎣0⎦ ⎣ 39.65 ⎦

⎡0⎤ ⎡ 20 ⎤
⎢0⎥ ⎢ 72 ⎥
⎢ ⎥ ⎢ ⎥
0 -7.85 ⎥
q ≔ stack ⎛⎝qnc , qc⎞⎠ = ⎢ ⎥ Q ≔ stack ⎛⎝Qc , Qnc⎞⎠ = ⎢
⎢0⎥ ⎢ 40.21 ⎥
⎢0⎥ ⎢ 86.59 ⎥
⎢⎣ 0 ⎥⎦ ⎢⎣ 39.65 ⎥⎦
5. CÁLCULO DE FUERZAS ELEMENTALES (SUMA PROBLEMA PRIMARIO MÁS SECUNDARIO)

Elemento 1 Elemento 2
⎡q ⎤ ⎡q ⎤
⎢ 6⎥ ⎡ ⎤ ⎢ 4⎥ ⎡ ⎤
⎢ q ⎥ ⎢ 3.65 ⎥ ⎢ q ⎥ ⎢ 7.85 ⎥
5 14.59 2 42.82
F1 ≔ k' ⎛E , I , L1⎞ ⋅ ⎢ ⎥ = ⎢ ⎥ F2 ≔ k' ⎛E , I , L2⎞ ⋅ ⎢ ⎥ = ⎢ ⎥
⎝ 1 1 ⎠ q4⎢ ⎥ ⎢ -3.65 ⎥ ⎝ 2 2 ⎠ q3⎢ ⎥ ⎢ -7.85 ⎥
⎢ ⎥ ⎣ 29.18 ⎦ ⎢ ⎥ ⎣ 20 ⎦
⎢q ⎥ ⎢q ⎥
⎣ 2⎦ ⎣ 1⎦