CAPITULO 15

ANALISIS ESTATICO CON CEINCI-LAB

RESUMEN

Se presenta la teoría y el uso del sistema de computación CEINCI-LAB para resolver:

armaduras planas y pórticos planos. Una buena parte de los programas que se utilizan para
resolver Armaduras Planas sirven para Pórticos Planos y Mallas Espaciales, por esta razón
cuando se estudie los pórticos planos se indicará con detalle el uso de estos programas, lo
propio se realizará en el Capítulo 18 dedicado a las Mallas.

Se presentan dos ejemplos para el caso de Armaduras Planas, el segundo ejemplo es

bastante práctico y tiene que ver con el reforzamiento de una estructura de acero de una
Piscina; estas estructuras son bastante rígidas en el sentido transversal porque se debe cubrir
una gran luz, pero en el sentido longitudinal son bastante flexibles por lo que se recomienda
reforzar con una celosía Tipo San Andrés. Se resuelve un ejemplo de ésta celosía pero
considerando únicamente tres elementos horizontales para no alargar la exposición y poder
resolver completamente en pocas páginas.

En el primer ejercicio que se resuelve con Pórticos Planos, se explica con mucho
detenimiento el uso de los diferentes programas que se requieren en su solución pero a más de
ello se presenta el marco teórico respectivos, para repasar lo estudiado en capítulos anteriores.

Luego se resuelven dos pórticos de un piso y un vano, en el primer caso se ve que no se

obtienen directamente las reacciones de los apoyos, en cambio en el segundo caso se
encuentra directamente las reacciones en los apoyos. Nótese que para el primer caso se tienen
6 grados de libertad mientras que para el segundo caso se tienen 12 grados de libertad ya que
se consideró tres grados de libertad en cada uno de los apoyos. En el primer caso la matriz de
rigidez admite inversa en cambio en el segundo caso la matriz de rigidez no admite inversa por
lo que se trabaja con la matriz de rigidez condensada.

Posteriormente se presenta la solución de un pórtico con cargas puntuales sobre los

elementos, esto con el objeto de ilustrar cómo se resuelve un pórtico con cualquier tipo de
carga. En este caso el usuario debe encontrar las acciones de empotramiento perfecto en
coordenadas locales y dar como dato al programa.

Una forma de reforzar una estructura de hormigón armado es mediante la colocación de

diagonales de acero, con lo que se convierte en una estructura mixta de hormigón y acero. La
solución de este tipo de estructuras se ilustra mediante la solución de un pórtico con vigas y
columnas de hormigón y diagonales de acero.
542 Roberto Aguiar Falconí
UFA - ESPE

Las estructuras de acero, si tienen conexiones simples, se calcula con dos grados de
libertad pero si tiene nudos “Precalificados” con capacidad de momento, se deben analizar con
tres grados de libertad, similar a un pórtico plano de hormigón armado. Para ilustrar lo expuesto
se presenta una conexión precalificada BFP, utilizada en la construcción del edificio de la
Unidad de Vigilancia Comunitaria, UVC, de la ciudad de Manta, que soportó el terremoto del 16
de abril de 2016, de magnitud 7.8; pero a más de ello interesa ilustrar la forma de calcular
estructuras de acero con CEINCI-LAB.

15.1 ARMADURAS PLANAS

Son estructuras que tienen dos grados de libertad por nudo, un desplazamiento
horizontal y un desplazamiento vertical; existen giros en una armadura plana pero estas son
coordenadas dependientes. Sus elementos solo trabajan axialmente, no trabajan a flexión
(Conexión Simple) por este motivo es que se los modela como A0 , I = ∞ .

15.1.1 Matriz de rigidez de un elemento

En la parte superior izquierda de la figura 15.1 se presenta el sistema de coordenadas

del elemento cuando no se consideran los desplazamientos como cuerpo rígido, de ahí la
presencia de los apoyos, a este sistema se denomina 1. En la parte superior izquierda se tiene
el sistema de coordenadas locales que se identifica con el 2, aquí ya se consideran los
desplazamientos como cuerpo rígido. Finalmente en la parte inferior se muestra el sistema 3
que corresponde al sistema de coordenadas globales de un elemento de una armadura plana.

Figura 15.1 Sistema de coordenadas de un elemento de una armadura plana.

La matriz de rigidez para el sistema uno se denomina k1 y es:

 EA 
k1 =  (15.1)
 L 

Donde E es el módulo de elasticidad del material; A es el área de la sección

transversal y L la longitud del elemento. La matriz de paso del sistema de coordenadas 1 al
sistema de coordenadas 2 es:

T1− 2 = [− 1 0 1 0] (15.2)
 ANALISIS ESTATICO CON CEINCI-LAB 543

La matriz de rigidez en el sistema 2, que se llama k2 , se halla del triple producto

matricial, siguiente.

k 2 = T1− 2 k1 T1− 2
t
(15.3)

k2 Es la matriz de rigidez en coordenadas locales (axial y transversal). Finalmente la

matriz de paso del sistema 2 al sistema 3, es:

 cos α senα 0 0 
− senα cos α 0 0 

T2−3 = 0 0 cos α senα  (15.4)
 
 0 0 − senα cos α 
 

Siendo α el ángulo que forma el eje del elemento con el eje de las X; así para una
barra horizontal α = 0 . Finalmente la matriz de rigidez en coordenadas globales (horizontal y
vertical) se halla con la siguiente expresión.

k 3 = T2−3 k 2 T2 −3
t
(15.5)

La subrutina que obtiene la matriz de rigidez de una armadura plana en la forma

indicada, en el programa CEINCI-LAB se denomina kdiagonal Los datos de entrada y forma
de uso de este programa, son:

• [k3]=kdiagonal (A,L,E,seno,coseno)

Donde A es el área de la sección transversal; L es la longitud del elemento; E es el

módulo de elasticidad del material; seno es el valor de 𝑠𝑠𝑠𝑠𝑠𝑠 𝛼𝛼; coseno es el valor de cos 𝛼𝛼.

15.1.2 Solución total de una armadura plana

Se presenta un resumen de los pasos a seguir para su solución, indicando los

programas del sistema de computación CEINCI-LAB que se utilizan.

i. Se seleccionan los grados de libertad de la estructura. Esto se lo hace con la subrutina

cg. Se destaca, una vez más, que hay dos grados de libertad por nudo.

ii. Se obtienen los vectores de colocación de cada uno de los elementos. El vector ce
colocación de un elemento está compuesto por los grados de libertad del nudo inicial y
del nudo final. La subrutina que obtiene los vectores de colocación de todos los
elementos se llama vc.

iii. Se halla la matriz de rigidez de la estructura K por ensamblaje directo, con la subrutina
krigidez . Esta subrutina a su vez llama a la subrutina kdiagonal para calcular la
matriz de rigidez de cada elemento.

iv. Se encuentra el vector de cargas generalizadas Q con la subrutina cargas_armadura.

Una armadura plana solo tiene cargas en las juntas de tal manera que el usuario puede
indicar directamente el vector de cargas Q .
544 Roberto Aguiar Falconí
UFA - ESPE

v. Se obtienen los desplazamientos de los nudos q , resolviendo el siguiente sistema de

ecuaciones lineales.
Q=K q (15.6)

vi. Se encuentran las fuerzas finales en cada uno de los elementos, para ello se determina
en primer lugar las deformaciones p en coordenadas globales, con la ayuda del vector
de colocación; luego las fuerzas 𝑷𝑷𝑮𝑮 en coordenadas globales y finalmente con la matriz
de paso T 2−3 las fuerzas 𝑷𝑷𝐿𝐿 en coordenadas locales. Todo esto se hace con la
subrutina fuerzas_armadura

𝑷𝑷𝑮𝑮 = 𝒌𝒌𝟑𝟑 𝒑𝒑 (15.7)

𝑷𝑷𝑳𝑳 = 𝑻𝑻𝟐𝟐−𝟑𝟑 𝑷𝑷𝑮𝑮 (15.8)

vii. Es fundamental que se halle el momento de volcamiento debido a que en algunas
estructuras de acero se han volteado durante un sismo o un fuerte viento, como se
aprecia en la figura 15.2

Figura 15.2 Volcamiento de una torre de transmisión. (Cortesía de Theofilos Toulkeridis)

• EJEMPLO 1

Resolver la armadura plana, con nudos simples, indicada en la figura 15.3, la longitud
L = 50 cm. el área de la sección transversal de los elementos es igual y vale 1.13 cm2. La
fuerza F aplicada en la parte superior vale 50 kg., y el módulo de elasticidad del acero
E = 2100000 kg / cm 2 .

• SOLUCIÓN

En la figura 15.4 se presenta, a la izquierda, la numeración de los nudos y elementos

de la armadura y a la derecha los grados de libertad de acuerdo a la numeración de los nudos.
Antes de resolver el problema, el lector debe leer la entrada de datos y salida de resultados de
algunos de los programas de CEINCI-LAB, que se indican más adelante, en el desarrollo del
 ANALISIS ESTATICO CON CEINCI-LAB 545

