EPN

ANÁLISIS ESTRUCTURAL AVANZADO

DEBER 13

PANELES A CORTE

POR:

DIEGO QUIZANGA
 ESCUELA POLITÉCNICA NACIONAL

MAESTRÍA DE INGENIERÍA ESTRUCTURAL

Por Diego Quizanga

PANEL A CORTE Y PANEL A FLEXION

Objetivo

El presente trabajo tiene como objetivo mostrar un ejercicio para elementos tipo pórtico
con un panel a corte

También se realizará un resumen del artículo “Deficiencies in the wide column analogy for
shearwall core analysis”

Resumen Artículo

(DEFICIENCIES IN THE WIDE COLUMN ANALOGY FOR SHEARWALL CORE

ANALYSIS)

Autor: B. Stafford Smith y Amal M. Girgis

La analogía de la columna ancha fue ideada para el análisis de paredes de corte

acopladas hace aproximadamente 3 décadas, se proporciona un elemento lineal que
representa bien la interacción de los muros de corte atráves de sus vigas de conexión y
permite el análisis mediante el uso de un computador usando un programa tipo pórtico
 muro representado
por una linea

brazos rígidos

columna a la cual se le asigna

un area e inercia del muro representado

Figura 1 Analogía de la columna ancha.

Recientemente se ha propuesto el uso de la analogía de la columna ancha, la que puede

extenderse al ensamblado de muros de corte conectados continuos, como diafragmas por
ejemplo diafragmas en elevadores

Como parte del estudio de comportamiento de diafragmas en elevadores los autores, han
testado la analogía del marco de columna ancha y han encontrado algunas deficiencias
graves en su desempeño.

Estas deficiencias se hicieron especialmente evidentes, en el análisis de muros con núcleos

cerrados total o parcialmente sometidos a torsión. Ver figura 2(a) y (2b).

(a) (b)

Figura 2 Núcleos de ascensores: (a) Núcleo cerrado; (b) Núcleo cerrado parcialmente.
Las deficiencias incluyen rotaciones excesivas, inexactitud de esfuerzos.

Los errores, han sido determinados producto de una deficiencia inherente en el modelo,
estos errores podrían en ciertas circunstancias tener magnitudes mayores, por lo que los
autores sugieren precaución en el uso del modelo

FUNCIONAMIENTO DEL MODELO DE COLUMNA ANCHA SIMPLE

El uso de la analogía de columna ancha para una estructura de pared acoplada se basa en
asumir que las deformaciones por corte de la pared son despreciables.

 Los valores del área de la sección y la inercia del muro son asignados a sus
correspondientes columnas del modelo de columna ancha de modo que las rigideces axial
y de deflexión de los muros sean representados.

La figura 3 muestra el modelo para la columna ancha para muros acoplados, la flecha
horizontal arriba indica que la altura del piso regresa flexionando las columnas
total
Flexión

Figura 3 Modelo de la columna ancha en paredes acopladas.

Si se usa el modelo de columna de ancha para representar dos paredes a corte conectadas
continuamente como se muestra en la figura 4(a) y 4(b), el corte continuamente
distribuido a lo largo del borde vertical común de las dos paredes puede ser agrupado en
el modelo de columna ancha dentro de la fuerza de corte vertical concentrada en el nivel
de las vigas.

Distribución
continua de corte a Fuerzas de corte
lo largo de la agrupadas a
conexión brazos rígidos

Figura 4 (a) Dos paredes a corte conectadas continuamente y (b) Modelo de columna ancha en
paredes.

 Este efecto de agrupamiento puede causar componentes de flexión inversa en las

columnas entre los niveles de las vigas con errores resultantes en los esfuerzos de flexión
del muro.

Los esfuerzos de flexión en el muro pueden ser corregidos con una precisión aceptable
promediando los esfuerzos de flexión superiores e inferiores de la columna de cada piso
eliminando el efecto de la inversión de momentos.

