UNIVERSIDAD NACIONAL DE CAJAMARCA

FACULTAD DE INGENIERÍA
ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA CIVIL

TRABAJO PRÁCTICO DE ANALISIS MATRICIAL DE UN PORTICO.

ASIGNATURA: ANALISIS ESTRUCTURAL I

GRUPO: B1

DOCENTE: Ing. VILCHEZ JARA JOSE MANUEL

ALUMNOS:

MEDINA GAMBOA, Neyder.

MOSQUEIRA CHALAN, Jhony
PAREDES RAMIREZ, Adrian

Cajamarca, 4 de diciembre del 2020 / PRESENTACIÓN VIRTUAL

Índice
Introducción ...................................................................................................................... 3
OBJETIVOS ..................................................................................................................... 4
MARCO TEÓRICO ......................................................................................................... 5
Matriz de rigidez de un elemento de pórtico ................................................................ 5
Matriz de rigidez de un elemento de pórtico ................................................................ 7
VECTOR DE CARGAS Q ............................................................................................. 10
Zona de los pórticos a trabajar .................................................................................... 13
CÁLCULOS DE LA MATRIZ DE RIGIDEZ DE CADA ELEMENTO...................... 18
COEFICIENTES DE LA MATRIZ DE RIGIDEZ GLOBAL DE LA ESTRUCTURA 24
MATRIZ DE RIGIDEZ GLOBAL DE LA ESTRUCTURA ........................................ 25
FUERZAS FICTICIAS POR EMPOTRAMIENTO ...................................................... 26
VECTOR DE CARGAS DE UN PORTICO ................................................................. 31
REFERENCIA ............................................................................................................... 34
Introducción
La formulación matricial del análisis estructural fue desarrollada a fines del siglo XIX
.Sin embargo, no se la desarrollo completamente hasta el avenimiento de las
computadoras capaces de realizar el trabajo numérico al que conduce un análisis con
matrices. El método matricial consiste en asignar a la estructura un objeto matemático
que relacione los desplazamientos de los nudos de la estructura con las fuerzas exteriores.
Si bien, tanto los métodos de flexibilidad como los de rigidez se pueden representar en
forma matricial, el método de rigidez es el más sistemático y adecuado para la
programación. El método matricial considera a la estructura como un montaje de
miembros conectados en sus extremos a nudos.
Con la creciente y amplio uso de la tecnología digital que se ha dado recientemente es
imposible que su influencia no llegue a la rama de la ingeniería civil y, de manera más
específica , al campo de la ingeniería estructural , tanto así que existen en la actualidad
infinidad de herramientas computacionales dirigidas a resolver problemas de análisis y
diseño estructural . El análisis estructural mediante tecnologías digitales es muy recorrido
ya que es de enorme utilidad al realizar cálculos laboriosos y extremadamente repetitivos,
evitando así cometer posibles errores. Sin embargo, la parte importante y crucial de
problemas de este tipo no es realizar las operaciones requeridas por el método de análisis,
sino plantear el problema en cuestión de forma correcta.
OBJETIVOS
 analizar la estructura de un pórtico y aplicar la matriz de rigidez
 Analizar los grados de libertad que tiene el pórtico
 Calcular el vector de cargas del pórtico
MARCO TEÓRICO
Matriz de rigidez de un elemento de pórtico

La figura 4.1 muestra un pórtico típico compuesto por vigas y columnas sometido
a la acción De cargas de gravedad (normalmente distribuidas de manera uniforme)
y sísmicas (normalmente crecientes Hacia arriba). Las uniones entre elementos
permiten la transmisión de fuerzas horizontales, Verticales y momentos
Con base en lo expuesto en el capítulo anterior, deduciremos a continuación la
matriz de unElemento de sección constante caracterizado por un módulo de
elasticidad E, momento de inerciaI y longitud l. Nótese que, con respecto al
elemento usado en el capítulo anterior para una viga, tenemos ahora la presencia
de fuerza axiales Ni y Nj , que se agregan a las fuerzas Vi,Mi, Vj y Mj . La
ecuación matricial correspondiente a las nuevas fuerzas Ni y Nj es:
Mientras que las propias fuerzas restantes es :

Al reunir ambas expresiones en una sola expresión tenemos:

Esto indica que para los vectores de desplazamiento y fuerzas dadas por

Matriz de rigidez de un elemento de pórtico

La matriz de rigidez que los relaciona es:

Donde se ha omitido el subíndice e en las propiedades de los elementos en aras de la