Ejemplo 3. Con esta acotación el programa que se ha elaborado para resolver el ejercicio es el
siguiente:

Figura 15.3 Armadura plana a resolver con CEINCI-LAB

Figura 15.4 Numeración de nudos, elementos y grados de libertad

function armadura_ej1
%
% Programa para resolver Armaduras Planas
%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
% Noviembre de 2009
%-------------------------------------------------------------
% armadura_ej1
%-------------------------------------------------------------
546 Roberto Aguiar Falconí
UFA - ESPE

nod=5; % Número de nudos de armadura plana

nr=3; % Número de nudos restringidos
L=50; %Dimensión
%
% MATRIZ DE RESTRICCIONES DE LOS NUDOS
RES=[ 1 1 1; % Nudo restringido Des X, Des Y
2 1 1;
3 1 1];
%
% Calculo de Coordenadas Generalizadas y grados de libertad
[CG,ngl]=cg(nod,nr,RES); % CG Matriz de coord generalizadas
% MATRIZ DE GENERACION DE NUDOS DE LOS ELEMENTOS
GEN=[1 1 2 1 1 1 1; %elem, ni, nf, elem a gene, inc en elem, inc en ni, inc en nf
3 1 4 0 0 0 0;
4 2 4 1 3 2 1;% Se genera elemento 7
5 3 4 0 0 0 0;
6 1 5 1 2 2 0];% Se genera elemento 8
% Generación de vector NI (Nudo Inicial) y NJ (Nudo Final) de elementos
[NI,NJ]=gn_portico(GEN);
% MATRIZ DE GENERACION DE COORDENADAS DE LOS NUDOS
NUDOS=[1 -L 0.0 0 0 0.0 0.0;% i, xi, yi, nudos a gener, incr num nudo, dx,dy
2 0.0 0.0 0 0 0.0 0.0;
3 L 0.0 0 0 0.0 0.0;
4 0 L 0 0 0.0 0.0;
5 0 2*L 0 0 0 0];
% Generación de las coordenadas de los nudos
[X,Y]=glinea_portico(NUDOS);
% Dibuja la armadura plana
dibujo (X,Y,NI,NJ);
% Vector de colocación
[VC]=vc(NI,NJ,CG);
%Longitud, Seno y Coseno de elementos
[L,seno,coseno]=longitud (X,Y,NI,NJ);
% Vector de cargas generalizadas
njc=1; % Número de juntas Cargadas
F=[5 50.0 0.0]; %Junta Cargada, FH, FV
[Q]=cargas_armadura(njc,ngl,CG,F)
% Areas de los elementos
[ELEM]=[1.12; %Area del elemento 1 en cm2
1.12; 1.12;1.12;1.12;1.12;1.12;1.12];
% Matriz de rigidez de la estructura
E=2100000; % Modulo de elasticidad del acero en kg/cm2
[K]=krigidez(ngl,ELEM,L,seno,coseno,VC,E)
% Desplazamientos y Giros
q=K\Q
% Fuerzas y momentos finales en los elementos
[FF]=fuerzas_armadura(ngl,ELEM,L,seno,coseno,VC,E,q)
% Cálculo de Momento de Volcamiento
punto=1; % Es el número de nudo con respecto al cual se calcula momento
[Momento]=mom_vol(ngl,ELEM,L,seno,coseno,VC,E,q,X,Y,NI,NJ,njc,F,punto)
% ---end---

El vector de colocación de un elemento está compuesto por los grados de libertad del
nudo inicial y final. La matriz que contiene a estos vectores se ha denominado VC, que se
indica a continuación.
 ANALISIS ESTATICO CON CEINCI-LAB 547

La primera fila de VC corresponde al vector de colocación del elemento 1, ver figura

izquierda de 15.4; la segunda fila al elemento 2, etc. Los elementos 1 y 2 no aportan a la matriz
de rigidez de la estructura ya que sus vectores de colocación son cero; los restantes elementos
si aportan a la matriz 𝑲𝑲

El programa que determina la matriz de rigidez de la armadura se denomina krigidez,

este programa llama al programa kdiagonal. A continuación se indica, a la izquierda, la matriz
de rigidez del elemento 3 en coordenadas globales; a la derecha se presenta la misma matriz
pero en la parte superior se ha anotado el vector de colocación, que indica la posición a la cual
van cada uno de los elementos, así por ejemplo el término de la primera fila y primera columna
de la matriz de rigidez va a la fila 0 y columna 0 pero como no existe esa fila simplemente no
contribuye a la matriz de rigidez, por este motivo se ha tachado las filas y columnas cuyo vector
de colocación es cero.

El tercer elemento si contribuye y va a la posición (1,1) de la matriz de rigidez de la

estructura que también es de 4 por 4, ya que se tienen 4 grados de libertad. A continuación se
presenta la matriz de rigidez del elemento 4 y su respectiva contribución a la matriz de rigidez
de la estructura, todo esto con el objeto de que el lector comprenda el ensamblaje directo.

La matriz de rigidez el elemento 5, tiene los mismos valores que la del elemento 3 pero
difieren en algunos signos por el ángulo 𝛼𝛼. La contribución del elemento 5 a la matriz de rigidez
se indica a continuación.

En forma similar se procede con los elementos 6, 7 y 8. La matriz de rigidez de la

estructura que se encuentra es la siguiente.
548 Roberto Aguiar Falconí
UFA - ESPE

El vector de cargas 𝑸𝑸 se pudo escribir directamente en el programa sin necesidad de

utilizar ningún programa ya que las carga solo actúa en el nudo 5. Este vector es:

0
𝑸𝑸 = � 0 �
50
0

De la solución del sistema de ecuaciones lineales, se halla el vector 𝒒𝒒 que contiene los
desplazamientos de los nudos 4 y 5.

𝑲𝑲 𝒒𝒒 = 𝑸𝑸

𝟎𝟎
𝒒𝒒 = � 𝟎𝟎 �
𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
𝟎𝟎

El desplazamiento horizontal del grado de libertad 3, es 5.9 mm. Las deformaciones de

los elementos 𝒑𝒑, se halla en base al vector 𝒒𝒒 y al vector de colocación; para los elementos 1 al
5 estas deformaciones son cero; para los elementos 6, 7 y 8 el vector 𝒑𝒑 es el mismo y es el
siguiente.

0
𝑝𝑝(6) = 𝑝𝑝(7) = 𝑝𝑝(8) = � 0 �
0.0059
0

La matriz de rigidez del elemento 7 es la misma que la del elemento 4, que fue ya se
indicó anteriormente, al multiplicar esta matriz por 𝒑𝒑(7) se halla que 𝑷𝑷𝐺𝐺 = 𝟎𝟎. Es cero vector; es
decir todos sus elementos son cero.

Las fuerzas en coordenadas globales para el elemento 6, se halla de la siguiente

manera:
(6) (6)
𝑃𝑃𝐺𝐺 = 𝑘𝑘3 𝑝𝑝(6)

De igual manera se procede con el elemento 8. En la figura 15.5 se indican las fuerzas
en coordenadas globales. Nótese que existe equilibrio en cada uno de los elementos y que
también existe equilibrio de juntas.
 ANALISIS ESTATICO CON CEINCI-LAB 549

Figura 15.5 Fuerzas en Coordenadas Globales en los elementos de armadura plana.

Aquí podría finalizar la solución de la armadura plana, para luego manualmente

encontrar la resultante de las dos fuerzas que actúan en cada nudo, descomponiendo las
fuerzas horizontales y verticales. En forma matricial esto se realiza por medio de la matriz de
paso de coordenadas locales a globales (Ecuación 15.4 ). Para el elemento (6) se indica la
forma de obtener el vector de cargas en coordenadas locales 𝑃𝑃𝐿𝐿
(6) (6) (6)
𝑷𝑷𝐿𝐿 = 𝑻𝑻2−3 𝑷𝑷𝐺𝐺

Los resultados que reporta el programa, a nivel de fuerzas en los elementos se indica
en la figura 15.6. El elemento 6 trabaja a tracción y el elemento 8 a compresión.

Figura 15.6 Fuerzas en Coordenadas Locales de la Armadura.

550 Roberto Aguiar Falconí
UFA - ESPE

Al multiplicar las fuerzas horizontales y verticales que aparecen en la figura 15.5, por su
distancia con respecto a la base, se obtiene el momento de volcamiento, con respecto al nudo
1, que vale M = 5000 kg * cm .

15.1.3 Reforzamiento de estructuras de acero

Las estructuras de acero que son utilizadas para cubrir Piscinas, Naves Industriales o
Coliseos, por lo regular se las construye pensando solo en la carga vertical, no en la sísmica,
es así que son muy rígidas en el sentido transversal, que tienen que cubrir grandes luces,
debido a que las columnas se las orienta de tal manera que tengan un gran momento de
inercia en ese sentido y las vigas son de dimensiones considerables. A más de ello hay una
gran cantidad de pórticos, todo esto hace que la estructura sea muy rígida en el sentido
transversal.

