ANÁLISIS MATRICIAL DE ESTRUCTURAS -
ARMADURAS 2D
Entrada de Datos:
Nodos (Coordenadas)
N X Y
1 0 0
2 6 8.4
3 18 0
4 12 0
5 9 3
6 6 0
Unidad de Longitud: UL ≔ m
Miembros (Continuidad y Propiedades)
M Ni Nj A E
⎛ tonnef ⎞
⎛⎝cm 2 ⎞⎠ ⎜――― ⎟
2
⎝ cm ⎠
1 1 6 1 1
2 6 5 1 1
3 5 4 1 1
4 5 3 1 1
5 3 2 1 1
6 4 2 1 1
7 2 6 1 1
8 2 1 1 1
Apoyos (Restricciones)
Nodo Ux Uy
Nota:
1 0 1
1 = Si restringido
3 0 1 0 = No restringido
5 1 1 Es necesario al menos definir
dos apoyos
Fuerzas
Fuerzas
NLF Fx Fy
((tonnef )) ((tonnef ))
2 -15 0
4 0 0
SOLUCIÓN
De la figura 2.13 tenemos que, conforme al sistema global de referencia:
1) Ensamblaje del Vector Fuerza { F }
⎡ 0⎤
⎢ 0⎥
⎢ ⎥ ⎡ 0⎤
⎢ -15 ⎥ ⎢ -15 ⎥
⎢ 0⎥ ⎢ ⎥
⎢ 0⎥ ⎢ 0⎥
⎢ 0⎥ 0⎥
F=⎢ ⎥ tonnef FR = ⎢ tonnef
0 ⎢ 0⎥
⎢ ⎥ ⎢ 0⎥
⎢ 0⎥ ⎢ ⎥
⎢ 0⎥ ⎢
0
⎥
⎢ 0⎥ ⎣ 0⎦
⎢ 0⎥
⎢ ⎥
⎣ 0⎦
μ'' T = [[ 0 1 0 0 0 1 0 0 1 1 0 0 ]]
2) Definición de las longitudes y cosenos directores de una barra
⎡ 6000 ⎤ ⎡ 1 ⎤ ⎡ 0 ⎤
⎢ 4242.641 ⎥ ⎢ 0.707 ⎥ ⎢ 0.707 ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢ 4242.641 ⎥ ⎢ 0.707 ⎥ ⎢ -0.707 ⎥
9486.833 ⎥ 0.949 ⎥ -0.316 ⎥
L=⎢ mm l = ⎢ m=⎢
⎢ 14647.867 ⎥ ⎢ -0.819 ⎥ ⎢ 0.573 ⎥
⎢ 10322.79 ⎥ ⎢ -0.581 ⎥ ⎢ 0.814 ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
8400 0 -1
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ 10322.79 ⎦ ⎣ -0.581 ⎦ ⎣ -0.814 ⎦
⎢ 10322.79 ⎥ ⎢ -0.581 ⎥ ⎢ 0.814 ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
8400 0 -1
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ 10322.79 ⎦ ⎣ -0.581 ⎦ ⎣ -0.814 ⎦
PASO 1: Definir las matrices de continuidad [a] y de rigideces de las barras [ k ]
⎡ ⎡ 1 ⎤ ⎡ 0.707 ⎤ ⎡ 0.707 ⎤ ⎡ 0.949 ⎤ ⎡ -0.819 ⎤ ⎡ -0.581 ⎤ ⎡ 0 ⎤ ⎡ -0.581 ⎤ ⎤
μ' T = ⎢ ⎢ ⎥ ⎢
⎣ ⎣ 0 ⎦ ⎣ 0.707 ⎥⎦ ⎢⎣ -0.707 ⎥⎦ ⎢⎣ -0.316 ⎥⎦ ⎢⎣ 0.573 ⎥⎦ ⎢⎣ 0.814 ⎥⎦ ⎢⎣ -1 ⎥⎦ ⎢⎣ -0.814 ⎥⎦ ⎥⎦
⎡ -1 0 0 0 0 0 0 0 0 0 1 0 ⎤
⎢ 0 0 0 0 0 0 0 0 0.707 0.707 -0.707 -0.707 ⎥
⎢ ⎥
⎢ 0 0 0 0 0 0 0.707 -0.707 -0.707 0.707 0 0 ⎥
⎢ 0 0 0 0 0.949 -0.316 0 0 -0.949 0.316 0 0 ⎥
a=
⎢ 0 0 -0.819 0.573 0.819 -0.573 0 0 0 0 0 0 ⎥
⎢ 0 0 -0.581 0.814 0 0 0.581 -0.814 0 0 0 0 ⎥
⎢ ⎥
0 0 0 1 0 0 0 0 0 0 0 -1
⎢ ⎥
⎣ -0.581 -0.814 0.581 0.814 0 0 0 0 0 0 0 0 ⎦
⎡ -1 0 0 0 0 0 1 0 ⎤
⎢ 0 0 0 0 0 0 -0.707 -0.707 ⎥
⎢ ⎥
⎢ 0 0 0 0 0.707 -0.707 0 0 ⎥
0 0 0 0.949 0 0 0 0
aR = ⎢ ⎥
⎢ 0 -0.819 0.573 0.819 0 0 0 0 ⎥
⎢ 0 -0.581 0.814 0 0.581 -0.814 0 0 ⎥
⎢ ⎥
0 0 1 0 0 0 0 -1
⎢ ⎥
⎣ -0.581 0.581 0.814 0 0 0 0 0 ⎦
⎡ 0 0 0 0 ⎤
⎢ 0 0 0.707 0.707 ⎥
⎢ ⎥
⎢ 0 0 -0.707 0.707 ⎥
0 -0.316 -0.949 0.316 ⎥
aH = ⎢
⎢ 0 -0.573 0 0 ⎥
⎢ 0 0 0 0 ⎥
⎢ ⎥
0 0 0 0 ⎡ 6000 ⎤
⎢ ⎥
⎣ -0.814 0 0 0 ⎦ ⎢ 4242.641 ⎥
⎢ ⎥
⎢ 4242.641 ⎥
9486.833 ⎥
L=⎢ mm
⎢ 14647.867 ⎥
⎢ 10322.79 ⎥
⎢ ⎥
8400
⎢ ⎥ μ'' T = [[ 0 1 0 0 0 1 0 0 1 1 0 0 ]]
⎣ 10322.79 ⎦
⎡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 ⎥ tonnef
k=⎢ ―――
⎢0 0 0 0 0 0 0 0 ⎥ mm
⎢0 0 0 0 0 0 0 0⎥
⎢ ⎥
0 0 0 0 0 0 0 0
⎢ ⎥
⎣0 0 0 0 0 0 0 0⎦
PASO 2: Obtener la matriz global de rigidez [ K ]
[ K ]=[ a ] T [ k ][ a ]
⎡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 ⎥ tonnef
K=⎢ ―――
⎢0 0 0 0 0 0 0 0 ⎥ mm
⎢0 0 0 0 0 0 0 0⎥
⎢ ⎥
0 0 0 0 0 0 0 0
⎢ ⎥
⎣0 0 0 0 0 0 0 0⎦
PASO 3: Obtener el vector de los desplazamientos de los nudos { μ }
⎡ -140.372 ⎤ { μ }=[ K ] -1 { F }
⎢ -331.866 ⎥
⎢ ⎥
⎢ 62.857 ⎥
-122.411 ⎥
μ=⎢ m
⎢ 1049.665 ⎥
⎢ 1049.665 ⎥
⎢ ⎥
-120.049
⎢ ⎥
⎣ 91.309 ⎦
PASO 4: Obtener el vector de las deformaciones de las barras { e }
{ e }=[ a ]{ μ }
⎡ 20.323 ⎤
⎢ 20.323 ⎥
⎢ ⎥
⎢ 0 ⎥
⎢ -116.129 ⎥
e= m
⎢ 207.639 ⎥
⎢ 0 ⎥
⎢ ⎥
-28.452
⎢ ⎥
⎣ -60.155 ⎦
PASO 5: Obtener el vector de las fuerzas en las barras { s }
{ s }=[ k ]{ e }
⎡ 3.387 ⎤
⎢ 4.79 ⎥
⎢ ⎥
⎢ 0 ⎥
⎢ -12.241 ⎥
s= tonnef
⎢ 14.175 ⎥
⎢ 0 ⎥
⎢ ⎥
-3.387
⎢ ⎥
⎣ -5.827 ⎦
PASO 6: Comprobar el vector fuerzas externas { F }
{ F }=[ a ] T { s }
⎡ 0⎤
⎢ -15 ⎥
⎢ ⎥
⎢ 0⎥
0⎥
F' = ⎢ tonnef
⎢ 0⎥
⎢ 0⎥
⎢ ⎥
0
⎢ ⎥
⎣ 0⎦
PASO 7: Obtener el vector de reacciones { H }
{ H }=[ aH ] T { s }
⎡ -1 0 0 0 0 0 0 0 0 0 1 0 ⎤
⎢ 0 0 0 0 0 0 0 0 0.707 0.707 -0.707 -0.707 ⎥
⎢ ⎥
⎢ 0 0 0 0 0 0 0.707 -0.707 -0.707 0.707 0 0 ⎥
⎢ 0 0 0 0 0.949 -0.316 0 0 -0.949 0.316 0 0 ⎥
a=
⎢ 0 0 -0.819 0.573 0.819 -0.573 0 0 0 0 0 0 ⎥
⎢ 0 0 -0.581 0.814 0 0 0.581 -0.814 0 0 0 0 ⎥
⎢ ⎥
0 0 0 1 0 0 0 0 0 0 0 -1
⎢ ⎥
⎣ -0.581 -0.814 0.581 0.814 0 0 0 0 0 0 0 0 ⎦
⎡ 0 0 0 0 ⎤
⎢ 0 0 0.707 0.707 ⎥
⎢ ⎥
⎢ 0 0 -0.707 0.707 ⎥
0 -0.316 -0.949 0.316 ⎥
aH = ⎢
⎢ 0 -0.573 0 0 ⎥
⎢ 0 0 0 0 ⎥
⎢ ⎥
0 0 0 0
⎢ ⎥
⎣ -0.814 0 0 0 ⎦
⎡ 4.742 ⎤
⎢ -4.258 ⎥ μ'' T = [[ 0 1 0 0 0 1 0 0 1 1 0 0 ]]
H=⎢ ⎥ tonnef
⎢ 15 ⎥
⎣ -0.484 ⎦
ΣFx