Los tipos de muros de corte acoplados para los cuales se pueden utilizar la metodología
con un error casi despreciable en los desplazamientos, y con esta simple corrección de los
esfuerzos de flexión son paredes de corte planas, no planas y con núcleos de sección
abierta
El análisis de secciones cerradas de paneles y paneles con aberturas parcialmente
cerradas con vigas o losas es ciertamente menos tolerante a las deficiencias que presenta
el método de la columna ancha

Mientras un panel con sección abierta resiste torque en el plano de flexión de los muros
de corte, un panel con sección abierta resiste torque mayormente al flujo de corte
circulando alrededor de la sección

Los núcleos de paneles parcialmente cerrados representan un caso intermedio en la

torsión y se da por la combinación de flexión en la pared y el flujo de corte, la proporción
tomada por cada uno depende de la rigidez relativa a la flexión de las paredes y de los
miembros de conexión. 

La rigidez relativa de vigas o losas a través de las aberturas determina el mayor grado de
proporción del torque resistido por el flujo de corte y una menor proporción resistida por
la flexión del muro, consecuentemente, si se modela un núcleo de sección cerrada o
parcialmente cerrada por un marco de columna ancha, la proporción reducida de la
deformación a flexión de la pared contribuye a la rotación, junto con el aumento en la
magnitud del esfuerzo del flujo de corte, causa deformaciones por corte y giros que son
ciertamente significativos

Un extremo caso de erróneos desplazamientos por corte ocurrirá en un panel de sección

poligonal regular cerrado y sujeto a torsión en un núcleo de sección poligonal y regular
cerrado sujeto a torsión, en este caso no habría flexión en las paredes y la rotación sería
enteramente atribuible a las actuales deformaciones por corte.
 Fig 5
Figura 5 Modelo de columna ancha en paredes de un núcleo de sección cuadrada sujeto a
torsión.

CONCLUSIONES

La analogía de la columna ancha no debe ser utilizada para el análisis de núcleos

de pared de corte con una sección transversal cerrada o parcialmente cerrada sin
embargo puede ser usado para el análisis de muros de corte planas y no planas
continuamente conectados

brazo rígido brazo rígido brazo rígido

columna columna

enlace columna enlace

arriostraminetos
arriostramineto

Fig 6a Fig 6b Fig 6c

Figura 6 Analogías de marcos: (a) Columna articulada, (b) Columna ancha arriostrada, y (c) Marco
arriostrado.
 EJERCICIO

CASO DE ESTUDIO.

Se considera un marco plano mostrado en la figura I, se asume que todos los elementos
tienen la misma sección transversal y el mismo módulo de elasticidad

Figura I
RESOLUCIÓN EN MATLAB

1) Geometría y materiales

Hormigón:

Se debe definir los datos del material como el módulo de elasticidad E y de corte G .

E
G=
E=13500 √ f ' c 2(1+μ )

El marco es de hormigón armado y todos los elementos tienen la siguiente sección:

Se define la inercia de esta sección, el área y el área de corte como se indica a

continuación para elementos tipo pórtico:
% HORMIGON
E=135000*((240)^.5); % MODULO DE ELASTICIDAD EN T/m2
G=E/(2*(1+0.2)); % MODULO DE CORTE EN T/m2
A=0.3*.4; % AREA DE LAS COLUMNAS Y VIGAS EN m2
Ac=A/(1.2); % AREA DE CORTE EN m2
Iz=(0.3*.4^3)/12; % INERCIA DE LAS COLUMNAS Y VIGAS EN m4

Mampostería:

Se debe definir los datos del panel de corte como el módulo de elasticidad E y de
corte G de la mampostería y fm

Fm =20 Kg/cm2

E=400 fm G=0 . 4∗E

Espesor de la pared=0.15 m
Se define la inercia de esta sección, el área y el área de corte como se indica a
continuación para elementos tipo panel:
% MAMPOSTERIA
e=0.15 ; %ESPESOR DE LA PARED EN METROS
fm=200; % FM DE LA MAMPOSTERIA EN T/m2
Em=400*fm; % MODULO DE ELASTICIDAD DE LA MAMPOSTERIA SEGUN NORMA MEXICANA
T/m2
Gm=0.4*Em; % MODULO DE CORTE MAMPOSTERIA T/m2

