Quispe AJC-SD
Quispe AJC-SD
AUTOR:
ASESOR:
LÍNEA DE INVESTIGACIÓN:
Diseño Sísmico y Estructural
LIMA – PERÚ
2022
Dedicatoria
ii
Agradecimientos
A Jesucristo por su
amor y su misericordia.
iii
Índice de contenidos
Dedicatoria .............................................................................................................. ii
Agradecimientos..................................................................................................... iii
Abstract .................................................................................................................. ix
I. INTRODUCCION ............................................................................................... 10
V. DISCUSIÓN ...................................................................................................... 79
REFERENCIAS..................................................................................................... 85
ANEXOS ............................................................................................................... 90
iv
Índice de tablas
Tabla 3. Detalle de los elementos de confinamiento, muro con ladrillo tipo III. .... 25
Tabla 4. Detalle de los elementos de confinamiento muro con ladrillo tipo IV. ..... 26
Tabla 11. Propiedades mecánicas del puntal con ladrillo tipo III. ......................... 54
Tabla 12. Propiedades mecánicas del puntal con ladrillo tipo IV.......................... 55
Tabla 13. Capacidad a compresión del puntal del muro N°1, ladrillo tipo V. ........ 56
Tabla 14. Propiedades al corte de los muros con ladrillo tipo III. ......................... 57
Tabla 15. Propiedades al corte de los muros con ladrillo tipo IV. ......................... 58
Tabla 16. Propiedades al corte de los muros con ladrillo tipo V. .......................... 60
Tabla 18. Propiedades a flexión de los elementos del muro N°1. ........................ 67
Tabla 19. Propiedades a flexión de los elementos del muro N°2 y 3. .................. 68
Tabla 21. Propiedades a compresión del puntal del muro N°1. ........................... 71
Tabla 22. Propiedades a compresión del puntal del muro N°2 y 3. ...................... 72
v
Índice de figuras
Figura 5. Ladrillo King Kong de 18 huecos Tipo III Unidad hueca. ....................... 24
Figura 10. Envolvente de comportamiento histerético muro ladrillo Tipo IV. ........ 28
Figura 24. Diagrama momento curvatura de las columnetas del muro N°1. ........ 48
Figura 25. Diagrama momento curvatura de la viga solera del muro N°1. ........... 49
vi
Figura 27. Diagrama momento curvatura de la viga solera del muro N°2 y 3. ..... 50
Figura 28. Curva de capacidad del pórtico de confinamiento del muro N°1. ........ 51
Figura 29. Curva de capacidad del pórtico de confinamiento, muro N°2 y 3. ....... 52
Figura 30. Diagrama carga – deformación del puntal con ladrillo tipo III. ............. 53
Figura 31. Diagrama carga – deformación del puntal con ladrillo tipo IV. ............ 54
Figura 32. Diagrama carga deformación del muro N°1 con ladrillo tipo V. ........... 55
Figura 33. Curva de capacidad de los muros con ladrillo tipo III. ......................... 57
Figura 34. Curva de capacidad de los muros con ladrillo tipo III. ......................... 58
Figura 35. Curva de capacidad del muro N°1 con ladrillo tipo V. ......................... 59
Figura 36. Comparación de la capacidad del muro N°1 con los antecedentes. ... 60
Figura 43. Momento curvatura de los elementos de confinamiento muro N°1. .... 67
Figura 46. Diagrama carga – deformación del puntal del muro N°1. .................... 71
Figura 47. Diagrama carga – deformación del puntal del muro N°2 y 3. .............. 72
Figura 48. Capacidad resistente del muro N°1 con ladrillos tipo III, IV y V. .......... 73
Figura 42. Capacidad resistente del muro N°2 con ladrillos tipo III, IV y V. .......... 75
Figura 43. Capacidad resistente del muro 3 con ladrillos tipo III, IV y V. .............. 76
vii
Resumen
viii
Abstract
ix
I. INTRODUCCIÓN
11
equivalente de los muros de albañilería confinada con ladrillos industriales, San
Ramón Junín 2021?.
12
Por lo cual la hipótesis general se plantea como tal: La capacidad resistente de
los muros de albañilería confinada con ladrillos industriales, es mayor de 5 T, en
San Ramón Junín 2021. Y como hipótesis específicas tenemos: La capacidad a
flexión de los elementos de confinamiento de los muros de albañilería confinada
con ladrillos industriales, es mayor a 50Tm en San Ramón Junín 2021, La
cortante de fluencia de los pórticos de confinamiento de los muros de albañilería
confinada con ladrillos industriales, es mayor a 10T en San Ramón Junín 2021,
y La capacidad a compresión del puntal equivalente de los muros de albañilería
confinada con ladrillos industriales, es mayor a 10T en San Ramón Junín 2021.
13
II. MARCO TEÓRICO
14
la máxima distorsión registrada en las pruebas (0.6% en la mayoría de los casos
aquí estudiados) para garantizar una resistencia residual del muro.
15
básicos y tendiendo a minimizar el daño estructural. Teniendo en cuenta los
resultados experimentales y analíticos para muros en contacto y aislados en el
plano, se establece que los muros en pleno contacto tendrán un daño leve para
una deriva de 0.20%, moderado para 0.35% y colapsarán para el 1.05 % de
deriva; y los muros aislados tendrán un daño leve para una deriva de 0.70%,
moderado para 1.20% y el 1.88 % de deriva. Estos resultados nos indican que
el sistema de muros aislados puede resultar eficiente ante cargas sísmica
actuando en el plano, siempre que se cumplan todas las especificaciones
constructivas del sistema; algo que en la práctica convencional no se cumple a
cabalidad. Con los niveles de daño obtenidos, se determina que la normativa
vigente está cumpliendo el criterio de diseño (para una deriva de 0.5%) de admitir
la posibilidad de daño, evitando las pérdidas de vidas y minimizando el daño
estructural.
16
En el ámbito nacional; (10), realizo la tesis de titulación: “Desarrollo de sistema
de computación para el cálculo de diagramas momento curvatura de secciones
de concreto armado en Huancayo – 2021”, que tuvo como objetivo principal
evaluar de qué manera la aplicación del Sistema de Computación influye en el
cálculo de diagramas momento curvatura de secciones de concreto armado. Fue
un estudio de tipo aplicada de diseño no experimental. Se tuvo como resultados
que la rigidez a flexión aumenta, respecto al estudio de referencia, de 11847 a
13167 T.m2 en un 11%; la rigidez de post fluencia aumenta de 2 a 129 T.m2 en
5076%, el porcentaje de rigidez efectiva aumenta de 0.19 a 0.21 en un 11% y la
ductilidad local por curvatura disminuye de 27.48 a 22.43 en un -18%. Se
concluye que la aplicación de un sistema de computación facilitará el cálculo de
diagramas momento curvatura de secciones de concreto armado.
17
estático de la misma edificación bajo la documentación FEMA 356, ATC 40 y
ASCE/SEI 31-03. La investigación concluye que el análisis no lineal estático fue
determinante en el cálculo de la capacidad estructural de la edificación de
albañilería confinada al que se investigó, ya que mediante ello se pudo
determinar el real comportamiento de la estructura y según los resultados
obtenidos la estructura tiene una capacidad para una fuerza cortante de Vxx =
175.75 tn un desplazamiento de 2.74 cm y Vyy = 153.19 tn un desplazamiento
de 2.90 cm, donde podemos determinar que las estructuras de albañilería
confinada poseen poca ductilidad a comparación de otros sistemas estructurales
como el caso de pórticos de concreto armado.
18
diseño y el programa Visual Basic Excel. La muestra fueron las secciones de
vigas y columnas de concreto armado, tiene como resultado los diagramas
momento curvatura con diversos modelos constitutivos del material. La
investigación concluye que la interfaz del programa permite la construcción de
diagramas de momento curvatura de manera sencilla y de licencia libre. La hoja
de cálculo es accesible al usuario de manera clara y su precisión ha sido
confirmada con el programa Section Designer y estudios analíticos por varios
autores. La hoja de cálculo tiene sus limitaciones, ya que es imposible realizar el
análisis de todas las posibles secciones de vigas.
En lo teórico tenemos:
• Muros reforzados.
Entre los muros reforzados tenemos los muros de albañilería confinada. Según
(16), estos muros de albañilería están confinados perimetralmente por una
19
cadena de elementos de concreto armado que son las vigas y columnas. Estas
vigas y columnas forman lo que se conoce como pórtico de confinamiento
Albañilería estructural
Pórticos de confinamiento
Columnas
Vigas soleras
LADRILLOS INDUSTRIALES
20
Los ladrillos industriales son unidades albañilería de arcilla cocida, fabricados de
manera industrial a máquina. Estas unidades se clasifican de acuerdo al área de
sus orificios en: King Kong sólido, King Kong de 18 huecos, Pandereta y
Caravista.
Unidades solidas
Son aquellas que pueden o no tener orificios. Si tienen orificios estos tienen que
ser necesariamente perpendiculares a la cara de asiento y el área que ocupan
no debe ser mayor del 30% del área bruta de la cara de asiento. Según (4), se
utilizan para construir muros portantes y no portantes y en el análisis estructural
no se consideran los orificios.
Unidades huecas
21
Figura 3. Unidad de albañilería hueca.
Fuente: Elaboración propia
Unidades tubulares
A diferencia de los casos anteriores, es estas unidades los orificios son paralelos
a la cara del asiento. El área que ocupan estos orificios no está normada. Según
(17), las unidades tubulares solamente se deben utilizar para la construcción de
muros no portantes. Dentro de estas unidades están los ladrillos “Pandereta” que
incorrectamente se utilizan para la construcción de muros portantes.
Unidades alveolares
22
Los ladrillos se clasifican según su calidad en 5 tipos según la variabilidad
dimensional en su ancho alto y largo, alabeo, absorción, porcentaje de vacíos y
su en especial su resistencia a compresión. Aunque la norma solo establece
valores mínimos y máximos, (16) estableció una descripción detallada de la
aplicación de cada tipo de ladrillo en la práctica como se muestra a continuación.
23
Tipo V: Estos ladrillos pueden utilizarse en construcciones en condiciones
exigentes de servicio, en contacto con agua, suelo y lluvias intensas. Tienen muy
alta durabilidad y resistencia.
Características del concreto y acero: Resistencia del concreto: f’c = 175 Kg/cm2,
resistencia de fluencia del acero: fy = 4200 Kg/cm2, Resistencia máxima del
acero: fu = 6330 Kg/cm2.
Dimensiones 24 x 13 x 9 cm
24
Figura 6. Muro de albañilería confinada con ladrillo tipo III.
Fuente: Coral, (5).
Tabla 3. Detalle de los elementos de confinamiento, muro con ladrillo tipo III.
ELEMENTO DIMENSIÓN LONGITUDINAL TRRANSVERSAL
Viga de cimentación 0.30 x 0.30m 4ø5/8” + 2ø1/2” [] 3/8” 1@5, 5@10, r@25
25
Figura 8. Envolvente de comportamiento histerético.
Fuente: Coral, (5).
Características del concreto y acero: Resistencia del concreto: f’c = 175 Kg/cm2,
Resistencia de fluencia del acero: fy = 4200 Kg/cm2, Resistencia máxima del
acero: fu = 6330 Kg/cm2.
Tabla 4. Detalle de los elementos de confinamiento muro con ladrillo tipo IV.
CONFINAMIENTO DIMENSIÓN LONGITUDINAL TRRANSVERSAL
26
Dimensiones 23 x 12.5 x 9 cm
27
Figura 10. Envolvente de comportamiento histerético muro ladrillo Tipo IV.
Fuente: Paredes, y otros, (18).
Características del concreto y acero: Resistencia del concreto: f’c = 175 Kg/cm2,
Resistencia de fluencia del acero: fy = 4200 Kg/cm2, Resistencia máxima del
acero: fu = 6330 Kg/cm2.
Tipo Solido
Fabricante REX
Porcentaje de huecos 32
28
Figura 11. Muro confinado con ladrillo Tipo V.
Fuente: Angles, (19 pág. 37).
29
Figura 13. Ladrillo de arcilla King Kong de 18 huecos.
Fuente: Angles, (19 pág. 4).
30
deformación de fluencia del acero, (𝜀𝑠ℎ ) es la deformación al inicio de la zona de
endurecimiento y (𝜀𝑠𝑢 ) es la deformación de rotura del acero.
𝑓𝑠 = 𝐸𝑠 𝜀𝑠 𝜀𝑠 ≤ 𝜀𝑦 𝐸𝑐. 1
𝑓𝑠 = 𝑓𝑦 𝜀𝑦 < 𝜀𝑠 ≤ 𝜀𝑠ℎ 𝐸𝑐. 2
𝑓𝑠 = 𝑓𝑦 + 𝐸𝑠ℎ (𝜀𝑠 − 𝜀𝑠ℎ ) 𝜀𝑠 > 𝜀𝑠ℎ 𝐸𝑐. 3
𝑓𝑦
𝐸𝑠 = 𝐸𝑐. 4
𝜀𝑦
𝑓𝑠𝑢 − 𝑓𝑦
𝐸𝑠ℎ = 𝐸𝑐. 5
𝜀𝑠𝑢 − 𝜀𝑠ℎ
31
Modelo de Park et al. para el concreto
El modelo de Park et al (1982) es una versión mejorada del Modelo de Kent y
Park (1971) y contempla el incremento de resistencia y ductilidad por el acero
transversal de confinamiento estribos con ayuda de un facto “k” (20). El modelo
de Park consta de las ecuaciones 6-13.
2 𝜀𝑐 𝜀𝑐 2
𝑓𝑐 = 𝑘 𝑓´𝑐 [ − ( ) ] 𝜀𝑐 ≤ 𝑘 𝜀𝑜 𝐸𝑐. 6
𝑘 𝜀𝑜 𝑘 𝜀𝑜
𝑓𝑐 = 𝑘 𝑓´𝑐 [1 − 𝑍(𝜀𝑐 − 𝑘 𝜀𝑜 )] 𝑘 𝜀𝑜 ≤ 𝜀𝑐 ≤ 𝜀𝑢 𝐸𝑐. 7
0,5
𝑍= 𝐸𝑐. 8
𝜀50𝑢 + 𝜀50ℎ + 𝑘 𝜀𝑜
3 + 0,29 𝑓´𝑐
𝜀50𝑢 = 𝐸𝑐. 9
145 𝑓´𝑐 − 1000
3 𝑏"
𝜀50ℎ = 𝜌𝑣 √ 𝐸𝑐. 10
4 𝑠
0,8
𝜀𝑢 = 𝑘 𝜀𝑜 + 𝐸𝑐. 11
𝑍
𝜌𝑣 𝑓𝑦𝑣
𝑘 = 1+ 𝐸𝑐. 12
𝑓´𝑐
2 (𝑏 ′ + ℎ′ )𝐴𝑠
𝜌𝑣 = 𝐸𝑐. 13
𝑏 ′′ ℎ′′ 𝑠
Donde:
32
𝑏 " es el ancho de la sección transversal del elemento, pero del concreto
confinado.
ℎ′′ es la distancia entre bordes externos del estribo en la dirección del peralte h.
33
Figura 17. Porción diferencial de un elemento sometido a flexión.
Fuente: Chang, 2015.
Pórticos de confinamiento
La capacidad de la estructura es representada por la curva Push Over que
relaciona la fuerza cortante en la base con el desplazamiento en el tope del
edificio. Se usa un patrón de distribución de fuerza lateral en la altura del edificio
que representa el primer modo de vibración. A medida que las cargas se
incrementan, el modelo matemático de la estructura se modifica para simular la
degradación de la rigidez de los elementos que empieza con la fluencia del
acero. Según (20), el incremento de cargas es continuado hasta que la estructura
alcance los límites preestablecidos.
34
Muros de albañilería confinada.
El traslado de las fuerzas laterales a través de los pórticos con albañilería
ocasiona una distribución de esfuerzos no uniforme al interior de la albañilería y
los elementos del pórtico (21). Como las fuerzas laterales son incrementadas, la
distribución de esfuerzos varia hasta que ocurre la falla de la albañilería. Cuando
la albañilería alcanza su resistencia al cortante o a la compresión ocurre la falla.
Procedimiento general
Usando este procedimiento general, se puede predecir la capacidad en el plano
de los pórticos con albañilería (21).
35
usando la Ecuación 16. Dónde: (ared) es el ancho del puntal equivalente
reducido por daño o aberturas, (a) es el ancho del puntal equivalente, (R1i) es el
factor de reducción debido a presencia de aberturas y (R2i) es el factor de
reducción debido al daño del relleno existente. Las propiedades del material que
deben ser asignadas al puntal consisten en la capacidad a compresión del puntal
(Rpuntal) y el módulo de elasticidad de la albañilería (Em), donde Rpuntal se
calcula usando la Ecuación 18. Dónde: (Rcomp) es la resistencia a compresión
de la albañilería simple, (Rcort) es la resistencia al cortante ó a la adherencia de
la albañilería simple, (tefec) es el espesor neto de la albañilería simple, (f’m) es
la resistencia a compresión de la albañilería simple, (tmort) es el espesor
transversal del mortero y (f’v) es la menor entre las resistencias al cortante y a la
adherencia mortero - ladrillo de la albañilería.
𝑎
𝐼𝑐𝑜𝑙 = 𝐸𝑐. 14
𝐶𝑜𝑠𝜃𝑐𝑜𝑙
ℎ − 𝐼𝑐𝑜𝑙
𝜃𝑐𝑜𝑙 = 𝐴𝑟𝑐𝑡𝑎𝑛 ( ) 𝐸𝑐. 15
𝑙
𝑎𝑝𝑢𝑛𝑡𝑎𝑙 = 3 𝑎𝑟𝑒𝑑 𝐸𝑐. 16
𝑎𝑟𝑒𝑑 = 𝑎 𝑅1𝑖 𝑅2𝑖 𝐸𝑐. 17
𝑅𝑐𝑜𝑚𝑝
𝑅𝑝𝑢𝑛𝑡𝑎𝑙 = 𝐸𝑙 𝑚𝑒𝑛𝑜𝑟 { 𝑅𝑐𝑜𝑟𝑡 } 𝐸𝑐. 18
𝐶𝑜𝑠𝜃𝑝𝑢𝑛𝑡𝑎𝑙
𝑅𝑐𝑜𝑚𝑝 = 𝑎𝑟𝑒𝑑 𝑡𝑒𝑓𝑒𝑐 𝑓′𝑚 𝐸𝑐. 19
𝑎 = 0.175 𝐷 (𝜆1 𝐻)−0.4 𝐸𝑐. 20
𝐸𝑚 𝑡 𝑆𝑒𝑛2𝜃 1/4
𝜆1 = ( ) 𝐸𝑐. 21
4 𝐸𝑐 𝐼𝑐𝑜𝑙 ℎ
𝑅𝑐𝑜𝑟𝑡 = 𝑡𝑚𝑜𝑟𝑡 𝑙 𝑓′𝑣 𝑅1𝑖 𝑅2𝑖 𝐸𝑐. 22
ℎ − 2 𝐼𝑐𝑜𝑙
𝜃𝑝𝑢𝑛𝑡𝑎𝑙 = 𝐴𝑟𝑐𝑡𝑎𝑛 ( ) 𝐸𝑐. 23
𝑙
36
Figura 19. Diagrama carga – deformación del puntal equivalente
Fuente: Carrillo, (21)
Asignar zonas rígidas en las uniones del pórtico para representar la intersección
real entre elementos viga y columna. Las zonas rígidas deben modelarse con un
factor de rigidez de 0.5, es decir la mitad de la zona rígida se considera efectiva.
37
inercia en el diseño sísmico. Las distribuciones de fuerza de inercia
recomendadas están dadas en la Sección 3.3.3.2 del (22) y del (23).
Rigidez elástica
La rigidez elástica es la propiedad mecánica del muro de albañilería cuando está
libre de agrietamiento y corresponde al rango elástico de respuesta del muro
(24). Es la pendiente de la línea inicial con mayor inclinación del modelo Trilineal
de la curva de capacidad.
Rigidez inelástica
La rigidez inelástica es la propiedad mecánica del muro de albañilería cuando
inicia el agrietamiento visible del muro y corresponde al nivel de servicio donde
se produce el término del rango elástico de respuesta del muro (24). Es la
pendiente de la segunda línea del modelo Trilineal de la curva de capacidad.
38
III. METODOLOGÍA
39
La curva de capacidad resistente se definirá operacionalmente en 3
dimensiones que son: La capacidad a flexión de los elementos de
confinamiento, la cortante de fluencia de los pórticos de confinamiento y la
Capacidad a compresión del puntal equivalente.
40
Según (28) estas fichas sirven para registrar información asociados a las
variables que está siendo estudiados por el investigador.
41
otra vez y producirá los mismos resultados (28). Por esta razón, los
instrumentos estándar mencionados representan especificaciones confiables
para el análisis sísmico de estructuras.
3.5. Procedimiento
42
Figura 21. Muro confinado N°2.
Fuente: Elaboración propia.
43
44
45
3.6. Método de análisis de datos
46
Tabla 10. Procedimiento de la investigación.
47
Paso 3. Construcción de diagramas momento – curvatura con el programa
“d_curvatura”, Mostrado en las figuras 15 y 16
Muro N°1
Figura 24. Diagrama momento curvatura de las columnetas del muro N°1.
Fuente: Elaboración propia.
Interpretación: En la figura 24, se muestra el diagrama momento curvatura
de las columnetas del muro de albañilería N°1. Se observa que la capacidad
a flexión (Mu) es de 2.90 Tm. La primera pendiente representa la rigidez
elástica a flexón (EIo) que es igual a 136.82 T/m, la segunda pendiente es la
rigidez inelástica a flexión (EIp) que es igual a 0.34 T/m y la ductilidad por
curvatura (u) es de 7.77. También se puede apreciar una degradación
sustancial de la rigidez debido a la fluencia del acero de refuerzo, la tercera
es una caída súbita de la resistencia que se equilibra en una resistencia
residual de 0.6 Tm. En esta última etapa se da la formación de una
articulación con rigidez nula.
48
Figura 25. Diagrama momento curvatura de la viga solera del muro N°1.
Fuente: Elaboración propia.
Interpretación: En la figura 25, se muestra el diagrama momento curvatura
de las vigas soleras del muro de albañilería N°1. Se observa que la
capacidad a flexión (Mu) es de 2.80 Tm. La primera recta representa la
rigidez elástica a flexión (EIo) que es igual a 136.82 T/m, la segunda es la
rigidez inelástica a flexión (EIp) que es igual a 0.56 T/m y la ductilidad por
curvatura (u) es de 14.77. También se puede apreciar una degradación
sustancial de la rigidez de la sección, debido a la fluencia del acero de
refuerzo, la tercera es una caída súbita de la resistencia que se equilibra en
una resistencia residual de 0.6 Tm. En esta última etapa se forma una
articulación con rigidez nula.
Muro N°2 y 3
49
Figura 26. Diagrama momento curvatura, columnetas del muro N° 2 y 3.
Fuente: Elaboración propia.
Interpretación: En la figura 26, se muestra el diagrama momento curvatura
de las columnetas del muro de albañilería N°2 y 3. Se observa una capacidad
a flexión (Mu) de 22.5 Tm, la primera pendiente es la rigidez elástica a flexión
(EIo) que es igual a 2285.71 T/m, la segunda pendiente es la rigidez
inelástica a flexión (EIp) que es igual a 63.11 T/m y la ductilidad por curvatura
(u) es de 15.71. También se puede apreciar una degradación sustancial de
la rigidez debido a la fluencia del acero de refuerzo, la tercera es una caída
súbita de la resistencia que se equilibra en una resistencia residual de 3.5
Tm. En esta última etapa se forma una articulación con rigidez nula.
Figura 27. Diagrama momento curvatura de la viga solera del muro N°2 y 3.
Fuente: Elaboración propia.
Interpretación: En la figura 27, se muestra el diagrama momento curvatura
de la viga solera del muro de albañilería N°2 y 3. Se observa una capacidad
a flexión (Mu) de 16.00 Tm, la primera pendiente es la rigidez elástica a
flexión (EIo) que es igual a 1357.14 T/m, la segunda pendiente es la rigidez
inelástica a flexión (EIp) que es igual a 48.87 T/m, y la ductilidad por
curvatura (u) es de 20.00. También se puede apreciar una degradación
sustancial de la rigidez debido a la fluencia del acero de refuerzo, la tercera
es una caída súbita de la resistencia que se equilibra en una resistencia
50
residual de 2 Tm. En esta última etapa se forma una articulación con rigidez
nula.
Muro N°1.
Figura 28. Curva de capacidad del pórtico de confinamiento del muro N°1.
Fuente: Elaboración propia.
51
Interpretación: En la figura 28, se muestra la curva de capacidad resistente
del pórtico de confinamiento del muro de albañilería N°1 de donde se
obtienen las propiedades al corte del sistema. Se observa una cortante de
fluencia (Vu) de 4.50 T, la primera pendiente es la rigidez elástica al corte
(Ke) que es igual a 150 T/m, la segunda pendiente es la rigidez inelástica al
corte (Ki) que es igual a 0.71 T/m y la ductilidad por corte (ud) es de 5.67.
Mientras más grande sea la ductilidad, la estructura disipara mayor energía
sísmica.
Muro N°2 y 3.
52
8.50. Mientras más grande sea la ductilidad, la estructura disipara mayor
energía sísmica.
Figura 30. Diagrama carga – deformación del puntal con ladrillo tipo III.
Fuente: Elaboración propia.
53
Interpretación: En la figura 30, se observa el diagrama carga deformación
del puntal de los muros con ladrillo industrial tipo III. El puntal del muro N°1
tendrá una capacidad a compresión (Ru) de 10.10 T y el muro N°2 y 3 una
capacidad a compresión (Ru) de 29.14 T, como también se puede observar
en la tabla 11.
Tabla 11. Propiedades mecánicas del puntal con ladrillo tipo III.
Figura 31. Diagrama carga – deformación del puntal con ladrillo tipo IV.
54
Fuente: Elaboración propia.
Interpretación: En la figura 31, se observa el diagrama carga deformación
del puntal de los muros con ladrillo industrial tipo IV. El puntal del muro N°1
tendrá una capacidad a compresión de 9.37 T y el puntal del muro N°2 y 3
una capacidad a compresión de 27.03 T.
Tabla 12. Propiedades mecánicas del puntal con ladrillo tipo IV.
PUNTO A B PROPIEDAD A LA COMPRESIÓN
dy Ry du Ru EAe EAi ud
MURO
(m) (T) (m) (T) (T/m) (T/m) (Adim)
1 0.002 9.5 0.01 9.37 4750.00 0.00 5.00
2Y3 0.004 27.00 0.013 27.03 6750.00 0.00 3.25
Figura 32. Diagrama carga deformación del muro N°1 con ladrillo tipo V.
55
Fuente: Elaboración propia.
Interpretación: En la figura 32, se observa el diagrama carga deformación
del puntal de los muros con ladrillo industrial tipo V. El puntal del muro N°1
tendrá una capacidad a compresión de 12.00 T y el puntal del muro N°2 y 3
una capacidad a compresión de 34.01 T, como también se puede observar
en la tabla 13.
Tabla 13. Capacidad a compresión del puntal del muro N°1, ladrillo tipo V.
PUNTO A B PROPIEDAD A LA COMPRESIÓN
dy Ry du Ru EAe EAi ud
MURO
(m) (T) (m) (T) (T/m) (T/m) (Adim)
1 0.002 12 0.01 12 6000.00 0.00 5.00
2Y3 0.004 34 0.013 34 8500.00 0.00 3.25
56
Figura 33. Curva de capacidad de los muros con ladrillo tipo III.
Fuente: Elaboración propia.
Interpretación: En la figura 33, se observa la curva de capacidad resistente
de los muros de albañilería confinada con ladrillo industrial tipo III. El muro
N°1,2 y 3 tienen una capacidad ultima (Vu) de 17.09, 136.00 y 136.00 T
respectivamente, como también se observa en la tabla 14.
Tabla 14. Propiedades al corte de los muros con ladrillo tipo III.
PUNTO A B PROPIEDAD AL CORTE
dy Vy du Vu Ke Ki ud
MURO
(m) (T) (m) (T) (T/m) (T/m) (Adim)
1 0.0025 16 0.0151 17.09 6400.00 79.37 6.04
2 0.0025 109 0.0145 136.00 43600.00 2250.00 5.80
3 0.0075 110 0.0281 136.00 14666.67 1262.14 3.75
57
inelástica al corte (Ki) de 1262.14 T/m y una ductilidad al desplazamiento
(ud) de 3.75.
Figura 34. Curva de capacidad de los muros con ladrillo tipo IV.
Fuente: Elaboración propia.
Interpretación: En la figura 34, se observa la curva de capacidad resistente
de los muros de albañilería confinada con ladrillo industrial tipo IV. El muro
N°1,2 y 3 tienen una capacidad ultima (Vu) de 15.79, 129.64 y 127.93 T
respectivamente, como también puede observar en la tabla 15.
Tabla 15. Propiedades al corte de los muros con ladrillo tipo IV.
PUNTO A B PROPIEDAD AL CORTE
dy Vy du Vu Ke Ki ud
MURO
(m) (T) (m) (T) (T/m) (T/m) (Adim)
1 0.003 14.23 0.015 15.79 4743.33 130.00 5.00
2 0.003 100.83 0.015 129.64 33610.00 2400.83 5.00
3 0.008 102.2 0.029 127.93 12775.00 1225.24 3.63
58
inelástica al corte (Ki) de 2400.83 T/m y una ductilidad al corte (ud) de 5.00.
El muro N°3 tiene una rigidez elástica al corte (Ke) de 12775.00 T/m, una
rigidez inelástica al corte (Ki) de 1225.24 T/m y una ductilidad al corte (ud)
de 3.63.
Figura 35. Curva de capacidad del muro N°1 con ladrillo tipo V.
Fuente: Elaboración propia.
Interpretación: En la figura 35, se observa la curva de capacidad resistente
de los muros de albañilería confinada con ladrillo industrial tipo V. El muro
N°1,2 y 3 tienen una capacidad ultima (Vu) de 19.42, 154.09 y 153.12 T
respectivamente, como también puede observarse en la tabla 16.
59
Figura 36. Comparación de la capacidad del muro N°1 con los antecedentes.
Fuente: Elaboración propia.
60
3.6.5. Aplicación de la curva de capacidad resistente en la evaluación del
desempeño sísmico de viviendas de albañilería confinada con ladrillos
industriales.
61
Figura 38. Curva de capacidad de los muros de la vivienda.
Fuente: Elaboración propia.
Interpretación: En la figura 38, se muestra la curva de capacidad resistente
de los muros de la vivienda de albañilería confinada con ladrillo industrial tipo
III. El muro de los ejes A, B, C, D, 1, 2 y 3 tienen una capacidad ultima (Vu)
de 66, 33, 33, 66, 29, 29 y 29 respectivamente, como también se muestra en
la tabla 17.
62
Interpretación: En la tabla 17, se muestran los puntos notables de la curva
de capacidad y las propiedades mecánicas al corte. El muro de los ejes A y
D tiene una rigidez elástica al corte (Ke) de 6058 T/m y una rigidez inelástica
al corte (Ki) de 273 T/m. El muro de los ejes B y C tienen una rigidez elástica
al corte (Ke) de 2897 T/m y una rigidez inelástica al corte (Ki) de 174 T/m. El
muro de los ejes 1, 2 y 3 tienen una rigidez elástica al corte (Ke) de 2430 T/m
y una rigidez inelástica al corte (Ki) de 175 T/m. (α) es la dirección del muro
respecto a la coordenada (1) y (r) es su posición respecto al centro de masas.
63
Punto de desempeño
dpi = 0.0073 m
PFR1 = 1.280 -
di = 0.0093 m
Hedif = 8.40 m
Deriva = 0.1113 %
64
Fuente: Elaboración propia.
Interpretación: En la figura 41, se muestra la aplicación de la curva de
capacidad resistente en la evaluación del desempeño sísmico de una
vivienda de albañilería confinada con ladrillo industrial tipo III en la dirección
Y. La curva de capacidad no alcanza el espectro reducido, dando a entender
que la estructura sufrirá el colapso en la dirección Y.
65
IV. RESULTADOS
Ubicación política
Muro N°1
66
Figura 43. Momento curvatura de los elementos de confinamiento muro N°1.
Fuente: Elaboración propia.
Interpretación: En la figura 43, se observa el diagrama momento curvatura
de los elementos de confinamiento del muro de albañilería N°1. Las
columnetas tienen una capacidad a flexión (Mu) de 2.90 Tm y la viga solera
tiene una capacidad a flexión de 2.80 Tm como también se observa en la
tabla 17. Esto significa que, a pesar de tener la misma sección, las tablas
ASCE 41- 13 otorga mayor ductilidad por curvatura a las columnas en un
90%.
67
(ud) de 7.77. La viga solera tienen una rigidez a flexión elástica (EIo) de
136.82 Tm2, una rigidez a flexión inelástica (EIp) de 0.56 T/m y una ductilidad
por desplazamiento (ud) de 14.77.
Muro N°2 y 3
A B PROPIEDAD A FLEXIÓN
PUNTO øy My øu Mu EIo EIp u
(1/m) (Tm) (1/m) (Tm) (Tm2) (Tm2) (Adim)
68
Columneta 1 0.007 16 0.11 22.5 2285.71 63.11 15.71
Columneta 2 0.007 16 0.11 22.5 2285.71 63.11 15.71
Viga solera 0.007 9.5 0.14 16 1357.14 48.87 20.00
Muro N°1,2 y 3
69
Interpretación: En la figura 45, se observa la curva de capacidad resistente
del pórtico de confinamiento de los muros de albañilería. El pórtico del muro
N°1 tiene una cortante de fluencia (My) de 4.50 T y el pórtico del muro N°2 y
3 tiene una cortante de fluencia de 39.00 T como también se puede observar
en la tabla 19. Esto significa que la resistencia de falla de los muros N° 2 y 3
son 8.67 veces mayor a la resistencia de falla del muro N° 1, debido a que
su longitud es 3.33 veces mayor, además de la robustez de los elementos
de confinamiento.
Muro N°1
70
Figura 46. Diagrama carga – deformación del puntal del muro N°1.
Fuente: Elaboración propia.
Interpretación: En la figura 46, se muestra la determinación de la capacidad
a compresión del puntal del muro N°1, se obtuvo como resultado una
capacidad a compresión (Ru) de 10.10, 9.37 y 11.79 T utilizando los ladrillos
industriales tipo III, IV y V. respectivamente. Esto significa que hay una
variación 17 % en la resistencia de rotura del muro de albañilería, debido a
la calidad del ladrillo, que pasa de ser un ladrillo hueco tipo III a un ladrillo
solido tipo V.
71
puntal con ladrillo tipo III, IV y V es distinta con valores de 5000, 4750 y 6000
T/m respectivamente. En cambio, la rigidez axial inelástica (EAi) de cada
puntal con ladrillo tipo III, IV y V son iguales a 0 T/m. Y la ductilidad por
desplazamiento (ud) de cada puntal con ladrillo tipo III, IV y V también son
iguales con un valor de 5.00.
Muro N°2 y 3
Figura 47. Diagrama carga – deformación del puntal del muro N°2 y 3.
Fuente: Elaboración propia.
Interpretación: En la figura 47, se observa el diagrama carga – deformación
del puntal del muro N°2 y 3, se obtuvo como resultado una capacidad a
compresión (Ru) de 29.14, 27.03 y 34.01 T utilizando ladrillos industriales
tipo III, IV y V respectivamente. Esto significa, como en el caso anterior, que
hay una variación 17 % en la resistencia de rotura del muro de albañilería,
debido a la calidad del ladrillo, que pasa de ser un ladrillo hueco tipo III a un
ladrillo solido tipo V.
72
IV 0.004 27 0.013 27.03 6750.00 0.00 3.25
V 0.004 34 0.013 34.01 8500.00 0.00 3.25
Muro N°1
Figura 48. Capacidad resistente del muro N°1 con ladrillos tipo III, IV y V.
Fuente: Elaboración propia.
Interpretación: En la figura 48, se observa la curva de capacidad resistente
de los muros de albañilería tipo N°1. Los ladrillos industriales tipo III, IV y V
dieron como resultado una cortante de fisuramiento (Vy) de 16.00, 14.23 y
73
17.79 T y una capacidad resistente (Vu) de 17.10, 15.79 y 19.42 T
respectivamente. Esto significa una variación de 11 % en la fuerza de corte
que origina el fisuramiento y de 14% en la fuerza de corte que provoca el
colapso parcial del muro de albañilería.
Muro N°2
74
Figura 49. Capacidad resistente del muro N°2 con ladrillos tipo III, IV y V.
Fuente: Elaboración propia.
Interpretación: En la figura 49, se observa la curva de capacidad resistente
de los muros de albañilería N° 2. Los ladrillos tipo III, IV y V dieron como
resultado una cortante de fisuramiento (Vy) de 109.00, 100.83 y 124.83 T y
una capacidad ultima (Vu) de 136.00, 129.64 y 154.09 T respectivamente.
Esto significa una variación de 15 % en la fuerza de corte que origina el
fisuramiento y de 13 % en la fuerza de corte que provoca el colapso parcial
del muro de albañilería.
75
ductilidad al desplazamiento (ud) es de 5.08. Con el ladrillo tipo IV la rigidez
elástica al corte (Ke) es de 33610 T/m, la rigidez inelástica al corte (Ki) es de
2401 T/m y la ductilidad al desplazamiento (ud) de 5.00. Con el ladrillo tipo
V la rigidez elástica al corte (Ke) es de 41610 T/m, la rigidez inelástica al
corte (Ki) es de 2438 T/m y la ductilidad al desplazamiento (ud) es de 5.00.
Muro N°3
Figura 50. Capacidad resistente del muro 3 con ladrillos tipo III, IV y V.
Fuente: Elaboración propia.
Interpretación: En la figura 40, se observa la curva de capacidad resistente
de los muros de albañilería tipo 3. Los ladrillos tipo III, IV y V dieron como
resultado una cortante de fisuramiento (Vy) de 110.00, 102.20 y 126.67 T y
una capacidad ultima (Vu) de 136.00, 127.93 y 153.12 T respectivamente.
Esto significa, como en el caso anterior, una variación de 15 % en la fuerza
de corte que origina el fisuramiento y de 13 % en la fuerza de corte que
provoca el colapso parcial del muro de albañilería.
76
III 0.0075 110.00 0.0281 136 14666.67 1262.14 3.75
IV 0.008 102.20 0.029 127.93 12775.00 1225.24 3.63
V 0.008 126.67 0.027 153.12 15833.75 1392.11 3.38
Límite en
Límite de Desplazamiento
Grado de desempeño desplazamiento (m)
distorsión (%) del tope (m)
Altura Edif. = 8.4 m
77
no corta el espectro del sismo catalogado como raro, por lo tanto, la
estructura alcanza el colapso.
78
V. DISCUSIÓN
Discusión 1.
Discusión 2.
79
estas cortantes se determina el límite de deformación de los muros de albañilería
según la tabla 7-9 del FEMA 256.
Discusión 3.
Como se puede observar, la capacidad a compresión del puntal del muro N°2 y
3 con ladrillos tipo III, IV y V es igual a 29, 27 y 34 T y están por encima del
resultado de (40) en 79, 67 y 110%, debido a que el ladrillo utilizado por Lema
80
es tipo I con una resistencia a compresión (fm) de 30 Kg/cm2. Cabe resaltar que,
en la presente investigación la capacidad a compresión del puntal con ladrillo
tipo IV es menor al valor del puntal con ladrillo tipo III, pero debería ser todo lo
contrario porque el primero es de mayor calidad. Esto revela errores en los
ensayos a compresión de pilas de albañilería con ladrillo tipo III, realizado por
(5), que dan como resultado una resistencia a compresión (fm) de 96.46 Kg/cm2
para una pila con ladrillos que tienen una resistencia a compresión (fb) de 110
kg/cm2. La diferencia entre estos valores (fb-fm) debería ser mayor a 60 kg/cm2
como lo muestran los demás autores.
Discusión 4.
Se puede observar que la capacidad ultima (Vu) del muro N°1 con ladrillo tipo V
es 19.42 T y es menor al resultado de (19) y (11) en un -8%.
81
VI. CONCLUSIONES
Conclusión 1.
Conclusión 2.
Conclusión 3.
82
especifica que la capacidad a compresión del puntal equivalente de los muros
de albañilería confinada con ladrillos industriales, sea mayor a 10T.
Conclusión 4.
83
VII. RECOMENDACIONES
84
REFERENCIAS
construidos con ladrillos King Kong de fabricación industrial. Lima : s.n., 2018.
https://bit.ly/3IszcCc.
lateral - Distorsión en muros de manpostería confinada. 2010. Vol. 83, págs. 25-42.
85
9. CANO LICONA, Jesús y TENA COLUNGA, Arturo. Diseño sísmico de
https://bit.ly/34HxXAV.
https://bit.ly/3taeDWj.
https://bit.ly/3q8714H.
https://repositorio.pucp.edu.pe/index/handle/123456789/145873.
86
17. ABANTO, Flavio. Análisis y diseño de edificaciones de albañilería. Lima :
un muro confinado con ladrillos de concreto y otro con ladrillos de arcilla. Lima :
https://repositorio.uniandes.edu.co/bitstream/handle/1992/10305/u250827.pdf?seq
uence=1.
22. FEMA 273. NEHRP guidelines for the seismic rehabilitation of buildings.
23. FEMA 356. Prestandard and Commentary for the Seismic Rehabilitation of
https://www.redalyc.org/pdf/618/61807003.pdf.
87
hormigon armado. Empleo del programa matlab en aplicaciones estructurales a
2019.
31. NTP E060. Norma técnica de concreto armado. Lima : Macro, 2009.
32. NEHRP NIST. Soil Structure Interaction for Building Structures. 2012.
33. ASCE/SEI 7-10. Minimum design loads for buildings and other structures.
34. ASCE/SEI 41. American Society of Civil Engineers. seismic evaluation and
36. ATC - 40. Applied Technology Council. Seismic evaluation and retrofit of
88
38. OSORIO, Juan y RAMÍREZ, José. Influencia de los muros de manpostería
https://bit.ly/3sV7Puf.
https://bit.ly/3MahLJh.
89
ANEXOS
90
Anexo N° 1. Matriz de Operacionalización de Variables
TÍTULO: GENERACIÓN DE CURVAS DE CAPACIDAD RESISTENTE DE MUROS DE ALBAÑILERÍA CONFINADA CON LADRILLOS INDUSTRIALES, SAN RAMÓN JUNÍN 2021.
DEFINICIÓN DEFINICIÓN
VARIABLE DIMENSIONES INDICADORES
CONCEPTUAL OPERACIONAL
92
I3: Ductilidad por
corte
93
Anexo N° 3. Programas de computación
a_nudos
function [np,nt,nn,ne,nc,nv,gl,glr,ij]=a_nudos(np,nt)
%--------------------------------------------------------------------------
% El programa "nudos1", asigna los nudos i-j a cada elemento de un sistema
% aporticado con elementos en posición horizontal o vertical.
%
% Por: Juan carlos Quispe Ango
% UNIVERSIDAD CESAR VALLEJO
% Febrero/2022
%--------------------------------------------------------------------------
% [np,nt,nn,ne,nc,nv,gl,glr,ij]=nudos1(np,nt)
%--------------------------------------------------------------------------
% ENTRADA
% np = Número de pisos
% nt = Número de tramos
% RESULTADOS
% np = Número de pisos
% nt = Número de tramos
% nn = Número de nudos
% ne = Número de elementos
% nc = Número de columnas
% nv = Número de vigas
% gl = Número de coordenadas
% glr = Número de coordendas restringidas
% ij = Matriz de nudos para cada elemento
%--------------------------------------------------------------------------
clc
% 1. ENUMERACIÓN DE NODOS COLUMNAS
nn=(nt+1)*(np+1);
nc=(nt+1)*np;
nv=nc-np;
ne=nc+nv;
gl=nn*3;
glr=(nt+1)*3;
gll=gl-glr;
% 2. ENUMERACIÓN DE NODOS
% Columnas
col=1:1:nc;
ic=col;
jc=ic+(nt+1);
% Vigas
vig=1:1:nv;
pisos=zeros(1,nv);
for i=1:nv
j=ceil(i/nt);
pisos(1,(j-1)*nt+1:j*nt)=j;
end
iv=nt+vig+pisos;
jv=iv+1;
% Total
ij=[ic',jc';iv',jv'];
a_puntal
function [ij,nd,ne]=a_puntal(ij,np,nt,ne)
%--------------------------------------------------------------------------
% El programa "a_puntal", incoorpora en las ultimas filas la matriz de
% nudos i-j del sistema aporticado, los nudos i-j de los puntales
% equivalentes.
%
% Por: Juan Carlos Quispe Ango
% UNIVERSIDAD CESAR VALLEJO
% Febrero/2021
%--------------------------------------------------------------------------
% [ij,nd,ne]=a_puntal(ij,np,nt,ne)
%--------------------------------------------------------------------------
% ENTRADA
% ij = Matriz de nudos para cada elemento
% np = Número de pisos
% nt = Número de tramos
% ne = Número de elementos
% RESULTADOS
% ij = Matriz de nudos para cada elemento
% nd = Número de diagonales
% ne = Número de elementos
%--------------------------------------------------------------------------
% 1. MATRIZ DE NUDOS PARA DIAGONALES (xy)
for n=1:np
x=(n-1)*nt+n+1:n*(nt+1);
y=n*(nt+1)+1:n*(nt+1)+nt;
xy((n-1)*nt+1:n*nt,1:2)=[x',y'];
end
% 2. INCOORPORACIÓN DE NUDOS DE LAS DIAGONALES (xy) A (ij)
nd=numel(xy(:,1));
ij(ne+1:ne+nd,1:2)=xy;
ne=ne+nd;
% 3. SALIDA
fprintf('INCOORPORACIÓN DE NUDOS DE LAS DIAGONALES A MATRIZ (ij) \n');
b_colocacion
function [xyz,VC]=b_colocacion(ne,nn,ij)
%--------------------------------------------------------------------------
% El programa "colocacion", construye el vector de colocación (VC) que
% sirve para el ensamblaje de la matriz de rigidez global de la estructura.
%
% Por: Juan Carlos Quispe Ango
% UNIVERSIDAD CESAR VALLEJO
% Febrero/2021
%--------------------------------------------------------------------------
% [xyz,VC]=colocacion2(ne,nn,ij)
%--------------------------------------------------------------------------
% ENTRADA
% ne = Número de elementos
% nn = Número de nudos
% ij = matriz de nudos para cada elemento
% RESULTADOS
% xyz = Coordenadas (gl) Horizontal, vertical y Momento para cada nudo
% VC = Vector de colocación
95
%--------------------------------------------------------------------------
% 1. CONSTRUCCIÓN DE COORDENADAS PARA NUDOS
xyz=[];
xyz(:,1)=1:3:nn*3-2;
xyz(:,2)=2:3:nn*3-1;
xyz(:,3)=3:3:nn*3;
% 2. ORDENAMIENTO DE COORDENADAS PARA CADA ELEMENTO
for i=1:ne
ei(i,:)=xyz(ij(i,1),[1,2,3]); % Recoge la fila ij(i,1) de xyz
ej(i,:)=xyz(ij(i,2),[1,2,3]); % Recoge la fila ij(i,2) de xyz
end
VC=[ei,ej];
% 3. SALIDA
fprintf('Grados de libertad en los nudos (xyz): \n');
fprintf('Vector de colocación (VC): \n');
c_cargas
function [x,c,m]=c_cargas(np,nt,ne,nc,nv,gl,glr,VC)
%--------------------------------------------------------------------------
% El programa "cargas3",
%
% Por: Juan carlos Quispe Ango
% UNIVERSIDAD CESAR VALLEJO
% Febrero/2022
%--------------------------------------------------------------------------
% [x,c,m]=cargas3(np,nt,ne,nc,nv,gl,glr,VC)
%--------------------------------------------------------------------------
% ENTRADA
% np = Número de pisos
% nt = Número de tramos
% ne = Número de elementos
% nc = Número de columnas
% nv = Número de vigas
% gl = Número de coordenadas o grados de libertad
% glr = Número de coordenadas restringidas
% VC = Vector de colocación
% RESULTADOS
% x = Matriz con F. AXIALES en los extremos de los elementos
% c = Matriz con CORTANTES en los extremos de los elemento
% m = Matriz con MOMENTOS en los extremos de los elementos
%--------------------------------------------------------------------------
% 1. ENTRADA
datos=input('Ingrese el archivo de datos: ');
L=xlsread(datos,'B3:B225');
c1=xlsread(datos,'C3:C225');
c2=xlsread(datos,'D3:D225');
an=xlsread(datos,'E3:E225');
EA=xlsread(datos,'L3:L225');
GA=xlsread(datos,'M3:M225');
EI=xlsread(datos,'N3:N225');
CARGAS=xlsread(datos,'O3:O225');
T=xlsread(datos,'P3:P225');
EIo=EI; % Rigidez inicial elasticos en los extremos
EIab=[EI,EI];
% 2. VECTOR DE CARGAS DE EMPOTRAMIENTO (R)
R=zeros(gl,1);
96
for i=1:ne
if i<=nc % Cargas de empotramiento en las columnas
F(:,i)=zeros(6,1);
else % Cargas de empotramiento en las vigas
w=CARGAS(i);
l=L(i);
F(:,i)=[0 w*l/2 w*l^2/12 0 w*l/2 -w*l^2/12]';
vc=VC(i,:);
iR=zeros(gl,1);
iR(vc)=F(:,i); % Se remplazo por las 7 instrucciones siguientes
% 06/01/22
% vc=VC(i,:);
% vci=vc;
% vc(vc==0)=[]; % Se eliminan las coordenadas iguales a (0).
% a=find(vci); % Se ubica los casilleros de las coordenadas...
% f=F(:,i); % distintas de (0).
% n=f(a);
% iR(vc)=n;
R=R+iR;
end
end
fprintf('Vector de CARGAS DE EMPOTRAMIENTO (R) \n');
R
% 3. VECTOR DE CARGAS DE GRAVEDAD EN LOS NUDOS (Q)
Q=zeros(gl,1);
for i=1:nc
c=2*i+(i-1);
Q(c)=-CARGAS(i);
end
fprintf('Vector de CARGAS DE GRAVEDAD en los nudos (Q) \n');
Q
% 4. ESTADO INICIAL (FUERZAS INTERNAS) DE LA ESTRUCTURA
% Matriz de rigidez de los elementos
for i=1:ne
k=rigidez1(L(i),an(i),EA(i),EIo(i),EIab(i,1),EIab(i,2),T(i));
ke(:,:,i)=k;
end
% Vector de colocación
% Matriz de rigidez de la estructura
K=zeros(gl,gl);
for i=1:ne
DK=zeros(gl,gl);
vc=VC(i,:);
vci=vc;
vc(vc==0)=[]; % Se eliminan las coordenadas iguales a (0).
a=find(vci); % Se ubica los casilleros de las coordenadas...
ki=ke(:,:,i); % distintas de (0).
n=ki(a,a);
DK(vc,vc)=n;
K=K+DK;
end
% Partición de la matriz de rigidez
a=glr+1:gl; % Coordenadas libres (a)
b=1:glr; % Coordendas restringidas (b)
Kaa=K(a,a);
Kab=K(a,b);
Kba=K(b,a);
Kbb=K(b,b);
% Vector de desplazamientos de la estructura (D)
Pa=Q(a); % Cargas nodales en coord. libres
97
Ra=R(a); % Cargas de empotramiento en coord. libres
Rb=R(b); % (Pa-Ra)= Matriz de fuerza de la estructura
Da=inv(Kaa)*(Pa-Ra); % Desplazamiento de coord. libres
Pb=Kba*Da+Rb; % Reacciones en coord. restringidas
D=zeros(gl,1);
D(a)=Da;
% Acumulación de fuerzas internas y desplazamiento
for i=1:ne
vc=VC(i,:); % Vector de colocación del elemento (i).
a=find(vc); % Posición de los valores distintos de (0).
vc(vc==0)=[]; % Vector de colocación sin (0).
di=zeros(6,1);
di(a)=D(vc); % Vector de desplazamientos del elemento.
p(:,i)=ke(:,:,i)*di+F(:,i); % Axial, cortante y momentos del elemento.
end
% Reorganización de fuerzas internas
for i=1:ne
if i<=nc
x(i,1:2)=p([2,5],i); % Solo axial
c(i,1:2)=p([1,4],i); % Solo Cortante
else
x(i,1:2)=p([1,4],i);
c(i,1:2)=p([2,5],i);
end
end
m=[p(3,:)',p(6,:)']; % Solo momentos
% 5. RESULTADOS
fprintf('Matriz con F. AXIALES en los extremos de los elementos (x) \n');
fprintf('Matriz con CORTANTES en los extremos de los elementos (c) \n');
fprintf('Matriz con MOMENTOS en los extremos de los elementos (m) \n');
d_puntal
function [vy,du,EAp]=d_puntal(np,nt,ne,nc,nv,nd,my)
%--------------------------------------------------------------------------
% El programa “d_puntal” dibuja el diagrama carga - deformación de
% un puntal equivalente con la metodología de carrillo (2007) y Osorio y
% Ramirez(2019).
%
% Por: Juan Carlos Quispe Ango
% UNIVERSIDAD CESAR VALLEJO
% Febrero/2022
%--------------------------------------------------------------------------
% [R]=deformacionpuntal(lm,hm,bc,hc,t,fc,fm,mp,Vpo)
%--------------------------------------------------------------------------
% ENTRADA
% np = Número de pisos
% nt = Número de tramos
% ne = Número de elementos
% nc = Número de columnas
% nv = Número de vigas
% nd = Número de puntales
% my = Momentos de fluencia
% RESULTADOS
% vy = Vector de capacidad a compresión
% du = Vector de desplazamiento ultimo
% EAp = Vector de rigidez elastica axial
%--------------------------------------------------------------------------
98
% 1. ENTRADAS
datos=input('Ingrese el archivo de datos: ');
datos1=xlsread(datos);
datos2=xlsread(datos);
datos2(1:(ne-nd),:)=[];
% 2. REORGANIZACIÓN
% Vector de colocación (x,y,z) para recoger secciones de columna
for i=1:np
n=(i-1)*nt+i:i*(nt+1)-1;
m=(i-1)*nt+i+1:i*(nt+1);
x((i-1)*nt+1:i*nt)=[n];
y((i-1)*nt+1:i*nt)=[m];
end
z=nc+1:nc+nv;
% Recogemos el promedio de las dimensiones de las columnas
hcp=zeros(2,nd); % Se promedia el peralte de las secciones
hcp(1,:)=datos1(x,7); % Se considera que el ancho es constante
hcp(2,:)=datos1(y,7);
hcp=mean(hcp); % Es vector fila
bca=zeros(2,nd);
bca(1,:)=datos1(x,6);
bca(2,:)=datos1(y,6);
bca=mean(bca); % Es vector fila
% Recogemos el menor (my) de los elementos que rodean al muro
me=min(my'); % Es vector fila
myp=zeros(3,nd);
myp(1,:)=me(x);
myp(2,:)=me(y);
myp(3,:)=me(z); % Menor (my=mp) de los confinamientos
myp=min(myp);
% Recogemos altura de las columnas
hmp=datos1(x,2);
% 2. PREPARACIÓN AL BUCLE
for j=1:nd
lm=datos1(nc+j,2)
hm=hmp(j)
bc=bca(j)
hc=hcp(j)
t=datos2(j,6)
fm=datos2(j,8)/500
mp=myp(j)
Vpo=datos2(j,17)
% 1. CÁLCULOS PREVIOS
Em=datos2(j,8) % T/m2
Ec=datos1(1,8) % T/m2
an=atan(hm/lm);
and=atand(hm/lm)
Ic=bc*hc^3/12 % m4
L=sqrt(lm^2+hm^2)
% 2. ANCHO DEL PUNTAL EQUIVALENTE
lam=(Em*t*sin(2*an)/(4*Ec*Ic*hm))^(1/4);
w=0.175*(lam*hm)^-0.4*L
% 3. TIPO DE FALLA DEL PUNTAL
m=8*mp/(fm*t*lm^2)
% 4. RESISTENCIA A COMPRESIÓN DEL PUNTAL
if m<1 % Falla por compresión
z=pi()/2*(4*Ec*Ic*hm/(Em*t*sin(2*an)))^(1/4);
Rc1=2*z*t*fm*sec(an)/3;
Rc2=fm*t*w;
R=min([Rc1,Rc2]);
99
else % Deslizamiento por cortante
tau=0.03*fm;
u=0.30;
R=tau*L*t/(1-u*(hm/lm));
end
% 5. RESISTENCIA A CORTANTE DE LA ALBAÑILERÍA
Vm=t*lm*sqrt(fm/102)*102/6
% 6. TABLA 7-9 FEMA 256
beta=Vpo/Vm
rel=lm/hm % Relación longitud/altura de la albañilería
if (beta<0.7) % (beta>=0.3)&&(beta<0.7)
if (rel>=0.5)&&(rel<1)
d=0.5-((0.1/0.5)*(rel-0.5));
elseif (rel>=1)&&(rel<2)
d=0.4-(((0.4-0.3)/(2-1))*(rel-1));
end
elseif (beta>=0.7)&&(beta<1.3)
if (rel>=0.5)&&(rel<1)
d=1-(((1-0.8)/(1-0.5))*(rel-0.5));
elseif (rel>=1)&&(rel<2)
d=0.8-(((0.8-0.6)/(2-1))*(rel-1));
end
elseif (beta>=1.3)
if (rel>=0.5)&&(rel<1)
d=1.5-(((1.5-1.2)/(1-0.5))*(rel-0.5));
elseif (rel>=1)&&(rel<2)
d=1.2-(((1.2-0.9)/(2-1))*(rel-1));
end
end
d
% 7. PUNTOS NOTABLES
wp=3*w % Rigidez elastica, intente calibrar con factor 3.
EA=Em*wp*t;
k=EA/L;
dy=R/(k*cos(an)) % Desplazamiento de fluencia horizontal
%du=d*hm*cos(an)/100
dup=(d+0.5*d)*hm/100; % (d+0.5d) expresión planteada por calibración
% 9. MATRIZ DE RESISTENCIA Y RIGIDEZ AXIAL SECANTE DE LOS PUNTALES
vy(j,1)=min([R/cos(an),Vm]); % Expresión planteada por calibración
EAp(j,1)=EA; % T/m
du(j,1)=dup;
% 8. DIAGRAMACIÓN ¡OK!
% figure(j)
% hold on % Sobreescribe
% plot([0,R/k,dup*cos(an)],[0,vy(j,1)*cos(an),vy(j,1)*cos(an)],'-
','LineWidth',3)
% title('CARGA - DEFORMACIÓN DEL PUNTAL')
% xlabel('Desplazamiento (m)')
% ylabel('Carga (T)')
% grid on
end
e_capacidad
function
[rotulas,Vy,Dx,dp]=e_capacidad(np,nt,ne,nc,nv,gl,glr,x,c,m,VC,my,mu,EIab,EIp,vy,
du,EAp)
100
%--------------------------------------------------------------------------
% El programa "e_capacidad", construye la curva de capacidad resistente de
% porticos planos de concreto armado y muros de albañilería confinada.
%
% Por: Juan Carlos Quispe Ango
% UNIVERSIDAD CESAR VALLEJO
% Febrero/2022
%--------------------------------------------------------------------------
% [rotulas]=resistencia(np,nt,ne,nc,gl,glr,x,c,m,VC,my,mu,EIab,EIp)
%--------------------------------------------------------------------------
% ENTRADA
%
% RESULTADOS
% Rotulas = Número de articulaciones generadas
%--------------------------------------------------------------------------
% 1. ENTRADA
datos=input('Ingrese el archivo de datos: ');
datos=xlsread(datos);
%H=input('¿Cual es la altura del edificio? (m): ');
L=datos(1:ne,2);
c1=datos(1:ne,3);
c2=datos(1:ne,4);
an=datos(1:ne,5);
EA=datos(1:ne,12);
GA=datos(1:ne,13);
T=datos(1:ne,16);
EIo=EIab(:,1); % Rigidez inicial elasticos en los extremos
% 2. REORGANIZACIÓN DE ENTRADAS
fprintf('PARAMETROS DE ENTRADA: \n');
fprintf('Momentos de fluencia (fisuramiento)(my). \n');
my % El momento de fluencia se reparte a los extremos
fprintf('Momentos de rotura (mu). \n');
mu % El momento de rotura se reparte a los extremos
if ne~=nc+nv
fprintf('Cortante maxima de los puntales (vu). \n');
vy
fprintf('Desplazamiento ultimo de los puntales (du). \n');
du
end
fprintf('Rigidez efectiva a flexión (EIab). \n');
EIab
fprintf('Rigidez en estado plastico (EIp). \n');
EIp % La rigidez inelastica se reparte a los extremos
if ne~=nc+nv
fprintf('Rigidez a carga axial(EA). \n');
EA(nc+nv+1:ne)=EAp % Se ingresa la rigidez plastica de los puntales
end
% 3. INCORPORACIÓN DE FUERZAS LATERALES AL VECTOR DE CARGAS EN LOS NUDOS
e=0.5; % Escala de fuerza lateral en el tope
Q=zeros(gl,1); % REVISAR EL TEMA DE INCORPORACIÓN
for i=1:np
vc=3*i*(nt+1)+1:3:3*(i+1)*(nt+1)-2;
P=i*e/np;
Q(vc)=P/(nt+1);
end
% 4. PREPARACIÓN A LA SECUENCIA REPETITIVA
fluencia=zeros(ne,2);
rotura=zeros(ne,2);
fisuramiento_albanileria=zeros(ne-nc-nv,1);
rotura_albanileria=zeros(ne-nc-nv,1);
101
rotulas=0;
dx=zeros(ne,2);
dc=zeros(ne,2);
dm=zeros(ne,2);
dr=zeros(ne-nc-nv,1);
v=0; % Para acumular cortante basal
d=0; % Para acumular desplazamiento en el tope.
j=0; % Contador para la captura de (v) y (d).
% 5. SECUENCIA REPETITIVA PUSHOVER HASTA EL COLAPSO
while rotulas<1
% Matriz de rigidez de los elementos
for i=1:ne
if (EIab(i,1)>0)&&(EIab(i,2)>0)
[k,t]=rigidez1(L(i),an(i),EA(i),EIo(i),EIab(i,1),EIab(i,2),T(i));
elseif (EIab(i,1)==0) && (EIab(i,2)>0)
[k]=rigidezb1(L(i),an(i),EA(i),EIo(i),EIab(i,1),EIab(i,2));
elseif (EIab(i,1)>0) && (EIab(i,2)==0)
[k]=rigidezb2(L(i),an(i),EA(i),EIo(i),EIab(i,1),EIab(i,2));
elseif (EIab(i,1)==0) && (EIab(i,2)==0)
[k]=rigidezb3(L(i),an(i),EA(i),EIo(i),EIab(i,1),EIab(i,2));
end
ke(:,:,i)=k;
te(:,:,i)=t;
end
% Vector de colocación
% Matriz de rigidez de la estructura
K=zeros(gl,gl);
for i=1:ne
DK=zeros(gl,gl);
vc=VC(i,:);
vci=vc;
vc(vc==0)=[]; % Se eliminan las coordenadas iguales a (0).
a=find(vci); % Se ubica los casilleros de las coordenadas
ki=ke(:,:,i); % distintas de (0).
n=ki(a,a);
DK(vc,vc)=n;
K=K+DK;
end
% Partición de la matriz de rigidez
a=glr+1:gl; % Coordenadas libres (a)
b=1:glr; % Coordendas restringidas (b)
Kaa=K(a,a);
Kab=K(a,b);
Kba=K(b,a);
Kbb=K(b,b);
% Vector de desplazamientos de la estructura (D)
Pa=Q(a); % Cargas nodales en coord. libres
Da=inv(Kaa)*Pa; % Desplazamiento de coord. libres.
Pb=Kba*Da; % Reacciones en coord. restringidas.
D=zeros(gl,1);
D(a)=Da;
% Acumulación de fuerzas internas y desplazamiento
for i=1:ne
vc=VC(i,:); % Vector de colocación del elemento (i).
a=find(vc); % Posición de los valores distintos de (0).
vc(vc==0)=[]; % Vector de colocación sin (0).
dg=zeros(6,1);
dg(a)=D(vc); % Vector de desplazamientos del elemento.
p(:,i)=ke(:,:,i)*dg; % Axial, cortante y momentos del elemento Glg
de(1:6,i)=dg; % Para el comportamiento de los puntales
102
end
% Reorganización de desplazamientos
for i=1:ne
if i<=nc
idx(i,1:2)=de([2,5],i); % Solo axial
idc(i,1:2)=de([1,4],i); % Solo Cortante
else
idx(i,1:2)=de([1,4],i);
idc(i,1:2)=de([2,5],i);
end
end
idm=[de(3,:)',de(6,:)']; % Solo momentos
dx=dx+idx; % Acumulación de axiales
dc=dc+idc; % Acumulación de cortantes
dm=dm+idm; % Acumulación de momentos
% Reorganización de fuerzas internas
for i=1:ne
if i<=nc
ix(i,1:2)=p([2,5],i); % Solo axial
ic(i,1:2)=p([1,4],i); % Solo Cortante
else
ix(i,1:2)=p([1,4],i);
ic(i,1:2)=p([2,5],i);
end
end
im=[p(3,:)',p(6,:)']; % Solo momentos
x=x+ix; % Acumulación de axiales
c=c+ic; % Acumulación de cortantes
m=m+im; % Acumulación de momentos
% Selección de axiales en los puntales
if ne~=nc+nv
xp=x(nc+nv+1:ne,2);
end
% Selección de desplazamientos horizontales de los puntales
if ne~=nc+nv
dp=dx(nc+nv+1:ne,2); % Desplazamiento absoluto
dp(nt+1:nt*np)=dp(nt+1:nt*np)-dp(1:nt*(np-1)); % Despl. relativo
end
% Acumulación de la cortante basal y el desplazamiento en el tope
for i=1:np
vc=3*i*(nt+1)+1:3:3*(i+1)*(nt+1)-2;
iv(i,:)=Q(vc);
end
v=v+sum(sum(iv));
d=d+D(gl-2);
% Pushover controlado OK!
% Captura de fallas
pos1=find(abs(m(:,:))>my(:,:)); % Se encuentra la posición de ...
fluencia(pos1)=1; % ... los m mayores a my.
pos2=find(abs(m(:,:))>mu(:,:));
rotura(pos2)=1;
if ne~=nc+nv % Para comportamiento de puntales
pos3=find(xp(:)>vy(:));
fisuramiento_albanileria(pos3)=1;
pos4=find(dp(:)>du(:));
rotura_albanileria(pos4)=1;
end
rotulas1=sum(sum(rotura));
rotulas2=sum(sum(rotura_albanileria));
rotulas=rotulas1+rotulas2;
103
% Captura de cortante y desplazamiento, cambio de estado de los
% elementos e impresión en pantalla.
if pos1~=0
j=j+1;
fprintf('ESTADO DE FALLA %2d (ne,[i,j]): \n',j);
fluencia
rotura
fprintf('Fuerza axial: \n');
x
dx
fprintf('Fuerzas cortantes: \n');
c
dc
fprintf('momentos flectores: \n');
m
dm
fprintf('Estado de la resistencia: \n');
my(pos1)=mu(pos1) % Se reemplazan los (my) por (mu).
%my(pos2)=0 % Se reemplazan los (mu) por (0).
fprintf('Estado de la rigidez a flexión: \n');
EIab(pos1)=EIp(pos1) % Se reemplazan las rigideces elas. por inel
%EIab(pos2)=0 % Se reemplazan las rigideces inel. por (0).
Vy(j+1)=v; % Captura de cortante y desplazamiento
Dx(j+1)=d;
end
% Captura de cortante y desplazamiento, cambio de estado de los
% puntales e impresión en pantalla.
if ne~=nc+nv % Para comportamiento de puntales
if pos3~=0
j=j+1;
fprintf('ESTADO DE FALLA %2d (ne,[i,j]): \n',j);
fisuramiento_albanileria
rotura_albanileria
fprintf('Fuerza axial: \n');
x
dx
fprintf('Fuerzas cortantes: \n');
c
dc
fprintf('momentos flectores: \n');
m
dm
fprintf('Estado de la resistencia: \n');
vy(pos3)=vy(pos3);
fprintf('Estado de la rigidez a axial: \n');
ea=EA(nc+nv+1:ne); % Se extrae EA de los puntales
ea(pos3)=0; % Se elimina la rigidez axial puntales
EA(nc+nv+1:ne)=ea;
Vy(j+1)=v; % Captura de cortante y desplazamiento
Dx(j+1)=d;
end
end
end
% 7. DIAGRAMACIÓN DE LA CURVA DE CAPACIDAD
hold on % Sobreescribe
plot(Dx,Vy,'-','LineWidth',3)
title('CURVA DE CAPACIDAD')
xlabel('d (m)')
ylabel('V (T)')
grid on
104
Anexo N° 4. Fichas de recolección de datos
105
106
107