El problema se presenta en el sentido longitudinal, debido en primer lugar al hecho de

que solo se tienen dos pórticos en dicha dirección, que son los exteriores; luego a que no
existe una verdadera viga entre las columnas de éstos pórticos, lo que existe es una correa que
está pensada para soportar el recubrimiento de la Piscina pero en el caso de un sismo debe
trabajar como viga. Finalmente en muy pocos proyectos se observa la presencia de tensores
que son muy adecuados para soportar las fuerzas sísmicas en sentido longitudinal.

En la figura 15.7 se presenta, parte de la cubierta de la Piscina del Colegio Fernández

Madrid, ubicada en el Centro Colonial de Quito, al fondo se observa el pórtico transversal
exterior con mampostería confinada en la parte interior y ventanas en la parte superior. Nótese
las dimensiones considerables que tienen las vigas transversales en contraste con las correas
que existen en sentido longitudinal.

Figura 15.7 Vista interior de la Piscina del Colegio Fernández Madrid de Quito.
Sin reforzamiento.

En la figura 15.8 se presenta el reforzamiento sísmico para la estructura de la Piscina

del Colegio Fernández Madrid, que consiste en encamisar con hormigón las columnas
existentes que confinan la mampostería en los pórticos exteriores transversales y prolongarles
hasta la cercha de la cubierta, a más de ello se colocarán tensores en los vanos exteriores
como se indica en la figura 15.8. Todo esto con el propósito de que la estructura resista un
sismo en sentido longitudinal.
 ANALISIS ESTATICO CON CEINCI-LAB 551

Figura 15.8 Reforzamiento sísmico de la estructura de la Cubierta del

Colegio Fernández Madrid. Aguiar (2012)

Los tensores, en los vanos exteriores, ante la acción sísmica en sentido longitudinal o
ante fuerzas de viento, trabajan como vigas de gran peralte sujetas a fuerzas horizontales. Con
este antecedente se presenta el siguiente ejemplo, que por motivos didácticos a la viga de
borde se la ha hecho muy pequeña para poder presentar los resultados, en forma fácil

• EJEMPLO 2

Resolver la estructura de acero, presentada a la izquierda de la figura 15.9, si está

compuesta por dos perfiles “G” de 100/50/15/3 mm. Las cargas verticales son las debidas a la
acción sísmica en sentido longitudinal. La altura de 3.74 m., de la cercha vendría a ser la
distancia entre los pórticos transversales. En realidad solo se coloca como correa un perfil “G”
de 100/50/15/3 mm., pero se refuerza con otro perfil de las mismas características. El módulo
de elasticidad del acero, es 𝐸𝐸𝑠𝑠 = 2100000 𝑘𝑘𝑘𝑘/𝑐𝑐𝑐𝑐2 .

Figura 15.9 Geometría de Cercha de Ejemplo 2.

552 Roberto Aguiar Falconí
UFA - ESPE

• SOLUCIÓN

Todos los elementos estructurales de la armadura tienen el mismo perfil, luego tienen
la misma sección transversal y es de 11.41 cm2. En realidad no se dibujo toda la viga de borde
en sentido transversal solo se dibujaron tres correas que están separadas 1.14 m.

Figura 15.10 Numeración de nudos, elementos y grados de libertad.

A la izquierda de la figura 15.10 se presenta la numeración de los nudos y elementos

de la armadura plana; a la derecha los grados de libertad considerados. Cada nudo de la
armadura tiene dos grados de libertad, la componente de desplazamiento horizontal y la de
desplazamiento vertical. Nótese que en el nudo 4, existen grados de libertad, ya que se
soldarán los perfiles que llegan a ese nudo.

function armadura_ej2
% Programa para resolver Armaduras Planas tipo San Andres
%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
% Diciembre de 2009
%-------------------------------------------------------------
% armadura_ej2
%-------------------------------------------------------------
nod=7; % Número de nudos de armadura plana
nr=2; % Número de nudos restringidos
F1=137; F2=274;% Se trabaja en kg y cm
%
% MATRIZ DE RESTRICCIONES DE LOS NUDOS
RES=[ 1 1 1; % Nudo restringido Des X, Des Y
2 1 1];
%
% Calculo de Coordenadas Generalizadas y grados de libertad
[CG,ngl]=cg(nod,nr,RES) % CG Matriz de coord generalizadas
% MATRIZ DE GENERACION DE NUDOS DE LOS ELEMENTOS
GEN=[1 1 3 0 0 0 0; %elem, ni, nf, elem a gene, inc en elem, inc en ni, inc en nf
2 3 2 0 0 0 0;
3 5 6 1 1 1 1;
5 1 5 1 3 1 2;
6 3 4 1 1 1 2;
9 1 4 1 3 3 3;
10 2 4 1 1 2 1];
% Generación de vector NI (Nudo Inicial) y NJ (Nudo Final) de elementos
[NI,NJ]=gn_portico(GEN);
 ANALISIS ESTATICO CON CEINCI-LAB 553

% MATRIZ DE GENERACION DE COORDENADAS DE LOS NUDOS

NUDOS=[1 0.0 0.0 1 4 0.0 374;% i, xi, yi, nudos a gener, incr num nudo, dx,dy
3 114 0.0 1 1 0.0 374/2;
2 114*2 0 1 5 0.0 374;
6 114 374 0 0 0.0 0.0];
% Generación de las coordenadas de los nudos
[X,Y]=glinea_portico(NUDOS);
% Dibuja la armadura plana
dibujo (X,Y,NI,NJ);
% Vector de colocación
[VC]=vc(NI,NJ,CG);
%Lngitud, Seno y Coseno de elementos
[L,seno,coseno]=longitud (X,Y,NI,NJ);
% Vector de cargas generalizadas
njc=3; % Número de juntas Cargadas
F=[5 0 -F1;
6 0 -F2;
7 0 -F1]; %Junta Cargada, FH, FV
[Q]=cargas_armadura(njc,ngl,CG,F);
% Areas de los elementos
area=11.41; %area de un elemento
SECCION=[1 area,11 1];
[ELEM]=gelem_armadura(SECCION);
% Matriz de rigidez de la estructura
E=2100000; % Modulo de elasticidad del acero en kg/cm2
[K]=krigidez(ngl,ELEM,L,seno,coseno,VC,E)
% Desplazamientos y Giros
q=K\Q
% Fuerzas y momentos finales en los elementos
[FF]=fuerzas_armadura(ngl,ELEM,L,seno,coseno,VC,E,q)
% Cálculo de Momento de Volcamiento
punto=1; % Es el número de nudo con respecto al cual se calcula momento
[Momento]=mom_vol(ngl,ELEM,L,seno,coseno,VC,E,q,X,Y,NI,NJ,njc,F,punto)
% ---end---

En la figura 15.11 se presenta la solución de la armadura. Los tensores en las cubiertas

de las estructuras de acero, son una excelente alternativa para soportar las fuerzas sísmicas o
de viento en sentido longitudinal.
554 Roberto Aguiar Falconí
UFA - ESPE

6.04 6.04 6.04 6.04

11.60 11.60
127.08 274 127.08

11.60 11.60
274

172.06 172.06

127.08 172.06 172.06 127.08

0
0 0 0 0
Figura 15.11 Fuerzas que actúan en la Armadura del Ejemplo 2.

15.2 PORTICOS PLANOS

Los pórticos planos tienen tres grados de libertad por nudo, un desplazamiento
horizontal, un desplazamiento vertical y una rotación. En este apartado se considera que todos
los elementos de un pórtico, son totalmente flexibles.

15.2.1 Matriz de rigidez de un elemento

En forma similar al caso de una armadura plana, en la figura 15.12 se muestran los tres
sistemas de un elemento de un pórtico plano. El sistema 1 corresponde al caso en que no se
consideran los desplazamientos como cuerpo rígido; el sistema 2 son las coordenadas locales
(axial y transversal) y el sistema 3 son las coordenadas globales (horizontal y vertical).

Figura 15.12 Sistema de coordenadas de un elemento de un pórtico plano.

 ANALISIS ESTATICO CON CEINCI-LAB 555

La matriz de rigidez en el sistema 1, es la siguiente:

 4 EI (1 + φ ) 2 EI (1 − 2φ ) 0

 L (1 + 4φ ) L (1 + 4φ )
 
 2 EI
k1 = 
(1 − 2φ ) 4 EI (1 + φ )
0
 (15.9)


L (1 + 4φ ) L (1 + 4φ ) 
 EA 
 0 0
 L 

Donde E es el módulo de elasticidad del material; I es el momento de inercia de la

sección transversal; A es el área de la sección transversal; L es la longitud del elemento; φ
es el factor que mide el efecto del corte.
3βEI
φ= (15.10)
GAL2

Siendo β el factor de forma que vale 1.2 para elementos rectangulares; G es el

módulo de corte.

La matriz de paso de los sistemas 1 al 2 y del 2 al 3, son:

 1 1 
 0 1 0 − 0
L L
 