Matrices de rigidez de los materiales en sistema reducido

La matriz de rigidez de un elemento marco plano considerando deformaciones por corte

es:

a1∗EI 2+ β b∗EI 1−β

[ ]
0
2 L 1+ β 2 L 1+ β
a 2∗EI 2+ β
k = b∗EI 1−β 0
2 L 1+ β 2L 1+ β
AE
0 0
L

Donde:

6 EI A
β= Ac=
A c GL 2 f

Con estos datos se puede utilizar la rutina [K_MP]=R_MP(A,Ac,Iz,CG,L,E,G) para encontrar

la matriz de rigidez de cada elemento.

En matlab para este ejercicio se tiene:

%FUNCION R_MP
%Matriz de rigidez de elementos marco plano
%[k1]=R_MP(A,Ac,Iz,CG,L,E,G)
[k1]=R_MP(A,Ac,Iz,[4,4,2],8,E,G);
[k2]=R_MP(A,Ac,Iz,[4,4,2],sqrt(32),E,G);
[k3]=R_MP(A,Ac,Iz,[4,4,2],sqrt(32),E,G);
[k4]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k5]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k6]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k7]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k8]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k9]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
%Matriz de rigidez de elemento panel
kp=Gm*(e*3.6*3.6);

2 Matrices de transformación de sistema reducido a coordenadas de nudo

Para pasar del sistema de coordenadas reducido al sistema de coordenadas de nudos se

emplea la matriz de trasformación que se muestra a continuación:

−S 1/L −S1 /L −C 1

[
T=¿ C 1/L C1 /L −S 1 ¿ ¿¿
(1−dx1/L) −dx 1/L −dy 1
¿
]
¿

Donde:

C1 =cos(α−β 1 ) S 1 =sen (α−β 1 )

C2 =cos(α−β 2 ) S 2 =sen(α−β 2 )

Para encontrar la matriz de transformación de cada elemento se emplea la rutina

MP_TTah

En matlab para este ejercicio se tiene:

%Matriz de transformación
[T1,Ta1,h1]=MP_TTah(8,90,0,0,[0 0 0 0]);
[T2,Ta2,h2]=MP_TTah(sqrt(32),45,0,0,[0 0 0 0]);
[T3,Ta3,h3]=MP_TTah(sqrt(32),315,0,0,[0 0 0 0]);
[T4,Ta4,h4]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T5,Ta5,h5]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T6,Ta6,h6]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T7,Ta7,h7]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T8,Ta8,h8]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T9,Ta9,h9]=MP_TTah(4,90,0,0,[0 0 0 0]);
[Tp]=TP_Corte(4,4); % VA HASTA LOS EJES
 Matrices de rigidez en coordenadas de nudo

Para encontrar la matriz de rigidez de los elementos en coordenadas de nudos se aplica el

triple producto mostrado a continuación:
T
K=T ×k×T

En matlab para este ejercicio se tiene:

% RIGIDECES EN COORDENADAS DE NUDOS
K1=T1*k1*T1';
K2=T2*k2*T2';
K3=T3*k2*T3';
K4=T4*k4*T4';
K5=T5*k5*T5';
K6=T6*k6*T6';
K7=T7*k7*T7';
K8=T8*k8*T8';
K9=T9*k9*T9';
Kp=Tp*kp*Tp';

Ensamble de las matrices de rigidez de los elementos para formar la matriz de

rigidez de la estructura

Los vectores de ensamble son muy importantes para el cálculo de la matriz de rigidez de la
estructura, corresponde a los grados de libertad del nudo inicial (1) y el nudo final (2) de
cada elemento. En este ejemplo por ser un marco plano cada elemento tendrá un vector
de ensamble de seis valores, tres en cada nudo.