claridad en la notación. La ecuación matricial es, en consecuencia
𝑝𝑒 = 𝑘𝑒 𝑑𝑒
Consideremos ahora la situación general en la que el elemento tiene un ángulo de
inclinación B
Con respecto a la horizontal (figura 4.2). Al tener en cuenta las deducciones de
transformación de fuerzas realizadas en el capítulo 2 y recordando que el vector de
momentos es libre (lo cual implica que Mi = Zi yMj = Zj), se tiene que la matriz de
transformación entre los sistemas de fuerzas local y global mostrados en la figura es

Donde, como de costumbre 𝜂 = cos 𝛽 𝑦 𝜇 = sin 𝛽 , con sus valores dados por:
𝑥𝑗 − 𝑥𝑖
𝜂=
𝑙𝑒
𝑦𝑗 − 𝑦𝑖
𝜂=
𝑙𝑒

Fuerzas internas en un elemento de pórtico. (a) Sistema local; (b) sistema global
En consecuencia, la matriz de rigidez del elemento en coordenadas globales se obtiene
por la operación ya conocida
𝑘𝑒 = 𝑇𝑒𝑇 𝑘𝑒 𝑇𝑒
Los elementos de esta matriz deben ser colocados en una matriz Ke, de tamaño n × n,

donde n es el número de grados de libertad, de acuerdo con el cuadro de correspondencias

entre las numeraciones local y global. La matriz de rigidez del pórtico será entonces
𝑀

𝐾 = ∑ ∆𝐾𝑒
𝑒=1

Donde m es el número de elementos. Esta matriz relaciona las fuerzas y

desplazamientos en un sistema global de coordenadas:
𝐾𝐷 = 𝑃

Una vez resuelto el problema de desplazamiento, las reacciones en los apoyos y las
fuerzas internas de cada elemento (es decir, axiales, cortantes y momentos) se calculan
por medio de las ecuaciones generales deducidas.

Figurara 10. Idealización del modelo estructural en pórtico de tres vanos

Pórticos planos con elementos flexibles

Se inicia el estudio calculando el número de grados de libertad de un pórtico plano
compuesto por elementos lineales que son totalmente flexibles, que no tienen restricción
para deformarse a los cuales se les ha identificado con las letras A , I . La primera letra
hace relación al área de la sección transversal y la segunda al momento de inercia
La configuración del sistema vendrá dada por la posición de los nudos. Por consiguiente,
la definición del número de grados de libertad no es la general enunciada en mecánica,
sino una particular limitada a describir la posición de las juntas. En consecuencia, el
número de grados de libertad, en estática, es el mínimo número de coordenadas que es
preciso determinar para definir la posición de las juntas o nudos

VECTOR DE CARGAS Q
Existen dos maneras de encontrar el Vector de Cargas Q en una estructura. La primera
Es mediante el Algoritmo denominado Problema Primario y Complementario. La
segunda mediante Trabajos Virtuales. En este apartado se va a presentar el cálculo del
vector Q de la primera forma, para el efecto se resuelve un pórtico interior del Colegio
mostrado en la figura 11. En la figura 12 se presenta a la izquierda el modelo numérico
de cálculo y a la derecha el sistema de Coordenadas Generalizadas
Zona de los pórticos a trabajar
Las medidas y propiedades de los materiales de los elementos se detallan a continuación:
Losas:
Dimensiones:
Peralte = 18 cm.
Longitud de losa, de eje a eje de columna = 6 m.
Ancho tributario = 3.6 m.
Columnas:
Dimensiones: Peralte = 30 cm.
Ancho = 30 cm.
Altura de entrepiso = 2.8 m.
Propiedades del material de losas y columnas:
Concreto armado: f´c = 41 MPa (418 kg/cm2).
Peso específico: W = 25000 N/m3 (2549 kg/m3).
Módulo de elasticidad: Ec = 35 GPa (3.56x105 kg/cm2) a los 28 días.
Puntal y repuntal:
Tipo: Ulma PN 1.75 - 3.10.
Altura: 2.8 m
Área de sección transversal = 1.54 cm2.
El número de puntales por piso en el ancho efectivo es de 15 unidades, distribuidas en 3
filas y 5 puntales en cada fila. De acuerdo al método, en el modelo se idealizarán sólo 3
puntales que tendrán una rigidez equivalente a la rigidez total de los puntales ubicados en
cada fila.
Propiedades del material:
Esfuerzo de fluencia = 248.2 MPa. (2532 kg/ cm2).
Esfuerzo último a la tensión = 399.8 MPa. (4077 kg/ cm2).
Módulo de elasticidad, Es = 210 GPa. (2.1x106 kg/cm2).
Peso específico = 76970 N/m3. (7850 kg/m3).
Datos
E = 2.00E+06 tonf/m² elasticidad

F1 = 10.00 tonf fuerza 1