= 0 1
1 1
T1−2 0 0 − (15.11)
 L L 
− 1 0 0 1 0 0 

 

 cos α senα 0 0 0 0
− senα cos α 0 0 0 0

 0 0 1 0 0 0
T2−3 =  (15.12)
 0 0 0 cos α senα 0
 0 0 0 − senα cos α 0
 
 0 0 0 0 0 1

Las matrices de rigidez en coordenadas locales k2 y en coordenadas globales k3 se

hallan con las ecuaciones (15.3) y (15.5) ya indicadas. Pero realizando los productos
matriciales con las matrices (15.9), (15.11) y (15.11). En CEINCI-LAB la matriz que determina
la matriz de rigidez de un elemento de un pórtico plano se denomina: kmiembro. Los datos de
entrada y forma de uso es la siguiente:

• [k3]=kmiembro (b,h,L,E,seno,coseno)

Donde 𝒌𝒌𝟑𝟑 es la matriz de rigidez del elemento en coordenadas globales; 𝑏𝑏, ℎ son la
base y la altura de la sección transversal del elemento; 𝐿𝐿 es la longitud; 𝐸𝐸 es el módulo de
elasticidad; 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 es el valor de 𝑠𝑠𝑠𝑠𝑠𝑠 𝛼𝛼; 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐es el valor del cos 𝛼𝛼.
556 Roberto Aguiar Falconí
UFA - ESPE

15.2.2 Solución de un Pórtico Plano

El Análisis Matricial de Estructuras, está orientado a resolver una estructura en el

computador, sin embargo de ello y con el propósito de conocer exactamente cómo se resuelve
una estructura por el Método de los Desplazamientos, se desarrolla el siguiente ejercicio,

• EJEMPLO 3

Se considera que todos los elementos de la estructura de la figura 15.13, son

totalmente flexibles, se pide:

i. Encontrar el Vector de Cargas Generalizadas 𝑸𝑸, empleando el Algoritmo denominado:

Problema Primario y Problema Complementario.
ii. Indicar como se obtiene el Vector de Cargas Generalizadas 𝑸𝑸 utilizando la librería de
programas de CEINCI-LAB.
iii. Describir la forma como se halla la matriz de rigidez de la estructura con CEINCI-LAB.
iv. Presentar la solución total con CEINCI-LAB

No se desea elaborar un Programa que resuelva la estructura en forma completa como

se hizo con el Ejemplo 1, lo que se quiere es ilustrar como se utilizan las diferentes
subrutinas de CEINCI-LAB.

Figura 15.13 Pórtico Plano correspondiente al Ejemplo 3.

• SOLUCIÓN

A la izquierda de la figura 15.14 se presenta la numeración de los nudos, que se ha

considerado en la solución, en concordancia con esta numeración de los nudos, a la derecha
se indican los grados de libertad. Por ser los elementos totalmente flexibles existente tres
grados de libertad en cada nudo.
 ANALISIS ESTATICO CON CEINCI-LAB 557

Figura 15.14 Numeración de nudos, elementos y sistema de coordenadas 𝑸𝑸 − 𝒒𝒒

• Problema Primario

Por definición el Problema Primario, es aquel en el cual actúan todas las cargas y la
condición es que el vector 𝒒𝒒 = 𝟎𝟎. Es decir que actúan todas las cargas y los nudos interiores no
se desplazan ni rotan, para lograr esto es necesario que existan cargas (fuerza o momento) de
fijación 𝑅𝑅𝑖𝑖 que impiden ese grado de libertad.

Estas cargas de fijación se colocan con sentido contrario al sistema de cargas 𝑸𝑸

indicado en la figura 15.14. El Problema Primario finaliza con la determinación de estas cargas
de fijación y ahí se obtendrá la verdadera dirección que tienen. En la figura 15.15 se presenta
el problema primario de la estructura analizada.

Figura 15.15 Problema Primario

El Problema Primario tiene 2 etapas, que son equilibrio de elementos y equilibrio de

juntas.

• Equilibrio de elementos

Para el equilibrio de elementos se considera que cada elemento está empotrado-

empotrado debido a la condición de 𝒒𝒒 = 𝟎𝟎. En la figura 15.16 se presentan los elementos
empotrados-empotrados y las acciones de empotramiento perfecto, se ha utilizado un
subíndice el mismo que corresponde a la identificación de los nudos que se ha dado en la
figura 15.13. Nótese que los momentos de empotramiento perfecto se han colocado con el
signo que salen.
558 Roberto Aguiar Falconí
UFA - ESPE

Figura 15.16 Equilibrio de Elementos

Para encontrar las acciones de empotramiento perfecto, la carga vertical 𝑊𝑊 debe

descomponerse en sus dos componentes 𝑊𝑊𝑥𝑥 que es una carga uniforme axial al eje del
elemento y en la componente 𝑊𝑊𝑦𝑦 que es una carga transversal uniforme; en la figura 15.17 se
presentan estas cargas pero en un elemento diferencial. Las ecuaciones con las que se hallan
estas cargas uniformes distribuidas son:

Figura 15.17 Acciones de empotramiento perfecto de elemento BC.

 ANALISIS ESTATICO CON CEINCI-LAB 559

Para el elemento CD se tiene algo similar y se presenta a continuación

• Equilibrio de Junta B

Las acciones de empotramiento perfecto actúan con sentido contrario en la Junta por
este motivo se han cambiado de sentido. Para el equilibrio de la Junta la fuerza axial y de corte
se descomponen en sus componentes horizontal y vertical.
560 Roberto Aguiar Falconí
UFA - ESPE

• Equilibrio de Junta C

• Equilibrio de Junta D

• Problema Complementario

Solo actúan las cargas de fijación con sentido contrario al del Problema Primario para
que generen el vector de coordenadas generalizadas 𝒒𝒒. En la figura 15.18 se presenta el
Problema Complementario, que es el que se resuelve con Análisis Matricial de Estructuras.
 ANALISIS ESTATICO CON CEINCI-LAB 561

Figura 15.18 Problema Complementario

En el Problema Complementario se tienen las cargas solo en las juntas y en la

dirección del sistema 𝒒𝒒. Por lo que se cumple con la segunda hipótesis de estructuras. La
primera hipótesis se refiere al hecho de que solo interesa la posición de las juntas y para ello
utilizamos coordenadas para definir la posición de las juntas o nudos.

De la figura 15.18 se desprende que el vector de cargas 𝑸𝑸 es:

0
⎡−2.22⎤
⎢ ⎥
⎢−1.11⎥
⎢ 0 ⎥
𝑸𝑸 = ⎢−4.11⎥
⎢ 0 ⎥
⎢ 0 ⎥
⎢−2.22⎥
⎣ 1.11 ⎦

La segunda pregunta del Ejemplo 3, se refiere a como se determina este vector de

cargas utilizando el sistema de computación CEINCI-LAB. Para ello se van a ir indicando cada
una de las subrutinas que se requieren, con la entrada de datos y resultados esperados.

• Descripción de algunos programas

• [CG,ngl]=cg (nod,nr,RES)

El programa cg sirve para encontrar el CG y ngl de estructuras que tienen tres grados
de libertad por nudo o dos grados de libertad por nudo. Sirve entonces por ahora, para
pórticos planos o armaduras planas. Si es armadura plana el usuario solo indicará en la
matriz RES tres datos por fila y son el número de nudo, código en X, y código en Y. Se
destaca que código vale 1 o 0; el número 1 es para cuando no existe ese grado de
libertad; el código 0 es para cuando existe ese grado de libertad; para pórticos planos
se tienen 4 datos por cada nudo restringido: el número del nudo restringido y los
códigos para los desplazamientos en X, Y, y la rotación.

Los datos de entrada, son: nod que es el número de nudos de la estructura; nr que es
el número de nudos restringidos; RES la matriz que contiene las restricciones de los
nudos restringidos y que ha sido descrito en el párrafo anterior.

El programa reporta una matriz que contiene las Coordenadas Generalizadas de la

estructura CG, y el número de grados de libertad ngl.

Para la estructura del Ejemplo 3, los datos son:

>> nod = 5; nr = 2; RES=[1 1 1 1; 2 1 1 1];

562 Roberto Aguiar Falconí
UFA - ESPE

>> [CG,ngl]=cg (nod,nr,RES)

0 0 0
⎡0 0 0⎤⎥
𝐶𝐶𝐶𝐶 = ⎢⎢1 2 3⎥ ngl=9
⎢4 5 6⎥
⎣7 8 9⎦

Es conveniente numerar primero los nudos restringidos y después todos los demás
como se lo ha realizado en la figura 15.14. La primera fila de CG son los grados de
libertad del nudo 1; la segunda del nudo 2; la tercera del nudo 3; etc.

• [NI,NJ]=gn_portico (GEN)