A continuación se muestran los vectores de ensamble para cada elemento:

S=zeros(ngl);
B=zeros(nglt,1);
% Vectores de esnsamble
lee1=[17 18 19 1 2 3];
lee2=[1 2 3 4 5 6];
lee3=[4 5 6 7 8 9];
lee4=[7 8 9 7 10 11];
lee5=[12 13 14 12 15 16];
lee6=[12 13 14 7 8 9];
lee7=[12 15 16 7 10 11];
lee8=[20 21 22 12 13 14];
lee9=[23 24 25 12 15 16];
leep=[8 7 10 13 12 15];
% FUNCION Ens2
[S,B]=Ens2(lee1,K1,ngl,S,B);
[S,B]=Ens2(lee2,K2,ngl,S,B);
[S,B]=Ens2(lee3,K3,ngl,S,B);
[S,B]=Ens2(lee4,K4,ngl,S,B);
[S,B]=Ens2(lee5,K5,ngl,S,B);
[S,B]=Ens2(lee6,K6,ngl,S,B);
[S,B]=Ens2(lee7,K7,ngl,S,B);
[S,B]=Ens2(lee8,K8,ngl,S,B);
[S,B]=Ens2(lee9,K9,ngl,S,B);
[S,B]=Ens2(leep,Kp,ngl,S,B);

Aplicamos las cargas para posteriormente obtener los desplazamientos en los

GDL

De la figura I tomamos el siguiente vector de cargas

B(1)=10; B(5)=-5; B(7)=5; B(12)=5;

El momento 20 es positivo porque gira en sentido contrario a las manecillas del reloj

Con la matriz de rigidez y el vector de cargas del sistema se pueden hallar los
desplazamientos:
−1
Δ L =S ( P L )

El vector de cargas PL está conformado por los primeros nueve valores del vector B. Por lo
tanto para encontrar los desplazamientos se aplica la siguiente orden:

Delta=S\B(1:ngl)
% SOLUCION DEL SISTEMA DE ECUACIONES
B(1)=10; B(5)=-5; B(7)=5; B(12)=5;
fprintf ('DESPLAZAMIENTOS EN LOS GDL\n ')
format long;Delta=S\B(1:ngl)

Para calcular las acciones en los elementos se debe tener una matriz que contenga los
desplazamientos de cada miembro.

Primero se tiene que ampliar la matriz de desplazamientos Δ de 16x1 a una de 25x1 para
considerar las restricciones. Lo que se hace es añadir ceros en las filas de la 17 a la 25 de la
siguiente manera:
Delta(25)=0

Con la ayuda de los vectores de ensamble se calculan los desplazamientos que

corresponden a cada elemento de la siguiente manera:

D1=Delta(lee1); D2=Delta(lee2); D3=Delta(lee3); D4=Delta(lee4); D5=Delta(lee5);

D6=Delta(lee6); D7=Delta(lee7); D8=Delta(lee8); D9=Delta(lee9);

Para calcular los desplazamientos de los elementos en el sistema de coordenadas

reducido se aplica la matriz de transformación T:

d1=T1'*D1; d2=T2'*D2; d3=T3'*D3; d4=T4'*D4; d5=T5'*D5; d6=T6'*D6; d7=T7'*D7;

d8=T8'*D8; d9=T9'*D9;

Finalmente se puede encontrar las acciones de los elementos aplicando la ecuación:

a=k ×d

a1=k1*d1; a2=k2*d2; a3=k3*d3; a4=k4*d4; a5=k5*d5; a6=k6*d6; a7=k7*d7; a8=k8*d8;

a9=k9*d9;

En matlab para este ejercicio se tiene:

%ACCIONES FINALES
Delta(nglt)=0;
D1=Delta(lee1);d1=T1'*D1;a1=k1*d1;
D2=Delta(lee2);d2=T2'*D2;a2=k2*d2;
D3=Delta(lee3);d3=T3'*D3;a3=k3*d3;
D4=Delta(lee4);d4=T4'*D4;a4=k4*d4;
D5=Delta(lee5);d5=T5'*D5;a5=k5*d5;
D6=Delta(lee6);d6=T6'*D6;a6=k6*d6;
D7=Delta(lee7);d7=T7'*D7;a7=k7*d7;
D8=Delta(lee8);d8=T8'*D8;a8=k8*d8;
D9=Delta(lee9);d9=T9'*D9;a9=k9*d9;
Dp=Delta(leep);Gama=Tp'*Dp;Tao=Gm*Gama;
En MATLAB más fácil trabajar en el editor que en el comand window se hizo este ejercicio
en un archivo .m
fprintf (' PANEL CORTE.m\n ')
fprintf ('PORTICO CON CARGAS PUNTUALES Y UN PANEL A CORTE \n ')
fprintf (' \n ')
fprintf ('ESCUELA POLITECNICA NACIONAL \n ')
fprintf ('MAESTRIA EN ESTRUCTURAS \n ')
fprintf (' \n ')
fprintf ('Por Diego Quizanga \n ')
clear
% MATERIALES
% HORMIGON
E=135000*((240)^.5); % MODULO DE ELASTICIDAD EN T/m2
G=E/(2*(1+0.2)); % MODULO DE CORTE EN T/m2
A=0.3*.4; % AREA DE LAS COLUMNAS Y VIGAS EN m2
Ac=A/(1.2); % AREA DE CORTE EN m2
Iz=(0.3*.4^3)/12; % INERCIA DE LAS COLUMNAS Y VIGAS EN m4
% MAMPOSTERIA
e=0.15 ; %ESPESOR DE LA PARED EN METROS
fm=200; % FM DE LA MAMPOSTERIA EN T/m2
Em=400*fm; % MODULO DE ELASTICIDAD DE LA MAMPOSTERIA SEGUN NORMA MEXICANA
T/m2
Gm=0.4*Em; % MODULO DE CORTE MAMPOSTERIA T/m2
% ESTRUCTURA
ngl=16; % NUMERO DE GRADOS DE LIBERTAD
nglt=25; % NUMERO DE GRADOS DE LIBERTAD MAS RESTRICCIONES

