1
EJERCICIOS PASO A PASO DEL METODO MATRICIAL DE RIGIDEZ PARA EL
CÁLCULO DE ESTRUCTURALES ESQUELETALES
EDICSON ALEXANDER ALVAREZ SANCHEZ
POLITECNICO COLOMBIANO JAIME ISAZA CADAVID
FACULTAD DE INGENIERIA
INGENIERIA CIVIL
MEDELLIN
2009
2
EJERCICIOS PASO A PASO DEL METODO MATRICIAL DE RIGIDEZ PARA EL
CÁLCULO DE ESTRUCTURALES ESQUELETALES
EDICSON ALEXANDER ALVAREZ SANCHEZ
Trabajo de grado para optar el título de INGENIERO CIVIL.
Asesor Temático:
GIOVANNI MARTÍNEZ MARTÍNEZ
Ingeniero Civil
Especialista en Análisis y Diseño Estructural
Magister en Ingeniería Sismo resistente
Asesor(a) metodológica:
MARTHA ELENA ZAPATA PEREZ
Ing. Civil. Especialista en gestión pública.
POLITECNICO COLOMBIANO JAIME ISAZA CADAVID
FACULTAD DE INGENIERIA
INGENIERIA CIVIL
MEDELLIN
2009
3
Nota de Aceptación.
__________________________________
Firma del Jurado.
___________________________________
Firma del Jurado.
___________________________________
Firma del Jurado.
____________________________________
Firma Presidente del Jurado.
Medellín, 09 de Diciembre del 2009.
4
DEDICATORIA
A mis padres, hermanos, hija, profesor Giovanni Martínez Martínez y todo el
grupo profesoral tan de buena calidad que ayudaron a mi formación como
ingeniero civil. Por ello les estoy muy agradecido el cual solo les ofrezco el
sentimiento tan invaluable que es el de la alegría humana por el de
acompañarme educacionalmente en un proceso tan vital en mi proyecto de vida
personal y profesional.
5
AGRADECIMIENTOS
Expreso mis más sinceros agradecimientos a:
Mis padres FERNANDO ALVAREZ Y NANCY SANCHEZ que estuvieron a mí lado
durante todo el proceso de formación profesional.
A todo ese grupo de profesionales tan competente del POLITECNICO
COLOMBIANO J.I.C. que hicieron que el proceso como ingeniero civil fuera un
total éxito en especial a:
Giovanni Martínez Martínez
Martha Elena Zapata Pérez.
Daniel Zapata.
Luis Guillermo Montoya Vivas.
Santiago Wilches.
Profesionales tan íntegros tanto personal como profesionalmente, de los cuales
aprendí tanto de lo mencionado anteriormente y obviamente de su cátedra,
espero que perduren mucho más para que sigan ayudando en la formación de
nuevos profesionales en el área de la Ing. Civil.
6
CONTENIDO
Pág.
GLOSARIO. 9
RESUMEN. 10
INTRODUCCION. 11
1. OBJETIVOS. 12
1.1. OBJETIVO GENERAL. 12
1.2. OBJETIVOS ESPECIFICOS. 12
2. RESEÑA HISTORICA. 13
3. METODO MATRICIAL DE RIGIDEZ. 14
4. CONVENCION DE SIGNOS POSITIVO. 19
5. NUMERACION DE GRADOS DE LIBERTAD. 19
6. SISTEMA LOCAL Y GLOBAL DE COORDENADAS. 20
7
6.1. SISTEMA GLOBAL. 21
6.2. SISTEMA LOCAL. 21
7. MATRIZ DE UN ELEMENTO TRIDIMENSIONAL EN COORDENADAS
LOCALES. 23
7.1. PROPIEDADES DE LA MATRIZ. 24
8. MATRIZ DE TRANSFORMACION DE COORDENADAS [ λ ]. 25
8.1. MATRIZ DE TRANSFORMACION PORTICO PLANO. 26
8.2. MATRIZ DE TRANSFORMACION CERCHA PLANA. 27
8.3. MATRIZ DE TRANSFORMACION CERCHA ESPACIAL. 27
9. DEMOSTRACION DE LA OBTENCION DE LA ECUACION GENERAL
{ F } = [ K ] * { U }. 28
10. DESGLOSE DE LA ECUACION GENERAL. 30
11. MATRIZ EN COORDENADAS GLOBALES CERCHA PLANA. 31
12. MATRIZ EN COORDENADAS GLOBALES PORTICO PLANO. 32
13. MATRIZ EN COORDENADAS GLOBALES ELEMENTO VIGA. 33
8
14. MODIFICACION DE LA ECUACON GENERAL CUANDO SE TIENEN CARGAS
EN LAS LUCES. 34
15. MOMENTOS DE EMPOTRAMIENTO Y REACCIONES. 35
16. CONDENSACION MATRICIAL. 37
17. EJEMPLOS DE APLICACIÓN. 39
17.1. EJEMPLO 1. VIGA. 39
17.2. EJEMPLO 2. VIGA. 47
17.3. EJEMPLO 3. PORTICO. 56
17.4. EJEMPLO 4. PORTICO. 69
17.5. EJEMPLO 5. PORTICO 3D. 79
18. CONCLUSIONES. 102
19. RECOMENDACIONES. 103
20. BIBLIOGRAFIA Y CIBERGRAFIA. 104
9
GLOSARIO
FLEXIBILIDAD: alargamiento o giro producido por una fuerza o par unidad.
GRADO DE LIBERTAD: es un posible movimiento de un nudo en una estructura.
PORTICO: Un pórtico es un espacio arquitectónico conformado por una galería
de columnas adosada a un edificio, abierta al aire libre, y situado generalmente
ante su acceso principal.
RIGIDEZ: fuerza o par, que aparece ante un alargamiento o giro unitario.
VIGA: En ingeniería y arquitectura se denomina viga a un elemento constructivo lineal que trabaja principalmente a flexión. En las vigas la longitud predomina sobre las otras dos dimensiones y suele ser horizontal.
10
RESUMEN
Los siguientes temas serán profundizados para la debida interpretación y
aplicación del método matricial de rigidez:
Se comenzara con una pequeña reseña histórica del método a tratar. Se ampliaran cada uno de los conceptos del tema como grados de
libertad, grados de libertad restringidos, grados de libertad libres, modulo de elasticidad, inercia, etc. Todos los conceptos implícitos en el método.
Se trataran las matrices para sistemas globales y locales, que son y cuando se emplean.
Matriz de rigidez y sus propiedades. Matrices de transformación de coordenadas.
Estos serán los temas que comprenderá el método, el cual será dividido en dos
partes para su mejor comprensión las cuales serán:
Con cargas en los nudos. Cargas en las luces.
Cada parte con sus respectivos ejemplos de pórtico, cercha y vigas.
Y el tema más importante que es la condensación matricial ya que mediante este
podemos hallar nuestros desplazamientos de piso cuando nos encontramos
modelando un pórtico muy grande.
11
INTRODUCCION
El método matricial de rigidez, es un método que evoluciono tanto, que en la
actualidad tiene una teoría ampliamente fundamentada con unas bases
definidas y estructuradas lo cual hace de este método un camino para la
implementación de software de modela miento estructural.
Mediante la aplicación de los siguientes ejercicios paso a paso, se pretende dar
una herramienta básica en la solución de pórticos, vigas y cerchas estructurales
mediante la aplicación del método matricial de rigidez. También será a la vez un
asesor a la hora de estudiar para un parcial ya que se enfocara detalladamente a
cada uno de los pasos que se deben efectuar a la hora de solucionar una
estructura, todo ello se llevara a cabo mediante la planteacion y solución de
cada tipo de estructuras como son cerchas, vigas y pórticos.
Todo lo anterior es de suma importancia ya que el método matricial de rigidez es
uno de los temas más importantes de los que abarca el análisis estructural ya
que el 100 % de modeladores de sistemas estructurales se basan en el método,
como por ejemplo el SAP2000.
Se llego a la determinación de los ejercicios paso a paso por la falta de un
manual en el que se explique paso a paso la solución o el cálculo de estructuras
esqueletales como pórticos, cerchas y vigas estructurales aplicando el método
matricial de rigidez, ya que en general solo se usa este método.
12
1. OBJETIVOS
1.1. OBJETIVO GENERAL:
Realizar una guía práctica para la debida interpretación, análisis y aplicación del
método matricial de rigidez para el cálculo de estructuras esqueletales como
pórticos, vigas y cerchas, mediante la aplicación paso a paso de método
matricial de rigidez.
1.2. OBJETIVOS ESPECIFICOS:
Plantear teoría de cómo calcular paso a paso estructuras esqueletales como pórticos, cerchas y vigas mediante el método matricial de rigidez.
Aplicar la teoría especificada al cálculo de ejercicios paso a paso a pórticos, y vigas con cargas distribuidas en las luces.
13
2. RESEÑA HISTORICA
Henry Manderla fue el primero en utilizar los desplazamientos (∆); y rotaciones
(Ѳ) en los nudos como incógnitas en el análisis de una estructura hiperestática.
En 1880 analizo un pórtico de nudos rígidos tomando en consideración las
deformaciones producidas en los elementos de la estructura por la acción de los
momentos flectores y las fuerzas axiales.
Esta técnica no resulto apropiada para la época por la complejidad del sistema
resultante de ecuaciones, expresado en términos de la translación y rotación
desconocidas de los nudos y que pretende describir el efecto de la flexión y de
la fuerza axial sobre cada elemento.
Posteriormente en 1892 Otto Mohr quien había contribuido al desarrollo del
método de flexibilidad para estructuras hiperestáticas, propuso un método
aproximado para el cálculo de los esfuerzos producidos por la flexión en un
pórtico de nudos rígidos. La técnica de Mohr requería la solución de un sistema
de ecuaciones expresado únicamente en términos de las rotaciones (Ѳ) de los
nudos.
En 1914 Alex Bendixen propuso el método pendiente-deflexión para el análisis
de estructuras que requieren la solución de un sistema de ecuaciones
expresado en términos de los desplazamientos (∆) y rotaciones (Ѳ) de los nudos.
En 1915 G. A. Maney dio a conocer el desarrollo formal de las ecuaciones
pendiente-deflexión. El método pendiente-deflexión propuesto por Bendixen y
Maney es semejante al método propuesto anteriormente por Mohr.
En 1930 Hardy Cross difundió el método de distribución de momentos, este
método aproxima progresivamente el valor de los momentos no equilibrados en
los nudos permitiendo de esta forma analizar estructuras planas con nudos
rígidos esta técnica tuvo gran aceptación por cuanto elimino la necesidad de
resolver el sistema de ecuaciones simultaneas lineales requerido en el método
pendiente –deflexión . E método pendiente –deflexión para el análisis de
estructuras hiperestáticas es el predecesor del método más generalizado de
análisis que se utiliza actualmente.
El advenimiento del computador digital para realizar operaciones matemáticas
elimino a la solución de ecuaciones simultáneas como una restricción u
obstáculo para el análisis estructural. Esta ha permitido la utilización de un
método muy general para el análisis de estructuras reticulares (formada por
barras esqueletales).
Las incógnitas de su formulación son los desplazamientos y las rotaciones de
los nudos. Este método de análisis se llama METODO MATRICIAL DE RGIDEZ.
14
3. METODO MATRICIAL DE RIGIDEZ
Comenzaremos conociendo que significa grado de libertad (GDL): Un grado de
libertad es un posible movimiento de un nudo en una estructura.
De este significado se desglosa:
Grados de libertad restringidos (GDR): Son aquellos que impiden el
movimiento de los nudos. Estos no lo dan por general los apoyos de la
estructura.
Grados de libertad libres (GDL): Son aquellos que se desplazan
libremente por lo general son los que no tienen apoyo.
CERCHA PLANA
Ejemplo 1.
La cercha sometida a cualquier tipo de carga.
No de nudos = 8
Un nudo de una cercha tiene 2 posibles movimientos que son Horizontal y
Vertical. Es decir que por nudo tiene 2 GDL.
GDLTotales = No * 2 = 8*2 =16 Ello para cualquier cercha plana.
Δy Δx
15
GDLLibres = 6*2 = 12.
GDLRestringidos = 2*2 = 4.
Ejemplo 2.
No de nudos = 13, ya que cuando especifican los cruces de barras con círculos
es porque es un nudo de lo contrario no sería nudo.
GDLTotales = 2*13 = 26.
GDLLibres = (2*11)+1 = 23.
GDLRestringidos = 3.
Sabemos que el triangulo es un apoyo que me restringe el desplazamiento
horizontal y vertical y el circulo es un apoyo simple que solo restringe el
desplazamiento vertical.
CERCHA ESPACIAL
Para una cercha en el espacio los GDLTotales = No de nudos * 3
Δx Δz
Δy
16
PORTICO PLANO (FRAME PLANE)
Los posibles movimientos que sufre un nudo de un pórtico plano son:
∆x = desplazamiento horizontal.
∆y = desplazamiento vertical.
Ѳz = la rotación con respecto al eje z.
GDLTotales = No de nudos * 3
Para el ejemplo GDLT=8*3 =24, GDLR=2*3=6 y los GDLL=6*2=18. Recordemos
que un empotramiento restringe todos los posibles movimientos de un nudo en
una estructura.
θZ
θz Δy
Δx
17
PORTICO ESPACIAL
Los posibles movimientos de un nudo en un pórtico espacial son 6 que son las
tres rotaciones y los tres desplazamientos:
∆x = desplazamiento horizontal.
∆y = desplazamiento vertical.
∆z = desplazamiento con respecto al eje z.
θy
Δy
θx
θz
Δz
Δx
18
Ѳx = rotación alrededor del eje x.
Ѳy = rotación alrededor del eje y.
Ѳz = rotación alrededor del eje z.
Los GDLTotales = No de nudos *6.
Para el ejemplo GDLT=36*6=216, GDLL=27*6=162 y los GDLR=9*3=27.
ELEMENTO VIGA
Las vigas se diferencian del pórtico porque siempre es recta por lo tanto como la
estructura debe ser estable no va a ver desplazamiento horizontal (∆x); es decir
que la axial se desprecia. Por lo tanto para el conteo de los GDL no se tiene en
cuenta el ∆x.
Los GDLTotales = No de nudos*2
Para el ejemplo GDLT=14, GDLR=5 y GDLL=9.
La clave es primero observar los tipos de apoyo y sacar de acuerdo a eso los
GDLRestringidos y por ultimo restarlos a los GDLT para así obtener los GDLL.
Δy
Θz
19
4. CONVENSION DE SIGNOS POSITIVA
La convención que manejaremos de ahora en adelante será:
5. NUMERACION DE LOS GRADOS DE LIBERTAD
Primero se numeran los grados de libertad libres y posteriormente los grados de
libertad restringidos y se hará en el siguiente orden.
Primero el horizontal, luego el vertical y por último la rotación.
+
3
2
1
20
6. SISTEMA LOCAL Y GLOBAL DE COORDENADAS ENUMERACION DE LOS
GRADOS DE LIBERTAD EN AL ESTRUCTURA
PORTICO PLANO
Primero numeramos los grados de libertad libres en la convención ya estipulada
anteriormente y posteriormente los restringidos. (Rojos son los grados de
libertad restringidos).
GDLT=12
GDLR=4
GDLL=8
La matriz total será de 12*12 y la matriz de incógnitas será de 8*8.
8
11
12
7
9
10
3
1
2 5 6
4
21
6.1. SISTEMA GLOBAL
Esta nos hace referencia a los ejes X, Y y Z del plano cartesiano. Este sistema es
utilizado para la ubicación de los desplazamientos en la estructura. Se denota
con letras mayúsculas. (F, U)
F = Fuerzas
U = Desplazamientos
6.2. SISTEMA LOCAL
Este nos hace referencia con el eje X paralelo al eje geométrico del elemento, es
decir que se realiza un giro al sistema global. Este sistema se utiliza para los
diagramas de fuerzas internas en las estructuras. Se denota con letras
minúsculas (f, u).
f = Fuerzas
u = Desplazamientos
Cuando se tienen elementos verticales u horizontales el sistema global coincide
con el sistema local.
f12 f8
f11
f5 f6
f4
F12 F8
F11
F6
F4
LOCAL GLOBAL
22
Lo contrario ocurre con los elementos inclinados.
Solo los momentos son los que son iguales. Es decir f7=f7 y f3=f3.
Lo que realiza después es un análisis de elemento por elemento para ir
analizando de acuerdo a la resistencia de materiales cada grado de libertad y así
aplicar principio de superposición para la extracción de la matriz de rigidez en
coordenadas locales de un elemento tridimensional.
Lo anterior no es necesario para el manual ya que este se enfocara directamente
a la aplicación directa de las matrices por lo tanto se irán dando a medida que
evolucionamos en el método.
F2 F3
F1
F5 F6
F4 GLOBAL
LOCAL
f2 f3
f1
f5 f6
f4
F
F10 F7
F9
F2
F3
f1 f3
f2
f9
f7
f10
GLOBAL LOCAL
23
7. MATRIZ DE RIGIDEZ DE UN ELEMENTO TRIDIMENSIONAL EN COORDENADAS LOCALES
1 2 3 4 5 6 7 8 9 10 11 12
f1 AE / L 0 0 0 0 0 -AE / L 0 0 0 0 0 u1
f2 0 12EIz / L3 0 0 0 6Eiz / L2 0
-12Eiz / L3 0 0 0 6Eiz / L2 u2
f3 0 0 12Eiy / L3 0 -6Eiy / L2 0 0 0 -12Eiy / L3 0 -6Eiy / L2 0 u3
f4 0 0 0 G Ix / L 0 0 0 0 0 -G Ix / L 0 0 u4
f5 0 0 -6Eiy / L2 0 4Eiy / L 0 0 0 6Eiy / L2 0 2Eiy / L 0 u5
f6 0 6Eiz / L2 0 0 0 4Eiz / L 0 -6Eiz / L2 0 0 0 2Eiz / L u6
f7 -AE / L 0 0 0 0 0 AE / L 0 0 0 0 0 u7
f8 0 -12Eiz / L3 0 0 0
-6Eiz / L2 0
12Eiz / L3 0 0 0 -6Eiz / L2 u8
f9 0 0 -12Eiy / L3 0 6Eiy / L2 0 0 0 12Eiy / L3 0 6Eiy / L2 0 u9
f10 0 0 0 -G Ix / L 0 0 0 0 0 G Ix / L 0 0 u10
f11 0 0 -6Eiy / L2 0 2Eiy / L 0 0 0 6Eiy / L2 0 4Eiy / L 0 u11
f12 0 6Eiz / L2 0 0 0 2Eiz / L 0 -6Eiz / L2 0 0 0 4Eiz / L u12
{ f } = [ k ] * { u }
12*1 12*1 12*12
24
{ f } = Vector de cargas en los nudos.
{ u } = Vector de desplazamientos.
[ k ] = Matriz de rigidez
7.1. PROPIEDADES DE LA MATRIZ DE RIGIDEZ
Simétrica.
Si aparece el modulo de elasticidad E. El material se comporta en el rango
elástico lineal, es decir que los desplazamientos y las rotaciones son
pequeñas.
No considera el efecto de la cortante.
No considera el efecto del pandeo.
No considera la rigidez de la unión viga-columna (nudos).
Todos los términos de la diagonal y tienden hacer los mayores valores de
cada fila.
Es una matriz singular. Es decir que no tiene inversa.
25
8. MATRIZ DE TRANSFORMACION [ λ ]
De la anterior matriz se extraen las matrices correspondientes de elemento
pórtico plano, elemento cercha plana y elemento viga en coordenadas locales.
Que para el manual tampoco es necesario plantearlas ya que las que se utilizan
son e coordenadas globales.
Ello se obtiene de transformar las matrices de coordenadas locales a globales.
Para lograr lo anterior se debe calcular primero una matriz que transforme las
coordenadas. Esta matriz se llama MATRIZDE TRANSFORMACION DE
COORDENADAS [ λ ].
PORTICO PLANO
F1 = f1*cosα – f2*senα
F2 = f1*senα + f2*cosα
F3 = f3
F4 = f4cosα – f5*senα
F5 = f4*senα – f5cosα
F6 = f6
Organizando matricialmente considerando cosα = cx y senα = cy obtenemos:
F2
f1
F1
F2
F5
f4
F4
f5
α
f3
F3
F6
f6
26
8.1. MATRIZ DE TRANSFORMACION DE COORDENADAS ELEMENTO PORTICO
PLANO
cx -cy 0 0 0 0
cy cx 0 0 0 0
0 0 1 0 0 0
λ = 0 0 0 cx -cy 0
0 0 0 cy cx 0
0 0 0 0 0 1
{ F } = [ λ ] * { f } (1)
{ f } = [ λ ] inv * { F } (2)
Las anteriores formulas son las que se utilizan para la transformación de las
coordenadas.
El ángulo α se mide de local a global (con respecto al eje x). Si es horario α es
negativo y si es anti horario α es positivo.
y
x
y
x
y
x
- α
+
- α
α
Xf - xi
Yf - yi
(xi,yi)
(xf,yf)
27
Xf = Coordenada x nudo final
Xi = Coordenada x nudo inicial
Yf = Coordenada y nudo final
Yi = Coordenada y nudo inicial
Cos α = (xf – xi)/L
Sen α = (yf – yi)/L
L = {(xf – xi) 2 + (yf – yi)2}^(1/2)
8.2. MATRIZ DE TRANSFORMACION DE COORDENADAS [ λ ] CERCHA PLANA
cx 0 0 0
λ = cy 0 0 0
0 0 cx 0
0 0 cy 0
8.3. MATRIZ DE TRANSFORMACION [ λ ] CERCHA ESPACIAL
CX 0 0 0 0 0
CY 0 0 0 0 0
λ = CZ 0 0 0 0 0
0 0 0 CX 0 0
0 0 0 CY 0 0
0 0 0 CZ 0 0
28
cx, cy ,cz son los cosenos directores
Lo que se procede a realizar es transformar cada una de las matrices de
elemento pórtico plano, elementó cercha plana y elemento viga, en matrices en
coordenadas globales. Cabe aclarar que las matrices en coordenadas locales
salen de la matriz de rigidez del elemento tridimensional dada anteriormente que
para el manual no es necesario dicha demostración ya que como se menciono
anteriormente el enfoque es dedicado a la aplicación directa del método, por lo
tanto a continuación se entregan las matrices de cada uno de los elementos en
coordenadas globales que son las que se utilizan a la hora de abarcar un
ejercicio.
29
9. DEMOSTRACION DE LA OBTENCION DE LA ECUACION GENERAL
{ F } = [ K ] * { U }
{ f } = [ k ] * { u } (1)
{ F } = [ λ ] * { f ] (2)
Como { F } es una cantidad vectorial, también se aplica para los desplazamientos
{ U } = [ λ ] * { u } (3)
Multiplicando por [ λ ]-1
[ λ ]-1 *{ U } = [ λ ]-1 * [ λ ] *{ u }
[ λ ]-1 * { U } = { u } (4)
4 en 1
{ f } = [ k ] * [ λ ]-1 * { U } (5)
5 en 2
{ F } = [ λ ] * [ k ] * [ λ ] * { U }
{ F } = [ K ] * { U } GLOBAL
Esta fórmula se desglosa en dos de la siguiente manera:
P0 KO K1 U0
P1 K2 K3 U1
Para la partición es de acuerdo al número de GDL libres.
K
30
10. DESGLOCE DE LA ECUACION GENERAL
10.1. ECUACION 1
Con esta calculamos los desplazamientos en los grados de libertad libres.
{ P0 } = [ K0 ] * { U0 } + [ K1 ] *{ U1 }
10.2. ECUACION 2
Con esta se calculan las reacciones
{ F1 } = [ K2 ] * { U0 } + [ K3 ] * { U1 }
DONDE:
{ F } = Vector de cargas en los grados de libertad libres.
{ F0 } = Vector de cargas en los grados de libertad libres y son conocidos.
{ F1 } = Vector de cargas en los grados de libertad restringidos que
corresponden a las reacciones.
{ U } = Vector de de desplazamientos nodales. Igual al número de GDL totales.
{ U0 ] = Vector de desplazamientos en los GDL libres y son desconocidos.
{ U1 } = Vector de desplazamientos en los GDL restringidos, son conocidos y
puede que sean cero o diferente de cero pero siempre conocidos.
31
11. MATRIZ DE RIGIDEZ EN COORDENADAS GLOBALES ELEMENTO CERCHA
PLANA
cos2 α cos α * sen α -cos2 α -cos α * sen α
cos α * sen α sen2 α
-cos α * sen α -sen2 α
K = A * E /L -cos2 α -cos α * sen α cos2 α cos α * sen α
-cos α * sen α -sen2 α cos α * sen α sen2 α
cos α = (xf - xi) / L
L = {(xf - xi ) 2 + (yf - yi)2 }^(1/2)
sen α = ( yf - yi) / L
32
12. MATRIZ DE RIGIDEZ EN COORDENADAS GLOBALES ELEMENTO PORTICO
PLANO
a b -c -a -b -c
b f h -b -f h
K = -c h k c -h k/2
-a -b c a b c
-b -f -h b f -h
-c h k/2 c -h k/2
cos α = ( xf - xi )/ L
sen α = ( yf - yi )/ L
L = { ( xf - xi )2 + ( yf - yi )2 ) }^(1/2)
a =( E*A/L)*cos2α + (12*E*I)/L3)*sen2α
b = { [(E*A)/L] - [(12*E*I)/L3] } * cosα * senα
c =[ (6*E*I)/L2 ] * senα
f =( E*A/L)*sen2α + (12*E*I)/L3)*cos2α
h =[ (6*E*I)/L2 ] * cosα
k = (4*E*I)/L
33
13. MATRIZ DE RIGIDEZ EN COORDENADAS GLOBALES ELEMENTO VIGA
12*E*Iz /L3 6*E*Iz /L2 -12*E*Iz /L3 6*E*Iz /L2
K = 6*E*Iz /L2 4*E*Iz /L -6*E*Iz /L2 2*E*Iz /L
-12*E*Iz /L3 -6*E*Iz /L2 12*E*Iz /L3 -6*E*Iz /L2
6*E*Iz /L2 2*E*Iz /L -6*E*Iz /L2 4*E*Iz /L
34
14. MODIFICACION DE LA ECUACION GENERAL CUANDO TENEMOS CARGAS
EN LAS LUCES
Mediante las ecuaciones obtenidas anteriormente solo se tenían cargas en los
nudos el cual se tenía lo siguiente:
{ F } = [ K ] * { U }
Ahora con cargas en las luces simplemente lo que se hace es agregarle el vector
de cargas en las luces el cual se denotara { R }, la formula general quedaría
entonces de la siguiente manera:
{ F } = { R } + [ K ] * { U }
{ F } – { R } = [ K ] * { U }
{ P } = [ K ] * { U }
Donde { P } = { F } – { R }
Para la de determinación o cálculo de las cargas en las luces { R }, se deben
emplear los momentos de empotramiento y reacciones.
35
15. MOMENTOS DE EMPOTRAMIENTO Y REACCIONES
15.1. CARGA RECTANGULAR
15.2. CARGA PUNTUAL
L
QL2/12
QL/2 QL/2
QL2/12
2
Q
Q
L
Qab2/L2
Qb2/L3*(3a+b) Qa2/L3*(a+3b)
Qa2b/L2
2
b a
36
15.3. CARGA TRIANGULAR
15.4. CARGA TRAPEZOIDAL
QL2/30
3QL/20L 7QL/20
QL2/20
L
L
(L2/60)*(3*Q1+2*Q2)
(L/20)*(7*Q1+3*Q2) (L/20)*(3*Q1+7*Q2)
(L2/60)*(2*Q1+3*Q2) Q1
Q2
37
16. CONDENSACION MATRICIAL
Condensar matricialmente, se refiere a reducir los grados de libertad, cuando se
considera que la diferencia entre los desplazamientos de los nudos de un mismo
nivel son muy parecidos.
Por lo general, lo que se condensa matricialmente son los desplazamientos
horizontales, es decir que solo habrá un grado de libertad por piso, cuando se
realiza manualmente ( y con la ayuda de una hoja de cálculo) es recomendable
enumerar primero los grados de libertad horizontales de los pisos.
Con esta numeración se garantiza tener los grados de libertad horizontales al
principio de la matriz principal o total.
Si solamente esos grados de libertad tienen cargas horizontales, se
tendría:
F0 KO' K1' U0'
0 K2' K3' U1'
{ U0’ } = Desplazamientos de piso.
[ K0’ ] = Igual número de pisos.
{ F0 } = [ K0’ ] * { U0’ } + [ K1’ ] * { U1’ ] (1)
{ 0 } = [ K2’ ] * { U0’ } + [ K3’ ] * { U1’ ] (2)
De (2)
[ K3’ ] * { U1’ } = - [ K2’ ] * { U0’ }
K0 ‘
K2 ‘ K3 ‘
K1 ‘
K =
38
{ U1’ } = -( [ K3’ ] -1) * [ K2’ ] * { U0’ } (3)
(3) en (1)
{ F0 } = [ K0’ ] * { U0’ } – [ K1’ ]* ([ K3’ ]-1) * [ K2’ ] * { U0’ }
{ F0 } =( [ K0’ ] – [ K1’ ]* ([ K3’ ]-1) * [ K2’ ] ) * { U0’ }
La matriz [ Kc ] es una de las matrices más utilizadas en la modelación
estructural. Esta forma de condensar o modelar es lo que se conoce como
“ DIAFRAGMA RIGIDO “, donde todos los puntos de un mismo nivel se
desplaza horizontalmente lo mismo.
[ KC ]
39
17. EJEMPLOS DE APLICACIÓN DEL METODO MATRICIAL DE RIGIDEZ
17.1. EJEMPLO 1.
Calcular los desplazamientos y las reacciones de la viga con rotación y
desplazamiento cuya sección es de 30 cm * 30 cm con un EI = cte. de 3200 Ton-
m2
PASO 1. Enumeración de nudos, sentido de análisis de las barras o tramos de
viga.
PASO 2. Enumeración de los grados de libertad, COLOR ROJO SON LOS
RESTRINGIDOS.
2,5 T/m 1 T/m
4 T/m
2 T/m
1 cm 3 cm
2 m 5 m 4 m 4 m
3 Ton 3 T-m
0,02 Rad
1 5 3 2 6
4 10 9 8
7
4 4
5 3
3 2
2
2 m 5 m 4 m 4 m
1 1
40
PASO 3. Calculo de las matrices elementales. Ki
ELEMENTO O TRAMO 1
L = 5,0 m
B = 0,3 H = 0,4 A = 0,1 m2
E = 2000000,0 I= 0,0 E*I = 3200,0 Ton-m2
6
7 8
1
307,20 768,00 -307,20 768,00 6
K 1 = 768,00 2560,00 -768,00 1280,00 7
-307,20 -768,00 307,20 -768,00 8
768,00 1280,00 -768,00 2560,00 1
ELEMENTO O TRAMO 2
L = 4,0 m
B = 0,3 H = 0,4 A = 0,1 m2
E = 2000000,0 I= 0,0 E*I = 3200,0 Ton-m2
8
1 9
2
600,00 1200,00 -600,00 1200,00 8
K 2 = 1200,00 3200,00 -1200,00 1600,00 1
-600,00 -1200,00 600,00 -1200,00 9
1200,00 1600,00 -1200,00 3200,00 2
ELEMENTO O TRAMO 3
L = 4,0 m
B = 0,3 H = 0,4 A = 0,1 m2
E = 2000000,0 I= 0,0 E*I = 3200,0 Ton-m2
9
2 10
3
600,00 1200,00 -600,00 1200,00 9
K 3 = 1200,00 3200,00 -1200,00 1600,00 2
-600,00 -1200,00 600,00 -1200,00 10
1200,00 1600,00 -1200,00 3200,00 3
41
ELEMENTO O TRAMO 4
L = 2,0 m
B = 0,3 H = 0,4 A = 0,1 m2
E = 2000000,0 I= 0,0 E*I = 3200,0 Ton-m2
10
3 4
5
4800,00 4800,00 -4800,00 4800,00 10
K 4 = 4800,00 6400,00 -4800,00 3200,00 3
-4800,00 -4800,00 4800,00 -4800,00 4
4800,00 3200,00 -4800,00 6400,00 5
PASO 4. Calculo de la matriz total de la viga. K
1
2 3 4 5 6 7 8 9 10
5760 1600 0 0 0 768 1280 432 -1200 0 1
1600 6400 1600 0 0 0 0 1200 0 -1200 2
0 1600 9600 -4800 3200 0 0 0 1200 3600 3
0 0 -4800 4800 -4800 0 0 0 0 -4800 4
0 0 3200 -4800 6400 0 0 0 0 4800 5
K = 768 0 0 0 0 307 768 -307 0 0 6
1280 0 0 0 0 768 2560 -768 0 0 7
432 1200 0 0 0 -307 -768 907 -600 0 8
-1200 0 1200 0 0 0 0 -600 1200 -600 9
0 -1200 3600 -4800 4800 0 0 0 -600 5400 10
42
PASO 5. Calculo de los vectores de carga y desplazamientos.
VECTOR DE CARGAS NODALES VECTOR DE CARGAS EN LAS LUCES
0
0,05
-3
1,60
0
-3,47
-3
2,60
F = F0 0
R =
R0 0,93
F1 F6
R1 2,50
F7
2,08
F8
5,40
F9
10,00
F10
10,50
Para el cálculo del vector de cargas en las luces se debe aplicar los momentos
de empotramiento y reacciones de acuerdo al tipo de carga que se tenga en la
estructura así:
L
(L2/60)*(3*Q1+2*Q2)
(L/20)*(7*Q1+3*Q2) (L/20)*(3*Q1+7*Q2)
(L2/60)*(2*Q1+3*Q2) Q1
Q2
L
QL2/12
QL/2 QL/2
QL2/12
2
Q
43
TRAMO 1
TRAMO 2
TRAMO 3
TRAMO 4
4m
R1=2,13
R8=2,90 R9=4,10
R2=2,53 1T/m
2,5T/m
5m
R7:2,08
R6=2,5 R8=2,5
R1:2,08
1T/m
2m
R10=3,40 R4=2,60
2T/m
4T/m
R3=1,07 R5=0,93
4T/m
4m
R2=4,13
R9=5,90 R10=7,10
R3=4,53 2,5T/m
44
VECTOR DE DESPLAZAMIENTOS VECTOR { P } = { F } – { R }
U1
-0,05
U2
-4,60
U3
3,47
U4
-5,60
U = U0 U5
P =
P0 -0,93
U1 0
P1 F6 - 2,50
0,02
F7 - 2,08
-0,03
F8 - 5,40
-0,01
F9 - 10
0
F10 - 10,5
PASO 6. Planteamiento general del sistema de ecuaciones.
{P} = [k] *{U}
P0
KO K1 U0
P1 K2 K3 U1
ECUACION 1 Cálculo de los desplazamientos.
{P0} = [K0]*{U0} + [K1]*{U1}
{P0}– {[K1]*{U1}*{K0INV}} = {U0}
45
24,64
-24,69
U1
-0,0063
Rad
-36
31,40
U2 0,0071 Rad
K1 * U1 -12
P0 - (K1*U1) 15,47
U3= -0,0025 Rad
0
-5,60
U4 -0,0103 m
0
-0,93
U5 -0,0066 Rad
ECUACION 2 Calculo de las reacciones.
{P1} = [K2]*{U0} *[K3]*{U1}
-4,8061
24,576
F6
22,27
Ton
-8,0102
74,24
F7 68,31
Ton-m
K2 * U0 5,8136
K3 * U1 -36,576
F8= -25,36 Ton
4,5011
6
F9 20,50 Ton
0,0915
6
F10 16,59 Ton
PASO 7. Dibujo de las reacciones.
20,50 Ton 68,31 Ton-m 22,27 Ton
25,36 Ton 16,59 Ton
46
DIAGRAMAS DE CORTANTE Y MOMENTO FLEXIONANTE
2,5 T/m 1 T/m
4 T/m
2 T/m
1 cm 3 cm
2 m 5 m 4 m 4 m
3 Ton 3 T-m
0,02 Rad
3
9
68,31
5,41
8,09
17,27
5,05
30,54
10,82
13,82
15,09 7,59
22,27
11,83
V (Ton)
M (T-m)
2,1m
47
17.2. EJEMPLO 2.
a). La viga mostrada tiene una sección de 30 cm * 30 cm, y un E = 2E7 kN/m2
calcular la flecha en la mitad de la viga.
b). Calcular la sección transversal de tal forma que la flecha en la mitad sea 2
cm.
PASO 1. Enumeración de nudos, sentido de análisis de las barras o tramos de
viga.
PASO 2. Enumeración de los grados de libertad, COLOR ROJO SON LOS
RESTRINGIDOS.
40 kN/m 40 kN/m
80 kN/m
2 m 4 m 4 m 2 m
2 1
3 6 9
4 5
7
8 10
2 4 3 1
2 m 4 m 4 m 2 m
4 3 2 1 5
48
PASO 3. Calculo de las matrices elementales. Ki
ELEMENTO O TRAMO 1
L = 2,0 m
B = 0,3 H = 0,3 A = 0,1 m2
E = 20000000,0 I= 0,0 E*I = 13500,0 kN-m2
1
2 9
3
20250,00 20250,00 -20250,00 20250,00 1
K 1 = 20250,00 27000,00 -20250,00 13500,00 2
-20250,00 -20250,00 20250,00 -20250,00 9
20250,00 13500,00 -20250,00 27000,00 3
ELEMENTO O TRAMO 2
L = 4,0 m
B = 0,3 H = 0,3 A = 0,1 m2
E = 20000000,0 I= 0,0 E*I = 13500,0 kN-m2
9
3 4
5
2531,25 5062,50 -2531,25 5062,50 9
K 2 = 5062,50 13500,00 -5062,50 6750,00 3
-2531,25 -5062,50 2531,25 -5062,50 4
5062,50 6750,00 -5062,50 13500,00 5
ELEMENTO O TRAMO 3
L = 4,0 m
B = 0,3 H = 0,3 A = 0,1 m2
E = 20000000,0 I= 0,0 E*I = 13500,0 kN-m2
4
5 10
6
2531,25 5062,50 -2531,25 5062,50 4
K 3 = 5062,50 13500,00 -5062,50 6750,00 5
-2531,25 -5062,50 2531,25 -5062,50 10
5062,50 6750,00 -5062,50 13500,00 6
49
ELEMENTO O TRAMO 4
L = 2,0 m
B = 0,3 H = 0,3 A = 0,1 m2
E = 20000000,0 I= 0,0 E*I = 13500,0 kN-m2
10
6 7
8
20250,00 20250,00 -20250,00 20250,00 10
K 4 = 20250,00 27000,00 -20250,00 13500,00 6
-20250,00 -20250,00 20250,00 -20250,00 7
20250,00 13500,00 -20250,00 27000,00 8
PASO 4. Calculo de la matriz total de la viga. K
1
2 3 4 5 6 7 8 9 10 20250 20250 20250 0 0 0 0 0
1
20250 27000 13500 0 0 0 0 0
2
20250 13500 40500 -5063 6750 0 0 0
3
0 0 -5063 5063 0 5063 0 0
4
0 0 6750 0 27000 6750 0 0
5
0 0 0 5063 6750 40500 -
20250 13500
6
0 0 0 0 0 -
20250 20250 -
20250
7
0 0 0 0 0 13500 -
20250 27000
8
-20250 -20250 -
15188 -2531 5063 0 0 0
9
0 0 0 -2531 -5063 15188 -
20250 20250
10
50
PASO 5. Calculo de los vectores de carga y desplazamientos.
VECTOR DE CARGAS NODALES VECTOR DE CARGAS EN LAS LUCES
0
44,00
0
15,11
0
69,33
0
288,00
0
0,00
0
-69,33
0
44,00
F = F0 0
R =
R0 -15,11
F1 F9
R1 172,00
F10
172,00
Para el cálculo del vector de cargas en las luces se debe aplicar los momentos
de empotramiento y reacciones de acuerdo al tipo de carga que se tenga en la
estructura así:
TRAMO 1
2m
R2=15,12
R1=44 R9=49,33
R3=16 40
53,33
L
(L2/60)*(3*Q1+2*Q2)
(L/20)*(7*Q1+3*Q2) (L/20)*(3*Q1+7*Q2)
(L2/60)*(2*Q1+3*Q2) Q1
Q2
51
TRAMO 2
TRAMO 3
TRAMO 4
4m
R3=85,33
R9=122,67 R4=144
R5=92,44 53,33
80
2m
R10=49,33 R7=44
40
53,33
R6=16 R8=15,11
4m
R4=144 R10=122,67
53,33
80
R5=92,44 R6=85,33
52
VECTOR {P } = { F } – { R } VECTOR DE DESPLAZAMIENTOS
-44,00
U1
-15,11
U2
-69,33
U3
-288,00
U4
0,00
U5
69,33
U6
-44,00
U7
P = P0 15,11
U =
U0 U8
P1 F9 - 172
U1 0
F10 - 172
0
PASO 6. Planteamiento general del sistema de ecuaciones.
{P} = [k] *{U}
P0
KO K1 U0
P1 K2 K3 U1
ECUACION 1 Cálculo de los desplazamientos.
{P0} = [K0]*{U0} + [K1]*{U1}
[P0]*[K0INV] = {U0}
53
U1
0,1621
m
U2 -0,0800 Rad
U3 -0,0843 Rad
U4 -0,2254 m
U5 0,0000 Rad
U6 0,0843 Rad
U7 0,1621 m
U8 0,0800 Rad
ECUACION 2 Calculo de las reacciones.
{P1} = [K2]*{U0}
F9 - 172
188
F9
360
kN
F10 - 172 188
F10 360 kN
PASO 7. Dibujo de las reacciones.
360 kN 360 kN
54
DIAGRAMAS DE CORTANTE Y MOMENTO FLEXIONANTE
40 kN/m 40 kN/m
80 kN/m
2 m 4 m 4 m 2 m
93,33
V (kN)
M (kN-m)
266,67
266,67
88,89
Mmax=328,89
88,89
93,33
55
b). Para hallar la sección basta cambiar los datos en las matrices de cada
elemento de la viga, ya que estos se encuentran montados en una hoja de
cálculo Excel esto funciona como una iteración.
Realizando lo anterior nos da una flecha de 2 cm, con una sección transversal de
40 cm * 55 cm.
56
17.3. EJEMPLO 3.
Modulo del material concreto = 2E6 Ton/m2
Secciones transversales:
Columnas 40cm * 60cm y vigas de 40cm * 40cm
Combinaciones de carga 1). 1,4cm + 1,7cv 2). Cm + cv 3). Cm + cv + ch
Para la solución de pórtico se procede a realizar primero que todas las
combinaciones el ejercicio solo se calculara con la 1 combinación de carga.
Como ya no dieron cuánto vale cm y cv simplemente se remplazan en cada una
de las combinaciones así.
CARGAS (Ton/m) CM CV CH 1 0,6 0 COMVINACIONES (Ton/m) CARGAS R
1,4CM + 174CV
CM + CH
CM + CV + CH qL/2 qL2/12
2,42 1 1,6 4,84 3,227 CON
1 2 3 2 1,333 2
L = 4 m 3,2 2,133 3
4,840 3,227 1
3 m
3 m
4 m 4 m 4 m
5 Ton
10 Ton
0,6 Ton/m = cv
1 Ton/m = cm
0,6 Ton/m = cv
1 Ton/m = cm
57
PASO 1. Numeración de cada elemento, nudo y sentido de los elementos.
PASO 2 Enumeración de los grados de libertad, COLOR ROJO SON LOS
RESTRINGIDOS.
3 m
3 m
4 m 4 m 4 m
1
5 4 6 3
10 7 9 8
2
1 2
10 8
6 5
9
7
4 3
1 2
3
4 5
14 13
15 12
11 10
9
7
23
21
22
19
8
24
20
26
27
30
28
29
16 17
18
25
6
58
PASO 3. Cálculo de matrices elementales. Ki
ELEMENTO 1
Xi = 0,0 Xf = 0,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,6 A = 0,2 E = 2000000,0 I= 0,0 a = 6400,0 b = 0,0 c = 9600,0 f = 160000,0 h = 0,0 k = 19200,0 k/2 = 9600,0
21
22 23 7 8
9
6400,0 0,0 -9600,0 -6400,0 0,0 -9600,0 21
0,0 160000,0 0,0 0,0 -160000,0 0,0 22
-9600,0 0,0 19200,0 9600,0 0,0 9600,0 23
-6400,0 0,0 9600,0 6400,0 0,0 9600,0 7
0,0 -160000,0 0,0 0,0 160000,0 0,0 8
-9600,0 0,0 9600,0 9600,0 0,0 19200,0 9
ELEMENTO 2
Xi = 4,0 Xf = 4,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,6 A = 0,2 E = 2000000,0 I= 0,0 a = 6400,0 b = 0,0 c = 9600,0 f = 160000,0 h = 0,0 k = 19200,0 k/2 = 9600,0
24
25 19 10 11
12
6400,0 0,0 -9600,0 -6400,0 0,0 -9600,0 24
0,0 160000,0 0,0 0,0 -160000,0 0,0 25
-9600,0 0,0 19200,0 9600,0 0,0 9600,0 19
-6400,0 0,0 9600,0 6400,0 0,0 9600,0 10
0,0 -160000,0 0,0 0,0 160000,0 0,0 11
-9600,0 0,0 9600,0 9600,0 0,0 19200,0 12
59
ELEMENTO 3
Xi = 8,0 Xf = 8,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,6 A = 0,2 E = 2000000,0 I= 0,0 a = 6400,0 b = 0,0 c = 9600,0 f = 160000,0 h = 0,0 k = 19200,0 k/2 = 9600,0
26
27 20 13 14
15
6400,0 0,0 -9600,0 -6400,0 0,0 -9600,0 26
0,0 160000,0 0,0 0,0 -160000,0 0,0 27
-9600,0 0,0 19200,0 9600,0 0,0 9600,0 20
-6400,0 0,0 9600,0 6400,0 0,0 9600,0 13
0,0 -160000,0 0,0 0,0 160000,0 0,0 14
-9600,0 0,0 9600,0 9600,0 0,0 19200,0 15
ELEMENTO 4
Xi = 12,0 Xf = 12,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,6 A = 0,2 E = 2000000,0 I= 0,0 a = 6400,0 b = 0,0 c = 9600,0 f = 160000,0 h = 0,0 k = 19200,0 k/2 = 9600,0
28
29 30 16 17
18
6400,0 0,0 -9600,0 -6400,0 0,0 -9600,0 28
0,0 160000,0 0,0 0,0 -160000,0 0,0 29
-9600,0 0,0 19200,0 9600,0 0,0 9600,0 30
-6400,0 0,0 9600,0 6400,0 0,0 9600,0 16
0,0 -160000,0 0,0 0,0 160000,0 0,0 17
-9600,0 0,0 9600,0 9600,0 0,0 19200,0 18
60
ELEMENTO 5
Xi = 4,0 Xf = 4,0 Yi = 3,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,6 A = 0,2 E = 2000000,0 I= 0,0 a = 6400,0 b = 0,0 c = 9600,0 f = 160000,0 h = 0,0 k = 19200,0 k/2 = 9600,0
10
11 12 1 2
3
6400,0 0,0 -9600,0 -6400,0 0,0 -9600,0 10
0,0 160000,0 0,0 0,0 -160000,0 0,0 11
-9600,0 0,0 19200,0 9600,0 0,0 9600,0 12
-6400,0 0,0 9600,0 6400,0 0,0 9600,0 1
0,0 -160000,0 0,0 0,0 160000,0 0,0 2
-9600,0 0,0 9600,0 9600,0 0,0 19200,0 3
ELEMENTO 6
Xi = 8,0 Xf = 8,0 Yi = 3,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,6 A = 0,2 E = 2000000,0 I= 0,0 a = 6400,0 b = 0,0 c = 9600,0 f = 160000,0 h = 0,0 k = 19200,0 k/2 = 9600,0
13
14 15 4 5
6
6400,0 0,0 -9600,0 -6400,0 0,0 -9600,0 13
0,0 160000,0 0,0 0,0 -160000,0 0,0 14
-9600,0 0,0 19200,0 9600,0 0,0 9600,0 15
-6400,0 0,0 9600,0 6400,0 0,0 9600,0 4
0,0 -160000,0 0,0 0,0 160000,0 0,0 5
-9600,0 0,0 9600,0 9600,0 0,0 19200,0 6
61
ELEMENTO 7
Xi = 4,0 Xf = 8,0 Yi = 6,0 Yf = 6,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 80000,0 b = 0,0 c = 0,0 f = 800,0 h = 1600,0 k = 4266,7 k/2 = 2133,3
1
2 3 4 5
6
80000,0 0,0 0,0 -80000,0 0,0 0,0 1
0,0 800,0 1600,0 0,0 -800,0 1600,0 2
0,0 1600,0 4266,7 0,0 -1600,0 2133,3 3
-80000,0 0,0 0,0 80000,0 0,0 0,0 4
0,0 -800,0 -1600,0 0,0 800,0 -1600,0 5
0,0 1600,0 2133,3 0,0 -1600,0 4266,7 6
ELEMENTO 8
Xi = 0,0 Xf = 4,0 Yi = 3,0 Yf = 3,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 80000,0 b = 0,0 c = 0,0 f = 800,0 h = 1600,0 k = 4266,7 k/2 = 2133,3
7
8 9 10 11
12
80000,0 0,0 0,0 -80000,0 0,0 0,0 7
0,0 800,0 1600,0 0,0 -800,0 1600,0 8
0,0 1600,0 4266,7 0,0 -1600,0 2133,3 9
-80000,0 0,0 0,0 80000,0 0,0 0,0 10
0,0 -800,0 -1600,0 0,0 800,0 -1600,0 11
0,0 1600,0 2133,3 0,0 -1600,0 4266,7 12
62
ELEMENTO 9
Xi = 4,0 Xf = 8,0 Yi = 3,0 Yf = 3,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 80000,0 b = 0,0 c = 0,0 f = 800,0 h = 1600,0 k = 4266,7 k/2 = 2133,3
10
11 12 13 14
15
80000,0 0,0 0,0 -80000,0 0,0 0,0 10
0,0 800,0 1600,0 0,0 -800,0 1600,0 11
0,0 1600,0 4266,7 0,0 -1600,0 2133,3 12
-80000,0 0,0 0,0 80000,0 0,0 0,0 13
0,0 -800,0 -1600,0 0,0 800,0 -1600,0 14
0,0 1600,0 2133,3 0,0 -1600,0 4266,7 15
ELEMENTO 10
Xi = 8,0 Xf = 12,0 Yi = 3,0 Yf = 3,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 80000,0 b = 0,0 c = 0,0 f = 800,0 h = 1600,0 k = 4266,7 k/2 = 2133,3
13
14 15 16 17
18
80000,0 0,0 0,0 -80000,0 0,0 0,0 13
0,0 800,0 1600,0 0,0 -800,0 1600,0 14
0,0 1600,0 4266,7 0,0 -1600,0 2133,3 15
-80000,0 0,0 0,0 80000,0 0,0 0,0 16
0,0 -800,0 -1600,0 0,0 800,0 -1600,0 17
0,0 1600,0 2133,3 0,0 -1600,0 4266,7 18
63
PASO 4. Calculo de la matriz total de la estructura. K.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
20
86400 0 9600 -80000 0 0 0 0 0 -6400 0 9600 0 0 0 0 0 0 0 0 0 160800 1600 0 -800 1600 0 0 0 0 -160000 0 0 0 0 0 0 0 0 0 9600 1600 23467 0 -1600 2133 0 0 0 -9600 0 9600 0 0 0 0 0 0 0 0 -80000 0 0 86400 0 9600 0 0 0 0 0 0 -6400 0 9600 0 0 0 0 0 0 -800 -1600 0 160800 -1600 0 0 0 0 0 0 0 -160000 0 0 0 0 0 0 0 1600 2133 9600 -1600 23467 0 0 0 0 0 0 -9600 0 9600 0 0 0 0 0 0 0 0 0 0 0 86400 0 9600 -80000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 160800 1600 0 -800 1600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9600 1600 23467 0 -1600 2133 0 0 0 0 0 0 0 0 -6400 0 -9600 0 0 0 -80000 0 0 172800 0 0 -80000 0 0 0 0 0 9600 0 0 -160000 0 0 0 0 0 -800 -1600 0 321600 0 0 -800 1600 0 0 0 0 0 9600 0 9600 0 0 0 0 1600 2133 0 0 46933 0 -1600 2133 0 0 0 9600 0 0 0 0 -6400 0 -9600 0 0 0 -80000 0 0 172800 0 0 -80000 0 0 0 9600 0 0 0 0 -160000 0 0 0 0 0 -800 -1600 0 321600 0 0 -800 1600 0 0 0 0 0 9600 0 9600 0 0 0 0 1600 2133 0 0 46933 0 -1600 2133 0 9600 0 0 0 0 0 0 0 0 0 0 0 0 -80000 0 0 86400 0 9600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -800 -1600 0 160800 -1600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1600 2133 9600 -1600 23467 0 0 0 0 0 0 0 0 0 0 0 9600 0 9600 0 0 0 0 0 0 19200 0 0 0 0 0 0 0 0 0 0 0 0 0 9600 0 9600 0 0 0 0 19200 0 0 0 0 0 0 -6400 0 -9600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -160000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9600 0 9600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6400 0 -9600 0 0 0 0 0 0 -9600 0 0 0 0 0 0 0 0 0 0 0 -160000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6400 0 -9600 0 0 0 0 -9600 0 0 0 0 0 0 0 0 0 0 0 0 0 -160000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6400 0 -9600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -160000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9600 0 9600 0 0
64
PASO5. Calculo de los vectores de carga y desplazamientos
VECTOR DE CARGAS NODALES VECTOR DE CARGAS EN LAS LUCES
10
0
0
4,84
0
3,227
0
0
0
4,84
0
-3,227
5
0
0
4,84
0
3,227
0
0
0
9,68
0
0
0
0
0
9,68
0
0
0
0
0
4,84
0
-3,227
0
0
F = FO 0
R =
R0 0
F1 F21
R1 0
F22
0
F23
0
F24
0
F25
0
F26
0
F27
0
F28
0
F29
0
F30
0
Para el cálculo del vector de cargas en las luces se debe aplicar los momentos
de empotramiento y reacciones de acuerdo al tipo de carga que se tenga en la
estructura así:
65
VIGA 7
VIGA 8
VIGA 9
VIGA 10
L
QL2/12
QL/2 QL/2
QL2/12
2
4 m
R3=3,23 R6=3,23
R2=4,84 R5=4,84
4 m
R14=4,84 R11=4,84
R12=3,23 R15=3,23
4 m
R2=4,84 R14=4,84
R15=3,23 R18=3,23
4 m
R11=4,84 R8=4,84
R9=3,23 R12=3,23
66
VECTOR DE DESPLAZAMIENTOS VECTOR { P } = { F }-{ R }
U1
10
U2
-4,84
U3
-3,23
U4
0
U5
-4,84
U6
3,23
U7
5
U8
-4,84
U9
-3,23
U10
0
U11
-9,68
U12
0,00
U13
0,00
U14
-9,68
U15
0
U16
0
U17
-4,84
U18
3,23
U19
0
U = U0 U20
P =
P0 0
U1 0
P1 F21
0
F22
0
F23
0
F24
0
F25
0
F26
0
F27
0
F28
0
F29
0
F30
67
PASO 6. Planteamiento general del sistema de ecuaciones.
{P} = [k] *{U}
P0
KO K1 U0
P1 K2 K3 U1
ECUACION 1 Cálculo de los desplazamientos.
{P0} = [K0]*{U0} + [K1]*{U1}
Como U1 es cero por lo tanto {P0}= [K0]*{U0}
{P0}*[K0INV] = {U0}
U1
0,0061
m
U2 -0,0001 m
U3 -0,0012 Rad
U4 0,0060 m
U5 -0,0002 m
U6 -0,0008 Rad
U7 0,0026 m
U8 0,0000 m
U9 -0,0011 Rad
U10 = 0,0026 m
U11 -0,0001 m
U12 -0,0007 Rad
U13 0,0026 m
U14 -0,0001 m
U15 -0,0008 Rad
U16 0,0025 m
U17 0,0000 m
U18 -0,0008 Rad
U19 -0,0010 Rad
U20 -0,0009 Rad
68
ECUACION 2 Calculo de las reacciones.
{P1} = [K2]*{U0} *[K3]*{U1}
Como U0 es cero por lo tanto {P1} = [K2]*{U0}
F21
-5,7457
Ton
F22 1,8977 Ton
F23 14,1261 Ton*m
F24 -0,7361 Ton
F25 = 11,8745 Ton
F26 -0,2804 Ton
F27 17,4758 Ton
F28 -8,2379 Ton
F29 7,4721 Ton
F30 16,2250 Ton*m
PASO 7. Dibujo de las reacciones.
7,47 Ton
8,24 Ton
16,23 Ton-m 17,48 Ton
0,28 Ton 0,74 Ton
11,88 Ton 1,898 Ton 14,13 Ton-m
5,75 Ton
69
17.4. EJEMPLO 4.
Aplicar el método matricial al pórtico propuesto cuyas columnas y vigas tienen
una sección transversal de 40 cm * 40 cm y un E = 2E7 kN/m.
5 kN/m 5 kN/m 10 kN/m
3 m
3 m
4 m 4 m
70
PASO 1. Numeración de cada elemento, nudo y sentido de los elementos.
PASO 2 Enumeración de los grados de libertad, COLOR ROJO SON LOS
RESTRINGIDOS.
1
2 4
3
5
4 m 4 m
3 m
3 m
2 3
4 1
15
13
14 12
10
11
7
9 8
6 4
5
2 3
1
71
PASO 3. Cálculo de matrices elementales. Ki
ELEMENTO 1
Xi = 0,0 Xf = 0,0 Yi = 0,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 6,0 B = 0,4 H = 0,4 A = 0,2
E = 20000000,
0 I= 0,0 a = 2370,4 b = 0,0 c = 7111,1 f = 533333,3 h = 0,0 k = 28444,4 k/2 = 14222,2
13
14 15 1 2
3
2370,4 0,0 -7111,1 -2370,4 0,0 -7111,1 13
0,0 533333,3 0,0 0,0 -533333,3 0,0 14
-7111,1 0,0 28444,4 7111,1 0,0 14222,2 15
-2370,4 0,0 7111,1 2370,4 0,0 7111,1 1
0,0 -533333,3 0,0 0,0 533333,3 0,0 2
-7111,1 0,0 14222,2 7111,1 0,0 28444,4 3
ELEMENTO 2
Xi = 0,0 Xf = 4,0 Yi = 6,0 Yf = 9,0 Cos α = 0,8 Sen α = 0,6 L = 5,0 B = 0,4 H = 0,4 A = 0,2
E = 20000000,
0 I= 0,0 a = 411074,6 b = 305233,9 c = 6144,0 f = 233021,4 h = 8192,0 k = 34133,3 k/2 = 17066,7
1
2 3 4 5
6
411074,6 305233,9 -6144,0 -411074,6 -305233,9 -6144,0 1
305233,9 233021,4 8192,0 -305233,9 -233021,4 8192,0 2
-6144,0 8192,0 34133,3 6144,0 -8192,0 17066,7 3
-411074,6 -305233,9 6144,0 411074,6 305233,9 6144,0 4
-305233,9 -233021,4 -8192,0 305233,9 233021,4 -8192,0 5
-6144,0 8192,0 17066,7 6144,0 -8192,0 34133,3 6
72
ELEMENTO 3
Xi = 4,0 Xf = 8,0 Yi = 9,0 Yf = 6,0 Cos α = 0,8 Sen α = -0,6 L = 5,0 B = 0,4 H = 0,4 A = 0,2
E = 20000000,
0 I= 0,0 a = 411074,6 b = -305233,9 c = -6144,0 f = 233021,4 h = 8192,0 k = 34133,3 k/2 = 17066,7
4
5 6 7 8
9
411074,6 -305233,9 6144,0 -411074,6 305233,9 6144,0 4
-305233,9 233021,4 8192,0 305233,9 -233021,4 8192,0 5
6144,0 8192,0 34133,3 -6144,0 -8192,0 17066,7 6
-411074,6 305233,9 -6144,0 411074,6 -305233,9 -6144,0 7
305233,9 -233021,4 -8192,0 -305233,9 233021,4 -8192,0 8
6144,0 8192,0 17066,7 -6144,0 -8192,0 34133,3 9
ELEMENTO 4
Xi = 8,0 Xf = 8,0 Yi = 0,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 6,0 B = 0,4 H = 0,4 A = 0,2
E = 20000000,
0 I= 0,0 a = 2370,4 b = 0,0 c = 7111,1 f = 533333,3 h = 0,0 k = 28444,4 k/2 = 14222,2
10
11 12 7 8
9
2370,4 0,0 -7111,1 -2370,4 0,0 -7111,1 10
0,0 533333,3 0,0 0,0 -533333,3 0,0 11
-7111,1 0,0 28444,4 7111,1 0,0 14222,2 12
-2370,4 0,0 7111,1 2370,4 0,0 7111,1 7
0,0 -533333,3 0,0 0,0 533333,3 0,0 8
-7111,1 0,0 14222,2 7111,1 0,0 28444,4 9
73
PASO 4. Calculo de la matriz total de la estructura. K.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
413445 305234 967
-411075
-305234 -6144 0 0 0
1
305234 766355 8192
-305234
-233021 8192 0 0 0
2
967 8192 62578 6144 -8192 17067 0 0 0
3
-411075
-305234 6144 822149 0 12288
-411075 305234 6144
4
-305234
-233021 -8192 0 466043 0 305234
-233021 8192
5
-6144 8192 17067 12288 0 68267 -6144 -8192 17067
6
0 0 0
-411075 305234 -6144 413445
-305234 967
7
0 0 0 305234
-233021 -8192
-305234 766355 -8192
8
K = 0 0 0 6144 8192 17067 967 -8192 62578
9
0 0 0 0 0 0 -2370 0 -7111
10
0 0 0 0 0 0 0
-533333 0
11
0 0 0 0 0 0 7111 0 14222
12
-2370 0 -7111 0 0 0 0 0 0
13
0
-533333 0 0 0 0 0 0 0
14
7111 0 14222 0 0 0 0 0 0
15
74
PASO 5. Calculo de los vectores de cargas y desplazamientos.
VECTOR DE CARGAS NODALES VECTOR DE CARGAS EN LAS LUCES
0
-33,00
0
10
0
4,42
0
-28,50
0
0
0
-6,50
0
-7,5
0
-10
F = F0 0
R = R0 -10,42
F1 F10 R1 0
F11
0
F12
0
F13
-6,00
F14
0
F15
8,00
Para el cálculo del vector de cargas en las luces se debe aplicar los momentos
de empotramiento y reacciones de acuerdo al tipo de carga que se tenga en la
estructura así:
L
QL2/12
QL/2 QL/2
QL2/12
2
QL2/30
3QL/20L 7QL/20
QL2/20
L
75
Calculo de cada elemento con su respectiva carga distribuida.
L
(L2/60)*(3*Q1+2*Q2)
(L/20)*(7*Q1+3*Q2) (L/20)*(3*Q1+7*Q2)
(L2/60)*(2*Q1+3*Q2) Q1
Q2
3 m
R3=6
R1=11,50 R4=13,50
R6=6,50 6,67
10
R15=8
R13=6 R1=14
R3=12
6 m
6,67 kN/m
R6=10,42
12,50 R5=10
R2=10
R1=7,5
12,50
R3=10,42
R4=7,5 R4=7,5
12,50 R5=10
R7=7,5
12,50
R8=10
R9=10,42
R6=10,42
76
VECTOR {P} = [F]-{R} VECTOR DE DESPLAZAMIENTOS
33,00
U1
-10,00
U2
-4,42
U3
28,50
U4
0,00
U5
6,50
U6
7,50
U7
10,00
U8
P = P0 10,42
U =
U0 U9
P1 F10
U1 0
F11
0
F12
0
F13 + 6
0
F14
0
F15 - 8
0
PASO 6. Planteamiento general del sistema de ecuaciones.
{P} = [k] *{U}
P0
KO K1 U0
P1 K2 K3 U1
77
ECUACION 1 Cálculo de los desplazamientos.
{P0} = [K0]*{U0} + [K1]*{U1}
Como U1 es cero por lo tanto {P0}= [K0]*{U0}
{P0}*[K0INV] = {U0}
U1
0,025
m
U2 0,000 m
U3 -0,003 Rad
U4 0,024 m
U5 0,001 m
U6 0,002 Rad
U7 0,023 m
U8 0,000 m
U9 -0,003 Rad
ECUACION 2 Calculo de las reacciones.
{P1} = [K2]*{U0} *[K3]*{U1}
Como U0 es cero por lo tanto {P1} = [K2]*{U0}
F10
-32,006
F10
-32,01
kN
F11 19,428
F11 19,43 kN
F12 118,531
F12 118,53 kN-m
F13 + 6 -36,994
F13 -42,99 kN
F14 -19,428
F14 -19,43 kN
F15 - 8 133,049
F15 141,05 kN-m
78
PASO 7. Dibujo de las reacciones.
118,53 kN-m 19,43 kN
32,01 kN
141,05 kN-m 19,43 kN
42,99 kN
79
17.5. EJEMPLO 5.
Aplicar el método matricial y condensar matricialmente para calcular los
desplazamientos por piso sabiendo que hay unas fuerzas por piso en dirección x
y dirección y, mostradas en la figura.
E = 2E6Ton/m2 sección de columnas 40 * 40 y sección de vigas 40 * 30.
4m 4m 4m
4m
4m
4m
4m
3m
3m
3m
3m
3m
30ton
40ton
5ton
20ton
10ton
40ton
30ton
20ton
10ton
5ton
80
PASO 1. Numeración de cada elemento, nudo y sentido de los elementos.
4m 4m 4m
2
8
9
21
1
3
32
23
18
7
13
14
26
22
24
19 12
10 11 20
35 33 34
5 6
4 17
16
29 28 27
30 31
25
15
11
12
3
2
1
14
13
2
2
23
24
4 9 16 21
15 10
19
5 20 8 17
6 7 18
81
PASO 2 Enumeración de los grados de libertad, COLOR ROJO SON LOS
RESTRINGIDOS.
PORTICO PLANO CONDENSADO
10 11
35 37
2
8 9
24 25 28 29 26 26
16 17 20 21 18 19
12 13
44 45
2
32 33 34 35
42 43
2 40
2
41
5
6 7
46
47 48
1
38 39
2
4
14 15
2
30 31
3
22 23
49
50 51
52
53 54
55
56 57
82
PASO 3. Cálculo de matrices elementales. Ki
ELEMENTO 1 columna Xi = 0,0 Xf = 0,0
Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
46
47 48 1 38
39
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 46
0,0 106666,7 0,0 0,0 -106666,7 0,0 47
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 48
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 1
0,0 -106666,7 0,0 0,0 106666,7 0,0 38
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 39
ELEMENTO 2 columna Xi = 0,0 Xf = 0,0
Yi = 3,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
1
38 39 2 30
31
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 1
0,0 106666,7 0,0 0,0 -106666,7 0,0 38
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 39
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 2
0,0 -106666,7 0,0 0,0 106666,7 0,0 30
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 31
83
ELEMENTO 3 columna
Xi = 0,0 Xf = 0,0 Yi = 6,0 Yf = 9,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
2
30 31 3 22
23
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 2
0,0 106666,7 0,0 0,0 -106666,7 0,0 30
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 31
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 3
0,0 -106666,7 0,0 0,0 106666,7 0,0 22
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 23
ELEMENTO 4 columna Xi = 0,0 Xf = 0,0
Yi = 9,0 Yf = 12,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
3
22 23 4 14
15
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 3
0,0 106666,7 0,0 0,0 -106666,7 0,0 22
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 23
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 4
0,0 -106666,7 0,0 0,0 106666,7 0,0 14
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 15
84
ELEMENTO 5 columna Xi = 0,0 Xf = 0,0
Yi = 12,0 Yf = 15,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
4
14 15 5 6
7
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 4
0,0 106666,7 0,0 0,0 -106666,7 0,0 14
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 15
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 5
0,0 -106666,7 0,0 0,0 106666,7 0,0 6
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 7
ELEMENTO 6 columna Xi = 4,0 Xf = 4,0
Yi = 12,0 Yf = 15,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
4
16 17 5 8
9
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 4
0,0 106666,7 0,0 0,0 -106666,7 0,0 16
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 17
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 5
0,0 -106666,7 0,0 0,0 106666,7 0,0 8
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 9
85
ELEMENTO 7 columna
Xi = 4,0 Xf = 4,0 Yi = 9,0 Yf = 12,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
3
24 25 4 16
17
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 3
0,0 106666,7 0,0 0,0 -106666,7 0,0 24
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 25
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 4
0,0 -106666,7 0,0 0,0 106666,7 0,0 16
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 17
ELEMENTO 8 columna Xi = 4,0 Xf = 4,0
Yi = 6,0 Yf = 9,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
2
32 33 3 24
25
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 2
0,0 106666,7 0,0 0,0 -106666,7 0,0 32
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 33
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 3
0,0 -106666,7 0,0 0,0 106666,7 0,0 24
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 25
86
ELEMENTO 9 columna
Xi = 4,0 Xf = 4,0 Yi = 3,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
1
40 41 2 32
33
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 1
0,0 106666,7 0,0 0,0 -106666,7 0,0 40
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 41
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 2
0,0 -106666,7 0,0 0,0 106666,7 0,0 32
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 33
ELEMENTO 10 columna
Xi = 4,0 Xf = 4,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
49
50 51 1 40
41
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 49
0,0 106666,7 0,0 0,0 -106666,7 0,0 50
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 51
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 1
0,0 -106666,7 0,0 0,0 106666,7 0,0 40
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 41
87
ELEMENTO 11 columna
Xi = 8,0 Xf = 8,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
52
53 54 1 42
43
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 52
0,0 106666,7 0,0 0,0 -106666,7 0,0 53
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 54
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 1
0,0 -106666,7 0,0 0,0 106666,7 0,0 42
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 43
ELEMENTO 12 columna
Xi = 8,0 Xf = 8,0 Yi = 3,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
1
42 43 2 34
35
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 1
0,0 106666,7 0,0 0,0 -106666,7 0,0 42
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 43
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 2
0,0 -106666,7 0,0 0,0 106666,7 0,0 34
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 35
88
ELEMENTO 13 columna
Xi = 8,0 Xf = 8,0 Yi = 6,0 Yf = 9,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
2
34 35 3 26
27
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 2
0,0 106666,7 0,0 0,0 -106666,7 0,0 34
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 35
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 3
0,0 -106666,7 0,0 0,0 106666,7 0,0 26
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 27
ELEMENTO 14 columna
Xi = 8,0 Xf = 8,0 Yi = 9,0 Yf = 12,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
3
26 27 4 18
19
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 3
0,0 106666,7 0,0 0,0 -106666,7 0,0 26
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 27
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 4
0,0 -106666,7 0,0 0,0 106666,7 0,0 18
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 19
89
ELEMENTO 15 columna
Xi = 8,0 Xf = 8,0 Yi = 12,0 Yf = 15,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
4
18 19 5 10
11
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 4
0,0 106666,7 0,0 0,0 -106666,7 0,0 18
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 19
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 5
0,0 -106666,7 0,0 0,0 106666,7 0,0 10
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 11
ELEMENTO 16 columna
Xi = 12,0 Xf = 12,0 Yi = 12,0 Yf = 15,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
4
20 21 5 12
13
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 4
0,0 106666,7 0,0 0,0 -106666,7 0,0 20
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 21
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 5
0,0 -106666,7 0,0 0,0 106666,7 0,0 12
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 13
90
ELEMENTO 17 columna
Xi = 12,0 Xf = 12,0 Yi = 9,0 Yf = 12,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
3
28 29 4 20
21
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 3
0,0 106666,7 0,0 0,0 -106666,7 0,0 28
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 29
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 4
0,0 -106666,7 0,0 0,0 106666,7 0,0 20
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 21
ELEMENTO 18 columna
Xi = 12,0 Xf = 12,0 Yi = 6,0 Yf = 9,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
2
36 37 3 28
29
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 2
0,0 106666,7 0,0 0,0 -106666,7 0,0 36
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 37
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 3
0,0 -106666,7 0,0 0,0 106666,7 0,0 28
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 29
91
ELEMENTO 19 columna
Xi = 12,0 Xf = 12,0 Yi = 3,0 Yf = 6,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
1
44 45 2 36
37
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 1
0,0 106666,7 0,0 0,0 -106666,7 0,0 44
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 45
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 2
0,0 -106666,7 0,0 0,0 106666,7 0,0 36
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 37
ELEMENTO 20 columna
Xi = 12,0 Xf = 12,0 Yi = 0,0 Yf = 3,0 Cos α = 0,0 Sen α = 1,0 L = 3,0 B = 0,4 H = 0,4 A = 0,2 E = 2000000,0 I= 0,0 a = 1896,3 b = 0,0 c = 2844,4 f = 106666,7 h = 0,0 k = 5688,9 k/2 = 2844,4
55
56 57 1 44
45
1896,3 0,0 -2844,4 -1896,3 0,0 -2844,4 55
0,0 106666,7 0,0 0,0 -106666,7 0,0 56
-2844,4 0,0 5688,9 2844,4 0,0 2844,4 57
-1896,3 0,0 2844,4 1896,3 0,0 2844,4 1
0,0 -106666,7 0,0 0,0 106666,7 0,0 44
-2844,4 0,0 2844,4 2844,4 0,0 5688,9 45
92
ELEMENTO 21 Viga
Xi = 0,0 Xf = 4,0 Yi = 3,0 Yf = 3,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
1
38 39 1 40
41
60000,0 0,0 0,0 -60000,0 0,0 0,0 1
0,0 337,5 675,0 0,0 -337,5 675,0 38
0,0 675,0 1800,0 0,0 -675,0 900,0 39
-60000,0 0,0 0,0 60000,0 0,0 0,0 1
0,0 -337,5 -675,0 0,0 337,5 -675,0 40
0,0 675,0 900,0 0,0 -675,0 1800,0 41
ELEMENTO 22 Viga Xi = 4,0 Xf = 8,0
Yi = 3,0 Yf = 3,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
1
40 41 1 42
43
60000,0 0,0 0,0 -60000,0 0,0 0,0 1
0,0 337,5 675,0 0,0 -337,5 675,0 40
0,0 675,0 1800,0 0,0 -675,0 900,0 41
-60000,0 0,0 0,0 60000,0 0,0 0,0 1
0,0 -337,5 -675,0 0,0 337,5 -675,0 42
0,0 675,0 900,0 0,0 -675,0 1800,0 43
93
ELEMENTO 23 Viga
Xi = 8,0 Xf = 12,0 Yi = 3,0 Yf = 3,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
1
42 43 1 44
45
60000,0 0,0 0,0 -60000,0 0,0 0,0 1
0,0 337,5 675,0 0,0 -337,5 675,0 42
0,0 675,0 1800,0 0,0 -675,0 900,0 43
-60000,0 0,0 0,0 60000,0 0,0 0,0 1
0,0 -337,5 -675,0 0,0 337,5 -675,0 44
0,0 675,0 900,0 0,0 -675,0 1800,0 45
ELEMENTO 24 Viga Xi = 8,0 Xf = 12,0
Yi = 6,0 Yf = 6,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
2
34 35 2 36
37
60000,0 0,0 0,0 -60000,0 0,0 0,0 2
0,0 337,5 675,0 0,0 -337,5 675,0 34
0,0 675,0 1800,0 0,0 -675,0 900,0 35
-60000,0 0,0 0,0 60000,0 0,0 0,0 2
0,0 -337,5 -675,0 0,0 337,5 -675,0 36
0,0 675,0 900,0 0,0 -675,0 1800,0 37
94
ELEMENTO 25 Viga Xi = 4,0 Xf = 8,0
Yi = 6,0 Yf = 6,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
2
32 33 2 34
35
60000,0 0,0 0,0 -60000,0 0,0 0,0 2
0,0 337,5 675,0 0,0 -337,5 675,0 32
0,0 675,0 1800,0 0,0 -675,0 900,0 33
-60000,0 0,0 0,0 60000,0 0,0 0,0 2
0,0 -337,5 -675,0 0,0 337,5 -675,0 34
0,0 675,0 900,0 0,0 -675,0 1800,0 35
ELEMENTO 26 Viga Xi = 0,0 Xf = 4,0
Yi = 6,0 Yf = 6,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
2
30 31 2 32
33
60000,0 0,0 0,0 -60000,0 0,0 0,0 2
0,0 337,5 675,0 0,0 -337,5 675,0 30
0,0 675,0 1800,0 0,0 -675,0 900,0 31
-60000,0 0,0 0,0 60000,0 0,0 0,0 2
0,0 -337,5 -675,0 0,0 337,5 -675,0 32
0,0 675,0 900,0 0,0 -675,0 1800,0 33
95
ELEMENTO 27 Viga
Xi = 0,0 Xf = 4,0 Yi = 9,0 Yf = 9,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
3
22 23 3 24
25
60000,0 0,0 0,0 -60000,0 0,0 0,0 3
0,0 337,5 675,0 0,0 -337,5 675,0 22
0,0 675,0 1800,0 0,0 -675,0 900,0 23
-60000,0 0,0 0,0 60000,0 0,0 0,0 3
0,0 -337,5 -675,0 0,0 337,5 -675,0 24
0,0 675,0 900,0 0,0 -675,0 1800,0 25
ELEMENTO 28 Viga Xi = 4,0 Xf = 8,0
Yi = 9,0 Yf = 9,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
3
24 25 3 26
27
60000,0 0,0 0,0 -60000,0 0,0 0,0 3
0,0 337,5 675,0 0,0 -337,5 675,0 24
0,0 675,0 1800,0 0,0 -675,0 900,0 25
-60000,0 0,0 0,0 60000,0 0,0 0,0 3
0,0 -337,5 -675,0 0,0 337,5 -675,0 26
0,0 675,0 900,0 0,0 -675,0 1800,0 27
96
ELEMENTO 29 Viga
Xi = 8,0 Xf = 12,0 Yi = 9,0 Yf = 9,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
3
26 27 3 28
29
60000,0 0,0 0,0 -60000,0 0,0 0,0 3
0,0 337,5 675,0 0,0 -337,5 675,0 26
0,0 675,0 1800,0 0,0 -675,0 900,0 27
-60000,0 0,0 0,0 60000,0 0,0 0,0 3
0,0 -337,5 -675,0 0,0 337,5 -675,0 28
0,0 675,0 900,0 0,0 -675,0 1800,0 29
ELEMENTO 30 Viga Xi = 8,0 Xf = 12,0
Yi = 12,0 Yf = 12,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
4
18 19 4 20
21
60000,0 0,0 0,0 -60000,0 0,0 0,0 4
0,0 337,5 675,0 0,0 -337,5 675,0 18
0,0 675,0 1800,0 0,0 -675,0 900,0 19
-60000,0 0,0 0,0 60000,0 0,0 0,0 4
0,0 -337,5 -675,0 0,0 337,5 -675,0 20
0,0 675,0 900,0 0,0 -675,0 1800,0 21
97
ELEMENTO 31 Viga
Xi = 4,0 Xf = 8,0 Yi = 12,0 Yf = 12,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
4
16 17 4 18
19
60000,0 0,0 0,0 -60000,0 0,0 0,0 4
0,0 337,5 675,0 0,0 -337,5 675,0 16
0,0 675,0 1800,0 0,0 -675,0 900,0 17
-60000,0 0,0 0,0 60000,0 0,0 0,0 4
0,0 -337,5 -675,0 0,0 337,5 -675,0 18
0,0 675,0 900,0 0,0 -675,0 1800,0 19
ELEMENTO 32 Viga Xi = 0,0 Xf = 4,0
Yi = 12,0 Yf = 12,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
4
14 15 4 16
17
60000,0 0,0 0,0 -60000,0 0,0 0,0 4
0,0 337,5 675,0 0,0 -337,5 675,0 14
0,0 675,0 1800,0 0,0 -675,0 900,0 15
-60000,0 0,0 0,0 60000,0 0,0 0,0 4
0,0 -337,5 -675,0 0,0 337,5 -675,0 16
0,0 675,0 900,0 0,0 -675,0 1800,0 17
98
ELEMENTO 33 Viga
Xi = 0,0 Xf = 4,0 Yi = 15,0 Yf = 15,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
5
6 7 5 8
9
60000,0 0,0 0,0 -60000,0 0,0 0,0 5
0,0 337,5 675,0 0,0 -337,5 675,0 6
0,0 675,0 1800,0 0,0 -675,0 900,0 7
-60000,0 0,0 0,0 60000,0 0,0 0,0 5
0,0 -337,5 -675,0 0,0 337,5 -675,0 8
0,0 675,0 900,0 0,0 -675,0 1800,0 9
ELEMENTO 34 Viga Xi = 4,0 Xf = 8,0
Yi = 15,0 Yf = 15,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
5
8 9 5 10
11
60000,0 0,0 0,0 -60000,0 0,0 0,0 5
0,0 337,5 675,0 0,0 -337,5 675,0 8
0,0 675,0 1800,0 0,0 -675,0 900,0 9
-60000,0 0,0 0,0 60000,0 0,0 0,0 5
0,0 -337,5 -675,0 0,0 337,5 -675,0 10
0,0 675,0 900,0 0,0 -675,0 1800,0 11
99
ELEMENTO 35 Viga
Xi = 8,0 Xf = 12,0 Yi = 15,0 Yf = 15,0 Cos α = 1,0 Sen α = 0,0 L = 4,0 B = 0,4 H = 0,3 A = 0,1 E = 2000000,0 I= 0,0 a = 60000,0 b = 0,0 c = 0,0 f = 337,5 h = 675,0 k = 1800,0 k/2 = 900,0
5
10 11 5 12
13
60000,0 0,0 0,0 -60000,0 0,0 0,0 5
0,0 337,5 675,0 0,0 -337,5 675,0 10
0,0 675,0 1800,0 0,0 -675,0 900,0 11
-60000,0 0,0 0,0 60000,0 0,0 0,0 5
0,0 -337,5 -675,0 0,0 337,5 -675,0 12
0,0 675,0 900,0 0,0 -675,0 1800,0 13
PASO 4. Calculo de la matriz total de la estructura. K.
Se procede a ensamblar la matriz K total de 70 * 70. No se adicionara por la
extensión de la matriz para que no quede el documento tan extenso. Solo
mostraremos [ Kc ]
122467 -4175 1563 -470 615
-4175 369477 -6900 2068 -470
1563 -6900 370674 -6900 1563
-470 2068 -6900 369477 -4175
615 -470 1563 -4175 362467
100
PASO 5. Calculo de los vectores de cargas y desplazamientos.
VECTOR DE CARGAS NODALES
F1 5 Ton
F2 10 Ton
F3 20 Ton
F4 30 Ton
F5 40 Ton
VECTOR DE DESPLAZAMIENTOS
{ U } = [ Kc ]Inv. * { F }
U1 0,00004 m
U2 0,00003 m
U3 0,00006 m
U4 0,00008 m
U5 0,00011 m
Estos son los desplazamientos por piso en dirección x.
Luego se procede a analizar el pórtico en dirección y de la misma manera que en
dirección x, por lo tanto estos son los desplazamientos por piso en dirección y,
con su respectiva matriz [ Kc ].
480495 -5954 2270 -600 761
-5954 492381 -9950 2558 -563
2270 -9950 490718 -7517 1827
-600 2558 -7517 492025 -5322
761 -563 1827 -5322 483171
101
VECTOR DE DESPLAZAMIENTOS
U1 0,00001 m
U2 0,00002 m
U3 0,00004 m
U4 0,00006 m
U5 0,00008 m
102
18. CONCLUSIONES
Se pudo observar mediante las matrices totales de los ejemplos, que son
la imagen fotográfica de la estructura.
Cuando en una estructura se encuentran acciones de desplazamientos en
los apoyos, la matriz total de la estructura se debe calcular completamente,
caso contrario ocurre cuando en los apoyos los desplazamientos son nulos.
Cuando en una estructura a porticada se desprecian las deformaciones
axiales de las vigas, todos los nudos de un mismo piso sufren el mismo
desplazamiento horizontal, esto fue lo de condensación matricial o también
conocido diafragma rígido.
Trabajar con las matrices en coordeadas globales.
103
19. RECOMENDACIONES
Plasmar la estructura en el primer cuadrante del plano cartesiano con el
fin de que todas las coordenadas sean positivas.
Llevar las unidades de fuerza, longitud y esfuerzos a un solo sistema.
Aplicar cada paso estipulado en los ejercicios de propuestos y
solucionados del manual.
Trabar todos los elementos de la estructura con un mismo sentido y
ubicar estos sentidos en la misma dirección.
104
20. BIBLIOGRAFIA Y CIBERGRAFIA
http://www.uhu.es/javier.pajon/apuntes/matricial.pdf
http://es.wikipedia.org/wiki/P%C3%B3rtico
ROCHEL AWAD. Roberto. Análisis matricial de estructuras.