El programa gn_portico sirve para obtener dos vectores que contienen al nudo inicial
NI y al nudo final NJ de cada uno de los elementos del pórtico plano. La forma de uso
del programa es: [NI,NJ]=gn_portico (GEN). Donde GEN es una matriz con la
siguiente información para cada uno de los elementos: número del elemento; nudo
inicial del elemento; nudo final del elemento; número de elementos a generar;
incremento en la numeración de los elementos; incremento en la numeración del nudo
inicial; incremento en la numeración del nudo final.

Para el Ejemplo 3 se indica el arreglo GEN sin generar los nudos para que sea más
fácil la explicación; a la derecha de la figura 15.14 se indica la numeración de los
elementos

>> GEN = [1 1 3 0 0 0 0;
2 3 4 0 0 0 0;
3 2 5 0 0 0 0;
4 4 5 0 0 0 0];

>> [NI,NJ]=gn_portico (GEN)

El programa reporta:

𝟏𝟏 𝟑𝟑
NI = �𝟑𝟑� NJ= �𝟒𝟒�
𝟐𝟐 𝟓𝟓
𝟒𝟒 𝟓𝟓

• [VC] = vc (NI,NJ,CG)

El programa vc determina los vectores de colocación de estructuras que tienen dos o

tres grados de libertad por nudo. La forma de uso es: [VC] = vc (NI,NJ,CG). Donde NI
es el vector que contiene a los nudos iniciales; NJ el vector con los nudos finales; CG
es una matriz con las coordenadas generalizadas. El programa reporta la matriz que
contiene a los vectores de colocación de cada uno de los elementos. El vector de
colocación está compuesto por los grados de libertad del nudo inicial y del nudo final.

A continuación se indica a la izquierda lo que reporta el programa, para el Ejemplo 3 y

a la derecha se recuerda la numeración de los elementos (número dentro de un círculo)
y los grados de liberad, para entender la matriz que contiene a los vectores de
colocación de cada uno de los elementos.
ANALISIS ESTATICO CON CEINCI-LAB 563

• [X,Y] = glinea_portico (Nudos)

El programa glinea_portico genera las coordenadas de los nudos, tanto en sentido X,

como en sentido Y. La forma de uso es: [X,Y] = glinea_portico (Nudos) Donde Nudos
es una matriz que contiene: Número de nudo; coordenada en X del nudo; coordenada
en Y del nudo; número de nudos a generar; incremento en la numeración del nudo;
incremento de longitud en sentido X; incremento de longitud en sentido Y. Reporta dos
vectores X,Y con las coordenadas en X y de cada nudo.

La entrada de datos, del Ejemplo 3, sin generar ningún nudo, es la siguiente.

El Programa reporta los vectores X, Y que tienen esta información.

𝟎𝟎 𝟎𝟎
⎡𝟔𝟔⎤ ⎡𝟎𝟎⎤
⎢ ⎥ ⎢ ⎥
X = ⎢𝟎𝟎⎥ Y = ⎢𝟑𝟑⎥
⎢𝟑𝟑⎥ ⎢𝟓𝟓⎥
⎣𝟔𝟔⎦ ⎣𝟑𝟑⎦

Con las coordenadas de los nudos se obtiene la longitud de los elementos; el seno y
coseno del ángulo 𝛼𝛼. Esto se lo realiza con el programa longitud que se indica más
adelante.

• dibujo (X,Y,NI,NJ)

El programa dibujo, presenta la estructura, es importante utilizar este programa para

verificar si se han ingresado bien los datos. La forma de uso es: dibujo (X,Y,NI,NJ)
Las variables de entrada de este programa ya fueron indicadas.

• [L,seno,coseno]=longitud (X,Y,NI,NJ)

El programa longitud determina la longitud; sen α y cos α de cada uno de los

elementos de una estructura plana. Como se indicó α es el ángulo que forma el eje
del elemento con el eje de las X. El programa reporta tres vectores, denominados L,
que contiene la longitud de cada uno de los elementos; seno con el valor de la función
seno de cada elemento y coseno con el valor del coseno de cada elemento.
564 Roberto Aguiar Falconí
UFA - ESPE

Para el Ejemplo 3, la ejecución del programa, reporta.

• [Q, Q2] = cargas (njc, nmc, ngl, L, seno, coseno, CG, VC, F, Fm, datos)

El programa cargas determina el vector de cargas generalizadas Q en pórticos con

cargas en las juntas y cargas en los elementos; la forma de uso de este programa es la
siguiente: [ Q, Q2]= cargas (njc, nmc, ngl, L,seno, coseno, CG, VC, F, Fm). Donde
njc es el número de juntas cargadas, nmc es el número de elementos cargados; ngl
es el número de grados de libertad; L vector con la longitud de cada uno de los
elementos; seno vector con el sen α de los elementos; coseno vector con el cos α
de los elementos; CG, matriz que contiene las coordenadas generalizadas de la
estructura; VC, matriz con los vectores de colocación de los elementos; F matriz con la
siguiente información: número del elemento cargado, fuerza horizontal (positiva hacia
la derecha), fuerza vertical (positiva hacia arriba) y el momento (anti horario positivo) de
cada una de las juntas cargadas; Fm matriz que contiene el número del elemento
cargado; la carga que actúa en el elemento; Código de Carga, actualmente solo trabaja
con carga uniforme distribuida, el código es 1; número de elementos a generar;
incremento en la numeración de los elementos; datos esta opción se utiliza cuando se
tiene cualquier tipo de carga en los elementos, en este caso, el usuario, determina
aparte las acciones de empotramiento perfecto y da como datos en la variable datos,
que tiene la siguiente información: número del elemento cargado; fuerzas axial y de
corte y momento en el nudo inicial; fuerza axial, de corte y momento en el nudo final,
son 7 datos por cada elemento cargado.

F = 0.
Si no hay cargas en las juntas se debe especificar njc = 0 y se debe indicar
Lo propio si no hay cargas en los elementos se debe indicar nmc = 0 y Fm = 0 .
Finalmente si no se desea ingresar las acciones de empotramiento perfecto, se
colocará datos = 0

Para el Ejemplo 3, los datos que todavía no se han indicado para usar el programa
cargas son:

>> njc=0; nmc=2; F=0; datos=0;

>> FM=[2 1.5 1 0 0 0;
4 1.5 1 0 0 0];

Los otros datos que se necesitan para ejecutar el programa cargas han sido ya
indicados. El programa reporta.
ANALISIS ESTATICO CON CEINCI-LAB 565

• [K] = krigidez (ngl, ELEM, L, seno, coseno, VC, E)

Programa que determina la matriz de rigidez de Pórticos Planos o de Armaduras

Planas. Las variables que han sido ya definidas son: ngl es el número de grados de
libertad de la estructura; L, es un vector con la longitud de cada uno de los elementos;
seno, coseno, vectores con los senos y cosenos de cada elemento; VC es la matriz
con los vectores de colocación de cada elemento. Las variables no definidas todavía,
son: E es el módulo de elasticidad del material; ELEM define la geometría de las
secciones transversales de los elementos.

Para armaduras planas, ELEM, está compuesto el área de la sección transversal de

cada uno de los elementos. Para pórticos planos, en ELEM se indica la base y la
altura de cada elemento.

• [FF]=fuerzas (ngl, ELEM, L, seno, coseno, VC, E, q, Q2)

De los datos de entrada, la variable no definida todavía es q que es el vector de

coordenadas generalizadas. Este vector se obtiene resolviendo el sistema de
ecuaciones lineales: Q = K q. En MATLAB la solución del sistema de ecuaciones se
realiza de la siguiente manera: q = K \ Q

Este programa determina las cargas finales en los elementos en coordenadas locales,
en una matriz que se ha denominado FF. Para armaduras planas la matriz FF tiene
cuatro cantidades que son las fuerzas axiales y de corte, del nudo inicial y final,
respectivamente; la fuerza de corte será cero en armaduras planas. Para pórticos
planos, FF contiene la fuerza axial, fuerza de corte, momento, del nudo inicial y del
nudo final. La primera fila de FF corresponde al primer elemento, la segunda fila al
segundo elemento, etc.

Se destaca que la matriz Q2 contiene las cargas de empotramiento perfecto de los

elementos en coordenadas locales y se encontró cuando se ejecutó el programa
cargas. Para una armadura plana la matriz Q2 es nula, debido a que en armaduras no
hay cargas en los elementos.
566 Roberto Aguiar Falconí
UFA - ESPE

15.2.3 Comentario de los Programas

La solución matricial de un pórtico plano, por el método de los desplazamientos, es

similar a la de una armadura plana, es por eso que se utilizan las mismas subrutinas. Así por
ejemplo la subrutina que obtiene la matriz CG que contiene los grados de libertad de un pórtico
plano es la subrutina cg, la misma que se utilizó para armadura plana. Evidentemente que en
esta subrutina existe un código con el cual se identifica a una armadura plana y a un pórtico
plano; ese código vale 2 para armaduras planas y 3 para pórticos planos.

En pórticos planos, existen cargas en los elementos por este motivo la solución total del
problema es igual a la solución del problema primario más la solución del problema
complementario. En el caso de armaduras planas la solución total es la solución del problema
complementario. Para el caso de pórticos planos se debe sumar la solución del problema
primario en que sus elementos se consideran que están empotrados. Aguiar (2004).

