Taller de Analisis Matricial
Taller de Analisis Matricial
Taller de Analisis Matricial
PROGRAMA
ESPECIALIZACION EN ESTRUCTURAS
DOCENTE
INGENIERO DONALDO BARRETO NUÑEZ
ALUMNO
INGENIERO VICTOR ALFONSO ROMERO DIAZ
FECHA
21/05/201
CARTAGENA - COLOMBIA
EJERCICIOS
En la siguiente viga, escoja un sistema de coordenadas, determine el vector Q de las dos formas posibles, halle las
deformaciones, las reacciones y dibuje los diagramas.
SOLUCION
Características de los materiales y encontramos el Vector Q haciendo todos los nudos empotrados cuando Qn=0 y
encontramos cada reacción
27523 7462
7462 29849
K
Hallando la inversa de la matriz por medio de Excel y multiplicando por el vector Q se obtiene los desplazamientos
de las estructura, en los puntos B y C.
-0,002517431
0,003141964
q= radianes
Multiplicando los desplazamientos reales por los valores de cada reacción y por los momentos generalizados se,
hallan los valores reales para encontrar las reacciones y graficar los diagramas de cortante y momento.
Ɵ1 -0,00251743
Ɵ2 0,00314196
SOLUCION
Características de los materiales y encontramos el Vector Q haciendo todos los nudos empotrados cuando Qn=0 y
encontramos cada reacción y momentos por medio de tablas
SECCION VIGA
B(m) 0,30
H(m) 0,36
AREA(m²) 0,108
SECCION COLUMNA
B(m) 0,24
H(m) 0,36
AREA(m²) 0,0864
-33,33
33,75
VECTOR Q=
-67,78
-45
Vector Q por Trabajo Virtual
Matriz de Rigidez
Se halla la matriz de rigidez , haciendo cada Qn=1 y los demás Qni=0 como ilustra la siguiente imagen se ensambla
cada columna para armar la matriz K
Reemplazando por las propiedad de los materiales como inercia, modulo y longitudes.
-0,001461788 Radianes
0,0001058551 mts
q -0,0001825923 mts
-0,001453985 Radianes
b) Condición que todo los elementos son axialmente rígidos An=0
Comprobación del Vector Q= T V y Matriz de Rigidez
Reemplazando las propiedades de los materiales en la matriz de rigidez obtenemos
K₁₁ K₂₁
K
K₁₂ K₂₂
4EI/Lv 2EI/Lv
K
2EI/Lv 4EI/Lv + 4EI/Lc
15552 5832
K
5832 23846
- 0,001580422
Radianes
q= - 0,001500561
c) Condición elemento Horizontal transversalmente rígido Ef=0
Comprobación del Vector Q por T.V y Matriz de Rigidez
Reemplazando las propiedades de los materiales en la matriz de rigidez obtenemos
Ɵ 0,16666667
V1 Para el Q1 -5184
V1 Para el Q2 2304
K₁₁ K₂₁
K
K₁₂ K₂₂
AE/Lv+12EI/LC^3 6EI/LC^2*Ɵ
K
V1 AE/Lc-V1
361555 518
K
5184 429696
0,000094
Radianes
q= - 0,000189
En la siguiente estructura, determine un sistema de coordenadas, el vector Q de las dos formas, y sus
deformaciones para las siguientes consideraciones: A) Todos los elementos elásticos. B) Todos los elementos
axialmente rígidos ( n = 0), y C) El elemento horizontal izquierdo, transversalmente rígido ( f = 0)
a) Todos los elementos elásticos
Comprobación del Vector Q por T.V
Matriz de Rigidez
Propiedades de los materiales
Matriz de Rigidez reemplazando los valores de las propiedades de lo materiales encontramos los Kn para
armar la matriz de rigidez de cada miembro o elemento haciendo las deformaciones cada una igual a 1.
K₁₁ K₂₁ K₃₁ K₄₁
K₁₂ K₂₂ K₃₂ K₄₂
K
K₁₃ K₂₃ K₃₃ K₄₃
K₁₄ K₂₄ K₃₄ K₄₄
0,017496663 Mts
-0,143730769 Mts
q 0,067802495 Radianes
-0,017559163 Mts
0,141133591 Mts
-0,026986376 Radianes
b) Todos los elementos En=0
Comprobación del Vector Q por T.V y Matriz de Rigidez
Matriz de Rigidez reemplazando los valores de las propiedades de los materiales encontramos los Kn para
armar la matriz de rigidez de cada miembro o elemento haciendo las deformaciones cada una igual a 1.
SECCION VIGA
B(m) 0,25
H(m) 0,40
AREA(m²) 0,10
SECCION COLUMNA A Y B
B(m) 0,25
H(m) 0,3
AREA(m²) 0,075
K 4EI/Lv+4EI/Lc 8EI/Lv+4EI/Lc
2EI/Lv 2EI/Lv
31750 181083
K 10667 10667
SECCION COLUMNA A Y B
B(m) 0,25
H(m) 0,3
AREA(m²) 0,075
EA/Lv+12EI/Lc^3 0 0 0
0 2EA/Lc+12EI/Lvb^3 60EI/Lvb^3+6EI/Lvb^2+5EA/Lc 0
6EI/Lc^2 -5EA/Lc+60EI/Lvb^3+
4EI/Lvb+36EI/Lvb^2+4EI/Lc+300EI/Lvb^3+25EA/Lc
6EI/Lvb^2 6EI/Lc^2
K 0 0 0 12EI/Lc^2
403472 0 0 0
0 1002560 2519200 0
5208 2519200 12634150 5208
K 0 0 0 3472
0,0000024785 0 0 0 400
0,00000000515 0,00000199903 -0,00000039860 0,00000059790 -1100
K^-1 -0,00000000205 -0,00000039860 0,00000015863 -0,00000023794
x -2781,25
0,00000000000 0,00000000000 0,00000000000 0,00028800000 25
0,00099139 mts
-0,00107333 mts
q -0,00000950 Radianes
0,00720000 mts
En la siguiente estructura, determine un sistema de coordenadas, el vector Q de las dos formas, y sus deformaciones
para las siguientes consideraciones: A) Todos los elementos elásticos. B) Todos los elementos axialmente rígidos ( n
= 0), y C) El elemento inclinado izquierdo, transversalmente rígido ( f = 0) .
a) Todos los elementos elásticos
Comprobación del Vector Q por Q=T.V
En la siguiente Pórtico , halle las deformaciones.
A) Todos los elementos elásticos
Comprobación del Vector Q por TV
Matriz de Rigidez
Reemplazando los valores de las propiedades de los materiales
SECCION VIGA
B(m) 0,24
H(m) 0,60
AREA(m²) 0,144
SECCION COLUMNA
B(m) 0,3
H(m) 0,24
AREA(m²) 0,072
360.600 0 1.800 0
0 242.025 8.100 8.100
K
1.800 8.100 50.400 21.600
0 8.100 21.600 43.200
Reemplazando los valores delas propiedades de los materiales hallamos las deformaciones por medio de las operaciones con
matrices
50400 21600
21600 43200
0,010195707
q = - 0,005387205
Radianes
AE/Lv+12EI/LC^3 6EI/LC^2*Ɵ
K
-6EI/LC^2*Ɵ AE/Lc-4EI/L*Ɵ
360600 225
K
-225 239100
0,000899771
q = 0,002410926
Radianes
En la siguiente viga, , halle las deformaciones, y las reacciones por medio de la matriz de rigidez directamente.
Hallamos el vector Q por medio del equilibrio estático y por medio de hacer los miembros empotrados
Tabulamos las propiedades de los materiales y hallamos las matrices de cada uno de los miembros por separado para realizar la
superposición de matrices y aplicar las fórmulas para hallar desplazamientos y reacciones.
3 4 5 1
60,47018596 755,8773245 -60,47018596 755,8773245 3
755,8773245 12597,95541 -755,8773245 6298,977704 4
K360 -60,47018596 -755,8773245 60,47018596 -755,8773245 5
755,8773245 6298,977704 -755,8773245 12597,95541 1
5 1 6 2
49,74914769 746,2372154 -49,74914769 746,2372154 5
K400a 746,2372154 14924,74431 -746,2372154 7462,372154 1
-49,74914769 -746,2372154 49,74914769 -746,2372154 6
746,2372154 7462,372154 -746,2372154 14924,74431 2
6 2 7 8
49,74914769 746,2372154 -49,74914769 746,2372154 6
746,2372154 14924,74431 -746,2372154 7462,372154 2
K400b -49,74914769 -746,2372154 49,74914769 -746,2372154 7
746,2372154 7462,372154 -746,2372154 14924,74431 8
Matriz de Superposicion de Estrucutura Viga
1 2 3 4 5 6 7 8
27522,69972 7462,372154 755,8773245 6298,977704 -9,640109056 -746,2372154 0 0 1
7462,372154 29849,48862 0 0 746,2372154 0 -746,2372154 7462,372154 2
755,8773245 0 60,47018596 755,8773245 -60,47018596 0 0 0 3
6298,977704 0 755,8773245 12597,95541 -755,8773245 0 0 0 4
k -9,640109056 746,2372154 -60,47018596 -755,8773245 110,2193337 -49,74914769 0 0 5
-746,2372154 0 0 0 -49,74914769 99,49829539 -49,74914769 746,2372154 6
0 -746,2372154 0 0 0 -49,74914769 49,74914769 -746,2372154 7
0 7462,372154 0 0 0 746,2372154 -746,2372154 14924,74431 8
Knn Kna
-45,84 27522,69972 7462,372154 755,8773245 6298,977704 -9,640109056 -746,237215 0 0 ƟB
75 Qn 7462,372154 29849,48862 0 0 746,2372154 0 -746,237215 7462,37215 ƟC qn
VA-25 755,8773245 0 60,47018596 755,8773245 -60,47018596 0 0 0 0
MA-140,16
= 6298,977704 0 755,8773245 12597,95541 -755,8773245 0 0 0 x 0
VB-55 -9,640109056 746,2372154 -60,47018596 -755,8773245 110,2193337 -49,7491477 0 0 0 qa
VC-40 -746,2372154 0 0 0 -49,74914769 99,4982954 -49,7491477 746,237215 0
VD-10 0 -746,2372154 0 0 0 -49,7491477 49,7491477 -746,237215 0
MD+75 Qa 0 7462,372154 0 0 0 746,237215 -746,237215 14924,7443 0
Kan Kaa
Haciendo las operaciones entre matrices y aplicando las fórmulas hallamos los desplazamientos.
qn=Knn^-1*Qn
VA -1,902869135 -25
MA -15,85724279 -104,16
VB 2,368918517 -55
VC
=
1,878600824 - -40
VD -2,344650206 -10
MD 23,44650206 75
VA 23,09713087
MA 88,30275721
VB 57,36891852
=
VC 41,87860082
VD 7,655349794
MD -51,55349794
Punto B cuando el apoyo baja 0.3 pulgadas halle los desplazamientos y reacciones
qn=(Knn^-1)*Qn-(Knn^(-1) )*Kna*qa
ƟB -0,002517431 0,000727125
ƟC
=
0,003141964 - -0,000181781
ƟB -0,003244556
=
ƟC 0,003323745
Reacciones
Qa=Kan*qn+Kaa*qa
VA -2,452486632 -25 0
MA -20,4373886 -104,16 0
VB 2,511580066 -55 1,243728692
VC
=
2,421208754 - -40 + -2,487457385
VD -2,480302188 -10 1,243728692
MD 24,80302188 75 -18,65593039
VA 22,54751337
MA 83,7226114
VB 58,75530876
=
VC 39,93375137
VD 8,763426504
MD -68,8529085
Para la siguiente parilla halle sus desplazamientos y sus reacciones y diagramas
DESCRIPCION PROPIEDADES DE LOS MATERIALES
MODULO DE
VIGA DESCRIPCION ELEMENTO AREA LONGITUD BASE ALTURA INERCIA
ELASTICIDAD
A 1 2-1 20.000.000 0,09000 4,5 0,3 0,3 0,000675
B 2 3-1 20.000.000 0,09000 7,5 0,3 0,3 0,000675
G 8.000.000
B 0,140833333
J 0,00114075
Ɵ 143,14
4 5 6 1 2 3
12000,000 0,000 4000,000 6000,000 0,000 -4000,000 4
Viga 1 0,000 2028,000 0,000 0,000 -2028,000 0,000 5
4000,000 0,000 1777,778 4000,000 0,000 -1777,778 6
6000,000 0,000 4000,000 12000,000 0,000 -4000,000 1
0,000 -2028,000 0,000 0,000 2028,000 0,000 2
-4000,000 0,000 -1777,778 -4000,000 0,000 1777,778 3
7 8 9 1 2 3
3369,75982 2871,64643 863,80098 516,44924 2311,83088 -863,80098 7
2871,64643 5047,04018 1152,14924 2311,83088 1866,75076 -1152,14924 8
Viga 2 863,80098 1152,14924 384,00000 863,80098 1152,14924 -384,00000 9
516,44924 2311,83088 863,80098 3369,75982 2871,64643 -863,80098 1
2311,83088 1866,75076 1152,14924 2871,64643 5047,04018 -1152,14924 2
-863,80098 -1152,14924 -384,00000 -863,80098 -1152,14924 384,00000 3
Knn Kna
1 2 3 4 5 6 7 8 9
15369,76 2871,65 -4863,80 6000,00 0,00 4000,00 516,44924 2311,8309 863,8009831 1
2871,65 7075,04 -1152,149236 0 -2028,00 0 2311,8309 1866,7508 1152,149236 2
K Rigidez de -4863,80 -1152,149236 2162 -4000,00 0 -1778 -863,801 -1152,149 -384 3
Superposición 6000,00 0 -4000,00 12000,00 0 4000,00 0 0 0 4
0,00 -2028,00 0 0 2028,00 0,00 0 0 0 5
Parilla 4000,00 0 -1778 4000,00 0,00 1778 0 0 0 6
516,4492383 2311,830878 -863,8009831 0 0 0 3369,7598 2871,6464 863,8009831 7
2311,830878 1866,750762 -1152,149236 0 0 0 2871,6464 5047,0402 1152,149236 8
863,8009831 1152,149236 -384 0 0 0 863,80098 1152,1492 384 9
Kan Kaa
Qn Knn qn
28,15 15369,75982 2871,646427 -4863,800983 Ɵ1CX
Ɵ1CX -0,04496107
Ɵ2CY = -0,00344834
∆CZ -0,14990194
Qa Kan
MAX-33,75 6000 0 -4000
MAY 0 -2028 0 qn
RAZ-45 4000 0 -1777,777778 -0,04496107
MXB-25,6 516,4492383 2311,830878 -863,8009831 x -0,00344834
=
MYB+19,2 2311,830878 1866,750762 -1152,149236 -0,14990194
RBZ+6,4 863,8009831 1152,149236 -384
MAX 363,591348
MAY 6,993225775
RAZ 131,6480673
MXB
= 123,8933609
MYB 43,12983642
RBZ 8,35193266
Para la siguiente Parilla hallar sus desplazamientos y reacciones
De acuerdo a la gráfica se sacan las áreas por medio de AutoCAD y se hace el avaluó de cargas , para cargar cada una de las
vigas del problema.
Se toma como base para el cálculo de las reacciones los parámetros del Nrs 10 como módulo de elasticidad para un concreto
de 4000 psi combinaciones de carga de acuerdo al uso del diseño.
En las siguientes tablas están los parámetros para encontrar los valores.
CARGA VIVA
VIGA DESCRIPCION CONECTIVIDAD CARGA VIVA KG AREA AFERENTE M2 ESPESOR MTS PESO PROPIO KG CARGA MUERTA KG CARGA TOTAL CARGA W (KN/ML)
NORMA KG/M2
1 VIGA A-C 500 2040 2,55 0,15 1296 2656,8 4696,8 11,742
2 VIGA C-D 500 3400 4,25 0,15 1749,6 3935,52 7335,52 13,584
3 VIGA B-D 500 6448 8,06 0,15 1511,3952 5295,59424 11743,59424 25,175
1 2 3 4 5 6
3521 0 0 -3521 0 0 1
0 41968 -11658 0 20984 11658 2
0 -11658 4318 0 -11658 -4318 3
VIGA C-D
-3521 0 0 3521 0 0 4
0 20984 -11658 0 41968 11658 5
0 11658 -4318 0 11658 4318 6
10 11 12 4 5 6
36802 -19635 13396 16782 -12515 -13396 10
-19635 15858 -8038 -12515 3433 8038 11
13396 -8038 6698 13396 -8038 -6698 12
VIGA B-D
16782 -12515 13396 36802 -19635 -13396 4
-12515 3433 -8038 -19635 15858 8038 5
-13396 8038 -6698 -13396 8038 6698 6
Knn Kna
1 2 3 4 5 6 7 8 9 10 11 12
60178,59 0,00 -21246,42 -3521,48 0,00 0,00 28328,56 0,00 21246,42 0,00 0,00 0,00 1
0,00 46722,23 -11657,84 0,00 20984,12 11657,84 0,00 -4754,00 0,00 0,00 0,00 0,00 2
-21246,42 -11657,84 14940,93 0,00 -11657,84 -4317,72 -21246,42 0,00 -10623,21 0,00 0,00 0,00 3
-3521,48 0,00 0,00 40323,10 -19635,08 -13395,85 0,00 0,00 0,00 16782,21 -12515,22 13395,85 4
0,00 20984,12 -11657,84 -19635,08 57825,77 19695,35 0,00 0,00 0,00 -12515,22 3432,64 -8037,51 5
0,00 11657,84 -4317,72 -13395,85 19695,35 11015,59 0,00 0,00 0,00 -13395,85 8037,51 -6697,87 6
MATRIZ K
28328,56 0,00 -21246,42 0,00 0,00 0,00 56657,11 0,00 21246,42 0,00 0,00 0,00 7
0,00 -4754,00 0,00 0,00 0,00 0,00 0,00 4754,00 0,00 0,00 0,00 0,00 8
21246,42 0,00 -10623,21 0,00 0,00 0,00 21246,42 0,00 10623,21 0,00 0,00 0,00 9
0,00 0,00 0,00 16782,21 -12515,22 -13395,85 0,00 0,00 0,00 36801,62 -19635,08 13395,85 10
0,00 0,00 0,00 -12515,22 3432,64 8037,51 0,00 0,00 0,00 -19635,08 15857,54 -8037,51 11
0,00 0,00 0,00 13395,85 -8037,51 -6697,87 0,00 0,00 0,00 13395,85 -8037,51 6697,87 12
Kan Kaa
qn = Qn x Knn-1
-0,00979078
0,007949038
-0,025344113
qn =
-0,018849959
0,008133068
-0,069010188
MxA-15.66
MyA
RvA-23.48
Qa = Kan X qn
MxB-39
MyB +23,49
RvB-58.72