F2 = 20.00 tonf fuerza 2

distrubuida
W1 = 2.00 tonf/m
1

distrubuida
W2 = 1.00 tonf/m
2

altura nivel
H1 = 3.00 m
1

altura nivel
H2 = 3.00 m
2

longitud de
L= 4.00 m
vano
 eje A
Columnas
b (m) t (m)

nivel 1 (AB) 0.30 0.40

nivel 2 (BC) 0.30 0.40

nivel 3(CD= 0.30 0.30

eje B
Columnas
b (m) t (m)

nivel 1 (DE) 0.30 0.40

nivel 2 (EF) 0.30 0.40

nivel 3(CD= 0.30 0.30

Vigas b (m) h (m)

nivel 1 (BE) 0.30 0.50

nivel 2 (CF) 0.30 0.50

nivel 3 (DG) 0.3 0.5

CÁLCULOS DE LA MATRIZ DE RIGIDEZ DE CADA ELEMENTO

 Para las columnas no se tomara en cuenta las deformaciones axiales , es por esto
que no se mustra grados de livertad verticales

 Para las vigas no se tomara en cuenta las deformaciones axiales ni cortantes . en

la resolución solo se usaran las deformaciones angulares ( solo para la matriz de
rigidez global )

Columna AB

E = 2.00 tonf/m²

I = 0.0016 m⁴

H = 3.00m

0 0 1 3
 GDL GDL
1 2 3 4
local global

1422 -2133 -1422 -2133 1 0

-2133 4267 2133 2133 2 0

KAB=
-1422 2133 1422 2133 3 1

-2133 2133 2133 4267 4 3

Coeficientes usados en la matriz de rigidez global.

k11 = 1422

k13 = 2133

k33 = 4267

Columna BC

E = 2.00E tonf/m²

I = 0.0016 m⁴

H = 3.00m
 1 3 2 5 GDL GDL

1 2 3 4 local global

1422 -2133 -1422 -2133 1 1

-2133 4267 2133 2133 2 3

KBC=
-1422 2133 1422 2133 3 2

-2133 2133 2133 4267 4 5

Coeficientes usados en la matriz de rigidez global.

k11 = 1422 k22 = 1422 k35 = 2133

k12 = -1422 k23 = 2133 k55 = 4267

k13 = -2133 k25 = 2133

k15 = -2133 k33 = 4267

Columna DE

E = 2.00E tonf/m²

I = 0.0016 m⁴

H = 3.00m
 0 0 1 4 GDL GDL

1 2 3 4 local global

1422 -2133 -1422 -2133 1 0

-2133 4267 2133 2133 2 0

KDE=
-1422 2133 1422 2133 3 1

-2133 2133 2133 4267 4 4

k11 = 1422

k14 = 2133

k44 = 4267

Columna EF

E= 2.00E+06 tonf/m²

I= 0.0016 m⁴

H= 3.00 m
 1 4 2 6 GDL GDL

1 2 3 4 local global

1422 -2133 -1422 -2133 1 1

-2133 4267 2133 2133 2 4

KEF=
-1422 2133 1422 2133 3 2

-2133 2133 2133 4267 4 6

k11 = 1422 k22 = 1422 k46 = 2133

k12 = -1422 k24 = 2133 k66 = 4267

k14 = -2133 k26 = 2133

k16 = -2133 k44 = 4267

Viga BE

E= 2.00E+06 tonf/m²

I= 0.003125 m⁴

L= 4.00 m

0 3 0 4 GDL GDL

1 2 3 4 local global

1172 2344 -1172 2344 1 0

2344 6250 -2344 3125 2 3

KBE=
-1172 -2344 1172 -2344 3 0

2344 3125 -2344 6250 4 4

 Coeficientes usados en la matriz de rigidez global.

k33 = 6250

k34 = 3125

k44 = 6250

Viga CF

E= 2.00E+06 tonf/m²

I= 0.003125 m⁴

L= 4.00 m

0 5 0 6 GDL GDL

1 2 3 4 local global

1172 2344 -1172 2344 1 0

2344 6250 -2344 3125 2 5

KCF=
-1172 -2344 1172 -2344 3 0

2344 3125 -2344 6250 4 6

Coeficientes usados en la matriz de rigidez global

k55 =6250

k56 =3125

k66 =3125
COEFICIENTES DE LA MATRIZ DE RIGIDEZ GLOBAL DE LA
ESTRUCTURA

k11 = 5689 k22 = 2844 k33 = 14783 k44 = 14783 k55 = 10517 k66 = 10517
k12 = -2844 k23 = 2133 k34 = 3125 k45 = 0 k56 = 3125
k13 = 0 k24 = 2133 k35 = 2133 k46 = 2133
k14 = 0 k25 = 2133 k36 = 0 coeficiente
k15 = -2133 k26 = 2133 coeficientes para la fila 6
k16 = -2133 coeficientes para la fila 5

coeficientes para la fila 4

coeficientes para la fila 3

coeficientes para la fila 2

para la fila 1
MATRIZ DE RIGIDEZ GLOBAL DE LA ESTRUCTURA

k11 k12 k13 k14 k15 k16

k22 k23 k24 k25 k26

k33 k34 k35 k36

K=
k44 k45 k46
simétrica
k55 k56

k66

Nota: La matriz de rigidez es una matriz simétrica, es decir:

k12=k21, k14=k41, k35=k53, etc.

1 2 3 4 5 6

5689 -2844 0 0 -2133 -2133 1

-2844 2844 2133 2133 2133 2133 2

0 2133 14783 3125 2133 0 3

K=
0 2133 3125 14783 0 2133 4

-2133 2133 2133 0 10517 3125 5

-2133 2133 0 2133 3125 10517 6

Grados de libertad

globales (GDL Globales)

FUERZAS FICTICIAS POR EMPOTRAMIENTO

Fuerzas ficticias para la viga BE:

0.00 1

0.00 2

2.67 3
Fo(AB) =
-2.67 4

0.00 5

0.00 6

W1 = 2.00 tonf/m
 L= 4.00 m

RB = 4.00 tonf

MB = 2.67 tonf-m

RE = 4.00 tonf

ME = 2.67 tonf-m

Fuerzas ficticias para la viga CF:

0.00 1

0.00 2

0.00 3
Fo(CF) =
0.00 4

1.33 5

-1.33 6

W2 = 1.00 tonf/m

L= 4.00 m

RC = 2.00 tonf
 MC = 1.33 tonf-m

RF = 2.00 tonf

MF = 1.33 tonf-m

Fuerzas puntuales en cada grado de libertad global

10.00 1

20.00 2

0.00 3
Fp =
0.00 4

0.00 5

0.00 6

Vector de fuerzas F=Fp-∑Fo

 10.00 0.00 0.00

20.00 0.00 0.00

0.00 2.67 0.00

F= - +
0.00 -2.67 0.00

0.00 0.00 1.33

0.00 0.00 -1.33

vector de fuerzas fuerzas

fuerzas ficticias ficticias

puntuales en BE en CF

10.00 1

20.00 2

-2.67 3
F=
2.67 4

-1.33 5

1.33 6
 Al no haber desplazamientos y giros en el Problema Primario, cada uno de los
elementos se considera empotrado-empotrado. Por este motivo en la figura 14 se presenta
las acciones de empotramiento perfecto que viene a ser las fuerzas y momentos con las
cuales se equilibran los elementos.
VECTOR DE CARGAS DE UN PORTICO

5689 -2844 0 0 -2133 -2133 u1 10.00

-2844 2844 2133 2133 2133 2133 u2 20.00

0 2133 14783 3125 2133 0 x u3 = -2.67

0 2133 3125 14783 0 2133 u4 2.67

-2133 2133 2133 0 10517 3125 u5 -1.33

-2133 2133 0 2133 3125 10517 u6 1.33

Matriz de rigidez global 6x6 vector de desplazamiento vector fuerzas

0.01568 1 (m) desplazamiento

0.03053 2 (m) desplazamiento

-0.00363 3 (rad) rotación

U=
-0.00322 4 (rad) rotación

-0.00191 5 (rad) rotación

-0.00167 6 (rad) rotación

Ensamblaje de Q
En cada Junta se debe verificar que la suma de fuerzas en sentido horizontal y vertical sea
cero; lo propio debe suceder con la sumatoria de momentos. En esta etapa se encuentran
los valores de las fuerzas de fijación con su respectivo signo. En el Problema
Complementario actúan únicamente las cargas de fijación pero con sentido contrario al
del Problema Primario y son las que generan los desplazamientos y giros en la estructura.

Por último es importante indicar que si se desea resolver un pórtico plano, para diferentes
estados de carga, se deberá obtener en primer lugar la matriz de rigidez de la estructura y
luego para cada uno de los estados se encuentra el vector de cargas generalizadas y se
hallan las fuerzas y momentos finales. Se hace incapié que solo se encuentra una sola vez
la matriz de rigidez
REFERENCIA
1. Aguiar Roberto (2004), Análisis Matricial de Estructuras, Centro de Investigaciones
Científicas, CEINCI. Escuela Politécnica del Ejército. Tercera Edición, 550 p. Quito,
Ecuador.