• EJEMPLO 4

Resolver, utilizando los programas del sistema de computación CEINCI-LAB, el pórtico

plano indicado en la figura 15.19, las columnas son de 30/40 y la viga de 30/30. Sobre la
estructura actúa una fuerza horizontal de 3 T., y una carga uniforme distribuida de 2 T/m. El
módulo de elasticidad E = 1800000 T / m 2 .
 ANALISIS ESTATICO CON CEINCI-LAB 567

Figura 15.19 Pórtico plano a resolver con CEINCI-LAB

• SOLUCIÓN

A la izquierda de la figura 15.20 se muestra la numeración de los nudos y elementos; a

la derecha se indican los grados de libertad. Las matrices CG (Coordenadas Generalizadas) y
VC (Vectores de colocación) son:

0 0 0
0 0 0 0 1 2 3 
0
VC = 0 0 0 4 5 6
0
CG = 
1 2 3
  1 2 3 4 5 6 
4 5 6

Figura 15.20 Numeración de nudos y elementos. Grados de libertad de la estructura.

El vector de cargas es:

Qt = [ 3.00 -4.50 -3.375 0.00 -4.50 3.375].

De la solución del sistema de ecuaciones K Q =q , se encuentra:

qt =[0.0026 -0.0001 -0.0017 0.0025 -0.0001 -0.0001].

568 Roberto Aguiar Falconí
UFA - ESPE

Tabla 15.1 Fuerzas y momentos finales en extremos de elementos

Elemento Nudo Inicial Nudo Final
Fuerza Fuerza de Momento Fuerza Fuerza de Momento
Axial Corte ( T.) ( T.m.) Axial Corte ( T.) ( T.m.)
( T. ) ( T. )
1 3.8762 0.0646 1.7255 -3.8762 -0.0646 -1.5316
2 5.1238 2.9354 4.4765 -5.1238 -2.9354 4.3386
3 2.9354 3.8762 1.5316 -2.9354 5.1238 -4.3386

En la tabla 15.1, se muestran las fuerzas y momentos en cada uno de los elementos,
en coordenadas locales. El archivo de datos del programa que resuelve el pórtico indicado es
el siguiente.

% Programa para resolver un Pórtico de 1 piso y 1 vano

%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
% 30 de agosto de 2012
%..............................................................
% portico_ej1
%..............................................................
nod=4; % Número de nudos del pórtico plano
nr=2; % Número de nudos restringidos
% MATRIZ DE RESTRICCIONES DE LOS NUDOS
RES=[ 1 1 1 1; % Nudo restringido Des X, Des Y, Rot
2 1 1 1];
[CG,ngl]=cg(nod,nr,RES) % CG Matriz de coord generalizadas
% MATRIZ DE GENERACION DE NUDOS DE LOS ELEMENTOS
GEN=[1 1 3 1 1 1 1; %elem, ni, nf, elem a gene, inc en elem, inc en ni, inc en nf
3 3 4 0 0 0 0];
[NI,NJ]=gn_portico(GEN);%Generacion de vectores NI (Nudo Inicial) y NJ (Nudo Final)
% MATRIZ DE GENERACION DE COORDENADAS DE LOS NUDOS
NUDOS=[1 0.0 0.0 1 1 4.5 0.0;% i, xi, yi, nudos a gener, incr num nudo, dx,dy
3 0.0 3.0 1 1 4.5 0.0];
[X,Y]=glinea_portico(NUDOS);% Generación de las coordenadas de los nudos
dibujo (X,Y,NI,NJ)% Dibuja el portico plano
[VC]=vc(NI,NJ,CG);% Vector de colocación
[L,seno,coseno]=longitud (X,Y,NI,NJ);% Longitud, Seno y Coseno de elementos
% Vector de cargas generalizadas
njc=1;F=[3 3 0 0];datos=0; % Cargas en las Juntas
nmc=1;Fm=[3 2.0 1 0 0]; %Elem carg, carga, código, elem a gener, incr numero elemento
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos)
% Secciones de los elementos
SECCION=[1 0.30 0.40 1 1; % Eleme, base, altura, elem a gener, increm en elem
3 0.30 0.30 0 0];
[ELEM]=gelem_portico(SECCION)
E=1800000; % Modulo de elasticidad del material
[K]=krigidez(ngl,ELEM,L,seno,coseno,VC,E);% Matriz de rigidez de la estructura
q=K\Q % Desplazamientos y Giros
[FF]=fuerzas(ngl,ELEM,L,seno,coseno,VC,E,q,Q2)% Fuerzas y momentos finales en los
elementos
% ---end---
 ANALISIS ESTATICO CON CEINCI-LAB 569

15.2.4 Cálculo de las Reacciones

En la forma como se resolvió el pórtico plano, del Ejemplo 4, no es posible calcular

directamente las reacciones de los apoyos, se debe hacer equilibrio en los apoyos y colocar las
fuerzas en los extremos de los elementos que se han obtenido. Si se desea encontrar
directamente, las reacciones se debe considerar como que no existe apoyo en la base y
numerar esos grados de libertad.

• EJEMPLO 5

En la figura 15.21, se presenta a la izquierda el pórtico resuelto en el apartado anterior

con los grados de libertad (gdl) indicados en la gráfica central. Ahora se va a resolver este
pórtico con los gdl., indicados a la derecha, con 12 gdl, considerando grados de libertad en los
empotramientos para hallar directamente las reacciones. Como se verá más adelante loa
grados de liberad del 1 al 6 que están relacionados con apoyos valen cero ya que es
empotramiento perfecto.

Figura 15.21 Modelo para calcular las reacciones de los apoyos directamente.

La ecuación básica de análisis matricial de estructuras, es la siguiente.

Q=Kq

Donde Q es el vector de cargas generalizadas; K es la matriz de rigidez de la

estructura y q es el vector de coordenadas generalizadas. Para encontrar directamente las
reacciones se debe reconocer dos tipos de coordenadas, las que están directamente
relacionadas con las reacciones, en el ejemplo sería de la 1 a la 6, a estas se denominan
coordenadas a, y los restantes grados de libertad que en este caso van del 7 al 12, a estas se
llaman coordenadas b. Luego la ecuación básica de análisis matricial, queda:

Q A + R   K AA K AB  q A 
 = K BB  q B 
(15.15)
 Q B   K BA

Siendo Q A el vector de cargas generalizadas asociadas a las coordenadas 1 a 6 pero

como no existen cargas en las columnas este vector vale cero. R es el vector que contiene las
reacciones de los apoyos; QB es el vector de cargas generalizadas asociado a las
coordenadas 7 a 12.

Con los grados de libertad considerados, la matriz de rigidez de la estructura es de 12

por 12 y la partición se realiza de acuerdo al número de coordenadas a y b ; se denominan
na al número de coordenadas a que en este caso es 6; nb al número de coordenadas b
570 Roberto Aguiar Falconí
UFA - ESPE

que también es 6. Luego las sub matrices K AA es de ( na, na ) ; K AB de ( na, nb) ; K BB de

(nb, nb) y K BA = K
t
AB . El vector qA contiene a los desplazamientos y giros de las
coordenadas a que vale cero ya que la estructura es empotrada y q B son los
desplazamientos y giros de las coordenadas b . Al desarrollar la ecuación (15.15) se tiene:

Q A + R = K AB * q B
Q B = K BB q B

De la segunda ecuación se obtiene qB resolviendo el sistema de ecuaciones indicado.

Al reemplazar qB en la primera ecuación se halla R .

R = K AB ∗ q B − Q A (15.16)

• SOLUCIÓN

El programa con el cual se obtiene directamente las reacciones de los apoyos de la

estructura indicada en la figura 15.20 se lista a continuación.

% Calculo Directo de las Reacciones de un Pórtico

%
% Dr. Roberto Aguiar
% 15 de diciembre de 2012
%----------------------------------------------
% portico_ej2
%----------------------------------------------
nod=4;nr=1;RES=[1 0 0 0];% Obligado colocar nr=1
[CG,ngl]=cg (nod,nr,RES);% CG con todos los grados de libertad
GEN=[1 1 3 1 1 1 1;3 3 4 0 0 0 0];
[NI,NJ]=gn_portico(GEN);
NUDOS=[1 0.0 0.0 1 1 4.5 0.0; 3 0.0 3.0 1 1 4.5 0.0];
[X,Y]=glinea_portico(NUDOS);
dibujo (X,Y,NI,NJ);
[L,seno,coseno]=longitud (X,Y,NI,NJ);
SECCION=[1 0.30 0.40 1 1; 3 0.30 0.30 0 0];
[ELEM]=gelem_portico(SECCION);
[VC]=vc(NI,NJ,CG);
E=1800000;
[K]=krigidez(ngl,ELEM,L,seno,coseno,VC,E);
Kaa=K(1:6,1:6);Kab=K(1:6,7:ngl);Kbb=K(7:ngl,7:ngl);
njc=1;F=[3 3 0 0];nmc=1;Fm=[3 2 1 0 0]; datos=0;
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos);
na=6;Qa=Q(1:na);
Qb=Q(na+1:ngl);
qb=Kbb\Qb;
qa=[0;0;0;0;0;0];
q=[qa;qb];
R=Kab*qb
[FF]=fuerzas(ngl,ELEM,L,seno,coseno,VC,E,q,Q2)
 ANALISIS ESTATICO CON CEINCI-LAB 571