%FUNCION R_MP
%Matriz de rigidez de elementos marco plano
%[k1]=R_MP(A,Ac,Iz,CG,L,E,G)
[k1]=R_MP(A,Ac,Iz,[4,4,2],8,E,G);
[k2]=R_MP(A,Ac,Iz,[4,4,2],sqrt(32),E,G);
[k3]=R_MP(A,Ac,Iz,[4,4,2],sqrt(32),E,G);
[k4]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k5]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k6]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k7]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k8]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
[k9]=R_MP(A,Ac,Iz,[4,4,2],4,E,G);
%Matriz de rigidez de elemento panel
kp=Gm*(e*3.6*3.6);
% FUNCION MP_TTaH
%Matriz de transformación
[T1,Ta1,h1]=MP_TTah(8,90,0,0,[0 0 0 0]);
[T2,Ta2,h2]=MP_TTah(sqrt(32),45,0,0,[0 0 0 0]);
[T3,Ta3,h3]=MP_TTah(sqrt(32),315,0,0,[0 0 0 0]);
[T4,Ta4,h4]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T5,Ta5,h5]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T6,Ta6,h6]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T7,Ta7,h7]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T8,Ta8,h8]=MP_TTah(4,90,0,0,[0 0 0 0]);
[T9,Ta9,h9]=MP_TTah(4,90,0,0,[0 0 0 0]);
[Tp]=TP_Corte(4,4); % VA HASTA LOS EJES
% RIGIDECES EN COORDENADAS DE NUDOS
K1=T1*k1*T1';
K2=T2*k2*T2';
K3=T3*k2*T3';
K4=T4*k4*T4';
K5=T5*k5*T5';
K6=T6*k6*T6';
K7=T7*k7*T7';
K8=T8*k8*T8';
K9=T9*k9*T9';
Kp=Tp*kp*Tp';
% B matriz que tiene P y delta R ( PRIMERA PARTE CARGAS SEGUNDO DES)
S=zeros(ngl);
B=zeros(nglt,1);
% Vectores de esnsamble
lee1=[17 18 19 1 2 3];
lee2=[1 2 3 4 5 6];
lee3=[4 5 6 7 8 9];
lee4=[7 8 9 7 10 11];
lee5=[12 13 14 12 15 16];
lee6=[12 13 14 7 8 9];
lee7=[12 15 16 7 10 11];
lee8=[20 21 22 12 13 14];
lee9=[23 24 25 12 15 16];
leep=[8 7 10 13 12 15];
% FUNCION Ens2
[S,B]=Ens2(lee1,K1,ngl,S,B);
[S,B]=Ens2(lee2,K2,ngl,S,B);
[S,B]=Ens2(lee3,K3,ngl,S,B);
[S,B]=Ens2(lee4,K4,ngl,S,B);
[S,B]=Ens2(lee5,K5,ngl,S,B);
[S,B]=Ens2(lee6,K6,ngl,S,B);
[S,B]=Ens2(lee7,K7,ngl,S,B);
[S,B]=Ens2(lee8,K8,ngl,S,B);
[S,B]=Ens2(lee9,K9,ngl,S,B);
[S,B]=Ens2(leep,Kp,ngl,S,B);
% SOLUCION DEL SISTEMA DE ECUACIONES
B(1)=10; B(5)=-5; B(7)=5; B(12)=5;
fprintf ('DESPLAZAMIENTOS EN LOS GDL\n ')
format long;Delta=S\B(1:ngl)
%ACCIONES FINALES
% OJO SE COMPLETA LOS DELTA A PARTIR DEL NGL+1
Delta(nglt)=0;
D1=Delta(lee1);d1=T1'*D1;a1=k1*d1;
D2=Delta(lee2);d2=T2'*D2;a2=k2*d2;
D3=Delta(lee3);d3=T3'*D3;a3=k3*d3;
D4=Delta(lee4);d4=T4'*D4;a4=k4*d4;
D5=Delta(lee5);d5=T5'*D5;a5=k5*d5;
D6=Delta(lee6);d6=T6'*D6;a6=k6*d6;
D7=Delta(lee7);d7=T7'*D7;a7=k7*d7;
D8=Delta(lee8);d8=T8'*D8;a8=k8*d8;
D9=Delta(lee9);d9=T9'*D9;a9=k9*d9;
Dp=Delta(leep);Gama=Tp'*Dp;Tao=Gm*Gama;
fprintf ('ACCIONES DE MIEMBROS EN SISTEMA REDUCIDO ')
format short; a=[a1 a2 a3 a4 a5 a6 a7 a8 a9]
Tao

RESULTADOS DE LA CORRIDA:
PANEL CORTE.m

PORTICO CON CARGAS PUNTUALES Y UN PANEL A CORTE

ESCUELA POLITECNICA NACIONAL

MAESTRIA EN ESTRUCTURAS

Por Diego Quizanga

DESPLAZAMIENTOS EN LOS GDL

Delta =

0.043256178502158

-0.000061665371246

-0.001285678997370

0.032794584772066

0.010197818209350

0.000638111522933

0.022363423735432

-0.000005740567420

-0.001261914512079

-0.000091976092042

0.000206673226296

0.018588604057513

0.000014950877956
 -0.002169183205757

-0.000063809207688

-0.002455374299348

ACCIONES DE MIEMBROS EN SISTEMA REDUCIDO

a=

Columns 1 through 8

12.4089 -11.3334 9.0574 -3.8028 -11.0419 -4.5247 -3.1083 19.1959

11.3334 -9.0574 6.8095 -1.3456 -11.5207 -3.0067 1.3456 15.5666

-1.9345 -6.3404 -7.1401 -5.4106 -4.9416 -1.2982 -1.7673 0.9381

Column 9

18.7372

14.6291

-4.0035

Tao =

29.5386