El vector que contiene las reacciones R que reporta el programa es:

R t = [−0.0646 3.8762 1.7255 − 2.9354 5.1238 4.4675]

Los restantes resultados son los ya indicados en el problema anterior.

15.2.5 Estructuras con cargas puntuales

Actualmente algunas estructuras se están reforzando o al menos están comprobando si

sus elementos estructurales están en capacidad de resistir la coloción de un Data Center, que
tienen un considerable número de equipos que son bastante pesados, ver figura 15.22

Los fabricantes de estos equipos proporcionan los pesos de cada uno de ellos y
conjuntamente con el Arquitecto se procede a ubicarlos. Estos equipos pueden ser modelados
como cargas puntuales o como carga repartida en un determinado sitio.

Sea como sea el Proyectista estructural deberá hacer su modelo de cálculo con cargas
puntuales o con cargas uniformes o una mezcla de las dos y encontrar las acciones de
empotramiento perfecto que deberá dar como dato al programa cargas del sistema de
computación CEINCI-LAB.

En el Ejemplo que se presenta a continuación se ha considerado solo dos cargas

puntuales en la viga, para ilustrar el cálculo con cualquier tipo de carga.

Figura 15.22 Equipos de un Data Center que se pueden modelar como masas puntuales.

• EJEMPLO 6

Resolver el pórtico plano de la figura 15.23, que tiene cargas puntuales sobre la viga.
Las dimensiones están en metros y la sección transversal de sus elementos en cm. EL módulo
de elasticidad 𝐸𝐸 = 1800000 𝑇𝑇/𝑚𝑚2
572 Roberto Aguiar Falconí
UFA - ESPE

Figura 15.23 Pórtico con cargas puntuales.

• SOLUCION

Se debe encontrar el vector de empotramiento perfecto de la viga con dos cargas

puntuales, esto se lo hace con funciones de forma o simplemente viendo en una tabla donde se
encuentren las fórmulas para esas cargas. Se denomina 𝑸𝑸𝟐𝟐 al vector de empotramiento
perfecto en coordenadas locales, para la viga este vecor es el siguiente.

𝑸𝑸𝟐𝟐 = [𝟎𝟎. 𝟎𝟎 𝟐𝟐. 𝟎𝟎𝟎𝟎 𝟏𝟏. 𝟓𝟓𝟓𝟓𝟓𝟓 𝟎𝟎. 𝟎𝟎 𝟐𝟐. 𝟎𝟎𝟎𝟎 − 𝟏𝟏. 𝟓𝟓𝟓𝟓𝟓𝟓]

La numeración de los grados de libertad se indica a la izquierda de la figura 15.24; a la

derecha se presenta la numeración de nudos y elementos. Se ha colocado un apoyo fijo y un
rodillo para ilustrar la entrada de datos en los nudos restringuidos.

Figura 15.24 Coordenadas generalizadas, numeración de nudos y elementos.

El programa con que se resuelve el Ejemplo 6 se lista a continuación.

% Solución de pórtico con cargas puntuales en un elemento

%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
% 27 de septiembre de 2013
%...............................................
% portico_ej3
%...............................................
 ANALISIS ESTATICO CON CEINCI-LAB 573

nod=4; % Número de nudos del pórtico plano

nr=2; % Número de nudos restringidos
RES=[ 1 1 1 0; % Nudo restringido Des X, Des Y, Rot
2 0 1 0];
[CG,ngl]=cg(nod,nr,RES) % CG Matriz de coord generalizadas
GEN=[1 1 3 1 1 1 1;
3 3 4 0 0 0 0];
[NI,NJ]=gn_portico(GEN) % Nudo Inicial NI y Nudo Final NJ
NUDOS=[1 0.0 0.0 1 1 4.5 0.0;% i, xi, yi, nudos a gener, incr num nudo, dx,dy
3 0.0 3.2 1 1 4.5 0.0];
[X,Y]=glinea_portico(NUDOS) % Generación de las coordenadas de los nudos
dibujo (X,Y,NI,NJ)% Dibuja el Pórtico
[VC]=vc(NI,NJ,CG) % Vector de colocación
[L,seno,coseno]=longitud (X,Y,NI,NJ) %Longitud, seno y coseno de los elementos
njc=0;F=0; % Número de juntas Cargadas
datos=[0 2 1.5556 0 2 -1.556]; % N V M N' V' M' de empotramiento perfecto
nmc=1; % Número de miembros cargados
Fm=[3 0 4 0 0]; %Elem carg, carga, código, elem a gener, incr numero elemento
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos)
SECCION=[1 0.30 0.40 1 1; % Eleme, base, altura, elem a gener, increm en elem
3 0.30 0.30 0 0];
[ELEM]=gelem_portico(SECCION)
E=1800000; % Modulo de elasticidad del material
[K]=krigidez(ngl,ELEM,L,seno,coseno,VC,E) % Matriz de rigidez de la estructura
q=K\Q % Desplazamientos y Giros
[FF]=fuerzas(ngl,ELEM,L,seno,coseno,VC,E,q,Q2) % Fuerzas y momentos finales en los
elementos
% ---end---

Se ha presentado el vector de cargas generalizadas 𝑄𝑄 y la matriz que contiene las

acciones de empotramiento perfecto 𝑸𝑸2 ; luego se indica la matriz de rigidez 𝑲𝑲; finalmente se
574 Roberto Aguiar Falconí
UFA - ESPE

muestra el vector de coordenadas generalizadas 𝒒𝒒 que contiene los desplazamientos y giros y

las fuerzas finales en coordenadas locales, en los elementos 𝑭𝑭𝑭𝑭.

15.2.6 Reforzamiento sísmico de estructuras

Una forma de reforzar sísmicamente una estructura es colocar diagonales de acero

entre las columnas, de esta manera se rigidiza notablemente el pórtico. Lo que se desea
ilustrar en este apartado es como se resuelve un pórtico de hormigón armado con diagonales
de acero, sometido a cargas verticales (carga muerta más carga viva) y cargas horizontales
(sismo).

Ese es el objetivo, razón por la cual no se profundiza en el reforzamiento pero vale la

pena mencionar que el Proyectista Estructural deberá tener en cuenta los problemas de torsión
en planta ya que cuando se colocan diagonales se tienen pórticos muy rígidos y el otro aspecto
que se debe chequear es el problema de pandeo cuando las diagonales trabajan a compresión,
debido a que la fuerza axial que resisten es muy alta.

En la figura 15.25 se presenta una estructura de hormigón armado reforzada con

diagonales de acero en forma de “V” invertida, ubicada en la ciudad de Viña del Mar, en Chile.
El ejemplo que se resuelve a continuación está inspirado en la estructura de la figura 15.25.

Figura 15.25 Estructura reforzada en Viña del Mar, con diagonales de acero.

• EJEMPLO 7

Resolver la estructura de la figura 15.26, considerando un módulo de elasticidad del

hormigón de 1600000 T/m2 y del acero 20000000 T/m2. Las columnas y vigas de hormig[on son
de 20/20 cm, y las diagonales de acero un perfil tubular de 100/100/4 mm.
 ANALISIS ESTATICO CON CEINCI-LAB 575

Figura 15.26 Estructura mixta de hormigón armado y diagonales de acero

• SOLUCIÓN

En la figura 15.27 se presentan los grados de libertad de la estructura analizada y en la

figura 15.28 la numeración de los nudos y elementos. Nótese que se ha considerado un nudo
en el punto medio de la viga central donde llegan las dos diagonales de hormigón.

Por otra parte, se ha numerado en primer lugar todos los elementos de hormigón y
después los elementos de acero; de esa manera se encontrará en primer lugar la contribución
de los elementos de hormigón a la matriz de rigidez y posteriormente la contribución de los
elementos de acero.

Figura 15.27 Grados de libertad

Figura 15.28 Numeración de nudos y elementos

Las diagonales de acero tienen conexiones simples, por lo tanto lo más fácil es que se
576 Roberto Aguiar Falconí
UFA - ESPE

sumistre al programa los vectores de colocación. A continuación se presenta el programa

desarrollado en MATLAB utilizando el sistema de computación CEINCI-LAB

% Solución de pórtico con vigas y columnas de hormigón y diagonales

% de acero.
%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
% 16 de mayo de 2014
%...............................................
% portico_ej4
%...............................................
nod=9; nr=4;
RES=[ 1 1 1 1; % Nudo restringido Des X, Des Y, Rot
2 1 1 1; 3 1 1 1; 4 1 1 1];
[CG,ngl]=cg(nod,nr,RES); % CG Matriz de coord generalizadas
GEN=[1 1 5 1 1 1 1;3 3 8 1 1 1 1;
5 5 6 3 1 1 1;9 2 7 1 1 1 0];
[NI,NJ]=gn_portico(GEN); % Nudo Inicial NI y Nudo Final NJ
NUDOS=[1 0.0 0.0 3 1 4.0 0.0;% i, xi, yi, nudos a gener, incr num nudo, dx,dy
5 0.0 3.0 1 1 4.0 0.0; 7 6.0 3.0 1 1 2.0 0.0;
9 12 3.0 0 0 0 0];
[X,Y]=glinea_portico(NUDOS); % Generación de las coordenadas de los nudos
dibujo (X,Y,NI,NJ)% Dibuja el Pórtico
[VC]=vc(NI,NJ,CG); % Vector de colocación
[L,seno,coseno]=longitud (X,Y,NI,NJ); %Longitud, seno y coseno de los elementos
njc=1;F=[9 0 0 0]; % Número de juntas Cargadas
datos=0;
nmc=4; % Número de miembros cargados
Fm=[5 1.5 1 1 3; %Elem carg, carga, código, elem a gener, incr num elemen
6 1.5 1 1 1];
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos);
SECCION=[1 0.20 0.20 7 1; % Eleme, base, altura, elem a gener, increm en elem
9 0.0001 0.0001 1 1];
[ELEM]=gelem_portico(SECCION);
E=1600000; % Modulo de elasticidad del material
[KH]=krigidez(ngl,ELEM,L,seno,coseno,VC,E); % Matriz de rigidez solo de elementos de acero
% Contribucion de diagonales de acero a la matriz de rigidez
LA=L(9:10);senA=seno(9:10);cosA=coseno(9:10);area=0.001536;
VCA=[0 0 7 8;0 0 7 8]; Es=20000000; ELEMA=[area;area];
[KA]=krigidez(ngl,ELEMA,LA,senA,cosA,VCA,Es);
KT=KH+KA;
q=KT\Q; % Desplazamientos y Giros
[FH]=fuerzas(ngl,ELEM,L,seno,coseno,VC,E,q,Q2) % Fuerzas y mom en Hormigon
Q2A=zeros(2,4)
[FA]=fuerzas(ngl,ELEMA,LA,senA,cosA,VCA,Es,q,Q2) % Fuerzas y mom en Acero
% ---end---

En la parte final del progama portico_ej4 se observa que se ha ejecutado dos veces el
programa fuerzas, en la primera se obtienen fuerzas y momentos en coordenadas locales de
los elementos de hormigón y en la segunda solo lasfuerzas en coordenadas locales de los
elementos de acero. Nótese que en este último caso se debe dar como dato Q2A una matriz de
2 por 4 con ceros, ya que no existen cargas de empotramiento perfecto en los elementos de
acero pero es necesario para que pueda ejecutarse el programa fuerzas. Los resultados se
indican a continuación.
 ANALISIS ESTATICO CON CEINCI-LAB 577

15.3 ESTRUCTURAS DE ACERO CON NUDOS PRECALIFICADOS

Uno de los puntos críticos en el diseño de estructuras de acero, tiene que ver con las
conexiones viga-columna. Esto se puso en evidencia luego del sismo de Northridge de 1994 y
Kobe de 1995, donde fallaron los nudos, lo que obligó a investigar más en esta temática y es
así como el ANSI/AISC 358-16 presenta 9 conexiones “Precalificadas” que han sido ensayadas
y pueden utilizarse en el diseño sismo resistente de estructuras.

Una de las conexiones precalificadas es la Bolted Flange Plate (BFP) Moment

Connection, que es una conexión con placas empernadas, como se observa en la figura 15.29,
que corresponde al edificio donde funciona la Unidad de Vigilancia Comunitaria UVC, de la
ciudad de Manta. Esta estructura fue probada durante el terremoto del 16 de abril de 2016, que
tuvo una magnitud de 7.8 y ninguno de sus nudos falló se comportaron muy bien.

Con este tipo de conexiones precalificadas, el análisis de una estructura de acero se

realiza en forma similar al de una estructura de hormigón armado en el sentido de considerar
tres grados de libertad por nudo; es así como se resuelve el ejemplo 8.

Figura 15.29 Estructura de acero con un nudo “Precalificado” que tiene capacidad a flexión,
corte y carga axial. Fuente: Aguiar et al. (2016).

En la figura 1.30 se presenta la planta baja del edificio de UVC de Manta, se aprecian
dos columnas en sentido transversal, tipo “H”, y una viga tipo “I” que tiene dos voladizos a los
extremos. La estructura sufrió un daño considerable en las paredes que se encontraban en la
parte inferior de los voladizos, simplemente porque no había elementos confinante.
578 Roberto Aguiar Falconí
UFA - ESPE

Figura 15.30 Orientación de las columnas tipo “H” y de la viga “I” que conforman un pórtico en
sentido transversal

En la figura 15.30 se aprecia que solo existen vigas secundarias perpendiculares al

pórtico transversal de tal manera que son estos pórticos los que soportan la carga transversal.
Las vigas secundarias limitan la longitud no arriostrada del ala superior de la viga “I”, con ello
se tiene una mayor capacidad a flexión. Nótese que el ala inferior no se encuentra arriostrada
de tal manera que no se está considerando la inversión de momentos durante un sismo.

• EJEMPLO 8

Todas las columnas de la estructura de acero indicada en la figura 15.31 son de un

perfil tipo “H” de 466/280/18/10 mm. Las vigas tipo “I” son de 530/200/15/10 mm, ver figura
15.32. A la izquierda de la figura 15.31 se indica la carga vertical que gravita en la estructura.
Se pide encontrar los desplazamientos y giros en los extremos del voladizo; y las fuerzas y
momentos en los elementos de la planta baja.

Figura 15.31 Carga Vertical que actúa en la estructura y numeración de nudo y elementos.
 ANALISIS ESTATICO CON CEINCI-LAB 579

Figura 15.32 Dimensiones de elementos estructurales y área e inercia de los mismos.

• SOLUCIÓN

No se debe pensar que la inercia de las columnas es el mayor valor de las dos inercias
que se tienen según la dirección X, Y. Se debe observar la forma como están orientadas. En el
presente caso si es el mayor valor para el pórtico analizado pero si se estaría analizando en el
sentido transversal es el menor valor. A continuación se lista el programa desarrollado en
CEINCI-LAB para resolver el ejercicio.

% Estructura de acero UVC de Manta

% Dr. Roberto Aguiar
% 31 de enero de 2016
%.....................................
nod=12;nr=2;RES=[1 1 1 1;2 1 1 1]; Es=21000000;
[CG,ngl]=cg (nod,nr,RES);
GEN=[1 1 4 1 1 1 1;3 4 7 1 1 1 1;5 7 9 3 1 1 1 ;
9 3 4 2 1 1 1; 12 7 8 2 1 2 2];
NUDOS=[1 2.65 0.0 1 3 0.0 4.18;2 11.75 0.0 1 3 0.0 4.18;
3 0.0 4.18 1 1 2.65 0.0; 5 11.75 4.18 1 1 1.5 0.0;
7 2.65 7.88 2 2 0.0 3.7; 8 11.75 7.88 2 2 0.0 3.7];
[NI,NJ]=gn_portico(GEN); % Vector de nudo inicial y final
[X,Y]=glinea_portico(NUDOS);% Coordenadas X, Y de los nudos
dibujo(X,Y,NI,NJ)
[L,seno,coseno]=longitud (X,Y,NI,NJ);% Longitud seno y cosen de
elemen.
SECCION=[1 0.0144 0.000572 7 1;9 0.0110 0.000502 5 1];
[ELEM]=gelem_portico(SECCION);
[VC]=vc(NI,NJ,CG); % Vector de colocacion de elementos
[K]=krigidez_acero(ngl,ELEM,L,seno,coseno,VC,Es);
njc=0;F=0;nmc=6; Fm=[9 3.57 1 0 0;10 3.57 1 0 0;
11 3.57 1 0 0;12 3.31 1 2 1];datos=0;
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos);
q=K\Q;
[FF]=fuerzas_acero(ngl,ELEM,L,seno,coseno,VC,Es,q,Q2)

En la tabla 15.2 se presentan los desplazamientos y giros de los voladizos y en las

figuras 15.33 los valores de las fuerzas y momentos en los elementos de la planta baja.
580 Roberto Aguiar Falconí
UFA - ESPE

Tabla 15.2 Desplazamientos y Giros en los voladizos

Nudo Desplazamiento Desplazamiento Rotación
Horizontal Vertical (rad)
(m) (m)
3 -0.0008 -0.0029 0.0009
6 -0.0007 0.0001 0.0007

Figura 15.33 Fuerzas y Momentos en los elementos de la planta baja.

Para verificar el equilibrio de los nudos, se necesita conocer las fuerzas y momentos de
las columnas del primer piso alto, que no se ha presentado.